Saltwater intrusion modeling in heterogeneous confined aquifers using two-relaxation-time lattice Boltzmann method     

Borja Servan-Camas  Frank T.-C. Tsai 《Advances in water resources》2009

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.  相似文献   

A temperature-dependent single-crystal Raman spectroscopic study of fayalite: evidence for phonon-magnetic excitation coupling     

B. A. ?Kolesov  C. A. ?GeigerEmail author 《Physics and Chemistry of Minerals》2004,31(3):155-161

The polarized single-crystal Raman spectrum of synthetic fayalite, Fe₂SiO₄, was recorded between 5 and 773 K in order to investigate its lattice dynamic behavior. A broad absorption envelope is observed at wavenumbers between 800 and 960 cm^–1 and it contains two intense bands at 816 and 840 cm^–1 at 293 K in the (cc) spectrum. The integral area of the envelope decreases upon cooling from 293 K and reaches a minimum around 55 K. It then increases again with a further decrease in temperature down to 5 K. It is proposed that the envelope in the (cc) spectra consists of seven different modes, some of which are symmetry-forbidden, that arise from combination scattering of nonsymmetric internal SiO₄-stretching modes of B_ig symmetry (i = 1, 2, 3) and low-energy excitations. The individual modes can be observed under different polarizations and agree in number and wavenumber with those obtained by fitting the broad envelope with Lorentzians. An analysis of the Raman spectrum as a function of temperature, using the known magnetic properties of fayalite, allows the assignment of the low-energy excitations to short-range magnetic interactions. Modulation of the Fe²⁺(1)–Fe²⁺(2) exchange energy leads to phonon-magnetic excitation coupling and the main role in the Fe²⁺(1)–Fe²⁺(2) magnetic interaction occurs via superexchange through the oxygens. The magnetic excitations are not magnons in the usual sense, that is as quasiparticles having a long wavelength in an ordered system. The degree of observed broadening of the SiO₄-stretching modes is consonant with a Fe²⁺(1)–Fe²⁺(2) exchange energy of 4.7 cm^–1 presented by Schmidt et al. (1992). At temperatures above 300 K the line width of the mode at 840 cm^–1 decreases slightly, whereas those of low energy lattice modes increase. This suggests that a decrease in mode broadening due to weakened magnetic interactions compensates any thermally related broadening. Complete Fe²⁺ spin disorder may not be reached until at least 530 K. Results from this study show that estimates of third-law entropies for silicates using simple crystal-chemical considerations that do not account for magnetic properties cannot give accurate values for many transition-metal-containing phases.  相似文献   

Laboratory validation of lattice Boltzmann method for modeling pore-scale flow in granular materials     

Muhammed E. Kutay  Ahmet H. Aydilek  Eyad Masad 《Computers and Geotechnics》2006,33(8):381-395

Characteristics of fluid flow through various engineering structures, such as granular filters and asphalt pavements, influence their design life. Numerical simulation of fluid flow is useful for evaluating the hydraulic characteristics of these materials. Among various techniques, the lattice Boltzmann (LB) method is widely accepted due to the ease of implementing boundary conditions and the numerical stability in a wide variety of flow conditions. It has proven to be extremely efficient in the simulation of fluid flow through the complex geometries of granular materials. In this study, two-dimensional and three-dimensional LB models were developed to represent pore-scale monophasic Newtonian incompressible fluid flow in granular materials. Three-dimensional geometries of compacted aggregates and asphalt specimens were generated from X-ray Computed Tomography technique and used as input for the LB model. The accuracy of the models was verified by comparing the results with analytical solutions of simple geometries and hydraulic conductivity measurements on the compacted aggregates and hot mix asphalt specimens. The results of LB simulations were in excellent agreement with those obtained from analytical calculations and laboratory measurements.  相似文献   

Pore-scale simulation of entrapped non-aqueous phase liquid dissolution     

C. Pan  E. Dalla  D. Franzosi  C.T. Miller 《Advances in water resources》2007

We investigated the dissolution of non-aqueous phase liquids (NAPLs) in a three-dimensional random sphere-pack medium using a pore-scale modeling approach to advance fundamental understanding and connect rigorously to microscale processes. Residual NAPL distributions were generated using a morphological approach and the entrapped non-wetting phase was quantitatively characterized by calculating volume, orientation, interfacial area, and shape of isolated NAPL regions. With a detailed aqueous-phase flow field obtained by a multiple-relaxation time lattice Boltzmann approach, we solved the advective–diffusive equation in the pore space using a high-resolution, adaptive-stencil finite-volume scheme and an operator-splitting algorithm. We show good agreement between the mass transfer rates predicted in the computational approach and previously published experimental observations. The pore-scale simulations presented in this work provide the first three-dimensional comparison to the considerable experimental work that has been performed to derive constitutive relations to quantify mass transfer from a residual NAPL to a flowing aqueous phase.  相似文献   

湘北杨家坪中新元古宇和下古生界的近变质作用与成岩作用   总被引：1，自引：1，他引：0  

王河锦  周钊  王玲  苑蕾 《岩石学报》2014,30(10):3013-3020

运用伊利石结晶度、K云母晶格常数b0、应变值τ、颗粒大小分布与矿物组合分析等低温变质指示数据对湘北杨家坪剖面中晚元古宇和下古生界的近变质作用和成岩作用进行了研究。Kisch国际标样用于校正伊利石结晶度的测定。中晚元古宇伊利石结晶度Kübler指数范围为:0.21~0.24°Δ2θ,而下古生界为0.28~0.67°Δ2θ。表明近变质作用影响了中晚元古宇和部分下古生界。按照Kisch国际标样的成岩/变质界限(0.21~0.38°Δ2θ,CuKα),将剖面划分为北部近变质带和中南部成岩带。由绿泥石化学成分地质温度计估计近变质峰期温度约为260℃。K云母晶格常数b0值变化范围0.9000~0.9045nm,平均为0.9017nm,表明近变质作用处于中压-中低压力范围。近变质带伊利石(主要是中晚元古宇)的多型为2M1型,而成岩带伊利石(主要是下古生界)的多型为2M1+1M混合类型。平均对数粒度与对数方差之积(GLR)称为对数总成熟效应具有不同的两组分布特征即近变质作用区(GLR1.88)与成岩作用区(GLR1.88),这是一个新的发现。后变质作用的构造应力不仅可产生断层,而且可引起矿物的应变,从而导致衍射峰的宽化,使得近变质带伊利石结晶度值次生变大到成岩值范畴。与剖面南部地区(长沙-澧陵-浏阳,黄土店-仙溪,沅古坪)比较,杨家坪剖面中晚元古宇和下古生界经历了中压-中低压近变质作用和成岩作用的影响。这一结果与多数前人认为的"元古宇板溪群处于绿片岩相或亚绿片岩相或低绿片岩相而震旦到下古生界为沉积盖层"的观点不同。  相似文献   

格构梁设计计算系统的开发研究     

曾勇辉  刘亚 《岩土工程技术》2010,24(1):12-15,43

格构梁作为一种新型的边坡防护体系,目前在工程实践中已广泛应用。然而工程设计中往往将格构梁拆分为单片梁计算,忽略了格构梁作为整体框架的特点,且设计主要采用手工计算和手工绘图,设计质量和设计速度都难以令人满意。基于以上考虑,利用十字交叉基础设计方法,对纵横梁的荷载分配情况进行讨论,采用Winkler弹性地基梁理论对格构梁的内力进行分析。在此基础上,将格构梁内力计算、配筋计算及施工图绘制集于一体,构建格构梁设计计算系统流程图,利用可视化程序设计语言Visual Basic结合CAD二次开发技术,开发相应的设计计算程序,实现了格构梁设计的智能化和一体化。  相似文献   

龙马溪组页岩数字岩心LSM-RVM数值建模方法研究及TOC含量影响分析          下载免费PDF全文

刘宁  符力耘  曹呈浩  刘建林 《地球物理学报》2020,63(7):2774-2785

页岩气储层中含有大量有机碳(TOC),其丰度与成熟度对页岩力学特性有重要影响.建立包含TOC的精细数值模型,将有助于探索页岩微结构与矿物组分含量对等效弹性模量的作用程度,是“甜点区”预测的重要理论基础.本文提出了一种离散数值建模方法,基于高精度成像技术,采用晶格弹簧-随机孔隙耦合模型(LSM-RVM)模拟包含多种矿物组分及不同成熟度干酪根的数字岩心,分析TOC成熟度及含量对弹性参数的影响.在该模型中,参数设置(数值阻尼与加载应变速率)至关重要,选取不当会对计算精度造成一定影响.研究结果表明,LSM-RVM能够生成符合TOC及多种矿物实际分布特征的数值模型,是一种精细数值建模方法.  相似文献   

地震波各向异性的研究进展   总被引：20，自引：10，他引：10  

黄晓葛  白武明 《地球物理学进展》1999,14(3):54-65

地壳和上地幔矿物的优势取向的改变是产生地震波速各向异性的主要原因。本文重点概述了在高温高压实验方面地震波各向异性的研究进展,总结了它的研究成果,并讨论了产生地震波各向异性的主要原因,阐述了地震波各向异性在地球动力学中的意义。  相似文献   

浅水动力学的格子Boltzmann模拟   总被引：6，自引：0，他引：6  

刘峰  胡非 《海洋与湖沼》2004,35(1):35-40

建立了一个格子Boltzmann模型,经过Chapman-Enskog展开可以恢复出守恒形式的浅水动力学方程组。理论分析和数值试验的结果都表明,本模型严格地满足质量守恒定律和动量定理。数值模拟了一维和二维溃坝流动,分别与理论解和文献值进行了对比,结果相当吻合。本模型具有很好的并行计算性能,在海洋、大气等大规模数值计算领域有着广阔的发展前景。  相似文献   

Numerical analysis of water exit for a sphere with constant velocity using the lattice Boltzmann method     

《Applied Ocean Research》2019

In this paper, the three-dimensional water exit of a sphere with different vertical velocities is investigated numerically using the lattice Boltzmann method (LBM). In this method, the liquid-gas two-phase flow is simplified as a single-phase free surface flow. To capture the free surface, a mass tracking algorithm is incorporated into the LBM. The gravity as a body force is introduced in the form of calculating the equilibrium distribution with an altered velocity, while the surface tension is neglected. Besides, the employed bounce-back boundary conditions are used for a moving sphere. What’s more, the Wall-Adapting Local Eddy (WALE) viscosity model is employed to capture the turbulent structures of the flow and stabilize the simulation. The accuracy of the numerical results is demonstrated through comparisons with the previous numerical and experimental results in the literature. The results show that the spike height is significantly influenced under the Froude number (Fr) below 4.12 and slightly affected under the Fr varying from 4.12 to 8.24. After the sphere exits water totally, the evolution of the free surface waterfall can be described as two phases and becomes more intense with the Froude number increasing. The non-uniform distribution of velocity results in the breaking of the free surface after the sphere completely exits the water. Moreover, the Reynolds number greatly affects the wake dynamics and hydrodynamics acting on the sphere when it moves beneath the water surface.  相似文献