Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

[Low-Frequency Flow Oscillation

The results of the research conducted under this grant are presented in detail in three Master theses, by Heinrich, Balow, and Broeren. Additional analysis of the experimental data can be found in two AIAA Journal articles and two conference papers. Citations for all of the studies' publications can be found in the bibliography which is attached. The objective of Heinrich's study was to document the low-frequency flow oscillation on the LRN-1007 airfoil, which had been previously observed at low Reynolds number, to determine its origin, and explore the phenomenon at higher Reynolds number. Heinrich performed detailed flow visualization on the airfoil using surface fluorescent oil and laser-sheet off-body visualization. A large leading-edge separation bubble and trailing-edge separation was identified on the airfoil just prior to the onset of the unsteady stall flow oscillation. From the laser-sheet data, the unsteady flow appeared as a massive boundary-layer separation followed by flow reattachment. Hot-wire data were taken in the wake to identify the presence of the flow oscillation and the dominant frequency. The oscillation was found in the flow from a Reynolds number of 0.3 to 1.3 x 10 exp 6. The Strouhal number based on airfoil projected height was nominally 0.02 and increased slightly with increasing Reynolds number and significantly with increasing airfoil angle of attack. Balow focused his research on the leading-edge separation bubble which was hypothesized to be the origin of the low-frequency oscillation. Initially, experimental measurements in the bubble at the onset of the low-frequency oscillation were attempted to study the characteristics of the bubble and explain possible relationships to the shear-layer-flapping phenomena. Unfortunately, the bubble proved to be extremely sensitive to the probe interference and it drastically reduced the size of the bubble. These detailed measurements were then abandoned by Balow. However, this led to a series of tests where the leading-edge bubble and trailing-edge separation were altered and the affect on the flow-oscillation studied. Balow found that by tripping the airfoil boundary-layer with "zigzag" tape ahead of bubble separation, the bubble was effectively eliminated mid the oscillation suppressed. Wake survey drag measurements showed a drastic reduction in airfoil drag when the bubble and oscillation were eliminated. Using the "zigzag" tape, the trailing-edge separation was moved downstream approximately 5 percent chord. This was found to reduce the amplitude of the oscillation, particularly in the onset stage at low angle of attack (around 14 degrees). Through detailed analysis of the wake behind the airfoil during the unsteady flow oscillation, Balow provided a better understanding of the wake flowfield. Broeren studied the oscillating flowfield in detail at Reynolds number equal 3 x 10 exp 5 and an angle of attack of 15 degrees using laser Doppler velocimetry (LDV). Two-dimensional LDV data were acquired at 687 grid points above the model upper surface while hot-wire data were taken simultaneously in the wake. Using the hot-wire signal, the LDV data were phase averaged into 24 bins to represent a single ensemble average of one oscillation cycle. The velocity data showed a flowfield oscillation that could be divided into three flow regimes. In the first regime, the flow over the airfoil was completely separated initially, the flowfield reattached from the leading edge and the reattachment point moved downstream with increasing time or phase. Broeren referred to this as the reattachment regime. The bubble development regime followed, where a leading-edge separation bubble formed at the leading edge and grew with increasing time. During the initial part of this regime the trailing-edge separation continued to move downstream. However, during the last 30 degrees of phase the trailing-edge separation moved rapidly forward and appeared to merge with the leading-edge bubble. During the third regime, the separation regime, the flow was segmented from the airfoil leading edge and did not reattach to the airfoil surface. The reverse flow was seen to grow in vertical extent up from the model surface as the phase increased. Next reattachment began again at the leading edge signaling the start of the reattachment regime, and so the cycle continued. From Broeren's work, the details of the unsteady flowfield over the airfoil were seen for the first time. From this research a great deal has been learned about the low-frequency flow oscillation which naturally occurs on the LRN-1007 airfoil near stall. The oscillation was seen to persist at higher Reynolds number, the dependence of the Strouhal number on angle of attack and Reynolds number were discovered, the critical role played by the laminar bubble was shown and the entire upper surface flowfield during a flow oscillation cycle was measured and analyzed. What still eludes understanding is the scaling of the flow oscillation and why certain airfoils, such as the LRN, have a very strong low-frequency mode and other airfoils exhibit no organized low-frequency oscillation at all.

Bragg, Michael B.↗
NASA NTRS2006

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan Brian G.↗
NASA NTRS2007

Criterion for Identifying Vortices in High-Pressure Flows

A study of four previously published computational criteria for identifying vortices in high-pressure flows has led to the selection of one of them as the best. This development can be expected to contribute to understanding of high-pressure flows, which occur in diverse settings, including diesel, gas turbine, and rocket engines and the atmospheres of Jupiter and other large gaseous planets. Information on the atmospheres of gaseous planets consists mainly of visual and thermal images of the flows over the planets. Also, validation of recently proposed computational models of high-pressure flows entails comparison with measurements, which are mainly of visual nature. Heretofore, the interpretation of images of high-pressure flows to identify vortices has been based on experience with low-pressure flows. However, high-pressure flows have features distinct from those of low-pressure flows, particularly in regions of high pressure gradient magnitude caused by dynamic turbulent effects and by thermodynamic mixing of chemical species. Therefore, interpretations based on low-pressure behavior may lead to misidentification of vortices and other flow structures in high-pressure flows. The study reported here was performed in recognition of the need for one or more quantitative criteria for identifying coherent flow structures - especially vortices - from previously generated flow-field data, to complement or supersede the determination of flow structures by visual inspection of instantaneous fields or flow animations. The focus in the study was on correlating visible images of flow features with various quantities computed from flow-field data.

Bellan, Josette↗
NASA NTRS2014

Spectroscopic Measurement Techniques for Aerospace Flows

The conditions that characterize aerospace flows are so varied, that a single diagnostic technique is not sufficient for its measurement. Fluid dynamists use knowledge of similarity to help categorize and focus on different flow conditions. For example, the Reynolds number represents the ratio of inertial to viscous forces in a flow. When the velocity scales, length scales, and gas density are large and the magnitude of the molecular viscosity is low, the Reynolds number becomes large. This corresponds to large scale vehicles (e.g Airbus A380), fast moving objects (e.g. artillery projectiles), vehicles in dense fluids (e.g. submarine in water), or flows with low dynamic viscosity (e.g. skydiver in air). In each of these cases, the inertial forces dominate viscous forces, and unsteady turbulent fluctuations in the flow variables are observed. In contrast, flows with small length scales (e.g. dispersion of micro-particles in a solid rocket nozzle), slow moving objects (e.g. micro aerial vehicles), flows with low density gases (e.g. atmospheric re-entry), or fluids with a large magnitude of viscosity (e.g. engine coolant flow), all have low Reynolds numbers. In these cases, viscous forces become very important and often the flows can be steady and laminar. The Mach number, which is the ratio of the velocity to the speed of sound in the medium, also helps to differentiate types of flows. At very low Mach numbers, acoustic waves travel much faster than the object, and the flow can be assumed to be incompressible (e.g. Cessna 172 aircraft). As the object speed approaches the speed of sound, the gas density can become variable (e.g. flow over wing of Learjet 85). When the object speed is higher than the speed of sound (Ma > 1), the presences of shock waves and other gas dynamic features can become important to the vehicle performance (e.g. SR-71 Blackbird). In the hypersonic flow regime (Ma > 5), large changes in temperature begin to affect flow properties, causing real-gas effects to occur (e.g. X-43 Scramjet). At even higher Mach numbers, chemistry and nonequilibrium effects come into play (e.g. Startdust re-entry capsule), further complicating the measurement. These limits can be predicted by calculating the ratio of chemical and thermal relaxation time to the flow time scales. Other non-dimensional numbers can be used to further differentiate types of aerospace flows.

Danehy, Paul M.↗
NASA NTRS2017

Field Detection of Chemical Assimilation in A Basaltic Lava Flow

Lava channels are features seen throughout the inner Solar System, including on Earth, the Moon, and Mars. Flow emplacement is therefore a crucial process in the shaping of planetary surfaces. Many studies, including some completed by members of this team at the December 1974 lava flow, have investigated the dynamics of lava flow emplacement, both on Earth and on the Moon and how pre-flow terrain can impact final channel morphology, but far fewer have focused on how the compositional characteristics of the substrate over which a flow was em-placed influenced its final flow morphology. Within the length of one flow, it is common for flows to change in morphology, a quality linked to rheology (a function of multiple factors including viscosi-ty, temperature, composition, etc.). The relationship between rheology and temperature has been well-studied but less is known about the relationship between an older flow's chemistry and how the interaction between this flow and the new flow might affect lava rheology and therefore emplacement dynamics. Lava erosion. Through visual observations of active terrestrial flows, mechanical erosion by flowing lava has been well-documented. Lava erosion by which flow composition is altered as the active lava melts and assimilates the pre-flow terrain over which it moves is also hypothesized to affect channel formation. However, there is only one previous field study that geochemically documents the process in recent basaltic flow systems.

Young, K. E.↗
NASA NTRSDate not supplied

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗
NASA NTRS2006

Numerical Modeling of Flow Control in a Boundary-Layer-Ingesting Offset Inlet Diffuser at Transonic Mach Numbers

This paper will investigate the validation of the NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as a baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet experiment conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a free-stream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the fanface diameter. The numerical simulations with and without tunnel walls are performed, quantifying tunnel wall effects on the BLI inlet flow. A comparison is made between the numerical simulations and the BLI inlet experiment for the baseline and VG vane cases at various inlet mass flow rates. A comparison is also made to a BLI inlet jet configuration for varying actuator mass flow rates at a fixed inlet mass flow rate. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCP avg, very well within the designed operating range of the BLI inlet. A comparison of the average total pressure recovery showed that the simulations were able to predict the trends but had a negative 0.01 offset when compared to the experimental levels. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion levels for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a misalignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery offset that was 0.01 lower than the experiment as was seen in the baseline. Comparisons of the flow features for the jet cases revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.

Allan, Brian G.↗
NASA NTRS2010

Minnowbrook VI: 2009 Workshop on Flow Physics and Control for Internal and External Aerodynamics

Topics covered include: Flow Physics and control for Internal and External Aerodynamics (not in TOC...starts on pg13); Breaking CFD Bottlenecks in Gas-Turbine Flow-Path Design; Streamwise Vortices on the Convex Surfaces of Circular Cylinders and Turbomachinery Blading; DNS and Embedded DNS as Tools for Investigating Unsteady Heat Transfer Phenomena in Turbines; Cavitation, Flow Structure and Turbulence in the Tip Region of a Rotor Blade; Development and Application of Plasma Actuators for Active Control of High-Speed and High Reynolds Number Flows; Active Flow Control of Lifting Surface With Flap-Current Activities and Future Directions; Closed-Loop Control of Vortex Formation in Separated Flows; Global Instability on Laminar Separation Bubbles-Revisited; Very Large-Scale Motions in Smooth and Rough Wall Boundary Layers; Instability of a Supersonic Boundary-Layer With Localized Roughness; Active Control of Open Cavities; Amplitude Scaling of Active Separation Control; U.S. Air Force Research Laboratory's Need for Flow Physics and Control With Applications Involving Aero-Optics and Weapon Bay Cavities; Some Issues Related to Integrating Active Flow Control With Flight Control; Active Flow Control Strategies Using Surface Pressure Measurements; Reduction of Unsteady Forcing in a Vaned, Contra-Rotating Transonic Turbine Configuration; Active Flow Control Stator With Coanda Surface; Controlling Separation in Turbomachines; Flow Control on Low-Pressure Turbine Airfoils Using Vortex Generator Jets; Reduced Order Modeling Incompressible Flows; Study and Control of Flow Past Disk, and Circular and Rectangular Cylinders Aligned in the Flow; Periodic Forcing of a Turbulent Axisymmetric Wake; Control of Vortex Breakdown in Critical Swirl Regime Using Azimuthal Forcing; External and Turbomachinery Flow Control Working Group; Boundary Layers, Transitions and Separation; Efficiency Considerations in Low Pressure Turbines; Summary of Conference; and Final Plenary Session Transcript.

LaGraff, John E.↗
DOE OSTI2023

Finding middle ground: Flow regimes designed for salmon and energy value

In regulated rivers, shaping seasonal flows to recover species at risk depends on understanding when to expect conflicts with competing water users and when their interests are aligned. Multi-objective optimization can be used to reveal such conflicts and commonalities. When species are involved, multi-objective optimization is challenged by the need to simulate complex species responses to flow regimes. Previously, we addressed that challenge by developing a simplified salmon model (Quantus) that defines cohorts of salmon by the river section and time in which they were spawned. Salmon in these space-time cohorts are tracked from the time redds (nests) are constructed until the cohort exits the tributary en route to the ocean. In this study, we modeled seasonal patterns in energy value and developed a Pareto-optimal frontier of seasonal flow patterns to maximize in-river salmon survival and hydropower value. Candidate flow regimes were characterized by two pulse flows varying in magnitude, timing, and duration and constrained by a total annual flow near the historical median. Our analysis revealed times when economic and salmon objectives were aligned and times when they differed. Pulse flows that favored higher energy value were timed to meet demand during extreme temperatures. Both salmon and hydropower objectives produced optimal flow regimes with pulse flows in early summer, but only solutions favoring hydropower value included high flows in mid-winter. Solutions favoring higher age-0 salmon survival provided an extended pulse flow in late winter/early spring, which suggests that access to productive floodplain habitat allowed faster growth and earlier out-migration and reduced the need for higher temperature-moderating flows later in spring. Minimum flows were also higher among solutions favoring salmon over energy. The tools used to produce these results can help to design simplified seasonal flow regimes by revealing compromise solutions that satisfy both fish and energy producers and highlighting when potential conflicts are likely.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗
DOE OSTI2020

Annular Flow Simulation Supported by Iterative In-Memory Mesh Adaptation

Various flow regimes exist in a boiling water reactor (BWR) as the steam quality increases in the uprising coolant flow, from bubbly flow, slug/churn flow, to annular flow. The annular flow is characterized by the presence of a fast-moving gas core and the surrounding liquid film flowing on the conduit wall. Additionally, entrained droplets can be observed in the gas core with ingested bubbles in the liquid film. The dynamics occurring on the wavy interface between the liquid film and gas core plays a crucial role in affecting the heat transfer rate and pressure drop within the BWR core. However, a fundamental understanding of annular flow is still lacking, partly due to the difficulty in obtaining detailed local data in annular flow experiments. In the current study, a novel simulation framework is developed for the annular flow by coupling a computational fluid dynamics flow solver with state-of-the-art meshing software. The gas-liquid interface is tracked with the level set method. Based on the computed flow solutions, the computational mesh is dynamically adapted in memory to meet the local mesh resolution requirement. This iterative simulation-adaptation framework can ensure the fine mesh resolution across the interface, which not only helps mitigate the mass conservation degradation known to level set methods but also improves the representation of dramatic interface topological changes such as wave breaking and droplet entrainment. The present investigation will shed light onto the complex interfacial processes involved in annular flow and generate much needed simulation data for annular flow modeling.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Overcoming persistent challenges in putting environmental flow policy into practice: a systematic review and bibliometric analysis

The implementation of environmental flows (e-flows) aims to reduce the negative impacts of hydrological alteration on freshwater ecosystems. Despite the growing attention to the importance of e-flows since the 1970s, actual implementation has lagged. Therefore, we explore the limitations in e-flows implementation, their systemic reasons, and solutions. We conducted a systematic review and a bibliometric analysis to identify peer-reviewed articles published on the topic of e-flows implementation research in the last two decades, resulting in 68 research and review papers. Co-occurrence of terms, and geographic and temporal trends were analyzed to identify the gaps in environmental water management and propose recommendations to address limitations on e-flows implementation. We identify the underlying causes and potential solutions to such challenges in environmental water management. The limitations to e-flow implementation identified were categorized into 21 classes. The most recognized limitation was the competing priorities of human uses of water (n = 29). Many secondary limitations, generally co-occurring in co-causation, were identified as limiting factors, especially for implementing more nuanced and sophisticated e-flows. The lack of adequate hydrological data (n = 24) and ecological data (n = 28) were among the most mentioned, and ultimately lead to difficulties in starting or continuing monitoring/adaptive management (n = 28) efforts. The lack of resource/capacity (n = 21), experimentation (n = 19), regulatory enforcement (n = 17), and differing authorities involved (n = 18) were also recurrent problems, driven by the deficiencies in the relative importance given to e-flows when facing other human priorities. In order to provide a clearer path for successful e-flow implementation, system mapping can be used as a starting point and general-purpose resource for understanding the sociohydrological problems, interactions, and inherited complexity of river systems. Secondly, we recommend a system analysis approach to address competing demands, especially with the use of coupled water-energy modeling tools to support decision-making when hydropower generation is involved. Such approaches can better assess the complex interactions among the hydrologic, ecological, socioeconomic, and engineering dimensions of water resource systems and their effective management. Lastly, given the complexities in environmental water allocation, implementation requires both scientific rigor and proven utility. Consequently, and where possible, we recommend a move from simplistic flow allocations to a more holistic approach informed by hydroecological principles. To ease conflicts between competing water demands, water managers can realize more 'pop per drop' by supporting key components of a flow regime that include functional attributes and processes that enhance biogeochemical cycling, structural habitat formation, and ecosystem maintenance.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1985

Correlation and analysis of oil flow data for an air-breathing missile model

This paper will present the results of an oil flow investigation on an airbreathing missile model. This oil flow study examined the flow around the model, which can be configured with both axisymmetric and two-dimensional inlets. Flow visualization analyses were conducted for both types of geometries by examining the surface flow patterns made visible by the oil flows for Mach numbers of 2.5 and 3.95. The analysis has shown the extent of flow spillage around the inlet which has helped explain the force and moment data collected during previous testing of the model. The oil flow data has also been used to develop guidelines for modeling the location of the crossflow separation line along inlet fairings. Finally, the oil flow analysis has been used to identify unique features of the boattail flow. These boattail flow characteristics have been correlated with previous oil flow analysis of noncircular body models. This paper demonstrates the use of this type of oil flow analysis in developing missile flow field analysis and aerodynamic predictions ranging from impact angle methods through Navier-Stokes methods.

Stoy, S. L.↗
NASA NTRS1988

Microgravity two-phase fluid flow pattern modeling

When gas and liquid mixtures flow in a pipe, the distribution of the two phases may take many forms. A flow pattern, or flow regime, is the characteristic spatial distribution of the phases of flow in a pipe. Because heat transfer and pressure drop are dependent on the characteristic distribution of phases, it is necessary to describe flow patterns in an appropriate manner so that a hydrodynamic or heat transfer theory applicable to that can be chosen. A theoretical two phase flow regime transition map under a microgravity environment was developed on physical concepts. These transitions use four basic flow patterns: dispersed flow, slug flow, stratified flow, and annular flow. The forces considered are body force, surface tension force, inertial force, friction, and the force from eddy turbulent fluctuation. Three dimensionless parameters were developed. Because these transition boundaries were developed based on physical concepts, they should be applicable to flow regimes occurring in various design conditions. Because the flow pattern data from KC-135 experiments are insufficient to verify these theoretical transition lines completely, an adiabatic experiment for flow regime analysis is recommended.

Lee, Doojeong↗
NASA NTRS1990

Navier-Stokes and Euler solutions for lee-side flows over supersonic delta wings. A correlation with experiment

An Euler flow solver and a thin layer Navier-Stokes flow solver were used to numerically simulate the supersonic leeside flow fields over delta wings which were observed experimentally. Three delta wings with 75, 67.5, and 60 deg leading edge sweeps were computed over an angle-of-attack range of 4 to 20 deg at a Mach number 2.8. The Euler code and Navier-Stokes code predict equally well the primary flow structure where the flow is expected to be separated or attached at the leading edge based on the Stanbrook-Squire boundary. The Navier-Stokes code is capable of predicting both the primary and the secondary flow features for the parameter range investigated. For those flow conditions where the Euler code did not predict the correct type of primary flow structure, the Navier-Stokes code illustrated that the flow structure is sensitive to boundary layer model. In general, the laminar Navier-Stokes solutions agreed better with the experimental data, especially for the lower sweep delta wings. The computational results and a detailed re-examination of the experimental data resulted in a refinement of the flow classifications. This refinement in the flow classification results in the separation bubble with the shock flow type as the intermediate flow pattern between separated and attached flows.

Mcmillin, S. Naomi↗
NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗
NASA NTRS1993

Impeller flow field characterization with a laser two-focus velocimeter

Use of Computational Fluid Dynamics (CFD) codes, prevalent in the rocket engine turbomachinery industry, necessitates data of sufficient quality and quantity to benchmark computational codes. Existing data bases for typical rocket engine configurations, in particular impellers, are limited. In addition, traditional data acquisition methods have several limitations: typically transducer uncertainties are 0.5% of transducer full scale and traditional pressure probes are unable to provide flow characteristics in the circumferential (blade-to-blade) direction. Laser velocimetry circumvents these limitations by providing +0.5% uncertainty in flow velocity and +0.5% uncertainty in flow angle. The percent of uncertainty in flow velocity is based on the measured value, not full range capability. The laser electronics multiple partitioning capability allows data acquired between blades as the impeller rotates, to be analyzed separately, thus providing blade-to-blade flow characterization. Unlike some probes, the non-intrusive measurements made with the laser velocimeter does not disturb the flow. To this end,, and under Contract (NAS8-38864) to the National Aeronautics and Space Administration (NASA) at Marshall Space Flight Center (MSFC), an extensive test program was undertaken at Rocketdyne. Impellers from two different generic rocket engine pump configurations were examined. The impellers represent different spectrums of pump design: the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) impeller was designed in the 1 1970's the Consortium for CFD application in Propulsion Technology Pump Stage Technology Team (Pump Consortium) optimized impeller was designed with the aid of modern computing techniques. The tester configuration for each of the impellers consisted of an axial inlet, an inducer, a diffuser, and a crossover discharge. While the tested configurations were carefully chosen to be representative of generic rocket engine pumps, several features of both testers were intentionally atypical. A crossover discharge, downstream of the impeller, rather than a volute discharge was used to minimize asymmetric flow conditions that might be reflected in the impeller discharge flow data. Impeller shroud wear ring radial clearances were purposely close to minimize leakage flow, thus increasing confidence in using the inlet data as an input to CFD programs. The empirical study extensively examined the flow fields of the two impellers via performance of laser two-focus velocimeter surveys in an axial plane upstream of the impellers and in multiple radial planes downstream of the impellers. Both studies were performed at the impeller design flow coefficients. Inlet laser surveys that provide CFD code inlet boundary conditions were performed in one axial plane, with ten radial locations surveyed. Three wall static pressures, positioned circumferentially around the impeller inlet, were used to identify asymmetrical pressure distributions in the inlet survey plane. The impeller discharge flow characterization consisted of three radial planes for the SSME HPFTP impeller and two radial planes for the Pump Consortium optimized impeller. Housing wall static pressures were placed to correspond to the radial locations surveyed with the laser velocimeter. Between five and thirteen axial stations across the discharge channel width were examined in each radial plane during the extensive flow mapping. The largely successful empirical flow characterization of two different impellers resulted in a substantial contribution to the limited existing data base, and yielded accurate data for CFD code benchmarking.

Brozowski, L. A.↗
NASA NTRS1998

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗