SEARCH · Engineering Papers

Results for “Fast Fluid Dynamics”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Abstract Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

Microbiology↗

NASA NTRS1999

A View of the Lunar Interior Through Lunar Laser Range Analysis

Laser ranges between observatories on the Earth and retroreflectors on the Moon started in 1969 and continue to the present. Recent range accuracies are 2 cm while earliest ranges are an C, order of magnitude less certain. Four retroreflectors are ranged: three located at the Apollo 11, 14, and 15 sites and one on the Lunakhod 2 rover. Accurate analysis of the range data determines a number of lunar science parameters. The lunar interior variables include a fluid core parameter. The Lunar Laser Ranging effort is reviewed elsewhere. Many parameters are detected through their influence on rotation. Also detected are solid-body tides and accurate selenocentric reflector locations. Determined through the rotation are moment-of-inertia differences, gravitational harmonics, potential Love number, and dissipation effects due to tides and molten core. The rotation of the Moon is not at its minimum energy state; some recently active process has caused free librations. The moment differences contributed to the recent improvement of the Moon's moment of inertia from the Lunar Prospector gravity field. The Love numbers provide bulk elastic properties. Future possibilities for measurement include oblateness of the core-mantle boundary and core moment. A study of dissipation signatures in the rotation determines tidal Q vs. frequency and concludes that the Moon has a molten core. At 1 month the tidal Q is 37 and at 1 yr it is 60. The core radius is < or = 352 km for Fe and < or = 374 km for the Fe-FeS eutectic. The core detection exceeds 3x its uncertainty. The spin of the core is not aligned with the spin of the mantle and torque arises from the velocity difference at the boundary. Yoder's turbulent boundary layer theory is used to compute the radii. The present heat generation from tides and core interaction is minor compared to radiogenic heating. The heating for ancient times is more interesting. Peale and Cassen investigated lunar tidal heating while the lunar orbit expanded due to tides on Earth. Their calculations predate the measurement of Q and should be multiplied by 3.45 to match the lunar-laser-determined Love number and monthly Q. Tidal-heating computations depend on how fast the lunar orbit evolved and whether the tidal dissipation is localized. Neither is known, but plausible assumptions lead to early central region temperature increases of several hundred degrees. Most of the energy is deposited early in the Moon's history. The turbulent boundary layer theory allows a prediction of energy dissipated at the core-mantle boundary during orbit evolution. Under the assumption that the properties of the early core are the same as at present, the energy dissipated by core-mantle interaction is about the same as for tidal dissipation, but it is deposited in a smaller volume. This source of energy is capable of promoting convection in an early fluid core and driving a dynamo. This is a transient phase with a duration depending on the rate of orbit evolution. Plausible assumptions lead to a duration of a few hundred million years. Thus, the remnant magnetization of many lunar rocks is compatible with a brief global magnetic field powered by dynamical energy dissipation. Analysis of the lunar laser ranges is providing information on lunar geophysics. Future data will improve accuracies of present solution parameters, and several more interior effects should be detectable.

Williams, J. G.↗

DOE OSTI2022

Theoretical Approach for the Fast Estimation of the Turbulent Kinematic Viscosity for Internal Flows

An analytical expression for turbulent kinematic viscosity (νt), based solely on the hydraulic Reynolds number (Re), was derived and evaluated. The analytical expression is valid for the fast estimation of νt for internal, isotropic, fully developed flows. The expression was compared with experimental and simulation data for air, water, and liquid sodium, and was shown to provide reasonable values for 2100 ≤ Re ≤ 3.6 × 106 and Prandtl number (Pr) range of 0.0107 ≤ Pr ≤ 5.65. In addition, new expressions suitable for the central portion of internal flows, away from the wall, were derived for the turbulent Reynolds number (ReT), showing its relationship to Re, as well as to the ratio of νt and the molecular kinematic viscosity (ν).

42 ENGINEERING↗

DOE OSTI2023

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Halo formation from Yukawa forces in the very early Universe

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modeling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modeling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. Here, the final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Technical Assessment of Off-Gas System Technologies for Potential Use in Molten Salt Reactors

Molten salt reactors (MSRs) are a class of Generation IV advanced reactor technologies aimed to enhance and improve the safety, fuel utilization and cost-effectiveness of nuclear power generation. MSR concepts are based on using a molten salt mixture as a primary nuclear reactor coolant, while the fuel can be either directly dissolved in the coolant (i.e., fluid-fueled) or can be in a separate solid form (i.e., salt-cooled solid-fueled). MSRs can be operated in either thermal or fast neutron spectra, as simple fissile convertors, or breeder reactors (e.g., utilizing both fissile and fertile fuels), as well as accelerator-driven sub-critical reactors or transuranic element burners. Operation of MSRs will necessitate management and treatment of the off-gas streams released to the headspace of the reactor where, in some designs, a cover gas can be circulated to remove certain fission products and maintain an inert atmosphere. An inert MSR cover gas swept over or sparged into the salt will confine radionuclides emerging from the free surface, including radioactive noble and non-noble gases, aerosols, and volatile species. The elements and compounds of these species will have a wide range of thermochemical and physical properties, some of which will decay to different products in the off-gas stream. Therefore, system components used to manage and treat MSR off-gas streams will need to be robust under dynamic reactor and chemistry conditions to support adequate reactor performance by managing the composition and species in the salt coolant and cover gas space. Off-gas system components would perform as a radionuclide boundary to the reactor vessel and a primary means for preventing radionuclide release; therefore, they are of high safety significance for MSR operations. This report reviews various technologies and materials for potential implementation in the design of an MSR off-gas system per a general conceptual framework. The reviewed technologies operate based on scrubbing, capture, delay, or separation of species in the off-gas stream. These include molten hydroxide scrubbers, solid sorbents, delay off-gas systems, particle traps, and cryogenic distillation. Since aspects of an off-gas system will be design-specific, the approach is to provide general discussions on attributes of these technologies and materials, including target off-gas stream species, operational parameters and conditions of importance to decontamination factors, implementation maturity, system component monitoring, potential off-normal conditions and considerations for waste products. The information aims to support the U.S. Nuclear Regulatory Commission in their safety evaluations for MSR designs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2026

Surrogates for Valve-Controlled Pipe Flow: Accelerating Nuclear Reactor Design

Neural surrogate models are developed to replace expensive steady-state RANS CFD simulations for valve-controlled pipe flow in nuclear reactor design. Using parametric CFD data generated with MOOSE Pronghorn across a range of valve geometry and flow conditions, three approaches are compared: a POD-based reduced-order model, a structured UNet on a cylindrical grid, and unstructured models (DeepONet and BiStride MeshGraphNet) on nondimensionalized point clouds. POD achieves the highest accuracy (99%) with fast inference but requires storing all solution snapshots, while the DeepONet and BSMS-GNN both achieve ~89% accuracy at sub-second inference, with the BSMS-GNN offering superior geometric generalizability. These surrogates enable rapid ranking of candidate valve designs and can warm-start CFD solvers to accelerate convergence, supporting agentic design iteration on the Prometheus platform.

42 - ENGINEERING↗

DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

NASA NTRS1993

Laplace's equation and the Dirichlet-Neumann map in multiply connected domains

A variety of problems in material science and fluid dynamics require the solution of Laplace's equation in multiply connected domains. Integral equation methods are natural candidates for such problems, since they discretize the boundary alone, require no special effort for free boundaries, and achieve superalgebraic convergence rates on sufficiently smooth domains in two space dimensions, regardless of shape. Current integral equation methods for the Dirichlet problem, however, require the solution of M independent problems of dimension N, where M is the number of boundary components and N is the total number of points in the discretization. In this paper, we present a new boundary integral equation approach, valid for both interior and exterior problems, which requires the solution of a single linear system of dimension N + M. We solve this system by making use of an iterative method (GMRES) combined with the last multipole method for the rapid calculation of the necessary matrix vector products. For a two-dimensional system with 200 components and 100 points on each boundary, we gain a speedup of a factor of 100 from the new analytic formulation and a factor of 50 from the fast multipole method. The resulting scheme brings large scale calculations in extremely complex domains within practical reach.

Greenbaum, A.↗

DOE OSTI2025

How fast ions mitigate turbulence and enhance confinement in tokamak fusion plasmas

Along with high temperature and density, magnetic fusion requires good confinement and a degree of transport control for thermal plasmas. Meanwhile, fast ions are generated by the external heating used to raise plasma temperature and by the fusion reactions. As a result, the fusion plasmas are effectively rendered into systems with two coexisting populations of main interest — namely, the fast ions and the thermal plasma. Interestingly, several recent experiments indicate that the fast-ion population can improve the confinement of the thermal plasmas by mitigating turbulence. Here, in this Review, we describe the physical mechanisms that underpin the improved confinement and discuss recent experimental results in terms of these mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Design and Production for Maximum Structural Efficiency With Respect to Fiber Orientation With Increased Understanding of Hybrid Fiber Flow Behavior

Discontinuous fiber-reinforced thermoplastic composites have gained considerable attention in automotive, aerospace, and other industries, due to their high-rate of production combined with their ability to attain complex and intricate shapes. Among other high-rate thermoplastic manufacturing processes, injection-molding is one of the most common manufacturing methods due to fast production and high surface finishing of complex geometries. Fiber orientation in discontinuous fiber composites plays a pivotal role in determining the mechanical, electrical, and thermomechanical properties, underscoring the necessity to comprehend fiber orientation in injection molded parts. Among different fiber types, glass and carbon fibers are most common in the composite industries. The recent trend of hybrid composites comprising both glass fiber (GF) and carbon fiber (CF) is also gaining importance in the automotive industry. Hybrid fiber options allow designers to optimize the balance between glass and carbon fibers by leveraging the high durability and low cost of GF while the strength and lightweight properties of CF. Consequently, comparing the fiber oriented distribution (FOD) of injection molded composites containing GF, CF, and a hybrid of GF/CF is critical to investigating the local mechanical properties of intricate structures for high-end applications. In Phase I of this project, FOD in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for GF, CF, and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varies with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Thermal Scattering Law Data Development for Paraffin Wax

Paraffin wax is often used as a nuclear moderator to slow down the fast neutrons in experimental critical assemblies [1]. It is a colorless and soft solid material that consists primarily of straight-chain alkanes (n-alkanes), which are hydrocarbons with the general formula CnH2n+2 [2-3]. The length of the hydrocarbon chain ranges from C20 to C30 and higher [2]. It is distinguished by its solid state at room temperature and begins to melt above approximately 310 K [4]. Paraffin wax is a commonly employed substance in the manufacture of shielding. One of its noteworthy characteristics is its ability to effectively absorb the neutrons. Also, it possesses a high macroscopic cross section, which enables it to efficiently moderate neutrons. As a result, paraffin wax is extensively utilized in various applications where moderation and shielding of neutrons are needed. For simulations, it is necessary to evaluate its thermal scattering law (TSL) and cross sections. Computationally, classical molecular dynamics (CMD) simulations provide the capability of simulating atomic details. For example, several unary, binary, and few multi component mixtures have been investigated of the paraffin model by using molecular dynamics simulations [5-12]. An assessment of thermal neutron scattering in a heavy paraffinic oil treated both as a solid and a viscous fluid containing 25% linear branched paraffin (C30H62), 35% one ring cycloalkane (C30H60), 15% two rings cycloalkane (C30H58), and 25% aromatic (C30H60) chains has been studied using CMD simulations for producing TSL data [13]. Nevertheless, there is lack of TSL and cross section data for paraffin wax as most of the reported analyses focus on the unary and binary mixture of n-alkanes, which is not consistent with actual paraffin wax [2]. In this work, we applied the equilibrium CMD simulations technique to explore the structure and dynamical properties of wax, which are fundamental input to calculate the TSL. A paraffin wax system was modeled using the CMD code LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) [14-15] with the semi-empirical COMPASS [16] force field. The density of state (DOS) was calculated from the normalized velocity autocorrelation function (VACF), which is the Fourier transform of the normalized VACF. The DOS was used for the calculation of the TSL and thermal scattering cross sections. The paraffin wax atomic system was constructed by using the MedeA material design platform [17], and was benchmarked using available properties (i.e., density, bond lengths, angles, diffusivity, and viscosity).

Nuclear Criticality Safety Program (NCSP)↗

DOE OSTI2023

Stream lamination and rapid mixing in a microfluidic jet for X-ray spectroscopy studies

Microfluidic mixers offer new possibilities for the study of fast reaction kinetics down to the microsecond time scale, and methods such as soft X-ray absorption spectroscopy are powerful analysis techniques. These systems impose challenging constraints on mixing time scales, sample volume, detection region size and component materials. The current work presents a novel micromixer and jet device which aims to address these limitations. The system uses a so-called ‘theta’ mixer consisting of two sintered and fused glass capillaries. Sample and carrier fluids are injected separately into the inlets of the adjacent capillaries. At the downstream end, the two streams exit two micron-scale adjoining nozzles and form a single free-standing jet. The flow-rate difference between the two streams results in the rapid acceleration and lamination of the sample stream. This creates a small transverse dimension and induces diffusive mixing of the sample and carrier stream solutions within a time scale of 0.9 microseconds. The reaction occurs at or very near a free surface so that reactants and products are more directly accessible to interrogation using soft X-ray. We use a simple diffusion model and quantitative measurements of fluorescence quenching (of fluorescein with potassium iodide) to characterize the mixing dynamics across flow-rate ratios.

Huyke, Diego A. (ORCID:0000000283356613)↗

DOE OSTI2026

Reduced-order CFD modeling of cryogenic hydrogen isotope extrusion for pellet fueling

This study presents a reduced-order model (ROM) for computational fluid dynamics (CFD) simulations of cryogenic hydrogen isotope extrusions, focusing on protium (H₂) and deuterium (D₂) piston extruders. Using a 2D axisymmetric ROM in ANSYS-Polyflow, significant computational savings were achieved (runtime reduced from 9∼24 h to 3∼5 min), with extrusion force discrepancies between the 2D ROM and 3D models being on the order of 1%. Parametric studies identified optimal cutoff shear rates in the viscosity model (0.01/s for H₂ and 0.001/s for D₂), providing recommendations for future simulations. Finally, a comprehensive comparison of ROM results with experimental data was performed across varying geometries, cryogenic materials, temperatures, extrusion lengths, and piston velocities. Predictions at low extrusion temperatures met the objective of providing quick and efficient solutions with an acceptable extrusion force error of approximately 10% or less, validating the effectiveness of the 2D ROM approach. However, at high temperatures closer to the triple point, extrusion force error grows, which necessitates developing an improved model that accounts for temperature effects, e.g. melting. Nevertheless, the findings still represent a significant improvement in efficiency of CFD modeling of cryogenic hydrogenic extrusion. The ROM framework can also be extended to tritium (T2) and screw extruders, which will ultimately provide a fast and effective tool for optimizing pellet injector design for ITER and future reactor systems.

Fan, Joy [ORNL] (ORCID:0000000229751735)↗

NASA NTRS1999

Experimental Studies of the Interaction Between a Parallel Shear Flow and a Directionally-Solidifying Front

It has long been recognized that flow in the melt can have a profound influence on the dynamics of a solidifying interface and hence the quality of the solid material. In particular, flow affects the heat and mass transfer, and causes spatial and temporal variations in the flow and melt composition. This results in a crystal with nonuniform physical properties. Flow can be generated by buoyancy, expansion or contraction upon phase change, and thermo-soluto capillary effects. In general, these flows can not be avoided and can have an adverse effect on the stability of the crystal structures. This motivates crystal growth experiments in a microgravity environment, where buoyancy-driven convection is significantly suppressed. However, transient accelerations (g-jitter) caused by the acceleration of the spacecraft can affect the melt, while convection generated from the effects other than buoyancy remain important. Rather than bemoan the presence of convection as a source of interfacial instability, Hurle in the 1960s suggested that flow in the melt, either forced or natural convection, might be used to stabilize the interface. Delves considered the imposition of both a parabolic velocity profile and a Blasius boundary layer flow over the interface. He concluded that fast stirring could stabilize the interface to perturbations whose wave vector is in the direction of the fluid velocity. Forth and Wheeler considered the effect of the asymptotic suction boundary layer profile. They showed that the effect of the shear flow was to generate travelling waves parallel to the flow with a speed proportional to the Reynolds number. There have been few quantitative, experimental works reporting on the coupling effect of fluid flow and morphological instabilities. Huang studied plane Couette flow over cells and dendrites. It was found that this flow could greatly enhance the planar stability and even induce the cell-planar transition. A rotating impeller was buried inside the sample cell, driven by an outside rotating magnet, in order to generate the flow. However, it appears that this was not a well-controlled flow and may also have been unsteady. In the present experimental study, we want to study how a forced parallel shear flow in a Hele-Shaw cell interacts with the directionally solidifying crystal interface. The comparison of experimental data show that the parallel shear flow in a Hele-Shaw cell has a strong stabilizing effect on the planar interface by damping the existing initial perturbations. The flow also shows a stabilizing effect on the cellular interface by slightly reducing the exponential growth rate of cells. The left-right symmetry of cells is broken by the flow with cells tilting toward the incoming flow direction. The tilting angle increases with the velocity ratio. The experimental results are explained through the parallel flow effect on lateral solute transport. The phenomenon of cells tilting against the flow is consistent with the numerical result of Dantzig and Chao.

Zhang, Meng↗

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↗

NASA NTRS1994

Inducer analysis/pump model development

Current design of high performance turbopumps for rocket engines requires effective and robust analytical tools to provide design information in a productive manner. The main goal of this study was to develop a robust and effective computational fluid dynamics (CFD) pump model for general turbopump design and analysis applications. A finite difference Navier-Stokes flow solver, FDNS, which includes an extended k-epsilon turbulence model and appropriate moving zonal interface boundary conditions, was developed to analyze turbulent flows in turbomachinery devices. In the present study, three key components of the turbopump, the inducer, impeller, and diffuser, were investigated by the proposed pump model, and the numerical results were benchmarked by the experimental data provided by Rocketdyne. For the numerical calculation of inducer flows with tip clearance, the turbulence model and grid spacing are very important. Meanwhile, the development of the cross-stream secondary flow, generated by curved blade passage and the flow through tip leakage, has a strong effect on the inducer flow. Hence, the prediction of the inducer performance critically depends on whether the numerical scheme of the pump model can simulate the secondary flow pattern accurately or not. The impeller and diffuser, however, are dominated by pressure-driven flows such that the effects of turbulence model and grid spacing (except near leading and trailing edges of blades) are less sensitive. The present CFD pump model has been proved to be an efficient and robust analytical tool for pump design due to its very compact numerical structure (requiring small memory), fast turnaround computing time, and versatility for different geometries.

Cheng, Gary C.↗

DOE OSTI2020

Comparison of Model Predictions and Performance Test Data for a Prototype Thermal Energy Storage Module

Although model predictions of thermal energy storage (TES) performance have been explored in previous investigations, relevant test data that enable experimental validation of performance models have been limited. This is particularly true for high-performance TES designs that facilitate fast input and extraction of energy. In this paper, we present a summary of experimental tests of a high-performance TES unit using lithium nitrate trihydrate phase change material as a storage medium. Performance data are presented for complete dual-mode cycles consisting of extraction (melting) followed by charging (freezing). These tests simulate the cyclic operation of a TES unit for asynchronous cooling in a variety of applications. Finally, the model analysis is found to agree reasonably well, within 10%, with the experimental data except for conditions very near the initiation of freezing, a consequence of subcooling that is required to initiate solidification.

25 ENERGY STORAGE↗