Search NASA⌕ Search

SEARCH · Search NASA

Results for “initialization”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A MPET 2 -mPBPK model for subcutaneous injection of biotherapeutics with different molecular weights: From local scale to whole-body scale

Subcutaneous injection of biotherapeutics has attracted considerable attention in the pharmaceutical industry. However, there is limited understanding of the mechanisms underlying the absorption of drugs with different molecular weights and the delivery of drugs from the injection site to the targeted tissue. Here, we propose the MPET 2 -mPBPK model to address this issue. This multiscale model couples the MPET 2 model, which describes subcutaneous injection at the local tissue scale from a biomechanical view, with a post-injection absorption model at injection site and a minimal physiologically-based pharmacokinetic (mPBPK) model at whole-body scale. Utilizing the principles of tissue biomechanics and fluid dynamics, the local MPET 2 model provides solutions that account for tissue deformation and drug absorption in local blood vessels and initial lymphatic vessels during injection. Additionally, we introduce a model accounting for the molecular weight effect on the absorption by blood vessels, and a nonlinear model accounting for the absorption in lymphatic vessels. The post-injection model predicts drug absorption in local blood vessels and initial lymphatic vessels, which are integrated into the whole-body mPBPK model to describe the pharmacokinetic behaviors of the absorbed drug in the circulatory and lymphatic system. We establish a numerical model which links the biomechanical process of subcutaneous injection at local tissue scale and the pharmacokinetic behaviors of injected biotherapeutics at whole-body scale. With the help of the model, we propose an explicit relationship between the reflection coefficient and the molecular weight and predict the bioavalibility of biotherapeutics with varying molecular weights via subcutaneous injection. The considered drug absorption mechanisms enable us to study the differences in local drug absorption and whole-body drug distribution with varying molecular weights. This model enhances the understanding of drug absorption mechanisms and transport routes in the circulatory system for drugs of different molecular weights, and holds the potential to facilitate the application of computational modeling to drug formulation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Experimental and theoretical investigation into the high pressure deflagration products of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Diamond anvil cell (DAC) laser ignition experiments and reactive ab initio molecular dynamics (AIMD) simulations were performed on the high explosive (HE) LLM-105 to investigate its high pressure (HP) deflagration chemistry. Raman and optical spectroscopy measurements reveal LLM-105 reacts into an opaque carbonaceous product at 4–25 GPa. At pressures >~ 27 GPa, the reaction product consists of an amorphous optically transparent solid and nitrogen (N 2 ) in the solid phase. While not a one-to-one comparison due to the small time and length scales, the HP AIMD simulations show that some of the product is molecular N 2 , in qualitative agreement with experiment, while above 20 GPa most of the product consists of large amorphous C x H y N z O k clusters. Clustering is enhanced with pressure and reduces with temperature. In the experiments with initial sample pressure >~ 25 GPa, the pressure within the DAC decreases with minimal change in DAC cavity area. At initial sample pressures of 43.9 GPa, when quenched to 0 K, simulations predict a product experiencing a lower pressure consistent with the experimental measurement at lower load pressures. In conclusion, the results are important for understanding the HP deflagration chemistry of LLM-105.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mesoscale modeling of restructuring in high burnup UO 2 fuel

Here, this work aims to simulate the restructuring behavior observed in different regions of high burnup fuel, providing a first-of-its-kind restructuring model for the dark zone and rim region of high-burnup UO 2 fuel. We employed a grand-potential-based phase-field model to concurrently evaluate subgrain formation and the growth of fission gas bubbles within the fuel. An energy-based subgrain formation criterion was introduced to simulate the restructuring process. The effects of different initial conditions and different modeling parameters were systematically studied to capture how each of these parameters influences the characteristics of the restructured fuel. Subgrain formation was observed to begin around existing fission gas bubbles and proceed toward triple junctions, grain boundaries, and grain interiors. Restructuring was demonstrated to be influenced by a combination of initial dislocation densities, burnup rate, subgrain formation rate, and temperature. Under a given subgrain formation rate, the rate of restructuring increases with rising fuel temperature. A restructuring bias was observed within the microstructure, due to the variation in defect accumulation when comparing different grains. Microstructures corresponding to the dark zone and rim region can be obtained by parameterizing the model with the appropriate defect production rate, as determined based on the burnup rate and temperature. Furthermore, bubble size and distribution do not significantly affect the rate of restructuring. The predicted microstructures are consistent with experimental observations of the restructured regions. Finally, we present a correlation demonstrating the evolution of the restructuring volume fraction as a function of local burnup.

UO2↗
DOE OSTI2025

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗
DOE OSTI2025

Toward more-robust, AI-enabled subsurface seismic imaging for geotechnical applications

Non-invasive seismic imaging has the potential to cost-effectively evaluate large volumes of subsurface material to inform geotechnical site investigation. However, seismic imaging using full waveform inversion (FWI) requires significant computational time and is dependent on an initial starting model. As a result, FWI has not yet been widely adopted into geotechnical practice. Previous efforts, on relatively simple two-layered models, indicate that data-driven artificial intelligence (AI) models may be as effective as FWI at predicting 2D images of shear wave velocity (V s ). Furthermore, the AI model predictions can be made almost instantaneously after data acquisition and do not require an initial starting model. We examine the generality of these findings by developing a new AI model for subsurface seismic imaging, whereby we make several notable contributions. First, we architect a multimodal AI model that combines time- and frequency-domain representations of the seismic wavefield to predict a 50 m by 20 m subsurface image of V s . Second, we developed a new diverse dataset of 100,000 images with their corresponding seismic wavefields to train the AI model. Third, we propose four physics-informed data augmentations for data-driven seismic imaging. Fourth, we develop two prediction consistency tests to evaluate the model’s performance when the true subsurface is unknown. Our final model, which has been made publicly available, is capable of predicting a subsurface V s image from a single seismic wavefield with an average, mean absolute percent error (MAPE) of 24 %. The predictive model is applied to a field dataset and shown to be consistent with local geology and shear-wave refraction measurements from the same location.

Artificial intelligence↗
DOE OSTI2025

Delamination-informed lifecycle decisions: A dielectric and machine learning framework for composite sorting and recycling

Composite materials are widely used in aerospace, marine, and automotive sectors due to their high strength-to-weight ratio and durability. However, their long-term reliability can be compromised by damage accumulation. Specifically, delamination initiation serves as a precursor to structural failure, which is often difficult to detect during damage inspection. Identifying and sorting delamination initiation in samples not only increases operational safety while providing critical information for end-of-life decisions, which influences both the service life extension value and the efficiency of fiber extraction during recycling. This research addresses two challenges: (1) developing a nondestructive, ex-situ framework to sort composite materials based on damage severity, particularly delamination, and (2) understanding how damage in composites influences resin removal during pyrolysis. Both experimental work and finite element analysis were performed to predict critical stress levels that are associated with delamination onset. Based on these results, three loading levels 50 %, 75 %, and 90 % of maximum stress, were selected for controlled experiments, generating composite samples with varying extents of damage for machine learning model training. Microscopic imaging of these samples confirmed the damage progression from matrix cracking to delamination, validating the computational predictions. We explored supervised machine learning using dielectric measurements to classify damage states. Preliminary results show an artificial neural network can identify early delamination which is a potential precursor to failure, with 94.44 % accuracy on our dataset. A parallel investigation into the effect of damage severity on pyrolysis recycling showed that heavily delaminated samples required significantly less energy for comparable matrix removal than undamaged samples.

dielectric variables↗
DOE OSTI2025

Hot corrosion behavior of 304 & P91 graded composite transition joint under molten sulfate salts

A novel graded composite transition joint (GCTJ) between AISI 304 stainless steel and ASTM A335 P91 steel, has been demonstrated remarkablely superior creep performance compared to the conventional dissimilar metal weldment (DMW) under equivalent conditions. However, this advance is challenged by hot corrosion under salt deposition at elevated temperatures. Here, this study investigates the hot corrosion behavior of 304&P91 GCTJ exposed to sulfate salts at 700°C. Compared to the 304 steel, corrosion attacks initiate in the P91 triangle, where the dual-phase microstructure of ferrite and tempered martensite significantly influences pitting initiation. Anodic dissolution mainly occurs within the tempered martensite due to more vulnerable sites within the tempered martensite. With prolonged exposure, corrosion propagates across both phases, and the corrosion depth within the P91 triangle is related to the exposed surface area ratio between 304 and P91. Electrochemical analysis reveals the occurrence of galvanic corrosion between the 304 and P91, with a positive linear relationship between anodic dissolution current density (I a ) and the exposed surface area ratio between 304 and P91 in molten salts, further emphasizing the critical impact of this ratio on the corrosion severity in the P91 triangle. These findings underscore the importance of transition zone design optimization in mitigating localized corrosion.

Galvanic corrosion↗
DOE OSTI2025

The microscopic mechanisms of high temperature oxidation of Haynes 282

Nickel-based superalloy finds widespread applications in aerospace and extreme environments. They are known for their high temperature oxidation resistance under extended periods. However, the oxide formation and evolution which sets the stages for the later parabolic oxidation kinetics is not fully understood. This paper aims to provide new insights into the transient stage oxidation mechanism and kinetics of a typical Ni-based superalloy (Haynes 282). While tracking the chemical and microstructural evolution under micrometer scale at 800°C, we show that the oxide scale and its grain boundary species change significantly during the initial stages of oxidation and can have a profound impact on the oxidation kinetics. Cr 2 O 3 forms initially at the grain boundary along with minor amount of Al 2 O 3 . Then, the grain boundary region is enriched with copious amounts TiO 2 while Cr migrates away from the grain boundary. A noticeable change in the isothermal oxidation kinetics observed likely results from a mechanistic change from uniform surface oxidation to preferential outward diffusion of Ti 4+ ions through the grain boundary. Through first-principles calculations combined with energy dispersive spectroscopy, we confirmed the preferential outward diffusion of Ti 4+ ions as the diffusion barrier is lower for Ti than Cr along the grain boundaries. In conclusion, these findings highlight the critical role of early-stage grain boundary oxidation dynamics in dictating the long-term oxidation resistance of Ni-based superalloys and provide a foundation for future strategies to enhance their performance in extreme environments.

DFT↗
DOE OSTI2025

Stirring and peristaltic pumping alter flow electrode particle size and morphology

Flowable carbon-based electrodes are widely studied in electrochemical systems for energy storage and water treatment. Most laboratory- and pilot-scale studies employ active mixing for slurry preparation and handling and use peristaltic pumps for slurry circulation. Magnetic stirring is the most prevalent mixing method at the laboratory scale, and the combination of magnetic mixing and peristaltic pumping is favored for convenience and compatibility with particulate suspensions. However, these methods can strongly affect particle size and morphology (and therefore cell performance). In this study, we evaluate how mixing and pumping influence particle degradation in two commonly used carbon materials: spherical carbon black (CB) beads and powdered activated carbon (AC). We performed 24 h stirring-only experiments using magnetic, orbital, and overhead stirring. Separately, we conducted 32 h circulation experiments using a commercial peristaltic pump. Particle morphology and size were analyzed using optical microscopy and laser diffraction. The initially spherical CB beads (~580 μm, 1.5 wt%) exhibit rapid and extensive breakdown under both peristaltic pumping and magnetic stirring. By comparison, AC particles (initially ~31 μm) tested at 1.5 and 15 wt% showed measurable size reductions of up to ~25 % under either stirring or pumping. In conclusion, these findings highlight the need to account for and quantify morphological changes caused by both stirring and pumping in the design and associated evaluation of flowable electrode systems.

Chemistry↗
DOE OSTI2025

Cataloging US state policy patterns towards microgrid deployment

Frequent extreme weather events have called for rigorous and timely efforts for alternative non-wire solutions. These efforts are getting more widespread to offer a perfect alternative as the conventional grid becomes progressively less resilient. One of these solutions is microgrids that can disconnect from the grid and offer grid resilience during an outage. While this technology is still finding its footing in the industry, states across the US are employing policy patterns and forms of instruments to support its deployment. This study includes a systemic review of the US by conducting a binary analysis of all 50 states (including Washington D.C, excluding other US territories) using seven variables. The results show four major policy approaches to microgrids: i) supporting microgrids through a definitive legislative activity leading to further policy action; ii) direct efforts from the public utilities commissions without a concrete legislative push; iii) initiatives from institutions other than the commissions; and lastly, iv) self-initiated community and private consumer efforts. The results help understand what policy instruments are being used in each of these patterns to support this niche technology that still faces regulatory challenges.

Furqan, Maham↗
DOE OSTI2024

Rhenium Isotope Reconnaissance of Uranium Ore Concentrates

Exploration of natural isotopic variations of the element rhenium (Re) is in its infancy, with initial studies revealing isotopic fractionation in a variety of geological materials. Here, in this work, we investigate Re isotope variation as a new geochemical tool, given its redox-sensitive properties and affinity for organic matter and sulfides. In this work, Re abundance and isotope ratio data were collected from uranium ore concentrates (UOCs) across a variety of depositional ages, locations, geologic settings, and deposit types. Ore types from which the UOC were derived include sandstone, unconformity, and quartz-pebble (QP) conglomerate. To isolate Re from the U-rich matrix of UOCs, a new purification method utilizing DGA ion exchange resin was developed. We found that UOCs exhibit a wide range of Re isotope ratios, with sandstone ore-derived UOCs having the isotopically lightest values, QP conglomerate ore-derived UOCs having the heaviest, and unconformity ore-derived UOCs in between (with some overlap with sandstone UOCs). The Re isotope ratio range observed in UOCs extends previously reported values by more than a factor of two. Industrial processing (e.g., incomplete recovery of Re from ore, contamination, fractionation during processing) may play a role in the isotopic variability in the UOCs. However, systematic differences between ore types suggest that the depositional setting is a significant factor. For nuclear forensic investigations, Re isotopic compositions combined with data from other isotopic systems provide geochemical signatures that can aid in provenance assessment of UOCs. Regardless of the specific causes for the wide range of Re isotope ratios in UOCs, these initial data indicate Re is a promising tool for nuclear forensic investigations on samples from early in the nuclear fuel cycle.

58 GEOSCIENCES↗
DOE OSTI2024

Synthesis of hierarchical graphene coated porous Si anode for lithium-ion batteries

The ultra-high capacity and widespread availability of Si materials have resulted in them being the focus of extensive studies to replace the graphite anode. However, the main barriers preventing the Si anodes from large-scale applications are their huge volume change and severe interface reactions in the lithiation/delithiation process, which pulverizes its structure and subsequently deteriorates its cycle life. Here, micrometer-scale porous Si coated with graphene coating (mpSi@G) has been synthesized by using SiO 2 nanoparticles and novel coal-derived humic acid as feedstocks through a magnesiothermic reduction, followed by spray drying and calcination techniques. SEM, Raman, and X-ray absorption analysis demonstrate that the hierarchical graphene shell and micrometer-sized porous Si structure effectively release the Si anode's mechanical stress upon lithiation to achieve good structural stability. Here, the synthesized mpSi@G anode delivers a high initial lithiation capacity of 2974.9 mAh g –1 at 0.1 A g –1 with an initial coulombic efficiency of 70.2 %. Furthermore, the conductive hierarchical graphene network, along with the tight contacts of porous-Si and the graphene shell, contribute to a high capacity of 1109.5 mAh g –1 at a high current density of 5.0 A g –1 , showing excellent rate capability.

25 ENERGY STORAGE↗
DOE OSTI2025

Insights Into Thermal Runaway Mechanisms: Fast Tomography Analysis of Metal Agglomerates in Lithium-Ion Batteries

Thermal Runaway (TR) in lithium-ion batteries (LIB) is a critical technological and social concern. Whilst such events are rare, TR is characterized by uncontrollable heating leading to catastrophic failures. To deepen the understanding of the failure process and subsequently develop more accurate TR prediction models and as a result safer battery systems, we present in this work high-speed X-ray tomography for in-depth investigations of the copper current collector melting and agglomeration during TR. The melting process presents valuable real-time internal information about heat evolution during TR, previously challenging to access but crucially important for validating TR models. In this work, controlled failure studies combined with high-speed X-ray tomography were performed on two different commercial LIB models, subjecting them to both external heating and nail penetration to induce TR. Through real-time observation via high-speed tomography, followed by segmentation, rendering, and analysis, the formation of copper agglomerates was qualitatively and quantitatively characterized and visualized for the first time. Agglomerates tended to form either from the battery's outermost layers or centrally, depending on the method of TR initiation, and gives an indirect insight into the internal temperature evolution and distribution. Moreover, an initial comparative analysis between the battery models also revealed differences in agglomerate size, which has been linked to the thicker copper current collectors of one of the cell models. We further discuss the impact of larger copper agglomerates on heat distribution and safety. This study not only sheds light on the intricate dynamics of TR in LIBs but also underscores the pivotal role of 'gold-standard' imaging techniques in advancing battery safety, crucial for the robust modeling of TR and the future design of electric vehicle safety systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Mechanically induced thermal runaway severity analysis of Li-ion batteries and continuous energy release monitoring

The large-scale deployment of Li-ion batteries in stationary energy storage and electrical vehicle applications demands a strong focus on safety, particularly on the thermal runaway risk and severity evaluation. A standardized single-side mechanical indentation test protocol was developed to induce an internal short-circuit (ISC) and evaluate cells' thermal runaway severity at different state of charge (SOC). The observed hazard severity (OHS in five categories) and evaluated scores in this work have a comprehensive consideration of each cell's capacity, initial voltage, SOC, temperature and voltage change, allowing a better evaluation of the cells' thermal runaway potential. This method was applied to about 200 Li-ion batteries in order to build an extensive thermal runaway database covering various SOCs, capacities and chemistries. In this study, we monitored the transitions of stored electrochemical energy and applied mechanical energy into both thermal energy and acoustic emissions (AE). The surface temperature and mechanical failures were monitored by infrared imaging and AE to capture critical events within battery cells throughout the mechanical indentation tests. Furthermore, the initial temperature maps can predict two types of follow-up events: thermal runaway or gradual heat release via conduction. Analyzing each cell's severity, AEs, and leveraging the evolving database offer insights into predicting occurrences of thermal runaway. The test method, thermal runaway severity evaluation and prediction, and the corresponding database provide battery designers, manufacturers, and end-users a clear overview of Li-ion batteries' thermal runaway potential under mechanical abuse, advancing the safety design of Li-ion batteries.

Acoustic emission↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2024

A comprehensive numerical investigation on spray models for Direct-Injection Spark-Ignition engines

Gasoline direct-injection spark-ignition (DISI) engines generate a large portion of their unburned hydrocarbon (UHC) and soot emissions during the cold-start phase. A predictive computational fluid dynamics (CFD) modeling framework can be used to understand the physical processes that characterize fuel spray evolution and fuel-film formation at cold start conditions, which can help to reduce engine-out particulate emissions. This study systematically evaluated spray submodels and developed a set of simulation best practices for physical-numerical submodels with the goal of enabling accurate simulations of liquid spray behavior in a DISI engine. Three comprehensive experimental datasets containing free-spray projected liquid volume (PLV), liquid volume fraction (LVF), and near-field X-ray radiography data were used to validate the simulation results and evaluate the spray submodels. Systematic analysis delved into injected parcel distribution, droplet collision, spray breakup, and evaporation via a detailed assessment of the relevant spray submodels. Moreover, the effects of turbulence models and the initial turbulent flow properties on the liquid spray evolution were examined. Based on extensive calibration efforts, a set of simulation best practices for the free spray was developed and validated against the PLV/LVF data. Simulation results indicated that the uniform distribution for parcel initialization, coupled with appropriate droplet collision submodels, provides an improved spray morphology compared to the cluster distribution. The findings also underscored the importance of calibrating the Kelvin-Helmholtz Rayleigh-Taylor (KH-RT) breakup model constants and droplet heat transfer coefficient scaling factor to achieve favorable agreement regarding measured liquid penetration and spray widths. In conclusion, this study marks a substantial stride towards accurately predicting fuel film evolution and soot formation within DISI engine performance.

ECN Spray G↗
DOE OSTI2025

Vacuum vessel design with lofted toroidal surfaces for a QHS configuration

A vacuum vessel in between the hot plasma and the magnetic coils is crucial for a fusion reactor. In this study, a tool is presented to perform the initial design of such a vessel automatically. Starting from the last closed flux surface, an algorithm is developed which lofts this manifold out as much as possible. At the same time, the algorithm makes sure that no intersection with the coils takes place. In the presented work, this algorithm is applied to the HSX stellarator and it has been verified that the resulting vessel wall can be built from an engineering perspective. Using field line tracing, the strikepoints on the vessel wall indicate the initial locations for a divertor.

HSX↗