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,387 records · Page 77

Determination of Thermal State of Charge in Solar Heat Receivers

The research project at Howard University seeks to develop analytical and numerical capabilities to study heat transfer and fluid flow characteristics, and the prediction of the performance of solar heat receivers for space applications. Specifically, the study seeks to elucidate the effects of internal and external thermal radiation, geometrical and applicable dimensionless parameters on the overall heat transfer in space solar heat receivers. Over the last year, a procedure for the characterization of the state-of-charge (SOC) in solar heat receivers for space applications has been developed. By identifying the various factors that affect the SOC, a dimensional analysis is performed resulting in a number of dimensionless groups of parameters. Although not accomplished during the first phase of the research, data generated from a thermal simulation program can be used to determine values of the dimensionless parameters and the state-of-charge and thereby obtain a correlation for the SOC. The simulation program selected for the purpose is HOTTube, a thermal numerical computer code based on a transient time-explicit, axisymmetric model of the total solar heat receiver. Simulation results obtained with the computer program are presented the minimum and maximum insolation orbits. In the absence of any validation of the code with experimental data, results from HOTTube appear reasonable qualitatively in representing the physical situations modeled.

Glakpe, E. K.↗

The Sensitivity of African Easterly Waves to Eastern Tropical Atlantic Sea-Surface Temperatures

The results of two regional atmospheric model simulations are compared to assess the influence of the eastern tropical Atlantic sea-surface temperature maximum on local precipitation, transient easterly waves and the West African summer monsoon. Both model simulations were initialized with reanalysis 2 data (US National Center for Environmental Prediction and Department of Energy) on 15 May 2006 and extended through 6 October 2006, forced by synchronous reanalysis 2 lateral boundary conditions introduced four times daily. One simulation uses 2006 reanalysis 2 sea-surface temperatures, also updated four times daily, while the second simulation considers ocean forcing absent the sea-surface temperature maximum, achieved here by subtracting 3 K at every ocean grid point between 0 and 15 N during the entire simulation. The simulation with 2006 sea-surface temperature forcing produces a realistic distribution of June-September mean precipitation and realistic westward propagating swaths of maximum rainfall, based on validation against Tropical Rainfall Measuring Mission (TRMM) estimates. The simulation without the sea-surface temperature maximum produces only 57% of the control June-September total precipitation over the eastern tropical Atlantic and about 83% of the Sahel precipitation. The simulation with warmer ocean temperatures generates generally stronger circulation, which in turn enhances precipitation by increasing moisture convergence. Some local precipitation enhancement is also attributed to lower vertical thermal stability above the warm water. The study shows that the eastern tropical Atlantic sea-surface temperature maximum enhances the strength of transient easterly waves and broadens the spatial extent of associated precipitation. However, large-scale circulation and its interaction with the African continent, and not sea-surface temperatures, control the timing and trajectories of the waves.

Druyan, Leonard M.↗

Analysis of the transient response of nuclear spins in GaAs with/without nuclear magnetic resonance

As an alternative to studying the steady-state responses of nuclear spins in solid state systems, working within a transient-state framework can reveal interesting phenomena. The response of nuclear spins in GaAs to a changing magnetic field was analyzed based on the time evolution of nuclear spin temperature. Simulation results well reproduced our experimental results for the transient oblique Hanle signals observed in an all-electrical spin injection device. The analysis showed that the so called dynamic nuclear polarization can be treated as a cooling tool for the nuclear spins: It works as a provider to exchange spin angular momentum between polarized electron spins and nuclear spins through the hyperfine interaction, leading to an increase in the nuclear polarization. In addition, a time-delay of the nuclear spin temperature with a fast sweep of the external magnetic field produces a possible transient state for the nuclear spin polarization. On the other hand, the nuclear magnetic resonance acts as a heating tool for a nuclear spin system. This causes the nuclear spin temperature to jump to infinity: i.e., the average nuclear spins along with the nuclear field vanish at resonant fields of {sup 75}As, {sup 69}Ga and {sup 71}Ga, showing an interesting step-dip structure in the oblique Hanle signals. These analyses provide a quantitative understanding of nuclear spin dynamics in semiconductors for application in future computation processing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An Investigation Regarding the Effects of Volume Dynamics on Gas Turbine Engine Performance and Operability

This work explores the impact of adding volume dynamics to an electrified turbofan engine model. The zero-dimensional single-input single-output volume elements capture the physics of fluid compressibility of flow in the engine core gas path turbomachinery during transient operation with greater fidelity. Some background regarding the motivation for this effort is provided. A model of the volume element is developed from first principles. The baseline engine model is described and an explanation is given of the modifications made to incorporate the volume elements into the new model. Details of the simulated inputs, designed to elicit exaggerated dynamic responses and to facilitate comparison between the baseline and modified engine models, are explicated. The results of these simulations are presented and discussed, demonstrating that the transient performance and operability of the modified engine model, measured by parameters such as shaft speeds, thrust, and high-pressure compressor (HPC) stall margin, are impacted by the additional dynamics. Commentary is provided regarding potential explanations for the differences in behavior exhibited by the two models. Finally, suggestions are offered for future work that might improve the volume models developed in this effort and leverage the improved fidelity enabled by the modeling of turbomachinery volume dynamics.

Volume dynamics↗

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↗

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.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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.↗

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)↗