共查询到20条相似文献，搜索用时 15 毫秒
1.
A modification of the Roe scheme called L2Roe for low dissipation low Mach Roe is presented. It reduces the dissipation of kinetic energy at the highest resolved wave numbers in a low Mach number test case of decaying isotropic turbulence. This is achieved by scaling the jumps in all discrete velocity components within the numerical flux function. An asymptotic analysis is used to show the correct pressure scaling at low Mach numbers and to identify the reduced numerical dissipation in that regime. Furthermore, the analysis allows a comparison with two other schemes that employ different scaling of discrete velocity jumps, namely, LMRoe and a method of Thornber et al. To this end, we present for the first time an asymptotic analysis of the last method. Numerical tests on cases ranging from low Mach number (M∞=0.001) to hypersonic (M∞=5) viscous flows are used to illustrate the differences between the methods and to show the correct behavior of L2Roe. No conflict is observed between the reduced numerical dissipation and the accuracy or stability of the scheme in any of the investigated test cases. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
2.
Numerical simulations of a very small amplitude acoustic wave interacting with a shock wave in a quasi-ID convergent-divergent nozzle is performed using an unstructured finite volume algorithm with piece-wise linear, least square reconstruction, Roe flux difference splitting, and second-order MacCormack time marching. First, the spatial accuracy of the algorithm is evaluated for steady flows with and without the normal shock by running the simulation with a sequence of successively finer meshes. Then the accuracy of the Roe flux difference splitting near the sonic transition point is examined for different reconstruction schemes. Finally, the unsteady numerical solutions with the acoustic perturbation are presented and compared with linear theory results. 相似文献
3.
4.
The HLLEM scheme is a popular contact and shear preserving approximate Riemann solver that is known to be plagued by various forms of numerical shock instability. In this paper, we clarify that the shock instability exhibited by this scheme is primarily triggered by the spurious activation of the antidiffusive terms present in the first and third Riemann flux components on the transverse interfaces adjoining the shock front due to numerical perturbations. These erroneously activated terms are shown to counteract the favorable damping mechanism provided by its inherent HLL-type diffusive terms, causing an unphysical variation of the conserved quantity ρu both along and across the numerical shock. To prevent this, two distinct strategies are proposed termed as S elective W ave M odification and A nti D iffusion C ontrol. The former focuses on enhancing the quantity of the favorable HLL-type dissipation available on these critical flux components by carefully increasing the magnitudes of certain nonlinear wave speed estimates, while the latter focuses on directly controlling the magnitude of these critical antidiffusive terms. A linear perturbation analysis is performed to gauge the effectiveness of these cures and to estimate a von Neumann–type stability bounds on the CFL number associated with their use. Results from a variety of classic shock instability test cases show that the proposed strategies are able to provide excellent shock stable solutions even on grids that are highly elongated across the shock front without compromising the accuracy on inviscid contact or shear dominated viscous flows. 相似文献
5.
6.
7.
A very simple linearization of the solution to the Riemann problem for the steady supersonic Euler equations is presented. When used locally in conjunction with the Godunov method, computing savings by a factor of about four relative to the use of exact Riemann solvers can be achieved. For severe flow regimes, however, the linearization loses accuracy and robustness. We then propose the use of a Riemann solver adaptation procedure. This retains the accuracy and robustness of the exact Riemann solver and the computational efficiency of the cheap linearized Riemann solver. Numerical results for two- and three-dimensional test problems are presented. 相似文献
8.
This paper describes a central‐difference interface‐capturing scheme applied to the prediction of flows with cavitation. Compressible cavitation schemes based on standard central‐difference solvers have been previously described, but the current scheme uses an incompressible formulation only previously implemented with an upwind solver. The central‐difference solver offers significant advantages in computational time compared with upwind schemes. Regions of cavitation are captured rather than tracked. This means that there is no need for complex tracking and reconstruction procedures for the interface of the cavitation region. The use of such schemes on an arbitrarily unstructured mesh is no more complicated than on its structured counterpart. Results for a number of test cases are presented, with comparisons made with both experimental data and other numerical solutions. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
9.
Leo Postma Jan K. L. van Beek Henk F. P. van den Boogaard Guus S. Stelling 《国际流体数值方法杂志》2013,71(10):1226-1237
Particle‐tracking models are often used for near field short‐term subgrid transport of substances. The consistency demand at the discrete level does not show up so dominantly for these applications. This demand refers to the use of a numerical advection scheme for particles that is fully compatible with the local mass conserving advection properties of the underlying flow field at the discrete level of that field. A noncompatible scheme will produce both local convergence and local divergence of particles in different parts of the model area and thus erroneous advection results and erroneous concentration patterns. This compatibility in particle tracking is especially important if smooth distributions over larger areas are modelled for longer times. These applications did not occur that often in the past because they require many particles and thus much computation time. These applications occur more frequently nowadays especially for environmental assessment such as for the modelling of transport of fish larvae growing during their journey in the model to juvenile stages. The advection scheme that is developed in this paper is shown to be exactly compatible with hydrodynamic flow fields computed by mass conserving curvilinear grid models. It is not only exact, it is fortunately also very simple to implement and fast, allowing for modelling a huge amount of particles with moderate computation time. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
10.
针对复杂管系内可压缩流体,基于有限体积法,采用HLLC(Harten-Lax-van Leer Contact)格式和黎曼求解器构建了有限控制体数值离散方法,引入虚拟节点用于连接有限控制体,借助虚拟节点给出控制体之间数值通量的计算格式,发展了一种管道内一维流动数值建模方法.针对含有分支管路的管系,在管道连接部位构建了分支管路拟一维流动数值计算模型.基于所发展的一维流动数值方法,建立了变径管道和含60°分支管道内流动计算模型,验证了该方法的收敛性和有效性;基于虚拟节点的数值格式处理变径管激波问题具有一定精度优势.研究了变径管和分支管模型中可压缩流体激波、稀疏波等的传播机理,分析了管径对相邻支管压力的影响,为工程管路设计提供了参考. 相似文献
11.
Wave equation models currently discretize the generalized wave continuity equation with a three‐time‐level scheme centered at k and the momentum equation with a two‐time‐level scheme centered at k+1/2; non‐linear terms are evaluated explicitly. However in highly non‐linear applications, the algorithm becomes unstable at even moderate Courant numbers. This paper examines an implicit treatment of the non‐linear terms using an iterative time‐marching algorithm. Depending on the domain, results from one‐dimensional experiments show up to a tenfold increase in stability and temporal accuracy. The sensitivity of stability to variations in the G‐parameter (a numerical weighting parameter in the generalized wave continuity equation) was examined; results show that the greatest increase in stability occurs with G/τ=2–50. In the one‐dimensional experiments, three different types of node spacing techniques—constant, variable, and LTEA (Localized Truncation Error Analysis)—were examined; stability is positively correlated to the uniformity of the node spacing. Lastly, a scaling analysis demonstrates that the magnitudes of the non‐linear terms are positively correlated to those that most influence stability, particularly the term containing the G‐parameter. It is evident that the new algorithm improves stability and temporal accuracy in a cost‐effective manner. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
12.
A. ARABSHAHI D. L. WHITFIELD 《International Journal of Computational Fluid Dynamics》2013,27(3-4):307-321
SUMMARY An implicit, upwind numerical scheme is presented for computing the unsteady transonic flowfield around complex aircraft configurations. This scheme solves the time-dependent Euler equations with a finite volume method that incorporates a high resolution Riemann solver to define the interface fluxes. A multiblock domain decomposition appproach is used to partition the computational domain into a completely arbitrary arrangement of blocks. However, this work is restricted to C1-continuous grid lines across block boundaries. Consequently, block-block interfaces are treated such as to mimic interior block communication, thus introducing no additional spatial differencing error. Computational results have been obtained for a complete wing-pylon-store configuration with the store in the captive and vertical launch positions. The numerical results for both stationary and dynamic grids have shown favorable agreement with experimental data. 相似文献
13.
To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme. The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1) and the National Natural Sciences Foundation of China (40175024 and 40035010) 相似文献
14.
Ben R. Hodges Francisco J. Rueda 《International Journal of Computational Fluid Dynamics》2013,27(9):593-607
The unsteady shallow-water equations for barotropic/baroclinic (free-surface/density-stratified) flows with non-linear coupling of density transport and momentum are solved using a family of two-time-level, semi-implicit predictor–corrector methods (PC2). The PC2 methods are a general family that includes the popular TRIM method for hydrostatic flows. PC2 is characterised by four ‘θ’ parameters controlling the time ‘n’ and ‘n + 1’ weighting of (1) free surface gradient, (2) predictor step, (3) baroclinic gradient and (4) density temporal interpolation. Stability of the non-linear coupling between momentum and density transport for PC2 is examined in the inviscid limit. Central difference and quadratic (QUICK) spatial interpolation for density are compared. Second-order temporal accuracy for both barotropic and baroclinic flows is simultaneously obtained with appropriate θ parameters, which has previously been shown to be impractical for TRIM. The 2nd-order PC2 method has near-neutral non-linear stability (slightly positive amplification factor) where linear theory predicts exactly neutral stability. QUICK is shown to be preferable to central difference spatial discretisation to reduce the amplification factor. Adjusting the baroclinic weighting or adding small artificial viscosities can stabilise the model for non-linear internal wave simulations. 相似文献
15.
This paper deals with the introduction of a multiresolution strategy into the semi‐intrusive scheme, recently introduced by the authors, aiming to propagate uncertainties in unsteady compressible fluid applications. The mathematical framework of the multiresolution setting is presented for the cell‐average case, and the coupling with the semi‐intrusive scheme is described from both the theoretical and algorithmic point‐of‐view. Some reference test cases are performed to demonstrate the convergence properties and the efficiency of the overall scheme: the linear advection problem for both smooth and discontinuous initial conditions, the inviscid Burgers equation, and an uncertain shock tube problem obtained by modifying the well‐known Sod shock problem. For all the cases, the convergence curves are computed with respect to semi‐analytical (exact) solutions. In the case of the shock tube problem, an original technique to obtain a reference highly‐accurate numerical stochastic solution has also been developed. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
16.
A perturbation finite volume (PFV) method for the convective-diffusion integral equation is developed in this paper. The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations, with the least nodes similar to the standard three-point schemes, that is, the number of the nodes needed is equal to unity plus the face-number of the control volume. For instance, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D linear and nonlinear problems, 2-D and 3-D flow model equations. Comparing with other standard three-point schemes, the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme (UDS). Its numerical accuracies are also higher than the second-order central scheme (CDS), the power-law scheme (PLS) and QUICK scheme. The project supported by the National Natural Science Foundation of China (10272106, 10372106) 相似文献
17.
A new numerical method for particle tracking (Lagrangian particle advection) on 2‐D unstructured grids with triangular cells is presented and tested. This method combines key attributes of published methods, including streamline closure for steady flows and local mass conservation (uniformity preservation). The subgrid‐scale velocity reconstruction is linear, and this linear velocity field is integrated analytically to obtain particle trajectories. A complete analytic solution to the 2‐D system of ordinary differential equations (ODEs) governing particle trajectories within a grid cell is provided. The analytic solution to the linear system of locally mass‐conserving constraints that must be enforced to obtain the coefficients in the ODEs is also provided. Numerical experiments are performed to demonstrate that the new method has substantial advantages in accuracy over previously published methods and that it does not suffer from unphysical particle clustering. The method can be used not only in particle‐tracking applications but also as part of a semi‐Lagrangian advection scheme.Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
18.
Mo‐Hong Chou 《国际流体数值方法杂志》2000,32(5):545-567
A numerical study is made of the unsteady two‐dimensional, incompressible flow past an impulsively started translating and rotating circular cylinder. The Reynolds number (Re) and the rotating‐to‐translating speed ratio (α) are two controlled parameters, and the influence of their different combinations on vortex shedding from the cylinder is investigated by the numerical scheme sketched below. Associated with the streamfunction (ψ)–vorticity (ω) formulation of the Navier–Stokes equations, the Poisson equation for ψ is solved by a Fourier/finite‐analytic, separation of variable approach. This approach allows one to attenuate the artificial far‐field boundary, and also yields a global conditioning on the wall vorticity in response to the no‐slip condition. As for the vorticity transport equation, spatial discretization is done by means of finite difference in which the convection terms are handled with the aid of an ENO (essentially non‐oscillatory)‐like data reconstruction process. Finally, the interior vorticity is updated by an explicit, second‐order Runge–Kutta method. Present computations fall into two categories. One with Re=103 and α≤3; the other with Re=104 and α≤2. Comparisons with other numerical or physical experiments are included. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
19.
An upstream flux‐splitting finite‐volume (UFF) scheme is proposed for the solutions of the 2D shallow water equations. In the framework of the finite‐volume method, the artificially upstream flux vector splitting method is employed to establish the numerical flux function for the local Riemann problem. Based on this algorithm, an UFF scheme without Jacobian matrix operation is developed. The proposed scheme satisfying entropy condition is extended to be second‐order‐accurate using the MUSCL approach. The proposed UFF scheme and its second‐order extension are verified through the simulations of four shallow water problems, including the 1D idealized dam breaking, the oblique hydraulic jump, the circular dam breaking, and the dam‐break experiment with 45° bend channel. Meanwhile, the numerical performance of the UFF scheme is compared with those of three well‐known upwind schemes, namely the Osher, Roe, and HLL schemes. It is demonstrated that the proposed scheme performs remarkably well for shallow water flows. The simulated results also show that the UFF scheme has superior overall numerical performances among the schemes tested. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献