Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Properties”

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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Piston geometry and stroke optimization for high efficiency propane spark ignition engines

Propane has unique properties and offers interesting characteristics for high-efficiency spark ignition engines. Its high volatility reduces or completely eliminates fuel-wall wetting and facilitates fuel air mixing. Furthermore, propane has a research octane number of 112 and a high octane sensitivity of 15. Finally, its laminar flame speed is on the same order as that of conventional gasoline, and it exhibits high dilution tolerance. Modern spark ignition internal combustion engines rely on fast combustion rates and high dilution to achieve high brake thermal efficiencies. To accomplish this, high stroke-to-bore ratios and high geometric compression ratios have been used in new engine designs. Therefore, propane’s relatively high laminar flame speeds, high knock resistance, and dilution tolerance make it an excellent candidate fuel for modern spark ignition engines. The objective of this work is to co-optimize the piston geometry and the engine stroke to maximize the efficiency of a spark-ignition engine fueled with propane. 3D computational fluid dynamics (CFD) simulations employing the extended coherent flamelet model were used to study the parametric effects of piston shape and stroke length. A piston geometry based on high performing pistons was parameterized using four controlling parameters. The piston geometry and engine stroke design space was explored using deterministic and quasi-random sampling techniques. In conclusion, a Gaussian process regression model was built using the simulation data to explain the results observed.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Hydraulic Modeling of Induced and Propagated Fractures: Analysis of Flow and Pressure Data From Hydromechanical Experiments in the COSC-1 Deep Borehole in Crystalline Rock Near Åre, Sweden

To characterize the coupled hydromechanical behavior of rock fractures, the step-rate injection method for fracture in-situ properties (SIMFIP) was conducted with a specialized downhole probe developed by Guglielmi et.al. In June 2019, a field campaign was carried out near Åre, Sweden, where the SIMFIP probe was applied in the Collisional Orogeny in the Scandinavian Caledonides-1 scientific borehole to understand the dynamics of injection-induced fracture initiation, fracture opening, and shearing due to water injection-withdrawal in a borehole interval isolated by two packers. Three intervals were investigated at ~500 m depth: (a) an unfractured section (intact rock), (b) a section with non-conductive fractures, and (c) a section with hydraulically conductive fractures. Pressure, injection flow rate, and borehole wall displacement were simultaneously measured during the tests. In the present study, the geometry of the induced fracture and deformation of existing fractures at different time stages of the tests are determined based on a hydrologic model by using the measured pressure and flow data during each time stage of the experiment. A numerical model for the fluid flow within the fracture and the packed-off borehole interval is implemented within COMSOL Multiphysics. By matching model simulations with observed data for all three sections, estimates of the induced and propagated fractures' radius and aperture at successive time stages have been obtained in each case. We could also determine the non-linear relationship between fracture aperture and pressure for values above fracture opening pressures. The model results provide insights for the understanding of pressure-induced fracture initiation and propagation in crystalline rock.

58 GEOSCIENCES↗
DOE OSTI2025

Efficient machine learning interatomic potentials robust for liquid and multiple solid polymorphs of NaF and KF

Achieving atomic-level understanding of crystallization of molten salts is of importance to a wide range of technological applications. Recent work [Fan et al., Proc. Natl. Acad. Sci. USA 122, e2425702122 (2025)] revealed that crystal nucleation in molten LiF salt is a multistage process according to the molecular-dynamics (MD) simulations based on an atomic cluster expansion (ACE) machine-learning interatomic potential (MLIP). In order to understand the influence of increasing cation size on nucleation pathways and nucleation rates of molten fluoride salts, here we develop two new ACE MLIPs for NaF and KF. The two ACE MLIPs feature DFT-SCAN-level accuracy for liquid and multiple solid polymorphs over a wide temperature (0–2000 K) and pressure (0–100 GPa) range, and also reproduce well a number of experimental data for solid and liquid equilibrium properties. The efficiency of the two ACE MLIPs enable million-atom-scale or microsecond-scale MD simulations. The two general-purpose ACE MLIPs are expected to be useful for atomistic simulations for different purposes, in addition to studying crystallization of molten NaF and KF salts.

Crystal melting↗
DOE OSTI2023

Physics-Driven Convolutional Autoencoder Approach for CFD Data Compressions: Preprint

With the growing size and complexity of turbulent flow models, data compression approaches are of the utmost importance to analyze, visualize, or restart the simulations. Recently, in-situ autoencoder-based compression approaches have been proposed and shown to be effective at producing reduced representations of turbulent flow data. However, these approaches focus solely on training the model using point-wise sample reconstruction losses that do not take advantage of the physical properties of turbulent flows. In this paper, we show that training autoencoders with additional physics-informed regularizations, e.g., enforcing incompressibility and preserving enstrophy, improves the compression model in three ways: (i) the compressed data better conform to known physics for homogeneous isotropic turbulence without negatively impacting point-wise reconstruction quality, (ii) inspection of the gradients of the trained model uncovers changes to the learned compression mapping that can facilitate the use of explainability techniques, and (iii) as a performance byproduct, training losses are shown to converge up to 12x faster than the baseline model.

auto-encoders↗
DOE OSTI2023

Design optimization of an ethanol heavy-duty engine using design of experiments and bayesian optimization

Diesel-fueled engines still hold a large market share in the medium and heavy-duty transportation sector. However, the increase in fossil fuel prices and the strict emission regulations are leading engine manufacturers to seek cleaner alternatives without a compromise in performance. Alcohol-based fuels, such as ethanol, offer a promising alternative to diesel fuel in meeting regulatory demands. Ethanol provides cleaner combustion and lower levels of soot due to its chemical properties, in particular its lower level of carbon content. In addition, the stoichiometric operating conditions of alcohol fueled engines enable the mitigation of NOx emissions in aftertreatment stage. With the promise of retrofitting diesel engines to run on ethanol to reduce emissions, the thermal efficiency of these engines remains the primary optimization target. In order to find the optimal ethanol-fueled engine design that maximizes the thermal efficiency, a large design space needs to be investigated using engineering tools. In this study, previous research by the authors on optimizing the design of a single-cylinder ethanol-fueled engine was extended to explore the design space for a heavy-duty multi-cylinder engine configuration. A heavy-duty engine setup with multiple operating conditions at different engine speeds and loads were considered. A design optimization analysis was performed to identify the potential designs that maximize the indicated thermal efficiency in an ethanol-fueled compression ignition engine. First, a computational fluid dynamics (CFD) model of the engine was validated using experimental data for four drive cycle points. Using a design of experiments (DoE) approach and a parameterized piston bowl geometry, the model was then exercised to explore the relationship among geometric features of the piston bowl and spray targeting angle and indicated thermal efficiency across all tested operating conditions. After evaluating 165~candidate designs, a piston bowl geometry was identified that yielded an increase between 1.3 to 2.2 percentage points in indicated thermal efficiency for all tested conditions, while satisfying the operational design constraints for peak pressure and maximum pressure rise rate. The increased performance was attributed to enhanced mixing that led to the formation of a more homogeneous distribution of in-cylinder temperature and equivalence ratio, higher combustion temperatures, and shorter combustion duration. Finally, a Bayesian optimization (BOpt) analysis was employed to find the optimal piston bowl geometry with a fixed spray injector angle for one of the operating conditions. Using BOpt, a piston candidate was identified that resulted in a 1.9~percentage point increase in thermal efficiency from the baseline design, yet only required 65\% of the design samples investigated using the DoE approach.

Tekgul, Bulut↗
DOE OSTI2020

Gaseous jet through an outward opening injector: Details of mixing characteristic and turbulence scales

Direct injection (DI) strategy of natural gas (NG) into internal combustion engines (ICE) has led to higher thermal efficiency and lower exhaust emissions. In order to thoroughly understand the most relevant phenomena affecting the performances of such engines, computational fluid dynamics (CFD) plays a key role as an accurate description of the jet evolution and interaction within the combustion chamber is required to that aim. Accurate description of high-pressure gaseous jets is rather challenging at high Mach numbers, as the injected gas is strongly under-expanded once in the ambient, giving room to shocks due to compressibility effects. Also the interaction between shock waves and mixing layers needs to be carefully represented with a multi-dimensional model, calling for substantial computational resources requirements. In this paper a numerical investigation of the behavior of a gaseous jet (Argon) through an outward opening injector has been carried out. A Large Eddy Simulation (LES) approach has been used in order to track the structures derived by the interaction of the injected fuel with the surrounding ambient. Although already good results were obtained using a Reynolds Averaged Navier-Stokes (RANS) approach, the adoption of LES is required to characterize more accurately the jet properties in terms of vortex structures and mixing effectiveness. Finally, the effect of the Nozzle Pressure Ratio (NPR) on the jet evolution has been highlighted in the paper, showing how a higher NPR would give a faster injection process, compromising however the homogeneity of the mixture.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Electrical Conductivity of Clayey Rocks and Soils: A Non-Linear Model

It is well-accepted that Archie's law is only applicable to “clean” rocks and soils but fails in “dirty” ones where clay minerals possess an additional component of surface conductivity. Although several models, for example, Waxman-Smits model, were presented to account for this phenomenon, surface conductivity is always inappropriately treated as constant, which actually only holds at high salinities. The essential non-linear characteristic differing over fluid salinities has not been physically or mathematically explained well in those models. Here, we scrutinize the conduction mechanism of clayey rocks and soils and ascribe this non-linear feature to (a) variation of the electrical double layer and (b) the intrinsic clay-and-water conduction pattern. With effective medium theory, we develop an easy-to-use non-linear model that both reflects electrochemical theories and explains the measurement data well. Our model can be used to produce more accurate results for laboratory- and field-scale petrophysical parameter evaluations than the previous models.

58 GEOSCIENCES↗
DOE OSTI2023

Durable and Low-Cost Fractal Structured Multifunctional Coatings for Next-Generation CSP

Fractal coatings on various low-cost and high-temperature metallic substrates such as SS316, SS347, In800H, In740H, In625, Haynes 230, etc. were fabricated employing our patented novel electrodeposition process at various deposition potentials and deposition times, as well as chemical etching. The best combination of the materials and processing conditions were co-optimized for the reliable and repeatable fabrication of the multifunctional coatings. Surface morphologies were examined using SEM and AFM, and fractal characterization was done using a 3D optical surface profiler. Fractal dimension was used to identify successful fabrication and material parameters. The mechanical durability of the coatings was evaluated using several standardized tests, such as sand abrasion (ASTM D968) and adhesion (ASTM D3359), to fully characterize and demonstrate the mechanical durability of the fabricated surfaces. HTF immersion tests on the fractal multifunctional coatings were conducted based on a modification of the ASTM D870, which is established for static water immersion testing. To calculate the corrosion rates and descaled weight loss of the samples subjected to thermal exposure for different time durations (100 h, 300 h, 500 h, and optionally 750 h) to estimate the corrosion rate. While protecting the substrates from corrosion is an important goal, an allied goals is also to examine the post-corrosion temperature stability and properties of the molten salts themselves due to the dissolution of the constituents of the alloys and the corrosion products from the coatings that cause “fouling” of the molten salts. Fouling characteristics of the different molten salt HTFs exposed to coated and uncoated substrates were investigated. Finally, considering a prototype molten salt/sCO2 heat exchanger, cost/performance model is developed to demonstrate that the coatings on low-cost alloys such as stainless steel provide lower levelized cost of heat exchanger compared to high nickel content Ha230.

14 SOLAR ENERGY↗
DOE OSTI2020

Small-angle Neutron Scattering (SANS) Characterization of Clay- and Carbonate-rich Shale at Elevated Pressures

Unconventional oil and gas from shale formations have emerged as some of the fastest growing energy resources in the United States, providing both cleaner energy to consumers and reducing the nation’s reliance on energy imports. To properly harness these important natural resources, the nanopore structure of associated shales must be fully understood, particularly under hydraulic fracturing conditions, where they are exposed to both overburden compressive and hydrostatic fluid pressures. The current study uses small-angle neutron scattering (SANS) to characterize pore structure, including porosity, pore accessibility, and pore size distribution, in the 1–100 nm regime at elevated pressures for mineralogically distinct clay- and carbonate-rich shales from the Permian Basin. Unlike typical porosity measurement techniques, SANS is uniquely capable of characterizing both open and closed porosity, allowing measurement of how pore accessibility changes with pressure and determination of the size range of accessible versus inaccessible pores. The porosity of the clay-rich shale was 7.7%, compared to 0.51% for the carbonate-rich shale. However, only 2.6% of the nanopores in the carbonate-rich shale were inaccessible to water at 8 kPSI (55.1 MPa) compared to 7.8% for the clay-rich shale. Furthermore, our analyses indicated that the closed pores fall within distinct size ranges, likely corresponding with the chemical nature of the pore host material. These results provide valuable insight into the effects of shale petrophysical properties on hydrocarbon extraction from unconventional reservoirs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Numerical Investigation of the Spark Discharge Process in a Crossflow

The present study numerically investigates the spark discharge process under crossflow conditions using a thermal equilibrium plasma solver that fully couples the electromagnetic physics and fluid dynamics in a computational framework. Numerical results are validated by the comparison with experimental data. Here, the spark discharge experiment is performed in a constant volume vessel using an inductive coil ignition system for automotive applications, and the evolution of the spark channel is measured using high-speed imaging. The crossflow in the gap between the spark-plug electrodes is generated by a rotating fan with two different fan speeds, and the flow velocity across the gap is characterized by particle image velocimetry (PIV) measurement. A computational fluid dynamics (CFD) solver is employed to simulate the crossflow and provide the flow field variables (velocity, pressure, temperature) to the plasma solver. The crossflow velocity predicted in the flow simulation agrees well with the PIV data in that the non-uniform velocity profiles at monitoring points are reproduced by the CFD code. With the crossflow initialization in the plasma solver, the simulated spark discharge process from the breakdown to spark discharge matches the experimental data, including the voltage and circuit waveforms and the high-speed images of the spark channel evolution. The stretch of spark channel captured by plasma simulations agrees with the measured data. The plasma simulation reveals that the mean temperature of the spark channel is maintained at 5000 K during the discharge phase, and the temperature varies along the spark channel so that the highest value is obtained at the spark root on the center electrode. Overall, the results presented in this paper are meant to provide valuable information about the properties of the plasma generated by the spark discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Preliminary modeling of triply periodic minimal surface (TPMS) structures using RELAP5-3D

With the United States Department of Energy (DOE)’s goal of quadrupling the nation’s nuclear energy supply by 2050, and with the Advanced Fuels Campaign pushing for new types of advanced reactor fuels and geometries, the need has arisen for new nuclear fuel designs. One such design is to swap out current nuclear fuel geometries in exchange for another type of geometry, called a Triply Periodic Minimal Surface (TPMS). TPMSs are self-supporting, infinitely repeating lattices—attributes that lend themselves well to additive manufacturing. These surfaces also possess enhanced heat transfer properties thanks to their internal area changes and large surface-area-to-volume ratios. Their drawback, however, is an increased pressure drop. Given the small amount of correlations and data (Reynolds numbers in the 2,000–8,000 range), and the minimal amount of experience so far obtained by modeling TPMS structures using 1D systems codes such as the Reactor Excursion and Leak Analysis Program (RELAP5-3D), further research into this topic was needed. Using data from the University of Wisconsin - Madison (UW), curve fits were created for both a Heat Transfer Coefficient (HTC) correlation and a Darcy friction factor empirical coefficient correlation. The curves’ coefficients and multipliers were then output and utilized in RELAP5-3D models of two upcoming experiments—Flow Loop for INFLUX Pressure drop (FLIP) and Microreactor Agile Non-nuclear Experimental Test (MAGNET)—aimed at increasing the available data for Reynolds numbers to the 16,000–36,000 range for TPMS structures. The models were run under the conditions utilized by a Computational Fluid Dynamics (CFD) analysis performed by another group at Idaho National Laboratory. Only CFD pressure drop values were obtained from the FLIP test, and those values showed that the RELAP5-3D models had a lower rate of pressure increase in comparison to the CFD values. In addition, there seemed to be a vertical shift upward in the pressure drop for both models whenever the TPMS porosity decreased, and the RELAP5-3D models showed a higher vertical shift in comparison to the CFD values. The MAGNET results did not correspond to any CFD or experimental results against which they could be compared, so they were instead compared against the proposed CFD input conditions. These values were then compared with each other to make sure the model seemed to be performing as expected, paving the way for future tests that can be run for the purpose of further analyses and comparisons. The pressure drop increased with temperature and mass flow rate independently. The temperature change would decrease with increasing mass flow rate and temperature, which was just as we expected based on the fact that the lower viscosity and decreased density would result in higher friction and churning losses. The last metric that was assessed was the enthalpy flow change, which increased with increasing mass flow rate and decreasing temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2027

A multi-region approach for the analysis of porous materials and gas-solid reactions using USAXS-SAXS-WAXS: application to CaO carbonation

The reaction rate of gas-solid non-catalytic reactions is typically investigated using reactant conversion data over time and ex-situ measurements of the porous solid reactant textural properties; these data enable the experimental estimation of the initial intrinsic reaction rate, whereas the evolution of the textural properties and the reaction rate over time are evaluated theoretically using reaction models. In this work, a different methodology is presented, based on: a) in-situ time-resolved USAXS-SAXS-WAXS measurements of the solid sample, and b) a multi-region modeling approach; this methodology allows for the estimation of the textural properties and the intrinsic reaction rate at any time during the reaction. Referring to a reaction in which a solid product is obtained from a solid reactant and a gaseous reactant, the porous particle is described as consisting of several distinct regions with different microstructural properties, and both a gas-solid reaction model and a SAXS model are derived and applied to the carbonation of porous CaO. The proposed reaction model highlights the role of the inaccessible reactant and accurately predicts the solid reactant conversion versus time profile. The SAXS model accurately predicts the measured linear trends of the Porod invariant and of the pre-factor of the power law scattering profile versus CaO mass fraction; in the absence of chemical reactions, the proposed equations extend classical SAXS theory to porous materials containing macropores and nonporous solid phases in addition to standard nanoscale inhomogeneous regions.

CO2 solid sorbents↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Coupling Thermal-Hydraulics with Reactor Pressure Vessel Fracture Models

Because of the importance of maintaining the integrity of the reactor pressure vessel (RPV) in light-water reactors, significant effort has been devoted over the last several decades to understand the potential failure mechanisms, characterize how material properties evolve due to environmental conditions, and develop tools to assess the probability of fracture in aged RPVs during potential transient events. One important aspect of this problem that has not yet been fully addressed is the effects of nonuniform cooling on the RPV inner wall due to the way low-temperature water injected into the primary coolant system is distributed within the RPV during potential loss-of-coolant accident transients. Colder temperatures are expected in the regions near the inlets, causing a “cold-plume” effect. Rapid cooling of the RPV wall increases the likelihood of fracture initiation, and localized cooling due to cold-plume effects could potentially exacerbate this problem. The Grizzly code can perform a three-dimensional simulation of the RPV thermomechanical response as well as the probabilistic assessment of fracture for a population of pre-existing flaws, which is necessary to account for spatial variations in coolant temperature. This report documents a first effort to couple Grizzly RPV fracture models with computational fluid dynamics models performed using the Cardinal code, which allows for a realistic representation of the temperature distribution. This report demonstrates this coupling approach on two test cases: a simplified proof-of-concept scenario and a realistic small-break loss-of-coolant accident scenario. In both cases, the effects of the spatially varying temperature are evident; although, because the small-break loss-of-coolant accident scenario was far from challenging the RPV’s margins, it is still inconclusive how important cold-plume effects are in realistic accident scenarios. However, this effort represents an important step toward performing realistic analyses to determine to what extent spatially nonuniform cooling might affect RPV integrity under more aggressive accident scenarios.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

HYDRODYNAMIC ANALYSIS OF MACROALGAE LOCAL MODEL USING COMPUTATIONAL FLUID DYNAMICS

In this article, a local scale, fully nonlinear coupled fluid-structural interaction (FSI) sugar kelp model has been developed using a computational fluid dynamics (CFD) method. In this model, to be consistent with available experimental data, the sugar kelp is approximated as elongated rectangles with smoothed isosceles triangles at the ends and a single kelp model with one end fixed in a channel with constant current model is developed. Several different current speeds are simulated, and the resulting drag forces and calculated drag coefficients are validated by comparison with experimental data from the literature. In a previous study, a global scale model was developed using a computational structural dynamics (CSD) method to simulate macroalgae farming system and guide the system configuration design. In the global scale model, the hydrodynamic forces are calculated using Morison’s equation and the kinematics and dynamics of the sugar kelp are simplified and the group of kelps attached to the long line is modeled as a slender structure with the same length and an effective diameter such that the volumes are consistent with the real physical system. This simplified model matches the weight and buoyancy but adjusting the hydrodynamic properties when the general hydrodynamic coefficients are employed. Therefore, optimal hydrodynamic coefficients used in global scale model were determined to obtain the hydrodynamic force more accurately. The validated local scale model is then be applied to determine the hydrodynamic coefficients of the simplified sugar kelp model for global dynamic analysis.

hydrodynamic modeling, macroalgae, CFD model↗
DOE OSTI2022

Turbulent flow characteristics in an 84-pin rod bundle for typical and damaged spacer grids

Hexagonal rod bundles arranged in a tightly packed triangular lattice are extensively used for heat transfer and energy generation applications. Staggered spacer grids are used to maintain the structural integrity of gas-cooled fast reactor (GFR) fuel assemblies, while inducing localized turbulence in flow. Damage to these spacer grids results in a disruption of flow fields within these hexagonal fuel bundles. Experimental flow visualizations are critical to identify the differences in local flow properties that the structural damage may cause. This experimental research investigates the flow-field characteristics at a near-wall and center plane in a prototypical 84-pin GFR fuel assembly. Newly installed typical spacers and spacers subject to naturally occurring damage due to material degradation over prolonged experimentation were investigated. Velocity fields were acquired by utilizing the matched-index-of-refraction method to obtain time-resolved particle image velocimetry measurements for a Reynolds number of 12 000. Reynolds decomposition statistical results divulged differences in the time-averaged velocity, velocity fluctuations, flow anisotropy, and Reynolds stress distributions. Galilean decomposition demarcated the influence of spacer grid damage on the velocity fields. To extract turbulent structures and elucidate mechanisms of flow instabilities, proper orthogonal decomposition analysis was employed. Reduced order flow reconstructions enabled the application of vortex identification algorithms to determine the spatial and statistical characteristics of vortices generated. This research work provides unique experimental data on the spacer grid condition-dependent flow. The results offer a deeper understanding of fluid dynamics behavior to support GFR rod bundle design efforts and computational fluid dynamics model validation.

36 MATERIALS SCIENCE↗