首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
    
Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

2.
    
Dexun  Fu  Yanwen  Ma 《Acta Mechanica Sinica》1986,2(2):100-108
A new scheme for solving the compressible Navier-Stokes equations is developed. For the inviscid portion of the equations the single step scheme used by the authors is factored according to the sign of the eigenvalues of Jacobian matrix. For the viscous portion of the equations a scheme corrected with operator addition is factored too. The scheme obtained has second order accuracy in time and in space and is used to solve two-dimensional problem. The numerical results of 2-D shock wave-boundary layer interaction are compared with experimental data.  相似文献   

3.
    
This paper presents a simple and efficient procedure developed for tracing discontinuities in flow fields. Numerical experiments are carried out to test the new sensor coupled with the associated nonlinear WENO dissipation filter developed to suppress the numerical dissipation. The tests show that, for a problem containing shocks and vortices, the implementation of the new sensor and an optimized WENO scheme can obtain a stable and highly resolved solution. The numerical experiments demonstrated that the new filter scheme performs efficiently both in parallel and serial running for one‐dimensional inviscid flow problems. Direct numerical simulation of a Mach 5 turbulent boundary layer over a flat plate was carried out to demonstrate the applicability of the scheme to the DNS practices. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

4.
    
Considering the importance of high‐order schemes implementation for the simulation of shock‐containing turbulent flows, the present work involves the assessment of a shock‐detecting sensor for filtering of high‐order compact finite‐difference schemes for simulation of this type of flows. To accomplish this, a sensor that controls the amount of numerical dissipation is applied to a sixth‐order compact scheme as well as a fourth‐order two‐register Runge–Kutta method for numerical simulation of various cases including inviscid and viscous shock–vortex and shock–mixing‐layer interactions. Detailed study is performed to investigate the performance of the sensor, that is, the effect of control parameters employed in the sensor are investigated in the long‐time integration. In addition, the effects of nonlinear weighting factors controlling the value of the second‐order and high‐order filters in fine and coarse non‐uniform grids are investigated. The results indicate the accuracy of the nonlinear filter along with the promising performance of the shock‐detecting sensor, which would pave the way for future simulations of turbulent flows containing shocks. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

5.
    
A space and time third‐order discontinuous Galerkin method based on a Hermite weighted essentially non‐oscillatory reconstruction is presented for the unsteady compressible Euler and Navier–Stokes equations. At each time step, a lower‐upper symmetric Gauss–Seidel preconditioned generalized minimal residual solver is used to solve the systems of linear equations arising from an explicit first stage, single diagonal coefficient, diagonally implicit Runge–Kutta time integration scheme. The performance of the developed method is assessed through a variety of unsteady flow problems. Numerical results indicate that this method is able to deliver the designed third‐order accuracy of convergence in both space and time, while requiring remarkably less storage than the standard third‐order discontinous Galerkin methods, and less computing time than the lower‐order discontinous Galerkin methods to achieve the same level of temporal accuracy for computing unsteady flow problems. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

6.
The existence of shock–turbulent boundary layer interactions lead to very complicated flow phenomena and pose a challenge for numerical simulation. In this paper, two turbulence models, the Baldwin–Lomax (B–L) model and the Johnson–King (J–K) model, which were originally developed for simple external flow simulation, are modified to model complex high-speed internal separated flows. The full Navier–Stokes solver used in this paper is based on a cell-centered finite volume method and multistepping time marching scheme. Both implicit residual smoothing and local time stepping techniques are incorporated to accelerate the convergence rate. To ensure the numerical stability with the present explicit scheme, a point-implicit treatment to the source term in the ordinary differential equation (ODE) of the J–K model has been developed and has proved to be very effective in modeling such a complex flow. An arc-bump channel flow case has been studied. Comparisons of computed results with experimental data show that the present solver, with the modified turbulence models, predicts the shock and the flow separation very well. The J–K model is found to predict the size of the separation bubble with a higher accuracy. © 1998 John Wiley & Sons, Ltd.  相似文献   

7.
    
The interaction between a curved shock wave and a compressible vortex is numerically investigated. The investigation concentrates on the local deformation of the shock structure due to the shock–vortex interaction. The essentially non‐oscillatory (ENO) scheme is used to solve the unsteady two‐dimensional Euler equations. A curved shock wave is obtained by the diffraction of an initially planar shock wave around a right‐angled corner and then allowed to interact with a strong compressible vortex superimposed on the flow. The same vortex affects the shock wave differently depending on the placement of the vortex because of the varying strength of the shock wave. This effect could range from a non‐symmetric deformation of the shock wave to a local disruption in the shock structure depending on the strength of the shock wave in the interaction region. This process leading to a local disruption in the shock structure is analyzed in detail. It is shown that such a disruption in the shock structure can be predicted by simple one‐dimensional considerations. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

8.
This paper presents the results of an experimental study of the unsteady nature of a hypersonic separated turbulent flow. The nomimal test conditions were a freestream Mach number of 7.8 and a unit Reynolds number of 3.5×107/m. The separated flow was generated using finite span forward facing steps. An array of flush mounted high spatial resolution and fast response platinum film resistance thermometers was used to make multi-channel measurements of the fluctuating surface heat trtansfer within the separated flow. Conditional sampling analysis of the signals shows that the root of separation shock wave consists of a series of compression wave extending over a streamwise length about one half of the incoming boundary layer thickness. The compression waves converge into a single leading shock beyond the boundary layer. The shock structure is unsteady and undergoes large-scale motion in the streamwise direction. The length scale of the motion is about 22 percent of the upstream influence length of the separation shock wave. There exists a wide band of frequency of oscillations of the shock system. Most of the frequencies are in the range of 1–3 kHz. The heat transfer fluctuates intermittently between the undisturbed level and the disturbed level within the range of motion of the separation shock wave. This intermittent phenomenon is considered as the consequence of the large-scale shock system oscillations. Downstream of the range of shock wave motion there is a separated region where the flow experiences continuous compression and no intermittency phenomenon is observed. The project supported by National Natural Science Foundation of China  相似文献   

9.
王世芬  李清泉 《力学学报》1991,23(4):426-432
本文给出高超音速湍流分离不稳定特性的实验研究结果。试验条件是:自由流马赫数为 7.8,单位长度雷诺数为 3.5×10~7/米。分离流场由有限展长前向台阶产生,并用有高空间分辨率和快速响应的一列平齐安装的铂膜电阻温度计和多通道系统测量其表面热流率脉动。信号的条件采样分析结果表明:分离激波的根部由一束压缩波构成,流向展长约二分之一来流边界层厚度,在边界层外汇聚成单一主激波。这种激波结构极其不稳定,出现大尺度运动,流向运动的尺度约为分离激波上游影响区域长度的22%。激波振荡频率为一宽频带,主要集中在 1~3 千赫。在分离激波运动区域,热流脉动呈间歇性,在无扰动和激波扰动间跳跃。可以认为这种间歇性是分离激波系统大尺度振荡的结果。在激波运动区域的下游为分离区,流体继续压缩,热流脉动无间歇。  相似文献   

10.
    
An exact similarity solution of the compressible‐flow Navier–Stokes equations is presented, which embeds supersonic, transonic, and subsonic regions. Describing the viscous and heat‐conducting high‐gradient flow in a shock wave, the solution accommodates non‐linear temperature‐dependent viscosity as well as heat‐conduction coefficients and provides the variation of all the flow variables and their derivatives. Also presented are methods to obtain time‐dependent and/or multi‐dimensional solutions as well as verification benchmarks of increasing severity. Comparisons between the developed analytical solution and CFD solutions of the Navier–Stokes equations, with determination of convergence rates and orders of accuracy of these solutions, illustrate the utility of the developed exact solution for verification purposes. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

11.
Air‐breathing propulsion systems for high‐speed space travel applications are studied. Ramjets and scramjets have been identified as potential candidates. The flow inlets of such systems are modelled with a simulation that can predict all complex inlet flow features, including shock due to forebody, multiple shock reflections, normal shock, shock–boundary layer interaction and associated separation for two‐dimensional and axisymmetric inlets. Computed values are in good agreement with experimental data. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

12.
    
A dual‐time implicit mesh‐less scheme is presented for calculation of compressible inviscid flow equations. The Taylor series least‐square method is used for approximation of spatial derivatives at each node which leads to a central difference discretization. Several convergence acceleration techniques such as local time stepping, enthalpy damping and residual smoothing are adopted in this approach. The capabilities of the method are demonstrated by flow computations around single and multi‐element airfoils at subsonic, transonic and supersonic flow conditions. Results are presented which indicate good agreements with other reliable finite‐volume results. The computational time is considerably reduced when using the proposed mesh‐less method compared with the explicit mesh‐less and finite‐volume schemes using the same point distributions. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

13.
Turbulence and aeroacoustic noise high‐order accurate schemes are required, and preferred, for solving complex flow fields with multi‐scale structures. In this paper a super compact finite difference method (SCFDM) is presented, the accuracy is analysed and the method is compared with a sixth‐order traditional and compact finite difference approximation. The comparison shows that the sixth‐order accurate super compact method has higher resolving efficiency. The sixth‐order super compact method, with a three‐stage Runge–Kutta method for approximation of the compressible Navier–Stokes equations, is used to solve the complex flow structures induced by vortex–shock interactions. The basic nature of the near‐field sound generated by interaction is studied. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

14.
    
In order to solve compressible turbulent flow problems, this study focuses on incorporating the Spalart–Allmaras turbulence model into gas‐kinetic BGK (Bhatnagar–Gross–Krook) scheme. The Spalart–Allmaras turbulence model is solved using finite difference discretization. The variables on the cell interface are interpolated via the van Leer limiter in the reconstruction stage. Simulation of subsonic and transonic flow over a NACA0012 airfoil has been implemented using two‐dimensional body‐fitted grids. The numerical results obtained appear in good agreement with the AGARD results, demonstrating the effectiveness and usefulness of the strategy of coupling the Spalart–Allmaras turbulence model with the BGK scheme for compressible turbulent flow simulation. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

15.
    
Matrix‐free implicit treatments are now commonly used for computing compressible flow problems: a reduced cost per iteration and low‐memory requirements are their most attractive features. This paper explains how it is possible to preserve these features for all‐speed flows, in spite of the use of a low‐Mach preconditioning matrix. The proposed approach exploits a particular property of a widely used low‐Mach preconditioner proposed by Turkel. Its efficiency is demonstrated on some steady and unsteady applications. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

16.
    
A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.  相似文献   

17.
  总被引:1,自引:0,他引:1  
We describe the formulation of the gas dynamics and high‐temperature thermochemical modules of the Eilmer code, an open‐source Navier–Stokes solver for transient compressible flow in two and three dimensions. The core gas dynamics formulation is based on finite‐volume cells, and the thermochemical effects are handled with specialised updating schemes that are coupled into the overall time‐stepping scheme. Verification of the code is explored via a number of case studies that use analytic and semi‐analytic solutions as comparison. These include both smooth and shocked flows and are used to demonstrate the order of spatial accuracy of the code. Cases include manufactured solutions for rather abstract inviscid and viscous flow, an idealised detonation wave supported by a curved body, and the transient flow of an idealised but high‐performance shock tube. Validation of the inviscid gas dynamics and thermochemical models is then explored using data from a selection of experimental studies. These studies include ballistic range experiments with chemically‐inert noble gases and high‐temperature chemically‐reacting air. These comparisons show that the code performs well and they provide a lesson in considering a range of experimental data rather than relying upon isolated data points for validation. These verification and validation cases are described in full detail and will be useful for other code developers of high‐temperature compressible flow solvers. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

18.
    
Numerical oscillation has been an open problem for high‐order numerical methods with increased local degrees of freedom (DOFs). Current strategies mainly follow the limiting projections derived originally for conventional finite volume methods and thus are not able to make full use of the sub‐cell information available in the local high‐order reconstructions. This paper presents a novel algorithm that introduces a nodal value‐based weighted essentially non‐oscillatory limiter for constrained interpolation profile/multi‐moment finite volume method (CIP/MM FVM) (Ii and Xiao, J. Comput. Phys., 222 (2007), 849–871) as an effort to pursue a better suited formulation to implement the limiting projection in schemes with local DOFs. The new scheme, CIP‐CSL‐WENO4 scheme, extends the CIP/MM FVM method by limiting the slope constraint in the interpolation function using the weighted essentially non‐oscillatory (WENO) reconstruction that makes use of the sub‐cell information available from the local DOFs and is built from the point values at the solution points within three neighboring cells, thus resulting a more compact WENO stencil. The proposed WENO limiter matches well the original CIP/MM FVM, which leads to a new scheme of high accuracy, algorithmic simplicity, and computational efficiency. We present the numerical results of benchmark tests for both scalar and Euler conservation laws to manifest the fourth‐order accuracy and oscillation‐suppressing property of the proposed scheme. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

19.
    
A new vortex particle‐in‐cell (PIC) method is developed for the computation of three‐dimensional unsteady, incompressible viscous flow in an unbounded domain. The method combines the advantages of the Lagrangian particle methods for convection and the use of an Eulerian grid to compute the diffusion and vortex stretching. The velocity boundary conditions used in the method are of Dirichlet‐type, and can be calculated using the vorticity field on the grid by the Biot–Savart equation. The present results for the propagation speed of the single vortex ring are in good agreement with the Saffman's model. The applications of the method to the head‐on and head‐off collisions of the two vortex rings show good agreement with the experimental and numerical literature. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

20.
    
Stability is achieved in most approximate Riemann solvers through ‘flux upwinding’, where the flux at the interface is arrived at by adding a dissipative term to the average of the left and right flux. Motivated by the existence of a collapsed interface state in the gas‐kinetic Bhatnagar–Gross–Krook (BGK) method, an alternative approach to upwinding is attempted here; an interface state is arrived at by taking an upwinded average of left and right states, and then the flux is calculated as a function of this ‘collapsed’ interface state. This so called ‘state‐upwinding’ approach gives rise to a new scheme called the linearized Riemann solver for the Euler and Navier–Stokes equations. The scheme is shown to be closely associated with the Roe scheme. It is, however, computationally less expensive and gives qualitatively comparable results over a wide range of problems. Most importantly, this scheme is found to preserve stationary contacts while not exhibiting the carbuncle phenomenon which plagues the Roe and other contact‐preserving schemes. The scheme is therefore motivated as a new starting point to analyze the origin of the carbuncle phenomenon. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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

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