Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel performance 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,261 records · Page 70

Development of Instrumented Advanced Test Reactor Irradiation Capsule Experiment for In-situ Thermal Conductivity Measurements of High-Density Fuels

Idaho National Laboratory (INL) is developing a first-of-a-kind leadout instrumented capsule experiment design to enable in-situ measurement capabilities in the Advanced Test Reactor (ATR) core. The Ceramic Advanced Thermal Evolution Research (CRATER) experiment supports the aLEU program objective to accelerate fuel performance irradiation testing for identifying alternative high-assay, low enriched uranium (HALEU) fuel systems. CRATER is a fueled, instrumented capsule experiment to measure in-situ temperature and thermal conductivity of ceramic fuels. Two ceramic fuel types will be used, uranium mono-nitride (UN) and uranium mono-carbide (UC), with a third metallic fuel used for comparison (UMo). The three fuel specimens will use a stainless-steel cladding. Programmatic objectives include linear heat generation rates (LHGR) of 210 ± 25 Watts per cm. and an inner clad temperature of 300-450 °C. The evolution of fuel thermal conductivity during irradiation has never been successfully measured in-situ for these systems and this experiment is designed to use advances in measurement sciences to characterize how thermal transport properties evolve while in reactor. Neutronic simulations of the experiment and its surrounding reactor environment were conducted using the Monte Carlo N-Particle Transport code (MCNP) and result in optimized fuel enrichment to meet target linear heat generation rates (LHGRs) influencing fuel temperatures, and fuel burnup requirements. Fabrication research and development (R&D) efforts are underway to produce annular right cylinder UC and UN pellets using carbothermic reduction and nitridation (or hydride-dehydride-nitride) synthesis methods, followed double-action die cold isostatic pressing.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development & Validation of Low-Cost, Highly-Durable, Spinel-Based Materials for SOFC Cathode-Side Contact (Final Report)

A cathode-side contact layer is required to provide and maintain stable electrical conduction paths between the interconnect and cathode in a solid oxide fuel cell (SOFC) stack assembly and thus minimize the ohmic resistance and stack power loss. Current cathode-interconnect contact materials are based on noble metals, electrically-conductive perovskites, their composite materials, etc. These materials are either too expensive or do not possess the overall balanced performance required for the cathode-side contact application. To achieve the DOE SOFC system cost and performance stability goals, a new generation of low-cost, high-performance contact materials needs to be developed. In this project, spinel-based materials thermally converted from the Fe-Ni and Co-Mn based alloy precursors were developed and validated for the cathode-side contact application. The precursor alloy compositions were optimized via a combination of composition screening in the (Ni,Fe) 3 O 4 and (Mn,Co) 3 O 4 spinel system, alloy design using physical metallurgy principles, and cost considerations. The alloy powders with the desired composition and particle size were manufactured via gas atomization. The optimal process parameters for thermal conversion of these alloy precursor layers to a spinel-based layer were identified, i.e., 900°C x 2h in air, which is close to the initial stack firing condition. The area-specific resistances (ASRs) of the interconnect/contact/cathode test assemblies with the developed contact layer were determined for various durations (up to 5000 h) under simulated cathodic operation conditions. Some of the alloy-derived spinel contacts exhibited the lowest ASR and ASR degradation rate. The in-stack performance of the most promising alloy-derived contact layer is currently being evaluated via stack testing. To reduce the stack cost, the Co-Mn based alloy powders were utilized as the precursor for synthesis of dense spinel-based interconnect coating. By optimizing both the initial powder size/distribution and the alloy powder composition, a dense (Mn,Co) 3 O 4 -based spinel coating was achieved. Furthermore, co-sintering of the coating/contact dual-layer structure under the initial stack firing condition was realized by utilizing the tailored Co-Mn alloy precursors. Cost analysis of the developed technology indicated a total stack cost reduction of around 10.6% with the implementation of co-sintering of the interconnect coating and the contact layer during initial stack firing. Since low-cost processes such as screen printing is utilized in the precursor application and no reduction heat treatment is needed for the coating formation, the developed technology can be readily implemented at the industrial partner’s manufacturing facilities with no additional capital investment needed.

08 HYDROGEN↗

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly 4 increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co 2 O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO 2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres

01 COAL, LIGNITE, AND PEAT↗

Impacts of uncertain feedstock quality on the economic feasibility of fast pyrolysis biorefineries with blended feedstocks and decentralized preprocessing sites in the Southeastern United States

Abstract This study performs techno‐economic analysis and Monte Carlo simulations (MCS) to explore the effects that variations in biomass feedstock quality have on the economic feasibility of fast pyrolysis biorefineries using decentralized preprocessing sites (i.e., depots that produce pellets). Two biomass resources in the Southeastern United States, that is, pine residues and switchgrass, were examined as feedstocks. A scenario analysis was conducted for an array of different combinations, including different pellet ash control levels, feedstock blending ratios, different biorefinery capacities, and different biorefinery on‐stream capacities, followed by a comparison with the traditional centralized system. MCS results show that, with depot preprocessing, variations in the feedstock moisture and feedstock ash content can be significantly reduced compared with a traditional centralized system. For a biorefinery operating at 100% of its designed capacity, the minimum fuel selling price (MFSP) of the decentralized system is $3.97–$4.39 per gallon gasoline equivalent (GGE) based on the mean value across all scenarios, whereas the mean MFSP for the traditional centralized system was $3.79–$4.12/GGE. To understand the potential benefits of highly flowable pellets in decreasing biorefinery downtime due to feedstock handling and plugging problems, this study also compares the MFSP of the decentralized system at 90% of its designed capacity with a traditional system at 80%. The analysis illustrates that using low ash pellets mixed with switchgrass and pine residues generates a more competitive MFSP. Specifically, for a biorefinery designed for 2,000 oven dry metric ton per day, running a blended pellet made from 75% switchgrass and 25% pine residues with 2% ash level, and operating at 90% of designed capacity could make an MFSP between $4.49 and $4.71/GGE. In contrast, a traditional centralized biorefinery operating at 80% of designed capacity marks an MFSP between $4.72 and $5.28.

09 BIOMASS FUELS↗

Computational Fluid Dynamic Modeling of Dry Cask Simulator with Crosswind

The purpose of this study is to create a STAR-CCM+ model of a Belowground Vertical Dry Cask Simulator (BVDCS) at Sandia National Laboratories (SNL) and validate the model with SNL’s experimental results. The BVDCS consists of a single boiling water reactor assembly fitted with electric heaters encompassed by a containment vessel and shell to represent a belowground spent nuclear fuel (SNF) dry storage system. Blowers are located near the inlet and outlet of the BVDCS to simulate crosswind conditions. In addition to the experimental results, the STAR-CCM+ model developed for this study is compared with a previous computational fluid dynamics (CFD) model in a different software program, which is used as a software-to-software benchmark. The experimental results provide a dataset to compare the STAR-CCM+ model results for a variety of different conditions. The main objective is to validate and improve STAR-CCM+ CFD models for spent nuclear fuel storage systems with explicitly modeled external environments and “wind driven” crossflows. These CFD models aide in the study of external particle deposition in spent nuclear fuel storage systems, which is important to predicting the significance of chloride induced stress corrosion cracking (CISCC). In addition to experimental comparison, a sensitivity analysis study is performed using the STAR-CCM+ model. The sensitivity analysis provides a quantitative assessment of the sensitivity of various parameters. This helps provide information on various parameters that are of particular importance to constructing a model representative of real life systems. The STAR-CCM+ model compared well to the experimental results showing similar responses to changes in cross wind flow, and a number of parameters are identified for model improvement.

Jensen, Ben J.↗

Effect of the J-Shaped Wind Turbine Airfoil Opening Ratio and Thickness on the Performance of Symmetrical Airfoils

Abstract In this work, the effect of the inner opening ratio on the J-shaped airfoils aerodynamic performance was studied and documented for symmetrical airfoils. Three different airfoil thicknesses were investigated: small (NACA0008), medium (NACA0015), and large (NACA0024). For each airfoil thickness, effects of three inner opening ratios were analyzed: one-third, one-half, and two-thirds. The performance of each opening ratio was compared with the performance of the solid airfoil “zero opening ratio” for different angles of attack between 5 deg and 20 deg. All designs were simulated using the computational fluid dynamics (CFD) technology against experimental results for solid NACA4412 airfoil in the University of Wisconsin-Milwaukee (UWM) wind tunnel facility and other published experimental data. It was found that large eddy simulation yields accurate solutions with a smaller number of mesh cells compared to the k–ω turbulence model but with much longer computational time. The lift-to-drag ratio for all studied airfoils has a maximum value for solid airfoils compared to those equipped with openings. For airfoils equipped with 00.00% opening ratio “solid,” NACA0015 airfoil has the maximum lift-to-drag ratio. Furthermore, it was found that NACA0008 equipped with a 33.33% opening ratio has the best performance of all studied J-shaped airfoils.

Energy & Fuels↗

Interpretation of Hydraulic Fractures Based on Microseismic Response in the Marcellus Shale, Monongalia County, West Virginia, USA: Implications for Shale Gas Production

Summary Hydraulic fracturing is critical for extracting shale gas in the subsurface. The treatment technique of multistage hydraulic fracturing is widely used to maximize production. In multistage hydraulic fracturing, not all pumping stages make the same contribution for production, although the designed stimulation process is almost same in every pumping stage. In this study, we characterize the microseismic responses of hydraulic fracturing. Significant variations of microseismic characteristics are observed among the different pumping stages. Pre-existing natural fractures are examined along the horizontal well in selected stages. Combined with the Mohr’s circle analysis, results of the fracture study show that induced hydraulic fractures can be captured by pre-existing natural fractures. Induced hydraulic fractures are simulated by the unconventional fracture model (UFM), and the result reveals that the stress-shadow effect diverts the direction of hydraulic fractures. The diverted hydraulic fracture affects the development of the hydraulic fracture network, which has influence on production. The treatment of two-step pumping is investigated by application of hydraulic diffusivity. Hydraulic fracturing performance and production could be optimized by the treatment of two-step pumping in a single stage. The second pumping creates new fractures and fills the fractures with additional proppants to maintain production for a long duration.

Energy & Fuels↗

Highly-Active and Contaminant-Tolerant Cathodes for Durable Solid Oxide Fuel Cells (Final Report)

The objective of this project is to investigate the fundamental degradation methods occurring in solid oxide fuel cell (SOFC) cathodes when exposed to chromium and carbon dioxide contamination and to rationally design alkaline earth-based catalysts to increase stability and tolerance to contaminate poisoning. With a mechanistic and fundamental understanding of the degradation methods, advanced catalytic surface modifications can be applied to reduce degradation. The specific technical objectives are: (1) To identify/develop new catalysts (alkaline-earth based-) that are compatible chemically with the state-of-the-art cathode materials at high temperatures required for fabrication and with contaminates commonly encountered under operating conditions; (2) To improve the infiltration process for optimal control of the thickness, composition, and uniformity of the catalyst coatings; (3) To evaluate the electro-catalytic activity toward ORR of the chemically-stable materials when exposed to different types of contaminants using electrical conductivity relaxation measurements on bar samples and performance evaluation of catalyst-infiltrated cathodes; (4) To unravel the contamination-tolerant mechanisms of the new catalyst coatings under realistic environmental conditions (with different types of contaminants) using powerful in situ and in operando characterization techniques performed on model cells with thin-film/pattern electrodes, as guided by modeling and simulation; (5) To establish scientific basis for rational design of new catalysts of high tolerance to contaminants; (6) To validate the long term stability of modified LSCF cathodes in commercially available cells/stacks under ROC. Alkaline earth-based catalysts have been systematically explored and tested under various contamination conditions. BaCoO 3-δ (BCO) was shown to produce the best catalytic activity enhancement as well as stability in a variety of contaminating conditions, including CO 2 and chromium. The microstructural evolution was investigated with SEM, EDS, and Raman spectroscopy, showing the BCO catalyst prevents the formation of insulating SrCrO 4 by forming electrically conductive BaCrO 4 . Electrochemical relaxation measurements determined the BCO catalyst coatings increased the surface diffusion coefficient but did not significantly affect the diffusion coefficient and determined an optimum surface modification layer of 100 nm. A novel Ba 0.9 Co 0.7 Fe 0.2 Nb 0.1 O x (BCFN) catalyst was also shown to have excellent stability in chromium containing atmospheres. A novel surface sol-gel (SSG) surface modification, which offers superior thickness and compositional control, was applied to SOFCs to produce BaO catalyst coatings. SSG BaO coatings were shown to produce conformal coatings on the electrode surface, greatly increasing stability in chromium containing atmospheres. An atomic level mechanistic approach was applied to investigate the stability of PrBa 0.8 Ca 0.2 Co2O 5+δ (PBCC) with respect to a variety of common contaminants, demonstrating thermodynamically its superior stability in CO2 containing atmospheres. Finally, the best catalyst coatings demonstrated in this project, BCO and PBCC, were applied to full cells which demonstrated superior stability for over 300 hours in CO 2 and Cr atmospheres.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Novel Deep Reinforcement Learning Approach to Traffic Signal Control with Connected Vehicles

The advent of connected vehicle (CV) technology offers new possibilities for a revolution in future transportation systems. With the availability of real-time traffic data from CVs, it is possible to more effectively optimize traffic signals to reduce congestion, increase fuel efficiency, and enhance road safety. The success of CV-based signal control depends on an accurate and computationally efficient model that accounts for the stochastic and nonlinear nature of the traffic flow. Without the necessity of prior knowledge of the traffic system’s model architecture, reinforcement learning (RL) is a promising tool to acquire the control policy through observing the transition of the traffic states. In this paper, we propose a novel data-driven traffic signal control method that leverages the latest in deep learning and reinforcement learning techniques. By incorporating a compressed representation of the traffic states, the proposed method overcomes the limitations of the existing methods in defining the action space to include more practical and flexible signal phases. The simulation results demonstrate the convergence and robust performance of the proposed method against several existing benchmark methods in terms of average vehicle speeds, queue length, wait time, and traffic density.

42 ENGINEERING↗

Summary Report Of The FY25 Computational Fluid Dynamics Verification And Validation Exercises In The Advanced Reactor Technologies - Gas-cooled Reactor Program

Verification and Validation (V&V) of numerical tools is critical for ensuring reasonable predictions during design, safety analysis, and licensing. Recent work in the Advanced Reactor Technologies - Gas-cooled Reactor (ART-GCR) program has focused on V&V of common Computational Fluid Dynamics (CFD) tools that are used within the Untied States. This report presents an update on these CFD V&V activities. These Generation IV Forum (GIF) Very High Temperature Reactor (VHTR) Computational Methods, Validation, and Benchmarking (CMVB) is an international organization dedicated to the verification and validation of High Temperature Gas-Cooled Reactor (HTGR) simulation tools. Participation in the CMVB provides additional value to the V&V activities, as it allows for access to a wider range of data, and provides valuable benchmarking exercises. Three HTGR phenomena are targeted: Reactor Cavity Cooling System (RCCS) performance, core bypass flow, and lower plenum mixing. Simulations of the University of Wisconsin-Madison (UW-Madison) RCCS facilities are performed with Reynolds Averaged Navier-Stokes (RANS) in StarCCM+. Results are compared for both forced and natural convection conditions, with both exhibiting good agreement with experimental measurements. The Idaho National Laboratory (INL) matched index of refraction (MIR) and Korean Atomic Energy Research Institute (KAERI) bypass flow expeirments are used to validation CFD predictions of bypass flow. Simulations are performed with RANS in StarCCM+ and with Large Eddy Simulation (LES) in NekRS. Finally, preliminary simulations of the Institute of Nuclear and New Energy Technology (INET) lower plenum mixing facilities are presented. Initial work has developed models with LES, RANS, and porous media models. These preliminary models are presented and compared to each other to gauge differences in predictions with each of the three methods.

and Benchmarking (CMVB)↗

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling and optimization of steady flow of natural gas and hydrogen mixtures in pipeline networks

Here, we extend the canonical problems of simulation and optimization of steady-state gas flows in pipeline networks with compressors to the transport of mixtures of highly heterogeneous gases injected throughout a network. Our study is motivated by proposed projects to blend hydrogen generated using clean energy into existing natural gas pipeline systems as part of efforts to reduce the reliance of energy systems on fossil fuels. Flow in a pipe is related to endpoint pressures by a basic Weymouth equation model, with an ideal gas equation of state, where the wave speed depends on the hydrogen concentration. At vertices, in addition to mass balance, we also consider mixing of incoming flows of varying hydrogen concentrations. The problems of interest are the heterogeneous gas flow simulation (HGFS), which determines system pressures and flows given fixed boundary conditions and compressor settings, as well as the heterogeneous gas flow optimization (HGFO), which extremizes an objective by determining optimal boundary conditions and compressor settings. We examine conditions for uniqueness of solutions to the HGFS, as well as compare and contrast mixed-integer and continuous nonlinear programming formulations for the HGFO. We develop computational methods to solve both problems, and examine their performance using four test networks of increasing complexity.

08 HYDROGEN↗

Modeling and control of nuclear–renewable integrated energy systems: Dynamic system model for green electricity and hydrogen production

The need for decarbonization and diversification of energy resources has led to the development of integrated energy systems (IESs), where multiple resources supply more than one energy sector. Here, one such IES with small modular nuclear reactors and renewables (wind and solar) as generating resources, catering to the demand of the electric grid while producing hydrogen for industries, is modeled in this paper. The physics-based component models are represented using the Modelica language and interconnected to form the IES. The control and coordination of the overall system are ensured by designing a suitable control architecture composed of individual subsystem-level controls and supervisory control. The dynamic performance and the load-following capability of the IES are evaluated, while satisfying the safe operational limits of the components. Different configurations and modes of IES operation are considered, where the adaptability of the control system in the presence of varying demands and renewable generations is validated. The simulation results indicate that hydrogen as a flexible load facilitates the supply of varying grid demand. Additionally, the renewables are also accommodated into the IES owing to the flexibility of the balance of plant associated with the nuclear reactors.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DPC Direct Disposal Postclosure Thermal Modeling

Performance of geologic radioactive waste repositories depends on near-field and far-field processes, including km-scale flow and transport in engineered and natural barriers, that may require simulations of up to 1 M years of regulatory period. For a relatively short time span (less than 1000 years), the thermohydro-mechanical-chemical (THMC) coupled processes caused by heat from the waste package will influence near-field multiphase flow, chemical/reactive transport, and mechanical behaviors in the repository system. This study integrates the heat-driven perturbations in thermo-hydro-mechanical characteristics into thermo-hydro-chemical simulations using PFLOTRAN to reduce dimensionality and improve computational efficiency by implementing functions of stress-dependent permeability and saturation-temperature-dependent thermal conductivity. These process couplings are developed for spent nuclear fuel in dual-purpose canisters in two different hypothetical repositories: a shale repository and a salt repository.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗

EVALUATION OF PRESSURE DROP CORRELATIONS FOR THE WIRE-WRAPPED ROD BUNDLES

The accurate prediction of pressure drop in the fuel rod bundle is of importance for the design and thermal-hydraulic analysis of nuclear reactor system. Due to very complex geometry of wire-wrapped fuel rod bundle geometry for the Sodium Fast Reactor (SFR), a simple method using hydraulic (equivalent) diameter is not sufficient to predict the pressure drop of rod bundle accurately. For several decades, many experimental works have been performed to develop an accurate correlation for wire-wrapped fuel bundle friction factor. In these works, the friction factor correlation of wire-wrapped fuel bundle has been developed as a function of Reynolds number of the flow and the design parameters such as fuel pin diameter, wire diameter, rod pitch, and wire lead length. The Cheng and Todreas correlations are well-known and used in many subchannel analysis code, such as CADET, TRIO-U, MATRA-LMW, SE2-ANL, etc. Novendstern model, Rehme model, Baxi and Dalle Donne model are also well-known and widely referred. Computational Fluid Dynamics is another typical approach to predict the pressure drop of wire-wrapped fuel rod bundle, but it requires the model validation because the accuracy of pressure drop prediction by the CFD simulation depends on what turbulence model is employed. The present work aims at evaluating the existing pressure drop correlations for the wire-wrapped fuel rod bundle and evaluating the pressure drop of 217-pin rod bundle by the CFD analysis.

99 GENERAL AND MISCELLANEOUS↗

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗