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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Exoplanet Transit Spectroscopy Using WFC3: WASP-12b, WASP-17b, and WASP-19b

We report an analysis of transit spectroscopy of the extrasolar planets WASP-12 b, WASP-17 b, and WASP-19 b using the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST). We analyze the data for a single transit for each planet using a strategy similar, in certain aspects, to the techniques used by Berta et al., but we extend their methodology to allow us to correct for channel- or wavelength-dependent instrumental effects by utilizing the band-integrated time series and measurements of the drift of the spectrum on the detector over time. We achieve almost photon-limited results for individual spectral bins, but the uncertainties in the transit depth for the band-integrated data are exacerbated by the uneven sampling of the light curve imposed by the orbital phasing of HST's observations. Our final transit spectra for all three objects are consistent with the presence of a broad absorption feature at 1.4 microns most likely due to water. However, the amplitude of the absorption is less than that expected based on previous observations with Spitzer, possibly due to hazes absorbing in the NIR or non-solar compositions. The degeneracy of models with different compositions and temperature structures combined with the low amplitude of any features in the data preclude our ability to place unambiguous constraints on the atmospheric composition without additional observations with WFC3 to improve the signal-to-noise ratio and/or a comprehensive multi-wavelength analysis. Key words: planetary systems - techniques: photometric - techniques: spectroscopic

HST↗
DOE OSTI2023

Exploiting Kronecker structure in exponential integrators: Fast approximation of the action of φ $-$functions of matrices via quadrature

Here, in this article, we propose an algorithm for approximating the action of φ $-$ functions of matrices against vectors, which is a key operation in exponential time integrators. In particular, we consider matrices with Kronecker sum structure, which arise from problems admitting a tensor product representation. The method is based on quadrature approximations of the integral form of the φ $-$ functions combined with a scaling and modified squaring method. Owing to the Kronecker sum representation, only actions of 1D matrix exponentials are needed at each quadrature node and assembly of the full matrix can be avoided. Additionally, we derive a priori bounds for the quadrature error, which show that, as expected by classical theory, the rate of convergence of our method is supergeometric. Guided by our analysis, we construct a fast and robust method for estimating the optimal scaling factor and number of quadrature nodes that minimizes the total cost for a prescribed error tolerance. We investigate the performance of our algorithm by solving several linear and semilinear time-dependent problems in 2D and 3D. The results show that our method is accurate and orders of magnitude faster than the current state-of-the-art.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Late Paleozoic Evolution of the Anadarko Basin: Implications for Laurentian Tectonics and the Assembly of Pangea

Abstract Interplay between marginal and intraplate late Paleozoic tectonics in southeastern Laurentia is a subject of debate. In southern Laurentia, intraplate deformation has been attributed to combinations of compressive forces from orogenic belts in the west, southwest, and southeast. While left‐lateral transpression is generally an accepted mechanism for formation of uplifts in the Southern Oklahoma Aulacogen (SOA), to date no studies have explained the driving mechanisms of these kinematics. Here, we analyze the Southern Oklahoma Transpressional System (SOTS), a ∼500 km long ∼50 km wide fault zone that extends from the southeastern Laurentian margin into the plate interior along Cambrian lithospheric weaknesses inherited from the SOA. We use a compilation of published studies to present detailed, time‐integrated sedimentary thickness maps of SOTS‐related basins. We combine these maps with a new interpretation of three 2D seismic reflection lines to constrain timing and kinematics of all major late Paleozoic SOTS fault zones. We show that the SOTS kinematics evolved from contraction to left‐lateral strike‐slip during Mississippian‐Pennsylvanian time. Tectonic activity and fault kinematics in the SOTS are explained by a model in which a remnant ocean basin closed diachronously from northeast to southwest along the Southern Appalachian and Ouachita‐Marathon orogens. Left‐lateral transpression during diachronous closure of the Rheic oceanic basin along the Ouachita‐Marathon margin is related to rotation of the stress field due to changes in slab‐pull and eastward drift of Laurentia toward Gondwana. We speculate that the SOTS acted in part as a STEP (subduction‐transfer edge propagator) fault that followed the pre‐existing SOA.

Hobbs, Noah. F.↗
DOE OSTI2024

Testing the optical components for the National Ignition Facility time-resolved soft x-ray opacity spectrometer (OpSpecTR)

Opacity measurements are being carried out at the Z-facility at Sandia National Laboratories and at the National Ignition Facility (NIF) at Lawrence Livermore National Laboratory. The current soft x-ray Opacity Spectrometer (OpSpec) used on the NIF uses two elliptically bent crystals in time-integrated mode on either an image plate or a film. Plans are under way to expand these opacity measurements into a mode of time-resolved detection, called OpSpecTR. Previously, considerations for the available hCMOS detector size and photometrics led to a crystal geometry redesign and the use of a grazing angle x-ray mirror. The mirror acts as a low-pass x-ray energy filter, reducing the contribution of higher energy x rays. The first tests of the mirror and the crystal for OpSpecTR are presented here. The size of the mirror reflection and the reflectivity is tested using a Manson x-ray source. The mirror coupled with the new elliptical crystal shape demonstrates OpSpecTR’s spectral coverage. Finally, the results from the x-ray optics performance testing are shown along with the intended design.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

High order explicit local time stepping methods for hyperbolic conservation laws

In this paper we present and analyze a general framework for constructing high order explicit local time stepping (LTS) methods for hyperbolic conservation laws. In particular, we consider the model problem discretized by Runge-Kutta discontinuous Galerkin (RK-DG) methods and design LTS algorithms based on the strong stability preserving Runge-Kutta (SSP-RK) schemes, that allow spatially variable time step sizes to be used for time integration in different regions of the computational domain. The proposed algorithms are of predictor-corrector type, in which the interface information along the time direction is first predicted based on the SSP-RK approximations and Taylor expansions, and then the fluxes over the region of the interface are corrected to conserve mass exactly at each time step. Following the proposed framework, we detail the corresponding LTS schemes with accuracy up to the fourth order, and prove their conservation property and nonlinear stability for the scalar conservation laws. Numerical experiments are also presented to demonstrate excellent performance of the proposed LTS algorithms.

58 GEOSCIENCES↗
NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithms for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements have been developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration shceme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. The paper presents a description of the Euler solvers along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1990

Implicit flux-split Euler schemes for unsteady aerodynamic analysis involving unstructured dynamic meshes

Improved algorithm for the solution of the time-dependent Euler equations are presented for unsteady aerodynamic analysis involving unstructured dynamic meshes. The improvements were developed recently to the spatial and temporal discretizations used by unstructured grid flow solvers. The spatial discretization involves a flux-split approach which is naturally dissipative and captures shock waves sharply with at most one grid point within the shock structure. The temporal discretization involves an implicit time-integration scheme using a Gauss-Seidel relaxation procedure which is computationally efficient for either steady or unsteady flow problems. For example, very large time steps may be used for rapid convergence to steady state, and the step size for unsteady cases may be selected for temporal accuracy rather than for numerical stability. Steady and unsteady flow results are presented for the NACA 0012 airfoil to demonstrate applications of the new Euler solvers. The unsteady results were obtained for the airfoil pitching harmonically about the quarter chord. The resulting instantaneous pressure distributions and lift and moment coefficients during a cycle of motion compare well with experimental data. A description of the Euler solvers is presented along with results and comparisons which assess the capability.

Batina, John T.↗
NASA NTRS1994

An Approach for Dynamic Grids

An approach is presented for the generation of two-dimensional, structured, dynamic grids. The grid motion may be due to the motion of the boundaries of the computational domain or to the adaptation of the grid to the transient, physical solution. A time-dependent grid is computed through the time integration of the grid speeds which are computed from a system of grid speed equations. The grid speed equations are derived from the time-differentiation of the grid equations so as to ensure that the dynamic grid maintains the desired qualities of the static grid. The grid equations are the Euler-Lagrange equations derived from a variational statement for the grid. The dynamic grid method is demonstrated for a model problem involving boundary motion, an inviscid flow in a converging-diverging nozzle during startup, and a viscous flow over a flat plate with an impinging shock wave. It is shown that the approach is more accurate for transient flows than an approach in which the grid speeds are computed using a finite difference with respect to time of the grid. However, the approach requires significantly more computational effort.

Slater, John W.↗
DOE OSTI2026

Adaptive Control for Load-Following of Boiling Water Reactors Part I: Linear Systems and Fully-Observable Dynamics

Automation control is a key strategy to improve the economic competitiveness of nuclear power plants. Not only does it help reduce operational costs, but it also extends the value proposition of these plants to nontraditional markets, including unattended operations in remote villages and space. However, the dynamics of the operating environments of nuclear reactors are subject to changes over time, and there are no widely adopted methods to ensure that the automation strategy will remain effective over the extended durations required for these applications. Adaptive control is a discipline that offers the possibility to accommodate such changes online. However, it relies on mathematical assumptions that must be respected to ensure robustness and reliability. In this work, we derive an adaptive control formulation for linear systems in which all states are observable and apply it to an instance of load-follow operation for Boiling Water Reactors. We assumed uncertainty in two factors: the temperature coefficient and the control rod worth, both of which are affected over time by the evolution of the nuclear reactor core environment. With an arbitrary penalty factor of 5, we found that the mean absolute and integral time absolute errors can be reduced by more than 90%, underscoring the strength of adaptive control. To extend the application to more challenges, different uncertainties and load-follow trajectories, as well as new formulations that include non-linearity and partial observability, are currently being developed.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1980

On time discretizations for spectral methods

New methods are introduced for the time integration of the Fourier and Chebyshev methods of solution for dynamic differential equations. These methods are unconditionally stable, even though no matrix inversions are required. Time steps are chosen by accuracy requirements alone. For the Fourier method both leapfrog and Runge-Kutta methods are considered. For the Chebyshev method only Runge-Kutta schemes are tested. Numerical calculations are presented to verify the analytic results. Applications to the shallow water equations are presented.

Gottlieb, D.↗
NASA NTRS1991

The effects of rotation on initially anisotropic homogeneous flows

Rotation of initially anisotropic homogeneous flows is studied using a model spectral tensor. It is shown that the anisotropy changes because of the influence of rotation through phase scrambling. Phase scrambling causes the Reynolds stresses to develop with damped oscillations. The final Reynolds stress anisotropy is found to be proportional to the initial structural tensor anisotropy. Closure models for the rapid pressure strain terms should reflect this change in anisotropy, and should drive the anisotropy to reach its final predicted state. Finally, it is shown that long-time integration using direct numerical simulations should be treated with care because phase scrambling effects on a discrete wave space can cause loss of resolution when time becomes large.

Mansour, Nagi N.↗
NASA NTRS1987

A time domain technique for mechanism extraction

The properties of scattered fields from a structure can be better evaluated from the characteristics of the individual scatterers. Decomposition techniques can be classified either as a matrix or an integral formulation. With either formulation, aspect pattern of frequency information of a scattering center can be obtained. Emphasis is placed on an integral (time domain) isolation extraction technique to obtain the frequency characteristics of scattering mechanisms. This technique has its origins in the time domain interpretation of scattered fields.

Dominek, Allen K.↗
DOE OSTI2022

Diagnosing low-mode (ℓ < 6) and mid-mode (6 ≤ ℓ ≤ 60) asymmetries in the post-stagnation phase of laser-direct-drive deuterium–tritium cryogenic implosions on OMEGA

Low- and mid-mode perturbations are possible candidates for performance limitations in cryogenic direct-drive implosions on the OMEGA laser at the Laboratory of Laser Energetics. Simulations with a 3D hydrocode demonstrated that hotspot imagers do not show evidence of the shell breakup in the dense fuel. However, these same simulations revealed that the low- and mid-mode perturbations in the dense fuel could be diagnosed more easily in the post-stagnation phase of the implosion by analyzing the peak in the x-ray emission limb at the coronal–fuel interface than before or at the stagnation phase. In experiments, the asymmetries are inferred from gated images of the x-ray emission of the implosion by using a 16-pinhole array imager filtered to record x-ray energies >800 eV and an x-ray framing camera with 40-ps time integration and 20-μm spatial resolution. Furthermore, a modal analysis is applied to the spatial distribution of the x-ray emission from deuterium and tritium cryogenic implosions on OMEGA recorded after the bang time to diagnose the low- and mid-mode asymmetries, and to study the effect that the beam-to-target ratio (Rb/Rt) has on the shell integrity.

47 OTHER INSTRUMENTATION↗
NASA NTRS1985

Prospects for determining the cosmological helium-3 to helium-4 ratio via absorption-line studies of the local interstellar medium

The prospects for determining the interstellar (He-3)/(He-4) ratio via absorption-line studies in the extreme ultraviolet are considered. A high-resolution extreme-ultraviolet spectrometer fed by a 1-meter-class telescope similar to that discussed in the Far Ultraviolet Spectroscopic Explorer Report (1983) is assumed, and it is found that detection of He-3 may be possible with less than three days of observing time. To measure the (He-3)/(He-4) ratio with sufficient precision to determine the number of light neutrino species will require over 50 days of integration time for each object studied.

Hurwitz, M.↗
NASA NTRS1992

The gravitational wave experiment

Since the optimum size of a gravitational wave detector is the wave length, interplanetary dimensions are needed for the mHz band of interest. Doppler tracking of Ulysses will provide the most sensitive attempt to date at the detection of gravitational waves in the low frequency band. The driving noise source is the fluctuations in the refractive index of interplanetary plasma. This dictates the timing of the experiment to be near solar opposition and sets the target accuracy for the fractional frequency change at 3.0 x 10 exp -14 for integration times of the order of 1000 sec. The instrumentation utilized by the experiment is distributed between the radio systems on the spacecraft and the seven participating ground stations of the Deep Space Network and Medicina. Preliminary analysis is available of the measurements taken during the Ulysses first opposition test.

Bertotti, B.↗
NASA NTRS2004

The Case for Space-Borne Far-Infrared Line Surveys

The combination of sensitive direct detectors and a cooled aperture promises orders of magnitude improvement in the sensitivity and survey time for far-infrared and submillimeter spectroscopy compared to existing or planned capabilities. Continuing advances in direct detector technology enable spectroscopy that approaches the background limit available only from space at these wavelengths. Because the spectral confusion limit is significantly lower than the more familiar spatial confusion limit encountered in imaging applications, spectroscopy can be carried out to comparable depth with a significantly smaller aperture. We are developing a novel waveguide-coupled grating spectrometer that disperses radiation into a wide instantaneous bandwidth with moderate resolution (R ~ 1000) in a compact 2-dimensional format. A line survey instrument coupled to a modest cooled single aperture provides an attractive scientific application for spectroscopy with direct detectors. Using a suite of waveguide spectrometers, we can obtain complete coverage over the entire far-infrared and sub-millimeter. This concept requires no moving parts to modulate the optical signal. Such an instrument would be able to conduct a far-infrared line survey 10 6 times faster than planned capabilities, assuming existing detector technology. However, if historical improvements in bolometer sensitivity continue, so that photon-limited sensitivity is obtained, the integration time can be further reduced by 2 to 4 orders of magnitude, depending on wavelength. The line flux sensitivity would be comparable to ALMA, but at shorter wavelengths and with the continuous coverage needed to extract line fluxes for sources at unknown redshifts. For example, this capability would break the current spectroscopic bottleneck in the study of far-infrared galaxies, the recently discovered, rapidly evolving objects abundant at cosmological distances.

Bock, J. J.↗
NASA NTRS2004

X-Ray and TeV Gamma-Ray Observations of the Strong TeV Gamma-Ray Blazars Mrk 501 and 1ES 1959+650

In June 2003, the blazar 1ES 1959+650 showed a very high level of X-ray activity. The 2-12 keV flux measured with the RXTE All Sky Monitor exceeded 10 mCrab for a prolonged time. We therefore invoked RXTE multiwavelength observations of the source. It was observed with RXTE between May 2,2003 and June 7,2003 with a total integration time of 33 ksec. We complemented the RXTE observations with TeV gamma- ray observations with the Whipple 10 m Cherenkov telescope (Mt. Hopkins, &.). In addition, we invoked radio observations with the UMRAO radio telescope (Aller et al.) and with the Bordeaux optical observatory (P. Charlot et al.). The RXTE observations revealed several X-ray flares with large amplitude variations of the 3 keV - 15 keV spectral index. While the X-ray flares were of intermediate amplitude - some of them comparable to the high X-ray fluxes measured in 2002 - the source remained unexpectedly quiet at TeV gamma-rays. The analysis of the Whipple 10 m TeV gamma-ray data did not result in a firm detection.

Krawczynski, Henric↗