基于level set的Eulerian-Lagrangian耦合方法及其应用

提出了一种基于水平集的Eulerian-Lagrangian耦合方法,其中Lagrangian方法采用相容显式有限元拉氏方法,Eulerian方法采用基于近似Riemann解的有限体积Eulerian方法,多介质界面处理采用新的水平集和Ghost方法计算. 给出了若干数值算例,包括激波管问题以及金属和气体的运动界面及其大变形问题,并分别与精确解和相容显式有限元拉氏方法的计算结果进行了对比. 数值结果表明,该方法计算结果正确,精度较高,能够准确捕捉物质界面,适用于处理大变形问题.      

改进的PPM格式及其在Riemann问题中的应用

在Collela和Skora提出的PPM格式基础之上,为消除舍入误差带来的影响,发展了一种改进的PPM格式。将改进的PPM格式结合Riemann近似解算子应用于求解Riemann问题。选取双膨胀波和Rayleigh-Taylor不稳定性问题作为算例进行了数值验证,并将改进的PPM格式与原始PPM格式的求解结果进行了对比分析。结果表明:改进的PPM格式相较于原始PPM格式的求解精度有了明显提升;计算结果更加合理,气泡发展曲线与解析解更为接近。     

奇异摄动方法在输电线非线性振动问题中的应用

同一系统内部快变量和慢变量的同时存在往往引发相异于一般系统的特殊效应,比如输电线的松弛振荡.本文推导了架空输电线具有初始垂度的非线性动力学模型,发现该模型是具有快慢变量耦合的数学模型,应用求解周期运动的奇异摄动方法,得到系统的近似解析解,考察了快慢变量对系统周期运动的影响规律.结果表明解析解较数值解略微偏小,但仍有很好的吻合度,说明本文结果的有效性和正确性.进一步计算表明,随着摄动方法应用过程中近似次数的增加,两解逐次接近.     

二维拉氏辐射流体力学人为解构造方法

人为构造解方法是复杂多物理过程耦合程序正确性验证的重要方法之一，适用于二维拉氏大变形网格的流体、辐射耦合人为解模型较为少见。针对拉氏辐射流体力学程序正确性验证的需要，从二维拉氏辐射流体力学方程组出发，基于坐标变换技术，给出了拉氏空间到欧氏空间的物理变量导数关系式，开展了辐射流体耦合的人为解构造方法研究，构造了一类质量方程无源项的二维人为解模型，并应用于非结构拉氏程序LAD2D辐射流体力学计算的正确性考核，为流体运动网格上的辐射扩散计算提供了一种有效手段。数值结果显示观测到的数值模拟收敛阶与理论分析一致。     

TVD自相容流体力学拉氏格式研究

流体力学数值模拟格式总体上可分为Eulerian(欧氏)、Lagrangian(拉氏)和ALE(Arbitrary LagrangianEulerian),TVD广泛应用于Eulerian格式。本文利用具有TVD保持性质的Runge-Kutta型时间离散方法,构造了流体力学Lagrangian(拉氏)自相容格式,应用von Neumann小扰动技术分析了该格式的稳定性,并进行了相应的数值模拟,较好地抑制了激波波后非物理振荡。     

求一类非线性振动微分方程的近似解的新方法

本文利用非线性振动微分方程中的非线性项是小量的特点,并利用变量置换,把原微分方程近似地变换成为常系数线性微分方程,求得了一级近似解.算例表明,一级近似解具有颇高的精度,且计算过程十分简单.     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

匀低速曲柄滑块机构近似解与精确解比较

 叙述和比较了曲柄滑块机构滑块运动方程的近似解和精确解. 导出了理想状态下滑块的运动方程的精确表达式;使用Taylor公式将精确表达式展开获得近似表达式,并保留一阶和二阶;通过数值算例分别比较了滑块的位移、速度和加速度表达式的近似解和精确解;观察数值结果:只有当曲柄和连杆长度比值较小时,近似解才接近精确解;它们的比值较大时,近似解和精确解的差别较大.     

多尺度法的设解形式之探讨

当采用多尺度法研究非线性振动问题时，经常会遇到不同情形下系统响应的设解形式不同的问题，不同的设解形式得到的结果是否相同，用哪种设解形式更为好一些，在其他文献中尚未见到有关讨论．本文针对一类具有平方和立方非线性项的单自由度和多自由度非线性系统，得到不同设解形式下的一次近似解和二次近似解，并与数值解相比较，发现两种设解形式的近似解与数值仿真解的解曲线吻合的很好，表明传统的各种设解形式在研究弱非线性系统时都有很好的近似性。     

奇异边界法模拟二维弹性力学问题

奇异边界法是一种新的边界型无网格数值离散方法。该方法使用基本解作为插值基函数,在继承传统边界型方法优点的同时,不需要费时费力的网格划分和奇异积分,数学简单,编程容易,是一个真正的无网格方法。为避免配置点与插值源点重合时带来的基本解源点奇异性,该方法提出了源点强度因子的概念,从而将边界型强格式方法的核心归结为求解源点强度因子。本文首次将该方法应用于求解平面弹性力学问题。数值算例表明,本文算法稳定,效率高,并可达到很高的计算精度。     

Harten-Lax-van Leer-contact (HLLC) approximation Riemann solver with elastic waves for one-dimensional elastic-plastic problems

A Harten-Lax-van Leer-contact (HLLC) approximate Riemann solver is built with elastic waves (HLLCE) for one-dimensional elastic-plastic flows with a hypoelastic constitutive model and the von Mises’ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the presented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves (TRRSE).     

An arbitrary Lagrangian–Eulerian method based on the adaptive Riemann solvers for general equations of state

Approximate or exact Riemann solvers play a key role in Godunov‐type methods. In this paper, three approximate Riemann solvers, the MFCAV, DKWZ and weak wave approximation method schemes, are investigated through numerical experiments, and their numerical features, such as the resolution for shock and contact waves, are analyzed and compared. Based on the analysis, two new adaptive Riemann solvers for general equations of state are proposed, which can resolve both shock and contact waves well. As a result, an ALE method based on the adaptive Riemann solvers is formulated. A number of numerical experiments show good performance of the adaptive solvers in resolving both shock waves and contact discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

近似黎曼解对高超声速气动热计算的影响研究

高超声速流场计算一般采用TVD型格式,这些格式中,大多采用了不同形式的近似黎曼解. 通过分析和数值验证,论述了激波捕捉格式中近似黎曼解的耗散性质,说明其对高超声速热流计算的影响. 数值实验证明,采用低耗散格式可大大提高热流计算精度,降低热流计算对网格的依赖程度,从而获得精确的热流数值解.      

An exact Riemann solver for detonation products

Numerical methods based upon the Riemann Problem are considered for solving the general initial-value problem for the Euler equations applied to real gases. Most of such methods use an approximate solution of the Riemann problem when real gases are involved. These approximate Riemann solvers do not yield always a good resolution of the flow field, especially for contact surfaces and expansion waves. Moreover, approximate Riemann solvers cannot produce exact solutions for the boundary points. In order to overcome these shortcomings, an exact solution of the Riemann problem is developed, valid for real gases. The method is applied to detonation products obeying a 5th order virial equation of state, in the shock-tube test case. Comparisons between our solver, as implemented in Random Choice Method, and finite difference methods, which do not employ a Riemann solver, are given.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

A Free-Lagrange method for unsteady compressible flow: simulation of a confined cylindrical blast wave

A Free-Lagrange numerical procedure for the simulation of two-dimensional inviscid compressible flow is described in detail. The unsteady Euler equations are solved on an unstructured Lagrangian grid based on a density-weighted Voronoi mesh. The flow solver is of the Godunov type, utilising either the HLLE (2 wave) approximate Riemann solver or the more recent HLLC (3 wave) variant, each adapted to the Lagrangian frame. Within each mesh cell, conserved properties are treated as piece-wise linear, and a slope limiter of the MUSCL type is used to give non-oscillatory behaviour with nominal second order accuracy in space. The solver is first order accurate in time. Modifications to the slope limiter to minimise grid and coordinate dependent effects are described. The performances of the HLLE and HLLC solvers are compared for two test problems; a one-dimensional shock tube and a two-dimensional blast wave confined within a rigid cylinder. The blast wave is initiated by impulsive heating of a gas column whose centreline is parallel to, and one half of the cylinder radius from, the axis of the cylinder. For the shock tube problem, both solvers predict shock and expansion waves in good agreement with theory. For the HLLE solver, contact resolution is poor, especially in the blast wave problem. The HLLC solver achieves near-exact contact capture in both problems. Received May 25, 1995 / Accepted September 11, 1995     

Dissipation matrix and artificial heat conduction for Godunov‐type schemes of compressible fluid flows

This paper aims to reassess the Riemann solver for compressible fluid flows in Lagrangian frame from the viewpoint of modified equation approach and provides a theoretical insight into dissipation mechanism. It is observed that numerical dissipation vanishes uniformly for the Godunov‐type schemes in the sense that associated dissipation matrix has zero determinant if an exact or approximate Riemann solver is used to construct numerical fluxes in the Lagrangian frame. This fact connects to some numerical defects such as the wall‐heating phenomenon and start‐up errors. To cure these numerical defects, a traditional numerical viscosity is added, as well as the artificial heat conduction is introduced via a simple passage of the Lax–Friedrichs type discretization of internal energy. Copyright © 2016 John Wiley & Sons, Ltd.     

Direct Riemann solvers for the time-dependent Euler equations

Approaches for finding direct, approximate solutions to the Riemann problem are presented. These result in three approximate Riemann solvers. Here we discuss the time-dependent Euler equations but the ideas are applicable to other systems. The approximate solvers are (i) assessed on local Riemann problems with exact solutions and (ii) used in conjunction with the Weighted Average Flux (WAF) method to solve the two-dimensional Euler equations numerically. The resulting numerical technique is assessed on a shock reflection problem. Comparison with experimental observation is carried out.     

HLLC scheme for the preconditioned pseudo-compressibility Navier–Stokes equations for incompressible viscous flows

A new HLLC (Harten-Lax-van leer contact) approximate Riemann solver with the preconditioning technique based on the pseudo-compressibility formulation for numerical simulation of the incompressible viscous flows has been proposed, which follows the HLLC Riemann solver (Harten, Lax and van Leer solver with contact resolution modified by Toro) for the compressible flow system. In the authors' previous work, the preconditioned Roe's Riemann solver is applied to the finite difference discretisation of the inviscid flux for incompressible flows. Although the Roe's Riemann solver is found to be an accurate and robust scheme in various numerical computations, the HLLC Riemann solver is more suitable for the pseudo-compressible Navier--Stokes equations, in which the inviscid flux vector is a non-homogeneous function of degree one of the flow field vector, and however the Roe's solver is restricted to the homogeneous systems. Numerical investigations have been performed in order to demonstrate the efficiency and accuracy of the present procedure in both two- and three-dimensional cases. The present results are found to be in good agreement with the exact solutions, existing numerical results and experimental data.     

An extension of Osher's Riemann solver for chemical and vibrational non-equilibrium gas flows

In this paper we study an extension of Osher's Riemann solver to mixtures of perfect gases whose equation of state is of the form encountered in hypersonic applications. As classically, one needs to compute the Riemann invariants of the system to evaluate Osher's numerical flux. For the case of interest here it is impossible in general to derive simple enough expressions which can lead to an efficient calculation of fluxes. The key point here is the definition of approximate Riemann invariants to alleviate this difficulty. Some of the properties of this new numerical flux are discussed. We give 1D and 2D applications to illustrate the robustness and capability of this new solver. We show by numerical examples that the main properties of Osher's solver are preserved; in particular, no entropy fix is needed even for hypersonic applications.     

拉格朗日高阶中心型守恒气体动力学格式

提出了拉格朗日高阶中心型守恒气体动力学格式。用产生于当前时刻子网格密度和当前时刻网格声速的子网格压力构造了子网格力,用加权本质无震荡方法构造的高阶子网格力构造了高阶空间通量,借助时间中点通量的泰勒展开完成了高阶时间通量离散,利用动量守恒条件使得格点速度以与网格面的数值通量相容的方式计算。编制了拉格朗日高阶中心型守恒气体动力学格式,对Saltzman活塞问题进行了数值模拟,数值结果表明,拉格朗日高阶中心型守恒气体动力学格式的有效性和精确性.