共查询到20条相似文献,搜索用时 31 毫秒
1.
2.
流动数值模拟中一种并行自适应有限元算法 总被引:1,自引:0,他引:1
给出了一种流动数值模拟中的基于误差估算的并行网格自适应有限元算法.首先,以初网格上获得的当地事后误差估算值为权,应用递归谱对剖分方法划分初网格,使各子域上总体误差近似相等,以解决负载平衡问题.然后以误差值为判据对各子域内网格进行独立的自适应处理.最后应用基于粘接元的区域分裂法在非匹配的网格上求解N-S方程.区域分裂情形下N-S方程有限元解的误差估算则是广义Stokes问题误差估算方法的推广.为验证方法的可靠性,给出了不可压流经典算例的数值结果. 相似文献
3.
4.
In this work, two-level stabilized finite volume formulations for the
2D steady Navier-Stokes equations are considered.
These methods are based
on the local Gauss integration technique and the lowest equal-order
finite element pair. Moreover, the two-level
stabilized finite volume methods involve solving one small Navier-Stokes
problem on a coarse mesh with mesh size $H$, a large general Stokes problem for the Simple and
Oseen two-level stabilized finite volume methods on the fine mesh with mesh size $h$=$\mathcal{O}(H^2)$ or a large general Stokes equations for the Newton two-level stabilized finite
volume method on a fine mesh with mesh size $h$=$\mathcal{O}(|\log h|^{1/2}H^3)$.
These methods we studied provide an
approximate solution $(\widetilde{u}_h^v,\widetilde{p}_h^v)$ with the convergence rate of same order
as the standard stabilized finite volume method, which involve solving one large
nonlinear problem on a fine mesh with mesh size $h$. Hence, our methods
can save a large amount of computational time. 相似文献
5.
6.
7.
A. Craggs 《Journal of sound and vibration》1976,48(3):377-392
The approximate equations governing the forced harmonic motion of a damped acoustic system are set up by using a variational principle. Acoustic finite elements are then used in a computer program to study the transmission loss and insertion loss performance of some expansion chamber mufflers. The manner in which the equations are set up allows a number of input and output nodes, and two-dimensional effects involving the influence of transverse acoustic modes to be taken into account. Although only the simplest of elements and coarse mesh sizes are used the resulting accuracy of the solutions is extremely good; thus the method should be a viable one for studying the performance of more complicated mufflers, having variable cross-sections and internal energy dissipation. 相似文献
8.
从最小势能原理出发,使用变分-差分方法构造带有弯曲边梁的薄板的小挠度弯曲问题的差分格式,所得格式仅依赖板面网格结点,从而避免了由于引入虚拟网格结点而带来的问题;编制求解差分方程组的MATLAB程序,给出数值模拟结果. 相似文献
9.
Jan Nordström Sofia Eriksson Peter Eliasson 《Journal of computational physics》2012,231(14):4867-4884
We study the influence of different implementations of no-slip solid wall boundary conditions on the convergence to steady-state of the Navier–Stokes equations. The various approaches are investigated using the energy method and an eigenvalue analysis. It is shown that the weak implementation is superior and enhances the convergence to steady-state for coarse meshes. It is also demonstrated that all the stable approaches produce the same convergence rate as the mesh size goes to zero. The numerical results obtained by using a fully nonlinear finite volume solver support the theoretical findings from the linear analysis. 相似文献
10.
An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations
下载免费PDF全文
![点击此处可从《中国物理 B》网站下载免费的PDF全文](/ch/ext_images/free.gif)
The present paper deals with the numerical solution of the coupled Schrdinger-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. 相似文献
11.
High-Order Mesh Generation for Discontinuous Galerkin Methods Based on Elastic Deformation
下载免费PDF全文
![点击此处可从《advances in applied mathematics and mechanics.》网站下载免费的PDF全文](/ch/ext_images/free.gif)
Hongqiang Lu Kai Cao Lechao Bian & Yizhao Wu 《advances in applied mathematics and mechanics.》2016,8(4):693-702
In this paper, a high-order curved mesh generation method for Discontinuous
Galerkin methods is introduced. First, a regular mesh is generated. Second,
the solid surface is re-constructed using cubic polynomial. Third, the elastic governing
equations are solved using high-order finite element method to provide a fully
or partly curved grid. Numerical tests indicate that the intersection between element
boundaries can be avoided by carefully defining the elasticity modulus. 相似文献
12.
13.
A Method of Lines Based on Immersed Finite Elements for Parabolic Moving Interface Problems
下载免费PDF全文
![点击此处可从《advances in applied mathematics and mechanics.》网站下载免费的PDF全文](/ch/ext_images/free.gif)
This article extends the finite element method of lines to a parabolic
initial boundary value problem whose diffusion coefficient is discontinuous
across an interface that changes with respect to time. The method presented
here uses immersed finite element (IFE) functions for the discretization in
spatial variables that can be carried out over a fixed mesh (such as a
Cartesian mesh if desired), and this feature makes it possible to reduce
the parabolic equation to a system of ordinary differential equations (ODE)
through the usual semi-discretization procedure. Therefore, with a suitable
choice of the ODE solver, this method can reliably and efficiently solve a
parabolic moving interface problem over a fixed structured (Cartesian) mesh.
Numerical examples are presented to demonstrate features of this new method. 相似文献
14.
Jianhong Yang Lei Gang & Jianwei Yang 《advances in applied mathematics and mechanics.》2014,6(5):663-679
In this paper, we consider a two-scale stabilized finite volume method for the two-dimensional
stationary incompressible flow approximated by the lowest equal-order element pair $P_1-P_1$
which does not satisfy the inf-sup condition. The two-scale method consists of solving a small non-linear system
on the coarse mesh and then solving a linear Stokes equations on the fine mesh. Convergence of the optimal
order in the $H^1$-norm for velocity and the $L^2$-norm for pressure is obtained. The error analysis shows
there is the same convergence rate between the two-scale stabilized finite volume solution and the usual
stabilized finite volume solution on a fine mesh with relation $h =\mathcal{O}(H^2)$. Numerical experiments completely
confirm theoretic results. Therefore, this method presented in this paper is of practical importance in
scientific computation. 相似文献
15.
基于MOF界面重构的多物质ALE方法 总被引:2,自引:0,他引:2
提出一种基于MOF(Moment-of-Fluid)界面重构的多物质ALE(Arbitrary Lagrangian-Eulerian)方法.流体力学方程组采用相容有限元方法进行空间离散.提出一种新的二维子网格力学模型,用来计算混合网格中的物理量经过一个拉氏步后发生的变化,混合网格内的界面重构采用MOF方法.提出一种精确积分守恒重映方法.给出数值算例,如空气和水的Riemann问题,Dukowicz问题,水中强激波与空气泡相互作用问题等.结果表明,方法具有较高的精度,能够处理物质界面和网格的大变形问题. 相似文献
16.
Based on the integral form of the fluid dynamic equations, a finite volume kinetic scheme with arbitrary control volume and mesh velocity is developed. Different from the earlier unified moving mesh gas-kinetic method [C.Q. Jin, K. Xu, An unified moving grid gas-kinetic method in Eulerian space for viscous flow computation, J. Comput. Phys. 222 (2007) 155–175], the coupling of the fluid equations and geometrical conservation laws has been removed in order to make the scheme applicable for any quadrilateral or unstructured mesh rather than parallelogram in 2D case. Since a purely Lagrangian method is always associated with mesh entangling, in order to avoid computational collapsing in multidimensional flow simulation, the mesh velocity is constructed by considering both fluid velocity (Lagrangian methodology) and diffusive velocity (Regenerating Eulerian mesh function). Therefore, we obtain a generalized Arbitrary-Lagrangian–Eulerian (ALE) method by properly designing a mesh velocity instead of re-generating a new mesh after distortion. As a result, the remapping step to interpolate flow variables from old mesh to new mesh is avoided. The current method provides a general framework, which can be considered as a remapping-free ALE-type method. Since there is great freedom in choosing mesh velocity, in order to improve the accuracy and robustness of the method, the adaptive moving mesh method [H.Z. Tang, T. Tang, Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws, SIAM J. Numer. Anal. 41 (2003) 487–515] can be also used to construct a mesh velocity to concentrate mesh to regions with high flow gradients. 相似文献
17.
气相爆轰高阶中心差分-WENO组合格式自适应网格方法 总被引:1,自引:0,他引:1
研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点. 相似文献
18.
This paper is concerned with preconditioning the stiffness matrix resulting from finite element discretizations of Maxwell’s equations in the high frequency regime. The moving PML sweeping preconditioner, first introduced for the Helmholtz equation on a Cartesian finite difference grid, is generalized to an unstructured mesh with finite elements. The method dramatically reduces the number of GMRES iterations necessary for convergence, resulting in an almost linear complexity solver. Numerical examples including electromagnetic cloaking simulations are presented to demonstrate the efficiency of the proposed method. 相似文献
19.
给出了柱坐标下Euler方程数值边界条件的一种处理方法.径向第一个网格点设在距离中心半点位置上.根据相应物理量的特性,在中心附近进行边界延拓,使得内点的高精度差分格式可以同样应用在网格中心附近,从而无需单侧差分格式,保持了一致的高阶精度.对于周向边界,也建立了一种周期延拓方法,使得在周向所有节点处都能够采用同样的高精度格式离散,并进行了数值试验. 相似文献
20.
Axisymmetric laser welding of ceramics: comparison of experimental and finite element results 总被引:1,自引:0,他引:1
In this paper, we compare experimental data for a laser spot weld on a ceramic to the solution from an adaptive finite element model of the system. Our focus is on validating the finite element model, which necessarily includes numerous simplifications. We assume an axisymmetric geometry and flow profile, with a flat free surface. Buoyancy and surface tension drive the liquid motion in the molten ceramic pool beneath the laser, which is calculated using the axisymmetric forms of the continuity, momentum and energy equations. Latent heat, temperature-dependent material properties and radiation effects are all included in the formulation. These equations are solved with standard finite element techniques utilizing mesh relocation with a movement indicator based on solution gradients. Comparision with experimental data indicates that the numerical techniques used successfully predicted the depth and diameter of the actual ceramic weld pool. 相似文献