守恒无伪解麦克斯韦方程间断元研究:(Ⅰ)一维和二维情况

考察了电、磁场分量分别基于不同近似函数空间展开的一维和二维Maxwell方程间断元求解方法。结合中心数值通量和电、磁场分量近似函数空间的不同组合,构造了各种间断元算子。通过用这些算子在规则和不规则网格上编码分析一维和二维金属腔的谐振模式,详细考察了算子的收敛和伪解支持性,并据此对基函数进行了优选。算子在时域和频域对谐振模式的计算结果彼此符合良好。优选的Maxwell方程间断元算子不仅同时具备能量守恒和免于伪解的特性,且无需引入辅助变量,为设计高品质Maxwell方程间断元算法和研发相关电磁场模拟软件提供了支撑。     

应用有限元和边界元法计算方坯软接触结晶器的电磁场

给出了用有限元和边界元相结合的方法计算方坯软接触结晶器内钢液电磁场分布的全过程,并对4面体单元基础上的Whitney边元素,H-Φ方程及边界积分方程的离散方法作了重点解释.采用有限元和边界元相结合的方法来计算电磁场的分布可以大大减少计算工作量和计算时间.自行开发了三维电磁场计算程序,将数值模拟结果与物理实际进行了比较.     

非均匀网格时域伪谱算法在超宽带技术中的应用

 与传统时域有限差分算法相比，采用以伪谱方法离散Maxwell微分方程为基础的时域伪谱(PSTD)算法计算大的电尺度电磁场时域问题，将大大提高计算效率，降低内存需求。为了拓宽PSTD算法的应用，近年来，基于网格插值方法的非均匀时域伪谱算法得到了发展。研究的重点是算法中非均匀网格技术的实现及其在时域瞬态脉冲电磁场模拟和高功率超宽带脉冲技术方面的应用。以高斯脉冲为激励源，用该算法计算了多层介质的反射和透射，并通过超宽带脉冲穿墙实验对这一方法的应用进行了验证。模拟和实验结果具有较好的一致性。     

电磁场2.5维时域间断有限元方法

介质沿空间固定方向均匀分布的结构在电磁导波器件中有十分广泛的应用,对这类器件的分析通常被称为2.5D电磁问题。利用器件在固定方向介质分布均匀的特点,将电磁场量沿该方向进行空间傅里叶变换,可以把对三维问题的分析转化为两维问题求解,从而极大地减小计算开销。针对传统基于差分的2.5D电磁场算法在弯曲形状逼近上有阶梯误差的缺陷,本文提出了基于三角形网格的2.5D时域间断有限元方法（DGTD）,并用它模拟了电偶极子与光纤的耦合效率和光子晶体光纤的色散特性。与基于规则网格的2.5D差分方法进行对比。结果表明,文中建立的2.5D DGTD方法对弯曲形状的模拟更加逼真,计算内存占用最大减少10.4%,计算精度最大相差0.011%,计算时间缩短74.9%,计算效率提高。     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

基于空间广义相干系数加权的医学合成孔径超声成像

为解决广义相干系数用于合成孔径成像中所存在的运算量大,图像对比度改善有限等问题,提出空间广义相干系数加权成像方法。该算法根据单个孔径成像结果之间的相干性来计算相干系数,通过加权空间合成进行成像。采用对FieldⅡ仿真点目标和吸声斑目标的数据进行成像表明,算法不仅使运算量减少N(N为阵元数)倍,而且相对于传统的广义相干系数算法,算法对散射点成像可提高信噪比7 dB,对于吸声斑成像可提高对比度3.2 dB。可见算法可以有效地提高成像速度和成像质量。      

基于LU-SGS的非结构弹簧网格迭代算法

针对运动间断拟合中需频繁更新网格点位置的特点,提出一种基于LU-SGS(lower-upper symmetricGauss-Seidel)迭代方法的非结构弹簧网格运动算法.根据弹簧网格原理构建与网格拓扑关系相对应的稀疏系数矩阵,将LU-SGS思想成功引入动网格迭代算法,并辅以合理的网格运动管理策略,实现动网格的快速迭代.研究表明,在非结构网格下,LU-SGS算法可以满足运动间断拟合的需求,在流场隐式时间推进时,仍能保证获得稳定解;与传统的SOR方法相比,计算时耗减少20%以上.     

一种改进的快速N-FINDR端元提取算法

为了解决传统N-FINDR算法计算量大,提取结果对噪声和初始端元选取敏感,且容易将异常点作为端元而造成误提取的问题,提出一种改进的快速N-FINDR端元提取算法.该方法通过光谱距离提取并去除高光谱图像中的冗余信息,减少N-FINDR提取端元的搜索范围,平滑噪声影响,并自适应剔除异常点,通过最大化光谱距离选取N-FINDR的初始端元,避免了随机选择的盲目性.采用合成数据和真实高光谱数据进行仿真分析并与现有算法进行对比,结果表明,本文算法能在噪声与奇异点干扰下正确提取端元,其提取效率和鲁棒性均优于现有算法.     

高阶精度CE/SE算法及其应用

对时-空守恒元解元算法(CE/SE)的网格设置做较大改进,提出一种新的六面体解元和元定义;同时在解元中对物理量进行高阶Taylor展开,给出一种在时间和空间上均具有高阶精度CE/SE算法.在此基础上,把新型的高阶精度CE/SE算法推广应用于高速流动捕捉激波间断、气相化学反应流动、计及固体动态效应的流体-弹塑性流动和非稳态多相不可压缩粘性流动中.数值实践表明,提出的新型网格结构上的高阶精度CE/SE算法具有算法简单、计算精度高、计算效率和计算效果好的优点,并大大改进和拓展了CE/SE算法的应用范围.     

基于非结构变形网格的间断装配法原理

在激波捕捉法计算得到的流场基础上采用辨识算法得到初始间断位置, 从ALE方程出发, 考虑离散几何守恒律, 采用变形网格和网格重构技术解决计算过程中间断运动和变形, 新旧网格之间流场采用高精度信息传递方法保持时间精度, 建立了基于非结构动网格技术的间断装配方法.通过激波管问题的二维模拟, 模拟了初始间断分解为激波和接触间断激波遇到固壁反射后与接触间断相交的非定常流动过程, 对这种新方法的基本原理进行了介绍.      

A Discontinuous Galerkin Method Based on a BGK Scheme for the Navier-Stokes Equations on Arbitrary Grids

A discontinuous Galerkin Method based on a Bhatnagar-Gross-Krook (BGK) formulation is presented for the solution of the compressible Navier-Stokes equations on arbitrary grids. The idea behind this approach is to combine the robustness of the BGK scheme with the accuracy of the DG methods in an effort to develop a more accurate, efficient, and robust method for numerical simulations of viscous flows in a wide range of flow regimes. Unlike the traditional discontinuous Galerkin methods, where a Local Discontinuous Galerkin (LDG) formulation is usually used to discretize the viscous fluxes in the Navier-Stokes equations, this DG method uses a BGK scheme to compute the fluxes which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous representation in the flux evaluation at a cell interface through a simple hybrid gas distribution function. The developed method is used to compute a variety of viscous flow problems on arbitrary grids. The numerical results obtained by this BGKDG method are extremely promising and encouraging in terms of both accuracy and robustness, indicating its ability and potential to become not just a competitive but simply a superior approach than the current available numerical methods.     

Time-Domain Numerical Solutions of Maxwell Interface Problems with Discontinuous Electromagnetic Waves

This paper is devoted to time domain numerical solutions of two-dimensional (2D) material interface problems governed by the transverse magnetic (TM) and transverse electric (TE) Maxwell's equations with discontinuous electromagnetic solutions. Due to the discontinuity in wave solutions across the interface, the usual numerical methods will converge slowly or even fail to converge. This calls for the development of advanced interface treatments for popular Maxwell solvers. We will investigate such interface treatments by considering two typical Maxwell solvers – one based on collocation formulation and the other based on Galerkin formulation. To restore the accuracy reduction of the collocation finite-difference time-domain (FDTD) algorithm near an interface, the physical jump conditions relating discontinuous wave solutions on both sides of the interface must be rigorously enforced. For this purpose, a novel matched interface and boundary (MIB) scheme is proposed in this work, in which new jump conditions are derived so that the discontinuous and staggered features of electric and magnetic field components can be accommodated. The resulting MIB time-domain (MIBTD) scheme satisfies the jump conditions locally and suppresses the staircase approximation errors completely over the Yee lattices. In the discontinuous Galerkin time-domain (DGTD) algorithm – a popular Galerkin Maxwell solver, a proper numerical flux can be designed to accurately capture the jumps in the electromagnetic waves across the interface and automatically preserves the discontinuity in the explicit time integration. The DGTD solution to Maxwell interface problems is explored in this work, by considering a nodal based high order discontinuous Galerkin method. In benchmark TM and TE tests with analytical solutions, both MIBTD and DGTD schemes achieve the second order of accuracy in solving circular interfaces. In comparison, the numerical convergence of the MIBTD method is slightly more uniform, while the DGTD method is more flexible and robust.     

Finite volume evolution Galerkin method for hyperbolic conservation laws with spatially varying flux functions

We present a generalization of the finite volume evolution Galerkin scheme [M. Luká?ová-Medvid’ová, J. Saibertov’a, G. Warnecke, Finite volume evolution Galerkin methods for nonlinear hyperbolic systems, J. Comp. Phys. (2002) 183 533– 562; M. Luká?ová-Medvid’ová, K.W. Morton, G. Warnecke, Finite volume evolution Galerkin (FVEG) methods for hyperbolic problems, SIAM J. Sci. Comput. (2004) 26 1–30] for hyperbolic systems with spatially varying flux functions. Our goal is to develop a genuinely multi-dimensional numerical scheme for wave propagation problems in a heterogeneous media. We illustrate our methodology for acoustic waves in a heterogeneous medium but the results can be generalized to more complex systems. The finite volume evolution Galerkin (FVEG) method is a predictor–corrector method combining the finite volume corrector step with the evolutionary predictor step. In order to evolve fluxes along the cell interfaces we use multi-dimensional approximate evolution operator. The latter is constructed using the theory of bicharacteristics under the assumption of spatially dependent wave speeds. To approximate heterogeneous medium a staggered grid approach is used. Several numerical experiments for wave propagation with continuous as well as discontinuous wave speeds confirm the robustness and reliability of the new FVEG scheme.     

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.     

The Locally Conservative Galerkin (LCG) Method — a Discontinuous Methodology Applied to a Continuous Framework

This paper presents a comprehensive overview of the element-wise locally conservative Galerkin (LCG) method. The LCG method was developed to find a method that had the advantages of the discontinuous Galerkin methods, without the large computational and memory requirements. The initial application of the method is discussed, to the simple scalar transient convection-diffusion equation, along with its extension to the Navier-Stokes equations utilising the Characteristic Based Split (CBS) scheme. The element-by-element solution approach removes the standard finite element assembly necessity, with an face flux providing continuity between these elemental subdomains. This face flux provides explicit local conservation and can be determined via a simple small post-processing calculation. The LCG method obtains a unique solution from the elemental contributions through the use of simple averaging. It is shown within this paper that the LCG method provides equivalent solutions to the continuous (global) Galerkin method for both steady state and transient solutions. Several numerical examples are provided to demonstrate the abilities of the LCG method.     

Analysis of an Implicit Fully Discrete Local Discontinuous Galerkin Method for the Time-Fractional Kdv Equation

In this paper, we consider a fully discrete local discontinuous Galerkin (LDG) finite element method for a time-fractional Korteweg-de Vries (KdV) equation. The method is based on a finite difference scheme in time and local discontinuous Galerkin methods in space. We show that our scheme is unconditionally stable and convergent through analysis. Numerical examples are shown to illustrate the efficiency and accuracy of our scheme.     

Numerical methods for solving one-dimensional cochlear models in the time domain

In this article, a robust numerical solution method for one-dimensional (1-D) cochlear models in the time domain is presented. The method has been designed particularly for models with a cochlear partition having nonlinear and active mechanical properties. The model equations are discretized with respect to the spatial variable by means of the principle of Galerkin to yield a system of ordinary differential equations in the time variable. To solve this system, several numerical integration methods concerning stability and computational performance are compared. The selected algorithm is based on a variable step size fourth-order Runge-Kutta scheme; it is shown to be both more stable and much more efficient than previously published numerical solution techniques.     

Brownian motion and thermophoresis effects on unsteady stagnation point flow of Eyring-Powell nanofluid: a Galerkin approach

This article concerns the analysis of an unsteady stagnation point flow of Eyring-Powell nanofluid over a stretching sheet. The influence of thermophoresis and Brownian motion is also considered in transport equations. The nonlinear ODE set is obtained from the governing nonlinear equations via suitable transformations. The numerical experiments are performed using the Galerkin scheme. A tabular form comparison analysis of outcomes attained via the Galerkin approach and numerical scheme (RK-4) is available to show the credibility of the Galerkin method. The numerical exploration is carried out for various governing parameters, namely, Brownian motion, steadiness, thermophoresis, stretching ratio, velocity slip, concentration slip, thermal slip, and fluid parameters, and Hartmann, Prandtl and Schmidt numbers. The velocity of fluid enhances with an increase in fluid and magnetic parameters for the case of opposing, but the behavior is reversed for assisting cases. The Brownian motion and thermophoresis parameters cause an increase in temperature for both cases (assisting and opposing). The Brownian motion parameter provides a drop-in concentration while an increase is noticed for the thermophoresis parameter. All the outcomes and the behavior of emerging parameters are illustrated graphically. The comparison analysis and graphical plots endorse the appropriateness of the Galerkin method. It is concluded that said method could be extended to other problems of a complex nature.     

Space–time discontinuous Galerkin finite element method for two-fluid flows

A novel numerical method for two-fluid flow computations is presented, which combines the space–time discontinuous Galerkin finite element discretization with the level set method and cut-cell based interface tracking. The space–time discontinuous Galerkin (STDG) finite element method offers high accuracy, an inherent ability to handle discontinuities and a very local stencil, making it relatively easy to combine with local hp-refinement. The front tracking is incorporated via cut-cell mesh refinement to ensure a sharp interface between the fluids. To compute the interface dynamics the level set method (LSM) is used because of its ability to deal with merging and breakup. Also, the LSM is easy to extend to higher dimensions. Small cells arising from the cut-cell refinement are merged to improve the stability and performance. The interface conditions are incorporated in the numerical flux at the interface and the STDG discretization ensures that the scheme is conservative as long as the numerical fluxes are conservative. The numerical method is applied to one and two dimensional two-fluid test problems using the Euler equations.     

Galerkin methods applied to some model equations for non-linear dispersive waves

The application of a dissipative Galerkin scheme to the numerical solution of the Korteweg de Vries (KdV) and Regularised Long Wave (RLW) equations, is investigated. The accuracy and stability of the proposed schemes is derived using a localised Fourier analysis. With cubic splines as basis functions, the errors in the numerical solutions of the KdV equation for different mesh-sizes and different amounts of dissipation is determined. It is shown that the Galerkin scheme for the RLW equation gives rise to much smaller errors (for a given mesh-size), and allows larger steps to be taken for the integrations in time (for a specified error tolerance). Also, the interaction of two solitons is compared for the KdV and RLW equations, and several differences in their behaviour are found.