Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARABOLIC EQUATION”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Waveform synthesis of surface waves in a laterally heterogeneous earth by the Gaussian beam method

The present investigation is concerned with an application of the Gaussian beam method to surface waves in the laterally heterogeneous earth. The employed method has been developed for ray tracing and synthesizing seismograms of surface waves in cases involving the laterally heterogeneous earth. The procedure is based on formulations derived by Yomogida (1985). Vertical structure of the wave field is represented by the eigenfunctions of normal mode theory, while lateral variation is expressed by the parabolic equation as in two-dimensional acoustic waves or elastic body waves. It is demonstrated that a large-amplitude change can result from a slight perturbation in the phase velocity model.

Yomogida, K.↗
NASA NTRS1986

Multi-dimensional initial-boundary value problems with strong nonlinearities

Initial boundary value problems for the general scalar singularly perturbed parabolic equation in a cylindrical domain are studied. The one-dimensional reaction-diffusion equation is considered in the absence of convection. Then, the general one-dimensional and N-dimensional reaction-convection-diffusion equations are analyzed.

Howes, F. A.↗
NASA NTRS1986

Parameter identifiability under approximation

The problem of injectivity of the parameter-to-state map is discussed for Galerkin approximations of a linear parabolic equation. A necessary and sufficient condition is derived and illustrated by means of simple examples. Finally, output least squares identifiability of the Galerkin approximations is discussed.

Kunisch, K.↗
NASA NTRS1988

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗
NASA NTRS1988

An algorithm for unsteady flows with strong convection

An implicit ADI numerical method for the calculation of 2-D unsteady flows with strong convection effects is described. The method is based on the conventional Crank-Nicholson approach for parabolic equations but an upwind-downwind differencing is used for the first order spatial derivatives associated with convection. The differencing is carried out in the current and previous time plane in such a way that the algorithm is second order accurate in both space and time. The difference equations are factored into sequential operators, one in each independent spatial variable; the solution at each time step may then be computed as a sequence of tridiagonal matrix problems. The method may be used in a noniterative manner although iteration at each time step is recommended in situations where the effects of convection are strong.

Peridier, Vallorie↗
NASA NTRS1989

The solution of the Elrod algorithm for a dynamically loaded journal bearing using multigrid techniques

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed utilizing a multigrid iteration technique. The method is compared with a noniterative approach in terms of computational time and accuracy. The computational model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobsson-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via a smeared mass or striated flow extending to both surfaces in the film gap. The mixed nature of the equations (parabolic in the full film zone and hyperbolic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, C. M.↗
NASA NTRS1989

Acoustical effects of a large ridge on low-frequency sound propagation in stationary and moving atmospheres

The effects of a ridge on a low-frequency acoustic propagation in quiescent and windy atmospheres are investigated using a parabolic approximation. A logarithmic wind-speed profile, commonly employed to model atmospheric wind currents, is modified and used to model two-dimensional atmospheric flow over a triangularly-shaped hill. The parabolic equation is solved using an implicit finite-difference algorithm. Several examples are examined to determine the combined effects of source-ridge distance, ridge dimensions, wind-speed profile, and CW source frequency on the received acoustic field.

Robertson, J. S.↗
NASA NTRS1990

PTA en route noise measurements

A long-range advanced turboprop en route noise database was obtained with weather, tracking, and onboard measurements. In-flight noise directivity measurements were made. Data repeatability within a test day was excellent. Day-to-day variability existed and is not completely understood and therefore not predicted. Comparison of a two-dimensional ray-tracing propagation model with the ensemble average ground-measured data was good; however, as stated above, the day-to-day data variability was not completely predicted. Future research will include looking at alternative propagation models. Three-dimensional ray tracing, fast field program, and the parabolic equation are possibilities. The effect of turbulence needs to be accessed.

Willshire, William L., Jr.↗
NASA NTRS1990

Distorted-wave born approximation calculations for turbulence scattering in an upward-refracting atmosphere

Weiner and Keast observed that in an upward-refracting atmosphere, the relative sound pressure level versus range follows a characteristic 'step' function. The observed step function has recently been predicted qualitatively and quantitatively by including the effects of small-scale turbulence in a parabolic equation (PE) calculation. (Gilbert et al., J. Acoust. Soc. Am. 87, 2428-2437 (1990)). The PE results to single-scattering calculations based on the distorted-wave Born approximation (DWBA) are compared. The purpose is to obtain a better understanding of the physical mechanisms that produce the step-function. The PE calculations and DWBA calculations are compared to each other and to the data of Weiner and Keast for upwind propagation (strong upward refraction) and crosswind propagation (weak upward refraction) at frequencies of 424 Hz and 848 Hz. The DWBA calculations, which include only single scattering from turbulence, agree with the PE calculations and with the data in all cases except for upwind propagation at 848 Hz. Consequently, it appears that in all cases except one, the observed step function can be understood in terms of single scattering from an upward-refracted 'skywave' into the refractive shadow zone. For upwind propagation at 848 Hz, the DWBA calculation gives levels in the shadow zone that are much below both the PE and the data.

Gilbert, Kenneth E.↗
NASA NTRS1993

Comparisons of calculated and measured helicopter noise near instrument hill

The polar parabolic equation (POPE) method solves for the diffraction of sound by a curved surface including a realistic sound speed profile. POPE is outlined briefly to describe diffraction which propagates the field over a hill. Experimental data are compared with POPE predictions using the measured sound speed profile and ground impedance. Two trial cases are considered for the comparisons: the helicopter located at the base of the hill and far away from the base of the hill, respectively. The physical mechanisms for sound propagation over a hill are examined with and of POPE calculations and experimental data. The shedding of rays from the hillside gives an interference effect with a wave along the flat surface beyond the base of a hill.

Bass, Henry E.↗
NASA NTRS1993

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗
NASA NTRS1994

Microwave heating and joining of ceramic cylinders: A mathematical model

A thin cylindrical ceramic sample is placed in a single mode microwave applicator in such a way that the electric field strength is allowed to vary along its axis. The sample can either be a single rod or two rods butted together. We present a simple mathematical model which describes the microwave heating process. It is built on the assumption that the Biot number of the material is small, and that the electric field is known and uniform throughout the cylinder's cross-section. The model takes the form of a nonlinear parabolic equation of reaction-diffusion type, with a spatially varying reaction term that corresponds to the spatial variation of the electromagnetic field strength in the waveguide. The equation is analyzed and a solution is found which develops a hot spot near the center of the cylindrical sample and which then propagates outwards until it stabilizes. The propagation and stabilization phenomenon concentrates the microwave energy in a localized region about the center where elevated temperatures may be desirable.

Booty, Michael R.↗
NASA NTRS1994

Atmospheric propagation at larger lateral distances from the flight track

Sonic booms received on the ground tend to be restricted to a region of finite lateral extent below the flight track. This occurs because of refraction and because the effective speed of sound, even with winds taken into account, decreases with altitude in the lower atmosphere. Not all rays proceeding initially downwards from the flight track within an allowable range of initial directions will reach the ground. The restricted region which can be reached by rays impacting the ground is known as the primary carpet. However, weak rumbles are heard in the nominal shadow zone beyond the edge of this carpet. A full wave theory is necessary for explaining waveforms in that region, and the present paper gives a matched asymptotic expansion technique for a suitable approximate full wave theory that involves a relatively small number of parameters. The outer solution is derived from the structure of the system of rays that impact near the corridor edge; the inner solution involves a solution of the parabolic equation and results in the special functions encountered in the diffraction of sound over the tops of hills.

Pierce, Allan D.↗
NASA NTRS1995

Harmonic oscillators and resonance series generated by a periodic unstable classical orbit

The presence of an unstable periodic classical orbit allows one to introduce the decay time as a purely classical magnitude: inverse of the Lyapunov index which characterizes the orbit instability. The Uncertainty Relation gives the corresponding resonance width which is proportional to the Planck constant. The more elaborate analysis is based on the parabolic equation method where the problem is effectively reduced to the multidimensional harmonic oscillator with the time-dependent frequency. The resonances form series in the complex energy plane which is equidistant in the direction perpendicular to the real axis. The applications of the general approach to various problems in atomic physics are briefly exposed.

Kazansky, A. K.↗
NASA NTRS1995

qSmoothed particle interpolation

Smoothed particle hydrodynamics (SPH) discretization techniques are generalized to develop a method, smoothed particle interpolation (SPI), for solving initial value problems of systems of a nonhydrodynamical nature. Under this approach, SPH is viewed as strickly as interpolation scheme and, as such, suitable for solving general hyperbolic and parabolic equations. The SPI method is tested on (1) the wave equation with inhomogenous sound speed and (2) Burger's equation. The efficiency of SPI is studied by comparing SPI solutions to those obtained with standard finite-difference methods. It is shown that the power of SPI arises when the smoothing particles are free to move.

Laguna, Pablo↗
NASA NTRS1994

On the Dynamics of TVD Schemes

The dynamics of a class of TVD schemes for model hyperbolic and parabolic equations is studied numerically using a highly parallel supercomputer (CM-5). The objective is to utilize the highly parallel property of the CM-5 to reveal the reliable time step and entropy parameter ranges, and the degree of compressible flux limiters to avoid slow convergence and the production of nonphysical numerical solutions. We choose to study the nonlinear stability property of TVD schemes numerically since it is otherwise not amenable analytically.

Yee, H. C.↗
NASA NTRS2011

The Atmospheric Mutual Coherence Function From the First and Second Rytov Approximations and Its Comparison to That of Strong Fluctuation Theory

An expression for the mutual coherence function (MCF) of an electromagnetic beam wave propagating through atmospheric turbulence is derived within the confines of the Rytov approximation. It is shown that both the first and second Rytov approximations are required. The Rytov MCF is then compared to that which issues from the parabolic equation method of strong fluctuation theory. The agreement is found to be quite good in the weak fluctuation case. However, an instability is observed for the special case of beam wave intensities. The source of the instabilities is identified to be the characteristic way beam wave amplitudes are treated within the Rytov method.

Manning, Robert M.↗
NASA NTRSDate not supplied

Evaluation of Finite Impulse Response Filters for Turbulence Effects on Sonic Booms

Turbulence effects on sonic booms lead to random variability of sonic boom waveforms measured on the ground, complicating the prediction of such waveforms using sonic boom propagation codes that do not account for turbulence effects. The NASA PCBoom software is one such code shown to accurately predict sonic booms above the atmospheric boundary layer but not those on the ground. Analyses of measured sonic booms show that, on average, the turbulence effects result in sonic boom loudness reduction, on average, that increases with propagation distance and turbulence strength. To efficiently account for such effects, a signal processing-based approach has been developed at NASA using finite impulse response filters derived from predicted waveforms by solving a nonlinear parabolic equation. Analyses indicate that the mean loudness reduction obtained using the filters does not increase with propagation distance or with increasing turbulence strength, in disagreement with flight test data. An alternative approach to the filters is currently under development employing the multiple scattering theory (MST) of wave propagation. Available results from the physics-based MST approach are used to evaluate the predictive capability of the filters, suggesting that the filters can severely underestimate the prediction of the MST approach.

sonic booms↗