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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS2010

Evaluation of Grounding Impedance of a Complex Lightning Protective System Using Earth Ground Clamp Measurements and ATP Modeling

A new Lightning Protection System (LPS) was designed and built at Launch Complex 39B (LC39B), at the Kennedy Space Center (KSC), Florida, which consists of a catenary wire system (at a height of about 181 meters above ground level) supported by three insulators installed atop three towers in a triangular configuration. A total of nine downconductors (each about 250 meters long, on average) are connected to the catenary wire system. Each of the nine downconductors is connected to a 7.62-meter radius circular counterpoise conductor with six equally spaced 6-meter long vertical grounding rods. Grounding requirements at LC39B call for all underground and above ground metallic piping, enclosures, raceways, and cable trays, within 7.62 meters of the counterpoise, to be bounded to the counterpoise, which results in a complex interconnected grounding system, given the many metallic piping, raceways, and cable trays that run in multiple direction around LC39B. The complexity of this grounding system makes the fall of potential method, which uses multiple metallic rods or stakes, unsuitable for measuring the grounding impedances of the downconductors. To calculate the downconductors grounding impedance, an Earth Ground Clamp (a stakeless grounding resistance measuring device) and a LPS Alternative Transient Program (ATP) model are used. The Earth Ground Clamp is used to measure the loop impedance plus the grounding impedance of each downconductor and the ATP model is used to calculate the loop impedance of each downconductor circuit. The grounding impedance of the downconductors is then calculated by subtracting the ATP calculated loop impedances from the Earth Ground Clamp measurements.

Mata, Carlos T.↗

NASA NTRS2005

The Importance of Detailed Component Simulations in the Feedsystem Development for a Two-Stage-to Orbit Reusable Launch Vehicle

To meet the requirements for the 2nd Generation Reusable Launch Vehicle (RLV), a unique propulsion feed system concept was identified using crossfeed between the booster and orbiter stages that could reduce the Two-Stage-to-Orbit (TSTO) vehicle weight and development cost by approximately 25%. A Main Propulsion System (MPS) crossfeed water demonstration test program was configured to address all the activities required to reduce the risks for the MPS crossfeed system. A transient, one-dimensional system simulation was developed for the subscale crossfeed water flow tests. To ensure accurate representation of the crossfeed valve's dynamics in the system model, a high-fidelity, three-dimensional, computational fluid-dynamics (CFD) model was employed. The results from the CFD model were used to specify the valve's flow characteristics in the system simulation. This yielded a crossfeed system model that was anchored to the specific valve hardware and achieved good agreement with the measured test data. These results allowed the transient models to be correlated and validated and used for full scale mission predictions. The full scale model simulations indicate crossfeed is ' viable with the system pressure disturbances at the crossfeed transition being less than experienced by the propulsion system during engine start and shutdown transients.

Mazurkivich, Pete↗

DOE OSTI2020

A Multi-Function AAA Algorithm Applied to Frequency Dependent Line Modeling

Modeling of power distribution system components that are valid for a wide range of frequencies are crucial for highly accurate modeling of electromagnetic transient (EMT) events. This has recently become of interest due to the improvements needed for the resilient operation of distribution systems. Vector fitting (VF) is a very popular and commonly used algorithm for wide band representations of power system components in EMT simulations. In this research, we present a new multi-input rational approximation algorithm (MIAAA) and illustrate its advantages with respect to VF using examples of approximations of admittance matrices discussed in the literature. We show that MIAAA not only outperforms VF in terms of achieving better accuracy using lesser number of poles, but also has no numerical issues achieving convergence. In contrast to VF, MIAAA is not sensitive to the location of input sample points and it does not require good estimates for the location of the desired approximation poles. The novelty of this research work is the use of recent mathematical results to solve existing challenges in distribution system modeling and to develop rational approximations for power system models that intend to be optimal in terms of accuracy and performance.

41 EE - Solar Energy Technologies Office (EE-4S)↗

DOE OSTI2021

A MIMO AAA Agorithm for Frequency Dependent Line Modeling

AAA algorithm↗

DOE OSTI2020

FY20 SAM Code Developments and Validations for Transient Safety Analysis of Advanced non-LWRs

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report summarizes major progress in SAM code development, capability enhancements, demonstration, and validation to support transient safety analysis of advanced non-LWRs. Rapid developments continued in fiscal year 2020 (FY20) to support various needs of the advanced reactor community, especially the NRC and industry on the licensing safety analysis of advanced reactor designs. Significant code changes were made to provide various capability enhancements, bug fixes, and user friendliness improvements. Major code updates are summarized in Section 1, while four important enhancements are detailed in Sections 2-5, including a multi-dimension flow model; reactivity feedback and decay heat models; control and trip system modeling, and additional fluid and solid thermophysical property models. Code validation activities in FY20 include using test data from the Fast Flux Test Facility (FFTF), the High Temperature Test Facility (HTTF), and several separate effects test facilities for pebble-bed modeling.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS1993

Heat Transfer from Radiatively Heated Material in a Low Reynolds Number Microgravity Environment

A mathematical model of the transient three-dimensional heat transfer between a slowly moving ambient gas stream and a thermally thick or thin flat surface heated by external radiation in a microgravity environment is presented. The problem is motivated in part by fire safety issues in spacecraft. The gas phase is represented by variable property convection-diffusion energy and mass conservation equations valid at low Reynolds numbers. The absence of gravity and low Reynolds number together permit the flow to be represented by a self-consistent velocity potential determined by the ambient velocity and the thermal expansion in the gas. The solid exchanges energy with the gas by conduction/convection and with the surroundings by surface absorption and re-emission of radiation. Heat conduction in the solid is assumed to be one dimensional at each point on the surface as a consequence of the limited times (of order of 10 seconds) of interest in these simulations. Despite the apparent simplicity of the model, the results show a complex thermally induced flow near the heated surface. The thermal exchange between the gas and solid produces an outward sourcelike flow upstream of the center of the irradiated area and a sinklike flow downstream. The responses of the temperature fields and the associated flows to changes in the intensity of the external radiation and the ambient velocity are discussed.

Yamashita, H.↗

DOE OSTI2024

Pre-Test Feasibility Study Of Reactivity Transient Testing On TRISO Fuel In The TREAT Reactor

Uranium oxycarbide (UCO)-bearing tri-structural isotropic (TRISO) particle fuels are expected to be used in numerous US commercial reactor applications within the next decade. In this work, we reviewed historical particle fuel transient experiments to identify gaps in the TRISO fuel performance transient testing. A BISON--Griffin modeling framework was then developed to conduct preliminary TRISO transient analyses and begin to address these gaps. The framework was demonstrated using limiting-case transient conditions from a prototypic high-temperature gas-cooled reactor. It was then applied to develop a matrix of experiments that could be performed in the Transient Reactor Test Facility (TREAT) to (1) evaluate UCO-fueled particle performance at moderate and high heat rates, (2) assess whether historical testing involving UO$_2$-fueled particles is applicable to modern UCO-fueled particles, (3) deconvolute the impacts of temperature and heat rate on particle transient response, and (4) collect the data needed for fuel performance model validation and/or further development.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

High Time Resolution Modelled and Measured New Particle Formation Rates During the 2023 Pittsburgh, PA Campaign

Particle size distributions are used to calculate nucleation rates during new particle formation (NPF). However, the low time resolution and poor sensitivity of instruments deployed to observe particle size distributions can lead to inaccurate estimates of nucleation rates. We calculated nucleation rates in Pittsburgh, PA from particle size distributions measured with the new high time resolution Condensation Particle Counters For Atmospheric Rapid Measurements (CPC FARM) and from gas precursor concentrations measured with a chemical ionization mass spectrometer. The CPC FARM nucleation rates were calculated by examining the appearance rate of particles in the smallest size bin at 1.6 nm. The nucleation rate of 1.6 nm determined every second resulted in positive and negative rates, indicating the influence of inhomogeneous mixing and advection of precursor gases and clusters. These effects are not currently considered when estimating nucleation rates from measurements. The nucleation rates from the CPC FARM were compared to the nucleation rates determined from a transient acid‐base reaction model using measured sulfuric acid concentration, strong base concentrations, and CPC detection efficiencies. The transient form of the model allowed the advection of heterogeneously mixed precursors to influence the nucleation rate. Good agreement between particle and gas‐based nucleation rates is observed during sulfuric acid‐driven NPF events. This study demonstrates that high time resolution measurements of short‐lived and rapidly formed particles are needed to accurately capture nucleation rates, which are likely influenced by advection of nucleation precursors and clusters in heterogeneous environments such as urban areas.

Cheng, D. [Carnegie Mellon Univ., Pittsburgh, PA (↗

DOE OSTI2026

Simulation and Experimental Validation of an Integrated Heat Pump – Thermal Energy Storage Using a Room-Temperature Phase Change Material

As the dependence on electrical heat pumps (HPs) and intermittent renewables increases, grid strains are expected to grow. This necessitates an energy storage system to reduce the mismatch between energy supply and demand. Thus, a proposed dual-mode commercially available 14.1 kW HP was integrated with a single 22°C phase change material (PCM) thermal storage system (TES) to load-shift both cooling and heating loads. The HP-TES system was manufactured and experimentally tested using a novel test matrix based on AHRI 210/240 psychrometric conditions. Furthermore, transient dual-mode system-level HP-TES models were developed in Modelica and validated using the experimental test conditions. Base HP cooling and heating experimental tests at ambient temperatures of 35°C and −8.3°C show that the modified HP-TES maintained the rated system capacity and performance. The HP-TES discharge provided approximately 30% and 50% reductions in cooling and heating demand, respectively. The transient HP-TES models predicted system capacity and total power input for discharge and recharge operating modes within ±4% mean percentage error, and recharge power input within ±2%, with maximum errors occurring at the equipment startup. During system operation, the sources of model deviations are first-order polynomial fits of the PCM digital scanning calorimetry (DSC) data and unaccounted supercooling in the PCM during solidification. Nonetheless, the model predictions agree with the experimental tests, demonstrating the availability of robust, accurate, and validated transient models that can be used for further validation and the development of system controls.

25 ENERGY STORAGE↗

NASA NTRS1997

DRA/NASA/ONERA Collaboration on Icing Research: Prediction of Airfoil Ice Accretion - Part 2

This report presents results from a joint study by DRA, NASA, and ONERA for the purpose of comparing, improving, and validating the aircraft icing computer codes developed by each agency. These codes are of three kinds: (1) water droplet trajectory prediction, (2) ice accretion modeling, and (3) transient electrothermal deicer analysis. In this joint study, the agencies compared their code predictions with each other and with experimental results. These comparison exercises were published in three technical reports, each with joint authorship. DRA published and had first authorship of Part 1 - Droplet Trajectory Calculations, NASA of Part 2 - Ice Accretion Prediction, and ONERA of Part 3 - Electrothermal Deicer Analysis. The results cover work done during the period from August 1986 to late 1991. As a result, all of the information in this report is dated. Where necessary, current information is provided to show the direction of current research. In this present report on ice accretion, each agency predicted ice shapes on two dimensional airfoils under icing conditions for which experimental ice shapes were available. In general, all three codes did a reasonable job of predicting the measured ice shapes. For any given experimental condition, one of the three codes predicted the general ice features (i.e., shape, impingement limits, mass of ice) somewhat better than did the other two. However, no single code consistently did better than the other two over the full range of conditions examined, which included rime, mixed, and glaze ice conditions. In several of the cases, DRA showed that the user's knowledge of icing can significantly improve the accuracy of the code prediction. Rime ice predictions were reasonably accurate and consistent among the codes, because droplets freeze on impact and the freezing model is simple. Glaze ice predictions were less accurate and less consistent among the codes, because the freezing model is more complex and is critically dependent upon unsubstantiated heat transfer and surface roughness models. Thus, heat transfer prediction methods used in the codes became the subject for a separate study in this report to compare predicted heat transfer coefficients with a limited experimental database of heat transfer coefficients for cylinders with simulated glaze and rime ice shapes. The codes did a good job of predicting heat transfer coefficients near the stagnation region of the ice shapes. But in the region of the ice horns, all three codes predicted heat transfer coefficients considerably higher than the measured values. An important conclusion of this study is that further research is needed to understand the finer detail of of the glaze ice accretion process and to develop improved glaze ice accretion models.

Wright, William B.↗

DOE OSTI2023

Assessment of Sodium Thermal Stratification Models Utilizing the TSTF Benchmark

As a result of certain transient scenarios, a thermally stratified layer of liquid sodium can develop in the bulk coolant volumes of a sodium-cooled fast reactor (SFR). In addition to the effects a stratification layer has on the temperature of the heat transport system, a stratification layer can also influence the transition to and establishment of natural circulation flow, which plays an important role in passive cooling and the inherent safety of a pool-type SFR. Therefore, the ability to accurately capture thermal stratification phenomena is important when demonstrating the safety basis of a pool-type SFR during transient sequences. The present work assesses various computational models with different fidelities in their ability to predict thermal stratification in the upper plenum of an SFR. Each computational model will be assessed using the data generated at the Thermal Stratification Test Facility (TSTF) located at the University of Wisconsin-Madison. Using measured flow rate and inlet temperature data, the measured temperature distributions of the tests are compared to the predictions of the lumped volume-based models in SAS4A/SASSYS-1, a 1D-based model in SAM, and a 3-D computational fluid dynamics (CFD) model using STAR-CCM+. The relative performance of the various computational methods is assessed with respect to key metrics such as bulk coolant temperature distribution and plenum exit temperature. A total of eight tests are analyzed, covering different combinations of flow rates (3 and 10 GPM) and upper internal structure (UIS) configurations (none, solid, porous, and open) The perfect mixing model of SAS4A/SASSYS-1 provides the highest accuracy when the flow rate is high and there is no UIS in the test vessel, as high flow rate injection promotes thermal mixing of the sodium in the test vessel. For most of the analyzed tests, the stratified volume model of SAS4A/SASSYS-1 is able to predict the delay in the outlet temperature drop and temperature distribution in the test vessel by a small number of layers to represent thermal stratification. However, the stratified volume model can only simulate a maximum of three temperature layers within a volume and when a layer approaches the elevation of the outlet, the predicted outlet temperature can demonstrate rapid, non-physical changes. The 1-D axial mixing model of SAM provides results that agree reasonably well with the measured data in the prediction of the temporal evolution of the outlet temperature with the exception of the case with a high flow rate and no UIS. The SAM 1-D model has a similar level of accuracy to CFD results when it comes to predicting the outlet temperature. CFD shows overall good agreement in predicting the temperature distribution in the test vessel and outlet temperature. As CFD can model the test vessel geometry in detail, it performs well in the cases of complex geometries such as tests that included a UIS and internal flow through the UIS resulting in active mixing of the coolant in the test vessel. Each of the models discussed in the present work has the potential to be useful during the various stages of reactor design, analysis, and licensing. The lumped-volume approach can be applied for fast turnaround safety calculations to obtain overall reactor behavior during transients. The 1-D models provide improved accuracy when stratification is expected for a relatively low increase in the computational cost. The CFD model can be utilized for confirmatory analysis of the 1-D model, when experimental measurements are not available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2026

Ultrafast (1‐5 sec) Lamination of Perovskite Solar Cells With Self‐Encapsulation Using Rapid Joule Heating

Perovskite solar cells (PSCs) are traditionally fabricated using sequential layer‐by‐layer deposition, in which each layer of the device is processed on top of the preceding layer. This constrains the processing techniques and selection of transport layer materials that can be used in the solar cell. To overcome these challenges, two half‐cells can be processed independently and then diffusion‐bonded through a lamination process. However, current lamination processes for perovskite solar cells suffer from relatively long process times, which can limit throughput when moving toward high‐volume manufacturing. In this study, a custom platform was designed for rapid‐joule heating of perovskite materials and devices. This enabled more than a 99% reduction in lamination time from 26 min to 1 s. Perovskite samples that were laminated in 1 s exhibited comparable values of percent bonded area, interfacial toughness, grain domain size, and X‐ray diffraction spectra to those laminated in greater than 10 min. As a proof‐of‐concept, 18.3% efficient devices were successfully laminated in 5 s. A transient heat transfer model was developed to describe the relationship between the perovskite temperature and the electrical power supplied to the heaters, establishing a baseline for predicting processing conditions in large‐scale manufacturing systems. Ultra‐fast lamination provides a pathway toward scalable roll‐to‐roll or sheet‐to‐sheet manufacturing of PSCs.

heat transfer↗

DOE OSTI2024

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Local summer insolation and greenhouse gas forcing drove warming and glacier retreat in New Zealand during the Holocene

Geological climate archives from the Holocene Epoch provide baseline information concerning natural climate variability. Temperate mountain glacier extent is sensitive to summer air temperature, thus geological records of past glacier length changes are a useful proxy for this climatic variable. Here we present a new cosmogenic 10 Be chronology of glacier length changes at Dart Glacier in the Southern Alps, New Zealand. Prominent moraines deposited 321 ± 44 yr ago (n = 11) and 7.8 ± 0.3 ka (n = 5) show glaciers during the Little Ice Age were less extensive than during the early Holocene. This pattern of net Holocene glacier retreat is consistent with emerging data from other catchments in New Zealand and across the southern mid-latitudes. Using the physical framework of a transient global climate model simulation, we suggest that cool summers in the early Holocene were promoted by the local summer insolation minimum, together with low atmospheric greenhouse gas concentrations, causing an early Holocene austral glacial maximum. An insolation-driven reduction in seasonality at southern mid-latitudes may reconcile differences between early Holocene temperature reconstructions where climate proxies have different seasonal sensitivities. Furthermore, we suggest that rising greenhouse gas concentrations after 7 ka caused regional-scale glacier retreat and appear to be the dominant driver of multi-millennial summer temperature trends in the southern mid-latitudes during the present interglacial.

58 GEOSCIENCES↗

DOE OSTI2025

A Synoptic System for Capturing Ecosystem Control Points Across Terrestrial‐Aquatic Interfaces

Interconnected landscape features such as terrestrial‐aquatic interfaces play an outsized role in biogeochemical cycles as ecosystem control points, but it is notoriously challenging to characterize these. Here, we document a synoptic sensor network design that is (a) flexible to accommodate diverse ecosystem interfaces and gradients, (b) adaptable to monitoring and modeling needs of small and large projects alike, (c) standardized for intercomparability across sites and field experiments, and (d) adequately replicated to capture heterogeneity of each parameter monitored. This real‐time monitoring of surface water, groundwater, soil, and vegetation supports configuration and evaluation of models that span upland, wetland, open water strata, and transitions between them. We established the network at seven sites along the Chesapeake Bay and Lake Erie coastlines, including large‐scale flood manipulation experiments in both regions. A central design element is “one data logger program to rule them all”—a collection of sensor‐specific modules deployed on 40 loggers controlling ∼2,000 sensors, with the goal of streamlining maintenance, debugging, and reproducible data processing. The network generates ∼6 M observations per month, capturing system dynamics at the broad spatial and fine temporal scales needed to initialize and benchmark models; measurement frequency can be modified remotely to capture events. This network design has also revealed behaviors not represented in Earth system models, such as transient groundwater oxygen pulses. Completely documented and open source, this standardized, flexible, and efficient sensor network design can reduce barriers to understanding environmental changes and ecosystem responses across systems and scales.

Ward, Nicholas D. [Pacific Northwest National Labo↗

DOE OSTI2020

Spatial separation of triplet excitons drives endothermic singlet fission

Molecules that undergo singlet fission, converting singlet excitons into pairs of triplet excitons, have potential as photovoltaic materials. The possible advantages of endothermic singlet fission (enhanced use of photon energy and larger triplet energies for coupling with common absorbers) motivated us to assess the role of exciton delocalization in the activation of this process. Here we report the synthesis of a series of linear perylene oligomers that undergo endothermic singlet fission and have endothermicities in the range 5-10 kBT at room temperature in solution. We study these compounds using transient spectroscopy and modelling to unravel the singlet and triplet dynamics. We show that the minimal number of coupled chromophores needed to undergo endothermic singlet fission is three, which provides sufficient statistical space for triplet excitons to separate and avoid annihilation - and a subsequent fast return to the singlet state. Our data additionally suggest that torsional motion of chromophores about the molecular axis following triplet-pair separation contributes to the increase in entropy, thus lengthening the triplet lifetime in longer oligomers.

14 SOLAR ENERGY↗

DOE OSTI2025

A Practical Comparison of Data-Driven Prognostics Methods for Energy Systems

This study explores data-driven prognostics for nuclear power plant (NPP) condensers, focusing on tube fouling. We utilized the Asherah nuclear power plant simulator (ANS) to compare four methods: Random Forest (RF), Support Vector Regressor (SVR), Fully Connected Neural Network (FCNN), and Long Short-Term Memory Neural Network (LSTM). By simulating various fouling scenarios in the ANS, we generated data with different degradation rates under transient operations. The models were trained and tested on these data, with performance evaluated visually and numerically including uncertainty assessment. The LSTM model excelled, exhibiting minimal prediction noise and the most accurate remaining useful life estimates across all degradation levels. Its ability to capture long-term dependencies and produce cleaner outputs makes it a strong candidate, although accurate training data across the entire component lifespan are crucial. The RF model emerged as a robust alternative, providing reliable predictions with high confidence. The FCNN and SVR models, while less effective overall, showed potential under specific conditions. FCNN offers a less complex alternative to LSTM and might benefit from larger datasets. SVR excels in precision when the quality of the training data is high. Furthermore, this study highlights the operational benefits of advanced prognostics in the energy sector and emphasizes the need for further research in NPP condenser health management through real-life experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Electromagnetic Transient (EMT) Simulation Algorithm for Evaluation of Photovoltaic (PV) Generation Systems

Use of inverter-based resources facilitating renewable energy resources such as photovoltaic (PV) generation is increasing rapidly with decreasing costs and reduced emissions associated. To accommodate such rapid growth of inverter-based resources like PV systems, electromagnetic transient (EMT) simulation models of both PV systems and grids are required to analyze the interaction of PVs in the grid (like the post-event analysis). In addition, the EMT simulation would help with the planning of future power grid with a large number of PVs as well as other inverter-based distributed generation systems. In this paper, the EMT simulation models of PV systems and grids are developed based on the differential algebraic equations (DAEs) representing their EMT dynamics. Furthermore, advanced simulation algorithms including numerical stiffness-based hybrid discretization, DAE clustering and aggregation, multi-order integration, and matrix splitting approaches are applied to accelerate the EMT simulation. The proposed algorithm was applied to 125 PV inverters within 52-bus medium-voltage (MV) distribution grid.

Choi, Jongchan↗