Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient Modeling”

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 937 records · Page 52

Current status of droplet and liquid combustion

The present understanding of spray combustion in rocket engine, gas turbine, Diesel engine and industrial furnace applications is reviewed. In some cases, spray combustion can be modeled by ignoring the details of spray evaporation and treating the system as a gaseous diffusion flame; however, in many circumstances, this simplification is not adequate and turbulent two-phase flow must be considered. The behavior of individual droplets is a necessary component of two-phase models and recent work on transient droplet evaporation, ignition and combustion is considered, along with a discussion of important simplifying assumptions involved with modeling these processes. Methods of modeling spray evaporation and combustion processes are also discussed including: one-dimensional models for rocket engine and prevaporized combustion systems, lumped zone models (utilizing well-stirred reactor and plug flow regions) for gas turbine and furnace systems, locally homogeneous turbulent models, and two-phase models. The review highlights the need for improved injector characterization methods, more information of droplet transport characteristics in turbulent flow and continued development of more complete two-phase turbulent models

Faeth, G. M.↗

Validating Solar and Heliospheric Models at the CCMC

The Community Coordinated Modeling Center (CCMC) hosts a growing number of models of the ambient and transient corona and heliosphere which are ultimately intended for use in space weather forecasting. independent validation of these models is a critical step in their development as potential forecasting tools for the space weather operations community. In this poster we report on validation studies of these models, all of which are also available for use by the research community through our runs-on-request system.

MacNeice, P. J.↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Incorporation of Plasticity and Damage Into an Orthotropic Three-Dimensional Model with Tabulated Input Suitable for Use in Composite Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive industries. While there are several composite material models currently available within commercial transient dynamic finite element codes, several features have been identified as being lacking in the currently available material models that could substantially enhance the predictive capability of the impact simulations. A specific desired feature pertains to the incorporation of both plasticity and damage within the material model. Another desired feature relates to using experimentally based tabulated stress-strain input to define the evolution of plasticity and damage as opposed to specifying discrete input properties (such as modulus and strength) and employing analytical functions to track the response of the material. To begin to address these needs, a combined plasticity and damage model suitable for use with both solid and shell elements is being developed for implementation within the commercial code LS-DYNA. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain-hardening based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. The effective plastic strain is computed by using the non-associative flow rule in combination with appropriate numerical methods. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used, in which a load in one direction results in a stiffness reduction in multiple coordinate directions. A specific laminated composite is examined to demonstrate the process of characterizing and analyzing the response of a composite using the developed model.

Impact↗

Evaluation of Inelastic Constitutive Models for Nonlinear Structural Analysis

The influence of inelastic material models on computed stress-strain states, and therefore predicted lives, was studied for thermomechanically loaded structures. Nonlinear structural analyses were performed on a fatigue specimen which was subjected to thermal cycling in fluidized beds and on a mechanically load cycled benchmark notch specimen. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic-kinematic, combined plus transient creep) were exercised. Of the plasticity models, kinematic hardening gave results most consistent with experimental observations. Life predictions using the computed strain histories at the critical location with a Strainrange Partitioning approach considerably overpredicted the crack initiation life of the thermal fatigue specimen.

Kaufman, A.↗

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Multiscale Thermal-hydraulic analysis of the MARVEL microreactor using a coupled SAM and SCM simulation. (PPTX)

This is .pptx document presenting a summary of the paper of the same name that has been submitted and accepted to NURETH-21: MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor's primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL?s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Grid-Supportive Loads - A New Approach to Increasing Renewable Energy in Power Systems

This paper demonstrates the potential of inverter-based loads to support grid reliability during power system transients thereby enabling reliable integration of renewable energy in power systems. Such loads are referred to in this paper as grid-supportive loads (GSLs). A new GSL model is developed that simulates the transient response capabilities that can be programmed in electronic loads. The model’s design enables it to be easily integrated in widely used commercial power system transient analysis software. Theoretical expressions are derived that explain the workings of the GSL model. The performance, numerical stability, and impact of the GSL model is validated on 9-bus and 2000-bus synthetic power system models using generator tripping and bus fault disturbances. Results on the 2000 bus system show that in the absence of frequency support from wind/solar generation resources, just 20% of loads with grid-supportive capabilities can improve frequency response by up to 2000 MW/0.1 Hz and reduce deviation in frequency at nadir by up to 60% compared to the situation when GSLs are absent. Power system reliability also improves under fault events. Here, it is further shown that GSLs can aid in integrating more renewable generation without degrading the overall transient response of the power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Swift /UVOT discovery of Swift J221951−484240: a UV luminous ambiguous nuclear transient

We report the discovery of Swift J221951−484240 (hereafter: J221951), a luminous slow-evolving blue transient that was detected by the Neil Gehrels Swift Observatory Ultraviolet/Optical Telescope (Swift/UVOT) during the follow-up of gravitational wave alert S190930t, to which it is unrelated. Swift/UVOT photometry shows the UV spectral energy distribution of the transient to be well modelled by a slowly shrinking blackbody with an approximately constant temperature of T ∼ 2.5 × 10 4 K. At a redshift z = 0.5205, J221951 had a peak absolute magnitude of M u,AB = −23 mag, peak bolometric luminosity Lmax = 1.1 x 10 45 ergs -1 and a total radiated energy of E > 2.6 × 10 52 erg. The archival Wide-field Infrared Survey Explorer IR photometry shows a slow rise prior to a peak near the discovery date. Spectroscopic UV observations display broad absorption lines in N V and O VI, pointing towards an outflow at coronal temperatures. The lack of emission in the higher H α lines, N I and other neutral lines is consistent with a viewing angle close to the plane of the accretion or debris disc. The origin of J221951 cannot be determined with certainty but has properties consistent with a tidal disruption event and the turn-on of an active galactic nucleus.

79 ASTRONOMY AND ASTROPHYSICS↗

NUMERICAL MODELING OF A SOLID OXIDE FUEL CELL FOR USE IN REAL-TIME SIMULATION AND CYBER-PHYSICAL SYSTEMS

Cyber-physical systems provide a mechanism with which to investigate the physical phenomena and behavior of traditionally cost-prohibitive or otherwise fragile equipment. For the National Energy Technology Laboratory (NETL), this approach resulted in the Hybrid Performance (Hyper) facility which features a gas turbine-SOFC hybrid cycle utilizing real turbomachinery and a simulated SOFC stack. This allows for the investigation of combined cycle performance and control strategies, in an exhaustive manner, both without fear of destroying delicate state-of-the-art fuel cells, and with the full accuracy of real-world turbomachinery. Issues arose between the transient response of the SOFC model being limited to a sample time of 80 milliseconds, due to the calculation time of the SOFC model taking on average 40 milliseconds to calculate for a given timestep with spikes in calculation time reaching the 80 millisecond threshold. In order to be able to match the speed of transients from the turbomachinery and likewise better discern transient behavior, it was determined that the SOFC model must be optimized to operate at a sample time of 5 milliseconds. Therefore, it is necessary to optimize the SOFC model in order to decrease the calculation time from around 40 milliseconds, down to at the most 5 milliseconds. To do this, both the electrochemical algorithm and the thermal algorithm used to simulate the physical behavior of the SOFC are investigated to determine where improvements can be made. To this end the rootfinding numerical recipes of the electrochemical algorithm are investigated as the complex electrochemistry requires a highly iterative nested dual convergence loop to resolve the voltage-current relationship, and likewise the temporal discretization of the thermal algorithm is modified for the sake of higher accuracy and stability. Ultimately the new electrochemical algorithm featuring higher order rootfinding schemes proves to be efficient enough to reach the sub 5 millisecond target, signifying an order of magnitude reduction in calculation time, and when coupled with the new temporal discretization similar calculation time characteristics show that a fully implicit, higher order temporal discretization can also successfully be used if desired. Ultimately this result means that the cyber-physical simulation system can operate at higher sample rates, and resolve transient events at significantly higher resolution and fidelity.

Arias, Jesus↗

Hybrid MD-Nernst Planck Model of Alpha-hemolysin Conductance Properties

Motivated by experiments in which an applied electric field translocates polynucleotides through an alpha-hemolysin protein channel causing ionic current transient blockade, a hybrid simulation model is proposed to predict the conductance properties of the open channel. Time scales corresponding to ion permeation processes are reached using the Poisson-Nemst-Planck (PNP) electro-diffusion model in which both solvent and local ion concentrations are represented as a continuum. The diffusion coefficients of the ions (K(+) and Cl(-)) input in the PNP model are, however, calculated from all-atom molecular dynamics (MD). In the MD simulations, a reduced representation of the channel is used. The channel is solvated in a 1 M KCI solution, and an external electric field is applied. The pore specific diffusion coefficients for both ionic species are reduced 5-7 times in comparison to bulk values. Significant statistical variations (17-45%) of the pore-ions diffusivities are observed. Within the statistics, the ionic diffusivities remain invariable for a range of external applied voltages between 30 and 240mV. In the 2D-PNP calculations, the pore stem is approximated by a smooth cylinder of radius approx. 9A with two constriction blocks where the radius is reduced to approx. 6A. The electrostatic potential includes the contribution from the atomistic charges. The MD-PNP model shows that the atomic charges are responsible for the rectifying behaviour and for the slight anion selectivity of the a-hemolysin pore. Independent of the hierarchy between the anion and cation diffusivities, the anionic contribution to the total ionic current will dominate. The predictions of the MD-PNP model are in good agreement with experimental data and give confidence in the present approach of bridging time scales by combining a microscopic and macroscopic model.

Cozmuta, Ioana↗

Steady state and transient temperature distributions in the human thigh covered with a cooling pad

An analytical and experimental study was done on the performance of cooling pads attached to a human thigh. Each cooling pad consisted of a long, water cooled tube formed into a serpentine shape with uniform spacing between the parallel sections. The analytical work developed a cylindrical model for the human thigh. The transient times predicted by this model ranged from 25 to 80 minutes, which is reasonably close to the experimental results. Calculated and measured steady state temperature profiles were in fair agreement. The transient times associated with a change from a high metabolic rate of 1800 Btu/hr (528 w) to a low level of 300 Btu/hr (88 w), were found to be about 120 minutes. A change from 300 Btu/hr (264 w) to 300 Btu/hr (88 w) resulted in 90 to 100 minute transients. However, the transient times for a change in metabolic rate in the opposite direction from 300 Btu/hr (88 w) to 1800 Btu/hr (528 w) were 40 to 60 minutes.

Leo, R. J.↗

Detailed mathematical models of a radioisotope thermoelectric generator.

Two new models for the design and performance analysis of RTG's are outlined in this paper. The first model assumes a small-signal transient-type calculational sequence that permits the separation of steady-state operation of the generator from its dynamic behavior. The second model uses a numerical (finite difference) solution of the performance equations of the RTG. Both models enable the investigation of transient and steady-state performance of RTG's. Simplifying assumptions have been kept to a minimum in the new RTG models and these models enable the inclusion of generator end losses, axial temperature gradients and heat interchange between thermoelements and thermal insulation in RTG performance calculations in a self-consistent manner.

Dewinter, F.↗

Some atmospheric scattering considerations relevant to BATSE: A model calculation

The orbiting Burst and Transient Source Experiement (BATSE) will locate gamma ray burst sources by analysis of the relative numbers of photons coming directly from a source and entering its prescribed array of detectors. In order to accurately locate burst sources it is thus necessary to identify and correct for any counts contributed by events other than direct entry by a mainstream photon. An effort is described which estimates the photon numbers which might be scattered into the BATSE detectors from interactions with the Earth atmosphere. A model was developed which yielded analytical expressions for single-scatter photon contributions in terms of source and satellite locations.

Young, John H.↗

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗

Drift reduced Landau fluid model for magnetized plasma turbulence simulations in BOUT++ framework

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture $n$ = 0 axisymmetric mode evolution in realistic tokamak configurations. As we report here, these improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Understanding the interplay between pilot fuel mixing and auto-ignition chemistry in hydrogen-enriched environment

The diesel-piloted dual-fuel compression ignition combustion strategy is well-suited to accelerate the decarbonization of transportation by adopting hydrogen as a renewable energy carrier into the existing internal combustion engine with minimal engine modifications. Despite the simplicity of engine modification, many questions remain unanswered regarding the optimal pilot injection strategy for reliable ignition with minimum pilot fuel consumption. The present study uses a single-cylinder heavy-duty optical engine to explore the phenomenology and underlying mechanisms governing the pilot fuel ignition and the subsequent combustion of a premixed hydrogen-air charge. The engine is operated in a dual-fuel mode with hydrogen premixed into the engine intake charge with a direct pilot injection of n-heptane as a diesel pilot fuel surrogate. Optical diagnostics used to visualize in-cylinder combustion phenomena include high-speed IR imaging of the pilot fuel spray evolution as well as high-speed HCHO* and OH* chemiluminescence as indicators of low-temperature and high-temperature heat release, respectively. Three pilot injection strategies are compared to explore the effects of pilot fuel mass, injection pressure, and injection duration on the probability and repeatability of successful ignition. The thermodynamic and imaging data analysis supported by zero-dimensional chemical kinetics simulations revealed a complex interplay between the physical and chemical processes governing the pilot fuel ignition process in a hydrogen containing charge. Hydrogen strongly inhibits the ignition of pilot fuel mixtures and therefore requires longer injection duration to create zones with sufficiently high pilot fuel concentration for successful ignition. Results show that ignition typically tends to rely on stochastic pockets with high pilot fuel concentration, which results in poor repeatability of combustion and frequent misfiring. In conclusion, this work has improved the understanding on how the unique chemical properties of hydrogen pose a challenge for maximization of hydrogen’s energy share in hydrogen dual-fuel engines and highlights a potential mitigation pathway.

33 ADVANCED PROPULSION SYSTEMS↗

Hunting for gamma-ray emission from fast radio bursts

Context. Fast radio bursts (FRBs) are a recently discovered class of GHz-band, ms-duration, Jansky-level-flux astrophysical transients. Although hundreds of models have been proposed so far for FRB progenitors (the most popular ones involve magnetars), their physical origin and emission mechanism are still a mystery, making them one of the most compelling problems in astrophysics. Aims. FRBs are caused by astrophysical processes that are not yet understood. Exploring their high-energy counterpart is crucial for constraining their origin and emission mechanism. Methods. Thanks to more than 13 years of gamma-ray data collected by the Large Area Telescope on board the Fermi Gamma-ray Space Telescope, and to more than 1000 FRB events (from 561 non-repeating and 22 repeating sources), one of the largest samples created thus far, we performed the largest and deepest search for high-energy emission from FRB sources to date (between 100 MeV and 1 TeV). In addition to the analysis involving individual FRB events on different timescales (from a few seconds up to several years), we performed, for the first time, a stacking analysis on the full sample of FRB events as well as a search for triplet photons in coincidence with the radio event. Results. We do not detect significant emission, reporting the most stringent constraints, on short timescales, for the FRB-like emission from SGR 1935+2154 with E γ <10 41 erg, corresponding to a factor η < 10 7 with respect to the emitted radio energy. Similarly, for the stacked signal of steady emission from all repeaters, the obtained upper limit (UL) on the FRBs luminosity (L γ < 1.6 × 10 43 erg s -1 is more than two orders of magnitude lower than those derived from the individual sources. Finally, no individual or triplet photons have been significantly associated with FRB events. We derived the LAT ms-sensitivity to be ~0.3 ph cm -2 s -1 and constrained the gamma-ray energy E γ,δT = 1 ms ≲ 10 47 (D L /150 Mpc) 2 erg, ruling out a gamma-ray-to-radio energy ratio greater than 10 9 on ms timescales. Conclusions. The results reported here represent the most stringent UL reported so far on the high-energy emission from FRBs on short and long time scales, as well as on cumulative emission and individual photon searches. While the origin of FRBs is still unclear, our work provides important constraints for FRB modelling, which might shed light on their emission mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗