Search NASA⌕ Search

SEARCH · Search NASA

Results for “explicit time stepping”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Extension of multigrid methodology to supersonic/hypersonic 3-D viscous flows

A multigrid acceleration technique developed for solving 3-D Navier-Stokes equations for subsonic/transonic flows was extended to supersonic/hypersonic flows. An explicit multistage Runge-Kutta type of time stepping scheme is used as the basic algorithm in conjunction with the multigrid scheme. Solutions were obtained for a blunt conical frustum at Mach 6 to demonstrate the applicability of the multigrid scheme to high speed flows. Computations were performed for a generic High Speed Civil Transport configuration designed to cruise at Mach 3. These solutions show both the efficiency and accuracy of the present scheme for computing high speed viscous flows over configurations of practical interest.

Vatsa, Veer N.↗
NASA NTRS1990

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
NASA NTRS1991

The CFL condition for spectral approximations to hyperbolic initial-boundary value problems

The stability of spectral approximations to scalar hyperbolic initial-boundary value problems with variable coefficients are studied. Time is discretized by explicit multi-level or Runge-Kutta methods of order less than or equal to 3 (forward Euler time differencing is included), and spatial discretizations are studied by spectral and pseudospectral approximations associated with the general family of Jacobi polynomials. It is proved that these fully explicit spectral approximations are stable provided their time-step, delta t, is restricted by the CFL-like condition, delta t less than Const. N(exp-2), where N equals the spatial number of degrees of freedom. We give two independent proofs of this result, depending on two different choices of approximate L(exp 2)-weighted norms. In both approaches, the proofs hinge on a certain inverse inequality interesting for its own sake. The result confirms the commonly held belief that the above CFL stability restriction, which is extensively used in practical implementations, guarantees the stability (and hence the convergence) of fully-explicit spectral approximations in the nonperiodic case.

Gottlieb, David↗
DOE OSTI2024

Pressure stability in explicitly coupled simulations of poromechanics with application to CO 2 sequestration

We study in detail the pressure stabilizing effects of the non-iterated fixed-stress splitting in poromechanical problems which are nearly undrained and incompressible. When applied in conjunction with a spatial discretization which does not satisfy the discrete inf–sup condition, namely a mixed piecewise linear–piecewise constant spatial discretization, the explicit fixed-stress scheme can have a pressure stabilizing effect in transient problems. This effect disappears, however, upon time step refinement or the attainment of steady state. The interpretation of the scheme as an Augmented Lagrangian method similar to Uzawa iteration for incompressible flow helps explain these results. Moreover, due to the slowly evolving solution within undrained seal regions, we show that the explicit fixed-stress scheme requires very large time steps to reveal its pressure stabilizing effect in examples of geologic CO 2 sequestration. We note that large time steps can result in large errors in drained regions, such as the aquifer or reservoir regions of these examples, and can prevent convergence of nonlinear solvers in the case of multiphase flows, which can make the explicit scheme an unreliable source of pressure stabilization. We conclude by demonstrating that pressure jump stabilization is as effective in the explicit fixed-stress setting as in the fully implicit setting for undrained problems, while maintaining the stability and convergence of the fixed-stress split for drained problems.

58 GEOSCIENCES↗
NASA NTRS1989

The 3-D Euler and Navier-Stokes calculations for aircraft components

An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Vatsa, Veer N.↗
NASA NTRS1988

An explicit Runge-Kutta method for 3D turbulent incompressible flows

A computer code has been developed to solve for the steady-state solution of the 3D incompressible Reynolds-averaged Navier-Stokes equations. The approach is based on the cell-center, central-difference, finite-volume formulation and an explicit one-step, multistage Runge-Kutta time-stepping scheme. The Baldwin-Lomax turbulence model is used. Techniques to accelerate the rate of convergence to a steady-state solution include the preconditioned method, the local time stepping, and the implicit residual smoothing. Improvements in computational efficiency have been demonstrated in several areas. This numerical procedure has been used to simulate the turbulent horseshoe vortex flow around an airfoil/flat-plate juncture.

Sung, Chao-Ho↗
NASA NTRS1989

Development of an efficient multigrid code for 3-D Navier-Stokes equations

A multigrid acceleration technique has been developed to solve the three-dimensional Navier-Stokes equations efficiently. An explicit multistage Runge-Kutta type of time-stepping scheme is used as the basic algorithm in conjunction with the multigrid scheme. Solutions for flow over a finite wing have been obtained on extremely fine meshes in order to achieve grid convergence of the solutions. Present solutions indicate that the number of multigrid cycles required to achieve a given level of convergence does not increase with the number of mesh points employed, making it a very attractive scheme for fine meshes.

Vatsa, Veer N.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗
NASA NTRS1992

Evaluation of a multigrid-based Navier-Stokes solver for aerothermodynamic computations

A multigrid acceleration technique developed for solving the three-dimensional Navier-Stokes equations is used for computing high Mach number flows over configurations of practical interest. An explicit multistage Runge-Kutta type of time-stepping scheme is used as the basic algorithm. Solutions are presented for a spherically blunted cone at Mach 10 and a modified shuttle orbiter at Mach 6. The computed surface heat-transfer distributions are shown to compare favorably with the experimental data. Effect of grid-refinement on computed heat-transfer distributions is also examined to assess the numerical accuracy of the computed solutions. The rapid convergence rate associated with multigrid schemes in previous applications at transonic speeds is observed at the higher Mach number flows investigated here.

Vatsa, Veer N.↗
NASA NTRSDate not supplied

CFD Modeling & Validation of Single Bubble Growth During Cavitation in Zero Gravity

Cryogenic Fluid Management (CFM) in zero gravity represents a vital area of research as NASA pursues future long-duration space missions. Success of these missions relies heavily on the ability to store and transfer the cryogenic liquid fuels efficiently and safely. In the Zero Boil-Off Tank (ZBOT-1) Experiment in 2018, unexpected and intense phase change occurred at the screened Liquid Acquisition Device (LAD) during tank pressure control that was later attributed to cavitation. The cavitation ultimately resulted in pump failure. Since cryogen cavitation during storage and transfer of propellants in microgravity has been a source of serious concern for the aerospace community, a series of drop tower experiments with liquid methane were also performed at the University of Bremen to carefully study cavitation from a single bubble nucleus in an engineered micron-sized wall cavity. As the test cell is depressurized during the 10-second drop, cavitation causes bubble growth from the nucleus. Here, in a collaboration between the Case Western Reserve University and University of Bremen, we present a two-phase Computational Fluid Dynamics (CFD) model of bubble growth from a single wall cavity during depressurization in microgravity. Both 2D axisymmetric and 3D parametric simulation case studies are performed to show the evolution of the bubble growth. The model uses a Volume of Fluid (VOF) approach with and without adaptive mesh refinement and adaptive time stepping to capture the growth dynamics using an explicit formulation with algebraic reconstruction of the interface. CFD results show good agreement with corresponding experiment. Based on the velocity, temperature, and volume fraction fields generated by the validation case studies, detailed discussion and analysis of fluid flow and heat transfer during bubble growth process are presented to delineate the intricacies of the cavitation process. Strengths and shortcomings of the present numerical approach in capturing this phenomenon are also highlighted.

Two Phase Flow↗
NASA NTRS1987

Parallel processors and nonlinear structural dynamics algorithms and software

An explicit-explicit subcycling procedure for the finite element analysis of structural dynamics is developed. This procedure has relaxed the usual constraint of requiring integer time step ratios for adjacent nodal groups. This allows for greater advantage to be taken of local stability criteria, and thus improves the efficiency of the explicit time integrator. Example problems are included to demonstrate the accuracy and stability of the method.

Belytschko, Ted↗
NASA NTRS1990

Tailoring explicit time-marching schemes to improve convergence characteristics

Multi-stage time-stepping schemes, tailored to chosen spatial-differencing operators, are derived and tested. The schemes are constructed to give optimal damping of the high-frequency waves. They are ideal for use with multi-grid acceleration. The concept of characteristic time-stepping, necessary for the extension of the scalar analysis to systems of equations, is presented. The schemes show a marked improvement over Runge-Kutta schemes.

Powell, Kenneth G.↗
NASA NTRS1984

Stability and accuracy analysis of some fully-discrete algorithms for the one-dimensional second-order wave equation

The present investigation is concerned with some basic results for a predictor-multicorrector algorithm applied to the one-dimensional wave equation, giving particular attention to so-called 2-pass explicit schemes in which both lumped and coupled mass matrices are employed. In an assessment of the accuracy and stability properties of the algorithms, use is made of the one-dimensional, second-order wave equation. The maximum stable time step of the lumped right-hand-side mass, 2-pass explicit algorithm is twice that of the 1-pass explicit algorithm. Improved accuracy is obtained by employing higher-order, or consistent, right-hand-side, mass.

Hughes, T. J. R.↗
NASA NTRS1993

Application of an unstructured grid flow solver to planes, trains and automobiles

Rampant, an unstructured flow solver developed at Fluent Inc., is used to compute three-dimensional, viscous, turbulent, compressible flow fields within complex solution domains. Rampant is an explicit, finite-volume flow solver capable of computing flow fields using either triangular (2d) or tetrahedral (3d) unstructured grids. Local time stepping, implicit residual smoothing, and multigrid techniques are used to accelerate the convergence of the explicit scheme. The paper describes the Rampant flow solver and presents flow field solutions about a plane, train, and automobile.

Spragle, Gregory S.↗
DOE OSTI2025

Toward quantum tunneling from excited states: Recovering imaginary-time instantons from a real-time analysis

We revisit the path integral description of quantum tunneling and its generalization to excited states. For clarity, we focus on the simple toy model of a point particle in a double-well potential, for which we perform all steps explicitly. Instead of performing the familiar Wick rotation from physical to imaginary time—which is inconsistent with the requisite boundary conditions when treating tunneling from excited states—we regularize the path integral by adding an infinitesimal complex contribution to the Hamiltonian, while keeping time strictly real. We find that this gives rise to a complex stationary-phase solution, in agreement with recent insights from Picard-Lefshetz theory. We then show that there exists a class of analytic solutions for the corresponding equations of motion, which can be made to match the appropriate boundary conditions in the physically relevant limits of a vanishing regulator and an infinite physical time. We provide a detailed discussion of this nontrivial limit. We find that, for systems without an explicit time-dependence, our approach reproduces the picture of an instantonlike solution defined on a finite Euclidean-time interval. Lastly, we discuss the generalization of our approach to broader classes of systems, for which it serves as a reliable framework for high-precision calculations. Published by the American Physical Society 2025

Steingasser, Thomas (ORCID:0000000217262117)↗
NASA NTRS1996

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗
NASA NTRS1988

A Navier-Stokes solver for cascade flows

A computer code for solving the Reynolds averaged full Navier-Stokes equations has been developed and applied using sheared H-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady state calculations. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications in less than 100 multigrid cycles on the finest mesh.

Arnone, A.↗