Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time step”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

A numerical method for the exact calculation of airloads associated with impulsively started wings

A numerical method is developed to calculate three-dimensional potential flows due to the steady and impulsive motion of isolated wing and wing-wing interaction problems. The velocity potential is represented by a discrete set of constant-doublet quadrilaterals on wing and wake surfaces. The exact surface boundary condition is enforced, and the solution is obtained in a step-by-step fashion, configurations being impulsively started from rest. Free-wake geometries are generated for each time step with Rankine or Lamb viscous vortex segments used in wake-velocity calculations. Sample results include calculated performance to steady state for a thick wing and indicial lift of a wing-wing interaction problem.

Summa, J. M.↗
NASA NTRS1984

Nonlinear evolution of the lower-hybrid drift instability

The results of simulations of the lower-hybrid drift instability in a neutral sheet configuration are described. The simulations use an implict formulation to relax the usual time step limitations and thus extend previous explicit calculations to weaker gradients, larger mass ratios, and long times compared with the linear growth time. The numerical results give the scaling of the saturation level, heating rates, resistivity, and cross-field diffusion and a demonstration by comparison with a fluid electron model that dissipation in the lower-hybrid drift instability is caused by electron kinetic effects.

Brackbill, J. U.↗
NASA NTRS1989

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially non-oscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1991

Solution of the hydrodynamic device model using high-order non-oscillatory shock capturing algorithms

A micron n+ - n - n+ silicon diode is simulated via the hydrodynamic model for carrier transport. The numerical algorithms employed are for the non-steady case, and a limiting process is used to reach steady state. The simulation employs shock capturing algorithms, and indeed shocks, or very rapid transition regimes, are observed in the transient case for the coupled system, consisting of the potential equation and the conservation equations describing charge, momentum, and energy transfer for the electron carriers. These algorithms, termed essentially nonoscillatory, were successfully applied in other contexts to model the flow in gas dynamics, magnetohydrodynamics, and other physical situations involving the conservation laws in fluid mechanics. The method here is first order in time, but the use of small time steps allows for good accuracy. Runge-Kutta methods allow one to achieve higher accuracy in time if desired. The spatial accuracy is of high order in regions of smoothness.

Fatemi, Emad↗
NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗
NASA NTRS2021

A Recursive Multi-step Machine Learning Approach for Airport Configuration Prediction

Airport configuration selection is a complex decision-making process that involves several operational and human factors. In this paper we propose a novel recursive multi-step machine learning (ML) approach to predict airport configuration. The multi-step approach guarantees stability of the predicted configuration by taking as input the configuration predicted at the previous time step. The features of the proposed model include weather data, future arrival and departure counts and current configuration. Due to the importance of arrival and departure counts in predicting the airport configuration, arrival counts are calculated using landing time predictions selected from physics-based landing time predictions available in FAA System Wide Information Management data feeds for each flight. The selection rules were developed and refined to select the most accurate time for different phases of flight. The proposed model predicts the airport configurations up to 6 hours ahead. In this paper we show the predictive performance of the proposed model for six major US airports, including Charlotte Douglas International Airport (CLT), Dallas/Fort Worth International Airport (DFW), John F. Kennedy International Airport (JFK), Newark Liberty International Airport (EWR), LaGuardia Airport (LGA) and Dallas Love Field Airport (DAL). We trained and evaluated models on 2019 and 2020 data in order to study the effect of the pandemic and how changes in traffic patterns affected the performance of the proposed model. Results are compared with a baseline assuming no airport configuration changes. In our results for DFW, we obtained a prediction accuracy of 89.3% for 3 hours ahead prediction, and 82.8% for 6 hours ahead when applied on 2019 data.

machine learning↗
NASA NTRS2021

A Recursive Multi-step Machine Learning Approach for Airport Configuration Prediction

Airport configuration selection is a complex decision-making process that involves several operational and human factors. In this paper we propose a novel recursive multi-step machine learning (ML) approach to predict airport configuration. The multi-step approach guarantees stability of the predicted configuration by taking as input the configuration predicted at the previous time step. The features of the proposed model include weather data, future arrival and departure counts and current configuration. Due to the importance of arrival and departure counts in predicting the airport configuration, arrival counts are calculated using landing time predictions selected from physics-based landing time predictions available in FAA System Wide Information Management data feeds for each flight. The selection rules were developed and refined to select the most accurate time for different phases of flight. The proposed model predicts the airport configurations up to 6 hours ahead. In this paper we show the predictive performance of the proposed model for six major US airports, including Charlotte Douglas International Airport (CLT), Dallas/Fort Worth International Airport (DFW), John F. Kennedy International Airport (JFK), Newark Liberty International Airport (EWR), LaGuardia Airport (LGA) and Dallas Love Field Airport (DAL). We trained and evaluated models on 2019 and 2020 data in order to study the effect of the pandemic and how changes in traffic patterns affected the performance of the proposed model. Results are compared with a baseline assuming no airport configuration changes. In our results for DFW, we obtained a prediction accuracy of 89.3% for 3 hours ahead prediction, and 82.8% for 6 hours ahead when applied on 2019 data.

machine learning↗
NASA NTRS1973

Error analysis of multi-conic techniques

A general error analysis of three recently developed multi-conic methods of three-body trajectory integration has been carried out. Single-step error functions for position and velocity have been derived as Taylor series in powers of the time step and also in integral form. These error functions are used to investigate the relative accuracy of the three methods in various regions of the earth-moon space and to provide a method of variable step size control for the trajectory integration procedure. Numerical results are used to compare the multi-step performance of the methods for both large and small step sizes.

D'Amario, L. A.↗
NASA NTRS1975

Shock effects from a large impact on the moon

The paper calculates the shock and shear deformation-induced internal energy distribution associated with a major basin-forming hypervelocity meteorite impact on the moon. The Hageman and Walsh formulation of the axisymmetric two-dimensional conservation equations in finite difference form is used, and the flow field induced upon impact of an iron meteorite traveling at 15 km/sec with a gabbroic anorthosite lunar crust is calculated for sequential time steps over a grid of hoop-shaped zones fixed in space. If an energy of 5 times 10 to the 32nd power ergs for the projectile energy is required to excavate a basin of Imbrium size, than a flow field of approximately 210 km in radius approximately 19 seconds after impact is indicated. The meteorite residue, the melt, and the ejecta are discussed, and the results are compared with the Gault and Heitowit (1963) formulation.

Okeefe, J. D.↗
NASA NTRS1988

Prediction of time-dependent fuselage pressures in the wake of a helicopter rotor

A method is described for the analysis of the unsteady, incompressible potential flow associated with a helicopter. The fuselage, rotor, and rotor wake are considered together in forward flight. This method is particularly useful in low advance ratio flight due to the contribution of the deformed wake to the airloads imposed on the fuselage. The rotor geometry is prescribed and the unsteady wake geometry is computed from the local flow perturbation velocities. The fuselage is modeled as a non-lifting body of source panels and the rotor and its wake are modeled as a full vortex lattice. The rotor geometry is arbitrary and several rotor blades can be represented. The unsteady airloads are computed in the presence of the deformed rotor wake by a time-stepping technique. The wake is started impulsively from rest, allowing a natural convection of the wake with time.

Berry, John D.↗
NASA NTRS1988

Grid embedding technique using Cartesian grids for Euler solutions

Grid-embedding techniques are presently applied to the solutions of the Euler equations on Cartesian grids for the NACA 0012 airfoil and the multielement SKF1.1 airfoil at transonic speeds. Based on comparisons with solutions on fine grids, it is shown that the present scheme, when used in conjunction with the Runge-Kutta time-stepping scheme as well as Cartesian grids, yields reduced memory requirements and a substantial reduction in computational time, without loss of accuracy.

Micheltree, R. A.↗
NASA NTRS1993

Unsteady three-dimensional thermal field prediction in turbine blades using nonlinear BEM

A time-and-space accurate and computationally efficient fully three dimensional unsteady temperature field analysis computer code has been developed for truly arbitrary configurations. It uses boundary element method (BEM) formulation based on an unsteady Green's function approach, multi-point Gaussian quadrature spatial integration on each panel, and a highly clustered time-step integration. The code accepts either temperatures or heat fluxes as boundary conditions that can vary in time on a point-by-point basis. Comparisons of the BEM numerical results and known analytical unsteady results for simple shapes demonstrate very high accuracy and reliability of the algorithm. An example of computed three dimensional temperature and heat flux fields in a realistically shaped internally cooled turbine blade is also discussed.

Martin, Thomas J.↗
NASA NTRS2004

MSFC Stream Model Preliminary Results: Modeling Recent Leonid and Perseid Encounters

The cometary meteoroid ejection model of Jones and Brown (1996b) was used to simulate ejection from comets 55P/Tempel-Tuttle during the last 12 revolutions, and the last 9 apparitions of 109P/Swift-Tuttle. Using cometary ephemerides generated by the Jet Propulsion Laboratory s (JPL) HORIZONS Solar System Data and Ephemeris Computation Service, two independent ejection schemes were simulated. In the first case, ejection was simulated in 1 hour time steps along the comet s orbit while it was within 2.5 AU of the Sun. In the second case, ejection was simulated to occur at the hour the comet reached perihelion. A 4th order variable step-size Runge-Kutta integrator was then used to integrate meteoroid position and velocity forward in time, accounting for the effects of radiation pressure, Poynting-Robertson drag, and the gravitational forces of the planets, which were computed using JPL s DE406 planetary ephemerides. An impact parameter was computed for each particle approaching the Earth to create a flux profile, and the results compared to observations of the 1998 and 1999 Leonid showers, and the 1993 and 2004 Perseids.

Cooke, William J.↗
NASA NTRS2014

Proto-Examples of Data Access and Visualization Components of a Potential Cloud-Based GEOSS-AI System

Once a research or application problem has been identified, one logical next step is to search for available relevant data products. Thus, an early component of a potential GEOSS-AI system, in the continuum between observations and end point research, applications, and decision making, would be one that enables transparent data discovery and access by users. Such a component might be effected via the systems data agents. Presumably, some kind of data cataloging has already been implemented, e.g., in the GEOSS Common Infrastructure (GCI). Both the agents and cataloging could also leverage existing resources external to the system. The system would have some means to accept and integrate user-contributed agents. The need or desirability for some data format internal to the system should be evaluated. Another early component would be one that facilitates browsing visualization of the data, as well as some basic analyses.Three ongoing projects at the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) provide possible proto-examples of potential data access and visualization components of a cloud-based GEOSS-AI system. 1. Reorganizing data archived as time-step arrays to point-time series (data rods), as well as leveraging the NASA Simple Subset Wizard (SSW), to significantly increase the number of data products available, at multiple NASA data centers, for production as on-the-fly (virtual) data rods. SSWs data discovery is based on OpenSearch. Both pre-generated and virtual data rods are accessible via Web services. 2. Developing Web Feature Services to publish the metadata, and expose the locations, of pre-generated and virtual data rods in the GEOSS Portal and enable direct access of the data via Web services. SSW is also leveraged to increase the availability of both NASA and non-NASA data.3.Federating NASA Giovanni (Geospatial Interactive Online Visualization and Analysis Interface), for multi-sensor data exploration, that would allow each cooperating data center, currently the NASA Distributed Active Archive Centers (DAACs), to configure its own Giovanni deployment, while also allowing all the deployments to incorporate each others data. A federated Giovanni comprises Giovanni Virtual Machines, which can be run on local servers or in the cloud.

access↗
NASA NTRS1975

DYNGEN: A program for calculating steady-state and transient performance of turbojet and turbofan engines

The DYNGEN, a digital computer program for analyzing the steady state and transient performance of turbojet and turbofan engines, is described. The DYNGEN is based on earlier computer codes (SMOTE, GENENG, and GENENG 2) which are capable of calculating the steady state performance of turbojet and turbofan engines at design and off-design operating conditions. The DYNGEN has the combined capabilities of GENENG and GENENG 2 for calculating steady state performance; to these the further capability for calculating transient performance was added. The DYNGEN can be used to analyze one- and two-spool turbojet engines or two- and three-spool turbofan engines without modification to the basic program. A modified Euler method is used by DYNGEN to solve the differential equations which model the dynamics of the engine. This new method frees the programmer from having to minimize the number of equations which require iterative solution. As a result, some of the approximations normally used in transient engine simulations can be eliminated. This tends to produce better agreement when answers are compared with those from purely steady state simulations. The modified Euler method also permits the user to specify large time steps (about 0.10 sec) to be used in the solution of the differential equations. This saves computer execution time when long transients are run. Examples of the use of the program are included, and program results are compared with those from an existing hybrid-computer simulation of a two-spool turbofan.

Sellers, J. F.↗
NASA NTRS1975

Three-dimensional finite element analysis for high velocity impact

A finite element algorithm for solving unsteady, three-dimensional high velocity impact problems is presented. A computer program was developed based on the Eulerian hydroelasto-viscoplastic formulation and the utilization of the theorem of weak solutions. The equations solved consist of conservation of mass, momentum, and energy, equation of state, and appropriate constitutive equations. The solution technique is a time-dependent finite element analysis utilizing three-dimensional isoparametric elements, in conjunction with a generalized two-step time integration scheme. The developed code was demonstrated by solving one-dimensional as well as three-dimensional impact problems for both the inviscid hydrodynamic model and the hydroelasto-viscoplastic model.

Chan, S. T. K.↗
NASA NTRS1991

Steady and unsteady solutions of the incompressible Navier-Stokes equations

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The equations are solved with a line-relaxation scheme that allows the use of very large pseudotime steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. The steady-state solution of flow through a square duct with a 90-deg bend is computed, and the results are compared with experimental data. Good agreement is observed. Computations of unsteady flow over a circular cylinder are presented and compared to other experimental and computational results. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗
NASA NTRS1988

Finite volume computation of unsteady inviscid rotational transonic flows past airfoils in rigid body motion

Unsteady inviscid transonic flow over airfoils in arbitrary rigid body motion is analyzed numerically by solving the two-dimensional unsteady Euler equations in integral form using a finite volume scheme. The solution procedure is based on an explicit Runge-Kutta time-stepping scheme wherein the spatial terms are central-differenced and a combination of second- and fourth-differences in the flow variables are used to form the numerical dissipation terms to stabilize the scheme. Unsteady calculations are started from converged steady-state solutions as initial conditions. Nonreflective boundary conditions are imposed on the far-field boundaries. Results are presented and, where possible, validated against available numerical and experimental data for airfoils subjected to a step change in angle of attack, airfoils oscillating and plunging in transonic flow, and airfoils immersed in a time-varying free stream.

Damodaran, Murali↗