Search NASA⌕ Search

SEARCH · Search NASA

Results for “time domain simulation”

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 631 records · Page 35

Atmospheric Carbon Cycle Dynamics Over the Above Domain: an Integrated Analysis Using Aircraft Observations (Arctic-CA) and Model Simulations (GEOS)

The Arctic Carbon Atmospheric Profiles (Arctic-CAP) project conducted six airborne surveys of Alaska and northwestern Canada between April and November 2017 to capture the spatial and temporal gradients of northern high-latitude carbon dioxide (CO2), methane (CH4) and carbon monoxide (CO) as part of NASA's Arctic-Boreal Vulnerability Experiment (ABoVE). The Arctic-CAP sampling strategy involved acquiring vertical profiles of CO2, CH4 and CO from the surface to 5 km altitude at 25 sites around the ABoVE domain on a 4- to 6-week time interval. We observed vertical gradients of CO2, CH4 and CO that vary by eco-region and duration of the sampling period, which spanned the majority of the seasonal cycle. All Arctic-CAP measurements were compared to a global simulation using the Goddard Earth Observing System (GEOS) modeling system. Comparisons with GEOS simulations of atmospheric CO2, CH4 and CO highlight the potential of these multi-species observations to inform improvements in surface flux estimates and the representation of atmospheric transport. GEOS simulations provide estimates of the near surface average CO2 and CH4 enhancements that are well correlated with aircraft observations (R=0.74 and R=0.60 respectively), suggesting that GEOS has reasonable fidelity over this complex and heterogeneous region. This model-data comparison over the ABoVE domain reveals that while current state-of-the-art models and flux estimates are able to capture broadscale spatial and temporal patterns in near-surface CO2 and CH4 concentrations, more work is needed to resolve fine-scale flux features that are observed. The study also provides a framework for benchmarking a global model at regional scales, which is needed to use climate models as tools to investigate high-latitude carbon-climate feedbacks.

C Sweeney↗

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗

Learning a general model of single phase flow in complex 3D porous media

Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.

36 MATERIALS SCIENCE↗

Learning a General Model of Single Phase Flow in Complex 3D Porous Media

Modeling effective transport properties of 3D porous media, such as permeability, at multiple scales is challenging as a result of the combined complexity of the pore structures and fluid physics—in particular, confinement effects which vary across the nanoscale to the microscale. While numerical simulation is possible, the computational cost is prohibitive for realistic domains, which are large and complex. Although machine learning (ML) models have been proposed to circumvent simulation, none so far has simultaneously accounted for heterogeneous 3D structures, fluid confinement effects, and multiple simulation resolutions. By utilizing numerous computer science techniques to improve the scalability of training, we have for the first time developed a general flow model that accounts for the pore-structure and corresponding physical phenomena at scales from Angstrom to the micrometer. Using synthetic computational domains for training, our ML model exhibits strong performance (R 2 = 0.9) when tested on extremely diverse real domains at multiple scales.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

(abstract) 3D Electromagnetic Plasma Particle Simulations

A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.

electromagnetic plasma particles simulations 3D al↗

The Impact of Dry Midlevel Air on Hurricane Intensity in Idealized Simulations with No Mean Flow

This study examines the potential negative influences of dry midlevel air on the development of tropical cyclones (specifically, its role in enhancing cold downdraft activity and suppressing storm development). The Weather Research and Forecasting model is used to construct two sets of idealized simulations of hurricane development in environments with different configurations of dry air. The first set of simulations begins with dry air located north of the vortex center by distances ranging from 0 to 270 km, whereas the second set of simulations begins with dry air completely surrounding the vortex, but with moist envelopes in the vortex core ranging in size from 0 to 150 km in radius. No impact of the dry air is seen for dry layers located more than 270 km north of the initial vortex center (approximately 3 times the initial radius of maximum wind). When the dry air is initially closer to the vortex center, it suppresses convective development where it entrains into the storm circulation, leading to increasingly asymmetric convection and slower storm development. The presence of dry air throughout the domain, including the vortex center, substantially slows storm development. However, the presence of a moist envelope around the vortex center eliminates the deleterious impact on storm intensity. Instead, storm size is significantly reduced. The simulations suggest that dry air slows intensification only when it is located very close to the vortex core at early times. When it does slow storm development, it does so primarily by inducing outward- moving convective asymmetries that temporarily shift latent heating radially outward away from the high-vorticity inner core.

Weather Research and Forecasting model↗

On the use of a multigrid-reduction-in-time algorithm for multiscale convergence of turbulence simulations

Simulations of turbulent flow present challenges in terms of accuracy and affordability on modern highly-parallel computer architectures. A multigrid-reduction-in-time algorithm is used to provide a framework for separately evolving different scales of turbulence and for parallelizing the temporal domain, thereby increasing the concurrency. It is hypothesized that the space–time locality of the small scales of turbulence can be used to circumvent difficulties in applying temporal multigrid to flows dominated by inertial physics. For algorithms that fall well short of spectral accuracy (fourth-order is used in this work) attention must be paid to the accuracy of features on scales transferred between multigrid levels. Numerical experiments were performed using implicit large-eddy simulation. Results from applying the approach to an infinite-Reynolds number Taylor–Green flow and a double-shear flow at a Reynolds number of 11650 provide strong evidence that the approach has merit. The multigrid-reduction-in-time framework can be used to parallelize the temporal domain of a high-Reynolds-number turbulent flow and permit independent convergence of different scales. Establishing this foundation allows for future research in reducing the wall-clock time to solve turbulent flows while retaining the same accuracy as sequential solvers. In conclusion, current performance results from parallelizing the temporal domain are not competitive with those from sequential-in-time methods.

97 MATHEMATICS AND COMPUTING↗

GT-MSOCC - A domain for research on human-computer interaction and decision aiding in supervisory control systems

The Georgia Tech-Multisatellite Operations Control Center (GT-MSOCC), a real-time interactive simulation of the operator interface to a NASA ground control system for unmanned earth-orbiting satellites, is described. The GT-MSOCC program for investigating a range of modeling, decision aiding, and workstation design issues related to the human-computer interaction is discussed. A GT-MSOCC operator function model is described in which operator actions, both cognitive and manual, are represented as the lowest level discrete control network nodes, and operator action nodes are linked to information needs or system reconfiguration commands.

Mitchell, Christine M.↗

Cirrus cloud model parameterizations: Incorporating realistic ice particle generation

Recent cirrus cloud modeling studies have involved the application of a time-dependent, two dimensional Eulerian model, with generalized cloud microphysical parameterizations drawn from experimental findings. For computing the ice versus vapor phase changes, the ice mass content is linked to the maintenance of a relative humidity with respect to ice (RHI) of 105 percent; ice growth occurs both with regard to the introduction of new particles and the growth of existing particles. In a simplified cloud model designed to investigate the basic role of various physical processes in the growth and maintenance of cirrus clouds, these parametric relations are justifiable. In comparison, the one dimensional cloud microphysical model recently applied to evaluating the nucleation and growth of ice crystals in cirrus clouds explicitly treated populations of haze and cloud droplets, and ice crystals. Although these two modeling approaches are clearly incompatible, the goal of the present numerical study is to develop a parametric treatment of new ice particle generation, on the basis of detailed microphysical model findings, for incorporation into improved cirrus growth models. For example, the relation between temperature and the relative humidity required to generate ice crystals from ammonium sulfate haze droplets, whose probability of freezing through the homogeneous nucleation mode are a combined function of time and droplet molality, volume, and temperature. As an example of this approach, the results of cloud microphysical simulations are presented showing the rather narrow domain in the temperature/humidity field where new ice crystals can be generated. The microphysical simulations point out the need for detailed CCN studies at cirrus altitudes and haze droplet measurements within cirrus clouds, but also suggest that a relatively simple treatment of ice particle generation, which includes cloud chemistry, can be incorporated into cirrus cloud growth.

Sassen, Kenneth↗

Evaluation of Lightning Induced Effects in a Graphite Composite Fairing Structure

Defining the electromagnetic environment inside a graphite composite fairing due to near-by lightning strikes is of interest to spacecraft developers. This effort develops a transmission-line-matrix (TLM) model with a CST Microstripes to examine induced voltages. on interior wire loops in a composite fairing due to a simulated near-by lightning strike. A physical vehicle-like composite fairing test fixture is constructed to anchor a TLM model in the time domain and a FEKO method of moments model in the frequency domain. Results show that a typical graphite composite fairing provides adequate shielding resulting in a significant reduction in induced voltages on high impedance circuits despite minimal attenuation of peak magnetic fields propagating through space in near-by lightning strike conditions.

Trout, Dawn H.↗

Computational Bayesian Methods Applied to Complex Problems in Bio and Astro Statistics

In this dissertation we apply computational Bayesian methods to three distinct problems. In the first chapter, we address the issue of unrealistic covariance matrices used to estimate collision probabilities. We model covariance matrices with a Bayesian Normal-Inverse-Wishart model, which we fit with Gibbs sampling. In the second chapter, we are interested in determining the sample sizes necessary to achieve a particular interval width and establish non-inferiority in the analysis of prevalences using two fallible tests. To this end, we use a third order asymptotic approximation. In the third chapter, we wish to synthesize evidence across multiple domains in measurements taken longitudinally across time, featuring a substantial amount of structurally missing data, and fit the model with Hamiltonian Monte Carlo in a simulation to analyze how estimates of a parameter of interest change across sample sizes.

Elrod, Chris↗

Lamb wave defect detection and evaluation using a fully non-contact laser system

Traditional Lamb wave structural health monitoring (SHM)/nondestructive evaluation (NDE) system employs contact type transducers such as PZT, ultrasonic transducers, and optical fibers. In application, transducer attachment and maintenance can be time and labor consuming. In addition, the use of couplant and adhesives can introduce additional materials on structures, and the interface coupling is often not well understood. To overcome these limitations, we proposed a fully non-contact NDE system by employing pulsed laser (PL) for Lamb wave actuation and scanning laser Doppler vibrometer (SLDV) for Lamb wave sensing. The proposed system is implemented on aluminum plates. The PL Lamb wave excitation is calibrated, and the optimal parameters are obtained. Lamb wave modes are then characterized through 1D wavefield analysis. With the calibrated and characterized system, defect detection and evaluation are achieved on aluminum plates with simulated defects (surfaced-bonded quartz rod, and machine milled crack) through 1D and 2D inspection in both time-space and frequency-wavenumber domains.

Ma, Zhaoyun↗

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↗

(Paper) TRIPWIRE: Multi-modal Distributed Sensing For Repository Verification

Underground geological repositories for storing used nuclear fuel are being planned and developed around the world. Research is being performed to understand technical aspects of sealing entombed materials and underground repositories to ensure the radiological and nuclear materials remain contained. To support safeguards verification for these facilities research at Idaho National Laboratory (INL) is underway to develop and demonstrate a multi-modal sensor system, TRIPWIRE, for containment verification in inaccessible radiological and nuclear waste repositories. The TRIPWIRE system will continuously monitor ionizing radiation and electromagnetic fields in the vicinity of emplaced nuclear materials buried in a repository, reporting on disturbances with a real-time alarm control station. The system will use long-length scintillating fiber bundles (SFBs) to perform area radiation monitoring; these will be coupled to kilometer-scale multimodal optical communication fibers – all light sensors and electronic components used with this system will be located above ground. Electromagnetic fields, and changes in local dielectric conditions caused by intrusion and soil movement, will be monitored using commercial grade, ported "leaky" coaxial cables (PCCs), with control electronics also located above ground. Tamper-indicating self-diagnostic assessments will be done using optical and electronic time domain reflectometry in the SFBs and PCCs, respectively. The result will be a long, kilometer-scale multi-modal SFB-PCC system. Simulation and modeling are being used to inform the work and a demonstration of the system's utility is planned at a nuclear storage facility in the future.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗

A technique for simulating turbulence for aerospace vehicle flight simulation studies

An atmospheric turbulence model which accommodates variability of turbulence properties along an aerospace vehicle trajectory was developed. The technique involves the use of Dryden spectral forms in which the defining parameters are the standard deviations (sigma) and integral scales (L) of turbulence. These spectra are expressed as nondimensional functions of the nondimensional frequency Omega = omega L/V where omega is dimensional radian frequency and V is the true air speed of the aerospace vehicle. The nondimensional spectra are factored by standard techniques to obtain nondimensional linear recursive filters in the time domain whereby band-limited white-like noise can be operated upon to obtain nondimensional longitudinal, lateral, and vertical turbulence velocities, as functions of nondimensional time, tV/L, where t is time. Application of the technique to the simulation of the space shuttle orbiter entry flight phase is discussed.

Fichtl, G. H.↗

Validation of the dynamic response of a blade-element UH-60 simulation model in hovering flight

The dynamic fidelity of an operational blade-element simulation model of the UH-60A helicopter is assessed for the hovering and low-speed flight regimes. Nonparametric frequency-response identification techniques and time-history comparisons are used to determine the validity of vehicle responses. Frequency-domain methods are applied to the model to isolate deficiencies and verify refinements and corrections. A dedicated flight-test program was conducted to provide data used in the analyses. Model deficiencies are also identified by using side-by-side pilot assessments of a motion-based simulation and of a test aircraft. Model refinements are found to improve fidelity significantly in the frequency range of interest to handling-qualities research. The applicability of the model to high-bandwidth flight-control research is also discussed.

Ballin, Mark G.↗

Skyrmion lattice formation and destruction mechanisms probed with TR-SANS

Magnetic skyrmions are topologically protected, nanoscale whirls of the spin configuration that tend to form hexagonally ordered arrays. As a topologically non-trivial structure, the nucleation and annihilation of the skyrmion, as well as the interaction between skyrmions, varies from conventional magnetic systems. Recent works have suggested that the ordering kinetics in these materials occur over millisecond or longer timescales, which is unusually slow for magnetic dynamics. The current work investigates the skyrmion ordering kinetics, particularly during lattice formation and destruction, using time-resolved small angle neutron scattering (TR-SANS). Evaluating the time-resolved structure and intensity of the neutron diffraction pattern reveals the evolving real-space structure of the skyrmion lattice and the timeframe of the formation. Measurements were performed on three prototypical skyrmion materials: MnSi, (Fe,Co)Si, and Cu 2 OSeO 3 . To probe lattice formation and destruction kinetics, the systems were prepared in the stable skyrmion state, and then a square-wave magnetic field modulation was applied. The measurements show that the skyrmions quickly form ordered domains, with a significant distribution in lattice parameters, which then converge to the final structure; the results confirm the slow kinetics, with formation times between 10 ms and 99 ms. Comparisons are made between the measured formation times and the fundamental material properties, suggesting the ordering temperature, saturation magnetization and magnetocrystalline anisotropy may be driving the timeframes. Micromagnetic simulations were also performed and support a scaling of the kinetics with sample volume, a behavior which is caused by the reconciling of misaligned domains.

Liyanage, W. L. N. C.↗

A synthetic Roman Space Telescope High-Latitude Time-Domain Survey: supernovae in the deep field

ABSTRACT NASA will launch the Nancy Grace Roman Space Telescope (hereafter Roman) in the second half of this decade, which will allow for a generation-defining measurement of dark energy through multiple probes, including Type Ia supernovae (SNe Ia). To improve decisions on survey strategy, we have created the first simulations of realistic Roman images that include artificial SNe Ia injected as point sources in the images. Our analysis combines work done on Roman simulations for weak gravitational lensing studies as well as catalogue-level simulations of SN Ia samples. We have created a time series of images over 2 yr containing ∼1050 SNe Ia, covering a 1 deg2 subarea of a planned 5 deg2 deep survey. We have released these images publicly for community use along with input catalogues of all injected sources. We create secondary products from these images by generating coadded images and demonstrating recovery of transient sources using image subtraction. We perform first-use analyses on these images in order to measure galaxy detection efficiency, point source detection efficiency, and host-galaxy association biases. The simulated images can be found here at https://roman.ipac.caltech.edu/sims/SN_Survey_Image_sim.html.

Astronomy & Astrophysics↗