Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cross 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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Calculation of Friction Factors in Heat Pipes Using CFD in Support of the Sockeye Code

A heat pipe is two-phase heat transfer device which relies on surface tension and capillary pressure to provide a very efficient heat transfer mechanism. Currently, there is interest in designing micro nuclear power reactors using high temperature heat pipes to provide totally passive cooling. To support this, Sockeye is being developed as a heat pipe simulation capability. To support the development of Sockeye, individual phenomena must be tested. This work focuses on applying a CFD model to evaluate friction factors for use in Sockeye. A CFD model has been developed in Nek5000 which simulates each component fluid, i.e. sodium liquid and vapor, independently and captures the effects of the wick on friction factor. The model tests the effect of the wick permeability and the thickness of the wick occupied by each fluid component. The thickness of the wick occupied by each component can be directly related to the local vapor volume fraction from the Sockeye heat pipe model. Results from the CFD model indicate that friction factor in the liquid flow decreases with increasing wick permeability, while the friction factor in the vapor flow increases with increasing wick permeability. In the liquid, this is likely a direct result of the increased cross sectional area available for the flow. For the vapor, this is explained as the wick acting like a rough wall, while the increased area has only a negligible impact. Results for the tests of the thickness of the wick occupied by each fluid showed a thresholding behavior for both liquid and vapor. When the thickness of the wick occupied by each fluid was increased, it increased the effect of the wick, but only up to a certain point. Once a certain thickness was reached, the friction factors remained constant. When the fluid occupied only a very small fraction of the wick, friction factors returned to nominal values.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Large eddy simulation of particle dispersion analysis in a pressurized reactive fluid–particle system

Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new combustion technologies. Because of high pressure, the volumetric fraction of particles in a pressurized oxy-combustor could be higher than that in the conventional combustor, consequently, the particle dynamics in the pressurized vessel need to be evaluated accordingly. Specifically, the present computational work analyzes the pattern of the particle dispersion by extracting the particle mass flow rate and the particle size distribution in the pressurized combustor. A large eddy simulation (LES) with the environment pressure 15-bar and 100 kWth thermal input of fuel particle has been performed. The particle mass flow rates sampled in several cross-section planes strongly fluctuate, indicating the impacts of the flame as well as the gas-phase flow on the particle movement. Animations of the particle-gas flow show that the flow recirculation plays a major role in producing the pulsing phenomenon of the particle mass flow rate.

Li, Lei↗
DOE OSTI2024

Preliminary Assessment of Additively Manufactured Cooling Channel Performance for Helium-Cooled Blanket Concepts

Sufficient cooling of the blanket first wall remains a critical challenge for the design and deployment of fusion power plants. Helium has been targeted as a potential blanket coolant due to its inertness and low neutron interactivity, among other advantages. However, the low thermal mass of helium creates a need for heat transfer enhancements in coolant channels to provide adequate cooling to the blanket’s first wall. Toward this end, a series of ribbed flow channels of various rib cross sections and configurations has been produced via additive manufacturing (AM) to study the efficacy of AM for first wall heat transfer enhancement and the optimization of heat transfer geometries. Helium cooling performance is studied in AM test articles at 4 MPa operating pressure, Reynolds numbers up to 197 000, and outer surface heat fluxes up to 42 kW/m 2 in the recently commissioned helium flow loop experiment (HFLE). Preliminary results of this study are presented herein. Heat transfer performance of nominally smooth (i.e., featureless) AM channels is quantified via measured Nusselt numbers and friction factors and compared to off-the-shelf smooth pipe experiments. Results are compared to existing correlations and used to discuss the effects of the AM processes on thermal-hydraulic performance. It is seen that the inherent roughness of the AM channels leads to an increase in both heat transfer coefficient and pressure drop when compared to the conventional pipe. In conclusion, recommendations are made for future studies based on these findings and additional considerations for the deployment of AM blanket cooling components.

36 MATERIALS SCIENCE↗
DOE OSTI2020

On the competition between mixing rate and uniformity in a coaxial hydrodynamic focusing mixer

Fast microfluidic mixers for use with line-of-sight integrating detection schemes pose unique challenges. Such detectors typically cannot discriminate signal from slow moving (e.g. near internal walls) and fast-moving portions of the fluid stream. This convolves reaction rate dynamics with fluid flow residence time dynamics. Further, the small cross sections of typical three-dimensional hydrodynamic focusing devices lead to lower detection signals. The current study focuses on achieving both small time scales of mixing and homogenous residence times. This is achieved by injecting sample through a center capillary and hydrodynamically focusing using a sheath flow within a tapered second capillary. The current design also features a third, larger coaxial capillary. The mixed stream flows into the large cross-section of this third capillary to decelerate and expand the stream by up to 14-fold to improve line-of-sight signal strength of reaction products. Hydrodynamic focusing, mixing, and expansion are studied using analytical and numerical models and also studied experimentally using a fluorescein-iodide quenching reaction. The experimentally validated models are used to explore trade-offs between mixing rate and uniformity. For the first time, this work presents detailed analysis of the Lagrangian time history of species transport during mixing inside coaxial capillaries to measure mixing nonuniformity. The mixing region enables order 100 μs mixing times and residence time widths of the same order (140 μs).

42 ENGINEERING↗
DOE OSTI2023

Hydrodynamic irreversibility of non-Brownian suspensions in highly confined duct flow

The irreversible behaviour of a highly confined non-Brownian suspension of spherical particles at low Reynolds number in a Newtonian fluid is studied experimentally and numerically. In the experiment, the suspension is confined in a thin rectangular channel that prevents complete particle overlap in the narrow dimension and is subjected to an oscillatory pressure-driven flow. In the small cross-sectional dimension, particles rapidly separate to the walls, whereas in the large dimension, features reminiscent of shear-induced migration in bulk suspensions are recovered. Furthermore, as a consequence of the channel geometry and the development and application of a single-camera particle tracking method, three-dimensional particle trajectories are obtained that allow us to directly associate relative particle proximity with the observed migration. Companion simulations of a steadily flowing suspension highly confined between parallel plates are conducted using the force coupling method, which also show rapid migration to the walls as well as other salient features observed in the experiment. While we consider relatively low volume fractions compared to most prior work in the area, we nevertheless observe significant and rapid migration, which we attribute to the high degree of confinement.

42 ENGINEERING↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2021

Model, data, and code for paper "Modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x"

The data package includes data, model, and code that support the analyses and conclusions in the paper titled “modeling of streamflow in a 30-kilometer-long reach spanning 5 years using OpenFOAM 5.x”. The primary goal of this paper is to demonstrate that key streamflow properties such as water depth, flow velocity, and dynamic pressure in a natural river at 30-kilometer scale over 5 years can be reliably and efficiently modeled using the computational framework presented in this paper. To support the paper, various data types from remote sensing, field observations, and computational models are used. Specific details are described as follows. Firstly, the river bathymetry data was obtained from a Light Detection and Ranging (LiDAR) survey. This data is then converted to a triangulated surface format, STL, for mesh generation in OpenFOAM. The STL data can be found in Model_Setups/BaseCase_2013To2015/constant/triSurface. The OpenFOAM mesh generated using this STL file can be found in constant/polyMesh. Other model setups, boundary and initial conditions can be found in /system and /0.org under folder BaseCase_2013To2015. A similar data structure can also be found in BaseCase_2018To2019 for the simulations during 2018 and 2019. Secondly, the OpenFOAM simulations need the upstream discharge and water depth information at the upstream boundary to drive the model. These data are generated from a one-dimensional hydraulic model and the data can be found under the folder Model_Setups /1D model Mass1 data. The mass1_65.csv and mass1_191.csv files include the results of the 1D model at the model inlet and outlet, respectively. The Matlab source code Mass1ToOFBC20182019.m is used to convert these data into OpenFOAM boundary condition setups.With the above OpenFOAM model, it can generate data for water surface elevation, flow velocity, and dynamic pressure. In this paper, the water surface elevation was measured at 7 locations during different periods between 2011 and 2019. The exact survey locations (see Fig1_SurveyLocations.txt) can be found in folder Fig_1. The variation of water stage over time at the 7 locations can be found in folder /Observation_WSE. The data type include .txt, .csv, .xlsx, and .mat. The .mat data can be loaded by Matlab.We also measured the flow velocities at 12 cross-sections along the river. At each cross-section, we recorded the x, y locations, depth, three velocity components u,v,w. These data are saved to a Matlab format which can be found under folder /Observation_Velocity and /Fig_1. The relative locations of velocity survey locations to the river bathymetry can be found in Figure 1c.The water stage data at the 7 locations from OpenFOAM, 1D, and 2D hydraulic models are also provided to evaluate the long-term performance of 3D models vs 1D/2D models. The water stage data for the 7 locations from OpenFOAM have been saved to .mat format and can be found in /OpenFOAM_WSE. The water stage data from the 1D model are saved in .csv format and can be found in /Mass1_WSE. The water stage from the 2D model is saved as .mat format and can be found in / Mass2_WSEIn addition, the OpenFOAM model outputs the information of hydrostatic and hydrodynamic pressure. They are saved as .mat format under folder /Fig_11/2013_1. As the files are too large, we only uploaded the data for January 2013. The area of different ratio of dynamic pressure to static pressure for all simulation range, i.e., 2013-2015, are saved to .mat format. They can be found in /Fig_11/PA. Further, the data of wall clock time versus the solution time of the OpenFOAM modeling are also saved to .mat format under folder /Fig_13/LogsMat. In summary, the data package contains seven data types, including .txt, .csv, .xlsx, .dat, .stl, .m, and .mat. The former 4 types can be directly open using a text editor or Microsoft Office. The .mat format needs to be read by Matlab. The Matlab source code .m files need to be run with Matlab. The OpenFOAM setups can be visualized in ParaView. The .stl file can be opened in ParaView or Blender. The data in subfolders Fig_1 to Fig_10 and Fig_12 are copied from the aforementioned data folders to generate specific figures for the paper. A readME.txt file is included in each subfolder to further describe how the data in each folder are generated and used to support the paper.Please use the data package's DOI to cite the data package. Please contact yunxiang.chen@pnnl.gov if you need more data related to the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗
DOE OSTI2021

Magmatic, hydrothermal and ore element transfer processes of the southeastern Archean Superior Province implied from electrical resistivity structure

Although magmatic- and metamorphic- derived fluids are widely recognized mineralizing agents, the role of crustal architecture in defining source and sink zones within the middle to lower crust and upper mantle of ancient orogens remains enigmatic. The globally largest and best-preserved Archean greenstone belts lie in the Superior Province, Canada. They provide an ideal location to investigate the influence of igneous construction and subsequent syn-deformational plutonism and metamorphism on the localization of metal-rich melts and fluids throughout the crustal column. Integration of three-dimensional magnetotelluric modelling and seismic reflection sections across the Abitibi subprovince reveals details of a 'whole-of-crust' magmatic and hydrothermal system. East-west low resistivity structures broadly underlie the surface traces of the major deformation zones that are host to significant gold endowment (>200 Moz), while mid-crustal cross trends suggest mineralized fluids flowed along 'pipes' within fault planes. Most low resistivity structures are inferred to represent domains containing interconnected zones of graphite and/or sulfide. Furthermore, these delineate relict mantle source/transit domains and crustal pathways enriched by the flow of magmas or metamorphic fluids genetically related to a late-stage pulse of ore-bearing magmatism, possibly as a result of slab break-off or delamination. Using the combined MT and seismic data, we develop a 3-D crustal-scale model which highlights how evolving orogenic architecture-controlled mass transfer and metallogenic processes developed.

58 GEOSCIENCES↗
DOE OSTI2024

Development of a high-temperature Inconel 625 heat exchanger by model design and binder jetting additive manufacturing

A nickel-based Inconel 625 superalloy heat exchanger for high-temperature applications was developed via binder jetting additive manufacturing. The material properties were characterized first on printed and sintered parts. Two sintering temperatures were used to investigate the effects of the temperature on densification and microstructure. For channel geometry design, the heat transfer capabilities of the heat exchanger were optimized for the cross-section geometry of fluid flow channels in a counterflow configuration, and a stress analysis was conducted to investigate the effects of channel geometry. After headers were incorporated into the heat exchanger as a one-piece component, a prototype was printed. A complete depowdering was achieved via compressed air blowing, and sintering was performed with the developed profile. Ultimately, a one-piece heat exchanger with a nearly full density was obtained. Heat transfer tests were performed on this unit, and the results were compared with those from the simulations. In this study, systematic processes were developed for an additively manufactured Inconel-based heat exchanger for high-pressure, high-temperature heat transfer applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

One piece ceramic heat exchanger for concentrating solar power electric plants

Using additive manufacturing techniques, a ceramic heat exchanger was optimized for a concentrating solar power (CSP) electric power plant with a corrosive molten salt at atmospheric pressure as the solar field heat transfer fluid and supercritical carbon dioxide at 200 bar as the power block fluid. Temperatures ranged from 750 °C to 540 °C. The plant heat exchanger was optimized for heat transfer and stress using 1 m³ modules arranged in parallel. Heat transfer was optimized using COMSOL Multiphysics software resulting in a unique semi-elliptical cross section for the fluid flow channels in a counter-flow configuration. Counter-flow was made possible through additive manufacturing of the fluid headers as part of each modular heat exchanger thus producing the entire module as a single piece and minimizing cost. Using reasonable parametric values, a heat exchanger module was optimized to 0.5 MW/m³, and a section of it was manufactured including the fluid headers. Here, parametric studies showed the potential for a module to reach 3.5 MW/m³.

14 SOLAR ENERGY↗
DOE OSTI2023

Evolution of LIBS technology to mobile instrumentation for expediting firearm-related investigations at the laboratory and the crime scene

Gunshot residue (GSR) is one of the few forensic disciplines that lack accurate screening techniques. This study proposes using a mobile LIBS instrument to detect inorganic GSR and compares performance to a previously validated laboratory instrument. The mobile LIBS is designed with advanced configurations specifically for on-site GSR analysis, including a CMOS detector and a sampling chamber that holds up to six typical GSR collection devices with separate gas flow ports to prevent cross-contamination. A significant novelty of the portable instrument is its image magnification, which allows quick searching and visualization of GSR particle morphology. The single-particle imaging and elemental composition capability is one of a kind and offers superior confirmatory features for GSR. The mobile LIBS performance was evaluated for residues collected from the hands of shooters (100 samples) and non-shooters (200 background samples), analyzed sequentially by the mobile instrument and then the laboratory instrument. Accuracies better than 98.8% were obtained by both instruments, demonstrating their suitability for trace IGSR detection from skin specimens. Implementation of this methodology is anticipated to drastically speed up response times (i.e., from several hours per sample by standard SEM-EDS practice to a few minutes by LIBS). The screening methods can be easily incorporated into workflows to improve decision-making processes at the crime scene and laboratory settings, reduce backlogs, and improve case management.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Additive manufacturing and testing of a ceramic heat exchanger for high-temperature and high-pressure applications for concentrating solar power

Heat exchangers with excellent corrosion and oxidation resistances are essential for the next-generation concentrating solar power (CSP) plants using a molten salt heat transfer fluid and a supercritical CO 2 Brayton power cycle that can operate at temperatures > 700 °C for higher efficiencies. Techniques were developed for additively manufacturing ceramic materials for applications involving high temperature, high pressure, and high corrosion resistance as needed for the CSP application. Based on a previous ceramic heat exchanger design with its cross-section geometry of the flow channels optimized for heat transfer and mechanical strength, lab-scale prototype heat exchangers, with integrated headers and incorporated flow channel dimensional compensations, were fabricated by using the binder jetting process. The printed prototype heat exchangers were successfully densified through the processes of liquid polymer infiltration and pyrolysis. In addition, thermophysical properties of the densified silicon carbide parts were measured to provide necessary design information. Experimental heat transfer testing of the lab-scale prototype heat exchanger was conducted, and the experimental data agreed reasonably well with the simulated results. Furthermore, this agreement provides validation of the simulation models for their applicability to the development of a full-scale ceramic heat exchanger.

14 SOLAR ENERGY↗
DOE OSTI2024

On anomalous transport of multi-species plasma associated with the resistive ballooning and resistive drift waves driven turbulence

Anomalous transport of multi-species plasma related to the resistive ballooning and resistive drift wave turbulence is considered in a “cold” ion approximation. It is found that similar to the resistive drift wave turbulence [see A. R. Knyazev and S. I. Krasheninnikov, Phys. Plasmas 31, 012502 (2024); and S. I. Krasheninnikov and R. D. Smirnov, Phys. Plasmas (to be published)] the addition of the ballooning drive does not change the main features of anomalous transport of the multi-species plasma: (i) The transport of all ion species is described as a transport of the passive scalars in the turbulent field of the electrostatic potential and electron density perturbation; (ii) the density of ion species with a larger ratio of the mass to charge has the tendency to the accumulation/depletion in the vortices of plasma flow; and (iii) the cross-field transport of all plasma species (including electrons and ions) is described by the same anomalous transport coefficient.

Physics↗
DOE OSTI2021

Simple, low-cost and accurate data-driven geophysical forecasting with learned kernels

Modelling geophysical processes as low-dimensional dynamical systems and regressing their vector field from data is a promising approach for learning emulators of such systems. We show that when the kernel of these emulators is also learned from data (using kernel flows, a variant of cross-validation), then the resulting data-driven models are not only faster than equation-based models but are easier to train than neural networks such as the long short-term memory neural network. In addition, they are also more accurate and predictive than the latter. When trained on geophysical observational data, for example the weekly averaged global sea-surface temperature, considerable gains are also observed by the proposed technique in comparison with classical partial differential equation-based models in terms of forecast computational cost and accuracy. When trained on publicly available re-analysis data for the daily temperature of the North American continent, we observe significant improvements over classical baselines such as climatology and persistence-based forecast techniques. Although our experiments concern specific examples, the proposed approach is general, and our results support the viability of kernel methods (with learned kernels) for interpretable and computationally efficient geophysical forecasting for a large diversity of processes.

58 GEOSCIENCES↗
DOE OSTI2022

Nanostructured Composite Alloys for Extreme Environments

Accumulative Roll Bonding was used to process nanolayered Cu/Nb and Zr/Nb composites with varying individual layer thickness. Mechanical tests show remarkable strength as the average layer thickness decreases. Relatively high strength is preserved even at elevated temperatures. Despite differences in texture between rolling paths (longitudinally rolled-RL and cross rolled-CR) Cu/Nb, flow stresses between the resultant materials were similar. In addition, thermal conductivity does not change significantly as the interface density increases. The alloys show remarkable resistance to irradiation damage at very high doses due to stability of the interfaces. The results on Zr/Nb indicate that liquid phase chemistry resulting from thermal spikes play an important role and should be taken into account when designing alloys for extreme environments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗