Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

Three dimensional flow simulation with application to aeroelastic analysis

The three-dimensional flowfield about realistic launch vehicle configurations is simulated using the Reynolds-averaged Navier-Stokes equations. Turbulent mixing is accounted for by means of the two-layer Baldwin and Lomax (1978) algebraic eddy viscosity model. The Beam and Warming (1976) implicit approximate factorization algorithm is used for the solution of the finite difference equations. Applications include the study of the flowfield about a hemisphere-cylinder configuration both in the subsonic and supersonic flight regimes, and about two hammerhead payload configurations at transonic speeds. A method is also described which permits incorporating this flow solver into a complete algorithm to perform time-domain aeroelastic stability analyses. The vehicle is modeled as a free-free beam and modal superposition techniques are used for the structural-dynamic formulation. Aeroelastic analyses were performed for the hammerhead configurations.

Azevedo, Joao Luiz F.↗

SRM Internal Flow Tests and Computational Fluid Dynamic Analysis: RSRM Full Scale Motor Analyses - Volume 2

The purpose of the RSRM Nozzle Slag Ejection Precursor Test is to investigate the effect that slag ejection from the RSRM nozzle has on the chamber pressure and trust of the SRB's. In past firings of the Reusable Solid Rocket Motor (RSRM) both static test and flight motors have shown small pressure perturbations occurring primarily between 65 and 80 seconds. A joint NASA/Thiokol team investigation concluded that the cause of the pressure perturbations was the periodic ingestion and ejection of molten aluminum oxide slag from the cavity around the submerged nozzle nose which tends to trap and collect individual aluminum oxide droplets from the approach flow. The conclusions of the team were supported by numerous data and observations from special tests including high speed photographic films, real time radiography, plume calorimeters, accelerometers, strain gauges, nozzle TVC system force gauges, and motor pressure and thrust data. A simplistic slag ballistics model was formulated to relate a given pressure perturbation to a required slag quantity. Also, a cold flow model using air and water was developed to provide data on the relationship between the slag flow rate and the chamber pressure increase. Both the motor and the cold flow model exhibited low frequency oscillations in conjunction with periods of slag ejection. Motor and model frequencies were related to scaling parameters. The data indicate that there is a periodicity to the slag entrainment and ejection phenomena which is possibly related to organized oscillations from instabilities in the dividing streamline shear layer which impinges on the underneath surface of the nozzle.

Source record↗

Coherent wave-zonal mean flow interactions in the troposphere

A comprehensive analysis of wave-zonal mean flow variations in the troposphere is carried out using results of observations on the spatial structure and temporal evolution of synoptic time scale variations of the tropospheric zonal mean flow. Seven years of daily operational analysis are examined along with data from the NCAR GCM. Coherent wave-zonal mean flow interactions were revealed via cross-correlation analyses, including fluctuations in zonal mean temperature, three-dimensional winds, and quadratic wave quantities. It was found that tropical zonal wind tendency fluctuations are coherent with those anticipated from quasi-geostrophic wave driving, while tropical temperature tendencies are incoherent with analyzed zonal mean vertical velocities.

Randel, William J.↗

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗

Describing the Mechanism of Instability Suppression Using a Central Pilot Flame With Coupled Experiments and Simulations

Pilot flames are commonly used to extend combustor operability limits and suppress combustion oscillations in low-emissions gas turbines. Combustion oscillations, a coupling between heat release rate oscillations and combustor acoustics, can arise at the operability limits of low-emissions combustors where the flame is more susceptible to perturbations. While the use of pilot flames is common in land-based gas turbine combustors, the mechanism by which they suppress instability is still unclear. In this study, we consider the impact of a central jet pilot on the stability of a swirl-stabilized flame in a variable-length, single-nozzle combustor. Previously, the pilot flame was found to suppress the instability for a range of equivalence ratios and combustor lengths. We hypothesize that combustion oscillation suppression by the pilot occurs because the pilot provides hot gases to the vortex breakdown region of the flow that recirculate and improve the static, and hence dynamic, stability of the main flame. This hypothesis is based on a series of experimental results that show that pilot efficacy is a strong function of pilot equivalence ratio but not pilot flow rate, which would indicate that the temperature of the pilot products as well as the combustion intensity of the pilot flame play more of a role in oscillation stabilization than the length of the pilot flame relative to the main flame. Further, the pilot-flame efficacy increases with pilot-flame equivalence ratio until it matches the main-flame equivalence ratio; at pilot equivalence ratios greater than the main equivalence ratio, the pilot-flame efficacy does not change significantly with pilot equivalence ratio. To understand these results, we use large-eddy simulation (LES) to provide a detailed analysis of the flow in the region of the pilot flame and the transport of radical species in the region between the main flame and pilot flame. Furthermore, the simulation, using a flamelet/progress variable-based chemistry tabulation approach and standard eddy viscosity/diffusivity turbulence closure models, provides detailed information that is inaccessible through experimental measurements.

42 ENGINEERING↗

In situ feature analysis for large-scale multiphase flow simulations

The study of multiphase flow is essential for designing chemical reactors such as fluidized bed reactors (FBR), as a detailed understanding of hydrodynamics is critical for optimizing reactor performance and stability. An FBR allows scientists to conduct different types of chemical reactions involving multiphase materials, especially interaction between gas and solids. During such complex chemical processes, the formation of void regions in the reactor, generally termed as bubbles, is an important phenomenon. The study of these bubbles has a deep implication in predicting the reactor’s overall efficiency. But physical experiments needed to understand bubble dynamics are costly and non-trivial due to the technical difficulties involved and harsh working conditions of the reactors. Therefore, to study such chemical processes and bubble dynamics, a state-of-the-art computational simulation MFIX-Exa is being developed. Despite the proven accuracy of MFIX-Exa in modeling bubbling phenomena, the large-scale output data prohibits the use of traditional post hoc analysis capabilities in both storage and I/O time. Herein, to address these issues and allow the application scientists to explore the bubble dynamics in an efficient and timely manner, we have developed an end-to-end analytics pipeline that enables in situ detection of bubbles, followed by a flexible post hoc visual exploration methodology of bubble dynamics. The proposed method enables interactive analysis of bubbles, along with quantification of several bubble characteristics, enabling experts to understand the bubble interactions in detail. Positive feedback from the experts has indicated the efficacy of the proposed approach for exploring bubble dynamics in very-large-scale multiphase flow simulations.

97 MATHEMATICS AND COMPUTING↗

Two-phase flow measurements and observations in a cooling panel of the reactor cavity cooling system

The analysis of the flow behavior within a cooling panel of a scaled water-cooled reactor cavity cooling system has been performed supported by high-resolution measurements of the void fraction at the risers’ outlet. The results have shown a stable, symmetric distribution of the flow through the risers during subcooled boiling conditions, characterized by small, slow bubbles with average void fraction lower than 0.3 for all risers. As the temperature increased to saturation, flow distribution appeared strongly asymmetrical and unstable, with peaks of void fraction as high as 0.9. Large, faster bubbles were observed during these peaks, immediately followed by flow of subcooled liquid water. Flow stagnation and inversion has also been observed. The high resolution void fraction results, complemented by the measurements of the coolant flow rate and temperatures, not only have further described the behavior of the system under hypothetical accident conditions, but represent a unique set of data to support the validation of system-level codes, and more advanced computational tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstruct the acoustic signals at each microphone in the array. The mode structure of the lowest order mode is shown to be significantly impacted by shear flow refraction effects, and higher order mode structures are also effected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest order mode, but discrepancies arise when higher order modes are the dominant component of the acoustic field.

Duct acoustics↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstructs the acoustic signal at each microphone in the array. The mode structure of the lowest-order mode is shown to be significantly impacted by shear flow refraction effects, and higher-order mode structures are also affected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest-order mode, but discrepancies arise when higher-order modes are the dominant component of the acoustic field.

Duct acoustics↗

TRANSEP: A program for high lift separated flow about airfoils

A method and program called TRANSEP is presented that can be used for the analysis of the flow about a low speed airfoil under high lift, massive separation conditions. Since the present program is a modification of the direct-inverse TRANDES code, it can also be used for the design and analysis of transonic airfoils, including the effects of weak viscous interaction. Interactions on program usage, program modifications to convert TRANDES to TRANSEP, and sample cases and results are given.

Carlson, L. A.↗

Modeling study of divertor particle flow pattern and in–out plasma density asymmetry due to drifts with SOLPS and BOUT++

A study of the effects of drifts on the particle flow pattern and in–out divertor plasma density asymmetry for L-mode and H-mode plasmas is carried out for EAST discharges by the edge plasma transport codes SOLPS and BOUT++. The simulation of L-mode plasmas is done by SOLPS while the simulation of H-mode plasmas is done by BOUT++. The toroidal magnetic field direction for the simulated discharge is artificially reversed in the codes to study the effects of different drift directions on the divertor particle flow pattern and the in–out asymmetry of divertor plasma density. The divertor particle flows induced by diamagnetic and E × B drifts are found to have similar directions in the divertor region for the same discharge. The directions of the flows induced by drifts would be reversed with the reverse of toroidal magnetic field direction. The diamagnetic drift seems to have no effect on the in–out asymmetry of divertor plasma density due to its divergence-free nature. However, the E × B drift could result in a pronounced asymmetry of plasma density between the inner and outer divertor targets. The density in–out asymmetry caused by E × B drift is reversed with the reverse of E × B drift flow direction. Detailed analysis shows that the radial component of the E × B drift flow is the main cause of density asymmetry. The results from the simulation of H-mode plasmas with BOUT++ are similar to those of the L-mode plasmas with SOLPS except that the drift effects seem to be slightly larger in the H-mode plasmas compared to the L-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

CFD Analysis in Advance of the NASA Juncture Flow Experiment

NASA through its Transformational Tools and Technologies Project (TTT) under the Advanced Air Vehicle Program, is supporting a substantial effort to investigate the formation and origin of separation bubbles found on wing-body juncture zones. The flow behavior in these regions is highly complex, difficult to measure experimentally, and challenging to model numerically. Multiple wing configurations were designed and evaluated using Computational Fluid Dynamics (CFD), and a series of wind tunnel risk reduction tests were performed to further down-select the candidates for the final experiment. This paper documents the CFD analysis done in conjunction with the 6 percent scale risk reduction experiment performed in NASA Langley's 14- by 22-Foot Subsonic Tunnel. The combined CFD and wind tunnel results ultimately helped the Juncture Flow committee select the wing configurations for the final experiment.

Flow Experiment↗

Conjugate Heat Transfer Analysis of the Supercritical CO2 Based Counter Flow Compact 3D Heat Exchangers

Compact heat exchangers using supercritical fluids such as CO2 are preferred due to their high heat transfer capacity and smaller footprint. Three-dimensional conjugate forced convection heat transfer analysis was performed on several shell-and-tube counter-flow microchannel heat exchangers. Numerical simulations were conducted to test effect of change in mass flow rate, hydraulic diameter and various cross sections on the heat transfer. Increasing mass flow rate improved heat transfer up to a maximum value and then decreased downstream with increasing turbulence. Maximum heat transfer was obtained for the micro channel with the smallest hydraulic diameter. Amongst the cross sections analyzed (circular, square, circular with radial ribs, and square with radial ribs), the most uniform distribution of temperature and maximum heat transfer were obtained for circular cross section with radial ribs. An optimally efficient operation of such a heat exchanger can be attained by considering these factors during multi-objective constrained optimization of geometric parameters and requirements for additive manufacturing of such compact heat exchangers.

Chitale, Janhavi↗

Conjugate Heat Transfer Analysis of the Supercritical CO2 Based Counter Flow Compact 3D Heat Exchangers

Compact heat exchangers using supercritical fluids such as CO2 are preferred due to their high heat transfer capacity and smaller footprint. Three-dimensional conjugate forced convection heat transfer analysis was performed on several shell-and-tube counter-flow microchannel heat exchangers. Numerical simulations were conducted to test effect of change in mass flow rate, hydraulic diameter and various cross sections on the heat transfer. Increasing mass flow rate improved heat transfer up to a maximum value and then decreased downstream with increasing turbulence. Maximum heat transfer was obtained for the micro channel with the smallest hydraulic diameter. Amongst the cross sections analyzed (circular, square, circular with radial ribs, and square with radial ribs), the most uniform distribution of temperature and maximum heat transfer were obtained for circular cross section with radial ribs. An optimally efficient operation of such a heat exchanger can be attained by considering these factors during multi-objective constrained optimization of geometric parameters and requirements for additive manufacturing of such compact heat exchangers.

Chitale, Janhavi↗

Unraveling the discontinuous plastic flow of a Co-Cr-Fe-Ni-Mo multiprincipal-element alloy at deep cryogenic temperatures

We report an analysis of the discontinuous plastic flow of a multiprincipal-element alloy, Co 17.5 Cr 12.5 Fe 55 Ni 10 Mo 5 (atomic percent, at. %), in the temperature range of 0.5–4.2 K showing serrated deformation curves. Using the analytical techniques, we studied the statistics of the stress drops associated with the unstable plastic flow. The analysis showed that the complexity and heterogeneity of a discontinuous plastic flow were reduced when the temperature was lowered. This behavior was associated with the effects of dynamic recovery and adiabatic heating on the dislocation-density evolution.

36 MATERIALS SCIENCE↗

Analyzing Flow Fields in Axial-Compressor Rotors and Stators

Computer program, BLADE3D, developed for analysis of inviscid threedimensional flow fields in single blade passage of axial-compressor rotor or stator. Applicable to arbitrary axial-compressor hub, tip, and blade geometries, including blades with part-span dampers. Also used for axial turbines but with lowered resolution of flows about blunt leading and trailing edges. Analysis method valid for subsonic, transonic and supersonic flows, including choked flows.

Thompkins, W. J., Jr.↗