Search NASA⌕ Search

SEARCH · Search NASA

Results for “fast time 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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Investigating the role of edge neutrals in exciting tearing mode activity and achieving flat temperature profiles in LTX-β

Abstract We present observations, numerical simulations, and analysis from experiments in the Lithium Tokamak Experiment-Beta (LTX- β ) in which the electron temperature profile ( T e ( r )) shifts from flat to peaked and a tearing mode is also destabilized when the average density ( n e ave ) exceeds ∼10 19 m −3 . Flat T e ( r ) is obtained routinely in LTX- β , with a lithium coated, low-recycling first wall, once the external fueling is stopped and density decays [Boyle et al 2023 Nucl. Fusion 63 056020]. In the present experiment, flat T e profiles can be sustained while maintaining constant n e ave below a line averaged density threshold ( n e ave th ) of ∼10 19 m −3 . Above n e ave th , T e ( r ) shifts from flat to peaked and a tearing mode is destabilized. Due to low recycling, the achieved n e ave can be controlled precisely by external fueling and hence, a certain threshold of the edge neutral inventory from the external fueling is experimentally manifested through n e ave th . The goal of the present work is to investigate the role of edge neutrals in determining T e ( r ) and MHD stability in the unique low-recycling regime of LTX- β . Our hypothesis is that the peaking of T e ( r ) beyond n e ave th is due ultimately to the edge cooling by the cold neutrals beyond a critical fueling flux. At lower fueling flux, flat T e ( r ) results in broader pressure profile and lower resistivity, which in turn stabilizes the tearing mode. This hypothesis is supported by edge neutral density estimation by DEGAS 2 code. Mode analysis by singular value decomposition confirms the tearing mode structure to be m / n = 2/1 ( m and n being the poloidal and toroidal mode numbers). Linear tearing stability analysis with M3D-C1 predicts that plasmas with n e ave > 10 19 are highly susceptible to a n = 1 tearing mode. ORBIT simulations, however, confirmed that the tearing modes do not contribute to the loss of fast ions from neutral beam injection. This study shows for the first time that the neutral inventory at the edge could be one of the deciding factors for the achievability of the unique operation regime of flat T e ( r ) and the excitation of tearing activity that could be disruptive for the plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

First Experiments and Radiographs on the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) Dense Plasma Focus

We report a dense plasma focus (DPF) is a relatively compact coaxial plasma gun, which completes its discharge as a Z-pinch. These devices are designed to operate at a variety of scales to produce short (<100 ns) pulses of ions, X-rays, and/or neutrons. LLNL recently constructed and brought into operation a new device, the MegaJOuLe Neutron Imaging Radiography (MJOLNIR) DPF, which is designed for radiography and high-yield operations. This device has been commissioned and has achieved neutron yields of up to 3.8E11 neutrons/pulse at 2.5-MA peak current while operating at up to 1 MJ of stored energy in its original pulsed power configuration. MJOLNIR is equipped with a wide range of diagnostics, including nuclear activation detectors, neutron time-of-flight (nToF) detectors, a fast-framing camera, optical light gates, and a time-gated neutron and X-ray imager. LLNL also runs unique particle-in-cell (PIC) simulations of DPF discharges in the Chicago code and has gained significant insight into the various physical factors that influence neutron yield. MJOLNIR is one of the first DPFs whose design and continual upgrades are heavily influenced by model predictions. In this article, we describe insights from modeling, device operation, and recent results. Comparisons between modeling predictions and measurements, as well as X-ray and neutron images are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Deep learning-assisted modeling for χ (2) nonlinear optics

Modeling second-order (χ(2)) nonlinear optical processes remains computationally expensive due to the need to resolve fast field oscillations and simulate wave propagation using methods such as the split-step Fourier method (SSFM). This can become a bottleneck in real-time applications, such as high-repetition-rate laser systems requiring rapid feedback and control. We present a long short-term memory-based surrogate model trained on SSFM simulations generated from a start-to-end model of the photocathode drive laser at SLAC National Accelerator Laboratory’s Linac Coherent Light Source II. The model achieves over 250× speedup while maintaining high fidelity, enabling future real-time optimization and laying the foundation for data-integrated modeling frameworks and digital twins of laser systems.

Accelerator Physics (physics.acc-ph)↗
NASA NTRS2005

Role of Probabilistic Micromechanics Modeling in Establishing Design Allowables in Composites

One of the major challenges in designing with any new material, and particularly with advanced composite materials, is the fidelity of material design allowables. In the case of composite materials, the concern arises from the inherent nature of these materials, i.e., their heterogeneous make-up and the various factors that affect their properties in a specific design environment. Composites have various scales - micro, macro, laminate and structural, as well as numerous other fabrication related parameters. Many advanced composites in aerospace applications involve complex two- and three-dimensional fiber architectures and requires high-temperature processing. Since there are uncertainties associated with each of these, the observed behavior of composite materials shows scatter. Evaluating the effect of each of these variables on the observed scatter in composite properties solely by teSting is cost and time prohibitive. One alternative is to evaluate these effects by computational simulation. The authors have developed probabilistic composite micromechanics techniques by combining woven composite micromechanics and Fast Probability Integration (FPI) techniques to address these issues. In this paper these techniques will be described and demonstrated through selected examples. Results in the form of cumulative distribution functions (CDF) of the composite properties of a MI (melt-infiltrated) SiC/SiC (silicon carbide fiber in a silicon carbide matrix) Composite will be presented. A CDF is a relationship defined by the value of the property (the response variable) with respect to the cumulative probability of occurrence. Furthermore, input variables causing scatter are identified and ranked based upon their sensitivity magnitude. Sensitivity information is very valuable in quality control. How these results can be utilized to develop design allowables so that these materials may be used by structural analysts/designers will also be discussed.

Mital, Subodh K.↗
NASA NTRSDate not supplied

A Heat and Mass Transfer Model of the Orion European Service Module Propulsion Sub-System

A tool for simulating the thermodynamics of the fluid components of the Orion European Service Module propulsion subsystem has been developed. The 1-dimensional,time-stepping method tracks energy and mass movement through the propellant and pressurant tanks during burn and coast phases of a simulated mission. The code is fast and flexible, allowing easy modifications and incorporation into other computational tools. The results from the simulation have been used to support mission planning, development of flight rules and to guide experimental testing of components and systems.

Spacecraft propulsion↗
NASA NTRSDate not supplied

Probability of Obstacle Collision for UAVs in Presence of Wind

For incorporation of unmanned aerial vehicles into the National Airspace, ensuring safety of the airspace including the vehicles, people, and property on the ground is of utmost importance. One of the safety-critical factors for unmanned aviation flights is the risk of deviating from a planned trajectory resulting in a variety of hazards, including potential loss of separation between vehicle and obstacles or unexpected battery energy consumption. Off-nominal conditions introduced by component failures, degraded controllability and environmental disturbances such as wind gusts can lead to unacceptable unexpected deviations from the flight trajectory. It is essential to accurately model such effects on the flight trajectory while computing safety thresholds such as minimum separation from surrounding obstacles, available battery resource to complete the mission or determining delay in the expected time of arrival of flights. In this paper, a tool is presented based on Gaussian Process Regression for wind representation over a pre-defined trajectory for fast, yet approximated, in-time evaluation of possible trajectory deviations caused by wind gusts. The deviation in the planned trajectory caused by wind is further simulated utilizing a 6 degrees-of-freedom (DOF) UAV trajectory simulator comprising of a rotorcraft lumped-mass model with LQRI controller. Both steady-state wind and wind gust effects are investigated. The probability of collision with obstacle is computed and demonstrated on real flight data from experimental flights of an octocopter at NASA Langley Research Center in the presence of simulated obstacles and wind conditions. Effect of varying wind conditions and varying UAV airspeed is further demonstrated on experimental flights in the presence of wind measured by ground based weather service stations. The proposed approach would eventually benefit timely mitigation of current and future safety-critical events in autonomous systems by enabling risk-informed decision making.

Portia Banerjee↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2024

Energy Exascale Computational Fluid Dynamics Simulations With the Spectral Element Method

Development and application of the open-source GPU-based fluid-thermal simulation code, NekRS, are described. Time advancement is based on an efficient kth-order accurate timesplit formulation coupled with scalable iterative solvers. Spatial discretization is based on the high-order spectral element method (SEM), which affords the use of fast, low-memory, matrix-free operator evaluation. Further, recent developments include support for nonconforming meshes using overset grids and for GPU-based Lagrangian particle tracking. Results of large-eddy simulations of atmospheric boundary layers for wind-energy applications as well as extensive nuclear energy applications are presented.

42 ENGINEERING↗
NASA NTRS2011

TPSAS-NF1676L-12387-DND

Fast-time wake prediction models are semi-empirical, and are designed to predict aircraft wake vortex information behind a generating aircraft. They are developed using theoretical concepts and often calibrated with data obtained from field measurements, laboratory studies, or numerical simulations. They rely upon inputs representing the generating aircraft (such as weight, airspeed, wingspan, lateral position, and altitude) and atmospheric environment (such as ambient winds, turbulence, and temperature lapse rate). Other input may include the type of aircraft (e.g., B-747- 400). Output from these models may include time histories of the wake vortex strength, and position, as well as a prediction of the uncertainties in strength and position. Fast-time models should provide solutions in a fraction of a second, and hence is an enabling technology able to perform in a real-time operational system, such as dynamic spacing systems envisioned by the Aircraft Vortex Spacing System or other Super Density Operations concepts. Other applications of fast-time models include 1) guidance for wake separation standards for new aircraft, 2) optimization of separation standards for existing aircraft, 3) guidance for setting vertical separation standards during cruise, 4) accident reconstructions, and 5) testing of new concepts for safe increases in airport capacity for both single and multiple runway configurations. Fast-time models are based on a number of assumptions which limit their accuracy and range of application. For example, they do not predict the meandering vortex paths that can result from the wake interaction with atmospheric turbulence. However, these models should account for the mean vortex position and may predict the level of uncertainty. This presentation describes a study in which comparative evaluations of fast-time wake transport and decay models were conducted using data sets from previous field programs taken at MEM (1995), DFW (1997, 1999, 2000), and Denver (2005, 2006) airports. These data sets contain, meteorological data, lidar measurements of wake vortices, and the aircraft parameters associated with these vortices.

Fred H Proctor↗
DOE OSTI2024

Stellar wind impact on early atmospheres around unmagnetized Earth-like planets

ABSTRACT Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability.

Astronomy & Astrophysics↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

SNIa Cosmology Analysis Results from Simulated LSST Images: From Difference Imaging to Constraints on Dark Energy

Abstract The Vera Rubin Observatory Legacy Survey of Space and Time (LSST) is expected to process ∼10 6 transient detections per night. For precision measurements of cosmological parameters and rates, it is critical to understand the detection efficiency, magnitude limits, artifact contamination levels, and biases in the selection and photometry. Here we rigorously test the LSST Difference Image Analysis (DIA) pipeline using simulated images from the Rubin Observatory LSST Dark Energy Science Collaboration Data Challenge (DC2) simulation for the Wide-Fast-Deep survey area. DC2 is the first large-scale (300 deg 2 ) image simulation of a transient survey that includes realistic cadence, variable observing conditions, and CCD image artifacts. We analyze ∼15 deg 2 of DC2 over a 5 yr time span in which artificial point sources from Type Ia supernova (SNIa) light curves have been overlaid onto the images. The magnitude limits per filter are u = 23.66 mag, g = 24.69 mag, r = 24.06 mag, i = 23.45 mag, z = 22.54 mag, and y = 21.62 mag. The artifact contamination levels are ∼90% of all detections, corresponding to ∼1000 artifacts deg –2 in g band, and falling to 300 deg –2 in y band. The photometry has biases <1% for magnitudes 19.5 < m < 23. Our DIA performance on simulated images is similar to that of the Dark Energy Survey difference-imaging pipeline on real images. We also characterize DC2 image properties to produce catalog-level simulations needed for distance bias corrections. We find good agreement between DC2 data and simulations for distributions of signal-to-noise ratio, redshift, and fitted light-curve properties. Applying a realistic SNIa cosmology analysis for redshifts z < 1, we recover the input cosmology parameters to within statistical uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Development of Predictive Model for Accurate Rupture Time from Multi-Axial Creep in Alloy 709 with Physics-Based Simulations

A physics-based model is developed to predict multiaxial creep behavior in Alloy 709 (A709), an advanced austenitic stainless steel intended for high-temperature applications such as Sodium Fast Reactors (SFRs). Compared to conventional stainless steels like 316H, A709 offers superior high-temperature performance; however, comprehensive data on its multiaxial creep response remain limited. To address this gap, a crystal plasticity finite element (CPFE) framework is used to simulate the deformation and failure mechanisms of A709 under multiaxial loading conditions. The model incorporates an extended Hu-Cocks dislocation creep formulation that accounts for precipitation effects, along with the Sham–Needleman model to capture grain boundary cavitation-driven failure. These advanced constitutive models enable a detailed understanding of the interplay between microstructural evolution and macroscopic creep response. Furthermore, the study evaluates the predictive accuracy of various effective stress measures in estimating creep rupture life, leveraging simulated multiaxial creep data. The findings provide critical insights into the applicability of different stress measures for engineering design and life prediction of A709 components operating under complex loading conditions. This work contributes to improving the reliability of high-temperature structural components by advancing predictive modeling capabilities for advanced austenitic steels.

Alloy 709↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Fine structure of microwave spike bursts and associated cross-field energy transport

The characteristics of the maser emission from a driven system where energetic electrons continue to flow through the source region is investigated using electronic particle simulations. It is shown that, under appropriate conditions, the maser can efficiently radiate a significant portion of the energy of the fast electrons in a very short time. The radiation is emitted in pulses even though the flow of electrons through the system is at a constant rate. The mission of these pulses is proposed as the source of the fine structure. Under other conditions the dominant maser emission changes from fundamental x-mode to either fundamental z-mode or to electrostatic upper hybrid or Bernstein modes. The bulk of the emission from the maser instability cannot propagate across field lines in this regime, and hence strong local plasma heating is expected, with little energy transport across the magnetic field lines.

Winglee, R. M.↗
DOE OSTI2021

Examination of Semi-Analytical Solution Methods in the Coarse Operator of Parareal Algorithm for Power System Simulation

With continuing advances in high-performance parallel computing platforms, parallel algorithms have become powerful tools for development of faster than real-time power system dynamic simulations. In particular, it has been demonstrated in recent years that parallel-in-time (Parareal) algorithms have the potential to achieve such an ambitious goal. Here, the selection of a fast and reasonably accurate coarse operator of the Parareal algorithm is crucial for its effective utilization and performance. This paper examines semi-analytical solution (SAS) methods as the coarse operators of the Parareal algorithm and explores performance of the SAS methods to the standard numerical time integration methods. Two promising time-power series-based SAS methods were considered; Adomian decomposition method and Homotopy analysis method with a windowing approach for improving the convergence. Numerical performance case studies on 10-generator 39-bus system and 327-generator 2383-bus system were performed for these coarse operators over different disturbances, evaluating the number of Parareal iterations, computational time, and stability of convergence. All the coarse operators tested with different scenarios have converged to the same corresponding true solution (if they are convergent) and the SAS methods provide comparable computational speed, while having more stable convergence to the true solution in many cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Chatter detection in simulated machining data: a simple refined approach to vibration data

Vibration monitoring is a critical aspect of assessing the health and performance of machinery and industrial processes. This study explores the application of machine learning techniques, specifically the Random Forest (RF) classification model, to predict and classify chatter—a detrimental self-excited vibration phenomenon—during machining operations. While sophisticated methods have been employed to address chatter, this research investigates the efficacy of a novel approach to an RF model. The study leverages simulated vibration data, bypassing resource-intensive real-world data collection, to develop a versatile chatter detection model applicable across diverse machining configurations. The feature extraction process combines time-series features and Fast Fourier Transform (FFT) data features, streamlining the model while addressing challenges posed by feature selection. By focusing on the RF model’s simplicity and efficiency, this research advances chatter detection techniques, offering a practical tool with improved generalizability, computational efficiency, and ease of interpretation. The study demonstrates that innovation can reside in simplicity, opening avenues for wider applicability and accelerated progress in the machining industry.

42 ENGINEERING↗
DOE OSTI2023

Co-Simulation of Electric Power Distribution Systems and Buildings including Ultra-Fast HVAC Models and Optimal DER Control

Smart homes and virtual power plant (VPP) controls are growing fields of research with potential for improved electric power grid operation. A novel testbed for the co-simulation of electric power distribution systems and distributed energy resources (DERs) is employed to evaluate VPP scenarios and propose an optimization procedure. DERs of specific interest include behind-the-meter (BTM) solar photovoltaic (PV) systems as well as heating, ventilation, and air-conditioning (HVAC) systems. The simulation of HVAC systems is enabled by a machine learning procedure that produces ultra-fast models for electric power and indoor temperature of associated buildings that are up to 133 times faster than typical white-box implementations. Hundreds of these models, each with different properties, are randomly populated into a modified IEEE 123-bus test system to represent a typical U.S. community. Advanced VPP controls are developed based on the Consumer Technology Association (CTA) 2045 standard to leverage HVAC systems as generalized energy storage (GES) such that BTM solar PV is better utilized locally and occurrences of distribution system power peaks are reduced, while also maintaining occupant thermal comfort. An optimization is performed to determine the best control settings for targeted peak power and total daily energy increase minimization with example peak load reductions of 25+%.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗