首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到10条相似文献,搜索用时 125 毫秒
1.
We present a finite volume scheme for solving shallow water equations with source term due to the bottom topography. The scheme has the following properties: it is high-order accurate in smooth wet regions, it correctly solves situations where dry areas are present, and it is well-balanced. The scheme is developed within a general nonconservative framework, and it is based on hyperbolic reconstructions of states. The treatment of wet/dry fronts is carried out by solving specific nonlinear Riemann problems at the corresponding intercells.  相似文献   

2.
张荣培  蔚喜军  崔霞  冯涛 《计算物理》2012,29(5):647-653
提出一种求解二维非平衡辐射扩散方程的数值方法.空间离散上采用加权间断Galerkin有限元方法,其中数值流量的构造采用一种新的加权平均;时间离散上采用隐-显积分因子方法,将扩散系数线性化,然后用积分因子方法求解间断Galerkin方法离散后的非线性常微分方程组.数值试验中在非结构网格上求解了多介质的辐射扩散方程.结果表明:对于强非线性和强耦合的非线性扩散方程组,该方法是一种非常有效的数值算法.  相似文献   

3.
Unstructured adaptive grid flow simulation is applied to the calculation of high-speed compressible flows of inert and reactive gas mixtures. In the present case, the flowfield is simulated using the 2-D Euler equations, which are discretized in a cell-centered finite volume procedure on unstructured triangular meshes. Interface fluxes are calculated by a Liou flux vector splitting scheme which has been adapted to an unstructured grid context by the authors. Physicochemical properties are functions of the local mixture composition, temperature, and pressure, which are computed using the CHEMKIN-II subroutines. Computational results are presented for the case of premixed hydrogen–air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave fronts. A parametric analysis of the influence of the adaptation parameters on the computed solution is performed.  相似文献   

4.
In the recent years, there has been an increasing interest in discontinuous Galerkin time domain (DGTD) methods for the solution of the unsteady Maxwell equations modeling electromagnetic wave propagation. One of the main features of DGTD methods is their ability to deal with unstructured meshes which are particularly well suited to the discretization of the geometrical details and heterogeneous media that characterize realistic propagation problems. Such DGTD methods most often rely on explicit time integration schemes and lead to block diagonal mass matrices. However, explicit DGTD methods are also constrained by a stability condition that can be very restrictive on highly refined meshes and when the local approximation relies on high order polynomial interpolation. An implicit time integration scheme is a natural way to obtain a time domain method which is unconditionally stable but at the expense of the inversion of a global linear system at each time step. A more viable approach consists of applying an implicit time integration scheme locally in the refined regions of the mesh while preserving an explicit time scheme in the complementary part, resulting in an hybrid explicit–implicit (or locally implicit) time integration strategy. In this paper, we report on our recent efforts towards the development of such a hybrid explicit–implicit DGTD method for solving the time domain Maxwell equations on unstructured simplicial meshes. Numerical experiments for 3D propagation problems in homogeneous and heterogeneous media illustrate the possibilities of the method for simulations involving locally refined meshes.  相似文献   

5.
We study the effects of surfaces on the kinetics of phase changes in Ising-type systems. If the surface effects can be modelled by a field which couples linearly to the local order parameter, the growth of wetting or drying layers occurs. The numerical solution of the corresponding time-dependent Ginzburg-Landau equation yields a temporally logarithmic growth for the thickness of a wetting (drying) layer growing from an unstable dry (wet) state. On the other hand, if one starts off with a metastable state, the radius of a supercritical plug (wet or dry) grows linearly in time, in accordance with recent experimental results.  相似文献   

6.
An adaptive implicit–explicit scheme for Direct Numerical Simulation (DNS) and Large-Eddy Simulation (LES) of compressible turbulent flows on unstructured grids is developed. The method uses a node-based finite-volume discretization with Summation-by-Parts (SBP) property, which, in conjunction with Simultaneous Approximation Terms (SAT) for imposing boundary conditions, leads to a linearly stable semi-discrete scheme. The solution is marched in time using an Implicit–Explicit Runge–Kutta (IMEX-RK) time-advancement scheme. A novel adaptive algorithm for splitting the system into implicit and explicit sets is developed. The method is validated using several canonical laminar and turbulent flows. Load balance for the new scheme is achieved by a dual-constraint, domain decomposition algorithm. The scalability and computational efficiency of the method is investigated, and memory savings compared with a fully implicit method is demonstrated. A notable reduction of computational costs compared to both fully implicit and fully explicit schemes is observed.  相似文献   

7.
A new approach for the construction of implicit subgrid-scale models for large-eddy simulation based on adaptive local deconvolution is proposed. An approximation of the unfiltered solution is obtained from a quasi-linear combination of local interpolation polynomials. The physical flux function is modeled by a suitable numerical flux function. The effective subgrid-scale model can be determined by a modified-differential equation analysis. Discretization parameters which determine the behavior of the implicit model in regions of developed turbulence can be adjusted so that a given explicit subgrid-scale model is recovered to leading order in filter width. Alternatively, improved discretization parameters can be found directly by evolutionary optimization. Computational results for stochastically forced and decaying Burgers turbulence are provided. An assessment of the computational experiments shows that results for a given explicit subgrid-scale model can be matched by computations with an implicit representation. A considerable improvement can be achieved if instead of the parameters matching an explicit model discretization parameters determined by evolutionary optimization are used.  相似文献   

8.
The proposed work concerns the numerical approximations of the shallow-water equations with varying topography. The main objective is to introduce an easy and systematic technique to enforce the well-balance property and to make the scheme able to deal with dry areas. To access such an issue, the derived numerical method is obtained by involving the free surface instead of the water height and this produces the scheme well-balanced independently from the numerical flux function associated with the homogeneous problem. As a consequence, we obtain an easy well-balanced scheme which preserves non-negative water height. When compared with the well-known hydrostatic reconstruction, the presented topography discretization does not involve any max function known to introduce some numerical errors as soon as the topography admits very strong variations or discontinuities. A second-order MUSCL accurate reconstruction is adopted. The proposed hydrostatic upwind scheme is next extended for considering 2D simulations performed over unstructured meshes. Several 1D and 2D numerical experiments are performed to exhibit the relevance of the scheme.  相似文献   

9.
非结构网格在平面叶栅内湍流流动数值模拟中的应用   总被引:2,自引:1,他引:1  
1控制方程和离散方程对二维非定常可压缩粘性流体,雷诺时均Navierstokes方程可写成如下积分形式:式中,0为计算区域;*D为区域边界;Q为守恒变量;F为通量,包括无粘通量F’和粘性通量F”。计算域采用三角形单元离散,三角形的边作为控制体边界,物理量置于单元中心。对每个三角形控制体,方程(1)可写成:·其中,Q;表示单元的守恒变量值;aC和Vi分别表示控制体的边界和面积,为计算式(2中右端值,将单元所有边界通量相加,最后可得如下方程:式中k(i)为单元边界,只,j为通过单元界面的通量,凸马为界面长度。2无粘通量本文采…  相似文献   

10.
This paper addresses the convergence properties of implicit numerical solution algorithms for nonlinear hyperbolic transport problems. It is shown that the Newton–Raphson (NR) method converges for any time step size, if the flux function is convex, concave, or linear, which is, in general, the case for CFD problems. In some problems, e.g., multiphase flow in porous media, the nonlinear flux function is S-shaped (not uniformly convex or concave); as a result, a standard NR iteration can diverge for large time steps, even if an implicit discretization scheme is used to solve the nonlinear system of equations. In practice, when such convergence difficulties are encountered, the current time step is cut, previous iterations are discarded, a smaller time step size is tried, and the NR process is repeated. The criteria for time step cutting and selection are usually based on heuristics that limit the allowable change in the solution over a time step and/or NR iteration. Here, we propose a simple modification to the NR iteration scheme for conservation laws with S-shaped flux functions that converges for any time step size. The new scheme allows one to choose the time step size based on accuracy consideration only without worrying about the convergence behavior of the nonlinear solver. The proposed method can be implemented in an existing simulator, e.g., for CO2 sequestration or reservoir flow modeling, quite easily. The numerical analysis is confirmed with simulation studies using various test cases of nonlinear multiphase transport in porous media. The analysis and numerical experiments demonstrate that the modified scheme allows for the use of arbitrarily large time steps for this class of problems.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号