Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Oscillation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Spurious Numerical Oscillations in Numerical Simulation of Supersonic Flows Using Shock Capturing Schemes

The numerical simulation of transitional and turbulent processes in hypersonic boundary layers often involves a physical process of a shock-disturbance wave interaction in complex two-dimensional and three-dimensional flow fields. For such simulations, it is required that there be a high order of accuracy in capturing both the shock and the small disturbances. The purpose of this paper is to evaluate the viability of using high order shock capturing schemes to track small disturbances in a multi-dimensional steady hypersonic flow. The numerical methods that are to be studied are the Total Variation Diminishing (TVD) scheme, and Essentially Non-Oscillatory (ENO) scheme. This paper shows that the presence of numerical oscillations in the flow field solution may drastically hinder any attempt at tracking the propagation of any physical disturbances. It has been found that the numerical oscillations that exist for shock capturing methods may be significant enough to pollute a flow field containing small physical disturbances. The effects of the refinement of the grid do not reduce the oscillations, but rather they decrease the wavelength of the oscillations. It is shown that by aligning the shock with the grid, the amplitude of these spurious oscillations may be greatly reduced.

Lee, Theodore K.↗
NASA NTRS1972

Solar seeing and the spatial properties of the five-minute oscillations.

Numerical simulation of observations of the spatial properties of the five-minute oscillations, assuming that the oscillations are internal gravity waves excited by granular convection according to the theory of Thomas et al. (1971). The simulation includes the effects of seeing and finite aperture. The details of the simulation are chosen to model the observational method of Frazier (1968). The results show that the peak in the observed power spectrum of the oscillations can occur at a wavelength considerably longer than the true wavelength of the oscillations. In particular, the peak in Frazier's observed power spectra at wavelength of about 5000 km is consistent with the considerably shorter true wavelength of about 1500 km predicted by the gravity-wave theory.

Thomas, J. H.↗
NASA NTRS1977

Improved numerical control of oscillator frequency

Numerically-controlled oscillator is key element in phase-locked loop that generates exact frequency and phase of noise-corrupted input signal. Device is applicable to communication and tracking equipment.

Cellier, A.↗
NASA NTRS1980

Digital numerically controlled oscillator

The frequency and phase of an output signal from an oscillator circuit are controlled with accuracy by a digital input word. Positive and negative alterations in output frequency are both provided for by translating all values of input words so that they are positive. The oscillator reference frequency is corrected only in one direction, by adding phase to the output frequency of the oscillator. The input control word is translated to a single algebraic sign and the digital 1 is added thereto. The translated input control word is then accumulated. A reference clock signal having a frequency at an integer multiple of the desired frequency of the output signal is generated. The accumulated control word is then compared with a threshold level. The output signal is adjusted in a single direction by dividing the frequency of the reference clock signal by a first integer or by an integer different from the first integer.

Cellier, A.↗
NASA NTRSDate not supplied

Computationally Efficient Frequency Domain Method for Dynamic Control Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a computationally efficient method for dynamic control derivative estimation technique via control surface oscillation numerical experiments. Unsteady RANS CFD simulations in FUN3D are performed to simulate the control surface oscillations with a prescribed truncated square wave containing sufficient frequencies of interest. The truncated square wave oscillation offers the computational efficiency which reduces the computational cost by almost an order of magnitude compared to a sine wave oscillation. The nonlinear effect of large control surface oscillation amplitudes creates a spillover effect whereby the frequency response at the same input frequencies contains not only the linear aerodynamic response but also nonlinear aerodynamic response. A correction procedure is developed to remove the spillover effect from the linear aerodynamic response. A frequency-domain regression is performed to estimate the dynamic control derivatives after the correction. The results generally agree with the previous results obtained from the sine wave oscillation.

Stability and Control↗
NASA NTRS1990

High order filtering methods for approximating hyberbolic systems of conservation laws

In the computation of discontinuous solutions of hyperbolic systems of conservation laws, the recently developed essentially non-oscillatory (ENO) schemes appear to be very useful. However, they are computationally costly compared to simple central difference methods. A filtering method which is developed uses simple central differencing of arbitrarily high order accuracy, except when a novel local test indicates the development of spurious oscillations. At these points, the full ENO apparatus is used, maintaining the high order of accuracy, but removing spurious oscillations. Numerical results indicate the success of the method. High order of accuracy was obtained in regions of smooth flow without spurious oscillations for a wide range of problems and a significant speed up of generally a factor of almost three over the full ENO method.

Lafon, F.↗
NASA NTRS1993

Investigation of three-dimensional separation at wing/body junctions in supersonic flows using TVD MacCormack's scheme

A Total Variation Diminishing (TVD) MacCormack scheme is used to study the three-dimensional separation at wing/body junctions in supersonic flows. The new scheme is robust, resolves discontinuities well with-out any numerical oscillations present in the flow field. Moreover, this new scheme does not contain any problem dependent parameters to be adjusted and many production codes employing MacCormack algorithm can be easily updated to good effect. Numerical simulation carried out in laminar supersonic junction flows using the new scheme yields improved prediction for the separation location and the axial velocity profiles in the separated flow region.

Lakshmanan, B.↗
NASA NTRS1992

Multi-dimensional high order essentially non-oscillatory finite difference methods in generalized coordinates

The nonlinear stability of compact schemes for shock calculations is investigated. In recent years compact schemes were used in various numerical simulations including direct numerical simulation of turbulence. However to apply them to problems containing shocks, one has to resolve the problem of spurious numerical oscillation and nonlinear instability. A framework to apply nonlinear limiting to a local mean is introduced. The resulting scheme can be proven total variation (1D) or maximum norm (multi D) stable and produces nice numerical results in the test cases. The result is summarized in the preprint entitled 'Nonlinearly Stable Compact Schemes for Shock Calculations', which was submitted to SIAM Journal on Numerical Analysis. Research was continued on issues related to two and three dimensional essentially non-oscillatory (ENO) schemes. The main research topics include: parallel implementation of ENO schemes on Connection Machines; boundary conditions; shock interaction with hydrogen bubbles, a preparation for the full combustion simulation; and direct numerical simulation of compressible sheared turbulence.

Shu, Chi-Wang↗
NASA NTRS1992

Nonlinearly stable compact schemes for shock calculations

The applications of high-order, compact finite difference methods in shock calculations are discussed. The main concern is to define a local mean which will serve as a reference for introducing a local nonlinear limiting to control spurious numerical oscillations while maintaining the formal accuracy of the scheme. For scalar conservation laws, the resulting schemes can be proven total-variation stable in one space dimension and maximum-norm stable in multiple space dimensions. Numerical examples are shown to verify accuracy and stability of such schemes for problems containing shocks. These ideas can also be applied to other implicit schemes such as the continuous Galerkin finite element methods.

Cockburn, Bernardo↗
NASA NTRS1988

A numerical method for acoustic oscillations in tubes

A numerical method to obtain the neutral curve for the onset of acoustic oscillations in a helium-filled tube is described. Such oscillations can cause a serious heat loss in the plumbing associated with liquid helium dewars. The problem is modelled by a second-order, ordinary differential eigenvalue problem for the pressure perturbation. The numerical method to find the eigenvalues and track the resulting points along the neutral curve is tailored to this problem. The results show that a tube with a uniform temperature gradient along it is much more stable than one where the temperature suddenly jumps from the cold to the hot value in the middle of the tube.

Gary, John M.↗
NASA NTRS1987

A new algorithm for the Navier-Stokes equations applied to transonic flows over wings

A new algorithm has been developed for the Navier-Stokes equations. For the convective terms, this algorithm employs flux vector splitting together with a locally rotated coordinate system that is aligned with the streamwise direction. The algorithm has been implemented into a Navier-Stokes code, the NASA Ames Research Center's TNS code, and several cases of steady flow have been calculated. The results show improvements over the original method in the code. First, in a case of separated flow, there is better agreement with the experimental results and also there is an absence of numerical oscillations that occur with the original method. Second, in a case with a strong supersonic to subsonic shock wave, comparisons show an improvement in the resolution of the shock wave and the reexpansion singularity.

Goorjian, Peter M.↗
NASA NTRS1989

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗
NASA NTRS1990

Some Aspects of Essentially Nonoscillatory (ENO) Formulations for the Euler Equations, Part 3

An essentially nonoscillatory (ENO) formulation is described for hyperbolic systems of conservation laws. ENO approaches are based on smart interpolation to avoid spurious numerical oscillations. ENO schemes are a superset of Total Variation Diminishing (TVD) schemes. In the recent past, TVD formulations were used to construct shock capturing finite difference methods. At extremum points of the solution, TVD schemes automatically reduce to being first-order accurate discretizations locally, while away from extrema they can be constructed to be of higher order accuracy. The new framework helps construct essentially non-oscillatory finite difference methods without recourse to local reductions of accuracy to first order. Thus arbitrarily high orders of accuracy can be obtained. The basic general ideas of the new approach can be specialized in several ways and one specific implementation is described based on: (1) the integral form of the conservation laws; (2) reconstruction based on the primitive functions; (3) extension to multiple dimensions in a tensor product fashion; and (4) Runge-Kutta time integration. The resulting method is fourth-order accurate in time and space and is applicable to uniform Cartesian grids. The construction of such schemes for scalar equations and systems in one and two space dimensions is described along with several examples which illustrate interesting aspects of the new approach.

Chakravarthy, Sukumar R.↗
NASA NTRS1991

Artificial Dissipation In Computations Of Hypersonic Flows

Four artificial-dissipation mathematical models introduced to suppress spurious numerical oscillations in finite-difference computations of hypersonic external flows. Models tested for their effects in capturing details of shocks in hypersonic flows about bodies of various shapes.

Yoon, S.↗
NASA NTRS1991

Entropy-Based Approach To Nonlinear Stability

NASA technical memorandum suggests schemes for numerical solution of differential equations of flow made more accurate and robust by invoking second law of thermodynamics. Proposes instead of using artificial viscosity to suppress such unphysical solutions as spurious numerical oscillations and nonlinear instabilities, one should formulate equations so that rate of production of entropy within each cell of computational grid be nonnegative, as required by second law.

Merriam, Marshal L.↗
NASA NTRS1992

Finite element computation of incompressible flows

Among the solution techniques presented for FEM computation of incompressible flows are stabilized formulations designed to prevent the numerical oscillations that may be generated by (1) the presence of dominant advection terms in the governing equations, or (2) the failure to use an acceptable combination of interpolation functions to represent velocity and pressure fields. It is shown that the space-time formulation, in conjunction with the Galerkin-least squares methods, can be used to compute steady incompressible flows involving moving boundaries and interfaces; it is recommended that iteration methods be used to reduce the associated computational costs.

Tezduyar, T. E.↗
NASA NTRS1994

New Developments in the Method of Space-Time Conservation Element and Solution Element-Applications to Two-Dimensional Time-Marching Problems

A new numerical discretization method for solving conservation laws is being developed. This new approach differs substantially in both concept and methodology from the well-established methods, i.e., finite difference, finite volume, finite element, and spectral methods. It is motivated by several important physical/numerical considerations and designed to avoid several key limitations of the above traditional methods. As a result of the above considerations, a set of key principles for the design of numerical schemes was put forth in a previous report. These principles were used to construct several numerical schemes that model a 1-D time-dependent convection-diffusion equation. These schemes were then extended to solve the time-dependent Euler and Navier-Stokes equations of a perfect gas. It was shown that the above schemes compared favorably with the traditional schemes in simplicity, generality, and accuracy. In this report, the 2-D versions of the above schemes, except the Navier-Stokes solver, are constructed using the same set of design principles. Their constructions are simplified greatly by the use of a nontraditional space-time mesh. Its use results in the simplest stencil possible, i.e., a tetrahedron in a 3-D space-time with a vertex at the upper time level and other three at the lower time level. Because of the similarity in their design, each of the present 2-D solvers virtually shares with its 1-D counterpart the same fundamental characteristics. Moreover, it is shown that the present Euler solver is capable of generating highly accurate solutions for a famous 2-D shock reflection problem. Specifically, both the incident and the reflected shocks can be resolved by a single data point without the presence of numerical oscillations near the discontinuity.

Chang, Sin-Chung↗