Improved shock-capturing of Jameson's scheme for the Euler equations

It is known that Jameson's scheme is a pseudo-second-order-accurate scheme for solving discrete conservation laws. The scheme contains a non-linear artificial dissipative flux which is designed to capture shocks. In this paper, it is shown that the, shock-capturing of Jameson's scheme for the Euler equations can be improved by replacing the Lax-Friedrichs' type of dissipative flux with Roe's dissipative flux. This replacement is at a moderate expense of the calculation time.     

Damping efficiency of the Tchamwa–Wielgosz explicit dissipative scheme under instantaneous loading conditions

To deal with dynamic and wave propagation problems, dissipative methods are often used to reduce the effects of the spurious oscillations induced by the spatial and time discretization procedures. Among the many dissipative methods available, the Tchamwa–Wielgosz (TW) explicit scheme is particularly useful because it damps out the spurious oscillations occurring in the highest frequency domain. The theoretical study performed here shows that the TW scheme is decentered to the right, and that the damping can be attributed to a nodal displacement perturbation. The FEM study carried out using instantaneous 1-D and 3-D compression loads shows that it is useful to display the damping versus the number of time steps in order to obtain a constant damping efficiency whatever the size of element used for the regular meshing. A study on the responses obtained with irregular meshes shows that the TW scheme is only slightly sensitive to the spatial discretization procedure used. To cite this article: L. Mahéo et al., C. R. Mecanique 337 (2009).     

Construction of third-order WNND scheme and its application in complex flow

IntroductionWiththedevelopmentofaeronauticsandaerospacetechnology ,moreandmorerequirementsarearisingforCFD (computationalfluiddynamics) .Oneoftheproblemsistodevelophigherorderaccuracyschemes.Forexample ,whenapplyingLES (largeeddysimulation)orDNS(directnumericalsimulation)methodtosimulatingturbulenceproblem ,theschemesneedthirdorderaccuracyormoreinspace .Anotherquestionistheinfluenceofgrid’sscaletotopologicalstructureofflowfield .Inordertosimulatecomplicatedflowswithseparationorturbulencec…     

Viscous dissipative fluid flow past a semi‐infinite vertical plate with variable surface temperature

An analysis of the effect of viscous dissipative heat on two‐dimensional viscous incompressible fluid flow past a semi‐infinite vertical plate with variable surface temperature is carried out. The dimensionless governing equations are unsteady, two‐dimensional, coupled, and non‐linear governing equations. A most accurate, unconditionally stable and fast converging implicit finite‐difference scheme is used to solve the non‐dimensional governing equations. Velocity and temperature of the flow have been presented graphically for various parameters occurring in the problem. The local and average skin friction and Nusselt number are also shown graphically. It is observed that greater viscous dissipative heat causes a rise in the temperature. Copyright © 2007 John Wiley & Sons, Ltd.     

Dissipative control of switched nonlinear systems with sliding mode: an event-triggered control scheme

This paper introduces the event-triggered scheme to study the problems of dissipative control and dissipativity-based sliding mode control (SMC) for switched nonlinear systems. Firstly, based on a designed event-triggered scheme, a piecewise continuous state feedback controller is constructed; sufficient conditions for \(\left( Q,S,R\right) \)-\(\alpha \)-dissipativity of the resulting system are derived by the approaches of a switched Lyapunov function and an average dwell time. Secondly, these two methods are extended to the study of \(H_\infty \) performance and passive performance; the corresponding sufficient conditions are provided. Thirdly, the dissipative control problem is solved by introducing an event-triggered scheme into SMC method. Sufficient conditions for dissipativity of the resulting system are derived, and an SMC law is designed to guarantee the reachability. At the end of this paper, the correctness of the proposed methods is proved by two numerical examples.

     

A new NND difference scheme of second-order in time and space

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅＬａｘ＿ＷｅｎｄｒｏｆｆｓｃｈｅｍｅａｎｄＭａｃＣｏｒｍａｃｋｓｃｈｅｍｅａｒｅｔｗｏｗｅｌｌ＿ｋｎｏｗｎｄｉｆｆｅｒｅｎｃｅｓｃｈｅｍｅｓｗｈｉｃｈａｒｅｕｔｉｌｉｚｅｄｉｎｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｆｌｕｉｄｆｌｏｗｓｇｏｖｅｒｎｅｄｂｙＥｕｌｅｒｏｒＮａｖｉｅｒ＿Ｓｔｏｋｅｓｅｑｕａｔｉｏｎｓ.Ｔｈｅｌｉｍｉｔａｔｉｏｎｏｆｔｈｅｓｅｔｗｏｓｃｈｅｍｅｓａｒｅ:ｗｈｅｎｔｈｅｆｌｏｗｆｉｅｌｄｃｏｎｔａｉｎｓｓｈｏｃｋｓ,ｔｈｅｒｅｗｉｌｌ…     

Discretization of Asymptotically Stable Stationary Solutions of Delay Differential Equations with a Random Stationary Delay

We prove the existence of a stationary random solution to a delay random ordinary differential system, which attracts all other solutions in both pullback and forwards senses. The equation contains a one-sided dissipative Lipschitz term without delay, while the random delay appears in a globally Lipschitz one. The delay function only needs to be continuous in time. Moreover, we also prove that the split implicit Euler scheme associated to the random delay differential system generates a discrete time random dynamical system, which also possesses a stochastic stationary solution with the same attracting property, and which converges to the stationary solution of the delay random differential equation pathwise as the stepsize goes to zero.     

NND schemes and numerical simulation of axial symmetric free jet flows

Through a study on one-dimensional Navier-Stokes equations, it was found that the spurious oscillations occuring near shock waves with finite difference equations are related to the dispersion term in the corresponding modified differential equations. If the sign of dispersion coefficient is properly adjusted so that the sign changes across shock waves, the undesirable oscillations can be totally suppressed. Based on this finding, the non-oscillatory, containing no free parameters and dissipative shheme (NND scheme) is developed. This scheme is one of “TVD”. The axisymmetric free jet flows are simulated numerically using this scheme. The results obtained by the present scheme are compared with the experimental picture. It is shown that the agreement is very good, and that this scheme has advantages of high resolution for capturing shocks and contact discontinuities. Project supported by National Science Foundation of China     

A low-dissipation third-order weighted essentially nonoscillatory scheme with a new reference smoothness indicator

The classical third-order weighted essentially nonoscillatory (WENO) scheme is notoriously dissipative as it loses the optimal order of accuracy at critical points and its two-point finite difference in the smoothness indicators is unable to differentiate the critical point from the discontinuity. In recent years, modifications to the smoothness indicators and weights of the classical third-order WENO scheme have been reported to reduce numerical dissipation. This article presents a new reference smoothness indicator for constructing a low-dissipation third-order WENO scheme. The new reference smoothness indicator is a nonlinear combination of the local and global stencil smoothness indicators. The resulting WENO-Rp3 scheme with the power parameter p=1.5 achieves third-order accuracy in smooth regions including critical points and has low dissipation, but numerical results show this scheme cannot keep the ENO property near discontinuities. The recommended WENO-R3 scheme (p=1) keeps the ENO property and performs better than several recently developed third-order WENO schemes.     

Dissipative discrete breathers in a chain of Rayleigh oscillators

The dynamics of chains of interacting active particles with Rayleigh-type dissipation and coupled by a Morse potential were previously studied. In this work we introduce an on-site potential to transform this system into a chain of oscillators. We study the evolution of modes previously found in a chain of active particles as a consequence of the on-site force. Beside this, a new class of modes, dissipative discrete breathers, is found. These new modes appear due to the new time scales introduced by the on-site potential. Interaction of the dissipative discrete breathers is also investigated. We find different behaviors; for short chains dissipative discrete breathers can meet during the transitory formation reaching the stationary modes including the optical one, sufficiently long chains the dissipative discrete breathers reach its stationary state without interactions and for medium lengths, the dissipative discrete breathers interact by forming a standing wave in between.

     

A non‐iterative implicit algorithm for the solution of advection–diffusion equation on a sphere

A numerical algorithm for the solution of advection–diffusion equation on the surface of a sphere is suggested. The velocity field on a sphere is assumed to be known and non‐divergent. The discretization of advection–diffusion equation in space is carried out with the help of the finite volume method, and the Gauss theorem is applied to each grid cell. For the discretization in time, the symmetrized double‐cycle componentwise splitting method and the Crank–Nicolson scheme are used. The numerical scheme is of second order approximation in space and time, correctly describes the balance of mass of substance in the forced and dissipative discrete system and is unconditionally stable. In the absence of external forcing and dissipation, the total mass and L₂‐norm of solution of discrete system is conserved in time. The one‐dimensional periodic problems arising at splitting in the longitudinal direction are solved with Sherman–Morrison's formula and Thomas's algorithm. The one‐dimensional problems arising at splitting in the latitudinal direction are solved by the bordering method that requires a prior determination of the solution at the poles. The resulting linear systems have tridiagonal matrices and are solved by Thomas's algorithm. The suggested method is direct (without iterations) and rapid in realization. It can also be applied to linear and nonlinear diffusion problems, some elliptic problems and adjoint advection–diffusion problems on a sphere. Copyright © 2015 John Wiley & Sons, Ltd.     

A hybrid Padé ADI scheme of higher‐order for convection–diffusion problems

A high‐order Padé alternating direction implicit (ADI) scheme is proposed for solving unsteady convection–diffusion problems. The scheme employs standard high‐order Padé approximations for spatial first and second derivatives in the convection‐diffusion equation. Linear multistep (LM) methods combined with the approximate factorization introduced by Beam and Warming (J. Comput. Phys. 1976; 22 : 87–110) are applied for the time integration. The approximate factorization imposes a second‐order temporal accuracy limitation on the ADI scheme independent of the accuracy of the LM method chosen for the time integration. To achieve a higher‐order temporal accuracy, we introduce a correction term that reduces the splitting error. The resulting scheme is carried out by repeatedly solving a series of pentadiagonal linear systems producing a computationally cost effective solver. The effects of the approximate factorization and the correction term on the stability of the scheme are examined. A modified wave number analysis is performed to examine the dispersive and dissipative properties of the scheme. In contrast to the HOC‐based schemes in which the phase and amplitude characteristics of a solution are altered by the variation of cell Reynolds number, the present scheme retains the characteristics of the modified wave numbers for spatial derivatives regardless of the magnitude of cell Reynolds number. The superiority of the proposed scheme compared with other high‐order ADI schemes for solving unsteady convection‐diffusion problems is discussed. A comparison of different time discretizations based on LM methods is given. Copyright © 2009 John Wiley & Sons, Ltd.     

An analysis of the stability of the least square finite difference scheme and its shock capturing ability in inviscid flows

A meshless method – The Least Square Finite Difference scheme (LSFD) with diffusion is analyzed and applied to inviscid flows. The scheme is made second-order by using a modified difference in the formulation of LSFD. Several numerical experiments, namely the Sod shock tube and the shallow water problems, are carried out and, in the limelight of the results obtained, the ability of the scheme to resolve shock wave, rarefaction wave, and contact discontinuity is discussed. The conditional stability of the LSFD scheme is established. The LSFD uses weights to diagonalize the least square matrix resulting in the spatial discretization in order to gain computational time. We prove that there exists a unique weight for the resulting optimization problem. The weighted version of LSFD is used to solve the isentropic vortex problem numerically and the results are used to discuss the dissipative nature of the scheme. Five configurations of the two-dimensional Riemann problems are used in our numerical experiments. The capability of the scheme to capture the complex interaction of multiple planar waves is discussed in the limelight of the results on the Riemann problems. The result of the shock reflection problem shows that the scheme is minimally dissipative and leads to sharp and well-resolved shocks.     

Stabilized seventh-order dissipative compact scheme using simultaneous approximation terms

To ensure time stability of a seventh-order dissipative compact finite difference scheme, fourth-order boundary closures are used near domain boundaries previously.However, this would reduce the global convergence rate to fifth-order only. In this paper,we elevate the boundary closures to sixth-order to achieve seventh-order global accuracy.To keep the improved scheme time stable, the simultaneous approximation terms(SATs)are used to impose boundary conditions weakly. Eigenvalue analysis shows that the improved scheme is time stable. Numerical experiments for linear advection equations and one-dimensional Euler equations are implemented to validate the new scheme.     

Heat transfer in a flow past a semi-in-finite plate with oscillating temperature

The response of laminar boundary layer flow past a semi-infinite flat plate to harmonic oscillations in the plate temperature in the form of a travelling wave convected in the direction of the free-stream has been studied here. Series solutions in terms of the small amplitude and the small oscillations to the non-linear system have been derived and the resulting nonlinear ordinary equations due to usual similarity transformations are solved numerically. The function affecting the temperature is shown on a graph. Due to greater viscous dissipative heat the function K ₁, increases and it decreases with increasing Prandtl number. Also the time averaged heat flux function K ₁(0) increases with Prandtl number and decreases due to greater viscous dissipative heat.     

Stability of a Crank–Nicolson pressure correction scheme based on staggered discretizations

In the context of LES of turbulent flows, the control of kinetic energy seems to be an essential requirement for a numerical scheme. Designing such an algorithm, that is, as less dissipative as possible while being simple, for the resolution of variable density Navier–Stokes equations is the aim of the present work. The developed numerical scheme, based on a pressure correction technique, uses a Crank–Nicolson time discretization and a staggered space discretization relying on the Rannacher–Turek finite element. For the inertia term in the momentum balance equation, we propose a finite volume discretization, for which we derive a discrete analogue of the continuous kinetic energy local conservation identity. Contrary to what was obtained for the backward Euler discretization, the dissipation defect term associated to the Crank–Nicolson scheme is second order in time. This behavior is evidenced by numerical simulations. Copyright © 2013 John Wiley & Sons, Ltd.     

The multi‐stage centred‐scheme approach applied to a drift‐flux two‐phase flow model

For two‐phase flow models, upwind schemes are most often difficult do derive, and expensive to use. Centred schemes, on the other hand, are simple, but more dissipative. The recently proposed multi‐stage (MUSTA ) method is aimed at coming close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. So far, the MUSTA approach has been shown to work well for the Euler equations of inviscid, compressible single‐phase flow. In this work, we explore the MUSTA scheme for a more complex system of equations: the drift‐flux model, which describes one‐dimensional two‐phase flow where the motions of the phases are strongly coupled. As the number of stages is increased, the results of the MUSTA scheme approach those of the Roe method. The good results of the MUSTA scheme are dependent on the use of a large‐enough local grid. Hence, the main benefit of the MUSTA scheme is its simplicity, rather than CPU ‐time savings. Copyright © 2006 John Wiley & Sons, Ltd.     

Baroclinic stability for a family of two‐level,semi‐implicit numerical methods for the 3D shallow water equations

The baroclinic stability of a family of two time‐level, semi‐implicit schemes for the 3D hydrostatic, Boussinesq Navier–Stokes equations (i.e. the shallow water equations), which originate from the TRIM model of Casulli and Cheng (Int. J. Numer. Methods Fluids 1992; 15 :629–648), is examined in a simple 2D horizontal–vertical domain. It is demonstrated that existing mass‐conservative low‐dissipation semi‐implicit methods, which are unconditionally stable in the inviscid limit for barotropic flows, are unstable in the same limit for baroclinic flows. Such methods can be made baroclinically stable when the integrated continuity equation is discretized with a barotropically dissipative backwards Euler scheme. A general family of two‐step predictor‐corrector schemes is proposed that have better theoretical characteristics than existing single‐step schemes. Copyright © 2006 John Wiley & Sons, Ltd.     

结构动力响应数值算法耗散和超调特性设计

对结构动力响应数值计算问题提出引入多个自由参数来获得所希望的算法特性. 多参数的一个明显的好处就是在算法设计上有更大的自由空间. 利用这些自由参数获得了两个新的无条件稳定、有二阶精度的、有好的耗散和没有超调的单步时间直接积分算法. 在存在阻尼情况下基于有限差分分析理论证明了新算法的这些特性. 其中一个有高频渐进消去特性,且在有阻尼情况下与Houbolt方法相比对高频有更强的耗散. 另一个在低频极限无耗散,高频耗散可以用一自由参数控制. 超调分析结果显示两个新算法都不显示超调,而HHT方法不仅有速度超调,还有位移超调. 最后使用一些算例并通过与传统方法的比较数值地验证了理论分析结果.     

To the boundary and back—a numerical study

This study identifies the key parameters upon which energy absorption at artificial boundaries depends. A thorough numerical study is presented, of typical reflections from open computational boundaries, for problems governed by hyperbolic systems of equations. The emphasis is on systems, where it is often the combination of all boundary procedures that determine the quality of boundary treatment. We study dissipative numerical models which have so far not been analysed to the same extent as non-dissipative models and employ a Law-Wendroff-type scheme as a prototype. While it is widely accepted that dissipative models tend to give fewer problems than non-dissipative ones, we show a variety of cases where substantial reflections do occur even in ID and quasi-ID set-ups, where theory predicts best results. This can partly be explained by the vanishing of dissipation in the far field. Group velocity analysis, justifiable on the grounds of weak dissipation, predicts a pathological behaviour which is confirmed by numerical experiments. We demonstrate strong focusing of asymptotic errors generated at the artificial boundary. Internal reflections due to slowly expanding grids are shown for non-linear systems. The need for high-frequency boundary conditions naturally arises and combined low-high-frequency boundary recipes following Higdon, Vichnevetsky and Pariser are adapted to systems and tested. Partial cures are also discussed, mainly in terms of pointing out their theoretically limited potential.