Search NASA⌕ Search

SEARCH · Search NASA

Results for “Implicit integration”

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

LSENS, a general chemical kinetics and sensitivity analysis code for gas-phase reactions: User's guide

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS, are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include static system, steady, one-dimensional, inviscid flow, shock initiated reaction, and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method, which works efficiently for the extremes of very fast and very slow reaction, is used for solving the 'stiff' differential equation systems that arise in chemical kinetics. For static reactions, sensitivity coefficients of all dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters can be computed. This paper presents descriptions of the code and its usage, and includes several illustrative example problems.

Radhakrishnan, Krishnan↗
NASA NTRS1993

High accuracy computation of fluid-structure interaction in transonic cascades

A coupling strategy for simulating fluid-structure interaction phenomena is formulated and applied to the prediction of flutter in transonic cascades. The flow is governed by the Euler equations and discretized using a finite volume flux-splitting scheme. The structure is modeled using an isoparametric finite element formulation. The coupling strategy successfully reconciles these two formulations at the fluid-structure interface by enforcing both kinematic and kinetic boundary conditions. In particular, the conservation laws applicable to the combined fluid-structure system are preserved across the interface. Since the primary mechanism driving aeroelastic phenomena involves energy exchange occurring at the interface, this highly accurate coupling mechanism is believed to improve the predictive capability of the scheme. The coupled equations are advanced simultaneously in time using an implicit time integration method. Results obtained using the coupling method are presented for cascade geometries operating in transonic flow.

Boschitsch, Alexander↗
NASA NTRS1992

Chemical vapor deposition fluid flow simulation modelling tool

Accurate numerical simulation of chemical vapor deposition (CVD) processes requires a general purpose computational fluid dynamics package combined with specialized capabilities for high temperature chemistry. In this report, we describe the implementation of these specialized capabilities in the spectral element code NEKTON. The thermal expansion of the gases involved is shown to be accurately approximated by the low Mach number perturbation expansion of the incompressible Navier-Stokes equations. The radiative heat transfer between multiple interacting radiating surfaces is shown to be tractable using the method of Gebhart. The disparate rates of reaction and diffusion in CVD processes are calculated via a point-implicit time integration scheme. We demonstrate the use above capabilities on prototypical CVD applications.

Bullister, Edward T.↗
NASA NTRS2000

LSENS, The NASA Lewis Kinetics and Sensitivity Analysis Code

A general chemical kinetics and sensitivity analysis code for complex, homogeneous, gas-phase reactions is described. The main features of the code, LSENS (the NASA Lewis kinetics and sensitivity analysis code), are its flexibility, efficiency and convenience in treating many different chemical reaction models. The models include: static system; steady, one-dimensional, inviscid flow; incident-shock initiated reaction in a shock tube; and a perfectly stirred reactor. In addition, equilibrium computations can be performed for several assigned states. An implicit numerical integration method (LSODE, the Livermore Solver for Ordinary Differential Equations), which works efficiently for the extremes of very fast and very slow reactions, is used to solve the "stiff" ordinary differential equation systems that arise in chemical kinetics. For static reactions, the code uses the decoupled direct method to calculate sensitivity coefficients of the dependent variables and their temporal derivatives with respect to the initial values of dependent variables and/or the rate coefficient parameters. Solution methods for the equilibrium and post-shock conditions and for perfectly stirred reactor problems are either adapted from or based on the procedures built into the NASA code CEA (Chemical Equilibrium and Applications).

Radhakrishnan, K.↗
NASA NTRS2007

One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure

The development and verification of a one-dimensional material thermal response code with ablation is presented. The implicit time integrator, control volume finite element spatial discretization, and Newton's method for nonlinear iteration on the entire system of residual equations have been implemented and verified for the thermochemical ablation of internally decomposing materials. This study is a continuation of the work presented in "One-Dimensional Ablation with Pyrolysis Gas Flow Using a Full Newton's Method and Finite Control Volume Procedure" (AIAA-2006-2910), which described the derivation, implementation, and verification of the constant density solid energy equation terms and boundary conditions. The present study extends the model to decomposing materials including decomposition kinetics, pyrolysis gas flow through the porous char layer, and a mixture (solid and gas) energy equation. Verification results are presented for the thermochemical ablation of a carbon-phenolic ablator which involves the solution of the entire system of governing equations.

Amar, Adam J.↗
NASA NTRS2011

Development and Verification of the Charring, Ablating Thermal Protection Implicit System Simulator

The development and verification of the Charring Ablating Thermal Protection Implicit System Solver (CATPISS) is presented. This work concentrates on the derivation and verification of the stationary grid terms in the equations that govern three-dimensional heat and mass transfer for charring thermal protection systems including pyrolysis gas flow through the porous char layer. The governing equations are discretized according to the Galerkin finite element method (FEM) with first and second order fully implicit time integrators. The governing equations are fully coupled and are solved in parallel via Newton s method, while the linear system is solved via the Generalized Minimum Residual method (GMRES). Verification results from exact solutions and Method of Manufactured Solutions (MMS) are presented to show spatial and temporal orders of accuracy as well as nonlinear convergence rates.

Amar, Adam J.↗
NASA NTRS2016

Simulation of an Electromechanical Spin Motor System of a Control Moment Gyroscope

A two-phase brushless DC motor (BDCM) with pulse-width modulated (PWM) voltage drive is simulated to control the flywheel speed of a control moment gyroscope (CMG). An overview of a double-gimballed control moment gyroscope (DGCMG) assembly is presented along with the CMG torque effects on the spacecraft. The operating principles of a two-phase brushless DC motor are presented and the system's electro-mechanical equations of motion are developed for the root-mean-square (RMS) currents and wheel speed. It is shown that the system is an extremely "stiff" set of first-order equations for which an implicit Euler integrator is required for a stable solution. An adaptive proportional voltage controller is presented which adjusts the PWM voltages depending on several control modes for speed, current, and torque. The simulation results illustrate the interaction between the electrical system and the load dynamics and how these influence the overall performance of the system. As will be shown, the CMG spin motor model can directly provide electrical power use and thermal power output to spacecraft subsystems for effective (average) calculations of CMG power consumption.

Inampudi, Ravi↗
NASA NTRS2015

Development of the US3D Code for Advanced Compressible and Reacting Flow Simulations

Aerothermodynamics and hypersonic flows involve complex multi-disciplinary physics, including finite-rate gas-phase kinetics, finite-rate internal energy relaxation, gas-surface interactions with finite-rate oxidation and sublimation, transition to turbulence, large-scale unsteadiness, shock-boundary layer interactions, fluid-structure interactions, and thermal protection system ablation and thermal response. Many of the flows have a large range of length and time scales, requiring large computational grids, implicit time integration, and large solution run times. The University of Minnesota NASA US3D code was designed for the simulation of these complex, highly-coupled flows. It has many of the features of the well-established DPLR code, but uses unstructured grids and has many advanced numerical capabilities and physical models for multi-physics problems. The main capabilities of the code are described, the physical modeling approaches are discussed, the different types of numerical flux functions and time integration approaches are outlined, and the parallelization strategy is overviewed. Comparisons between US3D and the NASA DPLR code are presented, and several advanced simulations are presented to illustrate some of novel features of the code.

CFD↗
DOE OSTI2023

Simulating Catalysis with Realistic Pellet Geometries Using Mesoflow: A Case Study of Catalytic Propane Dehydrogenation

We present a case study of catalytic propane dehydrogenation with our open-source multiphysics solver, Mesoflow. The solver was developed to simulate reactive flow coupled to heterogeneous catalytic reactions and deactivation in the context of complex, mesoscale geometry. The method leverages cartesian block-structured adaptive mesh refinement to capture realistic catalyst microstructural features acquired directly from X-ray computed tomography data. A kinetic model for propane dehydrogenation and catalyst deactivation was developed based on temporal analysis of products (TAP) reactor experiments. The TAP reactor experiments allow for precise characterization of intrinsic kinetic reaction steps which are implemented into Mesoflow simulations to model the spatial and temporal evolution of reactants, products, and catalyst active sites. The short-term and long-term deactivation behavior is studied by using XCT data collected from fresh and aged catalyst pellets, which exhibit different microstructural features. This study employs time-splitting algorithms to connect disparate reaction and flow timescales, enabling the simulations to achieve realistic deactivation timescales on the order of minutes while the flow time-scales for small particles (100 microns) are several milliseconds. We also introduce a flexible automated python script that writes the necessary files to construct a Mesoflow simulation from user-created chemical mechanisms. We will also introduce a few new features that are added to Mesoflow such as higher order schemes, implicit chemistry integrators and the ability to run on AMD and NVIDIA graphics-processing-units.

AMReX↗
NASA NTRS1977

An ILLIAC program for the numerical simulation of homogeneous incompressible turbulence

An algorithm and ILLIAC computer program, developed for the simulation of homogeneous incompressible turbulence in the presence of an applied mean strain, are described. The turbulence field is represented spatially by a truncated triple Fourier series (spectral method) and followed in time using a fourth-order Runge-Kutta algorithm. These include: (1) transformation of variables suggested by Taylor's sudden-distortion theory; (2) implicit viscous diffusion by use of an integrating factor; (3) implicit pressure calculation suggested by Taylor's sudden-distortion theory, and (4) inexpensive control of aliasing by random and phased coordinate shifts.

Rogallo, R. S.↗
DOE OSTI2025

Moment-based adaptive time integration for thermal radiation transport

Here, in this paper we develop a framework for moment-based adaptive time integration of deterministic multifrequency thermal radiation transpot (TRT). We generalize our recent semi-implicit-explicit (IMEX) integration framework for gray TRT to multifrequency TRT, and also introduce a semi-implicit variation that facilitates higher-order integration of TRT, where each stage is implicit in all components except opacities. To appeal to the broad literature on adaptivity with Runge–Kutta methods, we derive new embedded methods for four asymptotic preserving IMEX Runge–Kutta schemes we have found to be robust in our previous work on TRT and radiation hydrodynamics. We then use a moment-based high-order-low-order representation of the transport equations. Due to the high dimensionality, memory is always a concern in simulating TRT. We form error estimates and adaptivity in time purely based on temperature and radiation energy, for a trivial overhead in computational cost and memory usage compared with the base second order integrators. We then test the adaptivity in time on the tophat and Larsen problem, demonstrating the ability of the adaptive algorithm to naturally vary the timestep across 4–5 orders of magnitude, ranging from the dynamical timescales of the streaming regime to the thick diffusion limit.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Semi-implicit transient analysis procedures for structural dynamics analysis

A semi-implicit direct time integration procedure is presented which avoids factorization of the implicit difference solution matrix. The procedure, if properly implemented, requires only vectorial calculations and hence needs the same computer core space as explicit integration procedures. Guidelines for splitting the stiffness matrix into upper and lower matrices are established, which among other things are designed to satisfy a correct transmission of rigid-body motions from element (or grid) to its adjacent elements.

Park, K. C.↗
NASA NTRS1978

On a new algorithm for time step integration of nonlinear systems

A new implicit algorithm for time step integration of finite element structural dynamic equations is presented. Convergence, stability and numerical damping properties are discussed. Due to the way nonlinear structural behavior is taken into account, the algorithm is expected to compare favorably with existing ones. Some simple numerical results are presented. A related explicit algorithm is also derived and shortly discussed.

Anderheggen, E.↗
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 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 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 NTRS1990

Extended implicit model following as applied to integrated flight and propulsion control

An extended model-following control-synthesis method, including loop-transfer recovery, is presented and applied to synthesize control laws for integrated flight and propulsion control (IFPC). The vehicle considered is representative of an unstable modern fighter aircraft; with a two-dimensional thrust-vectoring and thrust-reversing nozzle. The linearized design model includes both airframe and engine dynamics. A model-following formulation of the LQR problem is extended to handle this hybrid problem. Compensators are then obtained to realize an output-feedback control law, by using a loop-transfer-recovery procedure. The airframe and engine responses are decoupled, and perfect airframe response following is obtained. The loop transfers also reveal good stability robustness and reasonable loop cross-over frequencies that would not lead to excessive actuation requirements.

Schmidt, David K.↗