首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 234 毫秒
1.
The objectives of this study are to investigate the third order accuracy and linear stability of the lattice Boltzmann method (LBM) with the two-relaxation-time collision operator (LTRT) for the advection–diffusion equation (ADE) and compare the LTRT model with the single-relaxation-time (LBGK) model. While the LBGK has been used extensively, the LTRT appears to be a more flexible model because it uses two relaxation times. The extra relaxation time can be used to improve solution accuracy and/or stability. This study conducts a third order Chapman–Enskog expansion on the LTRT to recover the macroscopic differential equations up to the third order. The dependency of third order terms on the relaxation times is obtained for different types of equilibrium distribution functions (EDFs) and lattices. By selecting proper relaxation times, the numerical dispersion can be significantly reduced. Furthermore, to improve solution accuracy, this study introduces pseudo-velocities to develop new EDFs to reduce the second order numerical diffusion. This study also derives stability domains based on the lattice Peclet number and Courant number for different types of lattices, EDFs and different values of relaxation times, while conducting linear stability analysis on the LTRT. Numerical examples demonstrate the improvement of the LTRT solution accuracy and stability by selecting proper relaxation times, lattice Peclet number and Courant number.  相似文献   

2.
A mirror-image method is proposed in this paper to solve the boundary conditions in the lattice Boltzmann model proposed by Zhang et al. [Adv. Water Resour. 25 (2002) 1] for the advection and anisotropic dispersion of solute transport in porous media. Three types of boundary are considered: prescribed concentration boundary, prescribed flux boundary and prescribed concentration-gradient boundary. The accuracy of the proposed method is verified against benchmark problems and finite difference method.  相似文献   

3.
如何有效压制数值频散是有限差分正演模拟研究中的关键问题之一.近年来,许多学者对二阶声波方程的差分算子开展了大量的优化工作,在压制频散方面取得不错的效果.一阶压强-速度方程广泛用于研究地震波在地下变密度模型中传播规律,目前针对一阶方程的优化工作大多只是在空间差分算子上展开.本文在前人研究的基础上,推导出一阶声波方程中压强场与偏振速度场之间的解析关系,据此在传统交错网格基础上给出一种高精度的显式时间递推格式,该递推格式将时间差分与空间差分算子结合在一起,并采用共轭梯度法得到精确时间递推匹配系数,实现时空差分算子的同时优化.在编程实现算法的基础上,通过频散分析与三个典型模型测试表明:本文方法能够较为有效地压制时间频散与空间频散,提高数值计算精度;同时对复杂模型也有很好适用性.  相似文献   

4.
This study formulates and analyzes continuous time random walk (CTRW) models in radial flow geometries for the quantification of non-local solute transport induced by heterogeneous flow distributions and by mobile–immobile mass transfer processes. To this end we derive a general CTRW framework in radial coordinates starting from the random walk equations for radial particle positions and times. The particle density, or solute concentration is governed by a non-local radial advection–dispersion equation (ADE). Unlike in CTRWs for uniform flow scenarios, particle transition times here depend on the radial particle position, which renders the CTRW non-stationary. As a consequence, the memory kernel characterizing the non-local ADE, is radially dependent. Based on this general formulation, we derive radial CTRW implementations that (i) emulate non-local radial transport due to heterogeneous advection, (ii) model multirate mass transfer (MRMT) between mobile and immobile continua, and (iii) quantify both heterogeneous advection in a mobile region and mass transfer between mobile and immobile regions. The expected solute breakthrough behavior is studied using numerical random walk particle tracking simulations. This behavior is analyzed by explicit analytical expressions for the asymptotic solute breakthrough curves. We observe clear power-law tails of the solute breakthrough for broad (power-law) distributions of particle transit times (heterogeneous advection) and particle trapping times (MRMT model). The combined model displays two distinct time regimes. An intermediate regime, in which the solute breakthrough is dominated by the particle transit times in the mobile zones, and a late time regime that is governed by the distribution of particle trapping times in immobile zones. These radial CTRW formulations allow for the identification of heterogeneous advection and mobile-immobile processes as drivers of anomalous transport, under conditions relevant for field tracer tests.  相似文献   

5.
A three-dimensional stochastic Lagrangian particle tracking sediment transport model is developed to solve the discrete advection-dispersion equation using a combination of empirical dispersion equations.The performance of three widely-used longitudinal dispersion coefficient equations was examined to select one of them as the primary dispersion equation term in the developed model. Also, a conditional empirical equation was used to consider the effect of vertical dispersion term in top layers n...  相似文献   

6.
Time nonlocal transport models such as the time fractional advection‐dispersion equation (t‐fADE) were proposed to capture well‐documented non‐Fickian dynamics for conservative solutes transport in heterogeneous media, with the underlying assumption that the time nonlocality (which means that the current concentration change is affected by previous concentration load) embedded in the physical models can release the effective dispersion coefficient from scale dependency. This assumption, however, has never been systematically examined using real data. This study fills this historical knowledge gap by capturing non‐Fickian transport (likely due to solute retention) documented in the literature (Huang et al. 1995) and observed in our laboratory from small to intermediate spatial scale using the promising, tempered t‐fADE model. Fitting exercises show that the effective dispersion coefficient in the t‐fADE, although differing subtly from the dispersion coefficient in the standard advection‐dispersion equation, increases nonlinearly with the travel distance (varying from 0.5 to 12 m) for both heterogeneous and macroscopically homogeneous sand columns. Further analysis reveals that, while solute retention in relatively immobile zones can be efficiently captured by the time nonlocal parameters in the t‐fADE, the motion‐independent solute movement in the mobile zone is affected by the spatial evolution of local velocities in the host medium, resulting in a scale‐dependent dispersion coefficient. The same result may be found for the other standard time nonlocal transport models that separate solute retention and jumps (i.e., displacement). Therefore, the t‐fADE with a constant dispersion coefficient cannot capture scale‐dependent dispersion in saturated porous media, challenging the application for stochastic hydrogeology methods in quantifying real‐world, preasymptotic transport. Hence improvements on time nonlocal models using, for example, the novel subordination approach are necessary to incorporate the spatial evolution of local velocities without adding cumbersome parameters.  相似文献   

7.
Groundwater contaminant transport processes are usually simulated by the finite difference (FDM) or finite element methods (FEM). However, they are susceptible to numerical dispersion for advection‐dominated transport. In this study, a numerical dispersion‐free coupled flow and transport model is developed by combining the analytic element method (AEM) with random walk particle tracking (RWPT). As AEM produces continuous velocity distribution over the entire aquifer domain, it is more suitable for RWPT than FDM/finite element methods. Using the AEM solutions, RWPT tracks all the particles in a vectorized manner, thereby improving the computational efficiency. The present model performs a convolution integral of the response of an impulse contaminant injection to generate concentration distributions due to a permanent contaminant source. The RWPT model is validated with an available analytical solution and compared to an FDM solution, the RWPT model more accurately replicates the analytical solution. Further, the coupled AEM‐RWPT model has been applied to simulate the flow and transport in hypothetical and field aquifer problems. The results are compared with the FDM solutions and found to be satisfactory. The results demonstrate the efficacy of the proposed method.  相似文献   

8.
This study develops a lattice Boltzmann method (LBM) with a two-relaxation-time collision operator (LTRT) to solve saltwater intrusion problems. A directional-speed-of-sound (DSS) technique is introduced to take into account the hydraulic conductivity heterogeneity and discontinuity, as well as the velocity-dependent dispersion coefficient. The forcing terms in the LTRT model are customized in order to recover the density-dependent groundwater flow and mass transport equations. Using the LTRT with the squared DSS achieves at least second-order accuracy. The LTRT results are verified with Henry’s analytical solution as well as compared with several numerical examples and modified Henry problems that consider heterogeneous hydraulic conductivity and velocity-dependent dispersion. The numerical results show good agreement with the Henry analytical solution and with the numerical solutions obtained by other numerical methods.  相似文献   

9.
This article presents a new public domain tool for generalized Lagrangian particle tracking in rivers. The approach can be applied with a variety of two- and three-dimensional flow solvers. Particle advection by the flow is incorporated using flow fields from the chosen solver assuming particles follow the Reynolds-averaged flow, although some other simple passive and active particle behaviors are also treated. Turbulence effects are treated using a random walk algorithm with spatial step lengths randomly chosen from Gaussian distributions characterized by the diffusivity from the flow solver. Our work extends this concept to a general framework that is solver and coordinate system independent to allow easy comparisons between differing flow treatments. To better treat problems where detailed information is required in specific regions, the approach includes novel cloning and colligation algorithms which enhance local resolution at modest computational expense. We also provide tools for computing local concentrations and total exposure over a user-specified time interval. Several examples of predictions are provided to illustrate applications of the technique, including examination of the role of curvature-driven secondary flows, storage in lateral separation eddies, treatment of larval drift, treatment of fuel spill dispersion, river-floodplain connections, and sedimentation in floodplain ponds by tie channel connections. We also demonstrate that the model can reproduce analytically derived concentration profiles for simple diffusivities. These examples show that the Lagrangian particle tracking approach and the extensions proposed here are broadly applicable and viable for treating difficult river problems with multiple temporal and spatial scales. The examples also illustrate the utility of the cloning/colligation extensions and show how these can decrease the computational effort required on problems where high local resolution is required. Enhancement of the tools and even broader applicability can be achieved through the inclusion of multiple particle populations and particle–particle interactions.  相似文献   

10.
Soil vapour extraction (SVE) is a common remediation technique for cleaning up unsaturated soils contaminated by volatile organic compounds (VOCs). Analytical solutions, which result from simple mathematical models, can allow the fast approximation of the time‐dependent effluent concentration and the gaining of insight into the processes that take place during soil remediation. Deriving the analytical solutions to advection–dispersion equations that simultaneously take into account the mechanical dispersion and molecular diffusion is very difficult because of the variable dependence of governing equations' coefficients. In this study, we first present two simplified analytical solutions that only consider mechanical dispersion or molecular diffusion. The two developed analytical solutions are compared with the numerical solution that simultaneously considers both mechanical dispersion and molecular diffusion to examine the applicability of the two simplified analytical solutions and distinguishes the individual contribution of the mechanical dispersion and molecular diffusion to total VOCs transport in an SVE system. Results show that dispersion plays an important role during SVE decontamination and neither the diffusion‐dominated solution nor the dispersion‐dominated solution can agree well with the numerical solution when both mechanical dispersion and molecular diffusion have significant contributions to the total VOCs transport flux. A composite analytical solution that linearly couples the diffusion‐ and dispersion‐dominated analytical solutions, which is proposed herein to eliminate the discrepancy between the analytical solutions and the numerical solution. Results indicate that the proposed composite analytical solution agrees well with the numerical solution and is an effective tool for quickly and accurately evaluating the time‐dependent effluent concentration for parameters of the different ranges of interest in an SVE remedial system. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

11.
Understanding the transport of chemical components in porous media is fundamentally important to many reservoir processes such as contaminant transport and reactive flows involved in CO2 sequestration. Carbonate rocks in particular present difficulties for pore-scale simulations because they contain large amounts of sub-micron porosity. In this work, we introduce a new hybrid simulation model to calculate hydrodynamic dispersion in pore-scale images of real porous media and use this to elucidate the origins and behaviour of stagnant zones arising in transport simulations using micro-CT images of carbonates. For this purpose a stochastic particle model for simulating the transport of a solute is coupled to a Lattice-Boltzmann algorithm to calculate the flow field. The particle method incorporates second order spatial and temporal resolution to resolve finer features of the domain. We demonstrate how dispersion coefficients can be accurately obtained in capillaries, where corresponding analytical solutions are available, even when these are resolved to just a few lattice units. Then we compute molecular displacement distributions for pore-spaces of varying complexity: a pack of beads; a Bentheimer sandstone; and a Portland carbonate. Our calculated propagator distributions are compared directly with recent experimental PFG-NMR propagator distributions (Scheven et al., 2005; Mitchell et al., 2008), the latter excluding spin relaxation mechanisms. We observe that the calculated transport propagators can be quantitatively compared with the experimental distribution, provided that spin relaxations in the experiment are excluded, and good agreement is found for both the sandstone and the carbonate. However, due to the absence of explicit micro-porosity from the carbonate pore space image used for flow field simulations we note that there are fundamental differences in the physical origins of the stagnant zones for micro-porous rocks between simulation and experiment. We show that for a given micro-CT image of a carbonate, small variations in the parameters chosen for the segmentation process lead to different amounts of stagnancy which diffuse away at different rates. Finally, we use a filtering method to show that this is due to the presence of spurious isolated pores which arise from the segmentation process and suggest an approach to overcome this limitation.  相似文献   

12.
多孔岩石波传播的热弛豫模型修正   总被引：2，自引：1，他引：1       下载免费PDF全文
经研究发现热弛豫模型的衰减比BISQ模型大得较多,与地球介质相比衰减量也过大;反演与实验结果相比在虚模量的低频(或低温)端和高频(或高温)端相差较大,仅在峰值附近符合较好.针对上述不足,将Arrhenius关系直接引进Biot模型,替换原模型引进的峰值点的频率对数和温度关系,并重新调整了模型参数.这样既改善了原模型衰减量过大,又克服了模型反演中虚模量峰值曲线两侧差异较大的缺点.进行了P波和S波的波传播分析,仍然在频率谱和温度谱上获得热弛豫峰和Biot峰.分析显示热弛豫峰导致波速随频率升高而上升的普遍规律,Biot峰导致波速随温度升高而上升的异常现象.在相同条件下对Biot模型,BISQ模型和热弛豫模型的P波波速和衰减进行了对比.热弛豫模型得到的速度频散更强,频散范围更宽,所得的衰减峰值频率比BISQ模型要低,衰减幅度比BISQ模型稍大.这些结果与实验结果相近,更符合实际.  相似文献   

13.
A multilayer lattice Boltzmann (LB) model is introduced to solve three-dimensional wind-driven shallow water flow problems. The multilayer LB model avoids the expensive Navier–Stokes equations and obtains stratified horizontal flow velocities as vertical velocities are relatively small and the flow is still within the shallow water regime. A single relaxation time BGK method is used to solve each layer coupled by the vertical viscosity forcing term. To increase solution stability, an implicit step is suggested to obtain flow velocities. The main advantage of using the LBM is that after selecting appropriate equilibrium distribution functions, the LB algorithm is only slightly modified for each layer and retains all the simplicities of the LBM within the high performance computing (HPC) environment. The performance of the parallel LB model for the multilayer shallow water equations is investigated on CPU-based HPC environments using OpenMP. We found that the explicit loop control with cache optimization in LBM gives better performance on execution time, speedup and efficiency than the implicit loop control as the number of processors increases. Numerical examples are presented to verify the multilayer LB model against analytical solutions. We demonstrate the model’s capability of calculating lateral and vertical distributions of velocities for wind-driven circulation over non-uniform bathymetry.  相似文献   

14.
It has been known for many years that dispersivities increase with solute displacement distance in a subsurface. The increase of dispersivities with solute travel distance results from significant variation in hydraulic properties of porous media and was identified in the literature as scale‐dependent dispersion. In this study, Laplace‐transformed analytical solutions to advection‐dispersion equations in cylindrical coordinates are derived for interpreting a divergent flow tracer test with a constant dispersivity and with a linear scale‐dependent dispersivity. Breakthrough curves obtained using the scale‐dependent dispersivity model are compared to breakthrough curves obtained from the constant dispersivity model to illustrate the salient features of scale‐dependent dispersion in a divergent flow tracer test. The analytical results reveal that the breakthrough curves at the specific location for the constant dispersivity model can produce the same shape as those from the scale‐dependent dispersivity model. This correspondence in curve shape between these two models occurs when the local dispersivity at an observation well in the scale‐dependent dispersivity model is 1·3 times greater than the constant dispersivity in the constant dispersivity model. To confirm this finding, a set of previously reported data is interpreted using both the scale‐dependent dispersivity model and the constant dispersivity model to distinguish the differences in scale dependence of estimated dispersivity from these two models. The analytical result reveals that previously reported dispersivity/distance ratios from the constant dispersivity model should be revised by multiplying these values by a factor of 1·3 for the scale‐dependent dispersion model if the dispersion process is more accurately characterized by scale‐dependent dispersion. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

15.
In this paper, frictional strength of hard solids, such as rock–rock sliding surfaces, is studied as a function of waiting time and shearing velocity. A one dimensional spring–mass sliding system is numerically simulated under the quasistatic condition using the rate and state dependent friction model. It is established that frictional strength varies linearly with the logarithm of waiting time (also known as time of stationary contact or relaxation time, etc.) as well as logarithm of shearing velocity. Analytical expression developed for frictional strength is found to be valid only in the case of high stiffness of the connecting spring. In the steady relaxation simulation, a steadily sliding mass is suddenly brought to zero velocity and relaxation of the interfacial stress and corresponding velocity at the sliding interface is studied as a function of relaxation time in the velocity strengthening regime of friction. A mathematical relation is derived between state variable and waiting time using the concept of steady relaxation. The relaxation model is also compared with the experimental data from the literature. Finally, the present study enables one to unify the slide–hold–slide friction experiments.  相似文献   

16.
This paper proposes a simple lattice model for collapse analysis of RC bridges subjected to earthquakes by using the extended distinct element method (EDEM). In the model, a concrete element consists of lumped masses connected to one another by springs, and a reinforcement bar is represented by a discrete model or an integrated model. The proposed lattice model is simple but its parameters are reasonably defined. It has fewer element nodes and connecting springs, which will be of benefit by shortening the CPU time. The processes to determine the initial stiffness of concrete and steel springs, the parameters of the constitutive model and the fracture criteria for springs are described. A re‐contact spring model is also proposed to simulate the re‐contact of the concrete after fracture of springs; and a general grid searching method is used to decrease the CPU time for judging re‐contact after fracture. The lattice model is assessed by numerical simulations and experiments. As an application, a damaged single‐column pier subjected to the Kobe Earthquake in 1995 is analysed by EDEM with the proposed model. The simulation results indicate that the proposed model predicts well qualitatively the collapse process of RC bridges. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

17.
Fluid flow exerts a critical impact on the convection of thermal energy in geological media, whereas heat transport in turn affects fluid properties, including fluid dynamic viscosity and density. The interplay of flow and heat transport also affects solute transport. To unravel these complex coupled flow, heat, and solute transport processes, here, we present a theory for the idealized scale‐dependent Poiseuille flow model considering a constant temperature gradient (?T) along a single fracture, where fluid dynamic viscosity connects with temperature via an exponential function. The idealized scale‐dependent model is validated based on the solutions from direct numerical simulations. We find that the hydraulic conductivity (K) of the Poiseuille flow either increases or decreases with scales depending on ?T > 0°C/m or ?T < 0°C/m, respectively. Indeed, the degree of changes in K depends on the magnitude of ?T and fracture length. The scale‐dependent model provides an alternative explanation for the well‐known scale‐dependent transport problem, for example, the dispersion coefficient increases with travel distance when ?T > 0°C/m according to the Taylor dispersion theory, because K (or equivalently flux through fractures) scales with fracture length. The proposed theory unravels intertwined interactions between flow and transport processes, which might shed light on understanding many practical geophysical problems, for example, geothermal energy exploration.  相似文献   

18.
A diffusive process is said to be anomalous if in any given direction the average square of the separation a particle experiences from its origin grows nonlinearly with time. Any diffusive process is anomalous if viewed on a short enough time scale, but interestingly, many diffusive processes remain anomalous over longer times. As a canonical example we study one such process here, diffusion in a laterally-confined nano-film as a function of the strain induced critical point. For this example we motivate and illustrate how a simple but novel process, Brownian motion run with a nonlinear clock (Bm-nlc), statistically mimics trajectories generated via Newton’s force law. The model is easily generalized to more complicated random processes and has application in many fields, including but not limited to, random conductivity field or terrain generation, Richardson turbulence in the atmosphere, and time dependent dispersion in hydrology.  相似文献   

19.
We investigate the early time induced polarization (IP) phenomenon in frozen unconsolidated rocks and its association with transient electromagnetic (TEM) signals measured in northern regions. The distinguishing feature of these signals is the distortion of the monotony or sign reversals in the time range from a few tens to a few hundreds of microseconds. In simulating TEM data, the IP effects in frozen ground were attributed to the dielectric relaxation phenomenon rather than to the frequency‐dependent conductivity. This enabled us to use laboratory experimental data available in the literature on dielectric spectroscopy of frozen rocks. In our studies we focused on simulating the transient response of a coincident‐loop configuration in three simple models: (i) a homogeneous frozen earth (half‐space); (ii) a two‐layered earth with the upper layer frozen; (iii) a two‐layered earth with the upper layer unfrozen. The conductivities of both frozen and unfrozen ground were assumed to exhibit no frequency dispersion, whereas the dielectric permittivity of frozen ground was assumed to be described by the Debye model. To simplify the presentation and the comparison analysis of the synthetic data, the TEM response of a frozen polarizable earth was normalized to that of a non‐polarizable earth having the same structure and resistivities as the polarizable earth. The effect of the dielectric relaxation on a TEM signal is marked by a clearly defined minimum. Its time coordinate tmin is approximately three times larger than the dielectric relaxation time constant τ. This suggests the use of tmin for direct estimation of τ, which, in turn, is closely associated with the temperature of frozen unconsolidated rock. The ordinate of the minimum is directly proportional to the static dielectric permittivity of frozen earth. Increasing the resistivity of a frozen earth and/or decreasing the loop size results in a progressively stronger effect of the dielectric relaxation on the TEM signal. In the case of unfrozen earth, seasonal freezing is not likely to have an appreciable effect on the TEM signal. However, for the frozen earth, seasonal thawing of a near‐surface layer may result in a noticeable attenuation of the TEM signal features associated with dielectric relaxation in a frozen half‐space. Forward calculations show that the dielectric relaxation of frozen unconsolidated rocks may significantly affect the transient response of a horizontal loop laid on the ground. This conclusion is in agreement with a practical example of inversion of the TEM data measured over the permafrost.  相似文献   

20.
I. Alhama  A. Soto Meca  F. Alhama 《水文研究》2012,26(24):3725-3735
The governing equations of coupled density‐driven flow and solute transport problems in porous media, with velocity‐dependent dispersion coefficient, are strongly nonlinear and must be solved numerically. This contribution presents a network model, based on the network simulation method, capable of simulating the transient solution to this kind of problem efficiently and with a relatively low computational time. The mathematical model is formulated using the stream function and concentration variables. Simulation of the network model is carried out in the standard electric circuit simulation code, Pspice. The present model is first applied to simulate the original benchmark Henry problem, and the solution is compared with those obtained by other authors. A study of the grid size is also carried out. In addition, the modified version of Simpson and Clement of the Henry problem, as well as the groundwater flow in the closed desert basin of Pilot Valley, is studied using the proposed model. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号