球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

Energy law preserving C finite element schemes for phase field models in two-phase flow computations

We use the idea in [33] to develop the energy law preserving method and compute the diffusive interface (phase-field) models of Allen–Cahn and Cahn–Hilliard type, respectively, governing the motion of two-phase incompressible flows. We discretize these two models using a C⁰ finite element in space and a modified midpoint scheme in time. To increase the stability in the pressure variable we treat the divergence free condition by a penalty formulation, under which the discrete energy law can still be derived for these diffusive interface models. Through an example we demonstrate that the energy law preserving method is beneficial for computing these multi-phase flow models. We also demonstrate that when applying the energy law preserving method to the model of Cahn–Hilliard type, un-physical interfacial oscillations may occur. We examine the source of such oscillations and a remedy is presented to eliminate the oscillations. A few two-phase incompressible flow examples are computed to show the good performance of our method.     

非规则形状介质内辐射-导热耦合传热的间断有限元求解

采用间断有限元法(discontinuous finite element method,DFEM)求解非规则形状介质内的辐射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热问题的高精度数值结果.和传统连续型有限元方法不同,DFEM将计算区域划分成相互独立的离散单元,形函数的构造、未知量的加权近似以及控制方程的求解均在每一个离散单元上进行.通过在单元之间施加迎风格式的数值通量,DFEM保证了整个计算区域的连续性,因此这种方法兼具良好的几何灵活性和局部守恒性.推导了辐射传输方程和能量扩散方程的射导热耦合传热问题,得到了典型非规则形状介质内辐射导热耦合传热的高精度数值结果.     

面光源矩形匀光照明反射器数学建模分析

通常反射器照明器件结构对面光源出射光线难以有效控制,用透镜又存在不可避免的色散问题。采用面光源置于反射器杯口,出光面背对目标面的结构,提出一种分点、分区控制光线的自由曲面反射器的设计方法。利用有限元法划分面光源,将微面光源视为点光源,根据矩形匀光照明配光条件,建立面光源与目标面的能量对应关系。根据微分几何理论和Snell定律,建立反射器的数学模型,只要给定参数,无需进行反馈优化,即可用Matlab软件给出自由曲面反射器的数值解。以台灯照明为例,用光学设计软件对其性能进行了模拟应用分析,结果表明,在有效照明区域内均匀度可达到80%以上,光能利用率达到90%。     

模拟生物组织冻结过程实验与分析

冷冻外科中组织冻结过程的分析对手术实施十分重要。本文建立了低温冷刀实验台,在模拟生物组织中进行冷冻实验,测量了冰球内某点的温度变化;并用有限元方法求解了建立在治法基础上的模拟生物组织冻结过程多维数学模型,计算结果与实验值符合较好;用该模型和方法计算了实验条件下冰球内的温度梯度变化和冷刀所需理论冷量;分析了不同的冷刀直径对冻结过程的影响。     

Diffusion front capturing schemes for a class of Fokker–Planck equations: Application to the relativistic heat equation

In this research work we introduce and analyze an explicit conservative finite difference scheme to approximate the solution of initial-boundary value problems for a class of limited diffusion Fokker–Planck equations under homogeneous Neumann boundary conditions. We show stability and positivity preserving property under a Courant–Friedrichs–Lewy parabolic time step restriction. We focus on the relativistic heat equation as a model problem of the mentioned limited diffusion Fokker–Planck equations. We analyze its dynamics and observe the presence of a singular flux and an implicit combination of nonlinear effects that include anisotropic diffusion and hyperbolic transport. We present numerical approximations of the solution of the relativistic heat equation for a set of examples in one and two dimensions including continuous initial data that develops jump discontinuities in finite time. We perform the numerical experiments through a class of explicit high order accurate conservative and stable numerical schemes and a semi-implicit nonlinear Crank–Nicolson type scheme.     

A least-squares/finite element method for the numerical solution of the Navier–Stokes-Cahn–Hilliard system modeling the motion of the contact line

In this article we discuss the numerical solution of the Navier–Stokes-Cahn–Hilliard system modeling the motion of the contact line separating two immiscible incompressible viscous fluids near a solid wall. The method we employ combines a finite element space approximation with a time discretization by operator-splitting. To solve the Cahn–Hilliard part of the problem, we use a least-squares/conjugate gradient method. We also show that the scheme has the total energy decaying in time property under certain conditions. Our numerical experiments indicate that the method discussed here is accurate, stable and efficient.     

非线性波动方程的新数值迭代方法

提出了一种新的求解非线性波动方程的数值迭代法,它是一种半解析的方法.与完全的数值计算方法扰法相比,它能够考虑各阶谐波的相互作用,且能够满足能量守恒定律.用它研究了非线性声波在液体中的传播性质,结果表明,在微扰法适用的声强范围内迭代法也适用,在微扰法不适用的一个较宽的声强范围内迭代法依然适用.     

An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

A stable and high-order accurate conjugate heat transfer problem

This paper analyzes well-posedness and stability of a conjugate heat transfer problem in one space dimension. We study a model problem for heat transfer between a fluid and a solid. The energy method is used to derive boundary and interface conditions that make the continuous problem well-posed and the semi-discrete problem stable. The numerical scheme is implemented using 2nd-, 3rd- and 4th-order finite difference operators on Summation-By-Parts (SBP) form. The boundary and interface conditions are implemented weakly. We investigate the spectrum of the spatial discretization to determine which type of coupling that gives attractive convergence properties. The rate of convergence is verified using the method of manufactured solutions.     

面向杆系结构的离散元异步长计算方法

为了提升离散元法处理连续介质问题的计算效率, 提出一种基于重叠颗粒的离散元分区异步长计算方法来处理连续介质问题。该方法采用颗粒分割将连续介质划分为若干子分区, 各分区采用velocity-Verlet积分格式求解运动方程。相邻分区通过重叠颗粒构成局部耦合区域, 边界数据传递过程中不涉及插值和截断过程。数值算例表明, 该方法在保证高精度的同时有效地降低了计算时间。在实际应用中可以有针对性地将连续介质划分为不同尺度颗粒的分区, 根据问题规模及分区颗粒尺度特性采用不同时间步长, 节省存储空间且大幅提升计算效率。     

An anisotropic mesh adaptation method for the finite element solution of heterogeneous anisotropic diffusion problems

Heterogeneous anisotropic diffusion problems arise in the various areas of science and engineering including plasma physics, petroleum engineering, and image processing. Standard numerical methods can produce spurious oscillations when they are used to solve those problems. A common approach to avoid this difficulty is to design a proper numerical scheme and/or a proper mesh so that the numerical solution validates the discrete counterpart (DMP) of the maximum principle satisfied by the continuous solution. A well known mesh condition for the DMP satisfaction by the linear finite element solution of isotropic diffusion problems is the non-obtuse angle condition that requires the dihedral angles of mesh elements to be non-obtuse. In this paper, a generalization of the condition, the so-called anisotropic non-obtuse angle condition, is developed for the finite element solution of heterogeneous anisotropic diffusion problems. The new condition is essentially the same as the existing one except that the dihedral angles are now measured in a metric depending on the diffusion matrix of the underlying problem. Several variants of the new condition are obtained. Based on one of them, two metric tensors for use in anisotropic mesh generation are developed to account for DMP satisfaction and the combination of DMP satisfaction and mesh adaptivity. Numerical examples are given to demonstrate the features of the linear finite element method for anisotropic meshes generated with the metric tensors.     

有限差分格式的一个通用色散-耗散条件

通过分析显式有限差分格式的数值色散和数值耗散,导出一个适于有限差分格式的通用色散-耗散条件.根据群速度和耗散率之间的物理关系,确定了用以抑制数值解中伪高波数波所需要的适度耗散.在以往发展的低耗散加权基本无振荡格式WENO-CU6-M2上的应用表明,该条件可用作优化线性或非线性有限差分格式的色散和耗散的通用指导准则.此外,满足色散-耗散条件的改进WENO-CU6-M2格式还可选作低分辨率数值模拟,以三维Taylor-Green涡向湍流转捩和自相似能量衰减问题展现了它的这种能力.与经典的动态Smagorinsky亚网格尺度模型相比,在Reynolds数Re=400~3000条件下,无黏和黏性Taylor-Green涡的数值模拟结果均得到明显改善.在保持激波捕捉特性同时,与最新的隐式大涡模拟模型的计算效果相当.      

A finite element method for analysis of vibration induced by maglev trains

This paper developed a finite element method to perform the maglev train–bridge–soil interaction analysis with rail irregularities. An efficient proportional integral (PI) scheme with only a simple equation is used to control the force of the maglev wheel, which is modeled as a contact node moving along a number of target nodes. The moving maglev vehicles are modeled as a combination of spring-damper elements, lumped mass and rigid links. The Newmark method with the Newton–Raphson method is then used to solve the nonlinear dynamic equation. The major advantage is that all the proposed procedures are standard in the finite element method. The analytic solution of maglev vehicles passing a Timoshenko beam was used to validate the current finite element method with good agreements. Moreover, a very large-scale finite element analysis using the proposed scheme was also tested in this paper.     

An energy law preserving C finite element scheme for simulating the kinematic effects in liquid crystal dynamics

In this paper, we use finite element methods to simulate the hydrodynamical systems governing the motions of nematic liquid crystals in a bounded domain Ω. We reformulate the original model in the weak form which is consistent with the continuous dissipative energy law for the flow and director fields in W^1,2+σ(Ω) (σ > 0 is an arbitrarily small number). This enables us to use convenient conformal C⁰ finite elements in solving the problem. Moreover, a discrete energy law is derived for a modified midpoint time discretization scheme. A fixed iterative method is used to solve the resulted nonlinear system so that a matrix free time evolution may be achieved and velocity and director variables may be solved separately. A number of hydrodynamical liquid crystal examples are computed to demonstrate the effects of the parameters and the performance of the method.     

平板大攻角平面绕流问题中涡生成的数值模拟

本文在涡团法的基础上结合有限元法及无限相似单元法在整个无界的流动区域上构造了Navier-Stokes方程的数值解,成功地克服了区域边界上凸角点邻域内解的奇性给数值计算带来的困难,同时克服了无穷远点给数值计算带来的困难。本文对平板大攻角平面绕流问题进行数值模拟,给出了初步的计算结果。     

蒸汽-冷流体接触冷凝流动的数值模拟

介绍了关于蒸汽-冷流体直接接触冷凝流动与传热的数值计算模型与部分研究结果。用Level Set方法确定蒸汽-冷流体接触界面的位置和形状,建立了对蒸汽和冷流体普遍适用的动量、能量和质量守恒方程,在能量和质量寺恒方程中增加了部分项用于计算蒸汽冷凝所产生的影响。用有限差分法在交错网格上离散控制方程,用Runge-Kutta法-五阶WENO组合格式求解Level Set输运方程,用压力修正的迭代Projection方法求解动量方程,而用SIMPLE方法求解温度控制方程。对算例的计算结果表明,本文所建立的数值计算模型能反映物理现象的宏观特性。根据计算结果,分析了本文模型的优缺点,并指出了今后改进的方向。     

Bifurcation Instability of sheet metal during spring-back

In automotive and home appliance industries, there are many complex-shaped sheet metal components which need to be fabricated in multiple stamping operations. For example, the manufacturing of an outer case of washing machine consists of stamping followed by a bending operation. After the first stage of the stamping process, a large amount of spring-back takes place, and therefore, it is difficult to proceed to the next stage of the bending process. In the stamping process of that kind of sheet component with low geometric constraint, the forming area is large compared to the forming depth. Therefore, the formed part is in an unstable state and is less geometrically constrained, which causes a large amount of spring-back. To investigate this phenomenon, finite element analyses are carried out. During a spring-back analysis after forming, bifurcation takes place and the finite element solution procedure using the Newton–Raphson scheme becomes unstable. To get a stable post-bifurcation solution, a bifurcation algorithm is introduced at the bifurcation point. The deformed shapes obtained from finite element analyses are in good agreement with the experimental data. From this study, it is shown that the bifurcation behaviour enlarges the spring-back and the degree of dimensional error. To obtain additional possible post-bifurcation solutions, non-bifurcation analyses using initial guesses obtained in a modal analysis are carried. For the initial guesses, lowed four eigenmodes are utilized. Finally, the post-bifurcation behaviour and spring-back amount are investigated for various process parameters including the forming depth, punch width and corner radius.     

层内突进传质扩散数学模型及求解

针对窜流型油藏的特点,抽象出物理原型,同时考虑纵向和横向扩散,建立层内突进传质扩散数学模型.应用拉普拉斯变换,求得层内突进传质扩散数学模型的解析解,并得到小段塞情况下的解析解.应用通用有限元分析软件,建立层内突进传质扩散的几何模型,并求得层内突进传质扩散数学模型的数值解.绘制层内突进传质扩散数学模型的浓度分布二维剖面及不同时间步的浓度变化剖面;通过无因次距离和无因次浓度关系及孔隙体积与无因次浓度关系可以看出,贝克莱特(Pe)数越大,峰值浓度越高,见剂时间越晚.通过解析解及数值解结合的方法,可加深对传质扩散本质的理解.     

Free vibrations of a mono-coupled periodic system

Vibration problems of periodic systems can be analyzed efficiently by means of the transfer matrix method. The frequency equation for the whole system is shown to be obtained in terms of the eigenvalues, or their natural logarithms, which are often called “propagation constants”, of the transfer matrix for a single periodic subsystem. In case of a mono-coupled system this frequency equation may be solved graphically by using the propagation constant curve, thereby saving a great deal of computational effort. Two types of mono-coupled systems are considered as numerical examples: a spring-mass oscillating system and a continuous Timoshenko beam resting on regularly spaced knife-edge supports. Depending on whether the transfer matrix is derived by an analytical procedure or by the finite element method, the numerical solutions become either exact or approximate.