Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximations”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Jupiter's Great Red Spot and other vortices

A theoretical explanation of Jupiter's Great Red Spot (GRS) as the self-organization of vorticity in turbulence is presented. A number of properties of the GRS and other Jovian vortices that are unambiguous from the data are listed. The simplest possible model that explains these properties one at a time rather than in a difficult all-encompassing planetary global circulation model is presented. It is shown that Jovian vortices reflect the behavior of quasi-geostrophic (QG) vortices embedded in an east-west wind with bands of uniform potential vorticity. It is argued that most of the properties of the Jovian vortices can be easily explained and understood with QG theory. Many of the signatures of QG vortices are apparent on Voyager images. In numerical and laboratory experiments, QG vortices relax to approximately steady states like the Jovian vortices, rather than oscillating or rotating Kida ellipses.

Marcus, Philip S.↗
NASA NTRS1994

Mass and radius of cosmic balloons

Cosmic balloons are spherical domain walls with relativistic particles trapped inside. We derive the exact mass and radius relations for a static cosmic balloon using Gauss-Codazzi equations. The cosmic balloon mass as a function of its radius, M(R), is found to have a functional form similar to that of fermion soliton stars, with a fixed point at 2GM(R)/R approximately or equal to 0.486 which corresponds to the limit of infinite central density. We derive a simple analytical approximation for the mass density of a spherically symmetric relativistic gas star. When applied to the computation of the mass and radius of a cosmic balloon, the analytical approximation yields fairly good agreement with the exact numerical solutions.

Wang, Yun↗
NASA NTRS1994

X-ray clusters in a cold dark matter + lambda universe: A direct, large-scale, high-resolution, hydrodynamic simulation

A new, three-dimensional, shock-capturing, hydrodynamic code is utilized to determine the distribution of hot gas in a cold dark matter (CDM) + lambda model universe. Periodic boundary conditions are assumed: a box with size 85/h Mpc, having cell size 0.31/h Mpc, is followed in a simulation with 270(exp 3) = 10(exp 7.3) cells. We adopt omega = 0.45, lambda = 0.55, h identically equal to H/100 km/s/Mpc = 0.6, and then, from the cosmic background explorer (COBE) and light element nucleosynthesis, sigma(sub 8) = 0.77, omega(sub b) = 0.043. We identify the X-ray emitting clusters in the simulation box, compute the luminosity function at several wavelength bands, the temperature function and estimated sizes, as well as the evolution of these quantities with redshift. This open model succeeds in matching local observations of clusters in contrast to the standard omega = 1, CDM model, which fails. It predicts an order of magnitude decline in the number density of bright (h nu = 2-10 keV) clusters from z = 0 to z = 2 in contrast to a slight increase in the number density for standard omega = 1, CDM model. This COBE-normalized CDM + lambda model produces approximately the same number of X-ray clusters having L(sub x) greater than 10(exp 43) erg/s as observed. The background radiation field at 1 keV due to clusters is approximately the observed background which, after correction for numerical effects, again indicates that the model is consistent with observations.

Cen, Renyue↗
NASA NTRS1995

Ground-to-Flight Handling Qualities Comparisons for a High Performance Airplane

A flight test program was conducted in conjunction with a ground-based piloted simulation study to enable a comparison of handling qualities ratings for a variety of maneuvers between flight and simulation of a modern high performance airplane. Specific objectives included an evaluation of pilot-induced oscillation (PIO) tendencies and a determination of maneuver types which result in either good or poor ground-to-flight pilot handling qualities ratings. A General Dynamics F-16XL aircraft was used for the flight evaluations, and the NASA Langley Differential Maneuvering Simulator was employed for the ground based evaluations. Two NASA research pilots evaluated both the airplane and simulator characteristics using tasks developed in the simulator. Simulator and flight tests were all conducted within approximately a one month time frame. Maneuvers included numerous fine tracking evaluations at various angles of attack, load factors and speed ranges, gross acquisitions involving longitudinal and lateral maneuvering, roll angle captures, and an ILS task with a sidestep to landing. Overall results showed generally good correlation between ground and flight for PIO tendencies and general handling qualities comments. Differences in pilot technique used in simulator evaluations and effects of airplane accelerations and motions are illustrated.

Brandon, Jay M.↗
NASA NTRS1997

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗
NASA NTRS1996

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗
NASA NTRS2003

In General, the Total Voltammetric Current from a Mixture of Redox-Active Substances will Not be the Sum of the Currents that Each Substance would Produce Independently at the Same Concentration as in the Mixture

At the potential range where both decamethylferrocene (dMeFc) and ferrocene (Fc) are oxidized with rates controlled by linear diffusion, electrogenerated Fc(+) radicals diffusing outwards from the electrode react quantitatively (K23 C=5.8 x 10(exp 8) with dMeFc diffusing towards the electrode and produce Fc and dMeFc. That reaction replaces dMeFc with Fc, whose diffusion coefficient is higher than that of dMeFc(+), and the total mass-transfer limited current from the mixture is increased by approximately 10%. Analogous observations are made when mass-transfer is controlled by convective-diffusion as in RDE voltammetry. Similar results have been obtained with another, and for all practical purposes randomly selected pair of redox-active substances, [Co(bipy)3](2+) and N - methylphenothiazine (MePTZ); reaction of MePTZ(+) with [Co(bipy)3](2+) replaces the latter with MePTZ, which diffuses faster and the current increases by approximately 20%. The experimental voltammograms have been simulated numerically and the role of (a) the rate constant of the homogeneous reaction; (b) the relative concentrations; and, (c) the diffusion coefficients of all species involved have been studied in detail. Importantly, it was also identified that within any given redox system the dependence of the mass-transfer limited current on the bulk concentrations of the redox-active species is expected to be non-linear. These findings are discussed in terms of their electroanalytical implications.

Leventis, Nicholas↗
NASA NTRS2010

Multiple Exhaust Nozzle Effects on J-2X Gas Generator Outlet Impedance

The current test setup of the J-2X gas generator system uses a multiple nozzle configuration to exhaust hot gases to drive the propellant supply turbines. Combustion stability assessment of this gas generator design requires knowledge of the impedance effects the multiple nozzle configuration creates on the combustion chamber acoustic modes. Parallel work between NASA and Sierra Engineering is presented, showing two methods used to calculate the effective end impedance resulting from multiple nozzle configurations. The NASA method is a simple estimate of the effective impedance using the long wavelength approximation. Sierra Engineering has developed a more robust numerical integration method implemented in ROCCID to accommodate for multiple nozzles. Analysis using both methods are compared to J-2X gas generator test data collected over the past year.

Kenny, R. Jeremy↗
NASA NTRS1993

Linear and nonlinear PSE for compressible boundary layers

Compressible stability of growing boundary layers is studied by numerically solving the partial differential equations under a parabolizing approximation. The resulting parabolized stability equations (PSE) account for nonparallel as well as nonlinear effects. Evolution of disturbances in compressible flat-plate boundary layers are studied for freestream Mach numbers ranging from 0 to 4.5. Results indicate that the effect of boundary-layer growth is important for linear disturbances. Nonlinear calculations are performed for various Mach numbers. Two-dimensional nonlinear results using the PSE approach agree well with those from direct numerical simulations using the full Navier-Stokes equations while the required computational time is less by an order of magnitude. Spatial simulation using PSE were carried out for both the fundamental and subharmonic type breakdown for a Mach 1.6 boundary layer. The promising results obtained show that the PSE method is a powerful tool for studying boundary-layer instabilities and for predicting transition over a wide range of Mach numbers.

Chang, Chau-Lyan↗
DOE OSTI2024

Enhanced Estimation of Quantum Properties with Common Randomized Measurements

We present a technique for enhancing the estimation of quantum state properties by incorporating approximate prior knowledge about the quantum state. This consists in performing randomized measurements on a quantum processor and comparing the results with those obtained from a classical computer that stores an approximation of the quantum state. We provide unbiased estimators for expectation values of multicopy observables and present performance guarantees in terms of variance bounds that depend on the prior knowledge accuracy. We demonstrate the effectiveness of our approach through experimental and numerical examples detecting mixed-state entanglement, and estimating polynomial approximations of the von Neumann entropy and state fidelities. Published by the American Physical Society 2024

Physics↗
NASA NTRS1980

Numerical solution techniques for unsteady transonic aerodynamics problems

Basic concepts of finite difference solution techniques for unsteady transonic flows are presented. The hierarchy of mathematical forumulations that approximate the Navier-Stokes equations are reviewed. The basic concepts involved in constructing numerical algorthms to solve these formulations are given. Semi-implicit and implicit schemes are constructed and analyzed. The discussion focuses primarily on techniques for solving the low frequency transonic small disturbance equation. This is the simplest formulation that contains the essence of inviscid unsteady transonic flow physics. The low frequency formulation is emphasized here because codes based on this theory can be run in minutes of processor time on currently available computers. Furthermore, numerical techniques involved in solving this simple formulation also apply to the more complicated formulations. Extensions to these formulations are briefly described. An indication of the present capability for solving unsteady transonic flows is provided. Important areas of future research for the advancement of computational unsteady transonic aerodynamics are described.

Ballhaus, W. F.↗
NASA NTRS1978

Direct numerical solution of the transonic perturbation integral equation for lifting and nonlifting airfoils

The linear transonic perturbation integral equation previously derived for nonlifting airfoils is formulated for lifting cases. In order to treat shock wave motions, a strained coordinate system is used in which the shock location is invariant. The tangency boundary conditions are either formulated using the thin airfoil approximation or by using the analytic continuation concept. A direct numerical solution to this equation is derived in contrast to the iterative scheme initially used, and results of both lifting and nonlifting examples indicate that the method is satisfactory.

Nixon, D.↗
NASA NTRS1986

Mutual coupling effects in antenna arrays, volume 1

Mutual coupling between rectangular apertures in a finite antenna array, in an infinite ground plane, is analyzed using the vector potential approach. The method of moments is used to solve the equations that result from setting the tangential magnetic fields across each aperture equal. The approximation uses a set of vector potential model functions to solve for equivalent magnetic currents. A computer program was written to carry out this analysis and the resulting currents were used to determine the co- and cross-polarized far zone radiation patterns. Numerical results for various arrays using several modes in the approximation are presented. Results for one and two aperture arrays are compared against published data to check on the agreement of this model with previous work. Computer derived results are also compared against experimental results to test the accuracy of the model. These tests of the accuracy of the program showed that it yields valid data.

Collin, R. E.↗
DOE OSTI2022

Tri-Sectional Approximation of the Shortest Path to Long-Term Voltage Stability Boundary with Distributed Energy Resources

Ensuring long-term voltage stability is critical for reliable operations of power grids. High share of distributed energy resources (DERs) can create complicated system operation modes that may invalidate the traditional long-term voltage stability analysis based on typical operation modes. To address this challenge, this paper investigates how to compute the shortest path to the voltage stability boundary in the DER aggregated load space with large dispersion. Instead of working in the Euclidean space, we establish the analysis and computations on the algebraic power flow manifold to better capture the curvature change of the shortest path along the direction of losing stability. A modified optimal control framework is presented for obtaining the ground-truth of the smooth shortest path on the manifold. To efficiently and accurately solve for the shortest path, we further leverage the geometric features of the power flow manifold and propose a tri-sectional approximation model that is scalable for large-scale systems. Several numerical examples, up to the 1354-bus system, with different DER penetration levels and high dimensional renewable power injection variations are evaluated. The simulation results demonstrate that the tri-sectional approximation achieves high accuracy and efficiency to approximate the shortest path to the voltage stability boundary.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1986

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗
NASA NTRS1996

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗
NASA NTRS1987

The use of a numerical filter to correct airborne temperature measurements for the effects of sensor lag

A numerical filter for transforming measured temperature signals into a close approximation of the actual temperature signal is described. The filter is derived by minimizing the mean-square error of the system, and assuming a knowledge of the characteristics of the sensing element and its housing. The equation representing the frequency-response function of the numerical filter is given. Input and output spectra for a filter applied to a case with negligible noise and a noise level of 1.5 percent of the total power in the input spectrum are analyzed, and the numerical weights for these two cases are calculated. Phase angle and gain for the entire system are examined. It is noted that the filter can enhance spectral components as high as 8 Hz with little phase and gain degradation over the bandwidth.

Ritter, John A.↗
NASA NTRS1997

Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.

Barth, Timothy J.↗