保持球对称的相容流体算法的构造

重点研究如何在二维柱坐标系中保持流体运动的球对称问题.该问题需要对整个流体算法做细致的分析与设计,从而获得不同的方法.基于Caramana介绍的修正压强梯度算子构造了保持球对称的相容流体算法,该算法还包括了边人工粘性方法、子网格压力的计算方法,以及预估步与校正步相结合的策略.并且,考虑了总能量的守恒性,讨论了网格节点质量保持不变的必要性,以及基于相容格式证明了边人工粘性的耗散性.数值模拟结果与已知数据的比对验证了程序的正确性、稳定性.     

基于子网格边界近似Riemann解的人为粘性方法

在交错网格型Lagrange(拉格朗日)流体力学算法中,通常采用人为粘性捕捉激波,人为粘性的好坏对于计算结果至关重要.研究了一种基于子网格边界处近似Riemann解的新型人为粘性.新人为粘性能够满足动量守恒和熵不等式.利用子网格边界速度差中引入的限制器,新人为粘性能够区分激波和等熵压缩,并能满足球对称问题中的波面不变性.新人为粘性在典型模型数值模拟及惯性约束聚变黑腔整体数值模拟中取得了较好的结果.     

二维柱坐标系辐射输运方程保球对称的离散纵标格式

 离散纵标格式是计算辐射输运方程的常用格式之一. 但是, 传统的离散纵标格式求解二维柱坐标系辐射输运方程模拟一维球对称问题时, 会出现明显的非对称现象, 球对称性被破坏. 针对该问题, 本文分析了传统离散纵标格式不能够保持球对称性的原因, 提出了空间基于柱坐标系、方向基于球坐标系的辐射输运方程, 并对该方程设计了新的离散纵标格式, 从理论上证明了当空间网格取球对称剖分时该离散格式能够保持一维球对称性的充分必要条件. 通过对真空球区域辐射输运、与物质耦合辐射输运等球对称算例的数值模拟, 验证了该格式的保球对称性, 球对称误差能够达到机器精度. 非对称辐射驱动模型以及非对称网格剖分条件下的数值模拟等算例也取得了较好的结果.     

一种自由界面追踪的模板化VOF方法

发展了一种模板化的volume-of-fluid (VOF)方法．该方法根据自由界面的法向建立一个模板,然后由已知的网格单元上的流体体积比值确定出自由界面的准确位置,使得在二维情形下一个网格单元被自由界面切割的形式只有3种．另一方面,引入了单元边流体占有长度的概念,在此基础上建立了一个统一的流体占有面积模型,可以使得自由界面输运方程的求解有统一的算法．该方法不受网格单元形式的限制,并且容易推广到三维情形．算例表明,该方法能保证自由界面的跟踪精度．     

粘性机制下一类熵相容格式的对比研究

研究了一种人工和物理耗散机制下的离散熵相容格式,探讨数值粘性和物理粘性的大小以及它们所起的作用.所得结论是:在激波捕捉的过程中,粘性系数越大,则无需加入人工粘性项;粘性系数较小时,除了物理粘性项,还需要加入人工粘性项来得到熵相容格式.首先研究了一维粘性Burgers方程离散熵相容格式,再将其推广至Navier-Stokes方程.数值算例采用空间半离散格式,并结合显式三步三阶Runge-Kutta(RK3)方法进行时间推进.这两类方程的数值结果表明,最终选取的熵相容格式能够准确地捕捉到激波.     

Lagrange柱坐标高阶中心型守恒格式

提出Lagrange柱坐标高阶中心型守恒格式.基于用对守恒律的单调迎风算法(MUSCL)构造的高阶子网格压力,引入了柱坐标高阶体权子网格力和柱坐标高阶面权子网格力,构造了柱坐标高阶体权中心型守恒格式和柱坐标高阶面权中心型格式.柱坐标高阶体权中心型守恒格式满足动量守恒、能量守恒,但不能确定保持一维球对称性.柱坐标高阶面权中心型格式满足能量守恒,保持一维球对称性.两种格式里,格点速度以与网格面的数值通量相容的方式计算.对Saltzman活塞问题等进行了数值模拟,数值结果显示Lagrange柱坐标高阶中心型守恒格式的有效性和精确性.     

高可扩展区域分解算法及其在流体模拟和优化问题中的应用

本文介绍一套求解复杂流体模拟和优化控制问题的高可扩展并行算法.该算法基于非结构化网格,结合了加稳定化项的有限元空间离散方法、全隐的时间离散格式、多物理场全耦合的求解算法、区域分解算法及求解非线性系统的Newton-Krylov-Schwarz算法等多套先进算法.利用该算法,本文对多个实际工程应用中流体模拟和优化设计问题进行了测试,数值结果显示,该算法对本文研究的几类问题,具有很好的收敛性和并行可扩展性,当使用8192个处理器核求解规模超过两千万个网格单元的问题时,仍然具有超过40%的并行效率.     

非对称稀疏线性方程组的快速外存解法及其在无网格法计算中的应用

针对局部Petrov-Galerkin无网格法(MLPG)等无网格方法的计算所产生的大型非对称稀疏线性方程组,介绍了一种新的直接解法．与一般非对称求解过程不同,该解法从现有的对称正定解法中演变出来,其分解过程在矩阵的上、下三角阵中对称进行．新的矩阵分解算法可以通过修改对称矩阵分解算法的代码来实现,这提供了从对称解法到非对称解法的快捷转换．还针对MLGP法以及有限元法所产生的方程组开发了多块外存算法(multi-blocked out-of-core strategy)来扩大求解规模．测试结果证明该方法大幅度提高了大型非对称稀疏线性方程组的求解速度．     

求解Oseen流的交替线松弛多重网格方法

利用Riemann解的通量差分分裂法——Godunov方法对Oseen流控制方程进行离散,得到了基于一阶上迎风格式的离散方程,并给出了使用多重网格方法求解该离散方程的V-循环算法和W-循环算法的收敛性分析.通过局部Fourier分析方法,对获得的离散方程的聚对称交替线GaussSeidel松弛的光滑性质进行了研究.结果表明:使用多重网格的两层网格及三层网格算法求解具有不同Reynolds数的Oseen流,即便是在高Reynolds数情况下,聚对称交替线Gauss-Seidel松弛具有很好的光滑性质,多重网格W-循环算法收敛性比V-循环算法好.     

求解浅水波方程的熵相容格式

提出了一种求解浅水波方程组的熵相容格式．在熵稳定通量中添加特征速度差分绝对值的项来抵消解在跨过激波时所产生的熵增,从而实现熵相容．新的数值差分格式具有形式简单、计算效率高、无需添加任何的人工数值粘性的特点．数值算例充分说明了其显著的优点．利用新格式成功地模拟了不同类型溃坝问题的激波、稀疏波传播及溃坝两侧旋涡的形成,是求解浅水波方程组较为理想的方法.     

Global Existence of the Equilibrium Diffusion Model in Radiative Hydrodynamics

This paper is devoted to the analysis of the Cauchy problem for a system of PDEs arising in radiative hydrodynamics. This system, which comes from the so-called equilibrium diffusion regime, is a variant of the usual Euler equations, where the energy and pressure functionals are modified to take into account the effect of radiation and the energy balance containing a nonlinear diffusion term acting on the temperature. The problem is studied in the multi-dimensional framework. The authors identify the existence of a strictly convex entropy and a stability property of the system, and check that the Kawashima-Shizuta condition holds. Then, based on these structure properties, the well-posedness close to a constant state can be proved by using fine energy estimates. The asymptotic decay of the solutions are also investigated.     

The initial value problem for a multi‐dimensional radiation hydrodynamics model with viscosity

In this paper, we study the existence and time‐asymptotic behavior of solutions to the Cauchy problem for the equations of radiation hydrodynamics with viscosity in ?³. The global existence of the solutions is obtained by using the energy method. With more elaborate energy estimates, we also give some decay rates of the solutions. Copyright © 2010 John Wiley & Sons, Ltd.     

Decomposition algorithms with parametric Gomory cuts for two-stage stochastic integer programs

We consider a class of two-stage stochastic integer programs with binary variables in the first stage and general integer variables in the second stage. We develop decomposition algorithms akin to the $L$ -shaped or Benders’ methods by utilizing Gomory cuts to obtain iteratively tighter approximations of the second-stage integer programs. We show that the proposed methodology is flexible in that it allows several modes of implementation, all of which lead to finitely convergent algorithms. We illustrate our algorithms using examples from the literature. We report computational results using the stochastic server location problem instances which suggest that our decomposition-based approach scales better with increases in the number of scenarios than a state-of-the art solver which was used to solve the deterministic equivalent formulation.     

A computer-assisted proof of the existence of traveling wave solutions to the scalar Euler equations with artificial viscosity

We establish the existence and local uniqueness of traveling wave solutions to the one-dimensional Euler equations with artificial viscosity. The equations are expressed as a fixed-point problem, which is solved by a computer-assisted method based on Yamamoto’s application of the Banach fixed-point theorem.     

求解中大规模复杂凸二次整数规划问题的新型分枝定界算法

针对现有分枝定界算法在求解高维复杂二次整数规划问题时所存在的诸多不足，本文通过充分挖掘二次整数规划问题的结构特性来设计选择分枝变量与分枝方向的新方法，并将HNF算法与原问题松弛问题的求解相结合来寻求较好的初始整数可行解，由此导出可用于有效求解中大规模复杂二次整数规划问题的改进型分枝定界算法．数值试验结果表明所给算法大大改进了已有相关的分枝定界算法，并具有较好的稳定性与广泛的适用性．     

Asymptotic stability of viscous contact wave to a radiation hydrodynamic limit model

This paper is concerned with the large time behavior of the solutions for 1D radiation hydrodynamic limit model without viscosity and its asymptotic stability of the viscous contact discontinuity wave under the smallness assumption of the strength of the contact wave and initial perturbations. The present pressure includes a fourth-order term about the absolute temperature from radiation effect which brings the main difficulty. Furthermore, the dissipative of the system is weaker for the lack of viscosity. All these make the problem more challenging. The prove is mainly based on the energy method, including normal and radial directions energy estimates.     

Well-Posedness of a Problem of Propagation of Nonlinear Acoustic-Gravity Waves in the Atmosphere Generated by Surface Pressure Variations

Currently there are many international microbarograph networks for high-resolution recording of wave pressure variations on the Earth’s surface. This arouses interest in wave propagation in the atmosphere generated by atmospheric pressure variations. A full system of nonlinear hydrodynamic equations for atmospheric gases with lower boundary conditions in the form of wavelike pressure variations on the Earth’s surface is considered. Since the wave amplitudes near the Earth’s surface are small, linearized equations are used in the analysis of well-posedness of the problem. With the help of a wave energy functional method, it is shown that in the non-dissipative case the solution to the boundary value problem is uniquely determined by the variable pressure field on the Earth’s surface. The corresponding dissipative problem is well-posed if, in addition to the pressure field, appropriate conditions on the velocity and temperature on the Earth’s surface are given. In the case of an isothermal atmosphere, the problem admits analytical solutions that are harmonic in the variables x and t. A good agreement between the numerical and analytical solutions is obtained. The study shows that the temperature and density can rapidly vary at the lower boundary of the boundary value problem. An example of solving the three-dimensional problem with variable pressure on the Earth’s surface taken from experimental observations is given. The developed algorithms and computer programs can be used to simulate atmospheric waves generated by pressure variations on the Earth’s surface.     

On a modification of Marquardt's compromise: Rationale and applications

The method of Marquardt's compromise for handling the problem of near-singularity in nonlinear regression is modified by (a) utilizing methods from ridge regression for computing the Marquardt parameter k, and (b) utilizing a search method of interpolation-extrapolation for determining the step size. The performance of these modifications is evaluated on a number of standard test problems from the literature. It is shown that the Lawless-Wang estimate of k without the benefit of interpolation- extrapolation competes favorably with standard algorithms implemented in packaged programs. An application to a model of plasma insulin dynamics is also given. Comparison of the modified procedure with the Nelder-Mead simplex method is given.     

Effect of pressure-dependent viscosity on the exiting sheet thickness in the calendering of Newtonian fluids

An analysis of calendering for an incompressible Newtonian fluid flow, with pressure-dependent viscosity is studied theoretically under assumptions of isothermal conditions. We predict the influence of pressure-dependent viscosity on the exiting sheet thickness of the sheet of fluid from the gap. The dimensionless mass and momentum balance equations, which are based on lubrication theory, were solved for the velocity and pressure fields by using perturbation techniques, where the exiting sheet thickness represents an eigenvalue of the mathematical problem. When the above variables were obtained, the dimensionless exiting sheet thickness was determined by considering the influence of the pressure variations in the calendering process. Moreover, quantities of engineering interest are also calculated, which include the cylinder-separating force and power required to drive both cylinders in terms of the geometrical and kinematical parameters of the system. The results show that the inclusion of pressure-dependent viscosity effect increases the leave-off distance and consequently the dimensionless exiting sheet thickness in comparison with the case of pressure-independent viscosity.