Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A two-dimensional, time-dependent, near-earth magnetotail

A kinetic approach is presented for obtaining slowly time-dependent, self-consistent plasma, and magnetic field configurations for this region, using anisotropic particle distributions. At present, a 2D configuration is used. The time variation of the plasma/field configuration is assumed to be driven by a slowly varying, externally imposed electric field (representing, for example, a change in solar wind conditions). The numerical approach involves solving a 2D nonlinear Poisson equation for the magnetic vector potential at each time step. The local current density is calculated from particle velocity distributions under the assumption of adiabatic particle motion. The obtained equations of state are compared with those of double-adiabatic theory. Numerical examples of self-consistent responses of the near-earth plasma and field to an externally imposed electric field pulse are presented.

Whipple, Elden↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results are shown a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
NASA NTRS1993

Development of an explicit multiblock/multigrid flow solver for viscous flows in complex geometries

A new computer program is being developed for doing accurate simulations of compressible viscous flows in complex geometries. The code employs the full compressible Navier-Stokes equations. The eddy viscosity model of Baldwin and Lomax is used to model the effects of turbulence on the flow. A cell centered finite volume discretization is used for all terms in the governing equations. The Advection Upwind Splitting Method (AUSM) is used to compute the inviscid fluxes, while central differencing is used for the diffusive fluxes. A four-stage Runge-Kutta time integration scheme is used to march solutions to steady state, while convergence is enhanced by a multigrid scheme, local time-stepping, and implicit residual smoothing. To enable simulations of flows in complex geometries, the code uses composite structured grid systems where all grid lines are continuous at block boundaries (multiblock grids). Example results shown are a flow in a linear cascade, a flow around a circular pin extending between the main walls in a high aspect-ratio channel, and a flow of air in a radial turbine coolant passage.

Steinthorsson, E.↗
NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/sq cm for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/sq cm, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8-species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗
NASA NTRS2006

Viking Afterbody Heating Computations and Comparisons to Flight Data

Computational fluid dynamics predictions of Viking Lander 1 entry vehicle afterbody heating are compared to flight data. The analysis includes a derivation of heat flux from temperature data at two base cover locations, as well as a discussion of available reconstructed entry trajectories. Based on the raw temperature-time history data, convective heat flux is derived to be 0.63-1.10 W/cm2 for the aluminum base cover at the time of thermocouple failure. Peak heat flux at the fiberglass base cover thermocouple is estimated to be 0.54-0.76 W/cm2, occurring 16 seconds after peak stagnation point heat flux. Navier-Stokes computational solutions are obtained with two separate codes using an 8- species Mars gas model in chemical and thermal non-equilibrium. Flowfield solutions using local time-stepping did not result in converged heating at either thermocouple location. A global time-stepping approach improved the computational stability, but steady state heat flux was not reached for either base cover location. Both thermocouple locations lie within a separated flow region of the base cover that is likely unsteady. Heat flux computations averaged over the solution history are generally below the flight data and do not vary smoothly over time for both base cover locations. Possible reasons for the mismatch between flight data and flowfield solutions include underestimated conduction effects and limitations of the computational methods.

Edquist, Karl T.↗
NASA NTRS1993

Finite element/finite volume approaches with adaptive time stepping strategies for transient thermal problems

An adaptive time stepping strategy for transient thermal analysis of engineering systems is described which computes the time step based on the local truncation error with a good global error control and obtains optimal time steps to be used during the analysis. Combined mesh partitionings involving FEM/FVM meshes based on physical situations to obtain numerically improved physical representations are also proposed. Numerical test cases are described and comparative pros and cons are identified for practical situations.

Mohan, Ram V.↗
DOE OSTI2025

Real-time steerable frequency-stepped Doppler backscattering (DBS) system for local helicon wave electric field measurements on the DIII-D tokamak

A new frequency-stepped Doppler backscattering (DBS) system has been integrated into a real-time steerable electron cyclotron heating launcher system to simultaneously probe local background turbulence (f < 10 MHz) and high-frequency (20–550 MHz) density fluctuations in the DIII-D tokamak. The launcher allows for 2D steering (horizontally and vertically) over wide angular ranges to optimize probe location and wavenumber response. The vertical steering can be optimized during a discharge in real time. The new DBS system employs a programmable frequency synthesizer with adjustable dwell time as a source to launch either O or X-mode polarized millimeter waves. This system can step in real-time over the entire E-band frequency range (60–90 GHz). This combination of capabilities allows for the diagnosis of the complex internal spatial structure of high power (>200 kW) helicon waves (476 MHz) injected from an external antenna during helicon current drive experiments in DIII-D. Broadband density fluctuations around the helicon frequency are observed during real-time scans of measurement location and wavenumber during these experiments. Analysis indicates that these broadband high-frequency fluctuations are a result of backscattering of the DBS millimeter-wave probe beam from plasma turbulence modulated by the helicon wave. It is observed that background turbulence is effectively locally “tagged” with the helicon wave electric field, forming images of the turbulent spectrum in the overall density fluctuation spectrum that appear as high-frequency sidebands of the turbulence. These observations of background turbulence and high-frequency fluctuations open up the possibility of monitoring local helicon wave amplitude by comparing the high-frequency signal amplitude to the simultaneously measured background turbulence. In combination with the real-time measurement location and wavenumber scanning capabilities (offered by real-time frequency-stepping and steering), this allows rapid determination of the spatial distribution of the helicon wave power during steady-state plasma operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 NTRS2008

Oxidation Behavior of Carbon Fiber-Reinforced Composites

OXIMAP is a numerical (FEA-based) solution tool capable of calculating the carbon fiber and fiber coating oxidation patterns within any arbitrarily shaped carbon silicon carbide composite structure as a function of time, temperature, and the environmental oxygen partial pressure. The mathematical formulation is derived from the mechanics of the flow of ideal gases through a chemically reacting, porous solid. The result of the formulation is a set of two coupled, non-linear differential equations written in terms of the oxidant and oxide partial pressures. The differential equations are solved simultaneously to obtain the partial vapor pressures of the oxidant and oxides as a function of the spatial location and time. The local rate of carbon oxidation is determined at each time step using the map of the local oxidant partial vapor pressure along with the Arrhenius rate equation. The non-linear differential equations are cast into matrix equations by applying the Bubnov-Galerkin weighted residual finite element method, allowing for the solution of the differential equations numerically.

Sullivan, Roy M.↗
NASA NTRS1993

An adaptive mesh-moving and refinement procedure for one-dimensional conservation laws

We examine the performance of an adaptive mesh-moving and /or local mesh refinement procedure for the finite difference solution of one-dimensional hyperbolic systems of conservation laws. Adaptive motion of a base mesh is designed to isolate spatially distinct phenomena, and recursive local refinement of the time step and cells of the stationary or moving base mesh is performed in regions where a refinement indicator exceeds a prescribed tolerance. These adaptive procedures are incorporated into a computer code that includes a MacCormack finite difference scheme wih Davis' artificial viscosity model and a discretization error estimate based on Richardson's extrapolation. Experiments are conducted on three problems in order to qualify the advantages of adaptive techniques relative to uniform mesh computations and the relative benefits of mesh moving and refinement. Key results indicate that local mesh refinement, with and without mesh moving, can provide reliable solutions at much lower computational cost than possible on uniform meshes; that mesh motion can be used to improve the results of uniform mesh solutions for a modest computational effort; that the cost of managing the tree data structure associated with refinement is small; and that a combination of mesh motion and refinement reliably produces solutions for the least cost per unit accuracy.

Biswas, Rupak↗
NASA NTRS1981

Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes

A new combination of a finite volume discretization in conjunction with carefully designed dissipative terms of third order, and a Runge Kutta time stepping scheme, is shown to yield an effective method for solving the Euler equations in arbitrary geometric domains. The method has been used to determine the steady transonic flow past an airfoil using an O mesh. Convergence to a steady state is accelerated by the use of a variable time step determined by the local Courant member, and the introduction of a forcing term proportional to the difference between the local total enthalpy and its free stream value.

Jameson, A.↗
DOE OSTI2020

Grid generation, coupling strategies, and spatially-dependent time stepping for ocean-tidal/estuary systems and other ESM components (Final Report)

There are many multiple scale issues in the ocean modeling, including different mesh resolutions (local) and different physical processes (global). In order to gain efficiency and fidelity in modern ESMs, it is vital to develop high quality, highly nonuniform meshes and spatially-dependent time stepping strategies for coupling various climate system components. The USC/FSU team has worked on the design, analysis, and implementation of fast CVT grid generation and novel time- stepping schemes that could address all the requirements for global and regional ocean models on highly nonuniform meshes.

58 GEOSCIENCES↗
NASA NTRS2003

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗
DOE OSTI2020

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Localized Exponential Time Differencing Method for Shallow Water Equations: Algorithms and Numerical Study

Here, we explore the performance of the exponential time differencing (ETD) method applied to the rotating shallow water equations. Comparing with explicit time stepping of the same order accuracy in time, the ETD algorithms could reduce the computational time in many cases by allowing the use of large time step sizes while still maintaining numerical stability. To accelerate the ETD simulations, we propose a localized approach that synthesizes the ETD method and overlapping domain decomposition. By dividing the original problem into many subdomain problems of smaller sizes and solving them locally, the proposed approach could speed up the calculation of matrix exponential vector products. Several standard test cases for shallow water equations of one or multiple layers are considered. The results show great potential of the localized ETD method for high-performance computing because each subdomain problem can be naturally solved in parallel at every time step.

58 GEOSCIENCES↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

A semi-Lagrangian approach to the shallow water equation

We present a formulation of the shallow water equations that emphasizes the conservation of potential vorticity. A locally conservative semi-Lagrangian time-stepping scheme is developed, which leads to a system of three coupled PDE's to be solved at each time level. We describe a smoothing analysis of these equations, on which an effective multigrid solver is constructed. Some results from applying this solver to the static version of these equations are presented.

Bates, J. R.↗