基于求解Riemann问题的界面处理方法

通过在界面处构造Riemann问题,根据流体的法向速度和压力在界面(接触间断)处连续的特性,利用Riemann问题的解不仅定义了ghost流体的值,而且对真实流体中邻近界面的点值进行了更新,使得在界面处的流体的状态满足接触间断的性质,给出了更加精确的界面边界条件,守恒误差分析表明该方法在界面计算过程中引入较小的误差.数值试验表明该方法能准确地捕捉界面和激波的位置.     

虚拟流体方法的设计原则

研究模拟可压缩多介质流的虚拟流体方法,建立一般状态方程下定义虚拟流体状态的基本原则.根据波系结构和使用的自由变量分别推导虚拟流体状态的定义方式.结果表明在这些方式下求解多介质Riemann问题理论上完全精确.进一步总结几种简单有效的虚拟流体方法,这些定义方式不依赖于虚拟流体区域可能产生的波系结构.其中一种类似于反射边界条件,只是界面速度需要首先精确预测出来.数值算例验证了研究结果的合理性.     

虚拟流体方法中界面处Riemann问题定义方式的改进

对RGFM中定义Riemann问题的方式进行改进,取距离界面适当远处的插值点处的状态作为Riemann问题的初值.并用数值算例对改进前后的RGFM进行比较.     

三角形网格上多介质流动问题界面处理方法

多介质流动问题的求解一般是在结构网格上实现,而三角形网格对于复杂计算区域具有更好的适应性,本文结合rGFM方法,给出三角形网格上多介质流动问题界面处理方法.利用level-set方法跟踪界面,在界面处构造Riemann问题,得到界面处流体准确的流动状态.通过定义界面边界条件,将多介质流动问题转化为单介质流动问题,利用高精度RKDG方法求解.采用多个算例验证该方法的稳健性和有效性,结果表明该方法能准确捕捉界面和激波的位置,保持界面清晰.     

Numerical method of the Riemann problem for two-dimensional multi-fluid flows with general equation of state

Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.     

A ghost fluid,level set methodology for simulating multiphase electrohydrodynamic flows with application to liquid fuel injection

In this paper, we present the development of a sharp numerical scheme for multiphase electrohydrodynamic (EHD) flows for a high electric Reynolds number regime. The electric potential Poisson equation contains EHD interface boundary conditions, which are implemented using the ghost fluid method (GFM). The GFM is also used to solve the pressure Poisson equation. The methods detailed here are integrated with state-of-the-art interface transport techniques and coupled to a robust, high order fully conservative finite difference Navier–Stokes solver. Test cases with exact or approximate analytic solutions are used to assess the robustness and accuracy of the EHD numerical scheme. The method is then applied to simulate a charged liquid kerosene jet.     

A cell-centered Lagrangian scheme with the preservation of symmetry and conservation properties for compressible fluid flows in two-dimensional cylindrical geometry

We develop a new cell-centered control volume Lagrangian scheme for solving Euler equations of compressible gas dynamics in cylindrical coordinates. The scheme is designed to be able to preserve one-dimensional spherical symmetry in a two-dimensional cylindrical geometry when computed on an equal-angle-zoned initial grid. Unlike many previous area-weighted schemes that possess the spherical symmetry property, our scheme is discretized on the true volume and it can preserve the conservation property for all the conserved variables including density, momentum and total energy. Several two-dimensional numerical examples in cylindrical coordinates are presented to demonstrate the performance of the scheme in terms of symmetry, accuracy and non-oscillatory properties.     

道路交通的流体物理模型与粒子仿真方法

针对受多种因素影响的复杂道路交通系统问题,基于颗粒动力学理论,结合传统的Lighthill-WhithamRichards物理模型,建立道路交通系统的流体物理模型,采用无网格粒子与网格相结合的方法进行数值仿真,并应用于典型道路交通问题的求解.在新模型中,将车辆比拟为硬颗粒,车辆的跟车比拟为颗粒间的碰撞相互作用,已知道路情况对驾驶员操作车辆的影响比拟为流-粒两相系统中的外部流体驱动力作用,不同车道间车辆的影响比拟为颗粒间的黏性作用,从而在颗粒动力学理论的基础上,推导建立了道路交通系统拟流体模型;引入光滑离散颗粒流体动力学(SDPH)对车辆系统模型进行离散,建立"SDPH车辆"与真实车辆之间的一一对应关系,再结合有限体积方法,对道路交通构建的双流体模型进行求解,建立求解交通流体物理模型的新型仿真方法.最后,采用所建立的模型和方法对车辆汇入以及机非混合对交通系统的影响过程进行了数值仿真,所得结果与实测值符合较好,表明新的模型和方法有效性好、可靠性高,为道路交通问题的解决提供了一条全新的途径.     

工程技术中的宽区域物态方程问题

物态方程参数是许多工程理论设计不可缺少的重要参数.文章介绍了实际工程问题的物态方程研究对象、物理问题、研究方法、主要研究阶段,提出了关于物态方程研究与发展的进一步构想.     

Mathematical and numerical studies on meshless methods for exterior unbounded domain problems

The method of fundamental solutions (MFS) is an efficient meshless method for solving boundary value problems in an exterior unbounded domain. The numerical solution obtained by the MFS is accurate, while the corresponding matrix equation is ill-conditioned. A modified MFS (MMFS) with the proper basis functions is proposed by the introduction of the modified Trefftz method (MTM). The concrete expressions of the corresponding condition numbers are given in mathematical forms and the solvability by these methods is mathematically proven. Thereby, the optimal parameter minimizing the condition number is also mathematically given. Numerical experiments show that the condition numbers of the matrices corresponding to the MTM and the MMFS are reduced and that the numerical solution by the MMFS is more accurate than the one by the conventional method.     

声子气的状态方程和声子气运动的守恒方程

根据爱因斯坦的狭义相对论中质能的等效关系，把固体(本文指非导体)晶格(原子)的质量分为晶格(原子)的静质量和晶格热振动能量的等效质量两个部分，后者就是固体中声子气的等效质量.晶格(原子)热振动的能量则分为晶格(原子)静质量具有的热能以及声子气质量具有的热能.基于固体的状态方程，导得了晶格静质量热振动的状态方程和声子气的状态方程.声子气在固体介质中的宏观运动就是热量在固体中的传递过程.建立了声子气运动的守恒方程组，分析表明，忽略惯性力时声子气的动量守恒方程就退化为傅里叶导热定律，阐明了傅里叶导热定律的物理本质是声子气驱动力与阻力的平衡方程.当热流密度很大惯性力不能忽略时，傅里叶导热定律不再适用. 关键词： 非傅里叶导热 声子气 声子气质量 状态方程 守恒方程     

无限深方势阱本征值和本征态的三种求解方法

一维无限深方势阱模型是量子力学理想模型,经典教材中势阱的边界一般取得比较特殊.或关于坐标原点具有对称性,或势阱左边界位于坐标原点.本文首先展示了如何利用3种方法求解一维任意边界无限深方势阱能量本征值和对应的本征态,不同方法得到的结果彼此之间等价,讨论分析了这3种方法的推导结果,然后得到关于一维任意边界无限深方势阱能量本征值和本征态的通式,从中比较容易看出这两个物理量均与阱宽有关,并且本征波函数与边界值有关,最后将一维结果拓展到二维和三维任意边界无限深方势阱情况.     

Generic van der Waals equation of state for polymers, modified free volume theory, and the self-diffusion coefficient of polymeric liquids

In this paper, a molecular theory of self-diffusion coefficient is developed for polymeric liquids (melts) on the basis of the integral equation theory for site-site pair correlation functions, the generic van der Waals equation of state, and the modified free volume theory of diffusion. The integral equations supply the pair correlation functions necessary for the generic van der Waals equation of state, which in turn makes it possible to calculate the self-diffusion coefficient on the basis of the modified free volume theory of diffusion. A random distribution is assumed for minimum free volumes for monomers along the chain in the melt. More specifically, a stretched exponential is taken for the distribution function. If the exponents of the distribution function for minimum free volumes for monomers are chosen suitably for linear polymer melts of N monomers, the N dependence of the self-diffusion coefficient is N⁻¹ for the small values of N, an exponent predicted by the Rouse theory, whereas in the range of 2.3?lnN?4.5 the N dependence smoothly crosses over to N⁻², which is reminiscent of the exponent by the reptation theory. However, for lnN?4.5 the N dependence of the self-diffusion coefficient differs from N⁻², but gives an N dependence, N^−2−δ(0<δ<1), consistent with experiment on polymer melts in the range. For polyethylene δ≈0.48 for the parameters chosen for the stretched exponential. Because the stretched exponential function contains undetermined parameters, the N dependence of diffusion becomes semiempirical, but once the parameters are chosen such that the N dependence of D can be successfully given for a polymer melt, the temperature dependence of the self-diffusion coefficient can be well predicted in comparison with experiment. The theory is satisfactorily tested against experimental and simulation data on the temperature dependence of D for polyethylene and polystyrene melts.     

Numerical methods for computing ground states and dynamics of nonlinear relativistic Hartree equation for boson stars

Efficient and accurate numerical methods are presented for computing ground states and dynamics of the three-dimensional (3D) nonlinear relativistic Hartree equation both without and with an external potential. This equation was derived recently for describing the mean field dynamics of boson stars. In its numerics, due to the appearance of pseudodifferential operator which is defined in phase space via symbol, spectral method is more suitable for the discretization in space than other numerical methods such as finite difference method, etc. For computing ground states, a backward Euler sine pseudospectral (BESP) method is proposed based on a gradient flow with discrete normalization; and respectively, for computing dynamics, a time-splitting sine pseudospectral (TSSP) method is presented based on a splitting technique to decouple the nonlinearity. Both BESP and TSSP are efficient in computation via discrete sine transform, and are of spectral accuracy in spatial discretization. TSSP is of second-order accuracy in temporal discretization and conserves the normalization in discretized level. In addition, when the external potential and initial data for dynamics are spherically symmetric, the original 3D problem collapses to a quasi-1D problem, for which both BESP and TSSP methods are extended successfully with a proper change of variables. Finally, extensive numerical results are reported to demonstrate the spectral accuracy of the methods and to show very interesting and complicated phenomena in the mean field dynamics of boson stars.     

TNT炸药和乳化炸药驱动飞板的通用状态方程特征线法研究

 在爆炸气体动力学的基础上,利用特征线差分法,对飞板的运动规律进行研究。差分过程中,基于稳定爆轰的基本假设,根据二维定常流理论,推导了通用状态方程的特征线相容关系。利用爆轰气体密度取代特征线相容关系中的马赫数,导出了与物质物态无关的通用状态方程特征线法,编写了爆轰产物作用下飞板内部和边界的差分计算程序。利用该程序计算了飞板在TNT炸药和乳化炸药爆轰作用下的运动参数,研究了在不同质量比下TNT炸药和乳化炸药爆轰驱动飞板的抛掷姿态曲线,并与Richter公式计算结果进行了对比研究。