Search NASA⌕ Search

SEARCH · Search NASA

Results for “dispersive 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗
NASA NTRS2001

Computational Aeroacoustics by the Space-time CE/SE Method

In recent years, a new numerical methodology for conservation laws-the Space-Time Conservation Element and Solution Element Method (CE/SE), was developed by Dr. Chang of NASA Glenn Research Center and collaborators. In nature, the new method may be categorized as a finite volume method, where the conservation element (CE) is equivalent to a finite control volume (or cell) and the solution element (SE) can be understood as the cell interface. However, due to its rigorous treatment of the fluxes and geometry, it is different from the existing schemes. The CE/SE scheme features: (1) space and time treated on the same footing, the integral equations of conservation laws are solve( for with second order accuracy, (2) high resolution, low dispersion and low dissipation, (3) novel, truly multi-dimensional, simple but effective non-reflecting boundary condition, (4) effortless implementation of computation, no numerical fix or parameter choice is needed, an( (5) robust enough to cover a wide spectrum of compressible flow: from weak linear acoustic waves to strong, discontinuous waves (shocks) appropriate for linear and nonlinear aeroacoustics. Currently, the CE/SE scheme has been developed to such a stage that a 3-13 unstructured CE/SE Navier-Stokes solver is already available. However, in the present paper, as a general introduction to the CE/SE method, only the 2-D unstructured Euler CE/SE solver is chosen as a prototype and is sketched in Section 2. Then applications of the CE/SE scheme to linear, nonlinear aeroacoustics and airframe noise are depicted in Sections 3, 4, and 5 respectively to demonstrate its robustness and capability.

Loh, Ching Y.↗
NASA NTRS2013

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel fl ow problem.

Morris, C. I.↗
NASA NTRS2013

Studies of Inviscid Flux Schemes for Acoustics and Turbulence Problems

The last two decades have witnessed tremendous growth in computational power, the development of computational fluid dynamics (CFD) codes which scale well over thousands of processors, and the refinement of unstructured grid-generation tools which facilitate rapid surface and volume gridding of complex geometries. Thus, engineering calculations of 10(exp 7) - 10(exp 8) finite-volume cells have become routine for some types of problems. Although the Reynolds Averaged Navier Stokes (RANS) approach to modeling turbulence is still in extensive and wide use, increasingly large-eddy simulation (LES) and hybrid RANS-LES approaches are being applied to resolve the largest scales of turbulence in many engineering problems. However, it has also become evident that LES places different requirements on the numerical approaches for both the spatial and temporal discretization of the Navier Stokes equations than does RANS. In particular, LES requires high time accuracy and minimal intrinsic numerical dispersion and dissipation over a wide spectral range. In this paper, the performance of both central-difference and upwind-biased spatial discretizations is examined for a one-dimensional acoustic standing wave problem, the Taylor-Green vortex problem, and the turbulent channel ow problem.

Morris, Christopher I.↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
NASA NTRSDate not supplied

Mars Sample Return Earth Entry System Helicopter Drop Test Reconstruction

The Mars Sample Return campaign conducted four successful helicopter drop tests of the Earth Entry System at the Utah Test and Training Range on October, 2023. The tests acquired data for the 52.5 deg sphere cone geometry that can be used to develop models for the terminal descent aerodynamics and flight dynamics used to model ground impact conditions. The acquired test data included both sensor data outputs from an on-board inertial measurement unit, global positioning system, and video cameras as well as atmospheric measurements from weather balloons. This data was processed using a Kalman filter/smoother to reconstruct the capsule’s trajectory and aerodynamics. An equation-error method was used to reconcile the aerodynamics by solving for a set of dispersions in the aerodynamic database that form a best-fit to the reconstructed aerodynamics. The reconstructed trajectories were compared with flight simulations generated using Program to Optimize Simulated Trajectories II. The results indicate that the simulation underpredicted the vehicle oscillation amplitudes. The reconstruction exhibits a persistent low-amplitude oscillation that does not damp out as predicted.

Chris D Karlgaard↗
NASA NTRSDate not supplied

Mars Sample Return Earth Entry System Helicopter Drop Test Reconstruction

The Mars Sample Return campaign conducted four successful helicopter drop tests of the Earth Entry System at the Utah Test and Training Range on October, 2023. The tests acquired data for the 52.5 deg sphere cone geometry that can be used to develop models for the terminal descent aerodynamics and flight dynamics used to model ground impact conditions. The acquired test data included both sensor data outputs from an on-board inertial measurement unit, global positioning system, and video cameras as well as atmospheric measurements from weather balloons. This data was processed using a Kalman filter/smoother to reconstruct the capsule’s trajectory and aerodynamics. An equation-error method was used to reconcile the aerodynamics by solving for a set of dispersions in the aerodynamic database that form a best-fit to the reconstructed aerodynamics. The reconstructed trajectories were compared with flight simulations generated using Program to Optimize Simulated Trajectories II. The results indicate that the simulation underpredicted the vehicle oscillation amplitudes. The reconstruction exhibits a persistent low-amplitude oscillation that does not damp out as predicted.

Chris D Karlgaard↗
NASA NTRS1991

Computational dispersion properties of horizontal staggered grids for atmospheric and ocean models

The computational dispersion properties of horizontally and time-horizontally staggered grids utilizing corresponding centered-difference techniques for approximation of the adjustment, or gravity wave equations, are examined in terms of their group velocity characteristics. Results are acquired for oceanic and atmospheric models, the former being characterized by a much smaller Rossby radius of deformation. For all grids considered additional filtering is required to control and even eliminate waves with poor computational dispersion characteristics. Computational dispersion properties along with other computational characteristics and requirements give some guidance for an optimal selection of an appropriate grid for an ocean or atmospheric model.

Fox-Rabinovitz, Michael S.↗
NASA NTRS1976

Numerical problems in the advection of pollutants

Prediction of pollutant transport involves solution of the advection equation numerically. Truncation errors arising from integrating the advection equation by using finite difference schemes are discussed. Numerical errors manifest themselves in artificial dissipation and dispersion of pollutants. Criteria for measuring the dissipation and dispersion are presented, and minimization of those errors is discussed. An alternative procedure by integration in the wavenumber space is examined. Examples illustrating quantitative discrepancy in actual forecasting are also shown.

Liu, C. Y.↗
NASA NTRS2000

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Mohseni, Kamran↗
NASA NTRS1993

A numerical treatment of geodynamic viscous flow problems involving the advection of material interfaces

Effective numerical treatment of multicomponent viscous flow problems involving the advection of sharp interfaces between materials of differing physical properties requires correction techniques to prevent spurious diffusion and dispersion. We develop a particular algorithm, based on modern shock-capture techniques, employing a two-step nonlinear method. The first step involves the global application of a high-order upwind scheme to a hyperbolic advection equation used to model the distribution of distinct material components in a flow field. The second step is corrective and involves the application of a global filter designed to remove dispersion errors that result from the advection of discontinuities (e.g., material interfaces) by high-order, minimally dissipative schemes. The filter introduces no additional diffusion error. Nonuniform viscosity across a material interface is allowed for by the implementation of a compositionally weighted-inverse interface viscosity scheme. The combined method approaches the optimal accuracy of modern shock-capture techniques with a minimal increase in computational time and memory. A key advantage of this method is its simplicity to incorporate into preexisting codes be they finite difference, element, or volume of two or three dimensions.

Lenardic, A.↗
NASA NTRS1974

Scattering and absorption from poly-dispersed aerosols

Results from numerical computations of the diffuse sky radiances and degree of polarization are given. These results were obtained by solving the equation of radiative transfer in the formulation of Eschelbach. In these computations, multiple scattering and absorption by aerosol particles were considered. In addition, results from experimental measurements of the complex refractive index of aerosol samples at different relative humidities are presented. From the radiation flux divergences which were computed based on the determined properties of atmospheric aerosol particles, atmospheric heating rates were derived which were found to be comparable to the heating rates by water vapor.

Bullrich, K.↗
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 NTRS1994

A double-moment multiple-phase four-class bulk ice scheme, Part 1: Description

A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Ferrier, Brad Shoenberg↗
NASA NTRS1998

Effect of Helical Slow-Wave Circuit Variations on TWT Cold-Test Characteristics

Recent advances in the state of the art of computer modeling offer the possibility for the first time to evaluate the effect that slow-wave structure parameter variations, such'as manufacturing tolerances, have on the cold-test characteristics of helical traveling-wave tubes (TWT's). This will enable manufacturers to determine the cost effectiveness of controlling the dimensions of the component parts of the TWT, which is almost impossible to do experimentally without building a large number of tubes and controlling several parameters simultaneously. The computer code MAxwell's equations by the Finite Integration Algorithm (MAFIA) is used in this analysis to determine the effect on dispersion and on-axis interaction impedance of several helical slow-wave circuit parameter variations, including thickness and relative dielectric constant of the support rods, tape width, and height of the metallized films deposited on the dielectric rods. Previous computer analyzes required so many approximations that accurate determinations of the effect of many relevant dimensions on tube performance were practically impossible.

Kory, Carol L.↗
NASA NTRS2011

Progress Report on Optimizing X-ray Optical Prescriptions for Wide-Field Applications

We report on the present status of our continuing efforts to develop a method for optimizing wide-field nested x-ray telescope mirror prescriptions. Utilizing extensive Monte-Carlo ray trace simulations, we find an analytic form for the root-mean-square dispersion of rays from a Wolter I optic on the surface of a flat focal plane detector as a function of detector tilt away from the nominal focal plane and detector displacement along the optical axis. The configuration minimizing the ray dispersion from a nested array of Wolter I telescopes is found by solving a linear system of equations for tilt and individual mirror pair displacement. Finally we outline our initial efforts at expanding this method to include higher order polynomial terms in the mirror prescriptions.

Elsner, R. F.↗
NASA NTRS1997

Computations of Complex Three-Dimensional Turbulent Free Jets

Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Wilson, Robert V.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗