Engineering Papers⌕ Search

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Stability of mixed time integration schemes for transient thermal analysis

A current research topic in coupled-field problems is the development of effective transient algorithms that permit different time integration methods with different time steps to be used simultaneously in various regions of the problems. The implicit-explicit approach seems to be very successful in structural, fluid, and fluid-structure problems. This paper summarizes this research direction. A family of mixed time integration schemes, with the capabilities mentioned above, is also introduced for transient thermal analysis. A stability analysis and the computer implementation of this technique are also presented. In particular, it is shown that the mixed time implicit-explicit methods provide a natural framework for the further development of efficient, clean, modularized computer codes.

Liu, W. K.↗
DOE OSTI2022

Parallel-in-Time Integration for Nonlinear Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to diffusion-dominated PDEs, with some recent success in scalar linear hyperbolic PDEs. Given the ubiquity of numerical simulations of hyperbolic PDEs throughout the sciences, in particular, nonlinear hyperbolic systems, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems beyond simple scalar and linear cases, which is the main topic of this subcontract. The main focus of the work was to further develop and perfect coarse-grid operators for the Multigrid Reduction-inTime (MGRIT) method applied to hyperbolic PDEs that were recently proposed in PhD thesis, based on a modified semi-Lagrangian approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Parallel-In-Time Integration for Hyperbolic Problems (Final Report)

The work for the subcontract was situated in the area of parallel-in-time integration for advection-dominated and hyperbolic partial differential equations (PDEs). Parallel-in-time integration is an active area of research due to its ability to enable faster numerical simulations for applications throughout many areas of science. Over the past two decades, much progress has been made in this area; however, this progress has largely been limited to non-hyperbolic, diffusion-dominated PDEs. Given the ubiquity of numerical simulations of advection-dominated hyperbolic PDEs throughout the sciences, there is a strong need to develop efficient parallel-in-time techniques for hyperbolic problems, which was the main topic of this subcontract.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A conservative discontinuous-Galerkin-in-time (DGiT) multirate time integration framework for interface-coupled problems with applications to solid–solid interaction and air–sea models

In this paper we extend the DGiT multirate framework, developed in Connors and Sockwell (2022) for scalar transmission problems, to a solid–solid interaction (SSI) problem involving two coupled elastic solids and a coupled air–sea model with the rotating, thermal shallow water equations. In so doing we aim to demonstrate the broad applicability of the mathematical theory and governing principles established in Connors and Sockwell (2022) to coupled problems characterized by subproblems evolving at different temporal scales. Further, multirate time integration algorithms employing different time steps, optimized for the dynamics of each subproblem, can significantly improve simulation efficiency for such coupled problems. However, development of multirate algorithms is a highly non-trivial task due to the coupling, which can impact accuracy, stability or other desired properties such as preservation of system invariants. DGiT provides a general template for multirate time integration that can achieve these properties. To elucidate the manner in which DGiT accomplishes this task, we fully detail each step in the application of the framework to the SSI and air–sea coupled problems. Numerical examples illustrate key properties of the resulting multirate schemes for both problems.

42 ENGINEERING↗
NASA NTRS1992

Direct time integration of Maxwell's equations in nonlinear dispersive media for propagation and scattering of femtosecond electromagnetic solitons

The initial results for femtosecond electromagnetic soliton propagation and collision obtained from first principles, i.e., by a direct time integration of Maxwell's equations are reported. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit the modeling of 2D and 3D optical soliton propagation, scattering, and switching from the full-vector Maxwell's equations.

Goorjian, Peter M.↗
DOE OSTI2023

Efficient Integration of Algebraic Constraint for Exponential Time Integration

This work is a continuation of a previous project where the high efficiency of exponential integration methods for DRMHD systems was demonstrated. Often algebraic constraints must also be enforced for DRMHD systems of interest. Straightforward application of exponential methods to DRMHD equations with constraints leads to prohibitively computationally expensive methods. In this work, we propose new exponential schemes that allow the constraints to be removed from evaluation of exponential matrix functions which drastically reduces computational complexity by eliminating the need to perform computations enforcing constraints in exponential calculations while still preserving a high order of accuracy and allowing for a large time step even when the problem is stiff. This idea is similar to the W-methods and is achieved by carefully designing a method with the desired order of accuracy, even with an incomplete Jacobian matrix used as an argument of exponential-like functions. The constraints are accounted for by including them in the evaluation of the right-hand-side forcing function of the spatially discretized system. We study performance of the new methods on test problems and outline future research directions that this work opens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An assessment of the utility of multirate time integration for the modeling of laser powder bed fusion

Finite element simulation of the laser powder bed fusion process is made challenging by the inherently multiscale nature of the process. When using the typical global time stepping techniques, slowly-evolving regions of the domain receive the same numerical treatment as the regions with the highest temperature rates. The current work details the implementation and evaluation of an implicit multirate method which is able to advance different regions of the domain with distinct time steps, depending on their current solution rate. Previous work indicates that this representation of the temporal scales of the problem can lead to significant reductions in wall clock run time, and it is shown herein that multirate time integration, when used with a uniform process-scale mesh, can result in speedups of approximately 19, 42, and 87 times, for domains with edge lengths of approximately 1 mm, 2 mm, and 5 mm, respectively. When used in conjunction with $\ h$-refinement (limited to two levels of refinement), the resulting speedups (taken relative to the uniform mesh) are around 28, 85, and 665 times. The method is demonstrated to converge as indicated by the literature, and its use with an AM-Bench domain is demonstrated.

42 ENGINEERING↗
NASA NTRS1994

The Relationship Between Area-Time Integrals Determined from Satellite Infrared Data by Means of a Fixed-Threshold Approach and Convective Rainfall Volumes

The relationship of the rainfall from convective clouds to area-time integrals determined from satellite infrared data using a fixed infrared-temperature threshold is investigated. Concurrent radar and rapid-scan satellite data obtained during field projects in the northern High Plains and the southeastern United States were used in this study. The fixed IR threshold appropriate for each region was determined by an optimization procedure that identified the brightness threshold that yields the strongest relationship between estimated rainfall from a cloud cluster and its satellite Area-Time Integral (ATI) for each dataset. For the North Dakota-Montana area the optimization procedure indicated that the area enclosed by the -22.5 C isotherm provides satellite ATI values most closely related to the estimated rainfalls. For the southeastern United States project, the optimized tem- perature threshold was 8.5'C. The difference between the thresholds determined for the two geographic areas suggests that a different 'calibration' for each distinct area may be needed to make use of this relationship. Slopes of the two log-log rainfall-ATI regressions are less than unity, indicating that a relative horizontal expansion and/or increase in persistence of a cloud cluster exceeds the associated increase in precipitation. Implications for the Geostationary Operational Environmental Satellite precipitation index are discussed. New results concerning the rain volume-radar ATI relationship for the southeastern United States are also appended to the paper.

Johnson, L. Ronald↗
NASA NTRS1994

Multiparameter radar study of rainfall: Potential application to area-time integral studies

Multiparameter radars measure one or more additional parameters in addition to the coventional reflectivity factor. The combination of radar observations from a multiparameter radar is used to study the time evolution of rainstorms. A technique is presented to self-consistently compare the area-time integral (ATI) and rainfall volume estimates from convective storms, using two different measurements from a multiparameter radar. Rainfall volumes for the lifetime of individual storms are computed using the reflectivity at S band (10-cm wavelength) as well as one-way specific attenuation at X band (3-cm wavelength). Area-time integrals are computed by summing all areas in each radar snapshot having reflectivities (S band) in excess of a preselected threshold. The multiparameter radar data used in this study were acquired by the National Center for Atmospheric Research (NCAR) CP-2 radar during the Cooperative Huntsville Meteorological Experiment (COHMEX) and the Convection and Precipitation/Electrification Experiment (CaPE), respectively. ATI studies were accomplished in this work using multiparameter radar data acquired during the lifetime of six convective events that occurred in the COHMEX radar coverage area. A case study from the COMHEX field campaign (20 July 1986) was selected to depict the various stages in the evolution of a storm over which the ATI and rainfall volume computations were performed using multiparameter radar data. Another case study from the CaPE field campaign (12 August 1991) was used to demonstrate the evolution of a convective cell based on differential reflectivity observations.

Raghavan, R.↗
NASA NTRS1973

A portable millisecond-integration-time photoelectric photometer.

Portable equipment for recording millisecond-integration-time photoelectric photometric data is described. Digital data are reliably recorded on standard 6.35 mm audio grade magnetic tape via a quadradial audio grade tape deck. The system is designed specifically for recording lunar occulations of stars, but the data recording technique is independent of the data source. Recovery of the data is made via minicomputer.

Mcgraw, J. T.↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a different yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1983

Mixed time integration methods for transient thermal analysis of structures

The computational methods used to predict and optimize the thermal-structural behavior of aerospace vehicle structures are reviewed. In general, two classes of algorithms, implicit and explicit, are used in transient thermal analysis of structures. Each of these two methods has its own merits. Due to the different time scales of the mechanical and thermal responses, the selection of a time integration method can be a difficult yet critical factor in the efficient solution of such problems. Therefore mixed time integration methods for transient thermal analysis of structures are being developed. The computer implementation aspects and numerical evaluation of these mixed time implicit-explicit algorithms in thermal analysis of structures are presented. A computationally-useful method of estimating the critical time step for linear quadrilateral element is also given. Numerical tests confirm the stability criterion and accuracy characteristics of the methods. The superiority of these mixed time methods to the fully implicit method or the fully explicit method is also demonstrated.

Liu, W. K.↗
NASA NTRS1993

On the spectral stability of time integration algorithms for a class of constrained dynamics problems

Incomplete field formulations have recently been the subject of intense research because of their potential in coupled analysis of independently modeled substructures, adaptive refinement, domain decomposition, and parallel processing. This paper discusses the design and analysis of time-integration algorithms for these formulations and emphasizes the treatment of their inter-subdomain constraint equations. These constraints are shown to introduce a destabilizing effect in the dynamic system that can be analyzed by investigating the behavior of the time-integration algorithm at infinite and zero frequencies. Three different approaches for constructing penalty-free unconditionally stable second-order accurate solution procedures for this class of hybrid formulations are presented, discussed and illustrated with numerical examples. The theoretical results presented in this paper also apply to a large family of nonlinear multibody dynamics formulations. Some of the algorithms outlined herein are important alternatives to the popular technique consisting of transforming differential/algebraic equations into ordinary differential equations via the introduction of a stabilization term that depends on arbitrary constants and that influences the computed so1ution.

Farhat, Charbel↗
NASA NTRS2002

Higher Order Time Integration Schemes for the Unsteady Navier-Stokes Equations on Unstructured Meshes

The rapid increase in available computational power over the last decade has enabled higher resolution flow simulations and more widespread use of unstructured grid methods for complex geometries. While much of this effort has been focused on steady-state calculations in the aerodynamics community, the need to accurately predict off-design conditions, which may involve substantial amounts of flow separation, points to the need to efficiently simulate unsteady flow fields. Accurate unsteady flow simulations can easily require several orders of magnitude more computational effort than a corresponding steady-state simulation. For this reason, techniques for improving the efficiency of unsteady flow simulations are required in order to make such calculations feasible in the foreseeable future. The purpose of this work is to investigate possible reductions in computer time due to the choice of an efficient time-integration scheme from a series of schemes differing in the order of time-accuracy, and by the use of more efficient techniques to solve the nonlinear equations which arise while using implicit time-integration schemes. This investigation is carried out in the context of a two-dimensional unstructured mesh laminar Navier-Stokes solver.

Jothiprasad, Giridhar↗
NASA NTRS2003

Effect of Variations in IRU Integration Time Interval On Accuracy of Aqua Attitude Estimation

During Aqua launch support, attitude analysts noticed several anomalies in Onboard Computer (OBC) rates and in rates computed by the ground Attitude Determination System (ADS). These included: 1) periodic jumps in the OBC pitch rate every 2 minutes; 2) spikes in ADS pitch rate every 4 minutes; 3) close agreement between pitch rates computed by ADS and those derived from telemetered OBC quaternions (in contrast to the step-wise pattern observed for telemetered OBC rates); 4) spikes of +/- 10 milliseconds in telemetered IRU integration time every 4 minutes (despite the fact that telemetered time tags of any two sequential IRU measurements were always 1 second apart from each other). An analysis presented in the paper explains this anomalous behavior by a small average offset of about 0.5 +/- 0.05 microsec in the time interval between two sequential accumulated angle measurements. It is shown that errors in the estimated pitch angle due to neglecting the aforementioned variations in the integration time interval by the OBC is within +/- 2 arcseconds. Ground attitude solutions are found to be accurate enough to see the effect of the variations on the accuracy of the estimated pitch angle.

Natanson, G. A.↗
NASA NTRS1982

Mixed time integration methods for transient thermal analysis of structures, appendix 5

Mixed time integration methods for transient thermal analysis of structures are studied. An efficient solution procedure for predicting the thermal behavior of aerospace vehicle structures was developed. A 2D finite element computer program incorporating these methodologies is being implemented. The performance of these mixed time finite element algorithms can then be evaluated employing the proposed example problem.

Liu, W. K.↗
DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2012

Neutrino Analysis of the September 2010 Crab Nebula Flare and Time-Integrated Constraints on Neutrino Emission from the Crab Using IceCube

We present the results for a search of high-energy muon neutrinos with the IceCube detector in coincidence with the Crab nebula flare reported on September 2010 by various experiments. Due to the unusual flaring state of the otherwise steady source we performed a prompt analysis of the 79-string configuration data to search for neutrinos that might be emitted along with the observed gamma-rays. We performed two different and complementary data selections of neutrino events in the time window of 10 days around the flare. One event selection is optimized for discovery of E(sub v)(sup -2) neutrino spectrum typical of 1st order Fermi acceleration. A similar event selection has also been applied to the 40-string data to derive the time-integrated limits to the neutrino emission from the Crab [35]. The other event selection was optimized for discovery of neutrino spectra with softer spectral index and TeV energy cut-offs as observed for various galactic sources in gamma-rays. The 90% CL best upper limits on the Crab flux during the 10 day flare are 4.73 x 10(exp -11) per square centimeter per second TeV(sup -1) for an E(sub v)(sup -2) neutrino spectrum and 2.50 x 10(exp -10) per square centimeter per second TeV(sup -1) for a softer neutrino spectra of E(sub v)(sup -2.7), as indicated by Fermi measurements during the flare. IceCube has also set a time-integrated limit on the neutrino emission of the Crab using 375.5 days of livetime of the 40-string configuration data. This limit is compared to existing models of neutrino production from the Crab and its impact on astrophysical parameters is discussed. The most optimistic predictions of some models are already rejected by the IceCube neutrino telescope with more than 90% CL.

Abbasi, R.↗