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 487 records · Page 27

All sources NASA NTRS DOE OSTI

DOE OSTI2021

An Orthogonal Recursive Bisection (ORB) Based Time Advancement Algorithm for CFD-DEM Solvers

The time integration of the granular phase in coupled computational fluid dynamics (CFD) – discrete element method (DEM) simulations presents a unique computational challenge brought about by the large variations in particle collisional time scales. Particles in the dilute regions of the computational domain can be advanced with large time steps while dense regions require much smaller time increments. However, the time step size in most solvers is globally set as the limit for accuracy and stability imposed by the collisions and is typically orders of magnitude less than that required away from collisions. This work addresses this precise issue and provides a strategy to avoid the use of a global conservative small time step size for the entire set of particles.A novel time stepping algorithm for CFD-DEM solvers using a partitioning approach using orthogonal recursive bisection (ORB) that allows for variable time steps among particles is described and its computational performance is compared against baseline explicit methods, typically used in several CFD-DEM solvers. ORB has advantages of being relatively quick and easy to update incrementally and has the required heuristic behavior (i.e., it will split the region in half with a cluster on each side) when groups of particles are well separated (clustered). The algorithm presented in this work uses a local time stepping approach to resolve collisional time scales for subsets of particles that are present at the leaves of the ORB, thereby resulting in substantial reduction of computational cost. The parallel implementation of this method where a ``knapsack” algorithm is used in tandem with ORB for effective load-balancing is also presented, where a best possible partitioning is obtained based on number of particles and local time-stepping costs. The algorithm is tested against benchmark problems with varying particle distributions that include fluidized bed and riser flow scenarios. Preliminary results indicate that the approach is 2-3X faster than traditional explicit methods for problems that involve both dense and dilute regions, while maintaining the same level of accuracy.

adaptive timestepping↗

NASA NTRS1977

Evaluation of subgrid-scale turbulence models using a fully simulated turbulent flow

An exact turbulent flow field was calculated on a three-dimensional grid with 64 points on a side. The flow simulates grid-generated turbulence from wind tunnel experiments. In this simulation, the grid spacing is small enough to include essentially all of the viscous energy dissipation, and the box is large enough to contain the largest eddy in the flow. The method is limited to low-turbulence Reynolds numbers, in our case R sub lambda = 36.6. To complete the calculation using a reasonable amount of computer time with reasonable accuracy, a third-order time-integration scheme was developed which runs at about the same speed as a simple first-order scheme. It obtains this accuracy by saving the velocity field and its first-time derivative at each time step. Fourth-order accurate space-differencing is used.

Clark, R. A.↗

NASA NTRS1988

SUMER: Solar Ultraviolet Measurements of Emitted Radiation

The SUMER (solar ultraviolet measurements of emitted radiation) experiment is described. It will study flows, turbulent motions, waves, temperatures and densities of the plasma in the upper atmosphere of the Sun. Structures and events associated with solar magnetic activity will be observed on various spatial and temporal scales. This will contribute to the understanding of coronal heating processes and the solar wind expansion. The instrument will take images of the Sun in EUV (extreme ultra violet) light with high resolution in space, wavelength and time. The spatial resolution and spectral resolving power of the instrument are described. Spectral shifts can be determined with subpixel accuracy. The wavelength range extends from 500 to 1600 angstroms. The integration time can be as short as one second. Line profiles, shifts and broadenings are studied. Ratios of temperature and density sensitive EUV emission lines are established.

Wilhelm, K.↗

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.↗

DOE OSTI2023

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗

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.↗

NASA NTRS1993

Climate forcings and feedbacks

Global temperature has increased significantly during the past century. Understanding the causes of observed global temperature change is impossible in the absence of adequate monitoring of changes in global climate forcings and radiative feedbacks. Climate forcings are changes imposed on the planet's energy balance, such as change of incoming sunlight or a human-induced change of surface properties due to deforestation. Radiative feedbacks are radiative changes induced by climate change, such as alteration of cloud properties or the extent of sea ice. Monitoring of global climate forcings and feedbacks, if sufficiently precise and long-term, can provide a very strong constraint on interpretation of observed temperature change. Such monitoring is essential to eliminate uncertainties about the relative importance of various climate change mechanisms including tropospheric sulfate aerosols from burning of coal and oil smoke from slash and burn agriculture, changes of solar irradiance changes of several greenhouse gases, and many other mechanisms. The considerable variability of observed temperature, together with evidence that a substantial portion of this variability is unforced indicates that observations of climate forcings and feedbacks must be continued for decades. Since the climate system responds to the time integral of the forcing, a further requirement is that the observations be carried out continuously. However, precise observations of forcings and feedbacks will also be able to provide valuable conclusions on shorter time scales. For example, knowledge of the climate forcing by increasing CFC's relative to the forcing by changing ozone is important to policymakers, as is information on the forcing by CO2 relative to the forcing by sulfate aerosols. It will also be possible to obtain valuable tests of climate models on short time scales, if there is precise monitoring of all forcings and feedbacks during and after events such as a large volcanic eruption or an El Nino.

Hansen, James↗

NASA NTRS2008

High Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

Taylor series integration is implemented in a spacecraft trajectory analysis code-the Spacecraft N-body Analysis Program (SNAP) - and compared with the code s existing eighth-order Runge-Kutta Fehlberg time integration scheme. Nine trajectory problems, including near Earth, lunar, Mars and Europa missions, are analyzed. Head-to-head comparison at five different error tolerances shows that, on average, Taylor series is faster than Runge-Kutta Fehlberg by a factor of 15.8. Results further show that Taylor series has superior convergence properties. Taylor series integration proves that it can provide rapid, highly accurate solutions to spacecraft trajectory problems.

Scott, James R.↗

DOE OSTI2022

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS1971

Improved reversible coulometer cell

Cell operates either as timer or current-time integrating device in any physical orientation with better than 2 percent accuracy over temperature range of 283 K to 398 K, with input current range from few microamperes to approximately 1000 microamperes over time period of 1 sec to several hours.

Boettcher, G. E.↗

NASA NTRS1984

Inelastic and Dynamic Fracture and Stress Analyses

Large deformation inelastic stress analysis and inelastic and dynamic crack propagation research work is summarized. The salient topics of interest in engine structure analysis that are discussed herein include: (1) a path-independent integral (T) in inelastic fracture mechanics, (2) analysis of dynamic crack propagation, (3) generalization of constitutive relations of inelasticity for finite deformations , (4) complementary energy approaches in inelastic analyses, and (5) objectivity of time integration schemes in inelastic stress analysis.

Atluri, S. N.↗

NASA NTRS1992

A transient FETI methodology for large-scale parallel implicit computations in structural mechanics

Explicit codes are often used to simulate the nonlinear dynamics of large-scale structural systems, even for low frequency response, because the storage and CPU requirements entailed by the repeated factorizations traditionally found in implicit codes rapidly overwhelm the available computing resources. With the advent of parallel processing, this trend is accelerating because explicit schemes are also easier to parallelize than implicit ones. However, the time step restriction imposed by the Courant stability condition on all explicit schemes cannot yet -- and perhaps will never -- be offset by the speed of parallel hardware. Therefore, it is essential to develop efficient and robust alternatives to direct methods that are also amenable to massively parallel processing because implicit codes using unconditionally stable time-integration algorithms are computationally more efficient when simulating low-frequency dynamics. Here we present a domain decomposition method for implicit schemes that requires significantly less storage than factorization algorithms, that is several times faster than other popular direct and iterative methods, that can be easily implemented on both shared and local memory parallel processors, and that is both computationally and communication-wise efficient. The proposed transient domain decomposition method is an extension of the method of Finite Element Tearing and Interconnecting (FETI) developed by Farhat and Roux for the solution of static problems. Serial and parallel performance results on the CRAY Y-MP/8 and the iPSC-860/128 systems are reported and analyzed for realistic structural dynamics problems. These results establish the superiority of the FETI method over both the serial/parallel conjugate gradient algorithm with diagonal scaling and the serial/parallel direct method, and contrast the computational power of the iPSC-860/128 parallel processor with that of the CRAY Y-MP/8 system.

Farhat, Charbel↗

NASA NTRS2008

Operational Characteristics of a Low-Energy FARAD Thruster

Data from a 100 J per pulse electrodeless accelerator employing pulsed RF-preionization are presented to gain insight into the accelerator's operating characteristics. The data suggest that the propellant distribution is highly unoptimized, with most of the gas inaccessible to the discharge and the remainder mostly concentrated at the inner radius of the coil. The pulsed RF-preionization discharge produces a visible plasma, but like the gas distribution it mostly appears concentrated at the inner radius of the thruster. Magnetic field probes in the discharge point to a current sheet that is not magnetically impermeable. These data also exhibit signs of nonrepeatability, and time-integrated discharge photography shows signs of spatial nonuniformity in both the radial and azimuthal directions. Terminal voltage measurements on the two capacitor banks of the thruster do not exhibit the asymmetric nature (in time) typically associated with an efficient pulsed plasma accelerator. Based on the experimental evidence, the poor performance of the thruster is thought to be due to insufficient preionization, which at these low discharge energy levels severely limits the ability of the main current pulse to couple with and effectively accelerate the propellant.

Polzin, Kurt A.↗

NASA NTRS1985

Numerical solution of the Euler equations for high-speed, blended wing-body configurations

Simulation of high-speed three-dimensional flow about blended wing-body combinations is investigated. A finite-volume explicit scheme with Runge-Kutta time integration is used to solve the compressible Euler equations in order to simulate the flow. The method, augmented by carefully chosen dissipative terms and convergence accelerators such as enthalpy damping and maximum local time-stepping, has been found to be very efficient in solving high-speed flows involving strong shocks. An analytic method for generating the body geometry and an algebraic method for quasi-three-dimensional grid generation are described. Results are presented for various blended wing-body configurations at different Mach numbers and angles of attack.

Moitra, A.↗

NASA NTRS1985

Electronically Steered Thinned Array Radiometer (ESTAR) system design, calibration, and sensitivity

The ESTAR system, an imaging L-band radiometer designed to provide measurements of soil moisture from space, is described. The large aperture size dictated by the resolution requirement (10 km or less) suggests the use of aperture synthesis, and plans to have ESTAR operate in low earth orbit require that the synthesis be performed in the short data acquisition time permitted by satellite motion. Two forms of ESTAR are described. ESTAR-0 is to be an airborne prototype designed for proof of concept; it will image a line of eight adjacent pixels with 8 degree spatial resolution of 0.1 degree Kelvin temperature resolution, given 1 sec integration time. ESTAR-1 is the satellite version to operate from space. A tentative sensor design will image a line of 128 0.5 degree pixels with 2.0 degree K resolution assuming 0.7 sec integration at 700 km altitude.

Ruf, C.↗

DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗

NASA NTRS1975

An infrared heterodyne radiometer for high-resolution measurements of solar radiation and atmospheric transmission

A tunable, Dicke-switched, infrared heterodyne radiometer (IHR) has been designed, fabricated, tested, and used to observe solar radiation and determine the atmospheric transmissivity in the 9-11-micron spectral band. The IHR provides a spectral resolution of 0.0067 reciprocal cm, a minimum detectable power level of 2.9 times ten to the minus twenty-third power W/Hz, and a temperature resolution of less than 1 K for a source temperature of 1000 K, an IF predetection bandwidth of 100 MHz, and an integration time of 30 s. Detailed design equations and measured IHR performance are presented. The IHR was used to make solar and atmospheric transmission measurements, and a vertical-path atmospheric attenuation of 2.3 dB has been established under favorable weather conditions. The attenuation of solar radiation due to cloud cover and haze has also been investigated.

Peyton, B. J.↗

NASA NTRS1978

Passive imaging through the turbulent atmosphere - Fundamental limits on the spatial frequency resolution of a rotational shearing interferometer

The signal-to-noise (S/N) ratio to be expected when a 180 deg rotationally shearing interferometer is used for image recovery at the diffraction limit of a large telescope is computed. The variance and covariance of the irradiance fluctuations at the detector array are shown to yield measures of the high-frequency spatial spectrum of the source. Four fundamental sources of noise are considered: temporal fluctuations of the source, space-time fluctuations of the atmosphere, shot noise in the detected photocurrents, and the effects of finite sampling. S/N is found to be directly proportional to the angular resolution of the telescope, the single-frame integration time, the square root of the number of frames, the cube of the operating wavelength, the quantum efficiency of the detector, and the average spectral irradiance from the source on the pupil. It is inversely proportional to the cube of the field angle subtended by the source (or part thereof) under study.

Burke, J. J.↗