Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.     

A model for moisture and heat flow in fractured unsaturated porous media

The present study investigates propagation of a cohesive crack in non‐isothermal unsaturated porous medium under mode I conditions. Basic points of skeleton deformation, moisture, and heat transfer for unsaturated porous medium are presented. Boundary conditions on the crack surface that consist of mechanical interaction of the crack and the porous medium, water, and heat flows through the crack are taken into consideration. For spatial discretization, the extended finite element method is used. This method uses enriched shape functions in addition to ordinary shape functions for approximation of displacement, pressure, and temperature fields. The Heaviside step function and the distance function are exploited as enrichment functions for representing the crack surfaces displacement and the discontinuous vertical gradients of the pressure and temperature fields along the crack, respectively. For temporal discretization, backward finite difference scheme is applied. Problems solved from the literature show the validity of the model as well as the dependency of structural response on the material properties and loading. Copyright © 2016 John Wiley & Sons, Ltd.     

The effective stress concept and the evaluation of changes in pore air pressure under jacketed isotropic compression tests for dry sand based on the 2‐phase porous theory

The effective stress concept for solid‐fluid 2‐phase media was revisited in this work. In particular, the effects of the compressibility of both the pore fluid and the soil particles were studied under 3 different conditions, i.e., undrained, drained, and unjacketed conditions based on a Biot‐type theory for 2‐phase porous media. It was confirmed that Terzaghi effective stress holds at the moment when soil grains are assumed to be incompressible and when the compressibility of the pore fluid is small enough compared to that of the soil skeleton. Then, isotropic compression tests for dry sand under undrained conditions were conducted within the triaxial apparatus in which the changes in the pore air pressure could be measured. The ratio of the increment in the cell pressure to the increment in the pore air pressure, m, corresponds to the inverse of the B value by Bishop and was obtained during the step loading of the cell pressure. In addition, the m values were evaluated by comparing them with theoretically obtained values based on the solid‐fluid 2‐phase mixture theory. The experimental m values were close to the theoretical values, as they were in the range of approximately 40 to 185, depending on the cell pressure. Finally, it was found that the soil material with a highly compressible pore fluid, such as air, must be analyzed with the multi‐phase porous mixture theory. However, Terzaghi effective stress is practically applicable when the compressibilities of both the soil particles and the pore fluid are small enough compared to that of the soil skeleton.     

Why asymptotic limit of the acid dissolution capacity can lead to a sharp dissolution front in chemical dissolution of porous rocks?

The use of the asymptotic limit can greatly simplify the theoretical analysis of chemical dissolution front instabilities in fluid‐saturated rocks and therefore make it possible to obtain mathematical solutions, which often play a crucial role in understanding the propagation behavior of chemical dissolution fronts in chemical dissolution systems. However, there has been a debate in recent years that the asymptotic limit of the acid dissolution capacity (i.e., the acid dissolution capacity number approaching zero) alone cannot lead to a sharp dissolution front of the Stefan type in the acidization dissolution system, in which the dissolvable minerals of carbonate rocks are chemically dissolved by the injected acid flow. The acid dissolution capacity number is commonly defined as the ratio of the volume of the carbonate rock dissolved by an acid to that of the acid. In this paper, we use four different proof methods, including (i) direct use of the fundamental concepts; (ii) use of the mathematical governing equations of an acidization dissolution system; (iii) use of the different time scaling approach; and (iv) use of a moving coordinate system approach, to demonstrate that the asymptotic limit of the acid dissolution capacity can indeed lead to sharp dissolution fronts of the Stefan type in acidization dissolution systems on a much larger time scale (than the dissolution time scale). Our new finding is that on the reaction time scale, the condition of the conventional time derivative of porosity approaching zero alone can ensure that the acidization dissolution front has a sharp shape of the Stefan type. Copyright © 2017 John Wiley & Sons, Ltd.     

Analytical solutions for advective–dispersive solute transport in double‐layered finite porous media

Analytical solutions for advection and dispersion of a conservative solute in a one‐dimensional double‐layered finite porous media are presented. The solutions are applicable to five scenarios that have various combinations of fixed concentration, fixed flux and zero concentration gradient conditions at the inlet and outlet boundaries that provide a wide number of options. Arbitrary initial solute concentration distributions throughout the media can be considered via explicit formulations or numerical integration. The analytical solutions presented have been verified against numerical solutions from a finite‐element‐based approach and an existing closed‐form solution for double‐layered media with an excellent correlation being found in both cases. A practical application pertaining to advective transport induced by consolidation of underlying sediment layers on contaminant movement within a capped contaminated sediment system is presented. Comparison of the calculated concentrations and fluxes with alternative approaches clearly illustrates the need to consider advection processes. Consideration of the different features of contaminant transport due to varying pore‐water velocity fields in primary consolidation and secondary consolidation stages is achieved via the use of non‐uniform initial concentration distributions within the proposed analytical solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

BDS中长距离非差实时动态厘米级定位算法

针对BDS常规实时动态定位(RTK)中,随着流动站与参考站间的距离增加,大气延迟误差的空间相关性大大降低,影响了整周模糊度的快速解算和流动站位置信息的解算精度问题。该文研究了一种基于非差观测误差的BDS中长距离常规RTK定位算法,采用非差误差改正方法为流动站提供误差改正,利用参考站的非差误差改正数以单颗卫星为对象进行误差改正。对电离层延迟误差和相对天顶对流层延迟误差进行参数估计,处理电离层延迟误差和对流层延迟误差的影响。最后通过BDS实测数据对该算法进行了算法验证和结果分析。实验结果表明,该算法可以实现BDS中长距离常规RTK的快速定位,并获得厘米级定位精度。     

强壮前沟藻(Amphidinium carterae Hulbert)在多重逆境条件下的生长和生理响应

尽管单一逆境条件对浮游植物产生的影响已被广泛重视,然而对于多重逆境条件对藻细胞生长生理的研究仍有空白。本研究通过测定强壮前沟藻(Amphidinium carterae Hulbert)光密度(OD750)、叶绿素a含量(chl a)、光合作用效率(F_v/F_m)、碳氮比(C/N)等多个生理指标,分析其在9种环境下(常温光照、低温和低温黑暗三种物理环境,结合全营养、缺氮、缺磷三种营养状态)对多重逆境的生长和生理响应。研究结果表明,在多重逆境条件下,单一物理逆境因素(低温或黑暗)较氮限制或磷限制对藻细胞生长和生理的不利影响更为显著。营养限制并同低温环境双重作用对生物量和碳氮比产生显著性影响(P0.01)。此外,低温和黑暗条件耦合作用下,SYTOX Green染色强度处于较低水平,chl a稳定、F_v/F_m有所升高,强壮前沟藻在逆境环境下作为群体的衰亡得以缓解。     

Interaction of cnoidal waves with an array of vertical concentric porous cylinders

An array of large concentric porous cylinder arrays is mounted in shallow water exposed to cnoidal waves. The interactions between waves and cylinders are studied theoretically using an eigenfunction expansion approach. Semi-analytical solutions of hydrodynamic loads and wave run-up on each cylinder are obtained using first approximation to cnoidal waves. The square array configuration of four-legged identical concentric porous cylinder is investigated in present study. Numerical results reveal the variation of dimensionless wave force and wave run-up on individual cylinder with angle of incidence, porosity parameter, spacing between outer and inner cylinders, spacing between concentric porous cylinders and wave parameter. Different mechanism of wave force is found under different range of scattering parameter.     

基于五种油源的中重度风化溢油的常规鉴定指标的应用效果评估

本文利用采集于我国三大油田的五种原油样品开展了长达210天的溢油风化模拟实验,并依据相对偏差和重复性限数学分析法,进行溢油风化过程分析和诊断比值应用效果评估。研究结果表明:经过210天的风化,溢油鉴定诊断比值发生明显改变;其中来源于萜烷、甾烷和多环芳烃的诊断比值变化率要远低于正构烷烃,可用于中长期风化溢油鉴定。此外,研究发现,在这些有效诊断比值中有4个变化率较小,相对偏差低于5%,保持了较好的稳定性,更适合于重度风化溢油鉴定。     

TTI介质拟声波方程数值模拟

TTI介质qP波数值模拟方法因为考虑了倾角因素,可以比VTI介质qP波数值模拟方法更加准确地描述各向异性介质中地震波场的传播规律。文中用拟声波方程对TTI介质中的地震波场进行了高阶有限差分数值模拟,在改进衰减函数分布方式后,通过坐标变换,利用改进的完全匹配层(perfectly matched layer,PML)边界控制方程对波场边界进行吸收处理,取得了良好的效果;然后分析了拟声波方程数值模拟中的稳定性问题,并对波场中的伪横波进行压制。通过对不同模型的数值模拟,验证了文中使用的TTI介质拟声波波动方程的稳定性以及所采用的PML边界控制方程的可靠性和适用性。