SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Abstract for CRADA between Capture6 CORP and NETL

The National Energy Technology Laboratory (NETL) and Capture6 CORP (Participant) will collaborate to refine Capture6’s transformative air contactor design for optimal direct air capture (DAC) under an awarded project from the Department of Energy’s Direct Air Capture Pre-Commercialization Technology Prize. Current cooling tower designs require exponentially more resources as the concentration of CO 2 in the gas phase decreases, resulting in a low conversion rate of carbonates, a large physical footprint, and high evaporation of water, which consequently limits the energy and water benefits that a DAC process can create. Capture6’s innovative tower design AURA (Absorption Unit Reactor for Ambient air) seeks to improve the mass transfer characteristics of legacy cooling towers and reduces total energy demand by: (1) introducing laminar and turbulent flow to increase mass transfer, (2) generating a solvent gradient to improve CO 2 absorption, (3) leveraging venturi and packed towers to increase the surface for the mass-transfer process, (4) adding charged/cool plates to reduce overall water loss in the system, and (5) utilizing surfactant additives to increase mass transfer. The collaboration will facilitate the deployment of technologies that reduce the total energy required and carbon footprint of solvent-based carbon dioxide removal (CDR) technologies. The anticipated results support NETL’s goal of fostering rapid development and commercialization pathways for technologies that remove CO 2 from the atmosphere.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗

DOE OSTI2020

Characterization and Scaling of Forced Convective Swirl in Sinusoidal Wavy-Plate-Fin Cores of Compact Heat Exchangers

The characterization and scaling of the thermal-hydraulic performance in wavy plate-fin compact heat exchanger cores, based on the understanding of the physical phenomena and heat transfer enhancement mechanism is delineated. Experimental data are presented for forced convection in air (Pr = 0.71) with flow rates in the range 50 ≤ Re ≤ 4000. A variety of wavy-fin cores that span viable applications, with geometrical attributes described by the cross section aspect ratio α (ratio of fin spacing to height), fin corrugation aspect ratio γ (ratio of 2× corrugation amplitude to wave pitch), and fin spacing ratio ζ (ratio of fin spacing to wave pitch), are considered. To characterize and correlate the vortex-flow mixing in interfin spaces, a Swirl number Sw is introduced from the balance of viscous, inertial and centrifugal forces. It is shown that the laminar, transitional and turbulent flow regimes can be explicitly identified by this Swirl number. Based on the experimental results and extended analysis, new correlations for Fanning friction factor f and Colburn factor j are developed with Sw, α, γ, and ζ as scaling parameters. Requisite expressions are devised by a superposition of both enhancement components due to the corrugated surface area enlargement and induced swirl flow field, and they are combined to cover the laminar, transitional and turbulent regimes by the method of asymptotic matching. The resulting generalized correlations are further shown to predict all available experimental data for f and j factors to within ±20% and ±15%, respectively.

42 ENGINEERING↗

DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗

DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗

DOE OSTI2024

Enhancing Lattice Kinetic Schemes for Fluid Dynamics with Lattice-Equivariant Neural Networks

A new class of equivariant neural networks is presented, hereby dubbed lattice-equivariant neural networks (LENNs), designed to satisfy local symmetries of a lattice structure. The approach develops within a recently introduced framework aimed at learning neural network-based surrogate models’ lattice Boltzmann collision operators. Whenever neural networks are employed to model physical systems, respecting symmetries and equivariance properties has been shown to be key for accuracy, numerical stability, and performance. Here, hinging on ideas from group representation theory, trainable layers are defined whose algebraic structure is equivariant with respect to the symmetries of the lattice cell. In this work, the presented method naturally allows for efficient implementations, in terms of both memory usage and computational costs, supporting scalable training/testing for lattices in two spatial dimensions and higher (in which the size of symmetry group grows). The approach is validated and tested considering 2D and 3D flowing dynamics, both in laminar and turbulent regimes. It is compared with group-averaged-based symmetric networks and with plain, nonsymmetric, networks, showing how the presented approach unlocks the (a posteriori) accuracy and training stability of the former models and the train/inference speed of the latter networks. (LENNs are about one order of magnitude faster than group-averaged networks in 3D.) The work in this paper opens toward practical use of machine learning-augmented lattice Boltzmann CFD in real-world simulations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Large eddy simulation of lean blow-off in swirl-stabilized flame with the subgrid dissipation concept

The lean blow-off mechanism of the premixed swirl flame is numerically investigated by large eddy simulation (LES) with the subgrid dissipation concept (SDC) combustion model. Three simulated cases cover stable, near blow-off, and transient conditions. Compared with the experiment, the LES-SDC approach captures the flow and combustion features for stable and near blow-off conditions. More importantly, the predictions of the blowoff procedure and duration agree satisfactorily with the experiment, indicating that the LES-SDC approach is a promising tool for predicting strong, unsteady turbulent combustion processes. Further, the numerical results are used to investigate the blow-off mechanism. Two stages in the blow-off procedure are specified. The first is the necking and extinction of the downstream flame surface, and the second is the shrinking of the upstream flame surface. The blow-off mechanism is well explained by the theory of stretched flame extinction. At the end of the recirculation zone, the large negative radial velocity pushes the flame to the central line. The combustion process here can be abstracted as the stretched counter-flame of the reactant-reactant configuration. The excessive flame stretch dominates the flame extinction and triggers the blow-off event. The upstream flame resists the intense stretch with the help of hot product recirculation, and the combustion here can be idealized as the counter-flame of the reactant-product configuration. The alignment of the temperature gradient and flow velocity, together with the excessive stretch, clearly indicates the tendency of flame local extinction. A Damköhler number-based blow-off criterion is raised from the mean flow strain rate and laminar flame bulk extinction strain rate.

42 ENGINEERING↗

DOE OSTI2022

Numerical simulation of evaporating wavy falling liquid films in laminar gas streams

Numerical simulations are performed to investigate the interfacial heat and mass transfer from evaporating wavy falling liquid films in interaction with laminar gas streams. The OpenFOAM solver has been used to conduct the simulations where the liquid-gas interface is resolved using the Volume of Fluid method of solving for three phases, i.e., liquid, vapor, and air. The configuration considered is a falling liquid (water) film on a heated vertical plate with a confined laminar moist air (gas) flow that is either (a) co-current or (b) counter-current to the downward liquid flow. The evaporation at the liquid-gas interface is driven by the interfacial gradient of the vapor mass fraction. Interfacial waves are triggered using a monochromatic forcing disturbance that leads to sinusoidal or solitary waves forming at the liquid-gas interface under respective forcing frequencies. Further, the numerical model is validated with the available experimental data. The results show nearly a 15% enhancement in time-averaged Sherwood number (Sh) due to film waviness (sinusoidal or solitary) at the lower volumetric gas flow rate, Q g = +50 (co-current) and Q g = -50 (counter-current). This enhancement in the Sh for both the waves further increases by 11% with Q g = +800 and 196% with Q g = -800. A closer examination of the mass transfer process over a wave demonstrates that with Q g = +50, the concentration of the gas side streamlines at the trough locations of the wave leads to higher values of Sh at these locations. However, with Q g = +800, although the overall Sh increases, vortices appear at the wave trough locations, leading to a corresponding decrease in the local Sh values. Correlations are proposed for predicting Sh under co-current and counter-current gas flow effects.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Chemical influence of ammonia on suppressing soot formation pathways

There is interest in mixing ammonia (NH3) with hydrocarbons, as this decreases overall CO2 emissions from the fuel mixture when burned. Co-firing NH3 with hydrocarbons like methane (CH4) also enables stable combustion over a larger set of operating conditions in engines and gas turbines than firing only NH3 permits. When NH3 is co-fired with hydrocarbons, the influence of NH3 on soot formation chemistry becomes relevant. In this work, soot concentrations in laminar nonpremixed co-flow flames of NH3/CH4 and N2/CH4 were analyzed using color-ratio pyrometry. The concentrations of combustion species (CH4, CO2, H2O, O2, C2H2, C6H6, HCN) across the different flames were also determined using on-line electron-impact mass spectrometry. The flames were also interrogated with 2D CFD simulations, which employ a detailed chemical mechanism tailored to capturing fuel-nitrogen chemistry. Addition of NH3 to CH4 flames was found to decrease soot concentrations more significantly than N2-addition (Fig. 1). The NH3/CH4 flames were found to have higher flame temperatures than the N2/CH4 flames, implying the suppression of soot afforded by NH3 is due to a chemical influence of NH3, as opposed to thermal or dilution effects. Based on experimental and simulated soot volume fractions and species concentrations, soot formation pathways affected by NH3 are postulated.

10 SYNTHETIC FUELS↗

DOE OSTI2022

Noninvasive Imaging of a Liquid Jet

Liquid jets are found in many applications, from printing to manufacturing to entertainment. Herein this study uses three different noninvasive imaging modalities to compare resulting images of a liquid jet operating at three Reynolds numbers that cover laminar, transitional, and turbulent flow. Selected measurement quantities from each image type are also compared. High-speed backlit (BL) imaging is a simple imaging technique found in many laboratories, and this is compared to two high-speed X-ray imaging techniques, white beam (WB) imaging and focused beam (FB) radiography. BL imaging can provide a wide field of view and is easy to implement, but it only shows the presence or absence of liquid. WB imaging can show detailed contours on the surface of the liquid jet, but the imaging region is much smaller. FB radiography produces a point-source measurement and can provide the quantitative, instantaneous local liquid path length, termed the equivalent path length (EPL). All three techniques provide similar measures of jet thickness, with the FB measurements having less variation. FB measurements can also provide detailed cross sections of the average liquid jet thickness at high spatial resolutions.

42 ENGINEERING↗

DOE OSTI2023

FY23 Progress Report on Viscosity and Thermal Conductivity Measurements of Molten Salts

As presented in this report, thermal conductivity and viscosity measurements were performed on key chloride pseudo-binary molten salt systems of relevance to molten salt reactor developers. Thermal conductivity measurements were conducted with a variable gap technique, in which a known heat flux is driven across a molten salt specimen and the temperature difference is measured, allowing calculation of the thermal conductivity. This is achieved by establishing a small gap between the bottom of a cylindrical inner containment, which houses electrical heating elements, and an outer containment, which houses cooling channels; the gap size can be varied by compression of a formed bellows. A new calibration scheme was developed herein, involving a correction factor to the heat flux based on He measurements at various temperatures. Furthermore, the data processing methodology was improved to minimize the impact of radiative heat transfer in the calculation of salt specimen thermal conductivity from the temperature difference measurements. Viscosity measurements were conducted with a rolling ball viscometer, in which a ball rolls some known distance in an angled tubular crucible, and the terminal velocity can be used to calculate the viscosity of the salt. The measurement can be performed in a quartz crucible, with which a standard camera can be used to track the ball, or in a metal crucible, with which x-ray radiography is required to track the ball. A new custom x-ray system was made and dedicated to the rolling ball viscometer to enable high throughput automated measurements with salts which require containment with metal. Both systems have been integrated with a ventilation stack which allows for off-gassing of radioactive material, enabling future measurements with U-bearing salts. The thermal conductivity measurements performed herein were with NaCl-KCl (44 mol% NaCl). This salt system was measured in the previous fiscal year, however the thermal conductivity values obtained were comparatively low, and so the system was remeasured with the aforementioned calibration scheme and improved post-processing techniques. The newly obtained thermal conductivity values for NaCl-KCl (44 mol% NaCl) indicate good agreement with kinetic theory and ab-initio models (within 5–10 %). The viscosity measurements performed herein were with three different compositions of NaCl-KCl: 75, 50, and 25 mol% NaCl; a new calibration scheme was applied to account for variable flow effects in the laminar regime. The results show reasonable agreement with literature (5–20 %, depending on the temperature, composition, and study); however, literature values are likely higher than true values based on pure end-member measurements performed in the comparative studies. The measurements conducted herein do show a trend such that viscosity increases with increasing NaCl concentration, which agrees well with one of two comparative studies. In general, the measurements conducted herein gives confidence in the capability to use these systems to accurately measure thermal conductivity and viscosity of actinide-bearing salts within the next fiscal year.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Multi-fidelity information fusion with concatenated neural networks

Recently, computational modeling has shifted towards the use of statistical inference, deep learning, and other data-driven modeling frameworks. Although this shift in modeling holds promise in many applications like design optimization and real-time control by lowering the computational burden, training deep learning models needs a huge amount of data. This big data is not always available for scientific problems and leads to poorly generalizable data-driven models. This gap can be furnished by leveraging information from physics-based models. Exploiting prior knowledge about the problem at hand, this study puts forth a physics-guided machine learning (PGML) approach to build more tailored, effective, and efficient surrogate models. For our analysis, without losing its generalizability and modularity, we focus on the development of predictive models for laminar and turbulent boundary layer flows. In particular, we combine the self-similarity solution and power-law velocity profile (low-fidelity models) with the noisy data obtained either from experiments or computational fluid dynamics simulations (high-fidelity models) through a concatenated neural network. We illustrate how the knowledge from these simplified models results in reducing uncertainties associated with deep learning models applied to boundary layer flow prediction problems. The proposed multi-fidelity information fusion framework produces physically consistent models that attempt to achieve better generalization than data-driven models obtained purely based on data. While we demonstrate our framework for a problem relevant to fluid mechanics, its workflow and principles can be adopted for many scientific problems where empirical, analytical, or simplified models are prevalent. In line with grand demands in novel PGML principles, this work builds a bridge between extensive physics-based theories and data-driven modeling paradigms and paves the way for using hybrid physics and machine learning modeling approaches for next-generation digital twin technologies.

42 ENGINEERING↗

DOE OSTI2020

Characterization of flowing liquid films as a regenerating plasma mirror for high repetition-rate laser contrast enhancement

In this paper, we characterize a high repetition-rate regenerating plasma mirror produced by the thin film of liquid formed when two laminar streams collide. The use of a flowing liquid film is inexpensive and the interaction surface refreshes automatically, avoiding buildup of on-target debris. The composition of the liquid material and the relative angle of the film-generating nozzles was optimized for this application. Spectra measured in reflection from a water-based plasma mirror showed a blue shift but an optical reflectivity of up to 30%. The thickness of the film was found to be of the order of 2 μm, and the stability of the reflected spot was ≈1 mrad. The reflected beam profile was highly distorted but stable. Further optimization of the nozzles to affect the fluid flow should enable significant improvements in control of the fluid films and increase in the reflectivity of these mirrors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Evaluating Scaling Effects of Steam Condensation Physics-based Model for Small Modular Reactor

This study focused on assessing condensation heat transfer (CHT) models for enhancing the passive safety of small modular reactors (SMRs). It also encompassed a review of the theoretical and semi-theoretical CHT models and correlations. Condensation models such as the degradation factor, heat and mass transfer analogy (HMTA), and boundary layer models are used to study advanced nuclear reactor PCCS. The previous CHT models were developed for small tube geometries, which differ from the containment of SMR. This study examined the scaled relations of a physics-based CHT model for analyzing laminar and turbulent film condensation downward flow conditions inside scaled tube geometries.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

Presentation: Evaluating scaling effects of steam condensation physics-based model for small modular reactors safety

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

HTR solver: An open-source exascale-oriented task-based multi-GPU high-order code for hypersonic aerothermodynamics

In this study, the open-source Hypersonics Task-based Research (HTR) solver for hypersonic aerothermodynamics is described. The physical formulation of the code includes thermochemical effects induced by high temperatures (vibrational excitation and chemical dissociation). The HTR solver uses high-order TENO-based spatial discretization on structured grids and efficient time integrators for stiff systems, is highly scalable in GPU-based supercomputers as a result of its implementation in the Regent/Legion stack, and is designed for direct numerical simulations of canonical hypersonic flows at high Reynolds numbers. Additionally, the performance of the HTR solver is tested with benchmark cases including inviscid vortex advection, low- and high-speed laminar boundary layers, inviscid one-dimensional compressible flows in shock tubes, supersonic turbulent channel flows, and hypersonic transitional boundary layers of both calorically perfect gases and dissociating air.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

On the enhancement of boundary layer skin friction by turbulence: an angular momentum approach

Turbulence enhances the wall shear stress in boundary layers, significantly increasing the drag on streamlined bodies. Other flow features such as free stream pressure gradients and streamwise boundary layer growth also strongly influence the local skin friction. In this paper, an angular momentum integral (AMI) equation is introduced to quantify these effects by representing them as torques that alter the shape of the mean velocity profile. This approach uniquely isolates the skin friction of a Blasius boundary layer in a single term that depends only on the Reynolds number most relevant to the flow's engineering context, so that other torques are interpreted as augmentations relative to the laminar case having the same Reynolds number. The AMI equation for external flows shares this key property with the so-called FIK relation for internal flows (Fukagata et al. , Phys. Fluids , vol. 14, 2002, pp. L73–L76). Without a geometrically imposed boundary layer thickness, the length scale in the Reynolds number for the AMI equation may be chosen freely. After a brief demonstration using Falkner–Skan boundary layers, the AMI equation is applied as a diagnostic tool on four transitional and turbulent boundary layer direct numerical simulation datasets. Regions of negative wall-normal velocity are shown to play a key role in limiting the peak skin friction during the late stages of transition, and the relative strengths of terms in the AMI equation become independent of the transition mechanism a very short distance into the fully turbulent regime. The AMI equation establishes an intuitive, extensible framework for interpreting the impact of turbulence and flow control strategies on boundary layer skin friction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗