Experimental and computational validation and verification of the Stokes–Darcy and continuum pipe flow models for karst aquifers with dual porosity structure     

Xiaolong Hu  Xiaoming Wang  Max Gunzburger  Fei Hua  Yanzhao Cao 《水文研究》2012,26(13):2031-2040

In our previous study, we developed the Stokes–Darcy (SD) model was developed for flow in a karst aquifer with a conduit bedded in matrix, and the Beavers–Joseph (BJ) condition was used to describe the matrix–conduit interface. We also studied the mathematical well‐posedness of a coupled continuum pipe flow (CCPF) model as well as convergence rates of its finite element approximation. In this study, to compare the SD model with the CCPF model, we used numerical analyses to validate finite element discretisation methods for the two models. Using computational experiments, simulation codes implementing the finite element discretisations are then verified. Further model validation studies are based on the results of laboratory experiments. Comparing the results of computer simulations and experiments, we concluded that the SD model with the Beavers–Joseph interface condition is a valid model for conduit–matrix systems. On the other hand, the CCPF model with the value of the exchange parameter chosen within the range suggested in the literature perhaps does not result in good agreement with experimental observations. We then examined the sensitivity of the CCPF model with respect to the exchange parameter, concluding that, as has previously been noted, the model is highly sensitive for small values of the exchange parameter. However, for larger values, the model becomes less sensitive and, more important, also produces results that are in better agreement with experimental observations. This suggests that the CCPF model may also produce accurate simulation results, if one chooses larger values of the exchange parameter than those suggested in the literature. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary     

Mohammad Reza Maleki Javan  Asadollah Noorzad  Manouchehr Latifi Namin 《国际地质力学数值与分析法杂志》2008,32(6):681-699

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

Assessing uncertainties in dust emission in the Aral Sea region caused by meteorological fields predicted with a mesoscale model     

Kremena Darmenova  Irina N. Sokolik   《Global and Planetary Change》2007,56(3-4):297

This study investigates how the choice of the planetary boundary layer (PBL) parameterization and dust emission scheme affects the prediction of dust entrainment simulated with a regional mesoscale model. For this analysis we consider a representative dust episode which occurred on April 2001 in the Aral Sea region. The meteorological fields were simulated using the PSU/NCAR MM5 modeling system considering two different boundary layer parameterizations. In each case, emitted dust fluxes were computed off-line by incorporating MM5 meteorological fields into the dust module DuMo. Several dust emission schemes with a prescribed erodible fraction and fixed threshold wind speed were the subject of our analysis. Implications to assessment of the anthropogenic fraction of dust emitted in the Aral region were investigated by conducting the full, half, and no lake modeling experiments.Our results show that the discrepancies in dust fluxes between the two different PBLs are much higher compared to the discrepancy associated with the use of considered dust production schemes. Furthermore, the choice of the PBL affects the timing and duration of a modeled dust event. We demonstrate that different combinations of the PBL parameterization and wind- or friction velocity-driven dust emission schemes can result in up to about a 50% difference in predicted dust mass caused by the Aral Sea desiccation. We found that the drying-up of the Aral cannot only affect the dust emission by expanding the source area, but also by affecting atmospheric characteristics, especially winds. These competitive factors add further complexity to quantification of the anthropogenic dust fraction in the region.  相似文献   

一种以能量平衡为基础的干旱指数   总被引：5，自引：0，他引：5  

唐登银 《地理研究》1987,6(2):21-31

本文以能量平衡公式为基础,根据实际蒸发与潜在蒸发的关系依赖于土壤水分含量的事实,导出一种表达干湿状况的指标——土壤水分干旱指数。  相似文献   

Spatial scale effect on seasonal streamflows in permafrost catchments on the Qinghai–Tibet Plateau     

Genxu Wang  Guangsheng Liu  Lin'an Liu 《水文研究》2012,26(7):973-984

The scale issue is of central concern in hydrological processes to understand the potential upscaling or downscaling methodologies, and to develop models for scaling the dominant processes at different scales and in different environments. In this study, a typical permafrost watershed in the Qinghai‐Tibet Plateau was selected. Its hydrological processes were monitored for 4 years from 2004 to 2008, measuring the effects of freezing and thawing depth of active soil layers on runoff processes. To identify the nature and cause of variation in the runoff response in different size catchments, catchments ranging from 1·07 to 112 km² were identified in the watershed. The results indicated that the variation of runoff coefficients showed a ‘V’ shape with increasing catchment size during the spring and autumn seasons, when the active soil was subjected to thawing or freezing processes. A two‐stage method was proposed to create runoff scaling models to indicate the effects of scale on runoff processes. In summer, the scaling transition model followed an exponential function for mean daily discharge, whereas the scaling model for flood flow exhibited a linear function. In autumn, the runoff process transition across multiple scales followed an exponential function with air temperature as the driving factor. These scaling models demonstrate relatively high simulation efficiency and precision, and provide a practical way for upscaling or downscaling runoff processes in a medium‐size permafrost watershed. For permafrost catchments of this scale, the results show that the synergistic effect of scale and vegetation cover is an important driving factor in the runoff response. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

Tests and applications of an approach to absorbing reflected waves towards incident boundary     

张洪生  王炎  许春辉  商辉  于小伟 《中国海洋工程》2013,27(6):703-718

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.  相似文献   

Under‐ice salinity and stable isotope distribution of Saroma‐ko Lagoon,Hokkaido, northern Japan     

Maki Morimoto  Kyoko Kawanobe  Osamu Abe  Takayoshi Kawai  Toshiyuki Kawamura  Kunio Shirasawa 《水文研究》2010,24(7):904-916

In winter, lakes and lagoons at high altitudes or high latitudes have interesting hydrological cycles that differ from those in other seasons or in other regions, because water surfaces are covered with ice. Hydrological balances of lakes and lagoons are complex dynamic systems, and to elucidate them, isotopic tracers of water have been used as effective tools along with observations of precipitation, evaporation, inflows, and outflows. Here, to understand hydrological processes during freezing periods in the brackish Saroma‐ko Lagoon, Hokkaido, northern Japan, we examined horizontal and vertical distributions of salinity and isotope compositions of lagoon water and ice in 2005 and 2006. Horizontal and vertical gradients of salinity and isotope compositions were observed from the river mouth to the sea channel, and factors determining these distributions were considered. The mixing of freshwater and seawater and a freezing effect were presumed to be factors in relationships between salinity and isotopes and in relationships between surface waters and ice just above the water. A simple box model for water balance was constructed based on these putative factors to reproduce the distributions of salinity and isotope compositions of surface waters and ice. An evaluation of the model revealed that this hydrological system is controlled primarily by horizontal advection of the epilimnion, freshwater influx, and the ice growth rate. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Importance of footing–soil separation on dynamic stiffness of piled embedded footings     

L. A. Padrón  G. E. Mylonakis  D. E. Beskos 《国际地质力学数值与分析法杂志》2009,33(11):1439-1448

Different phenomena such as soil consolidation, erosion, and scour beneath an embedded footing supported on piles may lead to loss of contact between soil and the pile cap underside. The importance of this separation on the dynamic stiffness and damping of the foundation is assessed in this work. To this end, a numerical parametric analysis in the frequency domain is performed using a rigorous three‐dimensional elastodynamic boundary element–finite element coupling scheme. Dimensionless plots relating dynamic stiffness functions computed with and without separation effects are presented for different pile–soil configurations. Vertical, horizontal and rocking modes of oscillation are analyzed for a wide range of dimensionless frequencies. It is shown that the importance of separation is negligible for frequencies below those for which dynamic pile group effects start to become apparent. Redistribution of stiffness contributions between piles and footing is also addressed. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

Evaluating the stress intensity factors of anisotropic bimaterials using boundary element method     

Chao‐Shi Chen  Chien‐Chung Ke  Chia‐Huei Tu 《国际地质力学数值与分析法杂志》2008,32(11):1341-1364

This paper presents a boundary element method (BEM) procedure for a linear elastic fracture mechanics analysis in two‐dimensional anisotropic bimaterials. In this formulation, a displacement integral equation is only collocated on the uncracked boundary, and a traction integral equation is only collocated on one side of the crack surface. A fundamental solution (Green's function) for anisotropic bimaterials is also derived and implemented into the boundary integral formulation so that except for the interfacial crack part, the discretization along the interface can be avoided. A special crack‐tip element is introduced to capture the exact crack‐tip behavior. A computer program using FORTRAN has been developed to effectively calculate the stress intensity factors of an anisotropic bimaterial. This BEM program has been verified to have a good accuracy with previous studies. In addition, a central cracked bimaterial Brazilian specimen constituting cement and gypsum is prepared to conduct the Brazilian test under diametral loading. The result shows that the numerical analysis can predict relatively well the direction of crack initiation and the path of crack propagation. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

船行波对H-ADCP 影响分析研究     

《海洋技术学报》2024,(1)

作为当今世界上最先进的测流仪器，水平式声学多普勒流速剖面仪（Horizontal Acoustic Doppler Current Profiler，H-ADCP） 因其精度高、速度快等突出优势在河流中的应用越来越广泛，但在航运业发达地区，由于受船行波的影响，很难得到稳定可靠的流速数据，阻碍了H-ADCP在平原河网地区的应用和推广。目前技术上很难做到完全消除船行波的影响，本文将H-ADCP应用在平原河网地区对流速进行测量，并分析了船行波对测量精度的影响。从H-ADCP设备选型、测量和安装方式、测量预处理、数据后处理等多方面入手，可以有效消减船行波对H-ADCP 的影响，从而提高测流精度。同时，通过实验测量断面一天内平均流速变化情况，得到受船行波影响后的流速波形，然后采用数据处理手段，验证了所提方法的可行性。  相似文献