Search NASA⌕ Search

SEARCH · Search NASA

Results for “Burgers equations”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A new explicit method for the numerical solution of parabolic differential equations

A new method is derived for solving parabolic partial differential equations arising in transient heat conduction or in boundary-layer flows. The method is based on a combination of the modified differential quadrature (MDQ) method with the rational Runge-Kutta time-integration scheme. It is fully explicit, requires no matrix inversion, and is stable for any time-step for the heat equations. Burgers equation and the one- and two-dimensional heat equations are solved to demonstrate the accuracy and efficiency of the proposed algorithm. The present method is found to be very accurate and efficient when results are compared with analytic solutions.

Satofuka, N.↗
NASA NTRS1994

Fourth order difference methods for hyperbolic IBVP's

Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.

Gustafsson, Bertil↗
NASA NTRS1976

Analytic numerical solutions for shock waves

Study of weak solutions of simple wave equation, inviscid Burgers equation, and Euler equations has resulted in technique for accurate prediction of shock waves occurring in inviscid supersonic flows.

Maccormack, R. W.↗
NASA NTRS1979

Application of the method of integral relations to unsteady fluid flow problems with shocks

A mixed method using both the method of integral relations and the finite difference technique is developed for the solution of unsteady flow problems. The integral relations method is based upon a chosen interpolating function dependent only on the time domain. The resulting local semidiscrete finite element equations obtained are assembled into a global form. The spatial derivatives at the nodes are replaced by finite difference operators and the discretized nonlinear algebraic system is solved by an iterative scheme. Solutions are obtained for the one-dimensional gasdynamics equation, the one-dimensional wave equation and Burger's model of turbulence. Agreement with other numerical and analytical solutions is excellent for the gasdynamics problem and satisfactory for the Burger equation in cases of small viscous effects.

Wadia, A. R.↗
NASA NTRS2014

A Walsh Function Module Users' Manual

The solution of partial differential equations (PDEs) with Walsh functions offers new opportunities to simulate many challenging problems in mathematical physics. The approach was developed to better simulate hypersonic flows with shocks on unstructured grids. It is unique in that integrals and derivatives are computed using simple matrix multiplication of series representations of functions without the need for divided differences. The product of any two Walsh functions is another Walsh function - a feature that radically changes an algorithm for solving PDEs. A FORTRAN module for supporting Walsh function simulations is documented. A FORTRAN code is also documented with options for solving time-dependent problems: an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the usage of the Walsh function module including such features as operator overloading, Fast Walsh Transforms in multi-dimensions, and a Fast Walsh reciprocal.

Gnoffo, Peter A.↗
NASA NTRS2015

Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series

Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.

Gnoffo, Peter A.↗
NASA NTRS1994

Steepened structures in equatorial spread F. 2: Theory

This paper investigates the properties of a one-dimensional fluid model of plasma convection in the equatorial F region ionosphere. The model equations are similar in form to Burgers equation except for additional higher-order spatial derivatives. Like Burgers equation, solution to the model have the form of propagating, shocklike structures. Numerical simulations of the model closely resemble the steepened structures observed by sounding rocket plasma density probes within equatorial spread F. Simulated denstiy power spectra, like the spectra computed from in situ data, seem to possess power law forms with a break at wavelengths of about 100 m. The precise wavenumber of the spectral break is determined by the ambipolar diffusion coefficient. The model predicts that electric field fluctuations perpendicular to the direction of plasma steepening should be proportional to the plasma density fluctuations. Electric field fluctuations parallel to the steepening will be due primarily to the ambipolar field and have a Boltzmann relationship with density (square of the absolute value of delta E) approximately equal to (K(exp 2))(square of the absolute value of (delta n/n)). At wavelengths less than about 300 m, the ambipolar field should be the dominant component of the total field intensity.

Hysell, D. L↗
NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier- Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle

Miller, Steven A. E.↗
NASA NTRS2015

Toward a Nonlinear Acoustic Analogy: Turbulence as a Source of Sound and Nonlinear Propagation

An acoustic analogy is proposed that directly includes nonlinear propagation effects. We examine the Lighthill acoustic analogy and replace the Green's function of the wave equation with numerical solutions of the generalized Burgers' equation. This is justified mathematically by using similar arguments that are the basis of the solution of the Lighthill acoustic analogy. This approach is superior to alternatives because propagation is accounted for directly from the source to the far-field observer instead of from an arbitrary intermediate point. Validation of a numerical solver for the generalized Burgers' equation is performed by comparing solutions with the Blackstock bridging function and measurement data. Most importantly, the mathematical relationship between the Navier-Stokes equations, the acoustic analogy that describes the source, and canonical nonlinear propagation equations is shown. Example predictions are presented for nonlinear propagation of jet mixing noise at the sideline angle.

Miller, Steven A. E.↗
NASA NTRS1989

Calculation of shocked flows by mathematical programming

A framework for using mathematical programming to solve Burgers' equation is presented. The steady-state inviscid Burgers' equation with given boundary conditions is considered, and the physically relevant solution is discretized using a four-point difference scheme for the viscous term and a two-point difference scheme for the inviscid term. The framework is then used to solve the Euler equations for quasi-one-dimensional flows on grids with variable spacing.

Lavery, John E.↗
NASA NTRS1992

A multiple scales approach to sound generation by vibrating bodies

The problem of determining the acoustic field in an inviscid, isentropic fluid generated by a solid body whose surface executes prescribed vibrations is formulated and solved as a multiple scales perturbation problem, using the Mach number M based on the maximum surface velocity as the perturbation parameter. Following the idea of multiple scales, new 'slow' spacial scales are introduced, which are defined as the usual physical spacial scale multiplied by powers of M. The governing nonlinear differential equations lead to a sequence of linear problems for the perturbation coefficient functions. However, it is shown that the higher order perturbation functions obtained in this manner will dominate the lower order solutions unless their dependence on the slow spacial scales is chosen in a certain manner. In particular, it is shown that the perturbation functions must satisfy an equation similar to Burgers' equation, with a slow spacial scale playing the role of the time-like variable. The method is illustrated by a simple one-dimenstional example, as well as by three different cases of a vibrating sphere. The results are compared with solutions obtained by purely numerical methods and some insights provided by the perturbation approach are discussed.

Geer, James F.↗
NASA NTRS1995

On the modulational instability of large amplitude waves in supersonic boundary layers

The evolution of large amplitude Tollmien-Schlichting waves in a supersonic boundary layer is investigated. Disturbances which have their wavenumber and frequency slowly varying in time and space are described using a phase equation type of approach. Unlike the incompressible case we find that the initial bifurcation to a finite amplitude Tollmien-Schlichting wave is subcritical for most Mach numbers. In fact the bifurcation is only supercritical for a small range of Mach numbers and even then for only a finite range of wave propagation angles. The modulational instability of large amplitude wavetrains is considered and is shown to be governed by an equation similar to Burgers equation but with the viscous term replaced by a fractional derivative. A numerical investigation of the solution of this equation is described. It is shown that uniform wavetrains are unstable.

Hall, Philip↗
NASA NTRS1993

Progressive wave equations and algorithms for sonic boom propagation

Raspet, Bass, and others have carried out a number of studies of sonic boom propagation using algorithms that can be traced back to a 1973 doctoral thesis written by Pestorius at the University of Texas. The algorithm was originally set forth without a formal derivation from a governing set of partial differential equations. Perusal of Pestorius's report and of subsequent reports from the Texas nonlinear accoustics group suggests that such a derivation has not yet appeared in the archival literature. The author and his colleagues, on the other hand, have been working with an explicit set of approximate partial differential equations analogous to Burgers' equation. The present paper traces through the pertinent statements in the Pestorius reports and shows that the derived algorithm, at least in the limit of sufficiently short steps, is equivalent to the numerical solution of a similar (and nearly equivalent) explicit set of partial differential equations. The tying-together of the two approaches is important because it is possible that some of the basic ideas embodied in the Pestorius algorithm may be incorporated in other contexts where more nearly general models of sonic boom propagation through realistic atmospheres are to be solved numerically.

Pierce, Allan D.↗
NASA NTRS1990

Weakly multi-dimensional cosmic-ray-modified MHD shocks

The multi-dimensional structure of weak, energetic-particle-modified shocks is investigated by means of appropriate perturbation techniques. The time-dependent shock-structure equation is found to be a generalized form of the well-known one-dimensional Burgers equation, whose steady state, in the absence of cosmic rays, is shown to be related to an equation modeling steady transonic flow in several dimensions. The time-dependent (1 + 2)- and (1 + 3)-dimensional Burgers equations are integrated exactly by means of Hirota's technique for one-shock solutions. On the basis of the exact solutions, a discussion relating the various length scales associated with the shock is presented.

Zank, G. P.↗
NASA NTRS2023

Multiple Scattering Theory for Modeling Sonic Booms in Atmospheric Turbulence

Sonic booms are often modeled using Burgers equations accounting for dominant propagation effects. Scattering effects of turbulence, however, have not been incorporated into such equations, although these effects are ubiquitous in measured sonic booms. This paper formulates the mean scattering effects, including backscattering, using multiple scattering theory and ensemble averaging. It obtains an acousto-turbulence interaction energy representing the interaction of an acoustic wave field with a turbulence field. The interaction energy gives rise to a scattering wavenumber, a complex-valued correction to the free-space wavenumber. The scattering wavenumber leads to a dispersion and attenuation of the mean waveform due to backscattering, and can be obtained from a derived exact solution. An existing Burgers equation is extended to include the scattering effects of turbulence as an additional linear term. Numerical simulations of an N-wave and a low boom show that the mean scattering effects lead to shock thickening in the mean waveforms as well as reduce the peak amplitudes and total energy contents. The energy reduction is less severe in the low boom than in the N-wave and is dependent on the variance, correlation length, and thickness of the turbulence field.

sonic booms↗
NASA NTRS1996

The Evolution of Finite Amplitude Wavetrains in Plane Channel Flow

We consider a viscous incompressible fluid flow driven between two parallel plates by a constant pressure gradient. The flow is at a finite Reynolds number, with an 0(l) disturbance in the form of a traveling wave. A phase equation approach is used to discuss the evolution of slowly varying fully nonlinear two dimensional wavetrains. We consider uniform wavetrains in detail, showing that the development of a wavenumber perturbation is governed by Burgers equation in most cases. The wavenumber perturbation theory, constructed using the phase equation approach for a uniform wavetrain, is shown to be distinct from an amplitude perturbation expansion about the periodic flow. In fact we show that the amplitude equation contains only linear terms and is simply the heat equation. We review, briefly, the well known dynamics of Burgers equation, which imply that both shock structures and finite time singularities of the wavenumber perturbation can occur with respect to the slow scales. Numerical computations have been performed to identify areas of the (wavenumber, Reynolds number, energy) neutral surface for which each of these possibilities can occur. We note that the evolution equations will breakdown under certain circumstances, in particular for a weakly nonlinear secondary flow. Finally we extend the theory to three dimensions and discuss the limit of a weak spanwise dependence for uniform wavetrains, showing that two functions are required to describe the evolution. These unknowns are a phase and a pressure function which satisfy a pair of linearly coupled partial differential equations. The results obtained from applying the same analysis to the fully three dimensional problem are included as an appendix.

Hewitt, R. E.↗
NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗
NASA NTRS1994

Effect of stratification and geometrical spreading on sonic boom rise time

The purpose of our investigation is to determine the effect of unsteadiness (not associated with turbulence) on rise time. The unsteadiness considered here is due to (1) geometrical spreading, (2) stratification, which includes variation in density, temperature, and relative humidity, and (3) N shaped waveform. A very general Burgers equation, which includes all these effects, is the propagation model for our study. The equation is solved by a new computational algorithm in which all the calculations are done in the time domain. The present paper is a progress report in which some of the factors contributing to unsteadiness are studied, namely geometrical spreading and variation in relative humidity. The work of Pierce and Kang, which motivated our study, is first reviewed. We proceed with a discussion of the Burgers equation model and the algorithm for solving the equation. Some comparison tests to establish the validity of the algorithm are presented. The algorithm is then used to determine the distance required for a steady-state shock, on encountering an abrupt change in relative humidity, to reach a new steady state based on the new humidity. It is found that the transition distance for plane shocks of amplitude 70 Pa is about 4 km when the change in relative humidity is 10 percent. Shocks of amplitude 140 Pa require less distance. The effect of spherical and cylindrical spreading is also considered. We demonstrate that a spreading shock wave never reaches steady state and that its rise time will be less than the equivalent steady state shock. Finally we show that an N wave has a slightly shorter rise time than a step shock of the same amplitude.

Cleveland, Robin O.↗