Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Prediction of vortex breakdown on a delta wing

Recent studies of leading-edge vortex flows with computational fluid dynamics codes using Euler or Navier-Stokes formulations have shown fair agreement with experimental data. These studies have concentrated on simulating the flowfields associated with a sharp-edged flat plate 70 deg delta wing at angles of attack where vortex breakdown or burst is observed over the wing. There are, however, a number of discrepancies between the experimental data and the computed flowfields. The location of vortex breakdown in the computational solutions is seen to differ from the experimental data and to vary with changes in the computational grid and freestream Mach number. There also remain issues as to the validity of steady-state computations for cases which contain regions of unsteady flow, such as in the post-breakdown regions. As a partial response to these questions, a number of laminar Navier-Stokes solutions were examined for the 70 deg delta wing. The computed solutions are compared with an experimental database obtained at low subsonic speeds. The convergence of forces, moments and vortex breakdown locations are also analyzed to determine if the computed flowfields actually reach steady-state conditions.

Agrawal, S.↗
NASA NTRS2019

Total Temperature Measurements in Icing Cloud Flows Using a Rearward Facing Probe

This paper reports on temperature and humidity measurements from a series of ice-crystal icing tunnel experiments conducted in June 2018 at the Propulsion Systems Laboratory at the NASA Glenn Research Center. The tests were fundamental in nature and were aimed at investigating the icing processes on a two-dimensional NACA0012 airfoil subjected to artificially generated icing clouds. Prior to the tests on the airfoil, a suite of instruments, including total temperature and humidity probes, were used to characterize the thermodynamic flow and icing cloud conditions of the facility. Two different total temperature probes were used in these tests which included a custom designed rearward facing probe and a commercial self-heating total temperature probe. The rearward facing probe, the main total temperature probe, was designed to reduce and mitigate the contaminating effects of icing and ingestion of ice crystals and water droplets at the probe's inlet. The probe also serves as an air-sample inlet for a light absorption based humidity measurement. The paper includes a section which discusses total temperature and humidity measurement considerations, and another section which provides an analysis of the main probe's performance characteristics. A computational fluid dynamic model of the flow around the probe was also conducted to gain insight into the trajectory of the flow entering the probe inlet. The experiments included a series of tests in which the relative humidity of the facility flow was swept through with increasingly larger values. The data showed that the rearward facing probe can reasonably capture the flow's total temperature and humidity under mild to moderate icing conditions but produces anomalous results under more intense icing conditions. The experimental data was also compared to an in-house developed thermodynamic model which takes into account the interaction of the main flow with the icing cloud. Comparison to the thermodynamic model showed that the rearward facing probe measured the predicted trends.

Agui, Juan H.↗
NASA NTRS1993

ASRM Multi-Port Igniter Flow Field Analysis

The Advanced Solid Rocket Motor (ASRM) program was initiated by NASA in response to the need for a new generation rocket motor capable of providing increased thrust levels over the existing Redesigned Solid Rocket Motor (RSRM) and thus augment the lifting capacity of the space shuttle orbiter. To achieve these higher thrust levels and improve motor reliability, advanced motor design concepts were employed. In the head end of the motor, for instance, the propellent cast has been changed from the conventional annular configuration to a 'multi-slot' configuration in order to increase the burn surface area and guarantee rapid motor ignition. In addition, the igniter itself has been redesigned and currently features 12 exhaust ports in order to channel hot igniter combustion gases into the circumferential propellent slots. Due to the close proximity of the igniter ports to the propellent surfaces, new concerns over possible propellent deformation and erosive burning have arisen. The following documents the effort undertaken using computational fluid dynamics to perform a flow field analysis in the top end of the ASRM motor to determine flow field properties necessary to permit a subsequent propellent fin deformation analysis due to pressure loading and an assessment of the extent of erosive burning.

Kania, Lee↗
DOE OSTI2022

Local Lagrangian reduced-order modeling for the Rayleigh-Taylor instability by solution manifold decomposition

The Rayleigh-Taylor instability is a classical hydrodynamic instability of great interest in various disciplines of science and engineering, including astrophysics, atmospheric sciences and climate, geophysics, and fusion energy. Analytical methods cannot be applied to explain the long-time behavior of the Rayleigh-Taylor instability, and therefore, numerical simulation of the full problem is required. However, in order to capture the growth of amplitude of perturbations accurately, both the spatial and temporal discretizations need to be extremely fine for traditional numerical methods, and long-time simulation may become prohibitively expensive. In this paper, we propose efficient reduced order model techniques to accelerate the simulation of the Rayleigh-Taylor instability in compressible gas dynamics. Here, we introduce a general framework for decomposing the solution manifold to construct the temporal domain partition and temporally-local reduced order model construction with varying Atwood number. We propose two practical approaches in this framework, namely decomposition by physical time and by penetration distance. Numerical results are presented to examine the performance of the proposed approaches.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

A large eddy simulation scheme for turbulent reacting flows

The recent development of the dynamic subgrid-scale (SGS) model has provided a consistent method for generating localized turbulent mixing models and has opened up great possibilities for applying the large eddy simulation (LES) technique to real world problems. Given the fact that the direct numerical simulation (DNS) can not solve for engineering flow problems in the foreseeable future (Reynolds 1989), the LES is certainly an attractive alternative. It seems only natural to bring this new development in SGS modeling to bear on the reacting flows. The major stumbling block for introducing LES to reacting flow problems has been the proper modeling of the reaction source terms. Various models have been proposed, but none of them has a wide range of applicability. For example, some of the models in combustion have been based on the flamelet assumption which is only valid for relatively fast reactions. Some other models have neglected the effects of chemical reactions on the turbulent mixing time scale, which is certainly not valid for fast and non-isothermal reactions. The probability density function (PDF) method can be usefully employed to deal with the modeling of the reaction source terms. In order to fit into the framework of LES, a new PDF, the large eddy PDF (LEPDF), is introduced. This PDF provides an accurate representation for the filtered chemical source terms and can be readily calculated in the simulations. The details of this scheme are described.

Gao, Feng↗
DOE OSTI2021

Simulation and analysis of slip flow of water at hydrophobic silica surfaces of nanometer slit pores

Water flow in nanopores is controlled fundamentally by pore wettability and interfacial water features. In the research herein reported, Molecular dynamics (MD) was used to simulate water flow in nanometer slit pores of silica surfaces with variable wetting characteristics. Compared to Couette flow, the simulated slit-pores of the talc (001) and octadecyltrichlorosilane (OTS) monolayer surfaces exhibited water flow enhancement and a positive slip length at the interface, whereas water flow velocity at the quartz (001) and amorphous silica surfaces was zero, adhering to the no slip condition. More intensive interactions between the interfacial water and hydrophilic silica surfaces were confirmed by MD simulation and sum-frequency vibrational spectroscopy (SFVS) results. A water exclusion zone of approximately 3 Å was found at the hydrophobic silica surfaces by quantitative MD interfacial water analysis and validated by SFVS. This exclusion zone dimension was similar in magnitude to the slip length (4–5 Å) as determined from simulated water flow in hydrophobic silica slit pores. The simulated slip lengths at the hydrophobic surfaces (0.2–1 nm) were significantly less than the experimental slip lengths (>40 nm) reported in the literature. In the case of the experimental results, this difference appears to be due to nanobubbles expected to be present at hydrophobic surfaces.

42 ENGINEERING↗
DOE OSTI2023

Measurement of ion temperature and toroidal flow during magnetic reconnection with a large guide field

Here, ion temperature and toroidal flow along the guide field direction are measured using a new ion tomographic diagnostic on the Magnetic Reconnection eXperiment (MRX) during magnetic reconnection with a guide field strength of about 1.4 and 2.1 times the strength of the reconnecting component. Strong toroidal flows, beyond what has been measured in anti-parallel and lower guide field conditions on MRX, are observed. Sustained ion heating with no discernible structure within the measurement region is also observed. Probe measurements including Langmuir and Mach probe measurements are made to support the tomographic inversion of line-integrated measurements, as well as to provide local measurements of plasma parameters. Measurements of toroidal velocity and ion temperature are supported with time series data. Energy flow into and out of the X-line region is estimated using a guiding center framework and presented in the Appendix of this manuscript, suggesting an outsized role played by parallel electric field in energizing ions. The guiding center approximation is not well satisfied in the region of interest; however, the estimates provide a springboard for future, further experimentation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Flow Strength Measurements of Wrought and AM SS304L via Pressure Shear Plate Impact Experiments

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Self-generated electrokinetic flows from active-charged boundary patterns

We develop a hydrodynamic description of self-generated electrolyte flow in capillaries whose bounding walls feature nonuniform distributions of charge nonuniform active ionic fluxes. The hydrodynamic velocity arising in such a system has components that are forbidden by symmetry in the absence of charge and fluxes. However, when these two boundary mechanisms are simultaneously present, they can lead to a symmetry broken state where steady flows with both unidirectional and circulatory components emerge. We show that these flow states arise when modulated boundary patterns of charge and fluxes are offset by a flux-charge phase difference, which is associated with the separation between sites of their peak densities on the wall. Mismatch in diffusivity of cationic and anionic species can modify the flow states and becomes an enhancing factor when fluxes of both ion species are being produced together at the same site. We demonstrate that this mechanism can be realized with a microfluidic generator that is powered by enzyme-coated patches that catalyze reactants in the solution to produce fluxes of ions. The local ionic elevation or depletion, which disrupts a nonuniform double layer, promotes self-induced gradients yielding persistent body forces to generate bulk fluid motion. Our work quantifies a boundary-driven mechanism behind self-sustained electrolyte flow in confined environments that exists without any external bulk-imposed fields or gradients. It provides a theoretical framework for understanding the combined effect of active and charged boundaries that are relevant in biological or soft matter systems, and can be utilized in electrofluidic and iontronic applications.

active matter↗
NASA NTRS1997

Large-Eddy Simulation of a High Reynolds Number Flow Around a Cylinder Including Aeroacoustic Predictions

The dynamic subgrid-scale model is employed in large-eddy simulations of flow over a cylinder at a Reynolds number, based on the diameter of the cylinder, of 90,000. The Centric SPECTRUM(trademark) finite element solver is used for the analysis. The far field sound pressure is calculated from Lighthill-Curle's equation using the computed fluctuating pressure at the surface of the cylinder. The sound pressure level at a location 35 diameters away from the cylinder and at an angle of 90 deg with respect to the wake's downstream axis was found to have a peak value of approximately 110 db. Slightly smaller peak values were predicted at the 60 deg and 120 deg locations. A grid refinement study suggests that the dynamic model demands mesh refinement beyond that used here.

Spyropoulos, Evangelos T.↗
DOE OSTI2023

Tensor-based flow reconstruction from optimally located sensor measurements

Reconstructing high-resolution flow fields from sparse measurements is a major challenge in fluid dynamics. Existing methods often vectorize the flow by stacking different spatial directions on top of each other, hence confounding the information encoded in different dimensions. Here, we introduce a tensor-based sensor placement and flow reconstruction method which retains and exploits the inherent multidimensionality of the flow. We derive estimates for the flow reconstruction error, storage requirements and computational cost of our method. We show, with examples, that our tensor-based method is significantly more accurate than similar vectorized methods. Furthermore, the variance of the error is smaller when using our tensor-based method. While the computational cost of our method is comparable to similar vectorized methods, it reduces the storage cost by several orders of magnitude. The reduced storage cost becomes even more pronounced as the dimension of the flow increases. We demonstrate the efficacy of our method on three examples: a chaotic Kolmogorov flow, in situ and satellite measurements of the global sea surface temperature and three-dimensional unsteady simulated flow around a marine research vessel.

Mechanics↗
DOE OSTI2022

Flexible cryogenic air separation unit—An application for low-carbon fossil-fuel plants

In this study, the rapid integration of intermittent renewable sources into the electricity grid is driving the need for a flexible cryogenic air separation unit (ASU) coupled with a low-carbon fossil-fuel plant. However, the state-of-the-art ASU process is highly integrated and nonlinear, which can significantly restrict its ramping rate. In this work, we study the fundamental dynamics of a state-of-the-art double-column ASU, focusing on the dynamic characteristics of the highly integrated and nonlinear heat and mass transfers, and we propose basic control methods for achieving high process ramping rates. We found that the vapor–liquid countercurrent flows in the low-pressure column are critical to the cryogenic rectification of air, which governs the ramp rate of the ASU. This countercurrent flow structure is created through a complex heat integration process in the ASU. Here, this process is simplified as a countercurrent heat transfer to reduce the complexity of studying the ASU dynamics. To preserve this flow structure, the heat integration is maintained so that its heat duty follows the ASU load, using several basic controllers on the critical stream flowrates. A flow-driven dynamic ASU model is built in Aspen Plus Dynamics to capture the basic ramping dynamics. This model revealed a fundamental mismatch in the dynamics in the ASU column that causes a significant loss of O 2 product purity when ramping down the ASU and slightly increases the purity when ramping up. Based on these findings, we propose a basic control method for rapidly ramping down/up the ASU while maintaining O 2 product purity. Simulation results show the ASU basic dynamic process successfully ramps at a rate up to 10 %/min (40–100 % load) while maintaining the O 2 product purity at 95.2–95.6 mol%.

42 ENGINEERING↗
NASA NTRS2005

Benchmark Wall Heat Flux Data for a GO2/GH2 Single Element Combustor

Wall heat flux measurements in a 1.5 in. diameter circular cross-section rocket chamber for a uni-element shear coaxial injector element operating on gaseous oxygen (GOz)/gaseous hydrogen (GH,) propellants are presented. The wall heat flux measurements were made using arrays of Gardon type heat flux gauges and coaxial thermocouple instrumentation. Wall heat flux measurements were made for two cases. For the first case, GOZ/GHz oxidizer-rich (O/F=l65) and fuel-rich preburners (O/F=1.09) integrated with the main chamber were utilized to provide vitiated hot fuel and oxidizer to the study shear coaxial injector element. For the second case, the preburners were removed and ambient temperature gaseous oxygen/gaseous hydrogen propellants were supplied to the study injector. Experiments were conducted at four chamber pressures of 750, 600, 450 and 300psia for each case. The overall mixture ratio for the preburner case was 6.6, whereas for the ambient propellant case, the mixture ratio was 6.0. Total propellant flow was nominally 0.27-0.29 Ibm/s for the 750 psia case with flowrates scaled down linearly for lower chamber pressures. The axial heat flux profile results for both the preburner and ambient propellant cases show peak heat flux levels a t axial locations between 2.0 and 3.0 in. from the injector face. The maximum heat flux level was about two times greater for the preburner case. This is attributed to the higher injector fuel-to-oxidizer momentum flux ratio that promotes mixing and higher initial propellant temperature for the preburner case which results in a shorter reaction zone. The axial heat flux profiles were also scaled with respect to the chamber pressure to the power 0.8. The results at the four chamber pressures for both cases collapsed to a single profile indicating that at least to first approximation, the basic fluid dynamic structures in the flow field are pressure independent as long as the chamber/njector/nozzle geometry and injection velocities remain the same.

Marshall, William M.↗
DOE OSTI2024

Impact of breech geometry and propellant flow on the release of large pellets for the ITER disruption mitigation system

Studies have been performed on the release mechanism for large pellets using high pressure gas in a shattered pellet injector. Typically, pellets are dislodged from the cryogenic surface and accelerated down a barrel using high pressure gas delivered by a fast-acting propellant valve. The pellets impact an angled surface which shatters the pellet into many small fragments before entering the plasma. This technique was initially demonstrated on DIII-D (Commaux et al 2016 Nucl. Fusion 56 046007) and is now deployed on JET, KSTAR, ASDEX-Upgrade, and other tokamaks around the world in support of ITER's disruption mitigation system design and physics basis. The large hydrogen, 28.5 mm diameter, 2 length-to-diameter ratio, pellets foreseen for ITER SPI operation have low material strength and low heat of sublimation, which cause the pellets to be fragile and highly reactive to the impact of warm propellant gas. Due to the size of the pellets, significantly more propellant gas is required to dislodge and accelerate them. This creates a potentially significant propellant gas removal issue as 2–6 bar-L of gas is expected to be required for release and speed control. The research presented in this paper is an in-depth exploration of the parameters that are keys to reliable pellet release and speed control. Computational fluid dynamics (CFD) modeling of propellant flows through various breech designs was conducted to determine the force generated on the back surface of a pellet. These simulations assumed the use of the ORNL designed flyer plate valve. CFD modeling combined with experimental measurements provide adequate insight to determine a path to an optimal valve and breech design for ITER SPI pellet release and speed control while minimizing propellant gas usage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Direct numerical simulation of open-channel flow over a heterogeneous particle bed at low relative submergence

In this article, this study investigates turbulent open-channel flows over beds of irregularly arranged particles, using direct numerical simulations at a friction Reynolds number of Re T = 300⁠. Two distinct cases are examined: a polydisperse bed (P800) composed of multiple layers of randomly distributed spheres of varying sizes, and a monodisperse bed (M1015) formed by a random distribution of uniform sized spheres, with a bottommost single layer of varied-sized particles to introduce realistic randomness. Our investigation unveils a rich network of low- and high-speed streaks within the flow field, exhibiting distinctive behaviors in different bed configurations. The P800 case presents a poorly organized flow pattern induced by the varied particle sizes and arrangements, while the M1015 case shows a more regular flow pattern, marked by larger streaks. We also observe that total wall shear stress is substantially influenced by surface roughness-induced drag, extending beyond the effects documented in existing studies of open-channel flows. The present study reveals intricate secondary flow patterns over irregular particle beds. Large-scale circulations are discerned around particle crests in the P800 case and localized circulations with increased turbulence in the M1015 case. Furthermore, analysis of Reynolds stress tensor components indicates that roughness disrupts coherent turbulent eddies, consequently mitigating peak stress. We quantify correlations between drag force and local fluid velocity fluctuations. Notably, a larger deviation in drag is observed in the P800 case compared to M1015, accentuating the influence of particle size and distribution on fluid–particle interactions.

42 ENGINEERING↗
NASA NTRSDate not supplied

Cavitation Effects on the Structural Dynamics of Turbomachinery Components

The structural integrity of inducer and impeller blades in rocket engine turbo machinery must be evaluated in the face of complex excitation mechanisms including fluctuating pressures due to cavitation. Cavitation occurs when the local fluid pressure drops below the vapor pressure, causing the formation of vapor-filled bubbles. Cavitation can exist to various extents within the typical operating range of rocket engine turbo pumps. Despite recent progress towards reducing uncertainties in structural dynamic models of turbo machinery components, the extent to which pump cavitation affects the structural dynamic properties (i.e., natural frequencies, damping, and mode shapes) of inducer and impeller blades remains largely unknown. To study the structural dynamic effects of cavitating flows, experiments are conducted in a high-speed water tunnel. The test article is a low aspect ratio cantilevered plate, thus mimicing a single inducer blade that is unwrapped from its hub. The test article is oriented at zero angle of attack and has a triangular trip at its leading edge to induce cavitation. To change the extent of cavitation, the static pressure of the water tunnel is controlled. A high-speed camera is used to quantify the coverage of the cavitation sheet.

Structural Dynamics↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗
NASA NTRSDate not supplied

Cavitation Effects on the Structural Dynamics of Turbomachinery Components

The structural integrity of inducer and impeller blades in rocket engine turbomachinery must be evaluated in the face of complex excitation mechanisms including fluctuating pressures due to cavitation. Cavitation occurs when the local fluid pressure drops below the vapor pressure, causing the formation of vapor-filled bubbles. Cavitation can exist to various extents within the typical operating range of rocket engine turbopumps. Despite recent progress towards reducing uncertainties in structural dynamic models of turbomachinery components, the extent to which pump cavitation affects the structural dynamic properties (i.e., natural frequencies, damping, and mode shapes) of inducer and impeller blades remains largely unknown. To study the structural dynamic effects of cavitating flows, experiments are conducted in a high-speed water tunnel. The test article is a low aspect ratio cantilevered plate, thus mimicing a single inducer blade that is unwrapped from its hub. The test article is oriented at zero angle of attack and has a triangular trip at its leading edge to induce cavitation. To change the extent of cavitation, the static pressure of the water tunnel is controlled. A high-speed camera is used to quantify the coverage of the cavitation sheet.

Fluid-Structure Interaction↗