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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A Modeling Approach for Time-Dependent Geometry Applied to Transient Heat Transfer of Aluminum Electrolysis Cells

The thermal balance of aluminum electrolysis cells (AEC) have to be rigorously controlled in order to improve the efficiency and sustainability of this industrial process. A new modeling strategy is developed to consider the displacements of solid bodies and moving boundaries in finite element models. The transient thermal-electric modeling of the AEC demonstrates the effect of an increase in operating voltage on both the anode cover and the side ledge. With higher heat generation, the anode cover deteriorates and the side ledge thickness decreases. Since the anode cover is characterized by irreversible transformations, the top heat dissipation remains higher even when the operating voltage comes back to its typical value. For the first time, the transient temperature and electric fields throughout the anode life are simulated and validated by industrial measurements. The modeling predictions have been validated from instrumented anodes and manual measurements, all performed on operating AEC.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Conjugate heat transfer↗
DOE OSTI2021

Development of the reactor core isolation cooling system model for the extended station black-out accident

In order to accurately predict the behavior of the Reactor Core Isolation Cooling System (RCIC) in the context of the FLEX strategies under the extended station black-out accident (SBO), we present a mechanistic model of the RCIC in this work. The control volume formulation of angular momentum conservation is employed to model the Terry turbine in the RCIC system. The two-phase nozzle flow model developed in the previous work is utilized to simulate the convergent-divergent nozzles in the Terry turbine. A nozzle flow table is generated to cover the operation range of the RCIC during beyond design basis accidents (BDBAs). The mechanistic RCIC turbine-pump model and the nozzle flow table are incorporated in a TRACE BWR model. Simulations of short-term and long-term SBO transients are performed to demonstrate dynamic RCIC responses. A system analysis code with the proposed mechanistic model of the RCIC system will have a better prediction of system safety response during SBO accidents. With the addition of such more realistic models, the system code will be capable of providing expanded understanding of safety margin, which can help help to form the technical basis for the successful implementation of the FLEX and SAMG measures under Extended Loss of AC Power (ELAP) conditions for BWR or PWR reactor designs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Nowcasting, predictive control, and feedback control for temperature regulation in a novel hybrid solar-electric reactor for continuous solar-thermal chemical processing

A novel 15 kW reactor for the hybridization of concentrated solar and conventional electric heat was fabricated for renewable and continuous chemical processing at temperatures up to 1700 °C. Solar-electric controllers based on feedback or predictive linear models were used to regulate the device at 925 °C for the production of syngas via the gasification of carbon. The system was challenged with cloud transients programmed on a 45 kW high-flux solar simulator that approximated weather observed in the San Luis Valley (Colorado, USA). In experiment it was found that model predictive control with a 1 min ahead nowcast of incipient clouding best regulated solar-electric reactor temperatures, potentially averting thermal fatigue. Upon clouding, model predictive control with a nowcast yielded temperature disturbances of ±10 °C, whereas feedback control alone featured ±25 °C excursions. Overall, the performance of model predictive control with a nowcast was 75% better than feedback as indicated by the integral squared error of temperature excursion residuals. The importance of nowcast accuracy was explored, and although the model-based controller was robust to forecast amplitude inaccuracy, temporal forecast inaccuracy wholly negated the benefits of predictive control.

14 SOLAR ENERGY↗
DOE OSTI2025

Decoupled few-femtosecond phase transitions in vanadium dioxide

Abstract The nature of the insulator-to-metal phase transition in vanadium dioxide (VO 2 ) is one of the longest-standing problems in condensed-matter physics. Ultrafast spectroscopy has long promised to determine whether the transition is primarily driven by the electronic or structural degree of freedom, but measurements to date have been stymied by their sensitivity to only one of these components and/or their limited temporal resolution. Here we use ultra-broadband few-femtosecond pump-probe spectroscopy to resolve the electronic and structural phase transitions in VO 2 at their fundamental time scales. Our experiments show that the system transforms into a bad-metallic phase within 10 fs after photoexcitation, but requires another 100 fs to complete the transition, during which we observe electronic oscillations and a partial re-opening of the bandgap, signalling a transient semi-metallic state. Comparisons with tensor-network simulations and density-functional theory calculations show these features result from an unexpectedly fast structural transition, in which the vanadium dimers separate and untwist with two different timescales. Our results resolve the structural and electronic nature of the light-induced phase transition in VO 2 and establish ultra-broadband few-femtosecond spectroscopy as a powerful tool for studying quantum materials out of equilibrium.

Science & Technology - Other Topics↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Leveraging Optimal Sparse Sensor Placement to Aggregate a Network of Digital Twins for Nuclear Subsystems

Nuclear power plants (NPPs) require continuous monitoring of various systems, structures, and components to ensure safe and efficient operations. The critical safety testing of new fuel compositions and the analysis of the effects of power transients on core temperatures can be achieved through modeling and simulations. They capture the dynamics of the physical phenomenon associated with failure modes and facilitate the creation of digital twins (DTs). Accurate reconstruction of fields of interest (e.g., temperature, pressure, velocity) from sensor measurements is crucial to establish a two-way communication between physical experiments and models. Sensor placement is highly constrained in most nuclear subsystems due to challenging operating conditions and inherent spatial limitations. This study develops optimized data-driven sensor placements for full-field reconstruction within reactor and steam generator subsystems of NPPs. Optimized constrained sensors reconstruct field of interest within a tri-structural isotropic (TRISO) fuel irradiation experiment, a lumped parameter model of a nuclear fuel test rod and a steam generator. The optimization procedure leverages reduced-order models of flow physics to provide a highly accurate full-field reconstruction of responses of interest, noise-induced uncertainty quantification and physically feasible sensor locations. Accurate sensor-based reconstructions establish a foundation for the digital twinning of subsystems, culminating in a comprehensive DT aggregate of an NPP.

42 ENGINEERING↗
DOE OSTI2018

Artifacts in time-resolved Kelvin probe force microscopy

Kelvin probe force microscopy (KPFM) has been used for the characterization of metals, insulators, and semiconducting materials on the nanometer scale. Especially in semiconductors, the charge dynamics are of high interest. Recently, several techniques for time-resolved measurements with time resolution down to picoseconds have been developed, many times using a modulated excitation signal, e.g. light modulation or bias modulation that induces changes in the charge carrier distribution. For fast modulation frequencies, the KPFM controller measures an average surface potential, which contains information about the involved charge carrier dynamics. Here, we show that such measurements are prone to artifacts due to frequency mixing, by performing numerical dynamics simulations of the cantilever oscillation in KPFM subjected to a bias-modulated signal. For square bias pulses, the resulting time-dependent electrostatic forces are very complex and result in intricate mixing of frequencies that may, in some cases, have a component at the detection frequency, leading to falsified KPFM measurements. Additionally, we performed fast Fourier transform (FFT) analyses that match the results of the numerical dynamics simulations. Small differences are observed that can be attributed to transients and higher-order Fourier components, as a consequence of the intricate nature of the cantilever driving forces. These results are corroborated by experimental measurements on a model system. In the experimental case, additional artifacts are observed due to constructive or destructive interference of the bias modulation with the cantilever oscillation. Also, in the case of light modulation, we demonstrate artifacts due to unwanted illumination of the photodetector of the beam deflection detection system. Lastly, guidelines for avoiding such artifacts are given.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Electromagnetic Signatures from the Tidal Tail of a Black Hole—Neutron Star Merger

Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle because of two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of ~3 with orientation while the peak in the optical bands varies by ~3 magnitudes. The spectrum is blueshifted at viewing angles along the bulk motion, which increases the V-band peak magnitude to ~-14 despite the lanthanide-rich composition.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1982

A variable structure approach to robust control of VTOL aircraft

This paper examines the application of variable structure control theory to the design of a flight control system for the AV-8A Harrier in a hover mode. The objective in variable structure design is to confine the motion to a subspace of the total state space. The motion in this subspace is insensitive to system parameter variations and external disturbances that lie in the range space of the control. A switching type of control law results from the design procedure. The control system was designed to track a vector velocity command defined in the body frame. For comparison purposes, a proportional controller was designed using optimal linear regulator theory. Both control designs were first evaluated for transient response performance using a linearized model, then a nonlinear simulation study of a hovering approach to landing was conducted. Wind turbulence was modeled using a 1052 destroyer class air wake model.

Calise, A. J.↗
NASA NTRS1984

Nonlinear transient finite element analysis of rotor-bearing-stator systems

This paper extends the finite element scheme to handle the highly nonlinear interfacial fields generated in the fluid filled annulli of squeeze film and journal bearings so as to model the transient response of rotor-bearing-stator systems. Since such simulations are highly nonlinear, direct numerical integration schemes are employed to generate the overall response. In this context, the paper gives consideration to such items as (1) numerical efficiency/stability, (2) comparison of implicit and explicit schemes, (3) determines extent of response nonlinearity as well as (4) extensively benchmarks the overall concept/methodologies.

Padovan, J.↗
NASA NTRS1984

A variable structure approach to robust control of VTOL aircraft

This paper examines the application of variable structure control theory to the design of a flight control system for the AV-8A Harrier in a hover mode. The objective in variable structure design is to confine the state trajectories to a subspace of the total state space. The motion in this subspace is insensitive to system parameter variations and external disturbances that lie in the range space of the control. A switching type of control law results from the design procedure. The control system was designed to track a vector-valued velocity command. For comparison, a proportional controller was designed using optimal linear regulator theory. Both controllers were evaluated for their transient response performance using a linear model; then a nonlinear simulation study of a hovering approach to landing was conducted. The variable structure controller outperformed its linear counterpart in the presence of wind disturbances and plant parameter uncertainties afforded by the simulation.

Calise, A. J.↗
NASA NTRS1985

Flame flashback in a premixed dump combustor

A stainless steel, rectangular center-dump, premixed-prevaporized combustor was used at the NASA Lewis Research Center to study flashback in controlled conditions simulating gas turbine combustor operations under steady-state and transient conditions. The combustor allowed visual access to the combustion area, and was operated at inlet air temperatures from 600 to 850 K, premixer wall temperatures from 450 to 900 K, and average premixer velocity from 40 to 80 ft/s. A parametric study revealed a slight decrease in the flashback equivalence ratio as the inlet air temperatures increased. It also indicated that the premixer wall temperature and premixer velocity are not governing parameters of flashback. It is suggested that flashback occurs through the premixer wall boundary layer flow reversal caused by combustion instability.

Proctor, M. P.↗
NASA NTRS1988

Development of a two-dimensional skin friction balance nulling circuit using multivariable control theory

Measurement of planar skin friction forces in aerodynamic testing currently requires installation of two perpendicularly mounted, single-axis balances; consequently, force components must be sensed at two distinct locations. A two-axis instrument developed at the Langley Research Center to overcome this disadvantage allows measurement of a two-dimensional force at one location. This paper describes a feedback-controlled nulling circuit developed for the NASA two-axis balance which, without external compensation, is inherently unstable because of its low friction mechanical design. Linear multivariable control theory is applied to an experimentally validated mathematical model of the balance to synthesize a state-variable feedback control law. Pole placement techniques and computer simulation studies are employed to select eigenvalues which provide ideal transient response with decoupled sensing dynamics.

Tripp, John S.↗
NASA NTRS1990

Heat-Flux-Measuring Facility

Apparatus simulates conditions in turbine engines. Automated facility generates and measures transient and steady-state heat fluxes at flux densities from 0.3 to 6 MW/m(Sup2) and temperatures from 100 to 1,200 K. Positioning arm holds heat-flux gauge at focal point of arc lamp. Arm previously chilled gauge in liquid nitrogen in Dewar flask. Cooling water flows through lamp to heat exchanger. Used to develop heat-flux gauges for turbine blades and to test materials for durability under rapidly changing temperatures.

Liebert, Curt H.↗
NASA NTRS1990

Influence of residual gravity on crystal growth processes

The effects of steady and time-dependent residual accelerations on representative melt and solution crystal growth systems have been investigated using numerical models. For the melt growth system, it was found that the lateral nonuniformity in crystal composition, depending on the particular operating conditions, can be very sensitive to the orientation of the steady component of the residual gravity vector down to 10 to the -6th g. Oscillatory residual accelerations (g-jitter) were found to lead to transients in the dopant distribution. For the solution growth, longtime simulations showed that at 10 to the -5th g, and the cooling rates used in Spacelab 3 experiments, diffusive transport dominates only in the beginning. Convective transport becomes increasingly significant after 3 to 5 hours, depending on the orientation of the crystal in the residual gravity field.

Rosenberger, F.↗
NASA NTRS1991

A computer program to perform flow and thermal analysis during pressurization of the Space Shuttle Solid Rocket Motor field joint

A computational technique for prediction of the flow and thermal environment in the SRM field joint cavities is described. The SRM field joint hardware was tested with a defect in the insulation, and due to this defect, the O-ring gland cavities are pressurized during the early part of the ignition. A computer model is developed to predict the thermal environment and flow through the simulated flaw, during the pressurization of the field joint. The transient mass, momentum, and energy conservation equations in the flow passage together with the thermodynamic equation of state are solved by a fully implicit iterative numerical procedure.

Clayton, J. Louie↗