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.

SourceResearch connectionFrom yearThrough year

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

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Fourier-Legendre spectral methods for incompressible channel flow

An iterative collocation technique is described for modeling implicit viscosity in three-dimensional incompressible wall bounded shear flow. The viscosity can vary temporally and in the vertical direction. Channel flow is modeled with a Fourier-Legendre approximation and the mean streamwise advection is treated implicitly. Explicit terms are handled with an Adams-Bashforth method to increase the allowable time-step for calculation of the implicit terms. The algorithm is applied to low amplitude unstable waves in a plane Poiseuille flow at an Re of 7500. Comparisons are made between results using the Legendre method and with Chebyshev polynomials. Comparable accuracy is obtained for the perturbation kinetic energy predicted using both discretizations.

Zang, T. A.↗

NASA NTRS1991

Dynamical capture of field stars by interstellar clouds

Evidence suggesting the existence of old stars in young galactic clusters, has motivated an investigation of a model in which such stars are captured from the field during the gas cloud collapse phase of star formation. A hyperbolic field star cannot be gravitationally captured by a static cloud since the star's total energy is conserved. However, if the cloud's potential is an explicit function of time, the energy of a penetrating star can be decreased. As a first step in modeling the dynamical capture of field stars, capture by freely falling, homogeneous spherical cloud is considered. This simplified model is solved, and a prediction is obtained for number of captured field stars in terms of the cloud parameters and dispersion velocities of the field stars.

Whitman, P. G.↗

NASA NTRS1982

Ageostrophic winds in the severe strom environment

The period from 1200 GMT 10 April to 0000 GMT 11 April 1979, during which time several major tornadoes and severe thunderstorms, including the Wichita Falls tornado occurred was studied. A time adjusted, isentropic data set was used to analyze key parameters. Fourth order centered finite differences were used to compute the isallobaric, inertial advective, tendency, inertial advective geostrophic and ageostrophic winds. Explicit isentropic trajectories were computed through the isentropic, inviscid equations of motion using a 15 minute time step. Ageostrophic, geostrophic and total vertical motion fields were computed to judge the relative importance of ageostrophy in enhancing the vertical motion field. It is found that ageostrophy is symptomatic of those mass adjustments which take place during upper level jet streak propagation and can, in a favorable environment, act to increase and release potential instability over meso alpha time periods.

Moore, J. T.↗

NASA NTRS1984

Multiple steady states for characteristic initial value problems

The time dependent, isentropic, quasi-one-dimensional equations of gas dynamics and other model equations are considered under the constraint of characteristic boundary conditions. Analysis of the time evolution shows how different initial data may lead to different steady states and how seemingly anamolous behavior of the solution may be resolved. Numerical experimentation using time consistent explicit algorithms verifies the conclusions of the analysis. The use of implicit schemes with very large time steps leads to erroneous results.

Salas, M. D.↗

NASA NTRS1983

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.

Gunzburger, M. D.↗

NASA NTRS1986

Smoothing and the second law

The technique of obtaining second order, oscillation free, total variation diminishing (TVD), scalar difference schemes by adding a limited diffusion flux (smoothing) to a second order centered scheme is explored. It is shown that such schemes do not always converge to the correct physical answer. The approach presented here is to construct schemes that numerically satisfy the second law of thermodynamics on a cell by cell basis. Such schemes can only converge to the correct physical solution and in some cases can be shown to be TVD. An explicit scheme with this property and second order spatial accuracy was found to have an extremely restrictive time step limitation (Delta t less than Delta x squared). Switching to an implicit scheme removed the time step limitation.

Merriam, Marshal L.↗

NASA NTRS1992

Upwind scheme for solving the Euler equations on unstructured tetrahedral meshes

An upwind scheme is presented for solving the three-dimensional Euler equations on unstructured tetrahedral meshes. Spatial discretization is accomplished by a cell-centered finite-volume formulation using flux-difference splitting. Higher-order differences are formed by a multidimensional linear reconstruction process. The solution gradients required for the higher-order differenes are computed by a novel approach that yields highly resolved solutions in regions of smooth flow while avoiding oscillations across shocks without explicitly applying a limiter. Solutions are advanced in time by a three-stage Runge-Kutta time-stepping scheme with convergence accelerated to steady state by local time stepping and implicit residual smoothing. Transonic solutions are presented for two meshes around the ONERA M6 wing and demonstrate substantial accuracy and insensitivity to mesh size.

Frink, Neal T.↗

NASA NTRS1984

Spectral methods in time for hyperbolic equations

A pseudospectral numerical scheme for solving linear, periodic, hyperbolic problems is described. It has infinite accuracy both in time and in space. The high accuracy in time is achieved without increasing the computational work and memory space which is needed for a regular, one step explicit scheme. The algorithm is shown to be optimal in the sense that among all the explicit algorithms of a certain class it requires the least amount of work to achieve a certain given resolution. The class of algorithms referred to consists of all explicit schemes which may be represented as a polynomial in the spatial operator.

Tal-Ezer, H.↗

NASA NTRS1987

Smoothing and the second law

The technique of obtaining second-order oscillation-free total -variation-diminishing (TVD), scalar difference schemes by adding a limited diffusive flux ('smoothing') to a second-order centered scheme is explored. It is shown that such schemes do not always converge to the correct physical answer. The approach presented here is to construct schemes that numerically satisfy the second law of thermodynamics on a cell-by-cell basis. Such schemes can only converge to the correct physical solution and in some cases can be shown to be TVD. An explicit scheme with this property and second-order spatial accuracy was found to have extremely restrictive time-step limitation. Switching to an implicit scheme removed the time-step limitation.

Merriam, Marshal L.↗

NASA NTRS1990

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Mahalingam, S.↗

NASA NTRSDate not supplied

Implicit Preconditioning for Explicit Multigrid Solvers on Cut-Cell Cartesian Meshes

This work assesses the effectiveness of linearized implicit Euler preconditioning for multigrid solvers using an unpreconditioned, Jacobian-free Newton Krylov method to converge the linear system of equations. Multigrid convergence rates improve to approximately 0.75 across the cases tested including a Mach 2 supersonic wedge, transonic NACA 0012 airfoil, and ONERA M6 wing. While larger Krylov subspaces increase the convergence rate, they also increase the computational cost, such that 4-8 Krylov vectors often offers the fastest turnaround. Further reductions in computational cost are achieved with a sequential hybrid preconditioner that begins with the explicit multigrid solver before transitioning to the preconditioned algorithm later on. In addition, a novel implementation of dual time stepping is extended to include both common BDF methods as well as high-order implicit Runge-Kutta schemes. This particular formulation, which uses A −1 preconditioning, is amenable to matrix-free solvers, and the L-stable methods are especially suited for meshes with arbitrarily small cut-cells. Asymptotic order of convergence is demonstrated for BDF1, BDF2, SDIRK2, and 3rd-order Radau IIA time integration with unsteady 2D vortex simulations.

ARMD↗

NASA NTRS2000

A Finite-Element Approach for Modeling Inviscid and Viscous Compressible Flows using Prismatic Grids

The Galerkin finite-element method is used to solve the Euler and Navier-Stokes equations on prismatic meshes. It is shown that the prismatic grid is advantageous for correctly and efficiently capturing the boundary layers in high Reynolds number flows. It can be captured accurately because of the ability to cluster grid points normal to the body. The efficiency derives from the implicit treatment of the normal direction. To treat the normal direction implicitly, a semi-implicit Runge-Kutta time stepping scheme is developed. The semi-implicit algorithm is validated on simple geometries for inviscid and viscous flows and its convergence history is compared to that of the explicit Runge-Kutta scheme. The semi-implicit scheme is shown to be a factor of 3 to 4 faster in terms of CPU time to convergence.

Pandya, S. A.↗

NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗

NASA NTRS1981

Optimum design of hybrid phase locked loops

The design procedure of phase locked loops is described in which the analog loop filter is replaced by a digital computer. Specific design curves are given for the step and ramp input changes in phase. It is shown that the designed digital filter depends explicitly on the product of the sampling time and the noise bandwidth of the phase locked loop. This technique of optimization can be applied to the design of digital analog loops for other applications.

Lee, P.↗

NASA NTRS1989

Unsteady Navier-Stokes computations past oscillating delta wing at high incidence

The unsteady, thin-layer, compressible Navier-Stokes equations, written in the moving frame of reference for the flow relative motion, is solved for the steady and unsteady supersonic flow around a round-edged delta wing. For supersonic flow, local conical flow solution has been obtained from the three-dimensional equations. Pseudotime stepping is used for the steady flow, while time-accurate stepping is used for the unsteady flow. The computational scheme is an implicit approximately-factored finite volume scheme which uses explicit and implicit dissipation terms. The scheme is verified for the steady flow solution. The scheme is then applied to a delta wing undergoing rolling oscillation at a reduced frequency of 1.137 with 15- deg maximum amplitude about a mean angle of attack of 10 deg for a Mach number of 2 and a Reynolds number of 500,000.

Kandil, Osama A.↗

NASA NTRS1990

Mapping unstructured grid problems to the connection machine

We present a highly parallel graph mapping technique that enables one to solve unstructured grid problems on massively parallel computers. Many implicit and explicit methods for solving discretizated partial differential equations require each point in the discretization to exchange data with its neighboring points every time step or iteration. The time spent communicating can limit the high performance promised by massively parallel computing. To eliminate this bottleneck, we map the graph of the irregular problem to the graph representing the interconnection topology of the computer such that the sum of the distances that the messages travel is minimized. We show that, in comparison to a naive assignment of processors, our heuristic mapping algorithm significantly reduces the communication time on the Connection Machine, CM-2.

Hammond, Steven W.↗

NASA NTRS1990

Solving unstructured grid problems on massively parallel computers

A highly parallel graph mapping technique that enables one to efficiently solve unstructured grid problems on massively parallel computers is presented. Many implicit and explicit methods for solving discretized partial differential equations require each point in the discretization to exchange data with its neighboring points every time step or iteration. The cost of this communication can negate the high performance promised by massively parallel computing. To eliminate this bottleneck, the graph of the irregular problem is mapped into the graph representing the interconnection topology of the computer such that the sum of the distances that the messages travel is minimized. It is shown that using the heuristic mapping algorithm significantly reduces the communication time compared to a naive assignment of processes to processors.

Hammond, Steven W.↗

NASA NTRS1976

Multidimensional difference schemes with fourth-order accuracy

An explicit finite-difference algorithm is presented for the solution of quasilinear divergence free multidimensional hyperbolic systems. The method consists of four steps per time level. The resulting scheme is fourth-order accurate in both space and time, though the intermediate steps are only first-order accurate. The family of schemes introduced is dissipative, and hence, suitable for both smooth flows and flows containing shocks. This method is compared, in several numerical examples, with both second-order schemes and others that are fourth order in space, but second order in time.

Turkel, E.↗