柱几何Rayleigh-Taylor不稳定性的数值模拟

给出了柱几何中流体力学方程组及其在数值模拟中采用的计算方法。对二维柱几何Rayleigh-Taylor不稳定性进行数值模拟,在线性阶段与线性理论符合得很好;不稳定性增长进入非线性区域的阈值依赖于界面的位置,并且明显不同于平面情况。     

An anti‐diffusive volume of fluid method for interfacial fluid flows

A volume of fluid (VOF) method is developed combining a first‐order limited downwind scheme with higher order accurate schemes. The method is characterized by retaining a sharp fluid interface and a reduction in numerical diffusion near the interface, but avoids complicated geometrical reconstruction as occurs in most volume tracing algorithms. To demonstrate the accuracy and robustness of the method, a selection of numerical experiments are presented involving a pure advection problem, a water wave impact caused by a dam breaking and liquid sloshing in a partially filled tank. Copyright © 2011 John Wiley & Sons, Ltd.     

A simplified adaptive Cartesian grid system for solving the 2D shallow water equations

This paper presents a new simplified grid system that provides local refinement and dynamic adaptation for solving the 2D shallow water equations (SWEs). Local refinement is realized by simply specifying different subdivision levels to the cells on a background uniform coarse grid that covers the computational domain. On such a non‐uniform grid, the structured property of a regular Cartesian mesh is maintained and neighbor information is determined by simple algebraic relationships, i.e. data structure becomes unnecessary. Dynamic grid adaptation is achieved by changing the subdivision level of a background cell. Therefore, grid generation and adaptation is greatly simplified and straightforward to implement. The new adaptive grid‐based SWE solver is tested by applying it to simulate three idealized test cases and promising results are obtained. The new grid system offers a simplified alternative to the existing approaches for providing adaptive mesh refinement in computational fluid dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

CH4/空气简化动力学机理数值验证

选择绕圆柱预混燃烧算例，验证CH₄/空气三种简化动力学机理（16s41r、15s19r和53s325r）.考虑均匀来流，忽略湍流和湍流与燃烧相互作用以及燃料扩散效应，假设层流有限反应速率，采用保自由流5阶WENO格式求解多组分Euler方程组，得到CH₄/空气预混燃烧流场温度等值线、沿驻点线压力和温度及其CH₄、CO和CO₂质量百分数分布.结果表明：三种简化动力学机理给出的流场均出现弓形激波和火焰面，弓形激波和火焰驻点距离及其形状、诱导区宽度和简化动力学机理相关.当圆柱直径增大，弓形激波和火焰向圆柱上游移动，对应的驻点距离均增大，诱导区宽度变短，点火延时变小，但火焰和弓形激波位置次序未变化.53s325r模型要比16s41r模型和15s19r模型精度要高，点火延时覆盖的压力和温度范围也较宽，所有简化机理均未完全反应，在较大圆柱直径下游达到化学平衡.     

A three‐dimensional Cartesian cut cell method for incompressible viscous flow with irregular domains

A three‐dimensional Cartesion cut cell method is presented for the simulations of incompressible viscous flows with irregular domains. A new model (referred to as ‘6+N’ model) is proposed to describe arbitrarily shaped cut cells and treat all the cells as polyhedrons with 6+N faces. The finite volume discretization of the Navier–Stokes equation is then implemented by using the ‘6+N’ model to separate the surface flux integrals into two parts, that is, the fluxes through the basic face of the hexahedron and those through the cutting surfaces. The previously proposed Kitta Cube algorithm and volume computer‐aided design platform (J. Comput. Aided. Des. 2005; 37(4): 1509–1520. Doi:10.1016/j.cad.2005.03.006) are adopted to generate cut cells and provide shape data and physical attributes for the numerical analysis. A modified SIMPLE‐based smoothing pressure correction scheme is applied to suppress checkerboard pressure oscillations caused by the collocated arrangement of velocities and pressure. The calculation accuracy of the numerical method expressed by L₁ and L _∞ norm errors is first demonstrated by the simulation of a pipe flow. Then its feasibility, efficiency, and potential in engineering applications are verified by applying it to solve natural convections between concentric spheres and between eccentric spheres. The heat transfer patterns in eccentric spheres are also obtained by using the numerical method. Copyright © 2011 John Wiley & Sons, Ltd.     

Implementation of ADER scheme for a bore on an unsaturated permeable slope

The paper details the implementation of the Godunov‐type finite volume Arbitrary high order schemes using Derivatives (ADER) scheme for the case of a large source term in the continuity equation of the nonlinear shallow water equations. The particular application is the movement of a bore on a highly permeable slope. The large source term is caused by the infiltration into the initially unsaturated slope material. Infiltration is modelled as vertical downwards piston‐like flow with Forchheimer quadratic parameterisation of the resistance law. The corresponding ODE is solved using the fourth‐order Runge–Kutta method. The surface and subsurface flow models have been tested by comparison with analytical solutions. Example predictions of surface bore propagation and wetting front propagation are presented for a range of slope permeabilities. The effects of permeability on bore run‐up, water depths and velocities are illustrated. The ADER scheme is capable of handling the source term, including the extreme case when this term dominates the volume balance. Copyright © 2011 John Wiley & Sons, Ltd.     

Investigation of the solubilization of car-emitted Pt,Pd and Rh in street dust and spiked soil samples

A sequential extraction procedure (three-step), proposed by the Standards, Measurements and Testing Programme (formerly BCR) of the European Union, was applied to street dust and spiked soil samples for the determination of PGEs. Analyses were carried out using inductively coupled plasma-mass spectrometry (ICP-MS). The results indicate that up to 5% from Pt, 70% from Pd and 14% from Rh are in mobile forms in street dust. The results for the soil samples spiked with crushed catalytic converter are significantly lower indicating that PGEs are oxidised more efficiently in natural conditions. Additionally Pt and Pd bound to humic acids were investigated.     

Longitudinal Relaxation Enhancement in 1H NMR Spectroscopy of Tissue Metabolites via Spectrally Selective Excitation

Nuclear magnetic resonance spectroscopy is governed by longitudinal (T₁) relaxation. For protein and nucleic acid experiments in solutions, it is well established that apparent T₁ values can be enhanced by selective excitation of targeted resonances. The present study explores such longitudinal relaxation enhancement (LRE) effects for molecules residing in biological tissues. The longitudinal relaxation recovery of tissue resonances positioned both down‐ and upfield of the water peak were measured by spectrally selective excitation/refocusing pulses, and compared with conventional water‐suppressed, broadband‐excited counterparts at 9.4 T. Marked LRE effects with up to threefold reductions in apparent T₁ values were observed as expected for resonances in the 6–9 ppm region; remarkably, statistically significant LRE effects were also found for several non‐exchanging metabolite resonances in the 1–4 ppm region, encompassing 30–50 % decreases in apparent T₁ values. These LRE effects suggest a novel means of increasing the sensitivity of tissue‐oriented experiments, and open new vistas to investigate the nature of interactions among metabolites, water and macromolecules at a molecular level.     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.