Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow physics”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Wickability-optimized textured liquid-desiccant air dehumidifiers for independent moisture management in energy-efficient buildings

Liquid-desiccant-based air conditioning systems are envisioned to enable independent humidity management, thereby improving the energy efficiency of future buildings. Existing liquid-desiccant-based air conditioning concepts, however, suffer from a poor liquid flow distribution deteriorating moisture removal rate. They are consequently flooded with the liquid-desiccant solution, which significantly degrades the energy efficiency of the dehumidification process. Here, in this work, the capillary forces and wickability effect of textured air dehumidifier surfaces are altered to minimize the liquid-desiccant flow rate of the fully wetted state, thereby transforming the physics of interfacial desiccant flow distribution. Consequently, the wickability-optimized air dehumidifier surface maximizes both moisture removal rate and dehumidification energy efficiency. It was interestingly found that the length scale of a textured air dehumidifier surface concept is optimized at an intermediate pattern density. Dry solid-air menisci appear at length scales exceeding the optimum pattern distance while the effective liquid–air interfacial area is reduced at smaller length scales, both of which degrade the moisture removal rate. At the optimum pattern density, the effective liquid–air interfacial area increases with the solution flow rate, thereby increasing the dehumidification rate. At a water vapor pressure potential of 3 kPa and a solution flow rate of 2.8 g/s, experimental results indicated a moisture removal rate of 0.16 g/m2-s for a textured surface concept with a capillary length scale of 3 mm, a 28% improvement compared with that of smooth-plate dehumidifier surfaces. A high moisture removal rate of the textured surface at a low desiccant flow rate led to a high thermal efficiency of 0.75 at a water vapor pressure potential of 5.6 kPa and a LiBr flow rate of 2.8 g/s. The insights gained from the present study accelerate the development of advanced textured surface concepts for next-generation liquid-desiccant-based air dehumidification systems offering independent humidity management for future energy-efficient buildings.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Flow field characteristics of a 127-pin rod bundle with hexagonal spacer grids

A determination of nominal flow phenomena in liquid metal fast reactor (LMFR) fuel assemblies is critical toward generation-IV reactor development. Axially positioned spacer grids are used to maintain the geometry of hexagonal rod bundles and simultaneously introduce perturbations in the flow. Three-dimensional (3D) printed asymmetric honeycomb spacer grids were installed in a prototypical 127-pin LMFR fuel assembly model to study complex fluid dynamics interactions induced by the spacer grid and rods. To characterize flow dynamics in this intricate geometry, time-resolved particle image velocimetry (TR-PIV) using the matched-index-of-refraction method was employed to obtain non-intrusive velocity measurements for three axial planes (one near-wall and two interior planes) at a Reynolds number of 6000. The statistical TR-PIV results compared sub-channel-dependent normalized time-averaged velocity, velocity fluctuations, Reynolds stress, vorticity, and turbulence kinetic energy distributions. TR-PIV line profiles characterized downstream spacer grid flow dynamics. Two-point spatial and spatial–temporal cross-correlation fields revealed local coherent structures and quantified convection velocities of traveling vortices. Spatial–temporal decomposition using dynamic mode decomposition (DMD) applied to the near-wall vorticity fields extracted turbulent structures and flow instabilities in the wake region of the spacer grid, along with their decay and frequency rates. Reduced-order velocity fields from DMD reconstructions identified the most energy-containing coherent structures persistent in the near-wall region. This research provides experimental data sets and analyses of flow behavior in rod bundles with hexagonal spacer grids. In conclusion, the results are critical toward LMFR design and geometry optimization, crucial for the validation of computational fluid dynamics and reduced-order flow models.

42 ENGINEERING↗
DOE OSTI2018

Unified first-principles theory of thermal properties of insulators

The conventional first-principles theory for the thermal and thermodynamic properties of insulators is based on the perturbative treatment of the anharmonicity of crystal bonds. While this theory has been a successful predictive tool for strongly bonded solids such as diamond and silicon, here we show that it fails dramatically for strongly anharmonic (weakly bonded) materials, and that the conventional quasiparticle picture breaks down at relatively low temperatures. To address this failure, we present a unified first-principles theory of the thermodynamic and thermal properties of insulators that captures multiple thermal properties within the same framework across the full range of anharmonicity from strongly bonded to weakly bonded insulators. This theory features a new phonon renormalization approach derived from many-body physics that creates well-defined quasiparticles even at relatively high temperatures, and it accurately captures the effects of strongly anharmonic bonds on phonons and thermal transport. Using a prototypical strongly anharmonic material, sodium chloride (NaCl), as an example, we demonstrate that our new first-principles framework simultaneously captures the apparently contradictory experimental observations of large thermal expansion and low thermal conductivity of NaCl on the one hand, and anomalously weak temperature dependence of phonon modes on the other, while the conventional theory fails in all three cases. Here we demonstrate that four-phonon scattering due to higher-order anharmonicity significantly lowers the thermal conductivity of NaCl and is required for a proper comparison to experiment. Furthermore, we show that our renormalization framework, along with four-phonon scattering, also successfully predicts the measured phonon frequencies and thermal properties of a weakly anharmonic material, diamond, indicating universal applicability for thermal properties of insulators. Our work gives new insights into the physics of heat flow in solids, and presents a computationally efficient and rigorous framework that captures the thermal and thermodynamic properties of both weakly and strongly bonded insulators simultaneously.

36 MATERIALS SCIENCE↗
NASA NTRS1984

The computation of mesh configurations for three-dimensional flow analysis

A geometry modeling method was developed to define three-dimensional mesh configurations, which were expressed by the set of surface intersections with a coordinate mesh. These configurations define bodies for transonic potential-flow computation in simple physical coordinates. Various configurations were prepared for 3-D potential flow analysis, including many bare inlet models in cylindrical coordinates. Trimmed surface models were used to prepare nacelle installation configurations. A variety of ducts and mixers were represented. Analytically defined mesh/surface intersections were compared with similar data obtained from the corresponding surface model. Parameters for an analytic surface definition from mesh-surface intersections obtained from a model of the surface were produced. In both studies the modeling accuracy was shown to be compatible with the flow analysis requirements. This technique fulfills the requirements of the associated three-dimensional potential-flow analysis. It preserves the flow analysis method's superior flexibility for adapting to unusual geometries. The geometry modeling technique is supported by a sufficient set of tools for effective use by general users, and sufficient accuracy was demonstrated.

Gibson, S. G.↗
NASA NTRS1987

Application of advanced diagnostics to airblast injector flows

This effort is concerned with the application of both conventional laser velocimetry and phase Doppler anemometry to the flow produced by an airblast nozzle. The emphasis is placed on the acquisition of data using actual engine injector/swirler components at (noncombusting) conditions simulating those encountered in the engine. The objective of the effort was to test the applicability of the instrumentation to real injector flows, to develop information on the behavior of injectors at high flow, and to provide data useful in the development of physical models of injector flows.

Mcvey, John B.↗
NASA NTRS1995

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗
NASA NTRS1982

A new theory of sources of Birkeland currents

An approach to collisionless plasma shows the existence of current orthogonal to B along the low latitude boundary layer of the magnetosphere driven by electric field which is orthogonal to both B and the layer. In this case relationship P perpendicular + (sq B/8 pi) - (sq epsilon E perpendicular/8 pi) = constant holds on a line orthogonal to B and the layer, where epsilon is the dielectric constant of the plasma for electric fields orthogonal to B. Across the geomagnetic tail there flows a current in the direction of the dawn-dusk electric field, and in this case a relationship P perpendicular + (sq B/8 pi) + (sq epsilon E perpendicular/8 pi = constant, holds aong a line orthogonal to E and B. Divergence of both these currents is shown to be a source of Birkeland currents. Some of the boundary layer current is continuous with current across the tail. Electric currents of physically similar origin flow in interplanetary space, and when the magnetosphere interrupts them, additional Birkeland currents are driven.

Cole, K. D.↗
NASA NTRS1984

A new theory of sources of Birkeland currents

An approach to collisionless plasma shows the existence of current orthogonal to B along the low latitude boundary layer of the magnetosphere driven by electric field which is orthogonal to both B and the layer. In this case relationship P perpendicular + (sq B/8 pi) - (sq epsilon E perpendicular/8 pi) = constant holds on a line orthogonal to B and the layer, where epsilon is the dielectric constant of the plasma for electric fields orthogonal to B. Across the geomagnetic tail there flows a current in the direction of the dawn-dusk electric field, and in this case a relationship P perpendicular + (sq B/8 pi) + (sq epsilon E perpendicular/8 pi) = constant, holds along a line orthogonal to E and B. Divergence of both these currents is shown to be a source of Birkeland currents. Some of the boundary layer current is continuous with current across the tail. Electric currents of physically similar origin flow in interplanetary space, and when the magnetosphere interrupts them, additional Birkeland currents are driven.

Cole, K. D.↗
NASA NTRS1991

An evaluation of three two-dimensional computational fluid dynamics codes including low Reynolds numbers and transonic Mach numbers

Full-potential, Euler, and Navier-Stokes computational fluid dynamics (CFD) codes were evaluated for use in analyzing the flow field about airfoils sections operating at Mach numbers from 0.20 to 0.60 and Reynolds numbers from 500,000 to 2,000,000. The potential code (LBAUER) includes weakly coupled integral boundary layer equations for laminar and turbulent flow with simple transition and separation models. The Navier-Stokes code (ARC2D) uses the thin-layer formulation of the Reynolds-averaged equations with an algebraic turbulence model. The Euler code (ISES) includes strongly coupled integral boundary layer equations and advanced transition and separation calculations with the capability to model laminar separation bubbles and limited zones of turbulent separation. The best experiment/CFD correlation was obtained with the Euler code because its boundary layer equations model the physics of the flow better than the other two codes. An unusual reversal of boundary layer separation with increasing angle of attack, following initial shock formation on the upper surface of the airfoil, was found in the experiment data. This phenomenon was not predicted by the CFD codes evaluated.

Hicks, Raymond M.↗
DOE OSTI2020

A Method to Represent a Well in a Three‐Dimensional Discrete Fracture Network Model

Abstract In discrete fracture network (DFN) modeling, fractures are randomly generated and placed in the model domain. The rock matrix is considered impermeable. Small fractures and isolated fractures are often ignored to reduce computational expense. As a result, the rock matrix between fractures could be large and intersections may not be found between a well introduced in the model and the hydraulically connected fracture networks (fracture backbones). To overcome this issue, this study developed a method to conceptualize a well in a three‐dimensional (3D) DFN using two orthogonal rectangular fractures oriented along the well's axis. Six parameters were introduced to parameterize the well screen and skin zone, and to control the connectivity between the well and the fracture backbones. The two orthogonal fractures were discretized using a high‐resolution mesh to improve the quality of flow and transport simulations around and along the well. The method was successfully implemented within dfnWorks 2.0 (Hyman et al. 2015) to incorporate a well in a 3D DFN and to track particles leaving an injection well and migrating to a pumping well. Verification of the method against MODFLOW/MODPATH found a perfect match in simulated hydraulic head and particle tracking. Using three examples, the study showed that the method ensured the connectivity between wells and fracture backbones, and honored the physical processes of flow and transport along and around wells in DFNs. Recommendations are given for estimating the values of the six introduced well parameters in a real‐world case study.

Pham, Hai (ORCID:0000000310164394)↗
DOE OSTI2019

Studying the effect of phenyl group content on the physical properties of un-crosslinked PDMS (Progress Summary)

Over the last couple of years, the emphasis on our A&L efforts has been on developing a thermal-age-aware constitutive model for foam cushions. From our studies so far, it appears that an $N = 2$ Ogden hyperfoam model within the framework of a two-network Tobolsky scheme is appropriate to describe the age-aware stress-strain response of both stochastic and additively manufactured (AM) foams. The models have so far been fitted only to a limited number of silicone-based materials systems, i.e., M97 and M9763 stochastic foams and AM foams made of SE1700 rubber. As we explore new and different rubber formulations for improved structural as well as A&L properties, we find it important to develop a deeper understanding of the effect of functional groups on the physical, mechanical, and flow properties of the constituent elastomer. To this end, in Q2 FY19, we have begun a project on multiscale modeling of polymers that realistically incorporates from bond-level dynamics to terminal dynamics of whole chains (e.g., tube dynamics). In such mesoscale simulations the polymer tube dynamics is modeled through coarse-graining of monomers to mesoscale “beads”, with the bead-bead potential being systematically fitted to accurate atomic-level forcefields. The aim of this collaboration is to develop a fully physical constitutive model that predicts the mechanical response of polymers of our interest based on structure and modeled dynamics alone.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Pre-equilibrium dynamics and heavy-ion observables

To bracket the importance of the pre-equilibrium stage on relativistic heavy-ion collision observables, we compare simulations where it is modeled by either free-streaming partons or fluid dynamics. These cases implement the assumptions of extremely weak vs. extremely strong coupling in the initial collision stage. Accounting for flow generated in the pre-equilibrium stage, we study the sensitivity of radial, elliptic and triangular flow on the switching time when the hydrodynamic description becomes valid. Using the hybrid code iEBE-VISHNU [1] we perform a multi-parameter search, constrained by particle ratios, integrated elliptic and triangular charged hadron flow, the mean transverse momenta of pions, kaons and protons, and the second moment of the proton transverse momentum spectrum, to identify optimized values for the switching time \tau_s from pre-equilibrium to hydrodynamics, the specific shear viscosity eta/s, the normalization factor of the temperature-dependent specific bulk viscosity (zeta/s)(T), and the switching temperature T_sw from viscous hydrodynamics to the hadron cascade UrQMD. With the optimized parameters, we predict and compare with experiment the p_T -distributions of pi, K, p, Lambda, Xi and Omega yields and their elliptic flow coefficients, focusing specifically on the mass-ordering of the elliptic flow for protons and Lambda hyperons which is incorrectly described by VISHNU without pre-equilibrium flow.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Measuring differential flow angle fluctuations in relativistic nuclear collisions

Recently studies of the differential nature of the flow angle fluctuations, known as event plane angular decorrelation, indicated that measurements that assume a common symmetry plane may need to consider the flow angle fluctuations effect. Using the HIJING and AMPT models, it is shown that the flow angle fluctuations measurements, obtained with the two- and three-subevents correlation method, could have significant nonflow effects associated with long- and short-range nonflow correlations. The current study demonstrates the four-subevents cumulant methods' ability to reduce the nonflow effects. Finally, it is further argued that the measurements using the four-subevents correlation method can be used to provide accurate quantification of the differential flow angle fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRSDate not supplied

Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations

Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.

Boundary Layer Transition↗
DOE OSTI2016

Unsteady fluid flow in a slightly curved pipe: A comparative study of a matched asymptotic expansions solution with a single analytical solution

The present work is motivated by the fact that blood flow in the aorta and the main arteries is governed by large finite values of the Womersley number α and for such values of α there is not any analytical solution in the literature. The existing numerical solutions, although accurate, give limited information about the factors that affect the flow, whereas an analytical approach has an advantage in that it can provide physical insight to the flow mechanism. Having this in mind, we seek analytical solution to the equations of the fluid flow driven by a sinusoidal pressure gradient in a slightly curved pipe of circular cross section when the Womersley number varies from small finite to infinite values. Initially the equations of motion are expanded in terms of the curvature ratio δ and the resulting linearized equations are solved analytically in two ways. In the first, we match the solution for the main core to that for the Stokes boundary layer. This solution is valid for very large values of α. In the second, we derive a straightforward single solution valid to the entire flow region and for 8 ≤ α < ∞, a range which includes the values of α that refer to the physiological flows. Each solution contains expressions for the axial velocity, the stream function, and the wall stresses and is compared to the analogous forms presented in other studies. The two solutions give identical results to each other regarding the axial flow but differ in the secondary flow and the circumferential wall stress, due to the approximations employed in the matched asymptotic expansion process. The results on the stream function from the second solution are in agreement with analogous results from other numerical solutions. The second solution predicts that the atherosclerotic plaques may develop in any location around the cross section of the aortic wall unlike to the prescribed locations predicted by the first solution. In addition, it gives circumferential wall stresses augmented by approximately 100% with respect to the matched asymptotic expansions, a factor that may contribute jointly with other pathological factors to the faster aging of the arterial system and the possible malfunction of the aorta.

42 ENGINEERING↗
DOE OSTI2024

Predicting turbulent wake flow of marine hydrokinetic turbine arrays in large-scale waterways via physics-enhanced convolutional neural networks

We present a physics-enhanced convolutional neural network (PECNN) algorithm for reconstructing the mean flow and turbulence statistics in the wake of marine hydrokinetic (MHK) turbine arrays installed in large-scale meandering rivers. The algorithm embeds the mass and momentum conservation equations into the loss function of the PECNN algorithm to improve the physical realism of the reconstructed flow fields. The PECNN is trained using large eddy simulation (LES) results of the wake flow of a single row of turbines in a virtual meandering river. Subsequently, the trained PECNN is applied to predict the wake flow of MHK turbines with arrangements and positionings different than those considered during the training process. The PECNN predictions are validated using the results of separately performed LES. The results show that the PECNN algorithm can accurately predict the wake flow of MHK turbine farms at a small fraction of the cost of LES. The PECNN can improve the accuracy by around 1% and reduce the physical constraint indices by around 50% compared to the CNN without physical constraints. This work underscores the potential of PECNN to develop reduced-order models for control co-design and optimization of MHK turbine arrays in natural riverine environments.

Mechanics↗