Engineering Papers⌕ Search

Engineering topics

Shu, Chi-Wang

Publications and source records attributed to Shu, Chi-Wang.

54 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

On the Gibbs phenomenon 1: Recovering exponential accuracy from the Fourier partial sum of a non-periodic analytic function

It is well known that the Fourier series of an analytic or periodic function, truncated after 2N+1 terms, converges exponentially with N, even in the maximum norm, although the function is still analytic. This is known as the Gibbs phenomenon. Here, we show that the first 2N+1 Fourier coefficients contain enough information about the function, so that an exponentially convergent approximation (in the maximum norm) can be constructed.

Gottlieb, David↗
NASA NTRS1992

Effective equations and the inverse cascade theory for Kolmogorov flows

We study the two dimensional Kolmogorov flows in the limit as the forcing frequency goes to infinity. Direct numerical simulation indicates that the low frequency energy spectrum evolves to a universal kappa (exp -4) decay law. We derive effective equations governing the behavior of the large scale flow quantities. We then present numerical evidence that with smooth initial data, the solution to the effective equation develops a kappa (exp -4) type singularity at a finite time. This gives a convenient explanation for the kappa (exp -4) decay law exhibited by the original Kolmogorov flows.

Weinan, E.↗
NASA NTRS1991

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗
NASA NTRS1991

High-order ENO schemes applied to two- and three-dimensional compressible flow

High order essentially non-oscillatory (ENO) finite difference schemes are applied to the 2-D and 3-D compressible Euler and Navier-Stokes equations. Practical issues, such as vectorization, efficiency of coding, cost comparison with other numerical methods, and accuracy degeneracy effects, are discussed. Numerical examples are provided which are representative of computational problems of current interest in transition and turbulence physics. These require both nonoscillatory shock capturing and high resolution for detailed structures in the smooth regions and demonstrate the advantage of ENO schemes.

Shu, Chi-Wang↗
NASA NTRS1991

Uniform high order spectral methods for one and two dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗
NASA NTRS1991

The P1-RKDG method for two-dimensional Euler equations of gas dynamics

A class of nonlinearly stable Runge-Kutta local projection discontinuous Galerkin (RKDG) finite element methods for conservation laws is investigated. Two dimensional Euler equations for gas dynamics are solved using P1 elements. The generalization of the local projections, which for scalar nonlinear conservation laws was designed to satisfy a local maximum principle, to systems of conservation laws such as the Euler equations of gas dynamics using local characteristic decompositions is discussed. Numerical examples include the standard regular shock reflection problem, the forward facing step problem, and the double Mach reflection problem. These preliminary numerical examples are chosen to show the capacity of the approach to obtain nonlinearly stable results comparable with the modern nonoscillatory finite difference methods.

Cockburn, Bernardo↗
NASA NTRS1990

Numerical experiments on the accuracy of ENO and modified ENO schemes

Further numerical experiments are made assessing an accuracy degeneracy phenomena. A modified essentially non-oscillatory (ENO) scheme is proposed, which recovers the correct order of accuracy for all the test problems with smooth initial conditions and gives comparable results with the original ENO schemes for discontinuous problems.

Shu, Chi-Wang↗
NASA NTRS1990

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗
NASA NTRS1990

Numerical experiments on the accuracy of ENO and modified ENO schemes

Numerical experiments have been performed using different ENO schemes and different time discretizations, in order to assess accuracy-degeneracy phenomena of the type described by Rogerson and Meiburg (1990). A modified ENO scheme recovers the correct order of accuracy for all the test problems with smooth initial conditions and gives results comparable to the original ENO schemes for discontinuous problems. It is concluded that ENO schemes may lose the full high-order accuracy predicted by local truncation-error analysis for scalar linear conservation laws with smooth initial conditions. A modified ENO scheme can overcome this accuracy-degeneracy problem for the test problems without increasing the computational cost.

Shu, Chi-Wang↗
NASA NTRS1989

Efficient implementation of essentially non-oscillatory shock-capturing schemes. II

The present elaboration of essentially nonoscillatory (ENO) shock-capturing schemes proceeds with a novel, simplified expression for the ENO construction procedure having its basis in numerical fluxes rather than cell-averages. The ENO-local Lax-Friedrichs and ENO-Roe schemes introduced yield sharper shock transitions and greater overall accuracy than previous ENO-scheme implementations. Numerical fluxes and TVD Runge-Kutta time discretizations are used to apply Harten's subcell-resolution concept (1986) and Yang's artificial compression method to the current ENO schemes, in the interest of sharper contact discontinuities.

Shu, Chi-Wang↗
NASA NTRS1989

Essentially nonoscillatory spectral Fourier methods for shock wave calculations

An essentially nonoscillatory spectral Fourier method for the solution of hyperbolic partial differential equations is presented. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that essentially no oscillations develop in the solution.

Cai, Wei↗
NASA NTRS1989

Recent progress on essentially non-oscillatory shock capturing schemes

An account is given of the construction of efficient implementations of 'essentially nonoscillatory' (ENO) schemes that approximate systems of hyperbolic conservation laws. ENO schemes use a local adaptive stencil to automatically obtain information from regions of smoothness when the solution develops discontinuities. Approximations employing ENOs can thereby obtain uniformly high accuracy to the very onset of discontinuities, while retaining a sharp and essentially nonoscillatory shock transition. For ease of implementation, ENO schemes applying the adaptive stencil concept to the numerical fluxes and employing a TVD Runge-Kutta-type time discretization are constructed.

Osher, Stanley↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock-capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably - especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗
NASA NTRS1988

Non-oscillatory spectral Fourier methods for shock wave calculations

A non-oscillatory spectral Fourier method is presented for the solution of hyperbolic partial differential equations. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that no oscillations develop in the solution. Also, the accuracy of the spectral solution of the inviscid Burgers equation is shown to be higher than a fixed order.

Cai, Wei↗
NASA NTRS1988

Efficient implementation of essentially non-oscillatory shock capturing schemes, 2

Earlier work on the efficient implementation of ENO (essentially non-oscillatory) shock capturing schemes is continued. A new simplified expression is provided for the ENO construction procedure based again on numerical fluxes rather than cell averages. Also considered are two improvements which are labeled ENO-LLF (local Lax-Friedrichs) and ENO-Roe, which yield sharper shock transitions, improved overall efficiency, and lower computational cost than previous implementation of the ENO schemes. Two methods of sharpening contact discontinuities, i.e., the subcell resolution idea of Harten and the artificial compression idea of Yang, which those authors used originally in the cell-average framework, are supplied to the current ENO schemes using numerical fluxes and TVD Runge-Kutta time discretizations. The implementation for nonlinear systems and multi-dimensions is given. Finally, many numerical examples, including a compressible shock turbulence interaction flow calculation, are given.

Shu, Chi-Wang↗
NASA NTRS1988

Essentially non-oscillatory shock capturing methods applied to turbulence amplification in shock wave calculations

ENO (essentially non-oscillatory) schemes can provide uniformly high order accuracy right up to discontinuities while keeping sharp, essentially non-oscillatory shock transitions. Recently, an efficient implementation of ENO schemes was obtained based on fluxes and TVD Runge-Kutta time discretizations. The resulting code is very simple to program for multi-dimensions. ENO schemes are especially suitable for computing problems with both discontinuities and fine structures in smooth regions, such as shock interaction with turbulence, for which results for 1-D and 2-D Euler equations are presented. Much better resolution is observed by using third order ENO schemes than by using second order TVD schemes for such problems.

Osher, Stanley↗
NASA NTRS1987

Efficient implementation of essentially non-oscillatory shock capturing schemes

In the computation of discontinuous solutions of hyperbolic conservation laws, TVD (total-variation-diminishing), TVB (total-variation-bounded) and the recently developed ENO (essentially non-oscillatory) schemes have proven to be very useful. In this paper two improvements are discussed: a simple TVD Runge-Kutta type time discretization, and an ENO construction procedure based on fluxes rather than on cell averages. These improvements simplify the schemes considerably -- especially for multi-dimensional problems or problems with forcing terms. Preliminary numerical results are also given.

Shu, Chi-Wang↗