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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Liquid lithium as divertor material to mitigate severe damage of nearby components during plasma transients

The successful operation of thermonuclear fusion reactors such as ITER, DEMO, and future commercial plants is mainly determined by the optimum choice of materials for various components. The objective of this work is to accurately and comprehensively simulate the entire device in 3D to predict pros and cons of various materials, e.g., liquid lithium in comparison to tungsten and carbon to predict future ITER-like and DEMO divertor performances. We used our comprehensive HEIGHTS simulation package to investigate ITER-like components response during transient events in exact 3D geometry. Starting from the lost hot core plasma particles through SOL, deposition on the divertor surface, and the generation of secondary plasma of divertor materials. Our simulations predicted significant reduction in the heat loading and damage to the divertor nearby and internal components in the case when lithium is used on the divertor plates. While if tungsten or carbon are used on the divertor plate, significant melting areas and vaporization spots can occur (less for carbon) on the reflector, dome, and stainless steel tubes, and even parts of the first walls can melt due to the high radiation power of the secondary divertor plasma. Lithium photon radiation deposition into the divertor and nearby surfaces was decreased by two orders of magnitude compared to tungsten and by one order of magnitude compared to carbon. This analysis showed that using liquid lithium for ITER-like surfaces and future DEMO can lead to significant enhancement in components lifetime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1978

Transient dynamics of a flexible rotor with squeeze film dampers

A series of simulated blade loss tests are reported on a test rotor designed to operate above its second bending critical speed. A series of analyses were performed which predicted the transient behavior of the test rig for each of the blade loss tests. The scope of the program included the investigation of transient rotor dynamics of a flexible rotor system, similar to modern flexible jet engine rotors, both with and without squeeze film dampers. The results substantiate the effectiveness of squeeze film dampers and document the ability of available analytical methods to predict their effectiveness and behavior.

Buono, D. F.↗
DOE OSTI2025

An airfoil-based synthetic actuator disk model for wind turbine aerodynamic and structural analysis

Here, this study introduces an airfoil-based refinement technique to enhance the Actuator Disk Model (ADM) for improved wind turbine aerodynamic load prediction and structural simulation in conjunction with Large Eddy Simulations of the wind flow. While ADM offers higher computational efficiency than the more detailed but resource-intensive Actuator Line Model (ALM), it traditionally lacks the resolution needed to capture the localized blade forces accurately. To address this limitation, we introduce a refinement technique that uses airfoil-specific data and employs interpolation-based grid point refinement, achieving ALM-comparable accuracy while preserving ADM's efficiency. Unlike conventional ADM that provides only rotor-disk averaged forces, our synthetic method tracks transient aerodynamic load variations over multiple blade revolutions, allowing us to calculate the distributions of maximum and minimum loads during typical cycles. Applied to the NREL 5 MW reference turbine, our enhanced ADM accurately predicts key aerodynamic parameters (angle of attack, axial velocity, lift, drag, axial and tangential forces along the blades) as well as structural responses (blade tip deflection, maximum stress, and stress concentration). Our results show that the tip deflection ranges from 2.33m (3.69 % of blade length) to 4.28m (6.79 %), with maximum stress concentration occurring near the blade root. This research demonstrates that a refined synthetic ADM approach can serve as a computationally efficient alternative for both aerodynamic analysis and structural simulation of wind turbine blades subjected to realistic wind fields.

17 WIND ENERGY↗
NASA NTRS1986

Linear simulation of the stationary eddies in a general circulation model. I - The no-mountain model

A linear primitive-equation stationary-wave model is developed to describe Northern Hemisphere winter stationary eddies on a sphere, and simulation results are compared with the predictions of a flat-lower-boundary GCM in extensive graphs and maps. The linear model is shown to reproduce the general behavior of the GCM well, with 10-30-percent underestimation of amplitudes in the Northern Hemisphere extratropical upper troposphere. The response to global heating plus transient eddy flux convergences is then decomposed into responses to total tropical forcing (heating plus transients) and total extratropical forcing (upper and lower tropospheric transients and latent, sensible, and radiative heating). The quality of the simulation is found to deteriorate rapidly if transient forcing is omitted, although somewhat better results are obtained by substituting thermal damping for transient forcing. The importance of low-level-transient parameterization for theories of time-mean low-level flow is emphasized.

Nigam, Sumant↗
NASA NTRS2008

Transient Three-Dimensional Startup Side Load Analysis of a Regeneratively Cooled Nozzle

The objective of this effort is to develop a computational methodology to capture the startup side load physics and to anchor the computed aerodynamic side loads with the available data from a regeneratively cooled, high-aspect-ratio nozzle, hot-fired at sea level. The computational methodology is based on an unstructured-grid, pressure-based, reacting flow computational fluid dynamics and heat transfer formulation, a transient 5 s inlet history based on an engine system simulation, and a wall temperature distribution to reflect the effect of regenerative cooling. To understand the effect of regenerative wall cooling, two transient computations were performed using the boundary conditions of adiabatic and cooled walls, respectively. The results show that three types of shock evolution are responsible for side loads: generation of combustion wave; transitions among free-shock separation, restricted-shock separation, and simultaneous free-shock and restricted shock separations; along with the pulsation of shocks across the lip, although the combustion wave is commonly eliminated with the sparklers during actual test. The test measured two side load events: a secondary and lower side load, followed by a primary and peak side load. Results from both wall boundary conditions captured the free-shock separation to restricted-shock separation transition with computed side loads matching the measured secondary side load. For the primary side load, the cooled wall transient produced restricted-shock pulsation across the nozzle lip with peak side load matching that of the test, while the adiabatic wall transient captured shock transitions and free-shock pulsation across the lip with computed peak side load 50% lower than that of the measurement. The computed dominant pulsation frequency of the cooled wall nozzle agrees with that of a separate test, while that of the adiabatic wall nozzle is more than 50% lower than that of the measurement. The computed teepee-like formation and the tangential motion of the shocks during lip pulsation also qualitatively agree with those of test observations. Moreover, a third transient computation was performed with a proportionately shortened 1 s sequence, and lower side loads were obtained with the higher ramp rate.

Wang, Ten-See↗
DOE OSTI2022

Multi-scale simulation of high burnup UO2 nuclear fuel during loss-of-coolant accident conditions

To improve the economics of commercial nuclear energy generation, U.S. utilities are currently seeking licensing approval to operate UO2 fuel to higher burnups. One significant safety issue that must be addressed to obtain approval is the potential for fine fragmentation/pulverization of the fuel during a loss-of-coolant accident (LOCA). It has been hypothesized that pulverization is caused by the rapid increase of pressure in fission gas bubbles in the high burnup region of the fuel during a LOCA. To better understand this phenomenon, a novel phase-field model of the fission gas bubble microstructure in UO2 has been developed and implemented in Idaho National Laboratory (INL)'s Marmot application for phase-field simulation of nuclear materials. Simulations of the bubble response to steady-state and transient conditions were conducted, and the results were used to inform a mechanistic model of pulverization in BISON, INL’s fuel performance simulation code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Creep performance of transient liquid phase bonded haynes 230 alloy

In this study, creep performance of a potential heat exchanger alloy, Haynes 230 was investigated. To simulate mechanical performance of a heat exchanger, sheet material of Haynes 230 transient liquid phase (TLP) bonded together. Sheets were plated on both sides with a Ni–P compound, subsequently stacked and TLP bonded together. Duplicate creep tests show low scatter in creep results, meaning a good repeatable joint was formed by TLP bonding. Overall creep performance showed less ductility and shorter times to failure than specimens originating from plate material. Lastly, evolution of a phosphorous compound at the bond line was observed and thought to potentially reduce ductility of the TLP bonded alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulations of the Collision between a Debris Stream and an Outer Dusty Torus: A Possible Channel for Forming a Fast-rise and Long-delay Radio Outburst in Tidal Disruption Events

A geometrically thick dusty torus structure is believed to exist in the nuclear regions of galaxies (especially in active galactic nuclei). The debris stream from a tidal disruption event (TDE) will possibly collide with the dusty torus and produce a transient flare. We perform three-dimensional hydrodynamic simulations to model the dynamical evolution of the interaction between unbound debris and a dusty torus. During the continuous interaction, shocked material will be spilled out from the interaction region and form an outflow. We calculate the temporal evolution of synchrotron emission by assuming that the shock accelerates a fraction of electrons in the outflow into a nonthermal distribution. We find that radio emission from the debris–torus collision generates a steep-rise and slow-decline radio light curve due to the sharp edge and dense gas of the dusty torus, where the radio outburst is delayed relative to the main optical/X-ray outburst by several years or even several tens of years. We apply our model to a TDE that happened in a narrow-line Seyfert I galaxy (PS16dtm), where both the radio spectrum and the light curve can be roughly reproduced. Future high-sensitivity, wide-field-of-view radio surveys have the opportunity to detect more such radio flares.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1998

A Second Law Based Unstructured Finite Volume Procedure for Generalized Flow Simulation

An unstructured finite volume procedure has been developed for steady and transient thermo-fluid dynamic analysis of fluid systems and components. The procedure is applicable for a flow network consisting of pipes and various fittings where flow is assumed to be one dimensional. It can also be used to simulate flow in a component by modeling a multi-dimensional flow using the same numerical scheme. The flow domain is discretized into a number of interconnected control volumes located arbitrarily in space. The conservation equations for each control volume account for the transport of mass, momentum and entropy from the neighboring control volumes. In addition, they also include the sources of each conserved variable and time dependent terms. The source term of entropy equation contains entropy generation due to heat transfer and fluid friction. Thermodynamic properties are computed from the equation of state of a real fluid. The system of equations is solved by a hybrid numerical method which is a combination of simultaneous Newton-Raphson and successive substitution schemes. The paper also describes the application and verification of the procedure by comparing its predictions with the analytical and numerical solution of several benchmark problems.

Majumdar, Alok↗
DOE OSTI2021

BISON microstructure-based pulverization criterion in high burnup structure

To improve the economics of commercial nuclear power production, utilities are seeking to increase the allowable burnup limit of UO$_2$ fuel. One of the main factors that contributes to the current burnup limit of 62 GWd/MTU in commercial light water reactors (LWRs) is the risk of fine fragmentation or pulverization during a loss of coolant accident (LOCA). Pulverization primarily occurs at high burnups, especially when the high burnup structure (HBS) has formed. To allow the industry to pursue increased burnup and develop mitigation strategies, it is essential to have improved capability to predict the onset of pulverization. However, the mechanism of pulverization is not well understood, and the existing predictive capabilities implemented in the BISON fuel performance code are empirical in nature. In this report, mesoscale simulations are used to improve understanding of the formation mechanism of the HBS and how it responds during a LOCA transient, and inform development of a BISON pulverization criterion. A phase-field model was used to simulate the evolution of bubble pressure as a result of HBS formation. The simulations showed that gas atoms diffuse from grain interiors to the new grain boundaries created during HBS formation, and diffuse rapidly along these grain boundaries to reach existing bubbles. This causes an increase in bubble pressure in existing bubbles, leading to bubble growth during steady-state operation. To simulate the response of HBS bubbles to a LOCA transient, a newly developed phase-field model was used; in agreement with preliminary results from FY20, bubble size did not change significantly during the duration of the transient. A phase-field fracture model was used to study fragmentation patterns in the HBS, including using input from the phase-field model as initial conditions. Phase-field fracture simulations were used to determine a pulverization criterion for BISON. A function for the critical pressure for pulverization to occur was fit to data from the phase-field fracture simulations, and this function was implemented as a material property in BISON. For comparison, an analytical criterion for pulverization was developed and implemented within the same material property.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1991

FD-TD numerical simulation of an entire lightning strike on the C160 aircraft

Experimental transient electromagnetic field measurements were performed on a Transall C160 aircraft during in-flight lightning strikes. The data allow a test of the predictive capabilities of a three dimensional time domain finite difference code (ALICE) developed at ONERA in order to investigate lightning-aircraft interactions. Using a transfer function technique in the 3D code, it is shown that a bi-leader attached to an aircraft can be simulated by a linear model, and so the electromagnetic fields can be calculated anywhere on the vehicle. Comparison of experimental and numerical results were made for several lightning strikes. Skin current density and electromagnetic field distributions are discussed in detail.

Alliot, J. C.↗
NASA NTRS1983

Modeling of transient two-component flow using a four-point implicit method

The four-point, centered implicit scheme that is extensively used in open channel flow simulation is shown to be applicable to rapid and slow pressure transient problems in conduits with nearly single phase and two-phase flows. It is only necessary to choose the proper weighting factor value, theta, of the Courant number. For rapid pressure transients such as waterhammer, the implicit method can yield reasonable results with limited numerical dispersion and attenuation if theta is only slightly greater than the critical value of 0.5. For slower pressure gradients in single and two-phase flows, reasonable numerical solutions may be achieved for Courant number values as high as 20.

Wiggert, D. C.↗
DOE OSTI2024

Analysis of Grid-Forming Inverter Controls for Grid-Connected and Islanded Microgrid Integration

Autonomous grid-forming (GFM) inverter testbeds with scalable platforms have attracted interest recently. In this study, a self-synchronized universal droop controller (SUDC) was adopted, tested, and scaled in a small network and a test feeder using a real-time simulation tool to operate microgrids without synchronous generators. We presented a novel GFM inverter control adoption to better understand the dynamic behavior of the inverters and their scalability, which can impact the distribution system (DS). This paper provides a steady-state and transient analysis of the GFM power inverter controller via simulation to better understand voltage and frequency stabilization and ensure that the critical electric loads are not affected during a prolonged power outage. The controllers of the GFM inverter are simulated in HYPERSIM to examine voltage and frequency fluctuations. This analysis includes assessing the black start capability for photovoltaic microgrids, both grid-connected and islanded, during transient fault conditions. The high photovoltaic PV penetration levels open exciting opportunities and challenges for the DS. The GFM inverter control demonstrated appropriate response times for synchronization, connection, and disconnection to the grid. The DS has become more resilient and independent of fossil fuels by increasing the penetration of inverter-based distributed energy resources (DERs).

Ward, Laura (ORCID:0000000345501581)↗
DOE OSTI2021

Temperature Effects on Droplet Oscillation Decay with Application to Fuel Property Measurement

Observation of oscillation decay in droplets has been shown to be an effective approach in determining physical properties such as viscosity and surface tension for emerging biofuels using μl quantities. Herein this work extends the approach to higher temperatures relevant to fuel injection conditions for internal combustion engines. Experiments are conducted that use high-resolution strobed imaging of moving, heated, μm-sized, fuel droplets to capture shape oscillation decay through image processing and analysis. Two fuels are investigated, iso-butanol and a primary reference fuel (PRF 84), which is a mixture of iso-octane and n-heptane. A piezoelectric droplet generator is used to generate a continuous train of single droplets, which are given an initial perturbation and subsequently undergo damped oscillations, captured using strobed imaging. Surface tension and viscosity calculated at four different temperatures ranging from 30°C-55°C using the frequency and decay time associated with the fundamental mode are found to be within ~ 10% of reference values obtained from the literature. Complementary numerical simulations are performed that utilize a volume-of-fluid approach to track transient droplet oscillation phenomena along with heat and mass transfer in a 2D axisymmetric domain. Simulations, where droplet size and temperature can be independently varied, capture an expected decrease in oscillation frequency and increase in decay time, with increase in fuel temperature. The simulations are further used to study the relative contribution to deviations in surface tension and viscosity predictions due to heat and mass transfer from the droplet, as well as viscous effects violating the inviscid flow assumption in the droplet oscillation theory. Mass loss effects are found to be negligible. Temperature change due to heat transfer is found to have the next highest sensitivity, particularly for the more volatile PRF 84, which has a lower viscosity and associated Ohnesorge number. For isobutanol, which has a higher viscosity and Ohnesorge number, viscous effects contribute the most to deviation in fuel property predictions.

09 BIOMASS FUELS↗
DOE OSTI2024

Phonons reveal coupled cholesterol-lipid dynamics in ternary membranes

Experimental studies of collective dynamics in lipid bilayers have been challenging due to the energy resolution required to observe these low-energy phonon-like modes. However, inelastic X-ray scattering (IXS) measurements — a technique for probing vibrations in soft and biological materials — are now possible with sub-meV resolution, permitting direct observation of low energy, phonon-like modes in lipid membranes. Here, in this study, IXS measurements with sub-meV energy resolution reveal a low-energy optic-like phonon mode at roughly 3 meV in the liquid-ordered (L o ) and liquid-disordered (L d ) phases of a ternary lipid mixture. This mode is only observed experimentally at momentum transfers greater than 5 nm -1 in the L o system. A similar gapped mode is also observed in all-atom molecular dynamics (MD) simulations of the same mixture, indicating that the simulations accurately represnt the fast, collective dynamics in the L o phase. It’s optical nature and the Q range of the gap together suggest that the observed mode is due to the coupled motion of cholesterol-lipid pairs, separated by several hydrocarbon chains within the membrane plane. Analysis of the simulations provide molecular insight into the origin of the mode in transient, nanoscale substructures of hexagonally packed hydrocarbon chains. This nanoscale hexagonal packing was previously reported based on molecular dynamics simulations and later by NMR measurements. Here, however, the integration of IXS and MD simulations identifies a new signature of the L o substructure in the collective lipid dynamics, thanks to the recent confluence of IXS sensitivity and MD simulation capabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Simulations of graphite boronization: A molecular dynamics study of amorphization resulting from bombardment

The molecular dynamics code LAMMPS was used to simulate the bombardment of a graphite structure by atomic boron with impact energies ranging from 50–250 eV. The transient structural evolution, penetration depth, and amorphous layer thickness were analyzed. Simulations show that larger impact energies lead to a greater volume of amorphization and penetration of boron, but that the growth rate of the amorphous layer decreases with increasing fluence. Furthermore, the change in surface chemistry of the amorphized structures was studied using the ReaxFF formalism, which found that the amorphization process introduces dangling bonds thus increasing reactivity in the amorphous region.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Wellbore integrity assessment with casing-based advanced sensing

Wellbore integrity is of paramount importance to subsurface resource extraction, energy storage 1 and waste disposal. After installation, well casing and cement are subject to mechanical stress due to near-well pressure changes and fluid induced corrosion. This is exacerbated for geothermal wells where produced fluid is at high temperature and corrosive. The current state-of-the-art technologies for wellbore integrity assessments are an array of cased hole logging tools. Wireline deployed acoustic, electromagnetic and mechanical tools are all available to inspect steel casing corrosion and casing-cement bond and these tools can provide high-resolution assessment of borehole conditions. They are intrusive, however, in terms of borehole preparation and interruption to the normal operation of the wells, and not suitable for high temperature or highly deviated well deployments. In addition, these measurements are performed infrequently due to high cost, and are therefore incapable of providing frequent data to better predict borehole degradation trajectory, which can help provide early warning of potential borehole failures. For this project we are developing a suite of novel, non-invasive, casing based tools for wellbore integrity assessment, combining fast/low cost screening with higher-precision investigation. Our approach is based on monitoring the response of the casing when energized at the wellhead, thereby interrogating the casing without well intervention. Lab, field and numerical approaches are used in our study. During the early stage of the research, we focus on numerical simulations, which have shown the sensitivity of the low frequency electromagnetic (EM) signals to changes in borehole depths and have successfully tested the concept at a field site with different length well casings. Initial seismic modeling efforts have also demonstrated our capability to simulate seismic tube wave and seismic field alterations due to borehole breakage and associated fluid leakage. Further numerical, laboratory and field experiments are underway for additional technology sensitivity analysis, particularly the transient EM/Seismic reflectometry methods, data acquisition optimization, and numerical simulation improvements.

Wilt, Michael↗
DOE OSTI2020

Stereoscopic detection of hot spots in superfluid 4 He (He II) for accelerator-cavity diagnosis

Superconducting radio-frequency (SRF) cavities cooled by superfluid 4He (He II) are building blocks of many modern particle accelerators due to their high quality factor. However, Joule heating from sub-millimeter surface defects on cavities can lead to cavity quenching, which limits the maximum acceleration gradient of the accelerators. Developing a non-contacting detection technology to accurately locate these surface hot spots is the key to improve the performance of SRF cavities and hence the accelerators. In a recent proof-of-concept experiment, we demonstrated that a molecular tagging velocimetry (MTV) technique based on the tracking of a He molecular tracer line created nearby a surface hot spot in He II can be utilized to locate the hot spot. In order to make this technique practically useful, here we describe our further development of a stereoscopic MTV setup for tracking the tracer line’s motion in three-dimensional (3D) space. We simulate a quench spot by applying a transient voltage pulse to a small heater mounted on a substrate plate. Images of the drifted tracer line, taken with two cameras from orthogonal directions, are used to reconstruct the line profile in 3D space. Here, a new algorithm for analyzing the 3D line profile is developed, which incorporates the finite size effect of the heater. We show that the center location of the heater can be reproduced on the substrate surface with an uncertainty of only a few hundred microns, thereby proving the practicability of this method.

43 PARTICLE ACCELERATORS↗