共查询到20条相似文献,搜索用时 15 毫秒
1.
2.
In this paper we introduce a high-order discontinuous Galerkin method for two-dimensional incompressible flow in the vorticity stream-function formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method. The stream function is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability. The method is efficient for inviscid or high Reynolds number flows. Optimal error estimates are proved and verified by numerical experiments. 相似文献
3.
Adrian Doicu Thomas Trautmann Franz Schreier Michael Hess 《Journal of Quantitative Spectroscopy & Radiative Transfer》2005,91(3):347-361
The finite element method is applied to the solution of the two-dimensional atmospheric radiative transfer. The analysis is mainly focussed on the derivation of the cell or element equation. The Galerkin method and several hybrid methods using the integral and finite difference form of the radiative transfer equation are employed to obtain the cell equation. The assembled system of equations relating the radiances at the lower and upper boundary of the domain is solved by a direct method. 相似文献
4.
5.
《中国物理 B》2015,(6)
A novel two-dimensional(2D) simulation method of positive corona current pulses is proposed.A control-volumebased finite element method(CV-FEM) is used to solve continuity equations,and the Galerkin finite element method(FEM) is used to solve Poisson's equation.In the proposed method,photoionization is considered by adopting an exact Helmholtz photoionization model.Furthermore,fully implicit discretization and variable time step are used to ensure the time-efficiency of the present method.Finally,the method is applied to a positive rod-plane corona problem.The numerical results are in agreement with the experimental results,and the validity of the proposed method is verified. 相似文献
6.
波的传播往往在复杂的地质结构中进行,如何有效地求解非均匀介质中的波动方程一直是研究的热点.本文将局部间断Galekin(local discontinuous Galerkin, LDG)方法引入到数值求解波动方程中.首先引入辅助变量,将二阶波动方程写成一阶偏微分方程组,然后对相应的线性化波动方程和伴随方程构造间断Galerkin格式;为了保证离散格式满足能量守恒,在单元边界上选取广义交替数值通量,理论证明该方法满足能量守恒性.在时间离散上,采用指数积分因子方法,为了提高计算效率,应用Krylov子空间方法近似指数矩阵与向量的乘积.数值实验中给出了带有精确解的算例,验证了LDG方法的数值精度和能量守恒性;此外,也考虑了非均匀介质和复杂计算区域的计算,结果表明LDG方法适合模拟具有复杂结构和多尺度结构介质中的传播. 相似文献
7.
We present spectral element (SE) and discontinuous Galerkin (DG) solutions of the Euler and compressible Navier–Stokes (NS) equations for stratified fluid flow which are of importance in nonhydrostatic mesoscale atmospheric modeling. We study three different forms of the governing equations using seven test cases. Three test cases involve flow over mountains which require the implementation of non-reflecting boundary conditions, while one test requires viscous terms (density current). Including viscous stresses into finite difference, finite element, or spectral element models poses no additional challenges; however, including these terms to either finite volume or discontinuous Galerkin models requires the introduction of additional machinery because these methods were originally designed for first-order operators. We use the local discontinuous Galerkin method to overcome this obstacle. The seven test cases show that all of our models yield good results. The main conclusion is that equation set 1 (non-conservation form) does not perform as well as sets 2 and 3 (conservation forms). For the density current (viscous), the SE and DG models using set 3 (mass and total energy) give less dissipative results than the other equation sets; based on these results we recommend set 3 for the development of future multiscale research codes. In addition, the fact that set 3 conserves both mass and energy up to machine precision motives us to pursue this equation set for the development of future mesoscale models. For the bubble and mountain tests, the DG models performed better. Based on these results and due to its conservation properties we recommend the DG method. In the worst case scenario, the DG models are 50% slower than the non-conservative SE models. In the best case scenario, the DG models are just as efficient as the conservative SE models. 相似文献
8.
9.
The Locally Conservative Galerkin (LCG) Method — a Discontinuous Methodology Applied to a Continuous Framework
下载免费PDF全文
![点击此处可从《advances in applied mathematics and mechanics.》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Rhodri L. T. Bevan Raoul vanLoon & Perumal Nithiarasu 《advances in applied mathematics and mechanics.》2009,1(3):319-340
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. 相似文献
10.
二维多介质可压缩流的RKDG有限元方法 总被引:1,自引:0,他引:1
应用RKDG(Runge-Kutta Discontinuous Galerkin)有限元方法、Level Set方法和Ghost Fluid方法数值模拟二维多介质可压缩流,其中Euler方程组、Level Set方程和重新初始化方程的空间离散采用DG(Discontinuous Galerkin)有限元方法,时间离散采用Runge-Kutta方法.对二维的气-气和气-液两相流进行了数值计算,得到了分辨率较高的计算结果. 相似文献
11.
Based on the complex variable moving least-square (CVMLS) approximation, the complex variable element-free Galerkin (CVEFG) method for two-dimensional viscoelasticity problems under the creep condition is presented in this paper. The Galerkin weak form is employed to obtain the equation system, and the penalty method is used to apply the essential boundary conditions, then the corresponding formulae of the CVEFG method for two-dimensional viscoelasticity problems under the creep condition are obtained. Compared with the element-free Galerkin (EFG) method, with the same node distribution, the CVEFG method has higher precision, and to obtain the similar precision, the CVEFG method has greater computational efficiency. Some numerical examples are given to demonstrate the validity and the efficiency of the method. 相似文献
12.
13.
Analysis of elastoplasticity problems using an improved complex variable element-free Galerkin method
下载免费PDF全文
![点击此处可从《中国物理 B》网站下载免费的PDF全文](/ch/ext_images/free.gif)
《中国物理 B》2015,(10)
In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin(ICVEFG) method is presented for two-dimensional(2D) elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods. 相似文献
14.
This second paper in a two part series describes the implementation of the finite element method for the solution of the problem of acoustic transmission through a non-uniform duct carrying a high speed subsonic compressible flow. A finite element scheme based on both the Galerkin method and the residual least squares method and with eight noded isoparametric elements is described. Multi-modal propagation is investigated by coupling of the solution in the duct non-uniform section to modal expansions in uniform sections. The accuracy of the finite element results for both the eigenvalue and transmission problems is assessed by comparison with exact solutions and with results from the method of weighted residuals in the form of a modified Galerkin method as introduced in Part I of this pair of papers. The results of calculations show that modal interactions, particularly in transmitted modes, become increasingly important with increasing duct flow Mach number. Power transmission coefficient calculations for the geometries studied reveal no indication of a linear basis for the phenomenon of subsonic acoustic choking. 相似文献
15.
L.H. Liu 《Journal of Quantitative Spectroscopy & Radiative Transfer》2007,103(3):536-544
To avoid the complicated and time-consuming computation of curved ray trajectories, a least-squares finite element method based on discrete ordinate equation is extended to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Four cases of radiative heat transfer are examined to verify this least-squares finite element method. Linear and nonlinear graded index are considered. The predicted dimensionless net radiative heat fluxes are determined by the least-squares finite element method and compared with the results obtained by other methods. The results show that the least-squares finite element method is stable and has a good accuracy in solving the multi-dimensional radiative transfer problem in a semitransparent graded index medium, while the Galerkin finite element method sometimes suffers from nonphysical oscillations. 相似文献
16.
A novel methodology for the simulation of 2D thermohaline double diffusive processes, driven by heterogeneous temperature and concentration fields in variable-density saturated porous media, is presented. The stream function is used to describe the flow field and it is defined in terms of mass flux. The partial differential equations governing system is given by the mass conservation equation of the fluid phase written in terms of the mass-based stream function, as well as by the advection–diffusion transport equations of the contaminant concentration and of the heat. The unknown variables are the stream function, the contaminant concentration and the temperature. The governing equations system is solved using a fractional time step procedure, splitting the convective components from the diffusive ones. In the case of existing scalar potential of the flow field, the convective components are solved using a finite volume marching in space and time (MAST) procedure; this solves a sequence of small systems of ordinary differential equations, one for each computational cell, according to the decreasing value of the scalar potential. In the case of variable-density groundwater transport problem, where a scalar potential of the flow field does not exist, a second MAST procedure has to be applied to solve again the ODEs according to the increasing value of a new function, called approximated potential. The diffusive components are solved using a standard Galerkin finite element method. The numerical scheme is validated using literature tests. 相似文献
17.
给出数值求解二维定常不可压Navier-Stokes型方程的非线性Galerkin有限元算法,并分析了数值解的正则性和收敛性,当粗网格参数H和细网格参数h满足关系式H=O(h1/2)时,该算法具有和Galerkin有限元算法同阶的收敛精度,然而在计算上比Galerkin有限元算法更为简单,可以节省可观的计算量.最后给出了数值试验,验证了上述结果。 相似文献
18.
We derive a tensor artificial viscosity suitable for use in a 2D or 3D unstructured arbitrary Lagrangian–Eulerian (ALE) hydrodynamics code. This work is similar in nature to that of Campbell and Shashkov [1]; however, our approach is based on a finite element discretization that is fundamentally different from the mimetic finite difference framework. The finite element point of view leads to novel insights as well as improved numerical results. We begin with a generalized tensor version of the Von Neumann–Richtmyer artificial viscosity, then convert it to a variational formulation and apply a Galerkin discretization process using high order Gaussian quadrature to obtain a generalized nodal force term and corresponding zonal heating (or shock entropy) term. This technique is modular and is therefore suitable for coupling to a traditional staggered grid discretization of the momentum and energy conservation laws; however, we motivate the use of such finite element approaches for discretizing each term in the Euler equations. We review the key properties that any artificial viscosity must possess and use these to formulate specific constraints on the total artificial viscosity force term as well as the artificial viscosity coefficient. We also show, that under certain simplifying assumptions, the two-dimensional scheme from [1] can be viewed as an under-integrated version of our finite element method. This equivalence holds on general distorted quadrilateral grids. Finally, we present computational results on some standard shock hydro test problems, as well as some more challenging problems, indicating the advantages of the new approach with respect to symmetry preservation for shock wave propagation over general grids. 相似文献
19.
Matthew F. Barone Irina Kalashnikova Daniel J. Segalman Heidi K. Thornquist 《Journal of computational physics》2009,228(6):1932-1946
The Galerkin projection procedure for construction of reduced order models of compressible flow is examined as an alternative discretization of the governing differential equations. The numerical stability of Galerkin models is shown to depend on the choice of inner product for the projection. For the linearized Euler equations, a symmetry transformation leads to a stable formulation for the inner product. Boundary conditions for compressible flow that preserve stability of the reduced order model are constructed. Preservation of stability for the discrete implementation of the Galerkin projection is made possible using a piecewise-smooth finite element basis. Stability of the reduced order model using this approach is demonstrated on several model problems, where a suitable approximation basis is generated using proper orthogonal decomposition of a transient computational fluid dynamics simulation. 相似文献
20.
We present a sparse Bayesian reconstruction method based on multiple types of a priori information for multispectral bioluminescence tomography (BLT). In the Bayesian approach, five kinds of a priori information are incorporated, reducing the ill-posedness of BLT. Specifically, source sparsity characteristic is considered to promote reconstruction results. Considering the computational burden in the multispectral case, a series of strategies is adopted to improve computational efficiency, such as optimal permissible source region strategy and node model of the finite element method. The performance of the proposed algorithm is validated by a heterogeneous three-dimensional (3D) micron scale computed tomography atlas and a mouse-shaped phantom. Reconstructed results demonstrate the feasibility and effectiveness of the proposed algorithm. 相似文献