Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fluid drag”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Analysis of convection heat transfer on multiscale rough superhydrophobic and liquid infused surfaces

Multiscale rough superhydrophobic or slippery liquid infused porous surfaces have gained much interest in recent years for their improved transport phenomena properties. While there have been several studies on drag reduction and condensation on non-wetting surfaces, convection heat transfer that is important in many thermal and thermochemical applications has not been addressed systematically. Here, this article utilizes a fractal description of rough surface topographies to develop analytical models for the Nusselt number and the thermal hydraulic factor for fluid flow and heat transfer inside a cylinder with non-wetting surfaces. For air-infused superhydrophobic surfaces, the model considers the dynamic stability of the air/fluid interface in the asperities. Using the analytical formulations and the stability criteria, systematic studies are presented on the effects of the fractal surface parameters, cylinder radius and Reynolds number on the convective heat transfer characteristics, from which surface texture design maps are developed for maximizing the convection heat transfer. It is shown that multiscale non-wetting surfaces are most effective in the range of lower Reynolds number and small cylinder radius for achieving the best convective heat transfer and thermal hydraulic performance. Applying the models to actual non-wetting surface topographies fabricated using electrodeposition and chemical etching, it is shown that contrary to prevailing notion, superhydrophobicity, characterized by the highest contact angles, does not always lead to the maximum convective heat transfer performance, and that under certain fluid flow conditions, hydrophobic surfaces may offer a greater thermal performance.

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

Magnetohydrodynamics in free surface liquid metal flow relevant to plasma-facing components

Abstract While flowing Liquid Metal (LM) Plasma-Facing Components (PFCs) represent a potentially transformative technology to enable long-pulse operation with high-power exhaust for fusion reactors, Magnetohydrodynamic (MHD) drag in the conducting LM will reduce the flow speed. Experiments have been completed in the linear open-channel LMX-U device [Hvasta et al 2018 Nucl. Fusion 58 01602] for validation of MHD drag calculations with either insulating or conducting walls, with codes similar to those used to design flowing LM PFCs for a Fusion Nuclear Science Facility [Kessel et al 2019 Fusion Sci. Technol . 75 886]. We observe that the average channel flow speed decreased with the use of conducting walls and the strength of the applied transverse magnetic field. The MHD drag from the retarding Lorentz force resulted in an increase of the LM depth in the channel that ‘piled up’ near the inlet, but not the outlet. As reproduced by OpenFOAM and ANSYS CFX calculations, the magnitude and characteristics of the pileup in the flow direction increased with the applied traverse magnetic field by up to 120%, as compared to the case without an applied magnetic field, corresponding to an average velocity reduction of ∼45%. Particle tracking measurements confirmed a predicted shear in the flow speed, with the surface velocity increasing by 300%, despite the 45% drop in the average bulk speed. The MHD effect makes the bulk flow laminarized but keeps surface waves aligned along the magnetic field lines due to the anisotropy of MHD drag. The 3D fringe field and high surface velocity generate ripples around the outlet region. It was also confirmed that the MHD drag strongly depends on the conductivity of the channel walls, magnetic field, and volumetric flow rate, in agreement with the simulations and a developed analytical model. These validated models are now available to begin to determine the conditions under which the ideal LM channel design of a constant flow speed and fluid depth could be attained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Perfect Coulomb drag and exciton transport in an excitonic insulator

Strongly coupled electron-hole bilayers can host quantum states of interlayer excitons, such as high-temperature exciton condensates at zero magnetic field. This state is predicted to feature perfect Coulomb drag, where a current in one layer is accompanied by an equal but opposite current in the other. We used an optical technique to probe the electrical transport of correlated electron-hole bilayers based on MoSe2/hBN/WSe2 heterostructures. We observed perfect Coulomb drag in the excitonic insulator phase at low temperatures; the counterflow resistance of interlayer excitons remained finite. These results indicate the formation of an exciton gas that does not condense into a superfluid. Our work demonstrates that dynamic optical spectroscopy provides a powerful tool for probing exciton transport behavior in correlated electron-hole fluids.

Qi, Ruishi↗
DOE OSTI2024

Hard jet substructure in a multistage approach

We present predictions and postdictions for a wide variety of hard jet-substructure observables using a multistage model within the framework. The details of the multistage model and the various parameter choices are described in []. A novel feature of this model is the presence of two stages of jet modification: a high-virtuality phase [modeled using the modular all twist transverse-scattering elastic-drag and radiation model ()], where modified coherence effects diminish medium-induced radiation, and a lower virtuality phase [modeled using the linear Boltzmann transport model ()], where parton splits are fully resolved by the medium as they endure multiple scattering induced energy loss. Energy-loss calculations are carried out on event-by-event viscous fluid dynamic backgrounds constrained by experimental data. The uniform and consistent descriptions of multiple experimental observables demonstrate the essential role of modified coherence effects and the multistage modeling of jet evolution. Using the best choice of parameters from [], and with no further tuning, we present calculations for the medium modified jet fragmentation function, the groomed jet momentum fraction z g and angular separation r g distributions, as well as the nuclear modification factor of groomed jets. These calculations provide accurate descriptions of published data from experiments at the Large Hadron Collider. Furthermore, we provide predictions from the multistage model for future measurements at the BNL Relativistic Heavy Ion Collider. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Characterization of lift force and torque in prolate ellipsoid suspensions

The paper derives correlations for lift force and fluid torques acting on stationary prolate ellipsoid suspensions of aspect ratio (AR) 2.5, 5, and 10 subjected to uniform flow. Here, Particle Resolved Simulations (PRS) are conducted on a suspension of infinite extent in two directions for Reynolds number 10≤Re≤200 and solid fractions (φ) between 0.1≤φ≤0.3.The suspension-mean lift-to-drag ratio varies between 7% to 14% at Re=10 which increases to 14%~22% at Re=200. The torque-induced tip rotational acceleration can reach 38%~85% of drag-induced translational acceleration at Re=200. Single particle lift force and torque correlations of (Fröhlich et al., 2020) are modified and adapted to predict current angular-mean lift and torque data. The resulting lift correlation captures the PRS data within an average deviation below 7%. Torque exhibits a somewhat more complex dependency at AR=10 than the assumed sinθ ∙ cosθ variation but nevertheless the correlations predict angular-mean values with mean relative deviations of 16.6% at AR=10.

42 ENGINEERING↗
DOE OSTI2026

Competition between excitonic insulators and quantum Hall states in correlated electron–hole bilayers

Excitonic insulators represent a unique quantum phase of matter that enables the study of exotic quantum bosonic states. Strongly coupled electron–hole bilayers, which host stable dipolar exciton fluids with an exciton density that can be adjusted electrostatically, offer an ideal platform to investigate correlated excitonic insulators. On the basis of electron–hole bilayers made of MoSe2/hexagonal boron nitride/WSe2 heterostructures, here we study the behaviour of excitonic insulators in a perpendicular magnetic field. We report the observation of excitonic quantum oscillations in both Coulomb drag signals and electrical resistance at low to medium magnetic fields. Under a strong magnetic field, we identify multiple quantum phase transitions between the excitonic insulator phase and the bilayer quantum Hall insulator phase. These findings underscore the interplay between the electron–hole interactions and Landau-level quantization, and enable further exploration of quantum phenomena in composite bosonic insulators.

Qi, Ruishi↗
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 OSTI2020

High-Temperature Heat Transfer Fluid Circulator for Concentrating Solar Power Systems

Concentrating Solar Power (CSP) plants with Thermal Energy Storage (TES), require compressors capable of circulating supercritical carbon dioxide (sCO2) heat transfer fluids (HTFs) at temperatures from 500 to 700+°C and pressures from 8 to 25 MPa (80 to 250 Bar or 1160 to 3626 psia) in order to achieve thermal-to-electric efficiencies >50%. A review of papers relating to SCO2 compressors reveals that in almost every case the conditions being considered either allow for cooling of the SCO2 working fluid or are considering operating conditions close to the critical point at a temperature and pressure of approximately 31°C and 1070 psi). While the technology reported in these papers is relevant to most turbomachinery, the current CSP gas pathway circulator operating pressures are well above the CO2 critical point and do not allow for cooling of the main HTF flow. The purpose of this paper was to describe the key system issues affecting the design of a HTF circulator that could be scaled from a small size research and proof of concept circulator compressor to sizes supporting power plants producint upwards of 100 to 300 Mega Watts of electrical power. Due to the high temperatures and pressures the research scale HTF circulator was designed to be hermetically sealed to prevent leakage of the super critical carbon dioxide to the environment. In order to achieve hermetic sealing the drive motor and centrifugal impeller were fully integrated in a single shaft system configuration. From a preliminary assessment of compression power as a function of operating speed and impeller diameter it was determined that a smaller, higher speed impeller would result in a more efficient circulator. However, due to the high pressure and corresponding density of the sCO2 coupled with high operating speeds, large windage losses (i.e., viscous drag on the high speed rotating shaft) would occur which could result in a substantial reduction in efficiency of the system. Therefore tradoff studies were complerted to identify the operating speeds and system component sizes to achieve high efficiency. Based on the trade studies, the final design was a 55 mm diamter compressor impeller operating at speeds to 85,000 rpm. The unit will be built and tested in the future. Key Words: Super Critical Carbon Dioxide, Heat Transfer Fluid, Turbomachinery, Compressor, Circulator, Critical Point, Concentrating Solar Power, Solar Energy

Super Critical Carbon Dioxide, Heat Transfer Fluid↗
DOE OSTI2021

Settling of two-way momentum and energy coupled particles subject to Boussinesq and non-Boussinesq heating

This work establishes a procedure to accurately compute heat transfer between an Eulerian fluid and Lagrangian point-particles. Recent work has focused on accurately computing momentum transfer between fluid and particles. The coupling term for momentum involves the undisturbed fluid velocity at the particle location which is not directly accessible in the simulation. Analogously, in the context of thermal coupling, the undisturbed fluid temperature at the particle location is not directly accessible in simulations and must be estimated. In this paper, we develop a scheme to accurately estimate the undisturbed fluid temperature of a point-particle exchanging thermal energy with a surrounding fluid. Furthermore, the temperature disturbance is correlated with the enhanced temperature curvature in the vicinity of the particle and is formally valid in the low heating, low convection limit. We conduct extensive verification of the correction procedure for a settling particle subject to radiation. This setup allows the simultaneous testing of thermal and momentum corrections. By considering equations of drag and Nusselt number extended to finite Péclet and Boussinesq numbers, we establish a large range over which the correction procedure can be applied.

42 ENGINEERING↗
DOE OSTI2022

OC6 Phase Ia: CFD Simulations of the Free-Decay Motion of the DeepCwind Semisubmersible

Currently, the design of floating offshore wind systems is primarily based on mid-fidelity models with empirical drag forces. The tuning of the model coefficients requires data from either experiments or high-fidelity simulations. As part of the OC6 (Offshore Code Comparison Collaboration, Continued, with Correlation, and unCertainty (OC6) is a project under the International Energy Agency Wind Task 30 framework) project, the present investigation explores the latter option. A verification and validation study of computational fluid dynamics (CFD) models of the DeepCwind semisubmersible undergoing free-decay motion is performed. Several institutions provided CFD results for validation against the OC6 experimental campaign. The objective is to evaluate whether the CFD setups of the participants can provide valid estimates of the hydrodynamic damping coefficients needed by mid-fidelity models. The linear and quadratic damping coefficients and the equivalent damping ratio are chosen as metrics for validation. Large numerical uncertainties are estimated for the linear and quadratic damping coefficients; however, the equivalent damping ratios are more consistently predicted with lower uncertainty. Some difference is observed between the experimental and CFD surge-decay motion, which is caused by mechanical damping not considered in the simulations that likely originated from the mooring setup, including a Coulomb-friction-type force. Overall, the simulations and the experiment show reasonable agreement, thus demonstrating the feasibility of using CFD simulations to tune mid-fidelity models.

17 WIND ENERGY↗
DOE OSTI2023

Evaluation of Spray and Combustion Models for Simulating Dilute Combustion in a Direct-Injection Spark-Ignition Engine

Dilute combustion in spark-ignition engines has the potential to improve thermal efficiency by mitigating knock and by reducing throttling and wall heat losses. However, ignition and combustion processes can become unstable for dilute operation due to a lowered laminar flame speed, resulting in excessive cycle-to-cycle variability (CCV) of the combustion process. To compensate for the slower combustion in less reactive mixtures, a modified intake port geometry can be employed to generate a strong tumble flow in the cylinder and elevate turbulence levels around the spark plug, thereby promoting a faster transition to turbulent deflagration. Consequently, optimizing combustion chamber geometry and operating strategy is crucial to maximizing the benefits of using dilute combustion with enhanced in-cylinder turbulence across a wide range of operating conditions. Computational fluid dynamics (CFD) simulations can be utilized for virtual engine optimization tasks, but this would require the models to be truly predictive regarding the impact of changes to the engine design and operational parameters.In this study, multicycle large-eddy simulations (LES) are performed for a direct-injection spark-ignition engine to investigate the model performance in predicting engine combustion characteristics with respect to changes in the intake configuration. A tumble plate that blocks the lower part of the intake port inlet is used to vary the tumble. A set of CFD models that have been recently developed are employed, which takes into account the drag of nonspherical droplets, flash-boiling behavior of liquid sprays, spray-wall interaction, surrogate formulation of a research-grade E10 gasoline, and fast chemical kinetic solvers. Simulation results are compared to experimental engine data in terms of cylinder pressure, apparent heat release rate, mass fraction burned timing, and flame images. It is found that LES employing the state-of-the-art CFD models are capable of properly predicting the spray processes and reproducing the measured mean cylinder pressure for the case with the tumble plate. On the other hand, the LES over-predicts the combustion rate during the early combustion stage and under-estimates the CCV, and these discrepancies become larger when the tumble plate is removed.

computational fluid dynamics simulation↗
DOE OSTI2022

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

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

Energy & Fuels↗
DOE OSTI2024

Viscosity Measurements at High Pressures: A Critical Appraisal of Corrections to Stokes' Law

Abstract Fluids and melts in planetary interiors significantly influence geodynamic processes from volcanism to global‐scale differentiation. The roles of these geofluids depend on their viscosities ( η ). Constraining geofluid η at relevant pressures and temperatures relies on laboratory‐based measurements and is most widely done using Stokes' Law viscometry with falling spheres. Yet small sample chambers required by high‐pressure experiments introduce significant drag on the spheres. Several correction schemes are available for Stokes' Law but there is no consensus on the best scheme(s) for high‐pressure experiments. We completed high‐pressure experiments to test the effects of (a) the relative size of the sphere diameter to the chamber diameter and (b) the top and bottom of the chamber, that is, the ends, on the sphere velocities. We examined the influence of current correction schemes on the estimated viscosity using Monte Carlo simulations. We also compared previous viscometry work on various geofluids in different experimental setups/geometries. We find the common schemes for Stokes' Law produce statistically distinct values of η . When inertia of the sphere is negligible, the most appropriate scheme may be the Faxén correction for the chamber walls. Correction for drag due to the chamber ends depends on the precision in the sinking distance and may be ineffective with decreasing sphere size. Combining the wall and end corrections may overcorrect η . We also suggest the uncertainty in η is best captured by the correction rather than propagated errors from experimental parameters. We develop an overlying view of Stokes' Law viscometry at high pressures.

Geochemistry & Geophysics↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2022

Removal of Non-Specifically Bound Proteins Using Rayleigh Waves Generated on ST-Quartz Substrates

Label-free biosensors are plagued by the issue of non-specific protein binding which negatively affects sensing parameters such as sensitivity, selectivity, and limit-of-detection. In the current work, we explore the possibility of using the Rayleigh waves in ST-Quartz devices to efficiently remove non-specifically bound proteins via acoustic streaming. A coupled-field finite element (FE) fluid structure interaction (FSI) model of a surface acoustic wave (SAW) device based on ST-Quartz substrate in contact with a liquid loading was first used to predict trends in forces related to SAW-induced acoustic streaming. Based on model predictions, it is found that the computed SAW body force is sufficient to overcome adhesive forces between particles and a surface while lift and drag forces prevent reattachment for a range of SAW frequencies. We further performed experiments to validate the model predictions and observe that the excitation of Rayleigh SAWs removed non-specifically bound (NSB) antigens and antibodies from sensing and non-sensing regions, while rinsing and blocking agents were ineffective. An amplified RF signal applied to the device input disrupted the specific interactions between antigens and their capture antibody as well. ST-quartz allows propagation of Rayleigh and leaky SH-SAW waves in orthogonal directions. Thus, the results reported here could allow integration of three important biosensor functions on a single chip, i.e., removal of non-specific binding, mixing, and sensing in the liquid phase.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Application of a Distributed Element Roughness Model to Additively Manufactured Internal Cooling Channels

Abstract Design for cooling effectiveness in turbine blades relies on accurate models for dynamic losses and heat transfer of internal cooling passages. Metal additive manufacturing (AM) has expanded the design space for these configurations, but can give rise to large-scale roughness features. The range of roughness length scales in these systems makes morphology resolved computational fluid dynamics (CFD) impractical. However, volumetric roughness models can be leveraged, as they have computational costs orders of magnitude lower. In this work, a discrete element roughness model (DERM), based on the double-averaged Navier–Stokes equations, is presented and applied to additively manufactured rough channels, representative of gas turbine blade cooling passages. Unique to this formulation of DERM is a generalized sheltering-based treatment of drag, a two-layer model for spatially averaged Reynolds stresses, and explicit treatment of dispersion. Six different AM rough surface channel configurations are studied, with roughness trough to peak sizes ranging from 15% to 60% nominal channel passage half-width, and the roughness Reynolds number ranges from Rek = 60 to 300. DERM predictions for spatially and temporally averaged mean flow quantities are compared to previously reported direct numerical simulation results. Good agreement in the mean velocity profiles, stress balances, and drag partitions are observed. While DERM models are typically calibrated to specific deterministic roughness morphologies at comparatively small roughness Reynolds numbers, the present more generalized DERM formulation has wider applicability. Here, it is demonstrated that the model can accommodate random roughness of large scale, typical of AM.

Engineering↗
DOE OSTI2020

CTF Theory Manual (V.4.2)

Coolant-Boiling in Rod Arrays—Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e. fluid film, fluid drops, and vapor) modeling approach. Both sub-channel and 3D Cartesian forms of 9 conservation equations are available for LWR modeling. The code was originally developed by Pacific Northwest Laboratory in 1980 and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently re-branded as CTF. As part of the improvement process, it was necessary to generate sufficient documentation for the open-source code which had lacked such material upon being adopted by RDFMG. This document serves mainly as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code as well as the numerical solution process utilized. Coding of the models is also discussed, all with consideration for updates that have been made when transitioning from COBRA-TF to CTF. Further documentation outside of this manual is also available at RDFMG which focus on code input deck generation and source code global variable and module listings

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗