Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow matching”

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 55 records · Page 3

FUN3D and USM3D Analysis of the Propulsion Aerodynamic Workshop 2018 S-Duct Test Case

This work presents the results of Fun3D and USM3D analyses that were performed for the 4th AIAA Propulsion Aerodynamics Workshop (PAW). The PAW workshop is separated into three sections that focus on internal duct flows, nozzle flows and a special topic. This paper focuses on the internal duct flow section of PAW04 while an accompanying paper discusses the analyses performed for the nozzle portion. For the internal duct flow section, the PAW04 participants were provided with the two configurations consisting of an S-duct with and without aerodynamic interface plane (AIP) rake legs modeled. The participants were asked to perform a grid refinement study as well as a turbulence model study for the configuration with the rake legs. The analyses discussed here were performed on custom grids developed under the guidelines of the workshop. Additionally, the paper discusses the development and use of flow controllers for matching the desired flow characteristics. The results show that both solvers do well for predicting internal flow characteristics of the S-duct based on direct comparison with the experimental data. However, the CFD-to-CFD comparison proved to be more challenging due to the localized occurrence of supersonic flow near the rake legs when using the mass flow controller. A turbulence model study was performed to compare the two-equation SST model to the SA-QCR model. The results show that although the turbulence model does affect the solution, it makes a minimal impact on pressure recovery and inlet distortion intensity for this case. Suggestions for future workshops include gridding guidelines similar to those employed for the Drag Prediction Workshop series for the grid refinement study and a time accuracy study.

Carter, Melissa B.↗

Joule-Thomson Cryogenic Nitrogen Orifice Flows

The Joule-Thomson Effect involves the expansion of a fluid through a throttling device, often a simple orifice, resulting in a change in the fluid's temperature. It has been leveraged in the design and analysis of Thermodynamic Vent Systems and integrated into injector designs for cryogenic no-vent fill tank-to-tank transfer testing. There is a scarcity of comprehensive data available to model the phenomenon effectively using tools such as GFSSP for cryogenic nitrogen flows. This experiment aims to model the Joule-Thomson effect of liquid nitrogen through seven traditionally manufactured, sharp-edged orifices manufactured by O’Keefe and three additively manufactured (AM) orifices manufactured by Cumberland Additive. The traditional orifice sizes were 0.016, 0.02, 0.024, 0.032, 0.037, and 0.04 inch in diameter. The three additive orifices had diameters with target/nominal sizes equal to the 0.016, 0.024, and 0.04 inch O’Keefe orifices. However, the actual diameters of the AM orifices were found to be 0.014, 0.024, and 0.035 inches due to additive manufacturing restrictions. Water was flowed through the orifices to verify that their flow coefficients matched published values and to establish flow coefficients for the AM orifices. Liquid nitrogen was flowed through the orifices in a separate test apparatus. The liquid nitrogen started as a subcooled liquid upstream of the orifice and transitioned into a two-phase gas-liquid mixture after passing through each orifice. The upstream pressure was set to 29 psia, 34 psia, 39 psia, and/or 42 psia or 44 psia. The pressure, temperature, and flow rate were measured upstream and downstream of the orifice for three trials at each set point. The average volumetric flow rate and Joule-Thomson Coefficients were calculated from this data. The AM orifices and traditionally manufactured orifices displayed similar results. This data will refine the dataset GFSSP uses to model the J-T Effect more accurately.

Joule-Thomson↗

New results concerning the Galactic fountain

The time-dependent nonequilibrium thermal and ionization history is traced of hot optically-thin radiatively-cooling gas in a 1D planar steady-state flow model of the Galactic fountain. The calculation of such a flow which neglects photoionization is found to be unable to simultaneously reproduce the observed C IV, Si IV, and N V column densities through Galactic halo gas. It is shown that, when photoionization is taken into account, a fountain flow can match the observations of both UV absorption and emission lines for a range of ionizing fluxes which depends on the characteristic size of the cooling regions within the flow. What levels of the external photoionizing flux, such as that contributed by Galactic starlight, supernova remnants, and the metagalactic radiation background, can produce agreement with observational results are determined. The ionizing radiation emitted by the cooling gas itself is found to be sufficient to cause the flow to match the observed column densities and line emission.

Shapiro, Paul R.↗

The absorption and emission spectrum of radiative cooling galactic fountain gas

We have calculated the time-dependent, nonequilibrium thermal and ionization history of gas cooling radiatively from 10(exp 6) K in a one-dimensional, planar, steady-state flow model of the galactic fountain, including the effects of radiative transfer. Our previous optically thin calculations explored the effects of photoionization on such a flow and demonstrated that self-ionization was sufficient to cause the flow to match the observed galactic halo column densities of C 4, Si 4, and N 5 and UV emission from C 4 and O 3 in the constant density (isochoric) limit, which corresponded to cooling regions homogeneous on scales D less than or approximately equal to 1 kpc. Our new calculations which take full account of radiative transfer confirm the importance of self-ionization in enabling such a flow to match the data but allow a much larger range for cooling region sizes, i.e. D(sub 0) greater than or approximately equal to 15 pc. For an initial flow velocity v(sub 0) approximately equal to 100 km/s, comparable to the sound speed of a 10(exp 6) K gas, the initial density is found to be n(sub h,0) is approximately 2 x 10(exp -2) cm(exp -3), in reasonable agreement with other observation estimates, and D(sub 0) is approximately equal to 40 pc. We also compare predicted H(alpha) fluxes, UV line emission, and broadband x-ray fluxes with observed values. One dimensional numerical hydrodynamical calculations including the effects of radiative cooling are also presented.

Benjamin, Robert A.↗

A two phase Mach number description of the equilibrium flow of nitrogen in ducts

Some additional thermodynamic properties of the usual two-phase form which is linear in the moisture fraction are derived which are useful in the analysis of many kinds of duct flow. The method used is based on knowledge of the vapor pressure and Gibbs function as functions of temperature. With these, additional two-phase functions linear in moisture fraction are generated, which ultimately reveal that the squared ratio of mixture specific volume to mixture sound speed depends on liquid mass fraction and temperature in the same manner as do many weighted mean two-phase properties. This leads to a simple method of calculating two-phase Mach numbers for various duct flows. The matching of one- and two-phase flows at a saturated vapor point with discontinuous Mach number is also discussed.

Bursik, J. W.↗

Backward boundary layers in downward flame spread

The gas-phase combustion of a vertical fuel slab pyrolyzing at a finite rate under downward flame spread conditions is investigated theoretically. The backward boundary layer character of the flow is exploited, and a transformation that shifts the singularity associated with boundary-layer equations to a location, upstream of flame lift-off, where the surface temperature is nearly ambient, is employed. The boundary-layer flow is matched asymptotically for large Grashof number to a potential flow that allows the flow induced by the flame and along the upstream surface of the fuel bed to be calculated. Gas-phase temperature and velocity predictions are compared to experiment, and although analysis of the gas phase is emphasized, an approximate spread rate, based on a solid-phase energy balance, is calculated.

Vedha-Nayagam, M.↗

Development of the GPM Observatory Thermal Vacuum Test Model

A software-based thermal modeling process was documented for generating the thermal panel settings necessary to simulate worst-case on-orbit flight environments in an observatory-level thermal vacuum test setup. The method for creating such a thermal model involved four major steps: (1) determining the major thermal zones for test as indicated by the major dissipating components on the spacecraft, then mapping the major heat flows between these components; (2) finding the flight equivalent sink temperatures for these test thermal zones; (3) determining the thermal test ground support equipment (GSE) design and initial thermal panel settings based on the equivalent sink temperatures; and (4) adjusting the panel settings in the test model to match heat flows and temperatures with the flight model. The observatory test thermal model developed from this process allows quick predictions of the performance of the thermal vacuum test design. In this work, the method described above was applied to the Global Precipitation Measurement (GPM) core observatory spacecraft, a joint project between NASA and the Japanese Aerospace Exploration Agency (JAXA) which is currently being integrated at NASA Goddard Space Flight Center for launch in Early 2014. From preliminary results, the thermal test model generated from this process shows that the heat flows and temperatures match fairly well with the flight thermal model, indicating that the test model can simulate fairly accurately the conditions on-orbit. However, further analysis is needed to determine the best test configuration possible to validate the GPM thermal design before the start of environmental testing later this year. Also, while this analysis method has been applied solely to GPM, it should be emphasized that the same process can be applied to any mission to develop an effective test setup and panel settings which accurately simulate on-orbit thermal environments.

Yang, Kan↗

Theoretical versus experimental results for the rotordynamic coefficients of eccentric, smooth, gas annular seal annular gas seals

This viewgraph presentation presents the following results: (1) The analytical results overpredict the experimental results for the direct stiffness values and incorrectly predict increasing stiffness with decreasing pressure ratios. (2) Theory correctly predicts increasing cross-coupled stiffness, K(sub YX), with increasing eccentricity and inlet preswirl. (3) Direct damping, C(sub XX), underpredicts the experimental results, but the analytical results do correctly show that damping increases with increasing eccentricity. (4) The whirl frequency values predicted by theory are insensitive to changes in the static eccentricity ratio. Although these values match perfectly with the experimental results at 16,000 rpm, the results at the lower speed do not correspond. (5) Theoretical and experimental mass flow rates match at 5000 rpm, but at 16,000 rpm the theoretical results overpredict the experimental mass flow rates. (6) Theory correctly shows the linear pressure profiles and the associated entrance losses with the specified rotor positions.

Childs, Dara W.↗

Euler Flow Computations on Non-Matching Unstructured Meshes

Advanced fluid solvers to predict aerodynamic performance-coupled treatment of multiple fields are described. The interaction between the fluid and structural components in the bladed regions of the engine is investigated with respect to known blade failures caused by either flutter or forced vibrations. Methods are developed to describe aeroelastic phenomena for internal flows in turbomachinery by accounting for the increased geometric complexity, mutual interaction between adjacent structural components and presence of thermal and geometric loading. The computer code developed solves the full three dimensional aeroelastic problem of-stage. The results obtained show that flow computations can be performed on non-matching finite-volume unstructured meshes with second order spatial accuracy.

Gumaste, Udayan↗

Theory of discrete vortex noise

A solution is obtained for the noise induced by vortex shedding from the trailing edge of a semiinfinite flat plate. The method of solution for the potential flow is based on a procedure whereby an inner incompressible flow is matched to an outer acoustic field. Acoustic wave fronts are calculated and compared qualitatively with schlieren photographs of similar flows. The acoustic directivity pattern is obtained, and the extreme sensitivity of this parameter to the incompressible flow near the trailing edge of the plate is investigated.

Davis, S. S.↗

Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing

Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.

boundary-layer transition↗

Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing

Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.

computational modeling↗

Flowfield measurements in a separated and reattached flat plate turbulent boundary layer

The separation and reattachment of a large-scale, two-dimensional turbulent boundary layer at low subsonic speed on a flat plate has been studied experimentally. The separation bubble was 55 cm long and had a maximum bubble thickness, measured to the height of the mean dividing streamline, of 17 cm, which was twice the thickness of the inlet boundary layer. A combination of laser velocimetry, hot-wire anemometry, pneumatic probing techniques, and flow visualization were used as diagnostics. Principal findings were that an outer inviscid rotational flow was defined which essentially convected over the blockage associated with the inner, viscously dominated bubble recirculation region. A strong backflow region in which the flow moved upstream 100 percent of the time was measured near the test surface over the central 35 percent of the bubble. A laminar backflow boundary layer having pseudo-turbulent characteristics including a log-linear velocity profile was generated under the highly turbulent backflow. Velocity profile shapes in the reversed flow region matched a previously developed universal backflow profile at the upstream edge of the separation region but not in the steady backflow region downstream. A smoke flow visualization movie and hot-film measurements revealed low frequency nonperiodic flapping at reattachment. However, forward flow fraction data at reattachment and mean velocity profiles in the redeveloping boundary layer downstream of reattachment correlated with backward-facing step data when the axial dimension was scaled by the distance from the maximum bubble thickness to reattachment.

Patrick, William P.↗

A study of test section configuration for shock tube testing of transonic airfoils

Two methods are investigated for alleviating wall interference effects in a shock tube test section intended for testing two-dimensional transonic airfoils. The first method involves contouring the test section walls to match approximate streamlines in the flow. Contours are matched to each airfoil tested to produce results close to those obtained in a conventional wind tunnel. Data from a previous study and the present study for two different airfoils demonstrate that useful results are obtained in a shock tube using a test section with contoured walls. The second method involves use of a fixed-geometry slotted-wall test section to provide automatic flow compensation for various airfoils. The slotted-wall test section developed exhibited the desired performance characteristics in the approximate Mach number range 0.82 to 0.89, as evidenced by good agreement obtained between shock tube and wind tunnel results for several airfoil flows.

Cook, W. J.↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

P-107A hydrogen leak meter Final report

Flow improvement, mass flow calculations, matching cell factors, and temperature control of meter for hydrogen leakage measurements in liquid propulsion systems

MEASURING APPARATUS↗