共查询到20条相似文献,搜索用时 186 毫秒
1.
S. Adjerid J. E. Flaherty P. K. Moore Y. J. Wang 《Physica D: Nonlinear Phenomena》1992,60(1-4):94-111
We consider the adaptive solution of parabolic partial differential systems in one and two space dimensions by finite element procedures that automatically refine and coarsen computational meshes, vary the degree of the piecewise polynomial basis and, in one dimension, move the computational mesh. Two-dimensional meshes of triangular, quadrilateral, or a mixture of triangular and quadrilateral elements are generated using a finite quadtree procedure that is also used for data management. A posteriori estimates, used to control adaptive enrichment, are generated from the hierarchical polynomial basis. Temporal integration, within a method-of-lines framework, uses either backward difference methods or a variant of the singly implicit Runge-Kutta (SIRK) methods. A high-level user interface facilitates use of the adaptive software. 相似文献
2.
Zhiliang Xu Yingjie Liu Huijing Du Guang Lin Chi-Wang Shu 《Journal of computational physics》2011,230(17):6843-6865
We develop a new hierarchical reconstruction (HR) method and for limiting solutions of the discontinuous Galerkin and finite volume methods up to fourth order of accuracy without local characteristic decomposition for solving hyperbolic nonlinear conservation laws on triangular meshes. The new HR utilizes a set of point values when evaluating polynomials and remainders on neighboring cells, extending the technique introduced in Hu, Li and Tang [9]. The point-wise HR simplifies the implementation of the previous HR method which requires integration over neighboring cells and makes HR easier to extend to arbitrary meshes. We prove that the new point-wise HR method keeps the order of accuracy of the approximation polynomials. Numerical computations for scalar and system of nonlinear hyperbolic equations are performed on two-dimensional triangular meshes. We demonstrate that the new hierarchical reconstruction generates essentially non-oscillatory solutions for schemes up to fourth order on triangular meshes. 相似文献
3.
4.
在移动网格上构造一种反应流的动理学格式.首先利用BGK模型推导含化学反应的流体力学方程组,并利用其积分形式构造移动网格上离散格式,再利用自适应移动网格方法得到网格速度,最后利用时间精确的动理学数值方法构造数值通量,得到移动网格单元上新的物理量.一维与二维的数值实验表明这种格式同时具有高精度、高分辨率的特点. 相似文献
5.
A parallel adaptive mesh refinement (AMR) algorithm is proposed and applied to the prediction of steady turbulent non-premixed compressible combusting flows in three space dimensions. The parallel solution-adaptive algorithm solves the system of partial-differential equations governing turbulent compressible flows of reactive thermally perfect gaseous mixtures using a fully coupled finite-volume formulation on body-fitted multi-block hexahedral meshes. The compressible formulation adopted herein can readily accommodate large density variations and thermo-acoustic phenomena. A flexible block-based hierarchical data structure is used to maintain the connectivity of the solution blocks in the multi-block mesh and to facilitate automatic solution-directed mesh adaptation according to physics-based refinement criteria. For calculations of near-wall turbulence, an automatic near-wall treatment readily accommodates situations during adaptive mesh refinement where the mesh resolution may not be sufficient for directly calculating near-wall turbulence using the low-Reynolds-number formulation. Numerical results for turbulent diffusion flames, including cold- and hot-flow predictions for a bluff-body burner, are described and compared to available experimental data. The numerical results demonstrate the validity and potential of the parallel AMR approach for predicting fine-scale features of complex turbulent non-premixed flames. 相似文献
6.
Basil Bayati Philippe Chatelain Petros Koumoutsakos 《Journal of computational physics》2011,230(1):13-26
We present an algorithm for adaptive mesh refinement applied to mesoscopic stochastic simulations of spatially evolving reaction–diffusion processes. The transition rates for the diffusion process are derived on adaptive, locally refined structured meshes. Convergence of the diffusion process is presented and the fluctuations of the stochastic process are verified. Furthermore, a refinement criterion is proposed for the evolution of the adaptive mesh. The method is validated in simulations of reaction–diffusion processes as described by the Fisher–Kolmogorov and Gray–Scott equations. 相似文献
7.
We present a coupling procedure for two different types of particle methods for the Boltzmann and the Navier–Stokes equations. A variant of the DSMC method is applied to simulate the Boltzmann equation, whereas a meshfree Lagrangian particle method, similar to the SPH method, is used for simulations of the Navier–Stokes equations. An automatic domain decomposition approach is used with the help of a continuum breakdown criterion. We apply adaptive spatial and time meshes. The classical Sod’s 1D shock tube problem is solved for a large range of Knudsen numbers. Results from Boltzmann, Navier–Stokes and hybrid solvers are compared. The CPU time for the hybrid solver is 3–4 times faster than for the Boltzmann solver. 相似文献
8.
A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on hierarchical a posteriori error estimates. A global hierarchical error estimate is employed in this study to obtain reliable directional information of the solution. Instead of solving the global error problem exactly, which is costly in general, we solve it iteratively using the symmetric Gauß–Seidel method. Numerical results show that a few GS iterations are sufficient for obtaining a reasonably good approximation to the error for use in anisotropic mesh adaptation. The new method is compared with several strategies using local error estimators or recovered Hessians. Numerical results are presented for a selection of test examples and a mathematical model for heat conduction in a thermal battery with large orthotropic jumps in the material coefficients. 相似文献
9.
提出一种在自由重映移动网格下的广义黎曼问题方法模拟反应流.该方法基于显式的自由重映移动网格广义黎曼问题的解.为保证在时间和空间上的高精度,应用广义黎曼问题方法构造数值通量.为保证反应区的高分辨率,采用变分法生成自适应移动网格.该方法不仅能够保证网格质量,而且能有效地避免任意拉格朗日—欧拉方法中由于显式重映过程而带来的数值误差.包括CJ爆轰及不稳定爆轰的数值实验说明该格式的精确性和鲁棒性,证明这种移动网格下的二阶广义黎曼问题方法可以较好地捕捉反应流的间断与光滑结构. 相似文献
10.
11.
We propose self-adaptive finite element methods with error control for solving elliptic and electromagnetic problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. New error indicators are introduced to control the error due to non-body-fitted meshes. Flexible h-adaptive strategies are developed, which can be systematically extended to a large class of interface problems. Extensive numerical experiments are performed to support the theoretical results and to show the competitive behavior of the adaptive algorithm even for interfaces involving corner or tip singularities. 相似文献
12.
Yu Wang Songhe Song Zhijun Tan Desheng Wang 《Journal of computational physics》2009,228(17):6333-6348
This paper presents an adaptive method for variational curve smoothing based on level set implementation. A suitable cost functional is minimized via solving the derived Euler–Lagrangian equation, of which the discretization is conducted on unstructured triangular meshes by employing a simple and effective finite volume scheme. Through adaptive refinement of the mesh, the geometry features of the given curve can be well resolved in a cost-effective way. Various numerical experiments demonstrate the effectiveness and efficiency of the proposed approach. 相似文献
13.
研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量. 相似文献
14.
This work proposes a domain adaptive stochastic collocation approach for uncertainty quantification, suitable for effective handling of discontinuities or sharp variations in the random domain. The basic idea of the proposed methodology is to adaptively decompose the random domain into subdomains. Within each subdomain, a sparse grid interpolant is constructed using the classical Smolyak construction [S. Smolyak, Quadrature and interpolation formulas for tensor products of certain classes of functions, Soviet Math. Dokl. 4 (1963) 240–243], to approximate the stochastic solution locally. The adaptive strategy is governed by the hierarchical surpluses, which are computed as part of the interpolation procedure. These hierarchical surpluses then serve as an error indicator for each subdomain, and lead to subdivision whenever it becomes greater than a threshold value. The hierarchical surpluses also provide information about the more important dimensions, and accordingly the random elements can be split along those dimensions. The proposed adaptive approach is employed to quantify the effect of uncertainty in input parameters on the performance of micro-electromechanical systems (MEMS). Specifically, we study the effect of uncertain material properties and geometrical parameters on the pull-in behavior and actuation properties of a MEMS switch. Using the adaptive approach, we resolve the pull-in instability in MEMS switches. The results from the proposed approach are verified using Monte Carlo simulations and it is demonstrated that it computes the required statistics effectively. 相似文献
15.
《Combustion Theory and Modelling》2013,17(3):503-534
We present an adaptive finite-element method for the simulation of reactive flows in general domains including reliable error control for quantities of physical interest. On the basis of computable a posteriori error bounds the mesh is locally adjusted by a hierarchical feedback process yielding economical mesh-size distributions for a prescribed maximum number of cells. The key feature is the computational evaluation of the sensitivity of the error with respect to the local residuals by solving an associated global dual problem. This general approach is applied for computing three examples of flames in two dimensions with an increasing degree of complexity, namely a simple diffusion flame described by the flame-sheet model, an ozone decomposition flame and a methane flame in a complex burner geometry using a detailed reaction mechanism. 相似文献
16.
Lenka Dubcova Pavel Solin Glen Hansen HyeongKae Park 《Journal of computational physics》2011,230(4):1182-1197
Multiphysics solution challenges are legion within the field of nuclear reactor design and analysis. One major issue concerns the coupling between heat and neutron flow (neutronics) within the reactor assembly. These phenomena are usually very tightly interdependent, as large amounts of heat are quickly produced with an increase in fission events within the fuel, which raises the temperature that affects the neutron cross section of the fuel. Furthermore, there typically is a large diversity of time and spatial scales between mathematical models of heat and neutronics. Indeed, the different spatial resolution requirements often lead to the use of very different meshes for the two phenomena. As the equations are coupled, one must take care in exchanging solution data between them, or significant error can be introduced into the coupled problem. We propose a novel approach to the discretization of the coupled problem on different meshes based on an adaptive multimesh higher-order finite element method (hp-FEM), and compare it to popular interpolation and projection methods. We show that the multimesh hp-FEM method is significantly more accurate than the interpolation and projection approaches considered in this study. 相似文献
17.
将龙格库塔间断有限元方法(RDDG)与自适应方法相结合,求解三维欧拉方程.区域剖分采用非结构四面体网格,依据数值解的变化采用自适应技术对网格进行局部加密或粗化,减少总体网格数目,提高计算效率.给出四种自适应策略并分析不同自适应策略的优缺点.数值算例表明方法的有效性. 相似文献
18.
19.
任意网格重映的样条逼近算法 总被引:2,自引:1,他引:1
在大变形流体力学问题的数值模拟中,任何方法都必须考虑网格重分或网格自适应,只要改动网格就涉及重分,或自适应后从旧的、扭曲的网格到新网格的守恒量重映,包括质量、动量和能量.在研究样条函数逼近的基础上,给出一种物理量重映的对结构网格和非结构网格均适应的算法,并给出了数值结果. 相似文献
20.
九点格式中网格边上切向流的计算——节点自适应加权格式 总被引:3,自引:1,他引:2
针对网格扭曲的不同情形,直接考虑网格边上切向流的离散.基于扩散方程法向流连续的条件,给出离散法向流的构造,导出扭曲网格上九点计算格式中网格边上离散切向流的表达式,从而推导出加权系数的计算公式,适应于各种扭曲的网格.数值结果表明,与九点格式中节点量简单加权的方法相比,基于网格边离散切向流的节点自适应加权九点格式的精度有明显改进,迭代求解次数减少,计算效率明显提高. 相似文献