SEARCH · Engineering Papers

Results for “time integrator”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

NASA NTRS1992

On computational methods for crashworthiness

The evolution of computational methods for crashworthiness and related fields is described and linked with the decreasing cost of computational resources and with improvements in computation methodologies. The latter includes more effective time integration procedures and more efficient elements. Some recent developments in methodologies and future trends are also summarized. These include multi-time step integration (or subcycling), further improvements in elements, adaptive meshes, and the exploitation of parallel computers.

Belytschko, T.↗

DOE OSTI2022

CondiDiag1.0: a flexible online diagnostic tool for conditional sampling and budget analysis in the E3SM atmosphere model (EAM)

Abstract. Numerical models used in weather and climate prediction take into account a comprehensive set of atmospheric processes (i.e., phenomena) such as the resolved and unresolved fluid dynamics, radiative transfer, cloud and aerosol life cycles, and mass or energy exchanges with the Earth's surface. In order to identify model deficiencies and improve predictive skills, it is important to obtain process-level understanding of the interactions between different processes. Conditional sampling and budget analysis are powerful tools for process-oriented model evaluation, but they often require tedious ad hoc coding and large amounts of instantaneous model output, resulting in inefficient use of human and computing resources. This paper presents an online diagnostic tool that addresses this challenge by monitoring model variables in a generic manner as they evolve within the time integration cycle. The tool is convenient to use. It allows users to select sampling conditions and specify monitored variables at run time. Both the evolving values of the model variables and their increments caused by different atmospheric processes can be monitored and archived. Online calculation of vertical integrals is also supported. Multiple sampling conditions can be monitored in a single simulation in combination with unconditional sampling. The paper explains in detail the design and implementation of the tool in the Energy Exascale Earth System Model (E3SM) version 1. The usage is demonstrated through three examples: a global budget analysis of dust aerosol mass concentration, a composite analysis of sea salt emission and its dependency on surface wind speed, and a conditionally sampled relative humidity budget. The tool is expected to be easily portable to closely related atmospheric models that use the same or similar data structures and time integration methods.

58 GEOSCIENCES↗

DOE OSTI2020

Evidence of restricted heat transport in National Ignition Facility Hohlraums

We present experimental evidence of restricted electron thermal conduction in the high-Z coronal plasma regions of laser-driven Hohlraums on the National Ignition Facility. Four separate measurements, three of which are direct observations of Hohlraum dynamics, corroborate this finding. (1) The velocity of the coronal plasma ablated and heated by the outer-cone laser beams is determined by time-dependent imaging of the gold plasma plume, or “bubble.” The velocities of the incoming plume (perpendicular to the Hohlraum axis) are consistent with high-fidelity 2D radiation-hydrodynamic simulations using flux-limited thermal electron conduction with a flux multiplier f = 0.03. Simulations using f = 0.15, which is very nearly classical Spitzer–Härm transport, predict plume velocities slower than measured. (2) Specific features in time-resolved images of the Hohlraum wall at an angle of 19° are also more consistent with f = 0.03 simulations compared to f = 0.15. (3) Spectroscopic tracers were added to the Hohlraum wall in the outer-beam bubble region. The ratios of hydrogen-like to helium-like line emission are sensitive to the electron temperature of the bubble. The hydrogen-like to helium-like ratios extracted from the time-integrated spectra of manganese and cobalt tracers from two observation angles are consistent with f = 0.03 and not with f = 0.15. (4) The time of peak capsule emission, or “bang time,” an integrated measurement, is also more consistent with f = 0.03 than with f = 0.15. Overall, while these findings do not identify the causes of restricted thermal conduction in Hohlraums, they motivate future experiments to test specific hypotheses and focus on model development in the regions of the plasma exhibiting restricted transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗

DOE OSTI2023

Closure of Earth’s Global Seasonal Cycle of Energy Storage

Abstract The global seasonal cycle of energy in Earth’s climate system is quantified using observations and reanalyses. After removing long-term trends, net energy entering and exiting the climate system at the top of the atmosphere (TOA) should agree with the sum of energy entering and exiting the ocean, atmosphere, land, and ice over the course of an average year. Achieving such a balanced budget with observations has been challenging. Disagreements have been attributed previously to sparse observations in the high-latitude oceans. However, limiting the local vertical integration of new global ocean heat content estimates to the depth to which seasonal heat energy is stored, rather than integrating to 2000 m everywhere as done previously, allows closure of the global seasonal energy budget within statistical uncertainties. The seasonal cycle of energy storage is largest in the ocean, peaking in April because ocean area is largest in the Southern Hemisphere and the ocean’s thermal inertia causes a lag with respect to the austral summer solstice. Seasonal cycles in energy storage in the atmosphere and land are smaller, but peak in July and September, respectively, because there is more land in the Northern Hemisphere, and the land has more thermal inertia than the atmosphere. Global seasonal energy storage by ice is small, so the atmosphere and land partially offset ocean energy storage in the global integral, with their sum matching time-integrated net global TOA energy fluxes over the seasonal cycle within uncertainties, and both peaking in April.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Dynamic Response of a Semiactive Suspension System with Hysteretic Nonlinear Energy Sink Based on Random Excitation by means of Computer Simulation

This paper aims to investigate the property and behavior of the hysteretic nonlinear energy sink (HNES) coupled to a half vehicle system which is a nine-degree-of-freedom, nonlinear, and semiactive suspension system in order to improve the ride comfort and increase the stability in shock mitigation by using the computer simulation method. The HNES model is a semiactive suspension device, which comprises the famous Bouc–Wen (B-W) model employed to describe the force produced by both the purely hysteretic spring and linear elastic spring of potentially negative stiffness connected in parallel, for the half vehicle system. Nine nonlinear motion equations of the half vehicle system are derived in terms of the seven displacements and the two dimensionless hysteretic variables, which are integrated numerically by employing the direct time integration method for studying both the variables of vertical displacements, velocities, accelerations, chassis pitch angle, and the ride comfort and driver safety, respectively, based on the bump and random road inputs of the pseudoexcitation method as excitation signal. Simulation results show that, compared with the HNES model and the magnetorheological (MR) model coupled to the half vehicle system, the ride comfort and stability have been evidently improved. A successful validation process has been performed, which indicated that both the ride comfort and driver safety properties of the HNES model coupled to half vehicle significantly improved.

Chen, Hui↗

NASA NTRS1972

Turbulent transport measurements with a laser Doppler velocimeter.

The power spectrum of phototube current from a laser Doppler velocimeter operating in the heterodyne mode has been computed. The spectral width and shape predicted by the theory are in agreement with experiment. For normal operating parameters the time-average spectrum contains information only for times shorter than the Lagrangian-integral time scale of the turbulence. To examine the long-time behavior, one must use either extremely small scattering angles, much-longer-wavelength radiation, or a different mode of signal analysis, e.g., FM detection.

Edwards, R. V.↗

NASA NTRS1980

The application of Seasat-1 radar altimetry to continental shelf circulation modeling

Seasat-1 radar altimetry data have been used to verify a continental shelf circulation model subjected to tidal forcing on the open boundary. The model makes use of the semi-implicit mode of time integration, removing the surface gravity wave time step restriction. Both the altimetry and the model predictions are compared with empirically determined tidal fluctuations and generally good agreement is obtained.

Cornillon, P.↗

NASA NTRS1992

Computational modeling of femtosecond optical solitons from Maxwell's equations

An algorithm is developed that permits the direct time integration of full-vector nonlinear Maxwell's equations. This capability permits the modeling of both linear and nonlinear instantaneous and dispersive effects in the electric polarization in material media. The modeling of the optical carrier is retained. The fundamental innovation is to notice that it is possible to treat the linear and nonlinear convolution integrals, which describe the dispersion, as new dependent variables. A coupled system of nonlinear second-order ordinary differential equations can then be derived for the linear and nonlinear convolution integrals, by differentiating them in the time domain. These equations, together with Maxwell's equations, are solved to determine the electromagnetic fields in nonlinear dispersive media. Results are presented of calculations in one dimension of the propagation and collision of femtosecond electromagnetic solitons that retain the optical carrier, taking into account as the Kerr and Raman interactions.

Goorjian, Peter M.↗

NASA NTRS1993

Summary of AXAF calibration requirements

The following summarizes requirements on the High Resolution Mirror Assembly (HRMA) and HRMA/SI calibration. The lists of calibration measurements assume that the HRMA meets the CTT requirement that the XRCF test environment shall have not more than a 10 percent effect on the encircled energy in 1 arcsec and that it can be calibrated to 1 percent. This implies that the offloading scheme has been implemented in the HRMA. It should be remembered that there are additional calibrations needed for: aspect system; tracking; gyros; other parts of the PCAS; spacecraft timing and checks on timing accuracy for all SIs after integration in the spacecraft; ficucial lights and periscope; alignments; optical metrology data: lengths and diameters with errors on mirror elements and optical interferometer data on surface figure; and throughput and imaging stability test: an end-to-end test that can be used after the XRCF to verify that the x-ray throughput and imaging quality have not been degraded. The tables presented give a summary of the integration times for HRMA and HRMA/SI calibration.

Kellogg, E.↗

DOE OSTI2020

BISON Robustness and Performance Improvements

BISON is a modern finite-element based nuclear fuel performance code that has been under development at the Idaho National Laboratory (USA) since 2009 [1]. The code is applicable to both steady and transient fuel behavior and can be used to analyze 1D (spherically symmetric), 2D (axisymmetric and generalized plane strain) or 3D geometries. BISON is the fuel performance code used within CASL for LWR fuel under both normal operating and accident conditions. BISON is built using the INL Multiphysics ObjectOriented Simulation Environment, or MOOSE [2, 3]. MOOSE is a massively parallel, finite element-based framework to solve systems of coupled non-linear partial differential equations using the Jacobian-Free Newton Krylov (JFNK) method [4]. This enables investigation of computationally large problems, for example a full stack of discrete pellets in a LWR fuel rod, or every rod in a full reactor core. MOOSE supports the use of complex two and three-dimensional meshes and uses implicit time integration, important for the widely varied time scale in nuclear fuel simulation. An object-oriented architecture is employed which greatly minimizes the programming effort required to add new material and behavioral models. The flexibility of the implicit and fully coupled multiphysics approach comes with a need for constructing suitable approximations for the Jacobian matrix of the coupled system used for either preconditioning a Krylov solve or in a direct Newton solve. Preconditioning options for Bison problems need to be revisited with new preconditioning methods becoming available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1984

Integrating amplifiers using cooled JFETs

It is shown how a simple integrating amplifier based on commercially available JFET and MOSFET switches can be used to measure photocurrents from detectors with noise levels as low as 1.6 x 10 to the -18th A/root Hz (10 electrons/sec). A figure shows the basic circuit, along with the waveform at the output. The readout is completely nondestructive; the reset noise does not contribute since sampling of the accumulated charge occurs between resets which are required only when the stored charge has reached a very high level. The storage capacity ranges from 10 to the 6th to 10 to the 9th electrons, depending on detector parameters and linearity requirements. Data taken with an Si:Sb detector operated at 24 microns are presented. The responsivity agrees well with the value obtained by Young et al. (1981) in the transimpedance amplifier circuit. The data are seen as indicating that extremely low values of NEP can be obtained for integration times of 1 sec and that longer integrations continue to improve the SNR at a rate faster than the square root of time when background noise is not present.

Low, F. J.↗

NASA NTRS1990

Temporal-adaptive Euler/Navier-Stokes algorithm for unsteady aerodynamic analysis of airfoils using unstructured dynamic meshes

A temporal adaptive algorithm for the time-integration of the two-dimensional Euler or Navier-Stokes equations is presented. The flow solver involves an upwind flux-split spatial discretization for the convective terms and central differencing for the shear-stress and heat flux terms on an unstructured mesh of triangles. The temporal adaptive algorithm is a time-accurate integration procedure which allows flows with high spatial and temporal gradients to be computed efficiently by advancing each grid cell near its maximum allowable time step. Results indicate that an appreciable computational savings can be achieved for both inviscid and viscous unsteady airfoil problems using unstructured meshes without degrading spatial or temporal accuracy.

Kleb, William L.↗

NASA NTRS1990

Temporal-adaptive Euler/Navier-Stokes algorithm for unsteady aerodynamic analysis of airfoils using unstructured dynamic meshes

Kleb, William L.↗

NASA NTRS1983

Optimal trajectories for aeroassisted orbital transfer

Consideration is given to classical and minimax problems involved in aeroassisted transfer from high earth orbit (HEO) to low earth orbit (LEO). The transfer is restricted to coplanar operation, with trajectory control effected by means of lift modulation. The performance of the maneuver is indexed to the energy expenditure or, alternatively, the time integral of the heating rate. Firist-order optimality conditions are defined for the classical approach, as are a sequential gradient-restoration algorithm and a combined gradient-restoration algorithm. Minimization techniques are presented for the aeroassisted transfer energy consumption and time-delay integral of the heating rate, as well as minimization of the pressure. It is shown that the eigenvalues of the Jacobian matrix of the differential system is both stiff and unstable, implying that the sequential gradient restoration algorithm in its present version is unsuitable. A new method, involving a multipoint approach to the two-poing boundary value problem, is recommended.

Miele, A.↗

DOE OSTI2025

Immersion Freezing in Particle-Based Aerosol-Cloud Microphysics: A Probabilistic Perspective on Singular and Time-Dependent Models

Cloud droplets containing immersed ice-nucleating particles (INPs) may freeze at temperatures above the homogeneous freezing threshold temperature in a process referred to as immersion freezing. In modeling studies, immersion freezing is often described using either so-called “singular” or “time-dependent” parameterizations. Here, we compare both approaches and discuss them in the context of probabilistic particle-based (super-droplet) cloud microphysics modeling. First, using a box model, we contrast how both parameterizations respond to idealized ambient cooling rate profiles and quantify the impact of the polydispersity of the immersed surface spectrum on the frozen fraction evolution. Presented simulations highlight that the singular approach, constituting a time-integrated form of a more general time-dependent approach, is only accurate under a limited range of ambient cooling rates. The time-dependent approach is free from this limitation. Second, using a prescribed-flow two-dimensional cloud model, we illustrate the macroscopic differences in the evolution in time of ice particle concentrations in simulations with flow regimes relevant to ambient cloud conditions. The flow-coupled aerosol-budget-resolving simulations highlight the benefits and challenges of modeling cloud condensation nuclei activation and immersion freezing on insoluble ice nuclei with super-particle methods. The challenges stem, on the one hand, from heterogeneous ice nucleation being contingent on the presence of relatively sparse immersed INPs, and on the other hand, from the need to represent a vast population of particles with relatively few so-called super particles (each representing a multiplicity of real particles). We discuss the critical role of the sampling strategy for particle attributes, including the INP size, the freezing temperature (for singular scheme) and the multiplicity.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1985

On the impact of initially stressed composite laminates

The impact response behavior of initially stressed composite laminates is investigated using the finite element method. An experimentally established contact law is incorporated into the finite element program. The Newmark time integration algorithm is used for solving the time dependent equations of the plate and the impactor. Numerical results, including the contact force history, deflection, and strain in the plate, are presented. Effects of impact velocity, initial stress, and the mass and size of the impactor are discussed.

Sun, C. T.↗

NASA NTRS1991

New approach to radar rainfall measurement

This paper integrates individual studies of the Area-Time Integrals method and climatic tuning methods. The latter is extended to a new approach which generalizes the technique so that it is no longer restricted to power law relations between effective reflectivity and rain rate.

Atlas, David↗