discontinuous galerkin method for discontinuous temperature field problems

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允 许插值函数在单元边界上存在跳变. 在建立有限元方程时,通过在单元边界上引入数值通量 项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出. 数值算例 表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局 部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高. 该方法也为研 究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

水下声波传播过程的时域间断Galerkin有限元方法

本文构建了声压波动方程的改进时域间断Galerkin有限元方法.传统时域连续有限元方法在计算高梯度、强间断特征水中声波传播问题时往往会出现虚假数值振荡现象,这些数值振荡会影响正常波动的计算精度.为了解决这一问题,本文通过引入人工阻尼的方式构建了改进的时域间断Galerkin有限元方法,并针对具有高梯度、强间断特征的多障碍物复杂边界和层合液体介质声传播问题进行了计算.计算结果表明,与传统时域连续方法如N ew mark方法计算结果对比,所发展方法能较好地消除高梯度和强间断声压力波传播过程中虚假的数值振荡,具有较高的计算精度.问题的求解为进一步流固声耦合问题的研究奠定了基础.     

混凝土断裂的连续-非连续方法

采用有限元形函数作为单位分解函数,位移间断用富集节点的附加自由度表示,建立了允许在单元内部位移非连续的局部富集公式以表征混凝土的开裂区域．富集基函数由节点形函数和节点形函数与间断函数的乘积的并集构成．非连续位移的扩展路径完全与网格结构无关．不同于以非协调应变为基础的嵌入非连续模型,对单元的类型没有限制而且间断位移可以贯穿单元边界．局部富集思想与扩展有限元类似,但富集点自由度保持节点位移的物理意义不变,使相邻单元无需进行富集运算．在变分公式中引入混凝土粘结本构定律,推导了考虑断裂过程区非线性影响的基本方程．对混凝土粘结裂纹扩展的数值模拟说明了该计算方法的有效性．     

比例边界坐标插值方法在谱元法中的应用——无穷域Euler方程的数值模拟

将比例边界坐标插值方法引入谱元法, 构成比例边界谱单元, 对无穷域Euler方程进行数值模拟.阐述了比例边界谱单元的基本使用方法以及基于比例边界谱元的Runge-Kutta间断Galerkin方法求解Euler方程的过程;计算了无穷域圆柱和NACA0012翼型绕流问题, 并与已有结果进行了比较, 显示了计算结果的正确性.用基于比例边界谱元的间断Galerkin方法求解无穷域Euler方程时, 最多只需将求解域划分为2个子域, 避免了一般谱方法将求解域划分为9个或者27个子域的麻烦. 比例边界谱单元为无穷域Euler方程的直接求解提供了一个可供参考的方法.     

高精度线接触时变弹流润滑问题的数值模拟研究

针对一类简单的粗糙度模型,采用高阶间断有限元法对线接触时变弹流润滑问题进行高精度数值模拟,采用网格自适应方法减少计算量,用适合高阶格式的惩罚方法处理自由边界,对不同几何参数下该类粗糙度模型进行计算.结果表明,自适应高阶间断有限元方法可以捕捉到时变情况下相对微观的物理细节,如由粗糙度和时变效应所引起的额外压强峰,而采用传统低精度方法很难捕捉到.     

热冲击作用下层合皮肤组织热传导过程数值模拟

在急剧温度变化等强间断温度冲击作用下的生物层合组织非傅里叶热传导分析中,经典时域连续有限元方法(如Newmark等方法)会在波阵面以后的和层合组织界面附近的区域表现出强烈的数值振荡。这类数值振荡会影响问题求解精度,并带来较大不确定性。针对这类现象,本文发展了改进时域间断Galerkin有限元方法,进一步开展了相关问题的数值模拟。其控制方程的基本未知数(温度)及其时间导数在指定时间间隔内假设存在间断且独立插值。在有限元离散列式中引入比例刚度阵人工阻尼,以成功消除波前位置的虚假数值振荡行为。通过算例对比分析,相比Newmark方法和传统间断Galerkin方法,所提出的改进时域间断Galerkin有限元方法较好消除了波前、波后以及组织界面处的数值振荡,有效捕捉了波阵面的间断行为,提高了计算的精度。     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

各向异性位势平面问题的规则化边界元法

基于转化域方程为边界积分方程的极限定理及一个新颖的基本解分解技术, 建立间接变量规则化边界积分方程, 它有效地避免了奇异积分的直接计算. 与已有方法比,该方法不将问题变换为各向同性的问题去处理, 因而无需反演运算, 也有别于Galerkin方法, 无需计算重积分. 可计算任意边界位势梯度, 而不仅限于法向通量. 针对椭圆边界的边值问题, 提交一种精确单元来描述边界几何. 数值算例表明, 所提算法稳定且效率高, 所得数值结果与精确解吻合较好.      

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

改进时域间断Galerkin有限元方法在弹塑性波传播数值模拟中的应用

对于高频、强脉动荷载作用下的结构动力学波传播分析,对比于传统的时域算法,时域间断Galerkin方法能捕捉到波阵面的间断,有效得避免了由于间断引起的数值振荡。但时域间断方法却带来了波前面的虚假数值振荡。本文针对上述波前数值振荡的现象进行研究,通过引入人工阻尼的方法对时域间断Galerkin有限元方法进行进一步改进。数值结果表明,所发展的方法能够有效的滤掉强动荷载产生的波前数值振荡现象,同时降低了时域间断Galerkin方法的网格依赖性。     

Analysis of a discontinuous Galerkin approximation of the Stokes problem based on an artificial compressibility flux

In this work, we propose and analyse a discontinuous Galerkin (DG) method for the Stokes problem based on an artificial compressibility numerical flux. A crucial step in the definition of a DG method is the choice of the numerical fluxes, which affect both the accuracy and the order of convergence of the method. We propose here to treat the viscous and the inviscid terms separately. The former is discretized using the well‐known BRMPS method. For the latter, the problem is locally modified by adding an artificial compressibility term of the form (1/c²)(?p/?t) for the sole purpose of interface flux computation. The flux is obtained as the exact solution of a local Riemann problem. The analysis of the method extends the well‐established strategies for the DG discretization of the Laplacian to the resulting partially coercive problem. Copyright © 2007 John Wiley & Sons, Ltd.     

一类指数型间断函数的嵌入非连续等参有限元法

在嵌入非连续有限元的基本思想下,提出一类附加位移形函数———指数型间断函数,来模拟由于非连续结构,如裂纹和节理,所导致的位移不连续规律,该附加函数是以到间断处的垂直距离为自变量,且随距离的增大而呈指数衰减的函数.指数型间断函数具有在数学上的便于积分和求导的优点,且比阶梯间断函数更能反映实际破裂后的变形情况.本文用弱解形式推导了嵌入非连续有限元格式,编制了二维4节点和三维8节点的嵌入非连续等参有限元程序,并分别给出了算例.算例表明在模拟裂纹追踪时,指数型间断函数的嵌入非连续等参有限元法可行且有效.     

具有单元分裂功能的间断有限元方法

介绍一种能够模拟材料开裂过程的有限元方法,该方法引入单元分裂和界面分离技术,结合具体的破坏准则,模拟材料变形中的破坏,同时还可以方便地处理材料中任意分布的界面结构．以三点弯曲实验为例,通过数值模拟结果和实验数据的比较,验证了该方法的适用性．     

A kinetic energy preserving nodal discontinuous Galerkin spectral element method

In this work, we discuss the construction of a skew‐symmetric discontinuous Galerkin (DG) collocation spectral element approximation for the compressible Euler equations. Starting from the skew‐symmetric formulation of Morinishi, we mimic the continuous derivations on a discrete level to find a formulation for the conserved variables. In contrast to finite difference methods, DG formulations naturally have inter‐domain surface flux contributions due to the discontinuous nature of the approximation space. Thus, throughout the derivations we accurately track the influence of the surface fluxes to arrive at a consistent formulation also for the surface terms. The resulting novel skew‐symmetric method differs from the standard DG scheme by additional volume terms. Those volume terms have a special structure and basically represent the discretization error of the different product rules. We use the summation‐by‐parts (SBP) property of the Gauss–Lobatto‐based DG operator and show that the novel formulation is exactly conservative for the mass, momentum, and energy. Finally, an analysis of the kinetic energy balance of the standard DG discretization shows that because of aliasing errors, a nonzero transport source term in the evolution of the discrete kinetic energy mean value may lead to an inconsistent increase or decrease in contrast to the skew‐symmetric formulation. Furthermore, we derive a suitable interface flux that guarantees kinetic energy preservation in combination with the skew‐symmetric DG formulation. As all derivations require only the SBP property of the Gauss–Lobatto‐based DG collocation spectral element method operator and that the mass matrix is diagonal, all results for the surface terms can be directly applied in the context of multi‐domain diagonal norm SBP finite difference methods. Numerical experiments are conducted to demonstrate the theoretical findings. Copyright © 2014 John Wiley & Sons, Ltd.     

VISCOUS FLOWFIELD BASED ON DISCONTINUOUS BOUNDARY ELEMENT METHOD AND VORTEX METHOD

基于非协调边界元方法和涡方法的联合应用, 模拟了二维和三维黏性不可压缩流场. 计算中利用离散涡元对漩涡的产生、凝聚和输送过程进行模拟, 并将整体计算域分解为采用涡泡模拟的内部区域和用涡列模拟的数字边界层区域. 计算域中涡量场的拉伸和对流由Lagrangian涡方法模拟, 用随机走步模拟涡量场的扩散. 内部区域涡元涡量场速度由广义Biot-Savart公式计算, 势流场速度则采用非协调边界元方法计算. 非协调边界元将所有节点均取在光滑边界处, 从而避免了法向速度的不连续现象; 而对于系数矩阵不对称的大型边界元方程组,引入了非常高效的预处理循环型广义极小残余(the generalized minimum residual, GMRES)迭代算法, 使得边界元法的优势得到了充分发挥, 同时, 在内部涡元势流场计算中对近边界点采用了正则化算法, 该算法将奇异积分转化为沿单元围道上一系列线积分, 消除了势流计算中速度及速度梯度的奇异性. 二维、三维流场算例证明了所用方法的正确性, 也验证了该算法可以大幅度提高模拟精度和效率.     

A high-order discontinuous Galerkin method for simulating incompressible fluid-thermal-structural interaction problems

This paper presents a framework for the application of the discontinuous Galerkin(DG) finite element method to the multi-physics simulation of the solid thermal deformation interacting with incompressible flow problems in two-dimensions. Recent applications of the DG method are primarily for thermoelastic problems in a solid domain or fluid-structure interaction problems without heat transfer. Based on a recently published conjugate heat transfer solver, the incompressible Navier-Stokes equation, the fluid advection-diffusion equation, the Boussinesq term, the solid heat equation and the solid linear elastic equation are solved using an explicit DG formulation. A Dirichlet-Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux computed at interface quadrature points in the fluid-solid interface. Formal hp convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved for each solver. The algorithm is then further validated against several existing benchmark cases including the in-plane loaded square, the Timoshenko Beam, the laminated beam subject to thermal-loads and the lid-driven cavity with a flexible bottom wall. The computational effort demonstrates that for all cases the highest order accurate algorithm has several magnitudes lower error than the second-order schemes for a given computational effort. It is a strong justification for the development of such high order discretisations. The solver can be employed to predict thermal deformation of a structure due to convective and conductive heat transfer at low Mach, such as chip deformation on a printed circuit board, wave-guide structure optimization, thermoelectric cooler simulation, and optics mounting method verification.     

An internal penalty discontinuous Galerkin method for simulating conjugate heat transfer in a closed cavity

Using the discontinuous Galerkin (DG) method for conjugate heat transfer problems can provide improved accuracy close to the fluid‐solid interface, localizing the data exchange process, which may further enhance the convergence and stability of the entire computation. This paper presents a framework for the simulation of conjugate heat transfer problems using DG methods on unstructured grids. Based on an existing DG solver for the incompressible Navier‐Stokes equation, the fluid advection‐diffusion equation, Boussinesq term, and solid heat equation are introduced using an explicit DG formulation. A Dirichlet‐Neumann partitioning strategy has been implemented to achieve the data exchange process via the numerical flux of interface quadrature points in the fluid‐solid interface. Formal h and p convergence studies employing the method of manufactured solutions demonstrate that the expected order of accuracy is achieved. The algorithm is then further validated against 3 existing benchmark cases, including a thermally driven cavity, conjugate thermally driven cavity, and a thermally driven cavity with conducting solid, at Rayleigh numbers from 1000 to 100 000. The computational effort is documented in detail demonstrating clearly that, for all cases, the highest‐order accurate algorithm has several magnitudes lower error than first‐ or second‐order schemes for a given computational effort.     

A γ‐DGBGK scheme for compressible multi‐fluids

The paper presents a Discontinuous Galerkin γ‐BGK (γ‐DGBGK) method for compressible multicomponent flow simulations by coupling the discontinuous Galerkin method with a γ‐BGK scheme based on WENO limiters. In this γ‐DGBGK method, the construction of the flux in the DG method is based on the kinetic scheme which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous terms in the flux formulation at cell interfaces. WENO limiters are used to obtain uniform high‐order accuracy and sharp non‐oscillatory shock transition, and time accuracy obtained by integration for the flux function at the cell interface. Numerical examples in one and two space dimensions are presented to illustrate the robust and accuracy of the present scheme. Copyright © 2010 John Wiley & Sons, Ltd.     

HIGH-ORDER DISCONTINUOUS GALERKIN SOLUTION OF N-S EQUATIONS ON HYBRID MESH

针对层流NS方程发展了混合网格上的高阶间断有限元方法,给出了物面边界高阶近似的具体步骤以及近物面弯曲单元的处理方法。对数值离散产生的非线性方程组采用牛顿迭代进行求解,每个牛顿循环采用预处理广义最小余量法求解产生的大型稀疏线性系统。使用该方法得到了典型算例的数值结果,并跟前人的计算结果进行了比较。计算结果表明,混合网格上应用高阶间断有限元方法求解黏性流动具有很好的应用前景。