Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow Conditions”

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 469 records · Page 26

Compressible flows with periodic vortical disturbances around lifting airfoils

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method that is presented fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. The numerical scheme is validated through extensive comparisons with known solutions to unsteady vortical flow problems. In general, it is seen that the agreement between the numerical and analytical results is very good for reduced frequencies ranging from 0 to 4, and for Mach numbers ranging from .1 to .8. Numerical results are also presented for a wide variety of flow configurations for the purpose of determining the effects of airfoil thickness, angle of attack, camber, and Mach number on the unsteady lift and moment of airfoils subjected to periodic vortical gusts. It is seen that each of these parameters can have a significant effect on the unsteady airfoil response to the incident disturbances, and that the effect depends strongly upon the reduced frequency and the dimensionality of the gust. For a one-dimensional (transverse) or two-dimensional (transverse and longitudinal) gust, the results indicate that airfoil thickness increases the unsteady lift and moment at the low reduced frequencies but decreases it at the high reduced frequencies. The results show that an increase in airfoil Mach number leads to a significant increase in the unsteady lift and moment for the low reduced frequencies, but a significant decrease for the high reduced frequencies.

Scott, James R.↗

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗

Computation of supersonic jet noise under imperfectly expanded conditions

The turbulent mixing noise of supersonic jet under imperfectly expanded conditions is calculated for convergent and convergent-divergent (CD) axisymmetric nozzle geometries. The noise prediction incorporates CFD solution of Navier-Stokes equations. The effect of grid resolution on shock structure computation is demonstrated. Mixing noise spectra predicted from fine and coarse grid solutions exhibit little sensitivity to the grid resolution. A proper grid resolution, however, results in a significant improvement in shock capturing capability and helps predictions agree favorably with experimental data. Good agreement between predicted noise spectra and data shows that the CFD-incorporated noise prediction scheme, which was demonstrated for shock-free conditions, works as well for shock-containing flow conditions.

Kim, Chan M.↗

Computation of supersonic jet noise under imperfectly expanded conditions

The turbulent mixing noise of supersonic jet under imperfectly expanded conditions is calculated for convergent and convergent-divergent (CD) axisymmetric nozzle geometries. The noise prediction incorporates CFD solution of Navier-Stokes equations. The effect of grid resolution on shock structure computation is demonstrated. Mixing noise spectra predicted from fine and coarse grid solutions exhibit little sensitivity to the grid resolution. A proper grid resolution, however, results in a significant improvement in shock capturing capability and helps predictions agree favorably with experimental data. Good agreement between predicted noise spectra and data shows that the CFD-incorporated noise prediction scheme, which was demonstrated for shock-free conditions, works as well for shock-containing flow conditions.

Kim, Chan M.↗

Hypersonic Boundary Layer Instability Over a Corner

A boundary-layer transition study over a compression corner was conducted under a hypersonic flow condition. Due to the discontinuities in boundary layer flow, the full Navier-Stokes equations were solved to simulate the development of disturbance in the boundary layer. A linear stability analysis and PSE method were used to get the initial disturbance for parallel and non-parallel flow respectively. A 2-D code was developed to solve the full Navier-stokes by using WENO(weighted essentially non-oscillating) scheme. The given numerical results show the evolution of the linear disturbance for the most amplified disturbance in supersonic and hypersonic flow over a compression ramp. The nonlinear computations also determined the minimal amplitudes necessary to cause transition at a designed location.

Balakumar, Ponnampalam↗

Aerosciences Implications for Uranus Aerocapture

Exploration of Uranus has been indicated as the highest priority new flagship planetary science missions for NASA in the 2022 Decadal Survey. Due to Uranus’ location as the second furthest planet from the Sun, interplanetary flight times of a traditional propulsive-capture orbiter can expect to take 13-15 years to reach Uranus while requiring up to two-thirds of the launch mass budget to be allocated to propellant alone. Aerocapture presents an alternative means to capture orbit around Uranus by using aerodynamic forces generated from flight through the planetary atmosphere to decelerate to orbital insertion while reducing time in interplanetary transit time and increasing on-orbit payload mass. To date, only one spacecraft (Voyager 2) has visited Uranus, and there have been no missions which have entered its atmosphere. While knowledge of planetary entry at rocky terrestrial planets with Air (N2/O2) and CO2 dominated atmospheres has greatly advanced over the past decades with a variety of missions to learn from, the H2/He composition and much larger atmosphere of ice giants such as Uranus poses new aerosciences implications which are not fully understood. These implications must be better understood not only for aerocapture, but for planetary entry at Uranus as a whole (such as a Uranus probe, as proposed by the 2022 Decadal Survey), as well as other ice giants like Neptune. To conduct a proper analysis of aerocapture methods, aerodynamic and aerothermal analysis of flight trajectories in Uranus atmosphere are required. Previous work has investigated aerocapture at Neptune (with similar atmosphere to Uranus) using a Mid-L/D configuration. Additional previous work considers the convective heating experienced in H2/He atmospheres. The current work looks to use heritage Low-L/D configurations, adding technical depth and maturation to analysis previously performed by Girija. Atmospheric flight at ice giants presents unique challenges, including but not limited to: fidelity of atmospheric models, large entry velocities due to location in the solar system, extended flight in rarefied flow conditions due to atmospheric scale, and unusual aerodynamic flow properties specific to Uranus and Ice Giant atmospheres.

Eli Shellabarger↗

Effects of Buoyancy on Lean Premixed V-Flames Part I: Laminar and Turblent Flame Structure

Laser schlieren and planar laser-induced fluorescence techniques have been used to investigate laminar and turbulent v-flames in +g, -g, and micro g under flow conditions that span the regimes of momentum domination (Ri < 0. 1) and buoyancy domination (Ri > 0.1). Overall flame features shown by schlieren indicate that buoyancy dominates the entire flow field for conditions close to Ri = 1. With decreasing Ri, buoyancy effects are observed only in the far-field regions. Analyses of the mean flame angles demonstrate that laminar and turbulent flames do not have similar responses to buoyancy. Difference in the laminar +g and -g flame angles decrease with Ri (i.e., increasing Re) and converge to the microgravity flame angle at the momentum limit (Ri - 0). This is consistent with the notion that the effects of buoyancy diminish with increasing flow momentum. The +g and -g turbulent flame angles, however, do not converge at Ri = 0. As shown by OH-PLIF images, the inconsistency in +g and -g turbulent flame angles is associated with the differences in flame wrinkles. Turbulent flame wrinkles evolve more slowly in +g than in -g. The difference in flame wrinkle structures, however, cannot be explained in terms of buoyancy effects on flame instability mechanisms. It seems to be associated with the field effects of buoyancy that stretches the turbulent flame brushes in +g and compresses the flame brush in -g. Flame wrinkling offers a mechanism through which the flame responds to the field effects of buoyancy despite increasing flow momentum. These observations point to the need to include both upstream and downstream contributions in theoretical analysis of flame turbulence interactions.

Cheng, Robert K.↗

Bubble Generation in a Continuous Liquid Flow Under Reduced Gravity Conditions

The present work reports a study of bubble generation under reduced gravity conditions for both co-flow and cross-flow configurations. Experiments were performed aboard the DC-9 Reduced Gravity Aircraft at NASA Glenn Research Center, using an air-water system. Three different flow tube diameters were used: 1.27, 1.9, and 2.54 cm. Two different ratios of air injection nozzle to tube diameters were considered: 0.1 and 0.2. Gas and liquid volumetric flow rates were varied from 10 to 200 ml/s. It was experimentally observed that with increasing superficial liquid velocity, the bubbles generated decreased in size. The bubble diameter was shown to increase with increasing air injection nozzle diameters. As the tube diameter was increased, the size of the detached bubbles increased. Likewise, as the superficial liquid velocity was increased, the frequency of bubble formation increased and thus the time to detach forming bubbles decreased. Independent of the flow configuration (for either single nozzle or multiple nozzle gas injection), void fraction and hence flow regime transition can be controlled in a somewhat precise manner by solely varying the gas and liquid volumetric flow rates. On the other hand, it is observed that uniformity of bubble size can be controlled more accurately by using single nozzle gas injection than by using multiple port injection, since this latter system gives rise to unpredictable coalescence of adjacent bubbles. A theoretical model, based on an overall force balance, is employed to study single bubble generation in the dynamic and bubbly flow regime. Under conditions of reduced gravity, the gas momentum flux enhances bubble detachment; however, the surface tension forces at the nozzle tip inhibits bubble detachment. Liquid drag and inertia can act either as attaching or detaching force, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with performed experiments. However, at higher superficial,liquid velocities, the bubble neck length begins to significantly deviate from the value of the air injection nozzle diameter and thus the theory no longer predicts the experiment behavior. Effects of fluid properties, injection geometry and flow conditions on generated bubble size are investigated using the theoretical model. It is shown that bubble diameter is larger in a reduced gravity environment than in a normal gravity environment at similar flow condition and flow geometry.

Pais, Salvatore Cezar↗

Downstream boundary conditions for viscous flow problems

The problem of the specification of artificial outflow conditions in flow problems is studied. It is shown that for transport type equations incorrect outflow conditions will adversely affect the solution only in a small region near the outflow boundary, while for elliptic equations, e.g. those governing the streamfunction or pressure, a correct boundary specification is essential. In addition, integral outflow boundary conditions for fluid dynamical problems are considered. It is shown that such conditions are well posed, and their effect on the solutions of the Navier-Stokes equations is also considered.

Fix, G.↗

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↗

Laser velocimeter measurements of the flow field generated by a forward-swept propfan during flutter

Results are presented from an investigation to measure the flow field generated by a forward-swept propfan operating in flutter at a low forward velocity. For comparison to the flutter condition, flow field data are also presented for a slightly reduced rotational speed just below flutter. The forward-swept propfan was tested as the front rotor of a counterrotating pusher propeller. A laser Doppler velocimeter (LDV) was used to measure the velocity field in planes normal to the model centerline downstream of the rotor and in planes of constant radius within the blade passages at each operating condition. A comparison of the data taken at the two different operating conditions indicated that the mean, time-averaged flow about the blades did not change drastically as the propfan rotational speed was increased from the stable operating point to the flutter condition. No regions of flow separation could be identified in the data plots of the mean intrablade flow field. The data also indicate that the relative flow about the blades remained subsonic during flutter operation. The blades were found to have a higher than expected tip loading at both operating conditions. This is thought to have been caused by the outer blade sections twisting under load to higher than expected effective blade angles. This high tip loading resulted in strong vortices and a very nonuniform flow downstream of the tips of the forward-swept blades. This high tip loading may also have caused the blade flutter.

Podboy, Gary G.↗

A 1D–0D–3D coupled model for simulating blood flow and transport processes in breast tissue

In this work, we present mixed dimensional models for simulating blood flow and transport processes in breast tissue and the vascular tree supplying it. These processes are considered, to start from the aortic inlet to the capillaries and tissue of the breast. Large variations in biophysical properties and flow conditions exist in this system necessitating the use of different flow models for different geometries and flow regimes. In total, we consider four different model types. First, a system of 1D nonlinear hyperbolic partial differential equations (PDEs) is considered to simulate blood flow in larger arteries with highly elastic vessel walls. Second, we assign 1D linearized hyperbolic PDEs to model the smaller arteries with stiffer vessel walls. The third model type consists of ODE systems (0D models). It is used to model the arterioles and peripheral circulation. Finally, homogenized 3D porous media models are considered to simulate flow and transport in capillaries and tissue within the breast volume. Sink terms are used to account for the influence of the venous and lymphatic systems. Combining the four model types, we obtain two different 1D–0D–3D coupled models for simulating blood flow and transport processes: The first model results in a fully coupled 1D–0D–3D model covering the complete path from the aorta to the breast combining a generic arterial network with a patient specific breast network and geometry. The second model is a reduced one based on the separation of the generic and patient specific parts. The information from a calibrated fully coupled model is used as inflow condition for the patient specific sub-model allowing a significant computational cost reduction. Several numerical experiments are conducted to calibrate the generic model parameters and to demonstrate realistic flow simulations compared to existing data on blood flow in the human breast and vascular system. Moreover, we use two different breast vasculature and tissue data sets to illustrate the robustness of our reduced sub-model approach.

multiscale blood flow models↗

Liquid-phase turbulence measurements in air-water two-phase flows using particle image velocimetry

Liquid-phase turbulence plays a vital role in determining various gas-liquid two-phase flow parameters, such as void fraction distribution, bubble morphology, bubble-bubble interactions, and interfacial area concentration. In this study, a two-phase flow database including both the gas- and liquid-phases measurements was developed focusing on three bubbly flow conditions in an air-water two-phase flow loop with a vertical one-inch diameter circular pipe test section. A particle image velocimetry (PIV) system integrating an optical phase separation method, i.e., the planar laser-induced fluorescence (PLIF) technique using fluorescent particles and optical filtration, was applied to measure the liquid-phase turbulence information, including the time-averaged velocity, Reynolds stress, and turbulent kinetic energy for the liquid phase. The PIV measurements were taken at three ports along the test section at 14.5, 51.5, and 88.5 pipe inner diameters downstream of a bubble injector. In addition, a double-sensor conductivity probe was used to measure radial distributions of the local time-averaged void fraction and gas velocity. The measured liquid-phase turbulence was used to benchmark Sato's turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. Lastly, in the two bubbly flows tested that have low void fractions being less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction.

42 ENGINEERING↗

Pressure Drop with Moving Contact Lines and Dynamic Contact Angles in a Hydrophobic Round Minichannel: Visualization via Synchrotron X-ray Imaging and Verification of Experimental Correlations

In hydrophobic mini- and microchannels, slug flow with moving contact lines is typically generated under various two-phase flow conditions. There is a significant pressure drop in this flow pattern with moving contact lines, which is closely related to the dynamic contact angles. Researchers have investigated dynamic contact angles experimentally for decades, but due to the limitations of visualization techniques, these experiments have typically been conducted in low Weber number regions (We < 10 –3 ). In this study, we clearly visualized the dynamic contact angles of a liquid slug in high Weber number regions (10 –3 < We <1) via synchrotron X-ray imaging with high temporal (~1000 fps) and spatial (~ 2 μm/pixel) resolutions. We precisely measured the pressure drop with moving contact lines in a hydrophobic minichannel (inner diameter = 1.018 mm). On the basis of our experimental data, we verified previous correlations for dynamic contact angles and explored the relationship between pressure drop with moving contact lines and dynamic contact angles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding spectral dependance of laser-induced damage precursors in dielectric materials (Abbreviated report_23-ERD-006)

High peak and average power laser systems are typically limited by the handling fluence of the optical components. In particular, the multilayer dielectric coatings are known to be much lower operational fluence than the more ideal bulk materials. In this project we proposed and succeeded in probing different established classes of damage prone precursors at different wavelengths to understand their fundamental laser damage response as a function of wavelength. This study helped to shed light on the fundamental physics of the non-linear precursors that govern laser damage phenomena for ns-regime pulsed laser damage. In this study, we utilized the onsite coating capabilities (VPL, IBS coating lab) to purposefully generate laser damage-prone precursors in hafnia-based coatings (both single and multi-layer coatings). Specifically, we engineered coatings with craze lines initiated by nodules, generated coatings with our xenon-based coating process to suppress nanobubble formation and generated hafnia coatings under controlled oxygen flow conditions to study hafnia sub-oxides and oxygen flow dependance. In these studies we found that craze lines are rife with precursors that are sensitive to ultra-violet light but not to infrared light; we found that the removal of nanobubbles helps with all wavelengths tested, but is most impactful for ultra-violet light; we also found that ultra-violet laser damage performance of hafnia is closely matched to oxygen flow rate, while the infrared performance may be slightly better at lower flow rates. During this LDRD we also successfully stood up a new laser damage testing capability, namely a wavelength agile damage test station to study the spectral response of known laser-induced damage precursors. This is a unique and important capability for Lawrence Livermore National Lab, allowing us to understand the spectral response of materials under high intensity irradiation and damage. This is a vital tool to understand non-linear optical response at wavelengths that we have previously been unable to test at.

36 MATERIALS SCIENCE↗

Topics in structural dynamics: Nonlinear unsteady transonic flows and Monte Carlo methods in acoustics

The results are reported of two unrelated studies. The first was an investigation of the formulation of the equations for non-uniform unsteady flows, by perturbation of an irrotational flow to obtain the linear Green's equation. The resulting integral equation was found to contain a kernel which could be expressed as the solution of the adjoint flow equation, a linear equation for small perturbations, but with non-constant coefficients determined by the steady flow conditions. It is believed that the non-uniform flow effects may prove important in transonic flutter, and that in such cases, the use of doublet type solutions of the wave equation would then prove to be erroneous. The second task covered an initial investigation into the use of the Monte Carlo method for solution of acoustical field problems. Computed results are given for a rectangular room problem, and for a problem involving a circular duct with a source located at the closed end.

Haviland, J. K.↗