Search NASA⌕ Search

SEARCH · Search NASA

Results for “time integration”

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 253 records · Page 14

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↗

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

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↗

Double phase-lock loop with rapid transient response - A concept

Integration loop, with narrow bandwidth and long response time, integrates phase estimate of sampling loop, resulting in slower but less jittery response to phase step changes. Sampling loop has wide bandwidth and quick response to phase steps of burst sync interval.

Balcewicz, J. F.↗

Digital aspect clock

Digital clock precisely set and reset by pulses from a solar sensor, combined with a logic system, provides accurate time-sector division of spin-stabilized satellite. Integral times for viewing physical phenomena from various directions are equal and mean angles of viewing can be determined.

Bartley, W. C.↗

A small-scale turbulence model

A previously derived analytical model for the small-scale structure of turbulence is reformulated in such a way that the energy spectrum may be computed. The model is an ensemble of two-dimensional (2D) vortices with internal spiral structure, each stretched by an axially symmetric strain flow. Stretching and differential rotation produce an energy cascade to smaller scales in which the stretching represents the effect of instabilities and the spiral structure is the source of dissipation at the end of the cascade. The energy spectrum of the resulting flow may be expressed as a time integration involving only the enstrophy spectrum of the time evolving 2D cross section flow, which may be obtained numerically. Examples are given in which a k exp -5/3 spectrum is obtained by this method. The k exp -5/3 inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the 2D enstrophy spectrum. The results are found to be insensitive to time dependence of the strain rate, including even intermittent on-or-off strains.

Lundgren, T. S.↗

A small-scale turbulence model

A model for the small-scale structure of turbulence is reformulated in such a way that it may be conveniently computed. The model is an ensemble of randomly oriented structured two dimensional vortices stretched by an axially symmetric strain flow. The energy spectrum of the resulting flow may be expressed as a time integral involving only the enstrophy spectrum of the time evolving two-dimensional cross section flow, which may be obtained numerically. Examples are given in which a k(exp -5/3) spectrum is obtained by this method without using large wave number asymptotic analysis. The k(exp -5/3) inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the two-dimensional enstrophy spectrum. The results are insensitive to time dependence of the strain-rate, including even intermittent on-or-off strains.

Lundgren, T. S.↗

Local-in-Time Adjoint-Based Method for Optimal Control/Design Optimization of Unsteady Compressible Flows

.We study local-in-time adjoint-based methods for minimization of ow matching functionals subject to the 2-D unsteady compressible Euler equations. The key idea of the local-in-time method is to construct a very accurate approximation of the global-in-time adjoint equations and the corresponding sensitivity derivative by using only local information available on each time subinterval. In contrast to conventional time-dependent adjoint-based optimization methods which require backward-in-time integration of the adjoint equations over the entire time interval, the local-in-time method solves local adjoint equations sequentially over each time subinterval. Since each subinterval contains relatively few time steps, the storage cost of the local-in-time method is much lower than that of the global adjoint formulation, thus making the time-dependent optimization feasible for practical applications. The paper presents a detailed comparison of the local- and global-in-time adjoint-based methods for minimization of a tracking functional governed by the Euler equations describing the ow around a circular bump. Our numerical results show that the local-in-time method converges to the same optimal solution obtained with the global counterpart, while drastically reducing the memory cost as compared to the global-in-time adjoint formulation.

Yamaleev, N. K.↗

Visible-Infrared Hyperspectral Image Projector

The VisIR HIP generates spatially-spectrally complex scenes. The generated scenes simulate real-world targets viewed by various remote sensing instruments. The VisIR HIP consists of two subsystems: a spectral engine and a spatial engine. The spectral engine generates spectrally complex uniform illumination that spans the wavelength range between 380 nm and 1,600 nm. The spatial engine generates two-dimensional gray-scale scenes. When combined, the two engines are capable of producing two-dimensional scenes with a unique spectrum at each pixel. The VisIR HIP can be used to calibrate any spectrally sensitive remote-sensing instrument. Tests were conducted on the Wide-field Imaging Interferometer Testbed at NASA s Goddard Space Flight Center. The device is a variation of the calibrated hyperspectral image projector developed by the National Institute of Standards and Technology in Gaithersburg, MD. It uses Gooch & Housego Visible and Infrared OL490 Agile Light Sources to generate arbitrary spectra. The two light sources are coupled to a digital light processing (DLP(TradeMark)) digital mirror device (DMD) that serves as the spatial engine. Scenes are displayed on the DMD synchronously with desired spectrum. Scene/spectrum combinations are displayed in rapid succession, over time intervals that are short compared to the integration time of the system under test.

Bolcar, Matthew↗

High-Order Implicit-Explicit Multi-Block Time-stepping Method for Hyperbolic PDEs

This work seeks to explore and improve the current time-stepping schemes used in computational fluid dynamics (CFD) in order to reduce overall computational time. A high-order scheme has been developed using a combination of implicit and explicit (IMEX) time-stepping Runge-Kutta (RK) schemes which increases numerical stability with respect to the time step size, resulting in decreased computational time. The IMEX scheme alone does not yield the desired increase in numerical stability, but when used in conjunction with an overlapping partitioned (multi-block) domain significant increase in stability is observed. To show this, the Overlapping-Partition IMEX (OP IMEX) scheme is applied to both one-dimensional (1D) and two-dimensional (2D) problems, the nonlinear viscous Burger's equation and 2D advection equation, respectively. The method uses two different summation by parts (SBP) derivative approximations, second-order and fourth-order accurate. The Dirichlet boundary conditions are imposed using the Simultaneous Approximation Term (SAT) penalty method. The 6-stage additive Runge-Kutta IMEX time integration schemes are fourth-order accurate in time. An increase in numerical stability 65 times greater than the fully explicit scheme is demonstrated to be achievable with the OP IMEX method applied to 1D Burger's equation. Results from the 2D, purely convective, advection equation show stability increases on the order of 10 times the explicit scheme using the OP IMEX method. Also, the domain partitioning method in this work shows potential for breaking the computational domain into manageable sizes such that implicit solutions for full three-dimensional CFD simulations can be computed using direct solving methods rather than the standard iterative methods currently used.

Nielsen, Tanner B.↗

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solar System Studies with the Space Infrared Telescope Facility (SIRTF)

SIRTF (Space Infrared Telescope Facility) is the final element in NASA's 'Great Observatories' program. It consists of an 85-cm cryogenically-cooled observatory for infrared astronomy from space. SIRTF is scheduled for launch in late 2001 or early 2002 on a Delta rocket into a heliocentric orbit trailing the Earth. Data from SIRTF will be processed and disseminated to the community through the SIRTF Science Center (SSC) located at the Infrared Processing and Analysis Center (IPAC) at Caltech. Some 80/% of the total observing time (estimated at a minimum of 7500 hours of integration time per year for the mission lifetime of about 4 years) will be available to the scientific community at large through a system of refereed proposals. Three basic instruments are located in the SIRTF focal plane. The Multiband Imaging Photometer (MIPS), the Infrared Array Camera (IRAC), and the Infrared Spectrometer (IRS), taken together, provide imaging and spectroscopy from 3.5 to 160 microns. Among the solar system studies suited to SIRTF are the following: 1) spectroscopy and radiometry of small bodies from the asteroid main belt, through the Trojan clouds, to the Kuiper Disk; 2) dust distribution in the zodiacal cloud and the Earth's heliocentric dust ring; 3) spectroscopy and radiometry of comets; and 4) spectroscopy and radiometry of planets and their satellites. Searches for, and studies of dust disks around other stars, brown dwarfs, and superplanets will also be conducted with SIRTF. The SORTIE web site (http://ssc.ipac.caltech.edu/sirtf) contains important details and documentation on the project, the spacecraft, the telescope, instruments, and observing procedures. A community-wide workshop for solar system studies with SIRTF is in the planning stages by the author and Martha S. Hanner for the summer of 1999.

Cruikshank, Dale P.↗

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗

Acoustic-optic spectrometer. 1: Noise contributions and system consideration

An acousto-optic spectrometer (AOS) used as an IF spectrometer to a heterodyne receiver is modeled as a total power multi-channel integrating receiver. Systematic noise contributions common to all total power, time integrating receivers, as well as noise terms unique to the use of optical elements and photo-detectors in an AOS are identified and discussed. In addition, degradation of signal-to-noise ratio of an unbalanced Dicke receiver compared to a balanced Dicke receiver is found to be due to gain calibration processing and is not an instrumental effect.

Chin, G.↗

Development of a Complete Landsat Evapotranspiration and Energy Balance Archive to Support Agricultural Consumptive Water Use Reporting and Prediction in the Central Valley, CA

Mapping evapotranspiration (ET) from agricultural areas in Californias Central Valley is critical for understanding historical consumptive use of surface and groundwater. In addition, long histories of ET maps provide valuable training information for predictive studies of surface and groundwater demands. During times of drought, groundwater is commonly pumped to supplement reduced surface water supplies in the Central Valley. Due to the lack of extensive groundwater pumping records, mapping consumptive use using satellite imagery is an efficient and robust way for estimating agricultural consumptive use and assessing drought impacts. To this end, we have developed and implemented an algorithm for automated calibration of the METRIC remotely sensed surface energy balance model on NASAs Earth Exchange (NEX) to estimate ET at the field scale. Using automated calibration techniques on the NEX has allowed for the creation of spatially explicit historical ET estimates for the Landsat archive dating from 1984 to the near present. Further, our use of spatial NLDAS and CIMIS weather data, and spatial soil water balance simulations within the NEX METRIC workflow, has helped overcome challenges of time integration between satellite image dates. This historical and near present time archive of agricultural water consumption for the Central Valley will be an extremely useful dataset for water use and drought impact reporting, and predictive analyses of groundwater demands.

valleys↗

Public Health Data Applications Using the CDC Tracking Network: Augmenting Environmental Hazard Information with Lower-latency NASA Data

Exposure to environmental hazards is an important determinant of health, and the frequency and severity of exposures is expected to be impacted by climate change. Through a partnership with the U.S. National Aeronautics and Space Administration, the U.S. Centers for Disease Control and Prevention’s National Environmental Public Health Tracking Network is integrating timely observations and model data of priority environmental hazards into its publicly accessible Data Explorer (https://ephtracking.cdc.gov/DataExplorer/). Newly integrated datasets over the contiguous U.S. (CONUS) include: daily 5-day forecasts of air quality based on the Goddard Earth Observing System Composition Forecast (GEOS-CF), daily historical (1980-present) concentrations of speciated PM2.5 based on the Modern Era Retrospective analysis for Research and Applications, version 2 (MERRA-2), and Moderate Resolution Imaging Spectroradiometer (MODIS) daily near real-time maps of flooding (MCDWD). Data integrated into the CDC Tracking Network are broadly intended to improve community health through action by informing both research and early warning activities, including (1) describing temporal and spatial trends in disease and potential environmental exposures, (2) identifying populations most affected, (3) generating hypotheses about associations between health and environmental exposures, and (4) developing, guiding, and assessing environmental public health policies and interventions aimed at reducing or eliminating health outcomes associated with environmental factors.

air quality↗

Stabilization by modification of the Lagrangian

In order to reduce the error growth during a numerical integration, a method for stabilizing the differential equations of Keplerian motion is offered. It is characterized by the use of the eccentric anomaly as independent variable in such a way that the time transformation is given by a generalized Lagrange formalism. The control terms in the equations of motion obtained by this modified Lagrangian give immediately a completely Lyapunov-stable set of differential equations. The equation of time integration is modified by a control term which leads to an integral which defined the time element for the perturbed Keplerian motion.

Baumgarte, J. W.↗

A numerical transport scheme which avoids negative mixing ratios

A new scheme for numerically integrating the transport equation has been developed. This square root method avoids the problem of negative mixing ratios by using the square root of the concentration instead of the concentration itself as an advective time step. Conservation of total mass at every time step is achieved by extending the concept of quadratic conservation to the numerical time integration method. A few time schemes that fulfill the requirement of step-by-step quadratic conservation are described, the simplest of which is the modified Lax-Wendroff method.

Schneider, H.-R.↗