Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient 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 1,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

From core collapse to superluminous: the rates of massive stellar explosions from the Palomar Transient Factory

ABSTRACT We present measurements of the local core-collapse supernova (CCSN) rate using SN discoveries from the Palomar Transient Factory (PTF). We use a Monte Carlo simulation of hundreds of millions of SN light-curve realizations coupled with the detailed PTF survey detection efficiencies to forward model the SN rates in PTF. Using a sample of 86 CCSNe, including 26 stripped-envelope SNe (SESNe), we show that the overall CCSN volumetric rate is $r^\mathrm{CC}_v=9.10_{-1.27}^{+1.56}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$ at 〈z〉 = 0.028, and the SESN volumetric rate is $r^\mathrm{SE}_v=2.41_{-0.64}^{+0.81}\times 10^{-5}\, \text{SNe yr}^{-1}\, \text{Mpc}^{-3}\, h_{70}^{3}$. We further measure a volumetric rate for hydrogen-free superluminous SNe (SLSNe-I) using eight events at z ≤ 0.2 of $r^\mathrm{SLSN-I}_v=35_{-13}^{+25}\, \text{SNe yr}^{-1}\text{Gpc}^{-3}\, h_{70}^{3}$, which represents the most precise SLSN-I rate measurement to date. Using a simple cosmic star formation history to adjust these volumetric rate measurements to the same redshift, we measure a local ratio of SLSN-I to SESN of ${\sim}1/810^{+1500}_{-94}$, and of SLSN-I to all CCSN types of ${\sim}1/3500^{+2800}_{-720}$. However, using host galaxy stellar mass as a proxy for metallicity, we also show that this ratio is strongly metallicity dependent: in low-mass (logM* < 9.5 M⊙) galaxies, which are the only environments that host SLSN-I in our sample, we measure an SLSN-I to SESN fraction of $1/300^{+380}_{-170}$ and $1/1700^{+1800}_{-720}$ for all CCSN. We further investigate the SN rates a function of host galaxy stellar mass, and show that the specific rates of all CCSNe decrease with increasing stellar mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Black-Starting Microgrids Using Inverter-Based Distributed Energy Resources

Recent trends in power systems show increasing deployment of inverter-based distributed energy resources (DERs) at the distribution level. Efforts have also been made to use inverter-based DERs for system restoration of critical sections of distribution feeders. In the past few years, there has been active consideration of using local grid-forming (GFM) DERs, which can form their own voltage and frequency, to black-start microgrids instead of waiting for the bulk system for supply restoration. However, using current-limited inverter-based DERs for black-start operations with high inrush currents can pose some technical challenges. This paper presents detailed electromagnetic transient studies on a digital real-time simulator to assess the operational feasibility of using local inverter-based GFM DERs for microgrid black-start operations as a bottom-up approach instead of a top-down approach where the substation is used to restore critical sections of the feeder. The paper presents key modeling details for distribution system components for studying GFM DER-based microgrid black-start operations. The paper presents simulation results comparing different black-start approaches using GFM DERs for two test systems: i) a simple test microgrid and ii) a real-world utility microgrid currently being designed and deployed by Holy Cross Energy in Colorado.

black start↗
DOE OSTI2020

Grid-Following Inverters and Synchronous Condensers: A Grid-Forming Pair?

As power systems across the globe continue to achieve higher instantaneous penetrations of power electronic converter-interfaced renewable sources, the stability of these power systems is challenged because of the removal of synchronous machines, from which stability is traditionally derived. Although technological solutions to these stability challenges are on the horizon, such as the use of grid-forming inverters, they are not yet widely applied to larger power systems, which presents operational challenges for power systems achieving these high instantaneous penetrations today. A potential interim solution using existing technologies is to pair synchronous condensers with grid-following inverters, which might prolong the stability of an operating power system while synchronous generators are turned off during periods of high renewable energy availability. This work examines the transient stability of such a solution using PSCAD simulations of a two-bus system coupled to a varied-length transmission line, with a synchronous condenser at one bus and a grid-following inverter with grid-support functionality at the other. The system is exposed to load step, balanced fault, and unbalanced fault perturbations. It was found that this simple system returned to a steady state for transmission line lengths up to 125 km after 10% load steps, or a variety of fault types.

grid-following inverters↗
DOE OSTI2022

Network Uncertainty Quantification for Analysis of Multi-Component Systems

To impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

42 ENGINEERING↗
DOE OSTI2021

FIERRO V.X

FIERRO is a modern C++ code intended to simulate quasi-static solid mechanics problems and transient, compressible material dynamic problems with Lagrangian methods, which have meshes with constant mass elements that move with the material. FIERRO is designed to aid a) material model research that has historically been done using commercial implicit and explicit finite element codes, b) numerical methods research, and c) computer science research. FIERRO supports user developed material models that adhere to several industry standard formats by using a C++ to Fortran interface to couple the model to the numerical solvers. FIERRO is built on the ELEMENTS library that supports a diverse suite of element types, including high-order elements, and quadrature rules. The mesh class within the ELEMENTS library is designed for efficient calculations on unstructured meshes and to minimize memory usage. FIERRO is designed to readily accommodate a range of numerical methods including continuous finite element, finite volume, and discontinuous Galerkin methods. FIERRO is designed to support explicit and implicit time integration methods. FIERRO only supports a single material in an element. No physical data exists within the code.

Morgan, Nathaniel↗
DOE OSTI2021

Implementation of Monte Carlo Transport in the General Relativistic SpEC Code

Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the last is a subdominant source of error at the accuracy reached by current simulations and for the interactions currently included in our code. We also provide tests of the most important features of this code.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Opportunities in AI/ML for the Rubin LSST Dark Energy Science Collaboration

The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will produce unprecedented volumes of heterogeneous astronomical data (images, catalogs, and alerts) that challenge traditional analysis pipelines. The LSST Dark Energy Science Collaboration (DESC) aims to derive robust constraints on dark energy and dark matter from these data, requiring methods that are statistically powerful, scalable, and operationally reliable. Artificial intelligence and machine learning (AI/ML) are already embedded across DESC science workflows, from photometric redshifts and transient classification to weak lensing inference and cosmological simulations. Yet their utility for precision cosmology hinges on trustworthy uncertainty quantification, robustness to covariate shift and model misspecification, and reproducible integration within scientific pipelines. This white paper surveys the current landscape of AI/ML across DESC's primary cosmological probes and cross-cutting analyses, revealing that the same core methodologies and fundamental challenges recur across disparate science cases. Since progress on these cross-cutting challenges would benefit multiple probes simultaneously, we identify key methodological research priorities, including Bayesian inference at scale, physics-informed methods, validation frameworks, and active learning for discovery. With an eye on emerging techniques, we also explore the potential of the latest foundation model methodologies and LLM-driven agentic AI systems to reshape DESC workflows, provided their deployment is coupled with rigorous evaluation and governance. Finally, we discuss critical software, computing, data infrastructure, and human capital requirements for the successful deployment of these new methodologies, and consider associated risks and opportunities for broader coordination with external actors.

Aubourg, Eric [APC, Paris] (ORCID:000000025592023X↗
NASA NTRS1973

Capabilities of electrodynamic shakers when used for mechanical shock testing

The results of a research task to investigate the capabilities of electrodynamic vibrators (shakers) to perform mechanical shock tests are presented. The simulation method employed was that of developing a transient whose shock response spectrum matched the desired shock response spectrum. Areas investigated included the maximum amplitude capabilities of the shaker systems, the ability to control the shape of the resultant shock response spectrum, the response levels induced at frequencies outside the controlled bandwidth, and the nonlinearities in structural response induced by a change in test level.

Keegan, W. B.↗
NASA NTRS1973

Numerical simulation of small-scale thermal convection in the atmosphere

A Boussinesq system is integrated numerically in three dimensions and time in a study of nonhydrostatic convection in the atmosphere. Simulation of cloud convection is achieved by the inclusion of parametrized effects of latent heat and small-scale turbulence. The results are compared with the cell structure observed in Rayleigh-Benard laboratory conversion experiments in air. At a Rayleigh number of 4000, the numerical model adequately simulates the experimentally observed evolution, including some prominent transients of a flow from a randomly perturbed initial conductive state into the final state of steady large-amplitude two-dimensional rolls. At Rayleigh number 9000, the model reproduces the experimentally observed unsteady equilibrium of vertically coherent oscillatory waves superimposed on rolls.

Somerville, R. C. J.↗
NASA NTRS1976

Experimental study of transient dynamics of a flexible rotor

The results of an experimental program to investigate the transient response of a flexible rotor are presented. The program consisted of a series of tests conducted on a rig with a rotor designed to operate above its first bending critical speed. The purpose of the tests was to obtain experimental data on the transient behavior of a flexible rotor under conditions simulating those which might occur in a jet aircraft engine. The scope of the program included tests to measure the response of both balanced and unbalanced rotors during steady-state operation, acceleration, deceleration, and simulated blade loss.

Hibner, D. H.↗
NASA NTRS1977

Computer aided nonlinear electrical networks analysis

Techniques used in simulating an electrical circuit with nonlinear elements for use in computer-aided circuit analysis programs are described. Elements of the circuit include capacitors, resistors, inductors, transistors, diodes, and voltage and current sources (constant or time varying). Simulation features are discussed for dc, ac, and/or transient circuit analysis. Calculations are based on the model approach of formulating the circuit equations. A particular solution of transient analysis for nonlinear storage elements is described.

Slapnicar, P.↗
NASA NTRS1979

Two-dimensional, time-dependent MHD description of interplanetary disturbances Simulation of high speed solar wind interactions

A two-dimensional, time-dependent nonlinear MHD model is used to simulate the general features of multiply-interacting, transient solar wind flows. The major advantage of the two-dimensional model over the one-dimensional one is the proper consideration of azimuthal gradients. In addition, the new model can provide contour maps for disturbed plasma properties in the equatorial plane, thereby enabling one to perform essential tests of the physical assumptions through direct comparison with space probe data.

Wu, S. T.↗
NASA NTRS1980

Propulsion system mathematical model for a lift/cruise fan V/STOL aircraft

A propulsion system mathematical model is documented that allows calculation of internal engine parameters during transient operation. A non-realtime digital computer simulation of the model is presented. It is used to investigate thrust response and modulation requirements as well as the impact of duty cycle on engine life and design criteria. Comparison of simulation results with steady-state cycle deck calculations showed good agreement. The model was developed for a specific 3-fan subsonic V/STOL aircraft application, but it can be adapted for use with any similar lift/cruise V/STOL configuration.

Cole, G. L.↗
NASA NTRS1985

SSME fuel preburner phase 1/phase 2+ comparison

An analytical comparison of the predicted combustion flowfields for the Phase 1 and Phase 2+ fuel preburner configurations using an advanced CFD combustion code is currently underway at Rocketdyne. The Phase 2+ injector has a modified element design. The redesigned element provides for improved atomization through a shift in the injector element/face bleed flow split while maintaining an acceptable pressure drop. Cold-flow tests of the new design have been encouraging. As a complement to those tests, the Advanced Rocket Injector Combustor Code (ARICC) was selected to conduct a detailed combustion flowfield simulation. ARICC models two-dimensional (or axisymmetric), transient, turbulent, two-phase mixing and combustion flowfields. It is derived from a member of the Los Alamos ICEd-ALE family of codes. Unique features of ARICC include explicit representation of the coaxial LOX jet and hydrogen gas flows, distributed atomization processes, droplet breakup, and supercritical vaporation processes. Predictions of temperature and OH concentration profiles for the Phase 1 and Phase 2+ injector flowfields indicate a trend toward more uniform temperature distributions and shorter flame lengths with the Phase 2+ design.

Jensen, R. J.↗
NASA NTRS1987

Nonlinear dynamics of attractive magnetic bearings

The nonlinear dynamics of a ferromagnetic shaft suspended by the force of attraction of 1, 2, or 4 independent electromagnets is presented. Each model includes a state variable feedback controller which has been designed using the pole placement method. The constitutive relationships for the magnets are derived analytically from magnetic circuit theory, and the effects of induced eddy currents due to the rotation of the journal are included using Maxwell's field relations. A rotor suspended by four electro-magnets with closed loop feedback is shown to have nine equilibrium points within the bearing clearance space. As the rotor spin speed increases, the system is shown to pass through a Hopf bifurcation (a flutter instability). Using center manifold theory, this bifurcation can be shown to be of the subcritical type, indicating an unstable limit cycle below the critical speed. The bearing is very sensitive to initial conditions, and the equilibrium position is easily upset by transient excitation. The results are confirmed by numerical simulation.

Hebbale, K. V.↗