A continuum‐discrete model using Darcy's law: formulation and verification

This paper presents a numerical scheme for fluid‐particle coupled discrete element method (DEM), which is based on poro‐elasticity. The motion of the particles is resolved by means of DEM. While within the proposition of Darcian regime, the fluid is assumed as a continuum phase on a Eulerian mesh, and the continuity equation on the fluid mesh for a compressible fluid is solved using the FEM. Analytical solutions of traditional soil mechanics examples, such as the isotropic compression and one‐dimensional upward seepage flow, were used to validate the proposed algorithm quantitatively. The numerical results showed very good agreement with the analytical solutions, which show the correctness of this algorithm. Sensitivity studies on the effect of some influential factors of the coupling scheme such as pore fluid bulk modulus, volumetric strain calculation, and fluid mesh size were performed to display the accuracy, efficiency, and robustness of the numerical algorithm. It is revealed that the pore fluid bulk modulus is a critical parameter that can affect the accuracy of the results. Because of the iterative coupling scheme of these algorithms, high value of fluid bulk modulus can result in instability and consequently reduction in the maximum possible time‐step. Furthermore, the increase of the fluid mesh size reduces the accuracy of the calculated pore pressure. This study enhances our current understanding of the capacity of fluid‐particle coupled DEM to simulate the mechanical behavior of saturated granular materials. Copyright © 2014 John Wiley & Sons, Ltd.     

A computer method based on simulated annealing to identify aquifer parameters using pumping‐test data

Conventional graphical or computer methods for identifying aquifer parameters have their own inevitable limitations. This paper proposes a computer method based on a drawdown model and a heuristic approach of simulated annealing (SA) to determine the best‐fit aquifer parameters of the confined and unconfined aquifer systems. The drawdown model for the confined aquifer is the Theis solution and the unconfined aquifer is the Neuman solution. The estimated results of proposed method have better accuracy than those of the graphical methods and agree well with those of the computer methods based on the extended Kalman filter and Newton's method. Finally, the sensitivity analyses for the control parameters of SA indicate that the proposed method is very robust and stable in parameter identification procedures. Copyright © 2007 John Wiley & Sons, Ltd.     

Stochastic multi-site generation of daily rainfall occurrence in south Florida

This paper presents a stochastic model to generate daily rainfall occurrences at multiple gauging stations in south Florida. The model developed in this study is a space–time model that takes into account the spatial as well as temporal dependences of daily rainfall occurrence based on a chain-dependent process. In the model, a Markovian method was used to represent the temporal dependence of daily rainfall occurrence and a direct acyclic graph (DAG) method was introduced to encode the spatial dependence of daily rainfall occurrences among gauging stations. The DAG method provides an optimal sequence of generation by maximizing the spatial dependence index of daily rainfall occurrences over the region. The proposed space–time model shows more promising performance in generating rainfall occurrences in time and space than the conventional Markov type model. The space–time model well represents the temporal as well as the spatial dependence of daily rainfall occurrences, which can reduce the complexity in the generation of daily rainfall amounts.     

An analytical study on the second-order resonance system of tide in an ideal partially-enclosed bay

In the numerical studies of a real tide M4 resonance system, the Xiangshan Port which is a partially-closed bay, Dong et al. [1999. Acta Oceanologica Sinica, 21 （3）： 1-6] found the interesting phenomenon that the advection plays an important role in inhibiting the growth of the amplitude of the tidal second-order resonance response （M4）. This result is contrary to the general traditional ideas for a non-resonance system. How this phenomenon is interpreted and what internal mechanism is behind the phenomenon are the main focuses of this study. The followings are examined： （1） the dynamic features of a second-order resonance system of tide; （2） the dominating factors on the second-order resonance responses; （3） the effects of both the friction and the advection on the second-order resonance responses; and （4） their roles in dominating the second-order resonance response and internal mechanisms by using the analytical methods. The respective results show that： （1） Both the bottom friction and the advection play significant roles in dominating the magnitude of the amplitude of the second-order resonance responses; （2） the effect of the friction on the second-order resonance response depends on the distribution ratio of the work-done of the system to friction force exhausted into between the damping of the first-order system and the inner excitation of the second-order system; （3） the advection plays a positive role in increasing the amplitude of the second-order non-resonance response in the second order non-resonance of tide; （4） in a second-order resonance system of tide, the effect of the advection may be either to increase or to decrease the amplitudes of the second-order resonance responses of tide, which depends on the distribution ratio mentioned above.     

Weather conditions associated with the passage of precipitation type transition regions over eastern Newfoundland

Abstract

The passage of a winter storm is accompanied by changes in many surface and near‐surface parameters including temperature, humidity, wind, pressure, precipitation rate and type, cloud base height, visibility and accretion. These parameters were measured in association with the passage of precipitation‐type transitions over Newfoundland during the Canadian Atlantic Storms Program II field experiment. Three simple summaries of the observed weather events were developed. These summaries depend on the observed large‐scale synoptic conditions, which include warm fronts, a cold front and a trough.     

Modeling crack propagation path of anisotropic rocks using boundary element method

In a cracked material, the stress intensity factors (SIFs) at the crack tips, which govern the crack propagation and are associated with the strength of the material, are strongly affected by the crack inclination angle and the orientation with respect to the principal direction of anisotropy. In this paper, a formulation of the boundary element method (BEM), based on the relative displacements of the crack tip, is used to determine the mixed‐mode SIFs of isotropic and anisotropic rocks. Numerical examples of the application of the formulation for different crack inclination angles, crack lengths, and degree of material anisotropy are presented. Furthermore, the BEM formulation combined with the maximum circumferential stress criterion is adopted to predict the crack initiation angles and simulate the crack propagation paths. The propagation path in cracked straight through Brazilian disc specimen is numerically predicted and the results of numerical and experimental data compared with the actual laboratory observations. Good agreement is found between the two approaches. The proposed BEM formulation is therefore suitable to simulate the process of crack propagation. Copyright © 2008 John Wiley & Sons, Ltd.     

Investigation of the impacts of vegetation distribution and evaporative cooling on synthetic urban daytime climate using a coupled LES—LSM model

Local flow properties and regional weather or climate are strongly affected by land‐atmosphere interactions of momentum and scalars within the daytime convective boundary layer (CBL). In this study, we investigate the impact of green space scale on the daytime atmospheric boundary layer (ABL) over a synthetic urban domain using a recently developed large‐eddy simulation‐land surface model (LES–LSM) framework. With the use of realistic soundings as initial conditions, a series of numerical experiments over synthetic urban surfaces with varied scale of vegetated area is performed. Simulated micrometeorological properties, surface fluxes, basic CBL characteristics, and cloud distribution are analysed. The results show reference‐level air potential temperature and specific humidity as well as surface fluxes over green space are significantly affected by the scale of green space in the urban domain. The surface organization due to vegetated area scale also has impacts on horizontally averaged scalar and momentum profiles; however, the magnitude in this study is smaller than the results of a previous study using a set of offline surface fluxes as the lower boundary condition for LES. In addition, even though this study only performs a daytime diurnal cycle, the impact of green space scale on cloud distribution in simulations is significant. The cases with more organized green space yield lower‐elevated cumulus cloud and larger‐cloud cover fraction, which impacts the energy budget at the top of boundary layer and, in turn, could lead to additional surface cooling with respect to longer‐term weather and climate. Copyright © 2010 John Wiley & Sons, Ltd.     

Neural network modelling for mean velocity and turbulence intensities of steep channel flows

The main purpose of this study is to evaluate the potential of simulating the profiles of the mean velocity and turbulence intensities for the steep open channel flows over a smooth boundary using artificial neural networks. In a laboratory flume, turbulent flow conditions were measured using a fibre‐optic laser doppler velocimeter (FLDV). One thousand and sixty‐four data sets were collected for different slopes and aspect ratios at different locations. These data sets were randomly split into two subsets, i.e. training and validation sets. The multi‐layer functional link network (MFLN) was used to construct the simulation model based on the training data. The constructed MFLN models can almost perfectly simulate the velocity profile and turbulence intensity. The values of correlation coefficient (γ) are close to one and the values of root mean square error (RMSE) are close to zero in all conditions. The results demonstrate that the MFLN can precisely simulate the velocity profiles, while the log law and Reynolds stress model (RSM) are less effective when used to simulate the velocity profiles close to the side wall. The simulated longitudinal turbulence intensities yielded by the MFLN were also fairly consistent with the measured data, while the simulated vertical turbulence intensities by the RSM were not consistent with the measured data. Copyright © 2007 John Wiley & Sons, Ltd.     

Mesozoic high-Ba–Sr granitoids from North China: geochemical characteristics and geological implications

Abstract Mesozoic magmatism occurred extensively in the North China block (NCB) and the Dabie–Sulu orogen (DSO) post-dating the North–South China collision, resulting in abundant intrusive and volcanic rocks ranging from basic to acidic compositions. The intermediate-acidic intrusive rocks can be grouped into two types, namely high-Ba–Sr granitoids and low-Ba–Sr granitoids that both have distinct geochemical characteristics. The high-Ba–Sr granitoids are similar in most of the incompatible trace element systematics to the associated basic rocks, which probably originated from melting of subcontinental lithospheric mantle, indicating significant mantle contributions to them. Geochemical similarities are observed between the basic rocks from the NCB and DSO, implying a regional-scale magma-generating mechanism and that mantle enrichment beneath the DSO was independent from the Triassic deep continental subduction in the region. We therefore interpret that the Mesozoic magmatism resulted from delamination of the ancient lithospheric mantle beneath the eastern part of North China.