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.

SourceResearch connectionFrom yearThrough year

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

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Combustion Feature Characterization using Computer Vision Diagnostics within Rotating Detonation Combustors

Rotating detonation engines (RDEs) theoretically achieve pressure gain by continuously propagating detonation waves in a cylindrical annulus. Research aims to implement them in gas turbines for improved propulsion and power generation efficiency. Current diagnostic methods, such as point measurements and optical diagnostics, face limitations due to high data acquisition rates required to analyze rapid detonation interactions. In contrast, image classification and time series classification achieve real-time capabilities with rates of 4 and 5 Hz, respectively. Object detection offers the highest time-step resolution at 20 μsec, while real-time methods require longer resolutions of 80 msec, highlighting the advancements in diagnostic capabilities for RDEs.

May, Kristyn Johnson↗

DOE OSTI2026

TANTE: Time-adaptive operator learning via neural Taylor expansion

Operator learning for time-dependent partial differential equations (PDEs) has seen rapid progress in recent years, enabling efficient approximation of complex spatiotemporal dynamics. However, most existing methods rely on fixed time step sizes during rollout, which limits their ability to adapt to varying temporal complexity and often leads to error accumulation. In this work, we propose the Time-Adaptive Transformer with Neural Taylor Expansion (TANTE), a novel operator-learning framework that produces continuous-time predictions with adaptive step sizes. TANTE predicts future states by performing a Taylor expansion at the current state, where neural networks learn both the higher-order temporal derivatives and the local radius of convergence. This allows the model to dynamically adjust its rollout based on the local behavior of the solution, thereby reducing cumulative error and improving computational efficiency. We demonstrate the effectiveness of TANTE across a wide range of PDE benchmarks, achieving superior accuracy and adaptability compared to fixed-step baselines, delivering accuracy gains of 60-80 % and speed-ups of 30-40 % at inference time.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1978

Comparative study of the convergence rates of two numerical techniques

The paper examines the applicability of the three-step Stetter (1968) method to the problem of hypersonic viscous flow over a blunt axisymmetric body used for planetary entry probes at zero angle of attack. The flow-field results using the two-step finite-difference MacCormack (1969) method are reported by Kumar and Graves (1977). Only the computational efficiency of Stetter's method is compared with that of MacCormack's in terms of the iterative time steps and computing time required for the steady-state solution. Advantages of Stetter's method over MacCormack's are established.

Kumar, A.↗

NASA NTRS1980

Direct integration of transient rotor dynamics

An implicit method was developed for integrating the equations of motion for a lumped mass model of a rotor dynamics system. As an aside, a closed form solution to the short bearing theory was also developed for a damper with arbitrary motion. The major conclusions are that the method is numerically stable and that the computation time is proportional to the number of elements in the rotor dynamics model rather than to the cube of the number. This computer code allowed the simulation of a complex rotor bearing system experiencing nonlinear transient motion and displayed the vast amount of results in an easily understood motion picture format - a 10 minute, 16 millimeter, color, sound motion picture supplement. An example problem with 19 mass elements in the rotor dynamics model took 0.7 second of central processing unit time per time step on an IBM 360-67 computer in a time sharing mode.

Kascak, A. F.↗

NASA NTRS1989

Unsteady Euler algorithm with unstructured dynamic mesh for complex-aircraft aeroelastic analysis

A finite-volume unstructured-grid FEM scheme with multistage Runge-Kutta time stepping is applied to the three-dimensional time-dependent Euler equations for inviscid flows on complex aircraft configurations undergoing structural deformation. The derivation of the model, the solution procedure, and the computer implementation are described, and results are presented graphically for a NASA Langley supersonic fighter aircraft model in steady and unsteady (harmonic oscillation in complete-vehicle bending mode) flow regimes. Good agreement between FEM predictions and experimental data is demonstrated.

Batina, John T.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

Kleb, William L.↗

NASA NTRS2009

Local-in-Time Adjoint-Based Method for Optimal Control/Design Optimization of Unsteady Compressible Flows

.We study local-in-time adjoint-based methods for minimization of ow matching functionals subject to the 2-D unsteady compressible Euler equations. The key idea of the local-in-time method is to construct a very accurate approximation of the global-in-time adjoint equations and the corresponding sensitivity derivative by using only local information available on each time subinterval. In contrast to conventional time-dependent adjoint-based optimization methods which require backward-in-time integration of the adjoint equations over the entire time interval, the local-in-time method solves local adjoint equations sequentially over each time subinterval. Since each subinterval contains relatively few time steps, the storage cost of the local-in-time method is much lower than that of the global adjoint formulation, thus making the time-dependent optimization feasible for practical applications. The paper presents a detailed comparison of the local- and global-in-time adjoint-based methods for minimization of a tracking functional governed by the Euler equations describing the ow around a circular bump. Our numerical results show that the local-in-time method converges to the same optimal solution obtained with the global counterpart, while drastically reducing the memory cost as compared to the global-in-time adjoint formulation.

Yamaleev, N. K.↗

NASA NTRS2014

Applying Motion Constraints Based on Test Data

MSC ADAMS is a simulation software that is used to analyze multibody dynamics. Using user subroutines, it is possible to apply motion constraints to the rigid bodies so that they match the motion profile collected from test data. This presentation describes the process of taking test data and passing it to ADAMS using user subroutines, and uses the Morpheus free-flight 4 test as an example of motion data used for this purpose. Morpheus is the name of a prototype lander vehicle built by NASA that serves as a test bed for various experimental technologies (see backup slides for details) MSC.ADAMS"TM" is used to play back telemetry data (vehicle orientation and position) from each test as the inputs to a 6-DoF general motion constraint (details in backup slides) The MSC.ADAMS"TM" playback simulations allow engineers to examine and analyze flight trajectory as well as observe vehicle motion from any angle and at any playback speed. This facilitates the development of robust and stable control algorithms, increasing reliability and reducing development costs of this developmental engine The simulation also incorporates a 3D model of the artificial hazard field, allowing engineers to visualize and measure performance of the developmental autonomous landing and hazard avoidance technology ADAMS is a multibody dynamics solver. It uses forces, constraints, and mass properties to numerically integrate equations of motion. The ADAMS solver will ask the motion subroutine for position, velocity, and acceleration values at various time steps. Those values must be continuous over the whole time domain. Each degree of freedom in the telemetry data can be examined separately; however, linear interpolation of the telemetry data is invalid, since there will be discontinuities in velocity and acceleration.

Burlone, Michael↗

NASA NTRS2016

Addressing Thermal Model Run Time Concerns of the Wide Field Infrared Survey Telescope using Astrophysics Focused Telescope Assets (WFIRST-AFTA)

The Wide Field Infrared Survey Telescope using Astrophysics Focused Telescope Assets (WFIRST-AFTA) utilizes an existing 2.4 m diameter Hubble sized telescope donated from elsewhere in the federal government for near-infrared sky surveys and Exoplanet searches to answer crucial questions about the universe and dark energy. The WFIRST design continues to increase in maturity, detail, and complexity with each design cycle leading to a Mission Concept Review and entrance to the Mission Formulation Phase. Each cycle has required a Structural-Thermal-Optical-Performance (STOP) analysis to ensure the design can meet the stringent pointing and stability requirements. As such, the models have also grown in size and complexity leading to increased model run time. This paper addresses efforts to reduce the run time while still maintaining sufficient accuracy for STOP analyses. A technique was developed to identify slews between observing orientations that were sufficiently different to warrant recalculation of the environmental fluxes to reduce the total number of radiation calculation points. The inclusion of a cryocooler fluid loop in the model also forced smaller time-steps than desired, which greatly increases the overall run time. The analysis of this fluid model required mitigation to drive the run time down by solving portions of the model at different time scales. Lastly, investigations were made into the impact of the removal of small radiation couplings on run time and accuracy. Use of these techniques allowed the models to produce meaningful results within reasonable run times to meet project schedule deadlines.

Thermal Analysis↗

NASA NTRS2019

Satellite-Based Assessment of Grassland Conversion and Related Fire Disturbance in the Kenai Peninsula, Alaska

Spruce beetle-induced (Dendroctonus rufipennis (Kirby)) mortality on the Kenai Peninsula has been hypothesized by local ecologists to result in the conversion of forest to grassland and subsequent increased fire danger. This hypothesis stands in contrast to empirical studies in the continental US which suggested that beetle mortality has only a negligible effect on fire danger. In response, we conducted a study using Landsat data and modeling techniques to map land cover change in the Kenai Peninsula and to integrate change maps with other geospatial data to predictively map fire danger for the same region. We collected Landsat imagery to map land cover change at roughly five-year intervals following a severe, mid-1990s beetle infestation to the present. Landcover classification was performed at each time step and used to quantify grassland encroachment patterns over time. The maps of land cover change along with digital elevation models (DEMs),temperature, and historical fire data were used to map and assess wildfire danger across the study area. Results indicate the highest wildfire danger tended to occur in herbaceous and black spruce land cover types, suggesting that the relationship between spruce beetle damage and wildfire danger in costal Alaskan forested ecosystems differs from the relationship between the two in the forests of the coterminous United States. These change detection analyses and fire danger predictions provide the Kenai National Wildlife Refuge (KENWR) ecologists and other forest managers a better understanding of the extent and magnitude of grassland conversion and subsequent change in fire danger following the 1990s spruce beetle outbreak.

wildfire↗

NASA NTRSDate not supplied

High-Fidelity Modeling Approach for Coupling Thermal, Fluid, and Neutronic Analysis in an NTP Engine

Nuclear thermal propulsion (NTP), with its increased efficiency compared to traditional rocket engines, is a potentially enabling technology for space missions of higher complexity, such as delivering heavy machinery or crew to the moon or other planets. Thus, exploring modeling approaches that can effectively couple the thermal, fluid, and neutronic behaviors of NTP reactors is essential to mitigate cost and schedule challenges with hardware development. This paper details one such approach which utilizes Cardinal to unify reactor physics analysis performed using OpenMC, heat transfer analysis performed using the MOOSE framework, and thermal hydraulic analysis performed using the MOOSE Thermal Hydraulic Module. This approach is demonstrated with a proof of concept that focuses on a legacy reactor geometry from the historic NERVA program. This preliminary model found an overprediction in, when compared to the reference data, axial power peaking and the temperatures of both solids and fluids by roughly 20%. It is clear, however, that the model conserved a number of broad trends, and a few model refinements should reduce the observed error. This approach has definite merit and should enable the production of informative results for more modern NTP systems. This work is a first step in creating a transient reactor model that can be updated across a sequence of time steps based on continually fluctuating, time-dependent boundary conditions.

Nuclear Thermal Propulsion↗

DOE OSTI2024

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

Many complex systems can be accurately modeled as a set of coupled time-dependent partial differential equations (PDEs). However, solving such equations can be prohibitively expensive, easily taxing the world’s largest supercomputers. One pragmatic strategy for attacking such problems is to split the PDEs into components that can more easily be solved in isolation. This operator splitting approach is used ubiquitously across scientific domains, and in many cases leads to a set of ordinary differential equations (ODEs) that need to be solved as part of a larger “outer-loop” time-stepping approach. The SUNDIALS library provides a plethora of robust time integration algorithms for solving ODEs, and the U.S. Department of Energy Exascale Computing Project (ECP) has supported its extension to applications on exascale-capable computing hardware. In this paper, we highlight some SUNDIALS capabilities and its deployment in combustion and cosmology application codes (Pele and Nyx, respectively) where operator splitting gives rise to numerous, small ODE systems that must be solved concurrently.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

Vectorization of a particle simulation method for hypersonic rarefied flow

An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.

Mcdonald, Jeffrey D.↗

NASA NTRS1989

A high-order Lagrangian-decoupling method for the incompressible Navier-Stokes equations

A high-order Lagrangian-decoupling method is presented for the unsteady convection-diffusion and incompressible Navier-Stokes equations. The method is based upon: (1) Lagrangian variational forms that reduce the convection-diffusion equation to a symmetric initial value problem; (2) implicit high-order backward-differentiation finite-difference schemes for integration along characteristics; (3) finite element or spectral element spatial discretizations; and (4) mesh-invariance procedures and high-order explicit time-stepping schemes for deducing function values at convected space-time points. The method improves upon previous finite element characteristic methods through the systematic and efficient extension to high order accuracy, and the introduction of a simple structure-preserving characteristic-foot calculation procedure which is readily implemented on modern architectures. The new method is significantly more efficient than explicit-convection schemes for the Navier-Stokes equations due to the decoupling of the convection and Stokes operators and the attendant increase in temporal stability. Numerous numerical examples are given for the convection-diffusion and Navier-Stokes equations for the particular case of a spectral element spatial discretization.

Ho, Lee-Wing↗

NASA NTRS1990

Numerical simulation of flow through biofluid devices

The results of a numerical simulation of flow through an artificial heart and through an artificial tilting-disk heart valve are presented. The simulation involves solving the incompressible Navier-Stokes equations; the solution process is described. The details and difficulties of modeling these particular geometries are discussed. The artificial heart geometry uses a single moving grid, and the valve computation uses an overlaid-grid approach with one moving grid and one stationary grid. The equations must be solved iteratively for each discrete time step of the computations, requiring a significant amount of computing time. It is particularly difficult to analyze and present the fluid physics represented by these calculations because of the time-varying nature of the flow, and because the flows are internal. Three-dimensional graphics and scientific visualization techniques have become instrumental in solving these problems.

Rogers, Stuart E.↗

NASA NTRS1991

Functional models of power electronic components for system studies

A novel approach to model power electronic circuits has been developed to facilitate simulation studies of system-level issues. The underlying concept for this approach is to develop an equivalent circuit, the functional model, that performs the same functions as the actual circuit but whose operation can be simulated by using larger time step size and the reduction in model complexity, the computation time required by a functional model is significantly shorter than that required by alternative approaches. The authors present this novel modeling approach and discuss the functional models of two major power electronic components, the DC/DC converter unit and the load converter, that are being considered by NASA for use in the Space Station Freedom electric power system. The validity of these models is established by comparing the simulation results with available experimental data and other simulation results obtained by using a more established modeling approach. The usefulness of this approach is demonstrated by incorporating these models into a power system model and simulating the system responses and interactions between components under various conditions.

Tam, Kwa-Sur↗

NASA NTRS1993

Tracking flow features using overset grids

A method is proposed to use overset grid topology to track dynamic flow features. Features of interest such as moving shock waves and vortices are overset with relatively fine tracker grids. Solutions are computed on the various grids and information is exchanged at intergrid boundaries. A grid track-sensor variable such as pressure is used to track the position of the flow feature to be resolved. The tracker grid is moved to the position where the track-sensor variable has the desired value (generally a maximum or a minimum) and new interpolation coefficients are computed for information exchange across grid boundaries. Solutions are computed at the current location and time-step, and grid motion is brought into the solution via time metrics. The method is demonstrated by tracking a moving shock and vortices shed behind a circular cylinder. It is conjectured that the method would show significant benefits in resolving features such as wakes behind oscillating airfoils and trajectories of jets issuing from rotating nozzles as encountered during thrust-vectoring.

Chawla, Kalpana↗