SEARCH · Engineering Papers

Results for “time discretization”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Numerical Solution of Incompressible Navier-Stokes Equations Using a Fractional-Step Approach

A fractional step method for the solution of steady and unsteady incompressible Navier-Stokes equations is outlined. The method is based on a finite volume formulation and uses the pressure in the cell center and the mass fluxes across the faces of each cell as dependent variables. Implicit treatment of convective and viscous terms in the momentum equations enables the numerical stability restrictions to be relaxed. The linearization error in the implicit solution of momentum equations is reduced by using three subiterations in order to achieve second order temporal accuracy for time-accurate calculations. In spatial discretizations of the momentum equations, a high-order (3rd and 5th) flux-difference splitting for the convective terms and a second-order central difference for the viscous terms are used. The resulting algebraic equations are solved with a line-relaxation scheme which allows the use of large time step. A four color ZEBRA scheme is employed after the line-relaxation procedure in the solution of the Poisson equation for pressure. This procedure is applied to a Couette flow problem using a distorted computational grid to show that the method minimizes grid effects. Additional benchmark cases include the unsteady laminar flow over a circular cylinder for Reynolds Numbers of 200, and a 3-D, steady, turbulent wingtip vortex wake propagation study. The solution algorithm does a very good job in resolving the vortex core when 5th-order upwind differencing and a modified production term in the Baldwin-Barth one-equation turbulence model are used with adequate grid resolution.

Kiris, Cetin↗

NASA NTRSDate not supplied

The Essence of Reactivity

Reactive programming, functional reactive programming, event-based programming, stream programming, and temporal logic all share an underlying commonality: values can vary over time. These languages differ in multiple ways, including the nature of time itself (e.g., continuous or discrete, dense or sparse, implicit or explicit), on how much of the past and future can be referenced, on the kinds of values that can be represented, as well as the mechanisms used to evaluate expressions or formulas. This paper presents a series of abstractions that capture the essence of different forms of time variance. By separating the aspects that differentiate each family of formalisms, we can better express the commonalities and differences between them. We demonstrate our work with a prototype in Haskell that allows us to write programs in terms of a generic interface that can be later instantiated to different abstractions depending on the desired target.

reactive programming↗

NASA NTRS1972

Numerical methods for solving the Boltzmann equation (a review)

The methods are reviewed which are utilized in principal attempts to obtain the numerical solution or modeling of the Boltzmann equation over a broad range of Knudsen numbers. The primary methods considered are the Monte Carlo and the discrete velocities methods. The conculsions drawn from the analysis include the following: (1) The Monte Carlo methods are not well suited in the area of small Knudsen numbers. (2) Among the Monte Carlo methods, the Bird method appears to be the most attractive, since it is more directly related to the Boltzmann equation. (3) The deterministic methods, which include the discrete ordinate technique, offer great possibilities but require exceedingly large computer times. (4) The use of approximating equations in combination with the discrete velocities method will possibly improve computation time and reduce the required memory volume.

Limar, Y. F.↗

NASA NTRS1991

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗

NASA NTRS1992

Finite element Euler calculations of unsteady transonic cascade flows

A Galerkin finite element procedure incorporating an explicit Runge-Kutta time-stepping scheme has been developed in this work to solve unsteady transonic flow in cascades. The computational domain is discretized by a globally unstructured but locally structured blade-fitted deformable mesh. The Galerkin approximation is applied to the unsteady Euler equations based on a mixed Eulerian-Lagrangian description. The semi-discretized equations are integrated forward in time using a multistage Runge-Kutta scheme. An artificial dissipation operator of the type proposed by Jameson is adapted in the current scheme to capture shocks and suppress nonphysical oscillations. Phase-shifted boundary conditions are used to reduce the computational domain to a single reference passage. Results for both steady and unsteady transonic flows through cascades are presented and compared to existing finite volume solutions.

Hsiao, Chingteng↗

NASA NTRS1988

Singular perturbations and time scales in the design of digital flight control systems

The results are presented of application of the methodology of Singular Perturbations and Time Scales (SPATS) to the control of digital flight systems. A block diagonalization method is described to decouple a full order, two time (slow and fast) scale, discrete control system into reduced order slow and fast subsystems. Basic properties and numerical aspects of the method are discussed. A composite, closed-loop, suboptimal control system is constructed as the sum of the slow and fast optimal feedback controls. The application of this technique to an aircraft model shows close agreement between the exact solutions and the decoupled (or composite) solutions. The main advantage of the method is the considerable reduction in the overall computational requirements for the evaluation of optimal guidance and control laws. The significance of the results is that it can be used for real time, onboard simulation. A brief survey is also presented of digital flight systems.

Naidu, Desineni S.↗

NASA NTRS1981

Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

Jameson, A.↗

NASA NTRS1987

Anomalous transport in discrete arcs and simulation of double layers in a model auroral circuit

The evolution and long-time stability of a double layer in a discrete auroral arc requires that the parallel current in the arc, which may be considered uniform at the source, be diverted within the arc to charge the flanks of the U-shaped double-layer potential structure. A simple model is presented in which this current re-distribution is effected by anomalous transport based on electrostatic lower hybrid waves driven by the flank structure itself. This process provides the limiting constraint on the double-layer potential. The flank charging may be represented as that of a nonlinear transmission. A simplified model circuit, in which the transmission line is represented by a nonlinear impedance in parallel with a variable resistor, is incorporated in a 1-d simulation model to give the current density at the DL boundaries. Results are presented for the scaling of the DL potential as a function of the width of the arc and the saturation efficiency of the lower hybrid instability mechanism.

Smith, Robert A.↗

NASA NTRS1990

Three-dimensional flux-split Euler schemes involving unstructured dynamic meshes

Improved algorithms for the solution of the three-dimensional time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multi-stage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the ONERA M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗

NASA NTRS1991

Accuracy of an unstructured-grid upwind-Euler algorithm for the ONERA M6 wing

Improved algorithms for the solution of the three-dimensional, time-dependent Euler equations are presented for aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured-grid flow solvers. The spatial discretization involves a flux-split approach that is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves either an explicit time-integration scheme using a multistage Runge-Kutta procedure or an implicit time-integration scheme using a Gauss-Seidel relaxation procedure, which is computationally efficient for either steady or unsteady flow problems. With the implicit Gauss-Seidel procedure, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady flow results are presented for both the NACA 0012 airfoil and the Office National d'Etudes et de Recherches Aerospatiales M6 wing to demonstrate applications of the new Euler solvers. The paper presents a description of the Euler solvers along with results and comparisons that assess the capability.

Batina, John T.↗

NASA NTRS2017

Flight Deck Surface Trajectory-Based Operations

Surface Trajectory-Based Operations (STBO) is a future concept for surface operations where time requirements are incorporated into taxi operations to support surface planning and coordination. Pilot-in-the-loop flight deck simulations have been conducted to study flight deck displays algorithms to aid pilots in complying with the time requirements of time-based taxi operations (i.e., at discrete locations in 3 12 D operations or at all points along the route in 4DT operations). The results of these studies (conformance, time-of-arrival error, eye-tracking data, and safety ratings) are presented. Flight deck simulation work done in collaboration with DLR is described. Flight deck research issues in future auto-taxi operations are also introduced.

flight deck↗

NASA NTRS1991

Survival of discrete structures in the solar wind

Survival time scales are estimated for structures ejected from the sun into the solar wind, which are subject to decay processes during their transit from the sun to the earth. With the assumption that the structures are magnetic, MHD estimates are made to determine the time scales for reconnection into the ambient magnetic field. It is found that for structures smaller than a few tens of Mm at the sun, the probability of surviving the transit to the earth is low. The transition in size between structures that survive the sun-earth transit and those that do not (20-30 Mm at the sun) corresponds to scales of (4-6) x 10 to the 11th cm at the earth orbit (assuming radial expansion).

Mullan, D. J.↗

NASA NTRS2017

Global Observations of Magnetospheric High-m Poloidal Waves During the 22 June 2015 Magnetic Storm

We report global observations of high-m poloidal waves during the recovery phase of the 22 June 2015 magnetic storm from a constellation of widely spaced satellites of five missions including Magnetospheric Multiscale (MMS), Van Allen Probes, Time History of Events and Macroscale Interactions during Substorm (THEMIS), Cluster, and Geostationary Operational Environmental Satellites (GOES). The combined observations demonstrate the global spatial extent of storm time poloidal waves. MMS observations confirm high azimuthal wave numbers (m approximately 100). Mode identification indicates the waves are associated with the second harmonic of field line resonances. The wave frequencies exhibit a decreasing trend as L increases, distinguishing them from the single-frequency global poloidal modes normally observed during quiet times. Detailed examination of the instantaneous frequency reveals discrete spatial structures with step-like frequency changes along L. Each discrete L shell has a steady wave frequency and spans about 1 RE, suggesting that there exist a discrete number of drift-bounce resonance regions across L shells during storm times.

Le, G.↗

NASA NTRS1974

Lateral conduction effects on heat-transfer data obtained with the phase-change paint technique

A computerized tool, CAPE, (Conduction Analysis Program using Eigenvalues) has been developed to account for lateral heat conduction in wind tunnel models in the data reduction of the phase-change paint technique. The tool also accounts for the effects of finite thickness (thin wings) and surface curvature. A special reduction procedure using just one time of melt is also possible on leading edges. A novel iterative numerical scheme was used, with discretized spatial coordinates but analytic integration in time, to solve the inverse conduction problem involved in the data reduction. A yes-no chart is provided which tells the test engineer when various corrections are large enough so that CAPE should be used. The accuracy of the phase-change paint technique in the presence of finite thickness and lateral conduction is also investigated.

Maise, G.↗

NASA NTRS1971

A method for determining the response of space shuttle to atmospheric turbulence. Volume 1: Space shuttle turbulence response

A computer program has been developed and demonstrated that can analyze the response of space shuttle to atmospheric turbulence. The method developed accounts for propellant slosh, gimballed engine, and stability augmentation system (SAS) coupling with the elastic vehicle. Statistical outputs are generated that relate vehicle loads and accelerations to level of random turbulence. For discrete turbulence descriptions, response time-histories are computed. Preliminary turbulence response analyses of space shuttle were conducted for one ascent and one booster flyback subsonic flight condition. The ascent case considered both symmetric and antisymmetric boundary conditions, while only the symmetric analysis was conducted for booster flyback. The three conditions were analyzed with SAS active and inactive, but apparent instabilities with SAS active (arising from improper gains) invalidated the results. Load and acceleration responses in all cases, however, were well within vehicle design limits.

Source record↗

NASA NTRS1974

Computer programs for pressurization (RAMP) and pressurized expulsion from a cryogenic liquid propellant tank

An analysis to predict the pressurant gas requirements for the discharge of cryogenic liquid propellants from storage tanks is presented, along with an algorithm and two computer programs. One program deals with the pressurization (ramp) phase of bringing the propellant tank up to its operating pressure. The method of analysis involves a numerical solution of the temperature and velocity functions for the tank ullage at a discrete set of points in time and space. The input requirements of the program are the initial ullage conditions, the initial temperature and pressure of the pressurant gas, and the time for the expulsion or the ramp. Computations are performed which determine the heat transfer between the ullage gas and the tank wall. Heat transfer to the liquid interface and to the hardware components may be included in the analysis. The program output includes predictions of mass of pressurant required, total energy transfer, and wall and ullage temperatures. The analysis, the algorithm, a complete description of input and output, and the FORTRAN 4 program listings are presented. Sample cases are included to illustrate use of the programs.

Masters, P. A.↗

NASA NTRS1974

Interplanetary MeV electrons of Jovian origin

Observations of low energy electron increases observed in interplanetary space on Pioneer 10 are reported as it approached Jupiter. These discrete bursts were several hundred times the normal quiet-time electron flux, and became more frequent as one approached Jupiter resulting in the quasi-continuous presence of large fluxes of these electrons in interplanetary space. It is noted that the integrated flux from quiet-time electrons is comparable to the integrated ambient electron flux itself. In addition, the spectrum of electrons observed in Jupiter's magnetosphere, on Pioneer 10 in interplanetary space near Jupiter, for the quiet-time increases near the earth, and for the ambient electron spectrum are all remarkably similar. These two lines of evidence suggest the possibility that Jupiter could be the source of most of the ambient electrons at low energies.

Teegarden, B. J.↗

NASA NTRS1981

Medium energy gamma ray astronomy with transpacific balloon flights

Transpacific balloon flights with the University of California, Riverside (UCR) double scatter telescope are discussed. With flight durations from 5 days up to perhaps 15 days the long observation times necessary for medium energy (1-30 MeV) gamma ray astronomy can be obtained. These flights would be made under the auspices of the Joint U.S.-Japan Balloon Flight Program at NASA. It is proposed that flights can provide at least 30 hours of observation time per flight for many discrete source candidates and 120 hours for detecting low intensity cosmic gamma ray bursts.

Zych, A. D.↗