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 19 records

Comparison of exponential integrators and traditional time integration schemes for the shallow water equations

We report the time integration scheme is probably one of the most fundamental choices in the development of an ocean model. In this paper, we investigate several time integration schemes when applied to the shallow water equations. This set of equations is accurate enough for the modeling of a shallow ocean and is also relevant to study as it is the one solved for the barotropic (i.e. vertically averaged) component of a three dimensional ocean model. We analyze different time stepping algorithms for the linearized shallow water equations. High order explicit schemes are accurate but the time step is constrained by the Courant-Friedrichs-Lewy stability condition. Implicit schemes can be unconditionally stable but, in practice lack accuracy when used with large time steps. In this paper we propose a detailed comparison of such classical schemes with exponential integrators. The accuracy and the computational costs are analyzed in different configurations.

97 MATHEMATICS AND COMPUTING↗

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

Estimation of areal rainfall using the radar echo area time integral

The Area Time Integral (ATI) method of Doneaud et al. (1984) is extended to the measurement of cumulative areawide rainfall for periods up to 12 h. The extended ATI method is used to analyze data from the Florida Area Cumulus Experiment II. The relationship between radar estimated rain volume and radar-measured area covered with echoes is examined to test the possibility of obtaining values similar to conventional reflectivity-rainfall estimates of rainfall using only area measurements. The correlation between gage-estimated rain volume and radar estimated area covered with showers is also studied, focusing on the possible estimation of rain volume values using a small number of echo area observations.

Lopez, Raul E.↗

Simulating Nonlinear Deformations of Solar Sail Membranes Using Explicit Time Integration

In this study, the explicit time integration method is employed to predict deformation of highly flexible solar sail structural components. The nonlinear static analysis of a highly flexible ribbon structure is presented to demonstrate the need for having the explicit time integration method in the analysis toolbox for solar sail. Static analyses of the ribbon structure produce ambiquous results whereas the explicit time integration method determines the correct results. Extensive benchmarking examples are also presented to build confidence in the use of the explicit method. Previously determined nonlinear wrinkling deformations of solar sail membranes are found by the explicit method. As the explicit method is known to often require more computational time than nonlinear static methods, a study on mass scaling was also conducted. The computational times are reported for the nonlinear static and explicit time integration solutions to calibrate the advantage of using mass scaling for these problems.

Wang, John T.↗

A unified formulation of splitting-based implicit time integration schemes

Splitting-based time integration approaches such as fractional step, alternating direction implicit, operator splitting, and locally one dimensional methods partition the system of interest into components, and solve individual components implicitly in a cost-effective way. Here this work proposes a unified formulation of splitting time integration schemes in the framework of general-structure additive Runge–Kutta (GARK) methods. Specifically, we develop implicit-implicit (IMIM) GARK schemes, provide the order conditions for this class, and explain their application to partitioned systems of ordinary differential equations. We show that classical splitting methods belong to the IMIM GARK family, and therefore can be studied in this unified framework. New IMIM-GARK splitting methods are developed and tested using parabolic systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods

A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.

97 MATHEMATICS AND COMPUTING↗

A new explicit variable time-integration self-starting methodology for computational structural dynamics

A new explicit variable time-integration methodology and architecture which possesses self-starting attributes, eliminates the need to involve acceleration computations, and which has improved accuracy characteristics in comparison to the traditional central-difference-type formulations customarily advocated is described for applicability to computational structural dynamics. To sharpen the focus of the present study, an explicit variable time-integration architecture which is relatively simple, yet effective, is described. Unlike variable explicit time-integration formulations adopted in the past, the present self-starting variable time-integration architecture and implementation aspects facilitate a simplified representation and a straightforward and effective approach for combining finite element meshes requiring different time steps in a single analysis. Numerical test cases are provided which demonstrate the applicability of the proposed formulations.

Tamma, Kumar K.↗

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

Finite difference schemes for long-time integration

Finite difference schemes for the evaluation of first and second derivatives are presented. These second order compact schemes were designed for long-time integration of evolution equations by solving a quadratic constrained minimization problem. The quadratic cost function measures the global truncation error while taking into account the initial data. The resulting schemes are applicable for integration times fourfold, or more, longer than similar previously studied schemes. A similar approach was used to obtain improved integration schemes.

Haras, Zigo↗

GaAs-based photorefractive time-integrating correlator

A potential application of the photorefractive time-integrating correlator is the real-time radar jamming interference rejection system, using the adaptive filter method; a fast photorefractive crystal is needed for adapting a rapidly changing jamming signal. An effort is presently made to demonstrate and characterize a GaAs-based photorefractive time-integrating correlator, since GaAs crystals are 2-3 orders of magnitude faster than most other alternatives.

Liu, Duncan T. H.↗

Conservative High-Order Time Integration for Lagrangian Hydrodynamics

In tis study, we develop novel time integration methods for the compressible Euler equations in the Lagrangian frame that are of arbitrary high order and exactly preserve the mass, momentum, and total energy of the system. The equations are considered in nonconservative form, that is, common for staggered grid hydrodynamics (SGH) methods; namely, the evolved quantities are mass, momentum, and internal energy. A general family of time integration schemes is formulated, and practical pairs for orders three and four are derived. Numerical results on standard hydrodynamics benchmarks confirm the high-order convergence on smooth problems and the exact numerical preservation of all physically conserved quantities.

97 MATHEMATICS AND COMPUTING↗

Interferometric surface-wave acousto-optic time-integrating correlators

A structure for a coherent-interferometric acousto-optic (AO) time-integrating correlator was implemented by using a single surface acoustic wave (SAW) device with tilted transducers to reduce intermodulation terms. The SAW device was fabricated on Y-Z LiNbO3 with a center frequency of 175 MHz, a bandwidth of 60 MHz, and a time aperture of about 10 micros. The density of the photodetector array, with a potential of 120 MHz. Typical integration times are 30 to 40 ms, providing processing gains in excess of 10 to the 6th power. Such a device is useful in providing fast synchronization of communication links and in demodulating to base band and simultaneously acting as a synchronization lock monitor for moderate data rates. Where processing may be limited by Doppler shifts, a two dimensional architecture was implemented to allow full processing gain. Two one-dimensional, SAW AO time-integrating correlators and a two dimensional correlator are evaluated.

Berg, N. J.↗

Integration Time Required to Extract Accurate Data from Transonic Wind-Tunnel Tests

Because the forces and pressures on wind-tunnel models tested at transonic speeds are not steady, even for static aerodynamic tests, integration time is required to obtain data of acceptable accuracy. The integration time required for both static and dynamic tests is evaluated analytically and confirmed by experimental measurements. It is shown that, for static and dynamic tests, the accuracy obtained is a function of integration time, frequency of the signal, and the ratio of the dynamic amplitude to the full signal of interest. In addition, for the dynamic case, the frequency bandwidth used in analysis is important. Results of this study indicate that, for typical data accuracy desired from models in a large transonic wind tunnel (11- by 11-ft), up to the following integration times are required: static force and moment tests, 0.5 s; static pressure tests, 1 s; flutter tests, 30 to 60 s; and random-dynamic tests, 10 s.

Muhlstein, Lado, Jr.↗

Integration time required to extract accurate static and dynamic data from transonic wind-tunnel tests

Because the forces and pressures on wind-tunnel models tested at transonic speeds are not steady, even for 'static' aerodynamic tests, integration time is required to obtain data of acceptable accuracy. The integration time required for both static and dynamic tests is evaluated analytically and confirmed by experimental measurements. It is shown that, for static and dynamic tests, the accuracy obtained is a function of integration time, frequency content of the signal, and the ratio of the dynamic amplitude to the full signal of interest. In addition, for the dynamic case, the frequency band width used in analysis is important. Results of this study indicate that, for typical data accuracy desired from models in a large transonic wind tunnel (11- by 11-ft), up to the following integration times are required: static force and moment tests, 0.5 sec; static pressure tests, 0.1 sec; flutter tests, 30 to 60 sec; and random-dynamic tests, 10 sec.

Muhlstein, L., Jr.↗

The relation of radar to cloud area-time integrals and implications for rain measurements from space

The relationships between satellite-based and radar-measured area-time integrals (ATI) for convective storms are determined, and both are shown to depend on the climatological conditional mean rain rate and the ratio of the measured cloud area to the actual rain area of the storms. The GOES precipitation index of Arkin (1986) for convective storms, an area-time integral for satellite cloud areas, is shown to be related to the ATI for radar-observed rain areas. The quality of GPI-based rainfall estimates depends on how well the cloud area is related to the rain area and the size of the sampling domain. It is also noted that the use of a GOES cloud ATI in conjunction with the radar area-time integral will improve the accuracy of rainfall estimates and allow such estimates to be made in much smaller space-time domains than the 1-month and 5-deg boxes anticipated for the Tropical Rainfall Measuring Mission.

Atlas, David↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Virtual-pulse time integral methodology: A new explicit approach for computational dynamics - Theoretical developments for general nonlinear structural dynamics

The present paper describes a new explicit virtual-pulse time integral methodology for nonlinear structural dynamics problems. The purpose of the paper is to provide the theoretical basis of the methodology and to demonstrate applicability of the proposed formulations to nonlinear dynamic structures. Different from the existing numerical methods such as direct time integrations or mode superposition techniques, the proposed methodology offers new perspectives and methodology of development, and possesses several unique and attractive computational characteristics. The methodology is tested and compared with the implicit Newmark method (trapezoidal rule) through a nonlinear softening and hardening spring dynamic models. The numerical results indicate that the proposed explicit virtual-pulse time integral methodology is an excellent alternative for solving general nonlinear dynamic problems.

Chen, Xiaoqin↗