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

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

Contingent derivatives and regularization for noncoercive inverse problems

Abstract

We study the inverse problem of parameter identification in noncoercive variational problems that commonly appear in applied models. We examine the differentiability of the set-valued parameter-to-solution map using the first-order and the second-order contingent derivatives. We explore the inverse problem using the output least-squares and the modified output least-squares objectives. By regularizing the noncoercive variational problem, we obtain a single-valued regularized parameter-to-solution map and investigate its smoothness and boundedness. We also consider optimization problems using the output least-squares and the modified output least-squares objectives for the regularized variational problem. We give a complete convergence analysis showing that for the output least-squares and the modified output least-squares, the regularized minimization problems approximate the original optimization problems suitably. We also provide the first-order and the second-order adjoint method for the computation of the first-order and the second-order derivatives of the output least-squares objective. We provide discrete formulas for the gradient and the Hessian calculation and present numerical results.     

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.     

Approximant crystals of icosahedral Mg–(Ga,Al)–Zn alloys

The zero field cooled (ZFC) and field cooled (FC) low-field magnetic moment m of a dense frozen ferrofluid containing Fe₅₅Co₄₅ particles of size 4.6nm in hexane exhibits irreversibility at temperatures T?T _b≈ 30?K. FC in μ ₀ H ≤ 1?T gives rise to shifted minor hysteresis loops below T _b. At T _c≈ 10?K, sharp peaks of m ^ZFC and of the ac susceptibility χ ′, a kink of the thermoremanent magnetic moment m ^TRM, a sizeable reduction of the coercive field H _c, and the appearance of a spontaneous moment m ^SFM indicate a phase transition with near mean-field critical behaviour of both m ^SFM and χ ′ . These features are explained within a core-shell model of nanoparticles, whose strongly disordered shells gradually become blocked below T _b, while their soft ferromagnetic cores couple dipolarly and become superferromagnetic (SFM) below T _c.     

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.     

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

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

三维多面体网格上扩散方程的保正格式

 针对三维任意(星形)多面体网格, 本文构造了扩散方程的一种单元中心型非线性有限体积格式, 证明了该格式具有保正性. 在该格式设计中, 除引入网格中心量外, 还引入网格节点量和网格面中心量作为中间未知量, 它们将用网格中心未知量线性组合表示, 使得格式仅有网格中心未知量作为基本未知量. 在节点量计算中, 利用网格面上的调和平均点, 设计了一种适用于三维多面体网格的局部显式加权方法. 该格式适用于求解非平面的网格表面和间断扩散系数的问题. 数值例子验证了它对光滑解具有二阶精度和保正性.