Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PDE systems”

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 NTRS1994

Numerical Wave Propagation 211 Based on Wave Primitives

Compact higher order finite difference equations are applied to a sequence of problems in wave propagation and aeroacoustics. Systems of PDE's are reduced to a sequence of simple wave primitives using a local eigenvector decomposition. The wave primitives are first order PDE's in two independent variable and allow natural boundary conditions to be imposed for both single- and multidimensional problems. The method uses a "discrete dispersion relation" approach to obtain high order approximations to the wave primitives on a 3 spatial point / 2 time level computational molecule. The scheme is fourth order accurate for the class of system with constant coefficients, e.g., those that support exponential solutions. Weakly non-linear PDE's are solved in a similar manner using a variant of the "method of frozen coefficients." Experience with the new algorithm for linear, non-linear, and multi-dimensional test problems will be described.

Davis, Sanford S.↗
NASA NTRS1987

Parallels between control PDE's (Partial Differential Equations) and systems of ODE's (Ordinary Differential Equations)

System theorists understand that the same mathematical objects which determine controllability for nonlinear control systems of ordinary differential equations (ODEs) also determine hypoellipticity for linear partial differentail equations (PDEs). Moreover, almost any study of ODE systems begins with linear systems. It is remarkable that Hormander's paper on hypoellipticity of second order linear p.d.e.'s starts with equations due to Kolmogorov, which are shown to be analogous to the linear PDEs. Eigenvalue placement by state feedback for a controllable linear system can be paralleled for a Kolmogorov equation if an appropriate type of feedback is introduced. Results concerning transformations of nonlinear systems to linear systems are similar to results for transforming a linear PDE to a Kolmogorov equation.

Hunt, L. R.↗
DOE OSTI2026

Integral Kernel Methods for Nonlinear Parabolic-Elliptic Systems

Nonlinear parabolic-elliptic systems arise in many physical, biological, and chemical phenomena such as chemotaxis, ion transport, self-gravitating particles, and Brownian vortices. Existing methods struggle with the strong coupling and high nonlinearity and nonlocality of some of these systems, especially the ill-conditioned, convection-dominated problems. To overcome numerical difficulties, current approaches rely on initial guesses, preconditioning, or iterative techniques with no convergence guarantees. They might suffer from poor scalability, large memory usage, and difficulty to parallelize. Inspired by the connection of parabolic-elliptic systems to stochastic processes, we introduce a novel meshless, monolithic, and fully explicit method that naturally encapsulates the elliptic and parabolic operators into a single step which updates each node deterministically with global information. By being fully quadrature-based, it avoids solving systems of discretized equations and does not utilize initial guesses or preconditioning, while requiring little memory and being easy to parallelize. We first derive the method in an integral kernel formulation with quadratic complexity in the number of integration nodes and then leverage kernel-independent fast multipole methods (FMM) to present a scalable algorithm with linear complexity. We provide numerical examples for the Poisson-Nernst-Planck equations in one, two, and three dimensions, together with the derivation of the integral kernel for each case. Furthermore, the examples demonstrate the fast convergence and scalability of the FMM-accelerated algorithm, as well as its suitability for convection-dominated problems, making it competitive against traditional PDE solvers.

PDE systems↗
NASA NTRS1981

Application of adaptive grids to fluid-flow problems with asymptotic solutions

Coordinate system selection is an important consideration in the asymptotic numerical solution of any fluid-flow or heat transfer problem. This paper uses a new technique that provides a simple way of moving the mesh points in physical space in order to reduce the error in the computed asymptotic solution relative to that obtained using a fixed mesh. Applications to fluid-flow problems are presented, including boundary layer flow and inviscid supersonic flow over cylinders, and wedges with associated detached shocks. The treatment of curved boundaries, stationary and nonstationary boundaries, and systems of PDE's is discussed. Significant error reductions are demonstrated.

Rai, M. M.↗
NASA NTRS2003

Performance Enhancement of Unsteady Ejectors Investigated Using a Pulsejet Driver

Unsteady ejectors are currently under investigation for use in some pulse detonation engine (PDE) propulsion systems. This is due primarily to their potential high performance in comparison to steady ejectors of similar dimensions relative to the source or driver jet. Although some experimental work has been done in the past to study thrust augmentation with unsteady ejectors, there is no proven theory by which optimal design parameters can be selected and an effective ejector constructed for a given pulsed flow. Therefore, an experimental facility was developed at the NASA Glenn Research Center to study the correlation between ejector design and performance, and to get a better understanding of the flow phenomena that result in thrust augmentation. A commercially available pulsejet was used for the unsteady driving jet. This was paired with a basic, yet flexible, ejector design that allowed parametric evaluation of the effects that length, diameter, and inlet radius have on performance.

Paxson, Daniel E.↗
NASA NTRS2005

Studies of Operating Frequency Effects On Ejector-based Thrust Augmentation in a Pulse Detonation Engine

Studies were performed in order to characterize the thrust augmentation potential of an ejector in a Pulse Detonation Engine application. A 49-mm diameter tube of 0.914-m length was constructed with one open end and one closed end. Ethylene, oxygen, and nitrogen were introduced into the tube at the closed end through the implementation of a fast mixing injector. The tube was completely filled with a stoichiometric mixture containing a one to one molar ratio of nitrogen to oxygen. Ethylene was selected as the fuel due to its detonation sensitivity and the molar ratio of the oxidizer was chosen for heat transfer purposes. Detonations were initiated in the tube through the use of a spark ignition system. The PDE was operated in a multi-cycle mode at frequencies ranging from 20-Hz to 50-Hz. Baseline thrust measurements with no ejector present were performed while operating the engine at various frequencies and compared to theoretical estimates. The baseline values were observed to agree with the theoretical model at low operating frequencies and proved to be increasingly lower than the predicted values as the operating frequency was increased. The baseline thrust measurements were observed to agree within 15 percent of the model for all operating frequencies. A straight 152-mm diameter ejector was installed and thrust augmentation percentages were measured. The length of the ejector was varied while the overlap percentage (percent of the ejector length which overlapped the tube) was maintained at 25 percent for all tests. In addition, the effect of ejector inlet geometry was investigated by comparing results with a straight inlet to those of a 38-mm inlet diameter. The thrust augmentation of the straight inlet ejector proved to be independent of engine operating frequency, augmenting thrust by 40 percent for the 0.914-m length ejector. In contrast, the rounded lip ejector of the same length seemed to be highly dependent on the engine operating frequency. An optimum operating frequency observed with the rounded inlet occurred at an operating frequency of 30-Hz, resulting in thrust augmentation percentages greater than 100 percent. The effect that the engine operating frequency had on thrust augmentation levels attained with an ejector was characterized and optimum performance parameters were established. Insight into the frequency dependent nature of the ejector performance was pursued. Suggestions for future experiments which are needed to fully understand the means in which thrust augmentation is achieved in a PDE-ejector configuration were noted.

Landry, K.↗
NASA NTRS1996

Pulse Detonation Engine Air Induction System Analysis

A preliminary mixed-compression inlet design concept for potential pulse-detonation engine (PDE) powered supersonic aircraft was defined and analyzed. The objectives of this research were to conceptually design and integrate an inlet/PDE propulsion system into a supersonic aircraft, perform time-dependent CFD analysis of the inlet flowfield, and to estimate the installed PDE cycle performance. The study was baselined to a NASA Mach 5 Waverider study vehicle in which the baseline over/under turboramjet engines were replaced with a single flowpath PDE propulsion system. As much commonality as possible was maintained with the baseline configuration, including the engine location and forebody lines. Modifications were made to the inlet system's external ramp angles and a rotating cowl lip was incorporated to improve off-design inlet operability and performance. Engines were sized to match the baseline vehicle study's ascent trajectory thrust requirement at Mach 1.2. The majority of this study was focused on a flight Mach number of 3.0. The time-dependent Navier Stokes CFD analyses of a two-dimensional approximation of the inlet was conducted for the Mach 3.0 condition. The Lockheed Martin Tactical Aircraft Systems-developed FALCON CFD code with a two equation 'k-1' turbulence model was used. The downstream PDE was simulated by an array of four sonic nozzles in which the flow areas were rapidly varied in various opening/closing combinations. Results of the CFD study indicated that the inlet design concept operated successfully at the Mach 3.0 condition, satisfying mass capture, total pressure recovery, and operability requirements. Time-dependent analysis indicated that pressure and expansion waves from the simulated valve perturbations did not effect the inlet's operability or performance.

Pegg, R. J.↗
NASA NTRS1989

An investigation of a mathematical model of an optically pumped Ti(3+):Al2O3 laser system

During the last several years, solid state lasers were developed that have the potential for meeting rigorous performance requirements for space-based remote sensing of the atmosphere. In order to design a stable and efficient laser and to understand the effect on laser output of changes in the physical and design parameters, an understanding of the development of the dynamical processes of the laser is necessary. Typically, the dynamical processes in a laser system are investigated via rate equations describing the evolution of the occupancy in the electronic levels and of the photon density in the laser cavity. There are two approaches to this type of study. Most often, for the sake of simplicity, the spatial variations of the dynamic variables in the laser system are disregarded and the mathematical model consists of a system of first order nonlinear ordinary differential equations (ODE). The second approach is to take into account both spatial and temporal variations in the dynamic variables in the laser cavity. The resulting model consists of a first order semilinear system of partial differential equations (PDE). The model which was studied was studied was generic in the sense that it was a four-level laser system, but the parameters used in the numerical study were specific to Titanium-doped sapphire. For simplicity, a constant, spatially uniform pumping scheme was considered. In addition, a simplification of the model was made so that it treats a single lasing wavelength with a narrow bandwidth. The purpose was to investigate both versions of the mathematical model and to determine whether the numerical solutions are similar both qualitatively and quantitatively. The systems of ordinary differential equations were solved numerically using a Runge-Kutta-Fehlberg algorithm which was very efficient for typical values of the physical parameters. A numerical scheme, based on the Modified Euler method, for computing solutions to the system of partial differential equations was developed and implemented. The PDE model was solved numerically at the expense of greatly increased computer time.

Roberts, Lila F.↗
NASA NTRS2001

Discontinuous Galerkin Methods for NonLinear Differential Systems

This talk considers simplified finite element discretization techniques for first-order systems of conservation laws equipped with a convex (entropy) extension. Using newly developed techniques in entropy symmetrization theory, simplified forms of the discontinuous Galerkin (DG) finite element method have been developed and analyzed. The use of symmetrization variables yields numerical schemes which inherit global entropy stability properties of the PDE (partial differential equation) system. Central to the development of the simplified DG methods is the Eigenvalue Scaling Theorem which characterizes right symmetrizers of an arbitrary first-order hyperbolic system in terms of scaled eigenvectors of the corresponding flux Jacobian matrices. A constructive proof is provided for the Eigenvalue Scaling Theorem with detailed consideration given to the Euler equations of gas dynamics and extended conservation law systems derivable as moments of the Boltzmann equation. Using results from kinetic Boltzmann moment closure theory, we then derive and prove energy stability for several approximate DG fluxes which have practical and theoretical merit.

Barth, Timothy↗
NASA NTRS2003

Pulsed Ejector Wave Propogation Test Program

The development of, and initial test data from, a nondetonating Pulse Detonation Engine (PDE) simulator tested in the NASA Glenn 1 x 1 foot Supersonic Wind Tunnel (SWT) is presented in this paper. The concept is a pulsed ejector driven by the simulated exhaust of a PDE. This pro- gram is applicable to a PDE entombed in a ramjet flowpath, i.e., a PDE combined-cycle propulsion system. The ejector primary flow is a pulsed, uiiderexpanded, supersonic nozzle simulating the supersonic waves ema- nating from a PDE, while the ejector secondary flow is the 1 x 1 foot SWT test section operated at subsonic Mach numbers. The objective is not to study the detonation details, but the wave physics including t,he start- ing vortices, the extent of propagation of the wave front, the reflection of the wave from the secondary flowpath walls, and the timing of these events of a pulsed ejector, and correlate these with Computational Fluid Dynamics (CFD) code predictions. Pulsed ejectors have been shown to result in a 3 to 1 improvement in LID (length-to-diameter) and a near 2 to 1 improvement in thrust augmentation over a steady ejector. This program will also explore the extent of upstream interactions between an inlet and large, periodically applied, backpressures to the inlet as would be present due to combustion tube detonations in a PDE. These interactions could result in inlet unstart or buzz for a supersonic mixed compression inlet. The design of the present experiment entailed the use of an 2-t diagram characteristics code to study the nozzle filling and purging timescales as well as a series of CFD analyses conducted using the WIND code. The WIND code is a general purpose CFD code for solution of the Reynolds averaged Navier-Stokes equations and can be applied to both steady state and time-accurate calculations. The first, proof-of-concept, test entry (spring 2001) pressure distributions shown here indicate the simulation concept was successful and therefore the experimental approach is sound.

Fernandez, Rene↗
NASA NTRS1993

A semi-Lagrangian approach to the shallow water equation

We present a formulation of the shallow water equations that emphasizes the conservation of potential vorticity. A locally conservative semi-Lagrangian time-stepping scheme is developed, which leads to a system of three coupled PDE's to be solved at each time level. We describe a smoothing analysis of these equations, on which an effective multigrid solver is constructed. Some results from applying this solver to the static version of these equations are presented.

Bates, J. R.↗
NASA NTRS2002

Innovative Airbreathing Propulsion Concepts for High-speed Applications

The current cost to launch payloads to low earth orbit (LEO) is approximately loo00 U.S. dollars ($) per pound ($22000 per kilogram). This high cost limits our ability to pursue space science and hinders the development of new markets and a productive space enterprise. This enterprise includes NASA's space launch needs and those of industry, universities, the military, and other U.S. government agencies. NASA's Advanced Space Transportation Program (ASTP) proposes a vision of the future where space travel is as routine as in today's commercial air transportation systems. Dramatically lower launch costs will be required to make this vision a reality. In order to provide more affordable access to space, NASA has established new goals in its Aeronautics and Space Transportation plan. These goals target a reduction in the cost of launching payloads to LEO to $lo00 per pound ($2200 per kilogram) by 2007 and to $100'~ per pound by 2025 while increasing safety by orders of magnitude. Several programs within NASA are addressing innovative propulsion systems that offer potential for reducing launch costs. Various air-breathing propulsion systems currently are being investigated under these programs. The NASA Aerospace Propulsion and Power Base Research and Technology Program supports long-term fundamental research and is managed at GLenn Research Center. Currently funded areas relevant to space transportation include hybrid hyperspeed propulsion (HHP) and pulse detonation engine (PDE) research. The HHP Program currently is addressing rocket-based combined cycle and turbine-based combined cycle systems. The PDE research program has the goal of demonstrating the feasibility of PDE-based hybrid-cycle and combined cycle propulsion systems that meet NASA's aviation and access-to-space goals. The ASTP also is part of the Base Research and Technology Program and is managed at the Marshall Space Flight Center. As technologies developed under the Aerospace Propulsion and Power Base Research and Technology Program mature, they are incorporated into ASTP. One example of this is rocket-based combined cycle systems that are being considered as part of ASTP. The NASA Ultra Efficient Engine Technology (UEET) Program has the goal of developing propulsion system component technology that is relevant to a wide range of vehicle missions. In addition to subsonic and supersonic speed regimes, it includes the hypersonic speed regime. More specifically, component technologies for turbine-based combined cycle engines are being developed as part of UEET.

Whitlow, Woodrow, Jr.↗
NASA NTRS1999

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗
NASA NTRS1993

An automatic multigrid method for the solution of sparse linear systems

An automatic version of the multigrid method for the solution of linear systems arising from the discretization of elliptic PDE's is presented. This version is based on the structure of the algebraic system solely, and does not use the original partial differential operator. Numerical experiments show that for the Poisson equation the rate of convergence of our method is equal to that of classical multigrid methods. Moreover, the method is robust in the sense that its high rate of convergence is conserved for other classes of problems: non-symmetric, hyperbolic (even with closed characteristics) and problems on non-uniform grids. No double discretization or special treatment of sub-domains (e.g. boundaries) is needed. When supplemented with a vector extrapolation method, high rates of convergence are achieved also for anisotropic and discontinuous problems and also for indefinite Helmholtz equations. A new double discretization strategy is proposed for finite and spectral element schemes and is found better than known strategies.

Shapira, Yair↗
DOE OSTI2026

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Optical properties of materials at low temperature and their application to optical detection

A lumped model to represent the photodielectric effect is developed. An analog simulation for a sample in a microwave cavity with a static magnetic field is developed. A system to measure continuously the PDE is analyzed. A performance factor to compared PD detectors versus ac photoconductors is computed. The operating conditions are defined for the appropriate noise conditions. The detectivity of the detector is found to be limited by the semiconductor sample noise.

Hartwig, W. H.↗
NASA NTRS1995

The estimation of material and patch parameters in a PDE-based circular plate model

The estimation of material and patch parameters for a system involving a circular plate, to which piezoceramic patches are bonded, is considered. A partial differential equation (PDE) model for the thin circular plate is used with the passive and active contributions form the patches included in the internal and external bending moments. This model contains piecewise constant parameters describing the density, flexural rigidity, Poisson ratio, and Kelvin-Voigt damping for the system as well as patch constants and a coefficient for viscous air damping. Examples demonstrating the estimation of these parameters with experimental acceleration data and a variety of inputs to the experimental plate are presented. By using a physically-derived PDE model to describe the system, parameter sets consistent across experiments are obtained, even when phenomena such as damping due to electric circuits affect the system dynamics.

Banks, H. T.↗