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 829 records · Page 46

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ultrafast Relaxations in Ruthenium Polypyridyl Chromophores Determined by Stochastic Kinetics Simulations

Maximizing the efficiency of solar energy conversion using dye assemblies rests on understanding where the energy goes following absorption. Transient spectroscopies in solution are useful for this purpose, and the time-resolved data are usually analyzed with a sum of exponentials. This treatment assumes that dynamic events are well separated in time, and that the resulting exponential prefactors and phenomenological lifetimes are related directly to primary physical values. Such assumptions break down for coincident absorption, emission, and excited state relaxation that occur in transient absorption and photoluminescence of tris(2,2'-bipyridine)ruthenium(2+) derivatives, confounding the physical meaning of the reported lifetimes. Here, we use inductive modeling and stochastic chemical kinetics to develop a detailed description of the primary ultrafast photophysics in transient spectroscopies of a series of Ru dyes, as an alternative to sums of exponential analysis. Commonly invoked three-level schemes involving absorption, intersystem crossing (ISC), and slow nonradiative relaxation and incoherent emission to the ground state cannot reproduce the experimentally measured spectra. The kinetics simulations reveal that ultrafast decay from the singlet excited state manifold to the ground state competes with ISC to the triplet excited state, whose efficiency was determined to be less than unity. The populations predicted by the simulations are used to estimate the magnitudes of transition dipoles for excited state excitations and evaluate the influence of specific ligands. The mechanistic framework and methodology presented here are entirely general, applicable to other dye classes, and can be extended to include charge injection by molecules bound to semiconductor surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Power Layout Design of a GaN HEMTs-Based High-Power High-Efficiency Three-Level ANPC Inverter for 800 V DC Bus System

Multiple commutation paths exist for switching devices in a three-level active neutral point clamped (3L-ANPC) inverter operation based on the selected switching state and current direction. In addition, the capacitive coupling path of the nonswitching device is a key design aspect for enabling high voltage and high current operation of gallium nitride (GaN) switches in 3L-ANPC topology. A comprehensive study of the switching transient events of inner, outer, and clamping devices of 3L-ANPC is presented in this article. The commutation mechanisms for worst-case transient voltage overshoots (TVOs) are identified. A simplified equivalent circuit model is presented to determine the design criteria for the power layout structure's parasitic inductances. A power layout strategy satisfying the design criteria is then proposed using an insulated metal substrate power printed circuit board (PCB) to enable efficient high-power operation. The proposed design minimizes the commutation and capacitive coupling path inductances to 6 nH and 11.5 nH, respectively. This enables the fast switching operation of GaN HEMTs at 800 V dc, 36 A with a low TVO of 31% verified through experimental three-level double pulse test results. Experimental evaluation of a three-phase 3L-ANPC hardware prototype based on the proposed power layout shows 99% efficiency at 800 V, 9.5 kVA and 50 kHz switching frequency. The proposed design achieves a low case-to-ambient thermal resistance of 2.3 °C/W.

42 ENGINEERING↗

Increased Ly α Opacity in White Dwarf Photospheres from Transient H − Resonances

The structure and emergent flux of hydrogen atmosphere white dwarfs depend on the opacity of the Lyα and Lyβ spectral lines. The opacity here is set by the strength and broadening of these lines; the latter is dictated by the far line wing, which is in the “quasi-static” limit of electron broadening, placing it in the incomplete collision regime, and describes the transient parts of electron and ion collisions. These transient stages of the collision form resonances: In the case of ions, they manifest as molecular resonances, while for electrons they are H − resonances, both of which can only be captured quantum-mechanically. Quantum-mechanical calculations have historically preserved only a handful of broadening terms that are most important near the center of the line. However, in the wings of the line, the previously neglected terms that describe the transient stages of the collision need to be included. This requirement arises because, in the line wings, the broadening from the 1s ground state, which is generally assumed to be extremely small compared to the broadening of the upper state, is no longer negligible within a quantum-mechanical model that takes into account exchange interactions. The inclusion of all the transient terms results in asymmetries and extra broadening. The increased broadening of Lyα increases the opacity at the energy where most of the flux leaves the star. The broader Lyα lines also impact the visible flux, raising it by an amount that exceeds previously estimated errors.

74 ATOMIC AND MOLECULAR PHYSICS↗

A novel open-source cultivation system helps establish the first full cycle chemosynthetic symbiosis model system involving the giant ciliate Zoothamnium niveum

Symbiotic interactions drive species evolution, with nutritional symbioses playing vital roles across ecosystems. Chemosynthetic symbioses are globally distributed and ecologically significant, yet the lack of model systems has hindered research progress. The giant ciliate Zoothamnium niveum and its sulfur-oxidizing symbionts represent the only known chemosynthetic symbiosis with a short life span that has been transiently cultivated in the laboratory. While it is experimentally tractable and presents a promising model system, it currently lacks an open-source, simple, and standardized cultivation setup. Following the FABricated Ecosystems (EcoFABs) model, we leveraged 3D printing and polydimethylsiloxane (PDMS) casting to develop simple flow-through cultivation chambers that can be produced and adopted by any laboratory. The streamlined manufacturing process reduces production time by 86% and cuts cost by tenfold compared to the previous system. Benchmarking using previously established optimal growth conditions, the new open-source cultivation system proves stable, efficient, more autonomous, and promotes a more prolific growth of the symbiosis. For the first time, starting from single cells, we successfully cultivated the symbiosis in flow-through chambers for 20 days, spanning multiple generations of colonies that remained symbiotic. They were transferred from chamber to chamber enabling long-term cultivation and eliminating the need for continuous field sampling. The chambers, optimized for live imaging, allowed detailed observation of the synchronized growth between the host and symbiont. Highlighting the benefit of this new system, we here describe a new step in the first hours of development where the host pauses growth, expels a coat, before resuming growth, hinting at a putative symbiont selection mechanism early in the colony life cycle. With this simple, open-source, cultivation setup, Z. niveum holds promises for comparative studies, standardization of research and wide adoption by the symbiosis research community.

59 BASIC BIOLOGICAL SCIENCES↗

Impact of the Exciter and Governor Parameters on Forced Oscillations

In recent years, the frequency of forced oscillation events due to control system malfunctions or improper parameter settings has increased. Tuning the parameters of exciters and governor models is crucial for maintaining power system stability. Traditional simulation studies typically involve small transient disturbances or step changes to find optimal parameter sets, but existing optimization algorithms often fall short in fine-tuning for forced oscillations. Identifying the sensitive parameters within these control models is essential for ensuring stability during large, sustained disturbances. This study focuses on identifying these critical exciter and governor model parameters by analyzing their influence on sustained forced oscillations. Using Kundur’s two-area system, we analyze common exciter models such as SCRX, ESST1A, and AC7B, along with governor models like GAST, HYGOV, and GGOV1, utilizing PSS®E software version 34. Sustained forced oscillations are injected at generator-1 of area-1, with individual parameter changes dynamically simulated. By considering a local oscillation frequency of 1.4 Hz and an inter-area oscillation mode of 0.25 Hz, we analyze the impact of each parameter change on the magnitude and frequency of forced oscillations as well as on active and reactive power outputs. This novel approach highlights the most influential parameters of each tested model—such as exciter, governor, and turbine gains, as well as time constant parameters—on the impact of forced oscillations. Based on our findings, the sensitive parameters of each tested model are ranked. These would provide valuable insights for industry operators to fine-tune control settings during oscillation events, ultimately enhancing system stability.

42 ENGINEERING↗

Multichannel deconvolution of vibrational signals: A state-space inverse filtering approach

Deconvolution of noisy measurements, especially when they are multichannel, has always been a challenging problem. The processing techniques developed range from simple Fourier methods to more sophisticated model-based parametric methodologies based on the underlying acoustics of the problem at hand. Methods relying on multichannel mean-squared error processors (Wiener filters) have evolved over long periods from the seminal efforts in seismic processing. However, when more is known about the acoustics, then model-based state-space techniques incorporating the underlying process physics can improve the processing significantly. The problems of interest are the vibrational response of tightly coupled acoustic test objects excited by an out-of-the-ordinary transient, potentially impairing their operational performance. Further, employing a multiple input/multiple output structural model of the test objects under investigation enables the development of an inverse filter by applying subspace identification techniques during initial calibration measurements. Feasibility applications based on a mass transport experiment and test object calibration test demonstrate the ability of the processor to extract the excitations successfully.

42 ENGINEERING↗

Intracore Natural Circulation Study in the High Temperature Test Facility

The development of the Modular High-Temperature Gas-Cooled Reactor is a significant milestone in advanced nuclear reactor technology. One of the concerns for the reactor’s safe operation is the effects of a loss-of-flow accident (LOFA) where the coolant circulators are tripped, and forced coolant flow through the core is lost. Depending on the steam generator placement, loop or intracore natural circulation develops to help transfer heat from the core to the reactor cavity, cooling system. This paper investigates the fundamental physical phenomena associated with intracore coolant natural circulation flow in a one-sixth Computational Fluid Dynamics (CFD) model of the Oregon State University High Temperature Test Facility (OSU HTTF) following a loss-of-flow accident transient. This study employs conjugate heat transfer and steady-state flow along with an SST k-ω turbulence model to characterize the phenomenon of core channel-to-channel natural convection. Previous studies have revealed the importance of complex flow distribution in the inlet and outlet plenums with the potential to generate hot coolant jets. For this reason, complete upper and lower plenum volumes are included in the analyzed computational domain. CFD results also include parametric studies performed for a mesh sensitivity analysis, generated using the STAR-CCM+ software. The resulting channel axial velocities and flow directions support the test facility scaling analysis and similarity group distortions calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Sensitivity analysis of in-pile critical heat flux experiments in TREAT for characterization of RIA power-transient effects

A reactivity-initiated accident (RIA) is one type of postulated design basis accident (DBA) that can cause a departure from nucleate boiling (DNB) event in pressurized-water reactors (PWRs). A DNB occurrence and its consequences depend on the thermophysical properties of the fuel components and coolant, characteristics of the transient energy insertion into the fuel rod, and the onset of the critical heat flux (CHF) phenomenon. To leverage the restart of the Transient Reactor Test (TREAT) Facility, an effort is currently underway to better understand the cladding-to-coolant heat transfer mechanisms and the CHF phenomenon under fast-transient irradiation conditions. This paper characterizes the impact of power transients on the thermal-hydraulic behavior of a TREAT Facility reactor heater rodlet CHF experiment to provide the priority of parameters that need to be investigated for an improved CHF model. Sobol sensitivity analysis methods and the Reactor Excursion and Leak Analysis Program (RELAP5-3D) code were used to identify key input parameters on the uncertainty in the prediction of peak outer- and inner-surface temperatures of the heater tube, as well as the time of the DNB event. A series of sensitivity analyses revealed the total energy deposition on the tube and the transient effects of power pulse had large impacts on the maximum temperatures. The CHF multiplier had the largest impact on the time occurrence of CHF. The overall results show the energy deposition rate in the tube is the most influencing factor to the manifestation of CHF and the resulting thermal-hydraulic behaviors of the tube. The multiplier for the CHF, which is interpreted as the predicted CHF value, has the largest Sobol indices for the time of the CHF in all cases, since it directly determines the occurrence of CHF. It is inferred that the uncertainties in the thermal-hydraulic behaviors of fuels increase with respect to the key parameters as the power pulse becomes broader, and an accurate estimation of the energy deposition rate is required to reduce the uncertainty in the evaluation of the integrity of fuel if the CHF is expected to occur near the peak power. Therefore, the outputs are expected to provide rigorous interpretation of ongoing in-pile CHF experiments in the TREAT Facility reactor regarding thermal-hydraulic behavior of the fuel system aiming for a new transient in-pile CHF model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Atomic-scale polarization switching in wurtzite ferroelectrics

Ferroelectric wurtzites have the potential to revolutionize modern microelectronics because they are easily integrated with multiple mainstream semiconductor platforms. However, the electric fields required to reverse their polarization direction and unlock electronic and optical functions need substantial reduction for operational compatibility with complementary metal-oxide semiconductor (CMOS) electronics. To understand this process, we observed and quantified real-time polarization switching of a representative ferroelectric wurtzite (Al 0.94 B 0.06 N) at the atomic scale with scanning transmission electron microscopy. The analysis revealed a polarization reversal model in which puckered aluminum/boron nitride rings in the wurtzite basal planes gradually flatten and adopt a transient nonpolar geometry. Independent first-principles simulations reveal the details and energetics of the reversal process through an antipolar phase. Furthermore, this model and local mechanistic understanding are a critical initial step for property engineering efforts in this emerging material class.

36 MATERIALS SCIENCE↗

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↗

Aggregation of Leslie matrix models with application to ten diverse animal species

Some grouping is necessary when constructing a Leslie matrix model because it involves discretizing a continuous process of births and deaths. Here, the level of grouping is determined by the number of age classes and frequency of sampling. It is largely unknown what is lost or gained by using fewer age classes, and I address this question using aggregation theory. I derive an aggregator for a Leslie matrix model using weighted least squares, determine what properties an aggregated matrix inherits from the original matrix, evaluate aggregation error, and measure the influence of aggregation on asymptotic and transient behaviors. To gauge transient dynamics, I employ reactivity of the standardized Leslie matrix. I apply the aggregator to 10 Leslie models developed for animal populations drawn from a diverse set of species. Several properties are inherited by the aggregated matrix: (a) it is a Leslie matrix; (b) it is irreducible whenever the original matrix is irreducible; (c) it is primitive whenever the original matrix is primitive; and (d) its stable population growth rate and stable age distribution are consistent with those of the original matrix if the least squares weights are equal to the original stable age distribution. In the application, depending on the population modeled, when the least squares weights do not follow the stable age distribution, the stable population growth rate of the aggregated matrix may or may not be approximately consistent with that of the original matrix. Transient behavior is lost with high aggregation.

54 ENVIRONMENTAL SCIENCES↗

A Novel Physiologically Based Model of Creatinine Renal Disposition to Integrate Current Knowledge of Systems Parameters and Clinical Observations

Creatinine is the most common clinical biomarker of renal function. As a substrate for renal transporters, its secretion is susceptible to inhibition by drugs, resulting in transient increase in serum creatinine and false impression of damage to kidney. Novel physiologically based models for creatinine were developed here and (dis)qualified in a stepwise manner until consistency with clinical data. Data from a matrix of studies were integrated, including systems data (common to all models), proteomics‐informed in vitro – in vivo extrapolation of all relevant transporter clearances, exogenous administration of creatinine (to estimate endogenous synthesis rate), and inhibition of different renal transporters (11 perpetrator drugs considered for qualification during creatinine model development and verification on independent data sets). The proteomics‐informed bottom‐up approach resulted in the underprediction of creatinine renal secretion. Subsequently, creatinine‐trimethoprim clinical data were used to inform key model parameters in a reverse translation manner, highlighting best practices and challenges for middle‐out optimization of mechanistic models.

Scotcher, Daniel↗

Ring Model Development and Validation for Prismatic HTGR Core Thermal-hydraulics and Safety Analysis

Because of the complex core geometry, prismatic high temperature gas-cooled reactors (prismatic HTGRs) often exhibit complex thermal fluid behaviors during both normal operating and transient conditions. Most HTGR designs rely on passive safety system for decay heat removal, such as the reactor cavity cooling system (RCCS). During postulated accidents like Pressurized Conduction Cooldown (PCC) event, the decay heat is first radially transferred from the core region to the reactor vessel outer surface, then to the RCCS cooling panels. The peak fuel temperature is controlled by heat transfer mechanisms with two distinctive characteristic length scales, i.e., the core-wise effective heat conduction and local heat conduction in the fuel pellet scale. As both length scales are essential to determine the fuel temperature, from the modeling perspective, computer codes must be able to capture heat transfer in both scales. This is challenging for both computational fluid dynamics (CFD) tools because of extremely large amount of computation resources required, and for system analysis codes because of the challenge to model the complex core geometry. Under the support of DOE-NE’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, efforts have been pursued to support HTGR technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior during safety-related transients. In our previous studies, the ring model was adopted in SAM and further developed to simulate the normal operating condition and a PCC event using the MHTGR-350 design of General Atomics as the reference design. This current work represents a continuation of these previous efforts, and the focus is to critically review and examine simplifications and assumptions made to develop the ring model, and to perform code validation using experimental data from an integral-effect test facility, the High Temperature Test Facility (HTTF) at the Oregon State University. In this study, the test PG-27 from the HTTF test suite was selected for code validation purpose. The test PG-27 is a transient test designed to simulate the PCC event of the MHTGR design. Very good agreements between SAM prediction and experimental measurements were found in both coolant and solid structure temperatures during the transient.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Fast and Accurate Transient Stability Assessment Method Based on Deep Learning: WECC Case Study

Transient stability is one of the critical aspects of power system stability assessment. The increasing integration of inverter-based resources and the retirement of conventional synchronous generators result in the decreasing system inertia and growing complexity of system operating conditions. Using a few selected typical operating conditions cannot guarantee system transient stability in all operating conditions, and the time-domain simulation of all operating conditions requires tremendous time and is often infeasible. This paper proposes a more efficient transient stability assessment method based on deep learning. The binary search method is used to determine the critical clearing time (CCT) in creating training databased by time-domain simulation. This method is fast and accurate with 1 ms resolution. The buses whose CCTs are lower than 200 ms are considered critical buses. Buses close to each other are grouped based on their mutual admittance matrix to reduce the search space of the critical buses. This paper also proposes the generator feature normalization based on the physical model. Case study on the reduced 240-bus WECC system model demonstrates that the proposed method can predict CCT accurately and efficiently.

critical clearing time↗

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗