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 649 records · Page 36

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↗

Measurement-driven large-eddy simulations of a diurnal cycle during a wake-steering field campaign

Abstract. High-fidelity flow modeling with data assimilation enables accurate representation of the wind farm operating environment under realistic, nonstationary atmospheric conditions. Capturing the temporal evolution of the turbulent atmospheric boundary layer is critical to understanding the behavior of wind turbines under operating conditions with simultaneously varying inflow and control inputs. This paper has three parts: the identification of a case study during a field evaluation of wake steering; the development of a tailored mesoscale-to-microscale coupling strategy that resolved local flow conditions within a large-eddy simulation (LES), using observations that did not completely capture the wind and temperature fields throughout the simulation domain; and the application of this coupling strategy to validate high-fidelity aeroelastic predictions of turbine performance and wake interactions with and without wake steering. The case study spans 4.5 h after midnight local time, during which wake steering was toggled on and off five times, achieving yaw offset angles ranging from 0 to 17°. To resolve nonstationary nighttime conditions that exhibited shear instabilities, the turbulence field was evolved starting from the diurnal cycle of the previous day. These background conditions were then used to drive wind farm simulations with two different models: an LES with actuator disk turbines and a steady-state engineering wake model. Subsequent analysis identified two representative periods during which the up- and downstream turbines were most nearly aligned with the mean wind direction and had observed yaw offsets of 0 and 15°. Both periods corresponded to partial waking on the downstream turbine, which had errors in the LES-predicted power of 4 % and 6 %, with and without wake steering. The LES was also able to capture conditions during which an upstream turbine wake induced a speedup at a downstream turbine and increased power production by up to 13 %.

17 WIND ENERGY↗

High-power ferro-electric fast reactive tuner

We present a novel design of a FerroElectric Fast Reactive Tuner (FE-FRT) capable of modulating Mega VAR reactive power on a sub-microsecond time scale. We show detailed analytical estimates of the performance of this device and benchmark these estimates against finite element method eigenmode and frequency domain electromagnetic simulations.

43 PARTICLE ACCELERATORS↗

Electron and Electron-Nuclear Dynamics of Open Quantum Subsystems

Leveraging an open-subsystem formulation of Density Functional Theory (DFT) we aim at describing periodic and molecular systems alike, including their electronic and nuclear dynamics. Subsystem DFT enables first principles simulations to approach realistic time- and length-scales, and most importantly sheds light on the dynamical behavior of complex systems. Taking subsystem DFT to the time domain allows us to inspect the electron dynamics of condensed-phase systems in real time. In liquids and interfaces, we observe all the relevant regimes proper of non-Markovian open quantum system dynamics, such as electronic energy transfer, and screening. In addition, the ab-initio modeling of system-bath interactions brought us to observe and justify the holographic time-dependent electron density theorem. Contrary to interactions between molecular (finite) systems, when molecules interact with metal or semiconductor surfaces the electron dynamics is strongly non-Markovian with dramatic repercussions to the molecule’s response to external perturbations. Metals and semiconductors typically have large polarizabilities, and even in a regime of low coupling their effect on impinging molecular species is significant–line broadening, peak shift, and intensity borrowing are observed, characterized, and explained in terms of inter-subsystem dynamical interactions and a many-body decomposition of the system’s density-density response function in a way that transcends the canons of Fermi Golden Rule.

33 ADVANCED PROPULSION SYSTEMS↗

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid-flow of a row of jets in crossflow - A numerical study

A detailed computer-visualized flow field of a row of jets in a confined crossflow is presented. The Reynolds averaged Navier-Stokes equations are solved using a finite volume method that incorporates a partial differential equation for incremental pressure to obtain a divergence-free flow field. The turbulence is described by a multiple-time-scale turbulence model. The computational domain includes the upstream region of the circular jet so that the interaction between the jet and the crossflow is simulated accurately. It is shown that the row of jets in the crossflow is characterized by a highly complex flow field that includes a horse-shoe vortex and two helical vortices whose secondary velocity components are co-rotating in space. It is also shown that the horse-shoe vortex is a ring of reversed flows located along the circumference of the jet exit.

Kim, S.-W.↗

An Analysis of Performance Enhancement Techniques for Overset Grid Applications

The overset grid methodology has significantly reduced time-to-solution of high-fidelity computational fluid dynamics (CFD) simulations about complex aerospace configurations. The solution process resolves the geometrical complexity of the problem domain by using separately generated but overlapping structured discretization grids that periodically exchange information through interpolation. However, high performance computations of such large-scale realistic applications must be handled efficiently on state-of-the-art parallel supercomputers. This paper analyzes the effects of various performance enhancement techniques on the parallel efficiency of an overset grid Navier-Stokes CFD application running on an SGI Origin2000 machine. Specifically, the role of asynchronous communication, grid splitting, and grid grouping strategies are presented and discussed. Results indicate that performance depends critically on the level of latency hiding and the quality of load balancing across the processors.

Djomehri, J. J.↗