Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LAMINAR FLOW”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
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 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 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 OSTI2025

Comparison of Artificial Compressibility Methods for Coupled Laminar Fluid Flow and Heat Transfer

Artificial compressibility (AC) methods for incompressible flows perturb the mass continuity equation by the time derivative of the pressure field to achieve direct coupling between pressure and velocity without the need for a projection method. Several variations of the original AC method have been proposed (e.g., the entropically damped artificial compressibility [EDAC] method) that can, in some scenarios, improve on the original model. In this paper, we methodically compare the results of the original AC model and the EDAC model for both steady-state and transient problems with coupling between the momentum and energy equations. We demonstrate that the EDAC model and the AC model yield similar numerical results, and the EDAC model is 10%-20% faster than the AC model for a given problem.

Oz, Furkan [ORNL] (ORCID:0000000265831724)↗
DOE OSTI2024

Inertia-induced mixing and reaction maximization in laminar porous media flows

Solute transport and biogeochemical reactions in porous and fractured media flows are controlled by mixing, as are subsurface engineering operations such as contaminant remediation, geothermal energy production, and carbon sequestration. Porous media flows are generally regarded as slow, so the effects of fluid inertia on mixing and reaction are typically ignored. Here, we demonstrate through microfluidic experiments and numerical simulations of mixing-induced reaction that inertial recirculating flows readily emerge in laminar porous media flows and dramatically alter mixing and reaction dynamics. An optimal Reynolds number that maximizes the reaction rate is observed for individual pore throats of different sizes. This reaction maximization is attributed to the effects of recirculation flows on reactant availability, mixing, and reaction completion, which depend on the topology of recirculation relative to the boundary of the reactants or mixing interface. Recirculation enhances mixing and reactant availability, but a further increase in flow velocity reduces the residence time in recirculation, leading to a decrease in reaction rate. The reaction maximization is also confirmed in a flow channel with grain inclusions and randomized porous media. Interestingly, the domain-wide reaction rate shows a dramatic increase with increasing Re in the randomized porous media case. This is because fluid inertia induces complex three-dimensional flows in randomized porous media, which significantly increases transverse spreading and mixing. This study shows how inertial flows control reaction dynamics at the pore scale and beyond, thus having major implications for a wide range of environmental systems.

Chen, Michael A. (ORCID:0000000294614710)↗
DOE OSTI2025

Optimization of laser-induced breakdown spectroscopy for liquid sodium analysis

Next-generation nuclear power reactors based on liquid sodium coolant are in development to improve efficiency and safety. Fast breeder reactors could potentially achieve considerably higher fuel burn-up while producing less nuclear waste compared with current commercial technology. Monitoring of impurities in these types of reactors is important in order to prevent blockage of the cooling system or hazardous reactions due to the high reactivity of sodium with water. In this work, laser-induced breakdown spectroscopy is performed on a liquid sodium sample using a custom-built immersion probe system. The probe provides direct line of sight to the sample vessel in argon atmosphere for laser delivery and signal collection. Signatures of sodium, and impurities such as potassium and oxygen have been obtained. The probe height and immersion have a strong influence on the signal level and stability. Argon flow rates of 2 to 5 liters per minute produce no substantial change in signal, suggesting a steady laminar flowing regime. A high signal-to-noise ratio for impurity elements confirms the suitability of this technique for reliable impurity detection and monitoring in the reactor environment.

Gardette, V. [University of Michigan, Ann Arbor, M↗
DOE OSTI2025

Uniting Surface Properties With Hydrodynamic Roughness in Shallow Overland Flow Models

Describing flow resistance from the properties of an underlying surface is a challenge in surface hydrology. Runoff models must specify a resistance formulation or “roughness scheme”—describing the functional relationship between flow resistance and flow depth/velocity—and its parameters. Uncertainty in runoff predictions derives from both the selected roughness scheme (e.g., Darcy Weisbach, Manning's, or laminar flow equations), and its parameterization with a roughness coefficient (e.g., Manning's n). Both choices are informed by model calibration to data, usually discharge, and, if available, velocity. In this study, a Saint Venant Equation-based runoff model is calibrated to discharge and velocity data from 112 rainfall simulator experiments. The results are used to identify the optimal roughness scheme among four widely-used options for each experiment, and to explore whether surface properties can be used to select the optimal roughness scheme and its coefficient. Among the tested roughness schemes, a transitional flow equation provided the best fit to the plurality of experiments. The most suitable roughness scheme for a given experiment was not related to measured surface properties. Regression models predicted the calibrated roughness coefficients with adjusted values between 0.48 and 0.54, depending on the roughness scheme used. Litter cover was the best predictor of the roughness coefficient, followed by soil cover and average canopy gap size. The results suggest that selection of an optimal roughness scheme based on surface properties alone remains difficult, but that once a scheme is selected, roughness coefficients can be estimated from surface properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Structural uniformity and compositional homogeneity of solid-phase alloyed rod

Solid-phase processes have emerged as an alternative to fusion-based alloying to avoid coarse microstructures, undesirable phase formation, and high energy consumption. However, achieving uniform distribution of alloying elements during friction-based processing remains challenging due to highly heterogeneous thermomechanical conditions. This work evaluates the structural uniformity and compositional homogeneity of Al–Cu–Zn alloyed rods produced by friction extrusion (FE) and establishes the role of the rotational speed to feed rate ratio (N/V) on alloying effectiveness. A systematic matrix of FE experiments was conducted at constant extrusion ratio with N/V values ranging from 3.7 to 300. Compositional uniformity was assessed along the rod length (ICP-OES), in three dimensions (X-ray computed tomography), and at the microscale (SEM–EDS), supported by a gray-level co-occurrence matrix (GLCM)–based homogeneity metric. Smoothed particle hydrodynamics (SPH) simulations were used to reveal material flow and thermomechanical fields. Results show that N/V = 100 produces a high-shear mixing zone that eliminates the unmixed core and enables near-full dissolution and dispersion of Cu and Zn. At lower N/V, a laminar flow region persists at the rod center, causing segregation and large composition gradients. The combined experimental–computational analysis provides mechanistic insight into the transition from fragmented particle dispersion to thermomechanically assisted metallurgical mixing. This study establishes processing–structure relationships for solid-phase alloying and provides guidance for achieving homogenized compositions comparable to wrought alloys via rapid, scalable FE processing.

Aluminum↗
DOE OSTI2025

Flow Reactor Study of the Soot Precursors of Novel Cycloalkanes as Synthetic Jet Fuel Compounds: Octahydroindene, p -Menthane, and 1,4-Dimethylcyclooctane

Sustainable aviation fuels (SAFs) or Synthetic aviation turbine fuels (SATFs) derived from nonpetroleum sources are essential for energy security and a strong rural and agricultural economy. Airplanes operating on SAF can have lower particle emissions compared to those of conventional jet fuel, reducing air quality impacts near airports. Processing biobased isoprene or wood and agricultural waste can produce cycloalkane-rich fuels with properties meeting ASTM International’s SATF requirements. The unique structures of these cycloalkanes yield lower soot emissions because of their lack of aromatic rings. We measured the soot formation tendency as yield sooting index (YSI) and used laminar flow reactor experiments to evaluate soot precursors formed for isoprene-derived compounds p-menthane and 1,4-dimethylcyclooctane (DMCO), and octahydroindene (OHI)─ produced from woody biomass via catalytic fast pyrolysis. The combustion chemistry of the OHI and DMCO has not been previously studied. Experiments were conducted at 10 bar from 800 to 1200 K, equivalence ratios of 1.0 and 3.0, and residence times of 1.0 and 0.6 s, respectively. Experimentally detected species were used to elucidate the mechanisms of soot precursor formation. OHI exhibited the highest YSI (94.5) and formed a high concentration of benzene primarily by direct dehydrogenation of the six-membered ring. p-Menthane (YSI 92.0) and DMCO (YSI 85.0) oxidation products included fewer aromatic components but higher benzene precursors, including 1,3-butadiene, propyne, and allene. This suggests that the ring-opening pathway is dominant over the dehydrogenation pathway in the benzene formation for these compounds. This experimental speciation provides insight into the influence of the cycloalkane structure on the sooting tendencies of potential SAF blend components, thereby aiding in fuel design processes.

09 BIOMASS FUELS↗
DOE OSTI2026

Design and demonstration of a direct air capture system with moisture-driven CO 2 delivery into aqueous medium

Two direct air capture (DAC) systems were designed and demonstrated to passively capture CO 2 from ambient air and use moisture to release the CO2 into an alkaline medium. A bench-scale system delivering ∼1 g CO 2 d –1 was demonstrated in a laminar flow hood, and a small pilot-scale system that could deliver ∼100 g CO 2 d –1 was operated outdoors in a 4.2 m 2 raceway pond. Novel elongated mesh-tube packets containing anion-exchange resin (AER) beads were found to reduce drying and CO 2 loading time 4.3-fold compared to larger mesh bags. Technoeconomic analysis (TEA) estimates the cost of capturing CO 2 into an alkaline solution, suitable for cultivating photosynthetic microorganisms, to be $\$229$ per tonne for a practical scenario based on current results and $\$72$ per tonne for an aspirational scenario considering improvements to sorbent capacity, hydrophobicity, and sorbent lifetime. TEA further estimates an additional $\$110$ per tonne to extract CO 2 from solution, purify it, and compress it to 15 MPa, suitable for sequestration. Furthermore, moisture-driven processes have the potential to use up to 87% less energy than thermal and/or vacuum swing DAC by using energy from water evaporation.

09 BIOMASS FUELS↗
DOE OSTI2024

Limiting current in a collisional crossed-field gap

Crossed-field devices are often used in pulsed power and high-power microwave applications. Previous studies derived closed-form solutions for the limiting current of a vacuum crossed-field system, corresponding to the maximum permissible current for laminar flow, below and above the Hull cutoff B H for magnetic insulation. We extend these studies by introducing collision frequency into the electron force law as a friction term to derive the limiting current in a collisional crossed-field gap. The resulting solution recovers the vacuum crossed-field case in the limit of no collisions and the collisional space-charge limited current with general initial velocity for magnetic field B → 0. In the limit of infinite collisions, we obtain a crossed-field equivalent to the Mott–Gurney law for the maximum current permissible in a collisional, nonmagnetic diode. When the collision frequency ν is less than the electron cyclotron frequency Ω, increasing initial velocity makes the critical current nonmonotonic with increasing ν with the critical current higher at B = B H for ν = Ω. As for a misaligned crossed-field gap where a component of the magnetic field was introduced parallel to the electric field across the gap, magnetic insulation is eliminated and the discontinuity at B = B H for limiting current observed in a vacuum crossed-field gap vanishes. As B → ∞, the limiting current approaches a constant that depends on the initial velocity and the collision frequency.

Breen, Lorin I.↗
DOE OSTI2025

Mesoscale fractal whey protein particles derived from microscale linear-shaped protein assemblies (Part 1): Manufacturing method and particle characteristics

Whey protein isolates (WPI) are widely used in processed foods for their versatile functional properties. Modifying the structural properties of proteins by assembling them into mesoscale or microscale particles may improve their functionality and broaden their applications. This study aims to manufacture and characterize mesoscale whey protein particles (WPP) derived from WPI. Two types of WPP, WPP1 (0.05 mL/min) and WPP2 (0.25 mL/min), were prepared through a multistep approach involving liquid antisolvent (LAS) precipitation, heat treatment, and microfluidization. Liquid antisolvent precipitation was performed by injecting a 20% (wt/vol) WPI dispersion (pH 7) into an ethanol-glycerol mixture (75:25, vol/vol) under laminar flow, followed by heat treatment at 80°C for 20 min as a particle hardening step. This process produced stable fiber- and ribbon-shaped whey protein assemblies (WPA), which served as precursors to WPP. Subsequent microfluidization (150 MPa, 6 passages) reduced the size of WPA, yielding mesoscale WPP with irregular morphologies and a more uniform size distribution, as revealed by microscopy and dynamic light scattering. ζ-Potential and fluorescence labeling indicated higher surface charge and surface hydrophobicity of WPP compared with untreated WPI. The WPP showed internal mass fractal and surface fractal structures at larger length scales, analyzed using small-angle X-ray scattering. Fourier transform infrared spectroscopy demonstrated an increased fraction of intermolecular β-sheets in WPP, suggesting that hydrogen bonding contributed to their formation. Gel electrophoresis confirmed that disulfide bonds served as the primary cross-links stabilizing the WPP structure. Furthermore, turbidity measurements showed that WPP exhibited superior colloidal phase stability compared with untreated WPI and maintained high colloidal stability under both acidic and neutral pH conditions.

Antisolvent precipitation↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗
DOE OSTI2025

Openpronghorn

OpenPronghorn is a simulation tool specifically tailored for modeling thermal-hydraulic phenomena in advanced nuclear reactors. It is built on the Multiphysics Object-Oriented Simulation Environment (MOOSE), an open-source platform that facilitates the development of high-performance scientific computing applications. OpenPronghorn solves the Navier-Stokes equations, which describe the conservation of mass, momentum, and energy in fluid flows, using the finite volume numerical method. The code supports a wide range of fluid flow conditions that are applicable to nuclear reactors, including incompressible and weakly compressible flows, as well as single-phase and multiphase flows. It is capable of modeling diverse flow regimes, including laminar and turbulent flows, using various turbulence models such as the standard k-epsilon models, the v2f model, and the mixing length model. For multiphase flows, OpenPronghorn employs a mixture a Eulerian modeling approach with mixture, drift-flux, and full Eulerian models, and includes open-sourced interfacial transfer correlations for drag, exchange, and heat transfer coming from the scientific literature. OpenPronghorn's modular design allows it to handle multiscale simulations, ranging from detailed Reynolds-Averaged Navier Stokes (RANS) simulations to coarse-mesh and lumped parameter models. This flexibility enables users to perform high-fidelity simulations of specific reactor components as well as system-level analyses of entire reactor circuits. The code can be coupled with other MOOSE-based tools using the MultiApp system, allowing for the transfer of coupling quantities such as mass flow rates, heat fluxes, and boundary conditions between different simulation scales. One of the main features of OpenPronghorn is the it includes built-in validation cases from the open-source scientific literature and supports the implementation of user-defined models and correlations through MOOSE's FunctorMaterial system. OpenPronghorn is designed to be computationally efficient, leveraging the SIMPLE projection method for large-scale problems, and can be run on high-performance computing systems to handle the extensive computational demands of detailed reactor simulations. Overall, OpenPronghorn is a versatile and robust tool that provides critical insights into the thermal-hydraulic behavior of advanced nuclear reactors, supporting the design, safety, and optimization of next-generation nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2024

Overview of On-Line Optical Measurements at High Pressure for Flue Gases, Particulates and Acid-dew Point of Pressurized Oxy-Combustion

Optical flow cells are critical measurement interfaces, yet sampling under harsh conditions — high pressure, high temperature, particles, moisture, or corrosive gases — makes it difficult to maintain optical quality without perturbing the measurement. To address this challenge, a new flow cell was developed using a laminar coaxial flow field that separates the purge and sample flows. A dedicated test system was built to evaluate particle size distribution (PSD) measurements using a Malvern Panalytical Insitec analyzer. Results demonstrated that the sample flow alone defines the measurement zone, while the purge flow effectively shields the optical windows from deposition, eliminating sampling bias. The flow cell enables reliable PSD measurement under high pressure and temperature in moist, corrosive environments. As a key demonstration, the instrument was successfully deployed for on-line PSD measurement of flue gas from a 100 kWth pressurized oxy-coal combustor at 15 bara.

Cheng, Mao↗
DOE OSTI2025

A Novel Optical Instrument for On-Line Measurement of Particle Size Distribution—Application to Clean Coal Technologies

A flow cell is a critical measurement interface for many optical instruments. However, the flows are often sampled under harsh conditions, such as under high pressure and/or high temperature, in the presence of particles, moisture, vapors with high dew points or corrosive gases. Therefore, obtaining a high-optical-quality flow cell that does not perturb the measurement is a significant challenge. To address this challenge, we proposed a new flow cell that employs a unique laminar coaxial flow field (for the purge and sample flows). A test system was built to conduct particle size distribution (PSD) measurements with no sampling bias using a state-of-the-art analyzer (Malvern Panalytical Insitec). The results revealed that the measurement zone is well defined solely by the sample flow, and the optical windows are well protected by the purge flow, with minimal risk of any depositions from the sample flow. Using this flow cell, the Insitec can successfully measure PSD under high pressure and temperature under moist, corrosive conditions without generating any sampling bias. Importantly, we successfully applied this flow cell for on-line PSD measurement for the flue gas of a 100 kWth pressurized oxy-coal combustor operating at 15 bara.

Cheng, Mao (ORCID:0000000213273900)↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗