SEARCH · Engineering Papers

Results for “implicit time 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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2018

Solution of Ordinary Differential Equations in Gradient-Based Multidisciplinary Design Optimization

A gradient-based approach to multidisciplinary design optimization enables efficient scalability to large numbers of design variables. However, the need for derivatives causes difficulties when integrating ordinary differential equations (ODEs) in models. To simplify this, we propose the use of the general linear methods framework, which unifies all Runge-Kutta and linear multistep methods. This approach enables rapid implementation of integration methods without the need to differentiate each one, even in a gradient-based optimization context. We also develop a new parallel time integration algorithm that enables vectorization across time steps. We present a set of benchmarking results using a stiff ODE, a non-stiff nonlinear ODE, and an orbital dynamics ODE, and compare integration methods. In a modular gradient-based multidisciplinary design optimization context, we find that the new parallel time integration algorithm with high-order implicit methods, especially Gauss-Legendre collocation, is the best choice for a broad range of problems.

Hwang, John T.↗

DOE OSTI2024

A high-order finite difference method for moving immersed domain boundaries and material interfaces

Here, we present a high-order sharp treatment of immersed moving domain boundaries and material interfaces, and apply it to the advection-diffusion equation in two and three dimensions. The spatial discretization combines dimension-split finite difference schemes with an immersed boundary treatment based on a weighted least-squares reconstruction of the solution, providing stable discretizations with up to sixth order accuracy for diffusion terms and third order accuracy for advection terms. The temporal discretization relies on a novel strategy for maintaining high-order temporal accuracy in problems with moving boundaries that minimizes implementation complexity and allows arbitrary explicit or diagonally-implicit Runge-Kutta schemes. The approach is broadly compatible with popular PDE-specialized Runge-Kutta time integrators, including low-storage, strong stability preserving, and diagonally implicit schemes. Through numerical experiments we demonstrate that the full discretization maintains high-order spatial and temporal accuracy in the presence of complex 3D geometries and for a range of boundary conditions, including Dirichlet, Neumann, and flux conditions with large jumps in coefficients.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Axisymmetric gyrokinetic simulation of ASDEX-Upgrade scrape-off layer using a conservative implicit BGK collision operator

Collisions play an important role in turbulence and transport of fusion plasmas. For kinetic simulations, as the collisionality increases in the domain of interest, the size of the time step to resolve the collisional physics can become overly restrictive in an explicit time integration scheme, leading to high computational cost. With the aim of overcoming such restriction, we have implemented an implicit Bhatnagar–Gross–Krook (BGK) collision operator for use in the discontinuous Galerkin full-f gyrokinetic solver within the Gkeyll framework, which, when combined with Gkeyll's traditional explicit time integrator for collisionless advection, can significantly increase the time step in gyrokinetic simulations of highly collisional regimes. To ensure conservation of density, momentum, and energy, we utilize an iterative scheme to correct the discretized approximation to the equilibrium Maxwellian distribution to which the BGK collision operator relaxes. We have further generalized the BGK infrastructure, both the implicit scheme and the correction routine, to handle cross-species collisions. This improved implicit and conservative BGK operator is benchmarked against the more accurate but more computationally expensive Lenard–Bernstein–Dougherty (LBD) operator, which has been utilized in prior studies with Gkeyll. The implicit BGK operator enables 2D axisymmetric simulations of the ASDEX-Upgrade scrape-off layer to run 56 times faster to completion than the simulations with the LBD operator, because the BGK operator is more robust and converges at a lower resolution than is required by the LBD operator. Additionally, in this more collisional limit, we demonstrate that the results of our simulations utilizing the implicit BGK operator agreed well with simulations utilizing the more computationally expensive LBD operator.

Gyrokinetic simulations↗

NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗

NASA NTRS1983

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗

DOE OSTI2024

A collision-based hybrid method for the BGK equation

In this article, we apply the collision-based hybrid method introduced by Hauck and McClarren to the Boltzmann equation with the BGK operator and a hyperbolic scaling. An implicit treatment of the source term is used to handle stiffness associated with the BGK operator. Although it helps the numerical scheme become stable with a large time step size, it is still not obvious to achieve the desired order of accuracy due to the relationship between the size of the spatial cell and the mean free path. Without asymptotic preserving property, a very restricted grid size is required to resolve the mean free path, which is not practical. Our approaches are based on the noncollision-collision decomposition of the BGK equation. We introduce the arbitrary order of nodal discontinuous Galerkin (DG) discretization in space with a semi-implicit time-stepping method; we employ the backward Euler time integration for the uncollided equation and the 2nd order predictor-corrector scheme for the collided equation, i.e., both source terms in uncollided and collided equations are treated implicitly and only streaming term in the collided equation is solved explicitly. This improves the computational efficiency without the complexity of the numerical implementation. Numerical results are presented for various Knudsen numbers to present the effectiveness and accuracy of our hybrid method. Also, we compare the solutions of the hybrid and non-hybrid schemes.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

A Fourier collocation time domain method for numerically solving Maxwell's equations

A new method for solving Maxwell's equations in the time domain for arbitrary values of permittivity, conductivity, and permeability is presented. Spatial derivatives are found by a Fourier transform method and time integration is performed using a second order, semi-implicit procedure. Electric and magnetic fields are collocated on the same grid points, rather than on interleaved points, as in the Finite Difference Time Domain (FDTD) method. Numerical results are presented for the propagation of a 2-D Transverse Electromagnetic (TEM) mode out of a parallel plate waveguide and into a dielectric and conducting medium.

Shebalin, John V.↗

DOE OSTI2022

Exponential integrators for non-linear diffusion

The goal of this project is to compare the performance of exponential time integrators with traditional methods such as diagonally implicit Runge-Kutta methods in the context of solving the system of reduced magnetohydrodynamics (RMHD). In this report, we present initial results of a proof of concept study that shows that exponential integrators can be an efficient alternative to traditional integration schemes.

42 ENGINEERING↗

NASA NTRS1991

Combustion of liquid fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both sub- and super-critical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influences on the fluid transport, gas/liquid interface thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibits a significant variation near the critical burning pressure, mainly as a result of reduced mass-diffusion rate and latent heat of vaporization with increased pressure. The influence of droplet size on the burning characteristics is also noted.

Shuen, J. S.↗

NASA NTRS1991

Vaporization and combustion of fuel droplets at supercritical conditions

Vaporization and combustion liquid-fuel droplets in both sub- and super-critical environments have been examined. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates finite-rate chemical kinetics and a full treatment of liquid-vapor phase equilibrium at the droplet surface. The governing equations and the associated interface boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to the supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures of 5-140 atm. In addition, the dynamic responses of droplet vaporization and combustion to ambient-pressure oscillations are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

NASA NTRS1992

Combustion of liquid-fuel droplets in supercritical conditions

A comprehensive analysis of liquid-fuel droplet combustion in both subcritical and supercritical environments has been conducted. The formulation is based on the complete conservation equations for both gas and liquid phases, and accommodates variable thermophysical properties, finite-rate chemical kinetics, and a full treatment of liquid-vapor phase equilibrium at the drop surface. The governing equations and associated interfacial boundary conditions are solved numerically using a fully coupled, implicit scheme with the dual time-stepping integration technique. The model is capable of treating the entire droplet history, including the transition from the subcritical to supercritical state. As a specific example, the combustion of n-pentane fuel droplets in air is studied for pressures in the range of 5-140 atm. Results indicate that the ambient gas pressure exerts significant control of droplet gasification and burning processes through its influence on fluid transport, gas-liquid interfacial thermodynamics, and chemical reactions. The droplet gasification rate increases progressively with pressure. However, the data for the overall burnout time exhibit a considerable change in the combustion mechanism at the critical pressure, mainly as a result of reduced mass diffusivity and latent heat of vaporization with increased pressure.

Shuen, J. S.↗

NASA NTRS1995

A NASTRAN-based computer program for structural dynamic analysis of Horizontal Axis Wind Turbines

This paper describes a computer program developed for structural dynamic analysis of horizontal axis wind turbines (HAWT's). It is based on the finite element method through its reliance on NASTRAN for the development of mass, stiffness, and damping matrices of the tower end rotor, which are treated in NASTRAN as separate structures. The tower is modeled in a stationary frame and the rotor in one rotating at a constant angular velocity. The two structures are subsequently joined together (external to NASTRAN) using a time-dependent transformation consistent with the hub configuration. Aerodynamic loads are computed with an established flow model based on strip theory. Aeroelastic effects are included by incorporating the local velocity and twisting deformation of the blade in the load computation. The turbulent nature of the wind, both in space and time, is modeled by adding in stochastic wind increments. The resulting equations of motion are solved in the time domain using the implicit Newmark-Beta integrator. Preliminary comparisons with data from the Boeing/NASA MOD2 HAWT indicate that the code is capable of accurately and efficiently predicting the response of HAWT's driven by turbulent winds.

Lobitz, Don W.↗

NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗

NASA NTRS1989

A time accurate finite volume method for propulsion chamber flows

An implicit three-dimensional time-accurate method for propulsion chamber flows is proposed which uses line Gauss-Seidel relaxation and multiple axial sweeps for the convergence of each time step. The general time-integration algorithm employed includes such schemes as the Euler implicit method. The results of spatial and temporal accuracy tests reveal that Roe's (1981) flux difference splitting provides excellent tracking of acoustic wave speeds. In comparison with other methods, no low mean flow Mach number convergence limitation or Courant number stabilization restriction is observed.

Beddini, R. A.↗

DOE OSTI2020

Implicit fast sweeping method for hyperbolic systems of conservation laws

Implicit time-accurate methods are often used to integrate stiff problems where explicit schemes impose severe time step restrictions. This paper presents an efficient numerical framework based on the Fast Sweeping Method (FSM) for solving linear and nonlinear hyperbolic systems of conservation laws. The solution at each discrete location is computed by sweeping the numerical domain in several predetermined directions that follow the causality of the characteristic families. The use of a fractional step strategy eliminates the need for a solution selection criterion while one-sided stencils limit the number of sweeps to at most 2 d for d space dimensions. This work focuses on the first-order implicit upwind method since it constitutes the building block for high-order conservative schemes. For problems where the degree of stiffness evolves over time, implicit-explicit hybridization can be accomplished with the same algorithm by simply switching the stencil at each time level. As opposed to traditional implicit solvers, the sweeping method does not require a local time linearization of the fluxes thereby preserving the nonlinear stability properties of the original implicit scheme. It also avoids the large computational and memory requirements associated with solving large block-diagonal systems of equations. Here, a series of one- and two-dimensional test cases are presented for the inviscid Burgers' equation and the reactive Euler equations. The results indicate that the implicit FSM can allow a major reduction in the number of time steps even in the presence of discontinuous solution profiles.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRS1990

Characteristic-based algorithms for flows in thermo-chemical nonequilibrium

A generalized finite-rate chemistry algorithm with Steger-Warming, Van Leer, and Roe characteristic-based flux splittings is presented in three-dimensional generalized coordinates for the Navier-Stokes equations. Attention is placed on convergence to steady-state solutions with fully coupled chemistry. Time integration schemes including explicit m-stage Runge-Kutta, implicit approximate-factorization, relaxation and LU decomposition are investigated and compared in terms of residual reduction per unit of CPU time. Practical issues such as code vectorization and memory usage on modern supercomputers are discussed.

Walters, Robert W.↗

NASA NTRS1987

Numerical experiments on transition control in wall-bounded shear flows

Results are presented from a numerical simulation of transition control in plane channel and boundary layer flows. The analysis is based on a pseudo-spectral/finite difference semi-implicit solution procedure employed to numerically integrate the time-dependent, three-dimensional, incompressible Navier-Stokes equations in a doubly periodic domain. In the channel flow, it was found that the active periodic suction/blowing method was effective in controlling strongly three-dimensional disturbances. In the boundary layer, the preliminary analysis indicated that in the early stages, passive control by suction is as effective as active control to suppress instabilities. The current work is focused on a detailed comparison of active and passive control by suction/blowing in the boundary layer.

Biringen, S.↗

NASA NTRS2016

The Impacts of Dry Dynamic Cores on Asymmetric Hurricane Intensification

The fundamental pathways for tropical cyclone (TC) intensification are explored by considering axisymmetric and asymmetric impulsive thermal perturbations to balanced, TC-like vortices using the dynamic cores of three different nonlinear numerical models. Attempts at reproducing the results of previous work, which used the community WRF Model, revealed a discrepancy with the impacts of purely asymmetric thermal forcing. The current study finds that thermal asymmetries can have an important, largely positive role on the vortex intensification, whereas other studies find that asymmetric impacts are negligible. Analysis of the spectral energetics of each numerical model indicates that the vortex response to asymmetric thermal perturbations is significantly damped in WRF relative to the other models. Spectral kinetic energy budgets show that this anomalous damping is primarily due to the increased removal of kinetic energy from the vertical divergence of the vertical pressure flux, which is related to the flux of inertia-gravity wave energy. The increased kinetic energy in the other two models is shown to originate around the scales of the heating and propagate upscale with time from nonlinear effects. For very large thermal amplitudes (50 K), the anomalous removal of kinetic energy due to inertia-gravity wave activity is much smaller, resulting in good agreement between models. The results of this paper indicate that the numerical treatment of small-scale processes that project strongly onto inertia-gravity wave energy can lead to significant differences in asymmetric TC intensification. Sensitivity tests with different time integration schemes suggest that diffusion entering into the implicit solution procedure is partly responsible for the anomalous damping of energy.

Guimond, Stephen R.↗