Search NASA⌕ Search

SEARCH · Search NASA

Results for “Internal Flow”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Baroclinic Residual Circulation and Mass Transport Due to Internal Tides

Baroclinic (BC) tidal residual circulation due to internal tides is investigated around islands over a shallow ridge using a numerical ocean model. Internal tides enhance vertical mixing over shallow slopes, leading to horizontal density gradients that drive BC residual circulation along the main thermocline. For a strongly stratified summer case, the vertical diffusivity estimated by the Mellor and Yamada turbulence closure model exceeds 1 × 10 -2 m 2 s -1 , and the velocity of BC residual circulations reaches 0.2 m s -1 . The magnitude of BC residual circulation is larger than that of barotropic residual tidal circulation, implying that BC residual circulation due to internal tides plays an important role in coastal ocean circulation. Furthermore, BC residual flow accounts for an equal percentage (5%) of the total tidal kinetic energy as the barotropic residual flow under summer stratification conditions. Results from a coupled sediment resuspension and transport model show the growth of intermediate nepheloid layers formed by strong bottom shear stress and BC residual flow. The residual component contributes largely (23% of the total in summer) to the total suspended sediment flux. Seasonal variability is explored, with weaker winter stratification leading to reduced mixing, and thus weaker BC residual circulation and sediment flux. Lastly, the magnitude of the BC residual circulation is also shown to be proportional to the square of the tidal amplitude.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Durable and Highly Selective Ion Transport of a Sulfonated Diels Alder Poly(phenylene) for Vanadium Redox Flow Batteries

In order to improve the durability and ion selectivity of a hydrocarbon membrane for vanadium redox flow batteries (VRFBs) a polymer was rationally designed with an external hydrophobic shell and internal hydrophilic core. The polymer was designed to prevent hydrophilic polymer chain aggregation by functionalizing the external polymer shell with hydrophobic side chains and attaching acid moieties onto the polymer backbone. The hydrophobic shell is the result of pentafluorobenzoyl group functionalization on the pendent aryl rings in a Diels Alder poly(phenylene) backbone. The internal polymer chain contains sulfonic acid moieties to impart hydrophilic character. The physical and electrochemical characteristics of the membrane were systematically studied and compared with widely used Nafion-212. The synthesized membranes were found to be superior to Nafion-212 not only in terms of lower vanadium ion permeability (4.25 x 10 -8 vs. Nafion-212: 6.0 x 10 -6 cm 2 min -1 ) but also in higher VRFB performance (EE = 89% and vs. Nafion-212: 85.5% at 40 mA/cm 2 ). The membrane also showed over three times higher capacity retention than Nafion-212 (0.075 % vs. 0.25 % capacity loss per cycle respectively) and excellent ex situ and in situ chemical stability.

25 ENERGY STORAGE↗
DOE OSTI2024

Integrated Turbine Component Cooling Designs Facilitated by Additive Manufacturing and Optimization

This research program involved the development of improved gas turbine cooling by optimizing the combined internal and film cooling configurations used for turbine blades and vanes. Improved cooling of turbine airfoils allows the gas turbines to be operated at higher temperatures which improves efficiency for gas turbine-based power generation. This also facilitates operation with hydrogen-based fuels which inherently generate higher operating temperatures. Designs of the new cooling configurations were based on the added geometric flexibility that is available when using metal additive manufacturing techniques. When making these designs, adjustments were made to account for imperfections occurring in the metal additive manufacturing process. New designs for film cooling holes were developed using computational adjoint gradient optimization techniques which generated a unique film cooling hole with 70% greater film cooling effectiveness compared to conventional shaped film cooling holes. A major factor in improving the film cooling effectiveness were small external protrusions on either side of the new hole which generated vortices that pushed coolant towards the surface and increased the lateral spreading of the coolant. Our studies also included optimization of film cooling holes that can be built with conventional manufacturing techniques and optimizing the internal coolant channels that supply coolant to the film cooling holes. Various internal coolant channel shapes were investigated and designs that maximized the cooling of the airfoil while minimizing pressure losses in the channel flows were identified. Testing and verification of the performance of the final configurations was done with a transonic wind tunnel facility which allowed testing at engine realistic high velocities. Final testing of the combined internal and film cooling configurations was done by incorporating the configurations in Penn State’s National Experimental Turbine (NExT) vane. These tests confirmed the enhanced overall cooling effectiveness provided by the new cooling configuration designs.

03 NATURAL GAS↗
DOE OSTI2021

Model Reduction for Steady Hypersonic Aerodynamics via Conservative Manifold Least-Squares Petrov–Galerkin Projection

High-speed aerospace engineering applications rely heavily on computational fluid dynamics (CFD) models for design and analysis. This reliance on CFD models necessitates performing accurate and reliable uncertainty quantification (UQ) of the CFD models, which can be very expensive for hypersonic flows. Additionally, UQ approaches are many-query problems requiring many runs with a wide range of input parameters. One way to enable computationally expensive models to be used in such many-query problems is to employ projection-based reduced-order models (ROMs) in lieu of the (high-fidelity) full-order model (FOM). In particular, the least-squares Petrov–Galerkin (LSPG) ROM (equipped with hyper-reduction) has demonstrated the ability to significantly reduce simulation costs while retaining high levels of accuracy on a range of problems, including subsonic CFD applications. This allows LSPG ROM simulations to replace the FOM simulations in UQ studies, making UQ tractable even for large-scale CFD models. This work presents the first application of LSPG to a hypersonic CFD application, the Hypersonic International Flight Research Experimentation 1 (HIFiRE-1) in a three-dimensional, turbulent Mach 7.1 flow. This paper shows the ability of the ROM to significantly reduce computational costs while maintaining high levels of accuracy in computed quantities of interest.

42 ENGINEERING↗
DOE OSTI2021

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

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

Modeling active nematics via the nematic locking principle

Active nematic systems consist of rod-like internally driven subunits that interact with one another to form large-scale coherent flows. They are important examples of far-from-equilibrium fluids, which exhibit a wealth of nonlinear behavior. This includes active turbulence, in which topological defects in the nematic order braid around one another in a chaotic fashion. One of the most studied examples of active nematics consists of a dense two-dimensional layer of microtubules, crosslinked by kinesin molecular motors that inject extensile deformations into the fluid. Though numerous theoretical studies have modeled microtubule-based active nematics, no consensus has emerged on how to fully and quantitatively capture the features of the experimental system. Here, to better understand the theoretical foundations for modeling this system, we propose a fundamental principle we call the nematic locking principle—individual microtubules cannot rotate without all neighboring microtubules also rotating. Physically, this is justified by the high density of the microtubules, their elongated nature, and their corresponding steric interactions. We assert that the nematic locking principle holds throughout the majority of the material, but breaks down in the neighborhood of topological defects and other regions of low density. We derive the most general nematic transport equation consistent with this principle and also derive the most general term that violates it, introducing fracturing into the material. We then examine the standard Beris–Edwards approach, commonly used to model this system, and show that it violates the nematic locking principle throughout the majority of the material due to fracturing. We then propose a modification to the Beris–Edwards model that enforces nematic locking nearly everywhere. This modification shuts off fracturing except in regions where the order parameter (a proxy for density) is reduced. In these regions fracturing is turned on. The resulting simulations in turn show strong nematic locking throughout the bulk of the material, with narrow bands of fracturing, consistent with experimental observation. One additional advantage of enforcing nematic locking is that nontrivial stationary state solutions, common in Beris–Edwards simulations but not seen in experiments, are eliminated.

Mitchell, Kevin A. [Univ. of California, Merced, C↗
DOE OSTI2025

Design and Engineering of LUPIN: A Test-Bed Radio-Frequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is a radio-frequency (RF) inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D neutral beam injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa ⋅ L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Finally, results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗
DOE OSTI2020

Water and vapor transport in algal‐fungal lichen: Modeling constrained by laboratory experiments, an application for Flavoparmelia caperata

ABSTRACT Algal‐fungal symbionts share water, nutrients, and gases via an architecture unique to lichens. Because lichen activity is controlled by moisture dynamics, understanding water transport is prerequisite to understand their fundamental biology. We propose a model of water distributions within foliose lichens governed by laws of fluid motion. Our model differentiates between water stored in symbionts, on extracellular surfaces, and in distinct morphological layers. We parameterize our model with hydraulic properties inverted from laboratory measurements of Flavoparmelia caperata and validate for wetting and drying. We ask: (1) Where is the bottleneck to water transport? (2) How do hydration and dehydration dynamics differ? and (3) What causes these differences? Resistance to vapor flow is concentrated at thallus surfaces and acts as the bottleneck for equilibrium, while internal resistances are small. The model captures hysteresis in hydration and desiccation, which are shown to be controlled by nonlinearities in hydraulic capacitance. Muting existing nonlinearities slowed drying and accelerated wetting, while exaggerating nonlinearities accelerated drying and slowed wetting. The hydraulic nonlinearity of F. caperata is considerable, which may reflect its preference for humid and stable environments. The model establishes the physical foundation for future investigations of transport of water, gas, and sugar between symbionts.

Potkay, Aaron↗
DOE OSTI2020

Cold plate with combined inclined impingement and ribbed channels

Heat transfer devices and methods for making the same that include a first enclosure having at least one inlet port; a second enclosure having a bottom plate and one or more dividing walls to establish channels, at least one internal surface of each channel having rib structures to create turbulence in a fluid flow; and a jet plate connecting the first enclosure and the second enclosure having impinging jets that convey fluid from the first enclosure to the channels, said impinging jets being set at an angular deviation from normal to cause local acceleration of fluid and to increase a local heat transfer rate.

42 ENGINEERING↗
DOE OSTI2023

HIGH-LOW FIDELITY THERMAL HYDRAULIC COUPLING USING AI/MACHINE LEARNING ALGORITHMS

The primary goal of the US Department of Energy (DOE) office of Nuclear Energy Integrated Energy Systems (IES) program is to develop the tools and framework for coupling multi-scale and multi-physical thermal and electrical energy usage and storage systems. High- and low-fidelity (high–low) coupling is a key feature of multi-scale, multi-component systems and has been an important focus of research in the nuclear energy community for the past two decades. An essential feature of demonstrating the capability to couple high-fidelity and low-fidelity systems for real-time applications are surrogate/reduced order models (ROM). For the purposes of this study, surrogate models are essentially Blackbox models, typically developed using supervised Machine learning (ML) algorithms. The surrogate models can be used to mimic the response of high-fidelity models to represent large historical datasets and coupled with more general low-fidelity system models distributed as Functional Mock-up Interface (FMI) or Functional Mock-up Units (FMU) modules. The example is demonstrated with Spallation Neutron Source (SNS) First Target Station flow loop data. The flow loop is a liquid mercury loop with a pump, piping, heat exchange, and internal heat generation in the target window. This work elucidates some of the potential benefits and future needs of developing tools for high–low system coupling of energy systems.

Williams, Wesley↗
DOE OSTI2024

Design and Engineering of LUPIN: A Test-Bed Radiofrequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is an RF inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D Neutral Beam Injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10 s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa·L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗
DOE OSTI2019

Gas-Induced Segregation in Pt-Rh Alloy Nanoparticles Observed by In Situ Bragg Coherent Diffraction Imaging

Bimetallic catalysts can undergo segregation or redistribution of the metals driven by oxidizing and reducing environments. Bragg coherent diffraction imaging (BCDI) was used to relate displacement fields to compositional distributions in crystalline Pt-Rh alloy nanoparticles. Three-dimensional images of internal composition showed that the radial distribution of compositions reverses partially between the surface shell and the core when gas flow changes between O 2 and H 2 . Our observation suggests that the elemental segregation of nanoparticle catalysts should be highly active during heterogeneous catalysis and can be a controlling factor in synthesis of electrocatalysts. Additionally, our study exemplifies applications of BCDI for in situ 3D imaging of internal equilibrium compositions in other bimetallic alloy nanoparticles.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Validation of RELAP-7 for Forced Convective and Natural Circulation Reactor Flows

New reactor designs such as small modular reactors and newly proposed passive safety systems rely on natural circulation two-phase flow cooling of the reactor core under accident conditions. Without a pump to drive the flow, the circulation rate is a strong function of void fraction. Under these conditions where the flow rate, pressure drop, and heat transfer are all closely coupled through the mixture density, the accurate prediction of void fraction is pivotal. Furthermore, while these passive systems boast a more reliable cooling option, the sensitivity of the system to two-phase flow instabilities is more likely than in forced convective cooling and can result in dangerous density wave or flow excursion conditions. Improvement of one-dimensional system codes and further development of three-dimensional code capability are all important efforts in addressing these issues. The proposed work is a critical step towards providing confidence through validation of the capability of thermal-hydraulic tools to meet the emerging challenges in reactor simulation. Historically, system analysis codes develop a stiffness over time due to fine-tuning of adjustable code parameters that result in a reluctance or inability to adopt improved models. Therefore, a rigorous effort should be invested in preparation for RELAP-7 to identify required modeling improvements of constitutive relations and validation of the code applicability to the wide range of possible reactor flow conditions. Part of this validation involves the uncertainty quantification which is a critical step in the Code Scaling, Applicability, and Uncertainty (CSAU) methodology and Best-Estimate Plus Uncertainty (BEPU) licensing method. BEPU allows recovery of safety margins and potentially optimizing plant performance, resulting in economic benefit for the reactor operators. Beginning with the collection of past experimental data covering prototypic LWR conditions, this document reviews three critical datasets available for the validation of system codes. The experimental methods, procedures, facilities, instrumentation, and supplementary information are consolidated to aid the validation efforts of this work and future projects. These three datasets (824 conditions) are chosen due to their wide acceptance and demonstrated application for code validation. Two experimental datasets from the Nuclear Power Engineering Corporation (NUPEC) span prototypic conditions of Boiling Water Reactors (BWR) and Pressurized Water Reactors (PWR). These datasets are labeled as the BFBT (BWR Full-size Fine-mesh Bundle Test) and the PSBT (PWR Subchannel and Bundle Test) datasets respectively. The third dataset is taken in an internally heated annulus that is scaled based on a BWR subassembly and includes two-phase flows important to post-blowdown analysis. Together, the datasets consider the full range of pressure and mass flux expected in LWR forced convective flows from steady-state to long term cooling conditions. Recognizing the deficiency of available data under low-pressure, low-flow, a new natural circulation dataset consisting of 107 conditions is obtained. With system pressures ranging between 145 and 950 kPa, natural circulation is driven by a heat flux up to 275 kW/m2 in a 3 m long heated region of the annulus test section which is followed by a 2 m long unheated region before entering a condenser. Stable natural circulation flow rates in this facility were measured to be between 182 and 590 kg/m2-s. Pressure and temperature are measured at ten locations in the facility. In addition, measurements of void fraction, interfacial area concentration, and gas velocity are obtained at five axial locations in the annulus test section. The RELAP5-3D code demonstrates good accuracy for void fraction predictions in BFBT and PSBT benchmarks, but shows relatively large discrepancy in the annulus experiments at lower pressures. In two-phase natural circulation, large oscillations are observed in many system parameters during the simulation of 19 experiment cases which were conducted at low pressures, and the model shows good predictive performance with the remaining 88 experiment cases. However, some discrepancies in void fraction and temperature predictions are found which may require improved constitutive modeling under these low-pressure, low-flow condition. The quantified uncertainty in void fraction predictions in three forced convection benchmarks due to uncertainties in boundary conditions and geometries are relatively small, except for a few cases in the annulus experiment. The output uncertainty of the natural circulation model is quantified by considering the measurement error in boundary conditions, and it is shown that the output uncertainty ranges are relatively small for the range of conditions considered. Through this work, a clear pathway is provided for validation of RELAP7 once the code is available.

42 ENGINEERING↗
DOE OSTI2021

Topology optimization of additively manufactured fluidic components free of internal support structures

This report integrates projection-based approaches for implementing overhang constraints with fluid topology optimization to design additively manufactured fluidic components that do not require internal support structures. Internal support structures are challenging, potentially impossible to remove in complex fluid networks, yet would degrade fluidic performance if left within the part. This presents a major challenge in coupling topology optimization with additive manufacturing in the design of fluid components. The proposed approach leverages past work in overhang constraints for solid mechanics, including projection formulations and adjoint sensitivity analysis, and considers incompressible Navier–Stokes equations for laminar fluid flow. The approach is demonstrated on 2D and 3D fluid topology optimization problems to minimize pressure drop and constrain mass flow rates in pipes and manifolds. Resulting designs are crisp, logical, satisfy performance constraints, and satisfy overhang constraints, eliminating the need for internal support structures.

42 ENGINEERING↗
DOE OSTI2023

Lieb-Schultz-Mattis, Luttinger, and 't Hooft - anomaly matching in lattice systems

We analyze lattice Hamiltonian systems whose global symmetries have ’t Hooft anomalies. As is common in the study of anomalies, they are probed by coupling the system to classical background gauge fields. For flat fields (vanishing field strength), the nonzero spatial components of the gauge fields can be thought of as twisted boundary conditions, or equivalently, as topological defects. The symmetries of the twisted Hilbert space and their representations capture the anomalies. We demonstrate this approach with a number of examples. In some of them, the anomalous symmetries are internal symmetries of the lattice system, but they do not act on-site. (We clarify the notion of “on-site action.”) In other cases, the anomalous symmetries involve lattice translations. Using this approach we frame many known and new results in a unified fashion. In this work, we limit ourselves to 1+1d systems with a spatial lattice. In particular, we present a lattice system that flows to the c=1 c = 1 compact boson system with any radius (no BKT transition) with the full internal symmetry of the continuum theory, with its anomalies and its T-duality. As another application, we analyze various spin chain models and phrase their Lieb-Shultz-Mattis theorem as an ’t Hooft anomaly matching condition. We also show in what sense filling constraints like Luttinger theorem can and cannot be viewed as reflecting an anomaly. As a by-product, our understanding allows us to use information from the continuum theory to derive some exact results in lattice model of interest, such as the lattice momenta of the low-energy states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A computational parametric study of ducted fuel injection implementation in a heavy-duty diesel engine

Experiments have shown that ducted fuel injection (DFI) effectively reduces soot emissions from direct-injection diesel engines. Although many computational studies have evaluated DFI’s spray development and soot reduction mechanisms in constant volume chambers, only limited computational work on internal combustion engines exists. The DFI duct assembly changes the engine’s in-cylinder flow, spray, and combustion development. Therefore, current production engine designs might not be optimal for achieving the best engine performance with DFI. Here, this work conducted an extensive numerical study to evaluate how parameter changes affect DFI performance. The parameters include swirl ratio, piston geometry, compression ratio (CR), number of injector orifices, split injection strategy, and exhaust gas recirculation (EGR) in a heavy-duty diesel engine utilizing DFI. The combustion and soot emission data from the Sandia compression ignition optical research engine were used for model validation. Simulations showed that an increased swirl ratio resulted in more intense jet flame-piston interaction, slowing down the combustion heat release during the late combustion stage and leading to lower indicated thermal efficiency (ITE) due to higher exhaust losses. A piston-bowl design with a reentrant inner piston edge yielded the highest thermal efficiency, due to the reduced cylinder head heat transfer loss. Additional injector orifices led to higher efficiency owing to a more advanced combustion phasing. Nevertheless, the maximum pressure rise rate (MPRR) and oxides of nitrogen (NO x ) emissions also increased with the number of injector orifices due to more rapid heat release and higher combustion temperature. Implementation of a split injection strategy combined with a higher EGR rate effectively inhibited the excessive MPRR and NO x formation. In general, the study concluded that DFI is not sensitive to most parameter changes but will benefit from future parameter optimization.

33 ADVANCED PROPULSION SYSTEMS↗