Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLOW, LAMINAR”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Glovebox CFD Analysis: Low-Flow, Laminar CFD Analysis of Glovebox Ventilation and Humidity Transport in Ansys Fluent 2025 R1

There are gloveboxes that frequently exceed the 30% relative humidity level required for operations. These exceedances result in alternative operations needing to be utilized. This report documents a computational fluid dynamics (CFD) analysis performed to characterize airflow and humidity transport within a representative glovebox under nominal, low-flow conditions. The analysis focuses on establishing the flow regime, identifying transport-limited regions, and evaluating the implications of laminar flow on humidity removal. A single-phase, laminar flow model with a species transport model was used to represent air and water vapor mixing. Results demonstrate that, at the nominal flow rates, airflow is laminar with weak mixing, which will result in long residence times, tens of minutes long or more, in low-velocity regions. These conclusions suggests that humidity reduction is transport-limited. This analysis provides qualitative insight into airflow structure and humidity transport trends. In addition to characterization of glovebox airflow and humidity transport, results are intended to inform experimental validation and future transient modeling efforts.

42 ENGINEERING↗
DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radiative heat transfer in FLiBe molten salt participating medium in a vertical heated tube under forced and mixed convection laminar flows

The contribution of radiative heat transfer (RHT) to convective heat transfer in a heated vertical pipe with laminar flow and constant wall temperature, for FLiBe molten salt, is investigated computationally. The combined effects of conduction in the fluid, forced convection, buoyancy, temperature-dependent physical properties, and thermal radiation are investigated. The P1 approximation is employed in the discretization of the Radiative Transfer Equation (RTE). The COMSOL Multiphysics software is used to generate the numerical solutions. The theoretical analysis developed here demonstrates that only for intermediate values of the optical thickness, i.e. τ D ~ O 1 , the participating media effects are expected to be important. This analysis is confirmed computationally and a value of τ D ≈ 4 is shown to lead to the highest increase in overall heat transfer behavior. Under forced convection, the Nusselt ratio N u total / N u no - rad reaches a peak value of 1.76 at τ D = 4 and z/D = 200 , and is less than 1.1 for τ D < 0 . 1 and τ D > 60 . Under aiding-flow mixed convection, N u total / N u no - rad reaches a peak of 1.34 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 4 and τ D > 36 . Under opposing-flow mixed convection, Nu total / Nu no - rad reaches a peak of 1.81 at τ D = 4 . 2 and is less than 1.1 for τ D < 0 . 2 and τ D > 50 . RHT effects on the overall heat transfer are most pronounced in opposing mixed convection, where N u total / N u no - rad is 2.26 observed at z/D = 50 . The sensitivity to wall emissivity is evaluated; in forced convection, the peak N u total / N u no - rad is 1.66 at ε = 0 (reflective wall) and 2.04 at ε = 1 (absorptive wall). Finally, the sensitivity to pipe diameter is also discussed. Overall, for a vertical heated tube, radiative heat transfer effects lead to enhancement of heat transfer by as high as a factor of two, and they depend on the optical thickness of the flow, mixed convection environment ( Gr / R e 2 and direction of flow relative to the gravitational force), surface emissivity, and entrance effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Laminar Flow Heat Transfer In Helical Oval-Twisted Tube for Heat Exchanger Applications

The heat transfer performance of a novel tube configuration that combines the swirling velocity induced by oval-twisting and secondary flow generated by a helical geometrical flow path is presented. The Nusselt number (Nu) and friction factor (f) are compared for the laminar flow regime (Re = 250- 2000) under isothermal wall conditions. Under the same flow and boundary conditions, the oval-twisted helical tube increased the Nu and slightly increased the f over the circular helical tube. The best performance with the highest Nuand lowest f occurs at the coil curvature ratio (d h /D) of 0.17. The quantified enhancement performance factor (h) shows a 46 - 56% increase for the oval-twisted helical tube over the circular helical tube with d h /D = 0.17. Correlations for the Nu and f are developed and agreed with the numerical results: ± 6% and ± 5% for the Nu and f, respectively.

42 ENGINEERING↗
DOE OSTI2023

Applying Physics-Informed Neural Networks to Solve Navier–Stokes Equations for Laminar Flow around a Particle

In recent years, Physics-Informed Neural Networks (PINNs) have drawn great interest among researchers as a tool to solve computational physics problems. Unlike conventional neural networks, which are black-box models that “blindly” establish a correlation between input and output variables using a large quantity of labeled data, PINNs directly embed physical laws (primarily partial differential equations) within the loss function of neural networks. By minimizing the loss function, this approach allows the output variables to automatically satisfy physical equations without the need for labeled data. The Navier–Stokes equation is one of the most classic governing equations in thermal fluid engineering. This study constructs a PINN to solve the Navier–Stokes equations for a 2D incompressible laminar flow problem. Flows passing around a 2D circular particle are chosen as the benchmark case, and an elliptical particle is also examined to enrich the research. The velocity and pressure fields are predicted by the PINNs, and the results are compared with those derived from Computational Fluid Dynamics (CFD). Additionally, the particle drag force coefficient is calculated to quantify the discrepancy in the results of the PINNs as compared to CFD outcomes. The drag coefficient maintained an error within 10% across all test scenarios.

Hu, Beichao (ORCID:0009000163215151)↗
DOE OSTI2021

Taylor dispersion in osmotically driven laminar flows in phloem

Sucrose is among the main products of photosynthesis that are deemed necessary for plant growth and survival. It is produced in the mesophyll cells of leaves and translocated to different parts of the plant through the phloem. Progress in understanding this transport process remains fraught with experimental difficulties, thereby prompting interest in theoretical approaches and laboratory studies. The Münch pressure and mass flow model is one of the accepted hypotheses describing the physics of sucrose transport in the phloem. It is based on osmosis creating an energy potential difference between the source and the sink. The flow responding to this energy potential is assumed laminar and described by the Hagen–Poiseuille equation. This study revisits such osmotically driven flows in tubes with membrane walls by including the effects of Taylor dispersion on mass transport. This effect has been overlooked in phloem flow studies. Taylor dispersion can increase the effective transport of solutes by increasing the apparent diffusion coefficient. It is shown that, in addition to the conventional diffusive correction derived for impermeable tube walls, a new adjustment to the mean advective terms arises because of osmotic effects. Because the molecular Schmidt number is very large for sucrose in water, the sucrose front speed and travel times have a direct dependence on the Péclet number for different ranges of the Münch number. Fially, this study establishes upper limits on expected Taylor dispersion enhancement of sucrose transport.

42 ENGINEERING↗
DOE OSTI2022

Critical Shear Stress for Erosion Under Laminar Jet Flow

Particle erosion is an important process in both natural and manmade environments. For example, in natural environments particle or sediment erosion rates determine the course of rivers, depth of streams, and fate of deltas. The interplay between erosion rates and flow patterns can have dramatic impact on shorelines. Where hydrodynamic jet flows or rip currents change direction, sediments may accumulate or disperse altering the shape of coastlines and the ecosystems and economies that rely on them. In man-made systems, jet flows may be used to mobilize particles from surfaces or within industrial scale mixing tanks. For example, in the nuclear industry, impinging jet mixers may be used to mobilize and suspend sludges at the bottom of waste tanks. In each of these examples, the Shields diagram is useful to determine the conditions under which erosion occurs. This diagram plots the scaled critical shear stress for erosion as a function of the particle size as embedded in particle Reynolds numbers, Archimedes numbers, or their equivalents. Although the data, that carries a large spread, is largely in agreement with mathematical models, for the smallest particles, model predictions diverge from each other and from experimental data. Indeed, there is substantial uncertainty about the Shields diagrams at relatively small particle sizes and flow rates. Here we evaluate critical shear stresses for erosion under laminar jet flow. To date critical shear stresses for erosion, as embodied in Shields diagrams, have largely been developed for turbulent flow both experimental and mathematically. However, curves on the Shields diagrams diverge among various models at lower particle Reynolds numbers. Indeed, how critical shear stresses for erosion develop under fully laminar jet flow conditions (not simply laminar boundary layers under turbulent flow) remains unclear. Here we address this gap, which is particularly important for mobilization of the smallest particles. We find the Shields parameter to be inversely proportional to the particle Reynolds numbers under laminar jet flow conditions similar to turbulent conditions.

critical shear stresses, erosion rates, effective ↗
DOE OSTI2023

Gravitational vortex disruption in mesofluidic separators

Typically, bump array devices separate submicrometer particles where gravitational effects remain small. However, as particle size increases from a few micrometers to hundreds of micrometers, gravitational effects become important and settling may induce vortex formation. Here, this article shows under laminar flow that these vortices suffice to disrupt flow fields, frustrating particle separation. The limiting diameter, above which vortices do form, is determined, and the dependence of the limiting diameter on viscosity, particle properties, and separator orientation is quantified.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Characterization of the Finite Element Computational Fluid Dynamics Capabilities in the Multiphysics Object Oriented Simulation Environment

We report the multiphysics object-oriented simulation environment (moose) is a code package that couples a variety of physics modules, allowing for highly accessible multiphysics simulations. The physics modules include a finite element Navier–Stokes (N–S) module that is designed to solve laminar fluid dynamics problems. The usage of this module in multiple recent studies coupled with the growing interest in moose for usage in nonlight water reactor safety studies by the Nuclear Regulatory Commission (NRC) prompted the authors to investigate the computational fluid dynamics capabilities of moose. A two-dimensional laminar flow past a circular cylinder scenario is simulated in the moose framework to investigate the effectiveness of the N–S module. Simulations assumed an unsteady laminar flow with a Reynolds number of 200. To verify the results from moose, similar simulations were conducted using the well-utilized simulation of turbulent flow in arbitrary regions—computational continuum mechanics C++ (star-ccm + ) finite volume code. Results from both codes are also compared to some results from literature. Velocity and pressure profiles of both transient simulations were compared. The numerical and input errors in moose are also visualized with contour plots to qualitatively understand the evolution of the errors across time and space. The comparisons between moose and star-ccm + showed nearly perfect agreement between the codes for velocity and pressure, especially after the development of the vortex street in later time-steps. The force coefficients showed excellent agreement after the development of the vortex street, but demonstrated notable discrepancies prior to the vortex street development, which is likely due to how each code simulated the approach to the vortex street in earlier time-steps.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Effective Permeability of a Nuclear Fuel Assembly

This report aids in the development of models to perform characterization studies of aerosol dispersal and deposition within a spent fuel cask system. Due to the complex geometry in a spent-fuel canister, direct simulation of buoyancy-driven flow through the fuel assemblies to model aerosol deposition within the fuel canister is computationally expensive. Identification of an effective permeability as given in this work for a nuclear fuel assembly greatly simplifies the requirements for thermal hydraulic computations. The results of computations performed using OpenFOAM® to solve the Navier-Stokes Equations for laminar flow are used to determine an effective permeability by applying Darcy's Law. The computations are validated against an analytical solution for the special case of an infinite array of pins for which the numerical and analytical solutions have excellent agreement. The effective permeability of a 1717 PWR nuclear fuel assembly in a basket without spacer grids is numerically determined to be 1.85010 -6 m 2 for the range of fluid viscosities and pressure drops expected in a spent fuel storage canister. However, the flow is not uniform on the scale of multiple pins. Instead, significantly higher velocities are attained in the space between the assembly and the basket walls compared to the flow between the fuel pins within the assembly. Comparison with an analytical solution for fully developed flow through an infinite array of pins shows that the larger spacing near the basket walls results in about a 20% larger permeability compared to the analytical solution which does not include the enhanced flow in the space between the assembly and basket wall, or entrance and exit effects. A preliminary assessment of turbulence effects shows that with a k-epsilon model, significantly higher flow velocities are attained between the fuel pins within the assembly compared to the flow velocity in the space between the assembly and the basket walls. This is the opposite of what is determined for laminar flow.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Example of exponentially enhanced magnetic reconnection driven by a spatially bounded and laminar ideal flow

In laboratory and natural plasmas of practical interest, the spatial scale Δ d at which magnetic field lines lose distinguishability differs enormously from the scale a of magnetic reconnection across the field lines. In the solar corona, plasma resistivity gives a/Δ d ~ 10 12 , which is the magnetic Reynold number R m . The traditional resolution of the paradox of disparate scales is for the current density j associated with the reconnecting field Brec to be concentrated by a factor of Rm by the ideal evolution, so j ~ B rec /µ 0 Δ d . A second resolution is for the ideal evolution to increase the ratio of the maximum to minimum separation between pairs of arbitrarily chosen magnetic field lines, Δ max /Δ min , when calculated at various points in time. Reconnection becomes inevitable where Δ max /Δ min ~ R m . A simple model of the solar corona will be used for a numerical illustration that the natural rate of increase in time is linear for the current density but exponential for Δ max /Δ min . Reconnection occurs on a time scale and with a current density enhanced by only ln(a/Δ d ) from the ideal evolution time and from the current density B rec /µ 0 a. In both resolutions, once a sufficiently wide region, Δ r , has undergone reconnection, the magnetic field loses static force balance and evolves on an Alfv´enic time scale. The Alfv´enic evolution is intrinsically ideal but expands the region in which Δ max /Δ min is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Fluid inertia controls mineral precipitation and clogging in pore to network-scale flows

Mineral precipitation caused by fluid mixing presents complex control and predictability challenges in a variety of natural and engineering processes, including carbon mineralization, geothermal energy, and microfluidics. Precipitation dynamics, particularly under the influence of fluid flow, remain poorly understood. Combining microfluidic experiments and three-dimensional reactive transport simulations, we demonstrate that fluid inertia controls mineral precipitation and clogging at flow intersections, even in laminar flows. We observe distinct precipitation regimes as a function of Reynolds number (Re). At low Reynolds numbers (Re < 10), precipitates form a thin, dense layer along the mixing interface, which shuts precipitation off, while at high Reynolds numbers (Re > 50), strong three-dimensional flows significantly enhance precipitation over the entire intersection, resulting in rapid clogging. When injection rates from two inlets are uneven, flow symmetry-breaking leads to unexpected flow bifurcation phenomena, which result in enhanced concurrent precipitation in both downstream channels. Finally, we extend our findings to rough channel networks and demonstrate that the identified inertial effects on precipitation at the intersection scale are also present and even more dramatic at the network scale. This study sheds light on the fundamental mechanisms underlying mixing-induced mineral precipitation and provides a framework for designing and optimizing processes involving mineral precipitation.

Science & Technology - Other Topics↗
DOE OSTI2024

Gasification of Coal/Plastic Mixtures: Fundamental Studies in a Laminar Entrained Flow Reactor

There is interest in leveraging entrained-flow coal gasification for co-processing plastic waste. However interactions between plastic and coal, particularly possible competition between coal volatiles and plastic vapors for oxygen, are not well understood. This study explores co-gasification of pulverized coal and plastic granules in a lab-scale entrained-flow reactor. Conversion and CO/CO2 ratio are given special attention.

Kareem, Afeez↗