多介质流动数值计算中的界面处理方法

 求解Riemann问题得到界面接触间断的流动状态,并以此构造带状区域的虚拟流体状态,对于多维问题设计了一种方便有效的算法。同时求解耦合的守恒形式欧拉方程组和非守恒界面捕捉方程,并用Level-Set函数捕捉界面,数值计算采用高分辨率MWENO格式。最后对可压缩多介质流动问题进行了数值模拟。     

用改进的耦合型Level Set方法计算一维双介质可压缩流动

用带有虚拟流体(Ghost Fluid)修正的Level Set方法计算了一维可压缩双介质流动,把描述流动的Euler方程和描述流体界面运动的Level Set方程耦合起来,得到一个整体的守恒律系统,应用高分辨率差分格式求解;为了解决流体界面附近的数值跳动问题,在界面附近引入了虚拟流体方法的Isobaric修正,并给出了算例.     

可压缩多介质流体数值模拟中的Level-Set间断跟踪方法

针对可压缩多介质流体的数值模拟,发展了一种Level-Set间断追踪技术,用LS(Level-Set)函数追踪激波和捕捉界面,用Riemann问题解构造带状区域内的虚拟流体状态,对物理量的外推方法、间断附近虚拟流体的构造、间断推进速度的计算等问题进行了研究.最后对可压缩多介质流体一维和二维守恒律方程组进行数值模拟,数值计算采用通量重构的高精度WENO格式,计算结果令人满意.     

双曲型守恒律的一种高精度TVD差分格式

构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.     

广义非线性Schrdinger方程的多辛格式与模方守恒律

通过正则变换,构造出广义非线性Schr(o)dinger方程的多辛方程组.对此多辛方程组,导出了一个新的模方守恒多辛格式.数值实验结果表明,多辛格式具有长时间的数值行为,且在保持模方守恒律方面优于蛙跳格式和辛欧拉中点格式.     

Mie-Grüneisen状态方程下多介质守恒型欧拉方程组的数值计算

针对Mie-Grüneisen状态方程下的多介质流场,将Euler方程中的部分复杂参数作为守恒量引入到新的方程中,使原方程结构得到简化.同时引入质量分数以表示流场中不同介质的变化,并对通量作近似处理,再利用MUSCL-TVD格式结合改进过的Euler方程进行求解,模拟Mie-Grüneisen状态方程下多介质相互作用的情况.结果表明,通过本方法既能使方程组满足守恒,又能较好地处理多种介质形成的间断面.     

SRLW方程的多辛Fourier谱格式及其守恒律

通过引进正则动量,将对称正则长波方程（简称SRLW方程）转化成多辛形式的方程组,它具有多辛守恒律;介绍了空间方向满足周期边界条件的函数的Fourier谱方法;对SRLW方程的多辛方程组在空间方向利用Fourer谱方法,时间方向上应用Euler中点格式离散,得到其多辛Fourier拟谱格式;证明此格式的一些离散守恒律．用此格式模拟了SRLW方程的单个孤立波,还模拟了多个孤立波的追赶、碰撞和分离过程．     

多介质流体力学计算的守恒型高分辨率格式

应用Lagrange坐标系下的守恒型格式计算多介质流体力学问题,在物质交界面附近采用一阶格式的数值通量,而在其余部分采用高分辨率格式的数值通量,不仅保持了高分辨率的良好性质,而且消除了一般的守恒型格式在界面附近所产生的震荡.     

求解双曲守恒律方程的高分辨率熵相容格式

为提高熵相容格式的精度,利用限制器机制构造高分辨率格式,将构造的通量限制器插入熵相容格式,得到一类高分辨率熵相容格式.构造Euler方程高分辨率熵相容格式时,对熵相容格式中的几个参数做简单调整,提高了接触间断处的分辨率.将所得格式的数值结果与熵相容格式的数值结果比较表明,构造的高分辨率熵相容格式具有稳健和基本无振荡等特性.     

不规则对接网格的分区求解计算

采用一种保持通量守恒的不规则对接网格分区求解中交界面耦合条件的计算方法, 结合有限体积法求解了Euler 方程, 无粘通量取用Van Leer 分裂格式, 构造了一种限制器以实现格式的二阶精度和TVD 性质, 并给出数值算例。     

Positivity-preserving high-resolution schemes for systems of conservation laws

A new class of flux-limited schemes for systems of conservation laws is presented that is both high-resolution and positivity-preserving. The schemes are obtained by extending the Steger–Warming method to second-order accuracy through the use of component-wise TVD flux limiters while ensuring that the coefficients of the discretization equation are positive. A coefficient is considered positive if it has all-positive eigenvalues and has the same eigenvectors as those of the convective flux Jacobian evaluated at the corresponding node. For certain systems of conservation laws, such as the Euler equations for instance, this condition is sufficient to guarantee positivity-preservation. The method proposed is advantaged over previous positivity-preserving flux-limited schemes by being capable to capture with high resolution all wave types (including contact discontinuities, shocks, and expansion fans). Several test cases are considered in which the Euler equations in generalized curvilinear coordinates are solved in 1D, 2D, and 3D. The test cases confirm that the proposed schemes are positivity-preserving while not being significantly more dissipative than the conventional TVD methods. The schemes are written in general matrix form and can be used to solve other systems of conservation laws, as long as they are homogeneous of degree one.     

On the implementation of WENO schemes for a class of polydisperse sedimentation models

The sedimentation of a polydisperse suspension of small rigid spheres of the same density, but which belong to a finite number of species (size classes), can be described by a spatially one-dimensional system of first-order, nonlinear, strongly coupled conservation laws. The unknowns are the volume fractions (concentrations) of each species as functions of depth and time. Typical solutions, e.g. for batch settling in a column, include discontinuities (kinematic shocks) separating areas of different composition. The accurate numerical approximation of these solutions is a challenge since closed-form eigenvalues and eigenvectors of the flux Jacobian are usually not available, and the characteristic fields are neither genuinely nonlinear nor linearly degenerate. However, the flux vectors associated with the widely used models by Masliyah, Lockett and Bassoon (MLB model) and Höfler and Schwarzer (HS model) give rise to Jacobians that are low-rank perturbations of a diagonal matrix. This property allows to apply a convenient hyperbolicity criterion that has become known as the “secular equation” [J. Anderson, A secular equation for the eigenvalues of a diagonal matrix perturbation, Lin. Alg. Appl. 246 (1996) 49–70]. This criterion was recently applied [R. Bürger, R. Donat, P. Mulet, C.A. Vega, Hyperbolicity analysis of polydisperse sedimentation models via a secular equation for the flux Jacobian, SIAM J. Appl. Math. 70 (2010) 2186–2213] to prove that the MLB and HS models are strictly hyperbolic under easily verifiable conditions, that their eigenvalues interlace with the velocities of the species that form the flux vector (so the velocities are good starting values for a root finder), and that the corresponding eigenvectors can be calculated with acceptable effort. In the present work, the newly available characteristic information is exploited for the implementation of characteristic-wise (spectral) weighted essentially non-oscillatory (WENO) schemes for the MLB and HS models. Numerical examples illustrate that WENO schemes which use this spectral information are superior in resolution, and even in efficiency for the same overall resolution, to component-wise WENO schemes.     

Hybrid Flux-Splitting Schemes for a Two-Phase Flow Model

In this paper we deal with the construction of hybrid flux-vector-splitting (FVS) schemes and flux-difference-splitting (FDS) schemes for a two-phase model for one-dimensional flow. The model consists of two mass conservation equations (one for each phase) and a common momentum equation. The complexity of this model, as far as numerical computation is concerned, is related to the fact that the flux cannot be expressed in terms of its conservative variables. This is the motivation for studying numerical schemes which are not based on (approximate) Riemann solvers and/or calculations of Jacobian matrix. This work concerns the extension of an FVS type scheme, a Van Leer type scheme, and an advection upstream splitting method (AUSM) type scheme to the current two-phase model. Our schemes are obtained through natural extensions of corresponding schemes studied by Y. Wada and M.-S. Liou (1997, SIAM J. Sci. Comput.18, 633–657) for Euler equations. We explore the various schemes for flow cases which involve both fast and slow transients. In particular, we demonstrate that the FVS scheme is able to capture fast-propagating acoustic waves in a monotone way, while it introduces an excessive numerical dissipation at volume fraction contact (steady and moving) discontinuities. On the other hand, the AUSM scheme gives accurate resolution of contact discontinuities but produces oscillatory approximations of acoustic waves. This motivates us to propose other hybrid FVS/FDS schemes obtained by removing numerical dissipation at contact discontinuities in the FVS and Van Leer schemes.     

γ-and α-Ce phase diagram: First-principle calculation

Controversies about the phase diagram for the isostructural γ ? α phase transition of cerium have long been standing out for several decades. To seek insight into the problems, high-precision equations of state(EOS) for γ-and α-cerium are constructed based on first-principle calculation. Versus previous works, the strong anharmonic effects of ion vibration and the variation of magnetism of γ-cerium are stressed. The new EOS generally agrees well with experimental data regarding thermodynamics, phase diagrams, and phase transitions. However, new EOS predicts that another part of phase boundary in pressure-temperature space may exist except for the commonly known boundary. In addition, the well-known critical point seems to be a critical point for γ-cerium to translate from a stable state to an unstable state.     

Splitting based finite volume schemes for ideal MHD equations

We design finite volume schemes for the equations of ideal magnetohydrodynamics (MHD) and based on splitting these equations into a fluid part and a magnetic induction part. The fluid part leads to an extended Euler system with magnetic forces as source terms. This set of equations are approximated by suitable two- and three-wave HLL solvers. The magnetic part is modeled by the magnetic induction equations which are approximated using stable upwind schemes devised in a recent paper [F. Fuchs, K.H. Karlsen, S. Mishra, N.H. Risebro, Stable upwind schemes for the Magnetic Induction equation. Math. Model. Num. Anal., Available on conservation laws preprint server, submitted for publication, URL: <http://www.math.ntnu.no/conservation/2007/029.html>]. These two sets of schemes can be combined either component by component, or by using an operator splitting procedure to obtain a finite volume scheme for the MHD equations. The resulting schemes are simple to design and implement. These schemes are compared with existing HLL type and Roe type schemes for MHD equations in a series of numerical experiments. These tests reveal that the proposed schemes are robust and have a greater numerical resolution than HLL type solvers, particularly in several space dimensions. In fact, the numerical resolution is comparable to that of the Roe scheme on most test problems with the computational cost being at the level of a HLL type solver. Furthermore, the schemes are remarkably stable even at very fine mesh resolutions and handle the divergence constraint efficiently with low divergence errors.     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

色散方程的一类具任意稳定性的显格式

本文对色散方程u_r=au_xxx构造了一类中间层包含四个结点,带两个参数m和θ的三层显式差分格式。当m和θ满足一定的关系时,其稳定性条件为|γ|≤(m+1)/(4(m-1))(|m|>1),从而当取m充分接近1时,可得到任意大的稳定性条件,并且保持截断误差阶不变。数值例子验证了理论分析的结果。     

Development of extended field Doppler velocimetry for turbomachinery applications

The development of a portable, single component Doppler global velocimetry (DGV) head, based around a wavelength-stabilised argonion laser and a fast digital image-processing system, is described. The normalised two-dimensional DGV image, in which intensities are linearly related to velocities, can be displayed and updated at the 25 Hz camera frame rate, greatly easing the problem of system alignment. The effect of each individual system component upon the velocity resolution achieved for the system as a whole is discussed, and correction factors are calculated to account for the finite aperture and field of view of real systems and for divergence of the illuminating light sheet. Axial velocities of up to 100 m/s in a straight duct flow have been measured, demonstrating an rms velocity resolution of 2.5 m/s. The potential of the technique for gas turbine applications has been demonstrated by measuring the position of a shock in a transonic flow. At a Mach number of 2.3 and mass flow rate of 0.79 kg/s the velocity change across the shock was measured to be approximately 130 m/s.     

高阶Adams公式的系数

当k的值不大时(例如k≤4),诸系数γiγi^*,β_kρ,β_kρ^*易于手工计算。但对大的k值(例如k≥5),手工计算已非易事。利用著名的计算机代数软件:REDUCE^[a],我们可以很容易地计算出诸系数γi,γi^*,β_kρ,β_kρ^*。     

中性粒子在Ioffe阱中的经典运动方程的Floquet解

研究中性粒子在Ioffe阱中近原点区域的囚禁时,阱中的磁场可以呈现出一种简明的形式.磁矩μ反平行于磁场的中性粒子在阱中与磁场发生相互作用,借助相互作用势,可以获得粒子在阱中的经典运动方程.在一定的条件下,采用迭代近似的目的,将方程演化为马丢方程的形式,利用传统的WKBJ目的可实现方程的近似求解.研究阱中中性粒子的囚禁问题时,感兴趣的是马丢方程的Floquet解,即周期为π,2π的全周期和半周期解,欲获得这种周期解,马丢方程中的参数λ和q必须满足一定的关系,为此必须选择阱的特定参数和粒子的特定初始条件,对这一问题进行了探索性的研究.