Search NASA⌕ Search

SEARCH · Search NASA

Results for “engine 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,441 records · Page 80

Multi-scale Simulation, Calibration, and Optimization of Calcium Carbonate Precipitation in Microbial Communities

Ensuring the efficient engineering of microbially induced calcium carbonate precipitation (MICP) is crucial for a variety of environmental and civil engineering applications, such as soil stabilization and carbon sequestration. Addressing this need, we present a comprehensive multi-scale workflow that begins with the isolation of calcium carbonate-producing microbes from soil samples, followed by metagenomic sequencing and metabolic reconstruction. We then characterize microbial growth phenotypes under diverse nutrient conditions, compare observed growth with metabolic model predictions, and apply the Consistent Reproduction of Phenotype (CROP) algorithm to refine these models. Furthermore, we analyze metabolite consumption and production, and develop a consumer-resource model that is calibrated using time-series measurements of growth rates, pH levels, and calcium carbonate precipitation. The primary benefit of our approach lies in its ability to predict and control MICP outcomes, facilitated by a Bayesian methodology that incorporates priors on initial conditions and parameters. This allows us to compute posteriors by integrating experimental data, and to solve a risk optimization problem under uncertainty to identify nutrient conditions that maximize calcium carbonate production. In contrast to non-Bayesian methods, which fail to quantify uncertainty accurately, our approach provides a more reliable pathway to optimizing nutrient conditions, enhancing the likelihood of achieving desired MICP outcomes. This positions our method as a superior alternative in the quest to improve MICP through engineered microbial consortia.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Optimal Net Load Management in Microgrids Using Hardware-in-the-Loop Simulation: Preprint

This paper presents the performance evaluation of a net load management (NLM) engine that balances load and generation in an isolated community to power a critical facility after a grid interruption event (e.g., lost a large generation unit). The NLM engine is a central dispatch control system that provides high-speed cost-optimal coordination of a set of net loads (combination of generation and deferrable loads) connected to the same electrical network. The NLM algorithm cost-optimally dispatches the grid-following inverters and deferrable loads based on the demanded power and load priorities, and the grid-forming inverters use droop control to form system voltages and share active and reactive power. A controller-hardware-in-the-loop platform is developed to evaluate the control performance of the NLM algorithm. The experimental results indicate that the NLM engine can achieve the targeted system voltage and frequency and balance load and generation to serve the critical facility.

droop control↗

PCI analysis of a commercial PWR using BISON fuel performance code

Through the early 2000’s, nuclear fuel cladding failure resulting from pellet-cladding interaction in commercial pressurized water reactor fuel was managed successfully by incorporating power ramp restrictions. Within the past decade, however, cladding failures have again begun to appear, often correlated to missing pellet surface defects, and further aggravated by modern reactor operating strategies demanding higher fuel duty operations. In response, renewed efforts in fuel behavior modeling and simulation are being pursued to improve our understanding of the conditions that can lead to missing pellet surface related cladding failure. This paper describes a methodology used to model pellet-cladding interaction related failures observed in commercial pressurized water reactors using the BISON fuel performance code developed at Idaho National Laboratory. In order to accomplish this goal, a three-step process is followed: 1) A failure-stress threshold is developed by using BISON code to model ramp test with known pellet-clad interaction related failures, 2) BISON analysis is applied to a commercial pressurized water reactor cycle startup with one rod confirmed to contain a missing pellet surface defect, and 3) an missing pellet surface sensitivity study of a cycle startup of the same commercial pressurized water reactor which contained a failed fuel rod of un-identified type. The analysis results described in this paper serve to illustrate the BISON code capabilities for modeling pressurized water reactor fuel subjected to higher duty power events.

Nuclear Science & Technology↗

Ferroelectric Fractals: Switching Mechanism of Wurtzite AlN

The advent of wurtzite ferroelectrics is enabling new ferroelectric devices for computer memory that have the potential to bypass the von Neumann bottleneck due to their robust polarization and silicon compatibility. However, the atomistic switching mechanism of wurtzites is still undetermined due to the limitations of density functional theory simulation size and experimental temporal and spatial resolution. Thus, physics-informed materials engineering to reduce coercive field and breakdown in these devices has been limited. In this work, the atomistic mechanism of domain wall migration and domain growth in aluminum nitride-based wurtzites is uncovered using molecular dynamics and Monte Carlo simulations. We reveal the anomalous switching mechanism of fast 1D single columns of atoms propagating from a slow-moving 2D fractallike domain wall. We find that the critical nucleus is a single aluminum ion that breaks its bond with one nitrogen and bonds to another nitrogen; this creates a cascade that flips atoms directly only in the same column, due to the extreme locality (sharpness) of the domain walls in wurtzites. We further show how the fractallike shape of the domain wall in the 2D plane breaks assumptions in the Kolmogorov, Avrami, and Ishibashi (KAI) model and leads to the anomalously fast switching in wurtzite structured ferroelectrics.

36 MATERIALS SCIENCE↗

Large-Scale Simulation of Regional Demand Flexibility Implementation and Customer Economic Impact

The Distribution System Operator with Transactive (DSO+T) study conducted a large-scale simulation of over 60,000 customers in a region the size of Texas to demonstrate the effective coordination of distributed energy resources (DERs) in commercial and residential buildings. The integrated simulation included both the bulk (wholesale generation and transmission) and distribution systems. The DERs (including batteries, electric vehicles, air conditioning, and water heaters) participated in a transactive energy retail market that was integrated into an existing competitive wholesale market. The engineering and economic performance of the resulting demand flexibility was evaluated over annual simulations for both moderate and high renewable generation scenarios. A detailed parametric cost model was developed to enable detailed economic analysis of key stakeholders. In addition, fixed and dynamic customer tariffs were designed and applied to the customer population. This allowed the impact on annual customer bills to be analyzed for various building types (residential versus commercial; single- versus multi-family). This paper presents results showing the relative flexibility of batteries, electric vehicles, and building loads throughout the year and under different renewable scenarios. This feeds a detailed breakdown of the impact this flexibility has on the operating costs of the grid and the resulting net economic benefit. Finally, the study showed that practically all customer classes (including non-participating customers) save money under the proposed demand flexibility scheme. The study found overall net annual economic savings of $3.3-5.0B for a region the size of Texas equating to average customer bill savings of 10-16%.

Reeve, Hayden M.↗

Modeling and Simulation of the Off-gas in an Electric Arc Furnace

The following paper describes an approach to process modeling and simulation of the gas phase in an electric arc furnace (EAF). The work presented represents the continuation of research by Logar, Dovžan, and Škrjanc on modeling the heat and mass transfer and the thermochemistry in an EAF. Due to the lack of off-gas measurements, Logar et al. modeled a simplified gas phase under consideration of five gas components and simplified chemical reactions. The off-gas is one of the main continuously measurable EAF process values and the off-gas flow represents a heat loss up to 30 pct of the entire EAF energy input. Therefore, gas phase modeling offers further development opportunities for future EAF optimization. This paper presents the enhancement of the previous EAF gas phase modeling by the consideration of additional gas components and a more detailed heat and mass transfer modeling. In order to avoid the increase of simulation time due to more complex modeling, the EAF model has been newly implemented to use an efficient numerical solver for ordinary differential equations. Compared to the original model, the chemical components H{sub 2}, H{sub 2}O, and CH{sub 4} are included in the gas phase and equilibrium reactions are implemented. The results show high levels of similarity between the measured operational data from an industrial scale EAF and the theoretical data from the simulation within a reasonable simulation time. In the future, the dynamic EAF model will be applicable for on- and offline optimizations, e.g., to analyze alternative input materials and mode of operations.

42 ENGINEERING↗

Computational Study of Droplet Dynamics in Dispersed Flow Film Boiling Regime in PWR Geometries

Dispersed flow film boiling (DFFB) regime is exhibited during low-flooding rate conditions following the postulated Loss-of-Coolant Accident (LOCA) in nuclear reactor flow channels. LOCA is a type of design-basis accident scenario, which is crucial for reactor design and licensing. Accurate models to predict the heat transfer in DFFB regime is critical for the reactor to maintain safety limits. The nominal operating linear heat rate is also set based on the maximum linear heat rate determined for a LOCA transient. It is understood that the behavior of entrained droplets in DFFB regime greatly affects the heat transfer by evaporating and desuperheating the surrounding steam. Furthermore, interaction with flow structures also aids in increasing the evaporation rate of entrained droplets. In order to quantify this enhancement in heat transfer, the size of droplets downstream of flow structures (Spacer Grids and Mixing Vanes) is of crucial importance. Thermal hydraulic analysis codes like COBRA-TF, CATHARE-3 and adopt a three-field equation model that allows them to model the droplet field separately. However, to ensure the accuracy of these field equations, various correlations and physical models are required. This paper proposes a method for characterizing the mean Sauter Diameter of droplets across spacer-grids and mixing vanes in PWR sub-channels using PHASTA - a high fidelity flow solver with multiphase capabilities using the level-set method. The PHASTA code is a finite element based computational fluid dynamics code designed for solving compressible and incompressible turbulent flows. It has been demonstrated to be massively scalable (up to 786,432 cores), which is crucial for simulating large realistic domains. The code incorporates the level-set method for interface tracking of the gas-liquid interface.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A statistical analysis of developing knock intensity in a mixture with temperature inhomogeneities

Knock formation and its intensity for a stoichiometric ethanol/air mixture under a representative end-gas auto-ignition condition in IC engines with temperature inhomogeneities are investigated using multi-dimensional direct numerical simulations (DNS) with a 40-species skeletal mechanism of ethanol. Two- and three-dimensional simulations are performed by systematically varying temperature fluctuations and its most energetic length scale, l T . The volumetric fraction of the mixture regions that have the propensity to detonation development, F D , is proposed as a metric to predict the amplitude of knock intensity. It is found that with increasing l T , F D shows a good agreement with the heat release fraction of the mixture regions with pressure greater than equilibrium pressure, F H . The detonation peninsula is well captured by F D and F H when plotting them as a function of the volume-averaged ξ, ξ, (ξ= a/S sp is the ratio of the acoustic speed, a to the ignition front speed, S sp ) Decreasing l T is found to significantly reduce the super-knock intensity. The results suggest that decreasing l T , as in engines with tumble designs resulting in a smaller turbulence scale, will be effective in mitigating the the super-knock development.

42 ENGINEERING↗

Performance Characterization of a Natural Gas-Air Rotating Detonation Engine at Elevated Pressure

Performance of a rotating detonation engine was studied to understand the challenges of achieving and quantifying a net pressure gain at engine-relevant conditions. The combustor under study was designed for operation with natural gas and air at elevated chamber pressures and air preheat temperatures representative of land based power generation gas turbine engines. The thrust produced by the combustor was measured to characterize its performance while high-frequency pressure transducers and high speed imaging of broadband combustion chemiluminescence provided information about the wave structure and dynamics. The rotating detonation engine was tested with two injector configurations to explore the effect of injector design and combustor operating mode on device performance. Analysis using multiple metrics showed nearly identical performance between similar operating conditions that alternately produced a single detonation wave or two counter-propagating waves in the combustion chamber. Across the range of tested conditions, a maximum pressure gain of -42% was measured. Several tests were observed to produce up to 10% more thrust than an ideal, constant pressure combustor, though this comparison ignores injection losses. Comparison of measured pressure gain with idealized simulations showed that the net pressure gain of a rotating detonation engine is closely tied to the area ratio between the inlet and outlet of the combustor.

20 FOSSIL-FUELED POWER PLANTS↗

The Influence of the Instantaneous Collapse of Tailings Pond on Downstream Facilities

To evaluate the evolutionary processes guiding the formation of the tailings‐water mixtures produced by the instantaneous collapse of tailings ponds and the influence of these on downstream facilities, a 2D simulation model with reasonable boundary and working conditions derived from actual engineering practice was built in this study, and the relationship between dam‐break elevation and impact on downstream facilities was also analyzed to determine the relevant mechanism of influence. Computational results indicated that lowering the dam‐break elevation caused the maximum velocity and flooding depth, along with the flooded area at monitoring points, to gradually increase. The occurrence times of maximum velocity and flooding depth were also gradually moved forward as the breaking elevation was reduced; this effect is directly related to the increase in the total potential energy at the lower break elevations. Further simulations of sand‐prevent dams with different heights located downstream from a tailings pond were carried out to identify methods for mitigating the impact of dam failure. The results revealed that increasing the height of the sand‐prevent dam reduced the production of tailings mixtures. Based on the results, the construction of a sand‐prevent dam with a crest elevation equal to that of the starter dam was recommended.

Wang, Sha (ORCID:0000000324875445)↗

Computational Simulation of an Additively Manufactured Marine Component

Early work in finite-element analysis of additive manufacturing (AM) components was not able to include distortion and residual stress of the entire build, yielding uncertainty in the computational model predictions. To address this, current versions of recently developed computational tools are being used to understand their role in validation at the Naval Surface Warfare Center, Carderock Division, for laser powder-bed fusion (L-PBF) AM. A 316L stainless steel L-PBF build containing multiple parts, which exhibited the same local part failure over multiple AM fabrication runs, was selected for computational analysis. Single-part simulations analyzed optimization of computational speed and prediction accuracy. Full-build simulations exposed the smallest re-coater tolerances early in part fabrication as a result of part shrinkage and subsequent cupping. This cupping was predicted to change the powder spreading and result in the non-consolidation of the support structures, which paralleled the failures seen in the physical builds. Part redesign of the conventionally produced part to align with AM fabrication requirements was recommended to alleviate the build failures. This project is part of a larger program focused on validation of the computational tools to determine best practices for insertion of AM simulation into the greater Naval enterprise.

36 MATERIALS SCIENCE↗

Electromagnetic Transient (EMT) Simulation Algorithms for Evaluation of Large-Scale Extreme Fast Charging Systems (T&D Models)

Simulation of high-fidelity models of extreme fast charging (XFC) systems and large-area power grids with many XFCs can be time consuming in traditional simulators. Traditional simulators use a single method of discretization for all the components that results in imposing a large computational burden of inverting a large matrix as well as increased computations related to single method of discretization (that is typically a trapezoidal method). To overcome the problem of simulating large-area power grids with many XFCs, in this paper, advanced numerical simulation algorithms are applied for the first time together to reduce the dimension of matrix inversion. Here, the algorithms include numerical stiffness-based segregation, time constant-based segregation, clustering and aggregation on differential algebraic equations (DAEs), and multi-order integration approaches. These algorithms apply multiple discretization algorithms rather than a single discretization algorithm that further reduces the computational burden. The approaches mentioned here have resulted in speed-up of up to 18x in the simulation of a single distribution system with 15 XFCs and of up to 271x in the simulation of a transmission-distribution system with 300 XFCs in multiple distribution feeders with respect to conventional simulators (like power systems computer aided design [PSCAD]).

42 ENGINEERING↗

Retrieval and Management of TRU Sludge at Hanford K-Basin - 19641

At the US Department of Energy Hanford site, the K-East and K-West reactor fuel storage basins (KE or KW Basin) were used in the 1970's and 1980's to store more than 100,000 uranium fuel rods and rod fragments that had been irradiated in the Hanford N Reactor. Approximately 2,100 metric tons of spent nuclear fuel that had been stored in the KW Basin were removed in 2004. Over time a mixture of small fuel corrosion particles, fuel rod metal fragments, and other debris accumulated in the KW Basin from the degradation of the spent nuclear fuel during storage. This is the source of the TRU sludge stored in the KW Basin. The corrosion products mixed with windblown sand, concrete dust, and other debris in the basin to form what is known as sludge. In 2006/2007 KE Basin sludge of similar origin as the KW Basin was transferred to KW Basin. Approximately 27 m{sup 3} of highly radioactive sludge was collected and consolidated underwater in six specially designed engineered containers. The mission of the Hanford Sludge Treatment Project (STP) is to develop and deploy technologies to safely retrieve, transport, and interim store the approximately 27 m{sup 3} of remote-handled transuranic (RH-TRU) sludge that is stored in the KW Basin. After extensive technology development, testing (with a bounding range of simulants), design, and procurement a sludge removal and storage system called the Engineered Container Retrieval and Transfer System (ECRTS) completed acceptance testing in 2018. The ECRTS system is designed to remotely retrieve and transport high density, highly erosive, uranium metal pieces as well as small low-density particles. The ECRTS system includes a retrieval tool based on the XAGO{sup tm} design from the United Kingdom, a heavily modified positive displacement pump, a hose-in-hose transfer system, a Sludge Transport and Storage Container (STSC), a settling and decant system, and a shielded transport system. The first sludge transfer was successfully completed on June 12, 2018 proving that the ECRTS process equipment worked as intended on the actual sludge. Three batches (i.e. transfer and decant) were required to obtain a ∼ 100% fill of the first STSC. The STSC was then transported to the on-site storage facility. Over 20 STSCs will be filled with RH-TRU sludge over approximately 16 months. This paper will discuss the issues and successes associated with this project including the operational aspects of this new system. Experience gained from this successful project should assist anyone with the need to retrieve and manage remote handled sludge from fuel storage basins. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Accurate machine-learning predictions of coercivity in high-performance permanent magnets

Increased demand for high-performance permanent magnets in the electric vehicle and wind-turbine industries has prompted the search for cost-effective alternatives. Discovering magnetic materials with the desired intrinsic and extrinsic permanent magnet properties presents a significant challenge to researchers because of issues with the global supply of rare-earth elements, material stability, and a low maximum magnetic energy product BH max . While first-principles density functional theory (DFT) predicts materials’ magnetic moments, magnetocrystalline anisotropy constants, and exchange interactions, it cannot compute extrinsic properties such as coercivity (H c ). Although it is possible to calculate H c theoretically with micromagnetic simulations, the predicted value is larger than the experiment by almost an order of magnitude due to the Brown paradox. To circumvent these issues, we employ machine-learning (ML) methods on an extensive database obtained from experiments, DFT calculations, and micromagnetic modeling. The use of a large experimental dataset enables realistic H c predictions for materials such as Ce-doped Nd 2 ⁢Fe 14 ⁢B, comparing favorably against micromagnetically simulated coercivities. Remarkably, our ML model accurately identifies uniaxial magneto-crystalline anisotropy as the primary contributor to H c . With DFT calculations, we predict the Nd-site-dependent magnetic anisotropy behavior in Nd 2 ⁢Fe 14 ⁢B, confirming that Nd 4⁢g sites mainly contribute to uniaxial magnetocrystalline anisotropy, and also calculate the Curie temperature (T c ). Finally, both calculated results are in good agreement with the experiments. The coupled experimental dataset and ML modeling with DFT input predict H c with far greater accuracy and speed than was previously possible using micromagnetic modeling. Further, we reverse engineer the grain-boundary and intergrain exchange coupling with micromagnetic simulations by employing the ML predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Neutron Scattering Investigation of the Relationship between Molecular Structure, Morphology and Dynamics in Conjugated Polymers

Conjugated polymers (CP) have enabled significant technological advances such as inexpensive organic solar cells (OPVs), biomedical sensor devices and flexible displays incorporating organic light emitting diodes (OLEDs) among others. Some of the distinguishing properties of CPs include simple processing from solutions, synthesis from abundant raw materials (C, H, S, N), flexible device form-factors and a capacity to fine-tune optical and electronic properties. Moreover, conductive polymers can assemble into large organized structures, such as nanowires and networks, that can significantly boost their performance. However, it is observed that small changes in the structure of the molecules (e.g. addition of alkyl groups) results in drastic and unpredictable changes in organized structures and electronic properties. This project utilized a comprehensive suite of elastic (SANS and USANS) and quasi-elastic (NSE and backscattering spectroscopy) neutron scattering techniques along with simultaneous property evaluation to unravel the relationships between molecular structure, self-assembly, directed-assembly and the resulting properties. Completed project objectives were: 1) to evaluate the influence of chain architecture on solution conformation, 2) to determine of the influence of molecular structure on CP self and directed assembly and 3) to study the correlation between dynamic molecular relaxations and macroscopic properties of CPs in the solid state. The experimental components of the project were further complemented by advanced molecular simulation strategies (e.g. atomistic MD and DFT) to fully understand the relationship between chain structure, mesoscale morphology, dynamic relaxation processes and the macroscopic properties.

36 MATERIALS SCIENCE↗

Generalized 3D beam dynamics model for industrial traveling wave linacs design and simulations

Beam dynamics simulations in traveling wave (TW) accelerating structures with significant beam loading is a challenging problem. Some codes are capable of calculating the TW electromagnetic fields, and some can track particles through such fields, but most cannot treat both self-consistently. A few commercial codes can model the physics correctly at the expense of many processor-hour computations to obtain converged self-consistent solutions. However, simple, accurate equations of motion for intensive beam dynamics in TW accelerating structure analysis have been previously obtained by Masunov and subsequently implemented in the Hellweg code, which was developed at the Moscow Engineering-Physics Institute. Hellweg is based on equations that allow fast simulations of beam dynamics while taking into account such effects as beam loading, space charge, and external magnetic fields. In this paper, we describe in detail some recent improvements to the Hellweg physics kernel. We have generalized the Masunov results into a 3D set of equations of motion, which include all spatial components of the radiofrequency (RF) and external magnetic fields. We have also improved the Lapostolle space charge model to the general 3D ellipsoid form for any dimensions’ ratio, consideration for the particles outside the beam core, and the fields from the neighboring bunches that can exist in the real machine. These modifications allow approaching the Hellweg accuracy to the self-consistent commercial codes while keeping the simulation time short, which is essential during the linac design and optimization stage. The implementation of these new capabilities in Hellweg is carefully benchmarked against other codes and analytical calculations. The code is freely available with an open source license.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Structure of disordered TiO 2 phases from ab initio based deep neural network simulations

Amorphous TiO 2 (a-TiO 2 ) is widely used in many fields, ranging from photo-electrochemistry to bio-engineering, hence detailed knowledge of its atomic structure is of scientific and technological interest. Here we use an ab initio-based deep neural network potential (DP) to simulate large scale atomic models of crystalline and disordered TiO 2 with molecular dynamics. Our DP reproduces the structural properties of all (11) TiO 2 crystalline phases, predicts the densities and structure factors of molten and amorphous TiO 2 with only a few percent deviation from experiments, and describes the pressure dependence of the amorphous structure in agreement with recent observations. Furthermore, it can be extended to model additional structures and compositions and can be thus of great value in the study of TiO 2 -based (nano-)materials.

36 MATERIALS SCIENCE↗

PARALLEL PERTURBATION MODEL FOR CYCLE TO CYCLE VARIABILITY PPM4CCV

This code consists of a Fortran 90 implementation of the parallel perturbation model to compute cyclic variability in spark ignition (SI) engines. Cycle-to-cycle variability (CCV) is known to be detrimental to SI engine operation resulting in partial burn and knock, and result in an overall reduction in the reliability of the engine. Numerical prediction of cycle-to-cycle variability (CCV) in SI engines is extremely challenging for two key reasons: (i) high-fidelity methods such as large eddy simulation (LES) are required to accurately capture the in-cylinder turbulent flow field, and (ii) CCV is experienced over long timescales and hence the simulations need to be performed for hundreds of consecutive cycles. In the new technique, the strategy is to perform multiple parallel simulations, each of which encompasses 2-3 cycles, by effectively perturbing the simulation parameters such as the initial and boundary conditions. The PPM4CCV code is a pre-processing code and can be coupled with any engine CFD code. PPM4CCV was coupled with Converge CFD code and a 10-time speedup was demonstrated over the conventional multi-cycle LES in predicting the CCV for a motored engine. Recently, the model is also being applied to fired engines including port fuel injected (PFI) and direct injection spark ignition engines and the preliminary results are very encouraging.

Ameen, Muhsin Mohammed↗