Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conditional Normalizing 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 379 records · Page 21

Gravitational Effects on Flow Instability and Transition in Low Density Jets

Recent experiments have shown that low-density gas jets injected into a high-density ambient gas undergo an instability mode, leading to highly-periodic oscillations in the flow-field for certain conditions. The transition from laminar to turbulent flow in these jets is abrupt, without the gradual change in scales. Even the fine scale turbulent structure repeats itself with extreme regularity from cycle to cycle. Similar observations were obtained in buoyancy-dominated and momentum-dominated jets characterized by the Richardson numbers, Ri = [gD(rho(sub a)-rho(sub j))/rho(sub j)U(sub j)(exp 2) ] where g is the gravitational acceleration, D is the jet diameter, rho(sub a) and rho(sub a) are, respectively, the free-stream and jet densities, and U(sub j) is the mean jet exit velocity. At high Richardson numbers, the instability is presumably caused by buoyancy since the flow-oscillation frequency (f) or the Strouhal number, St = [fD/U(sub j)] scales with Ri. In momentum-dominated jets, however, the Strouhal number of the oscillating flow is relatively independent of the Ri. In this case, a local absolute instability is predicted in the potential core of low-density jets with S [= rho(sub j)/rho(sub a)] < 0.7, which agrees qualitatively with experiments. Although the instability in gas jets of high Richardson numbers is attributed to buoyancy, direct physical evidence has not been acquired in experiments. If the instability is indeed caused by buoyancy, the near-field flow structure of the jet will change significantly when the buoyancy is removed, for example, in the microgravity environment. Thus, quantitative data on the spatial and temporal evolutions of the instability, length and time scale of the oscillating mode and its effects on the mean flow and breakdown of the potential core are needed in normal and microgravity to delineate gravitational effects in buoyant jets. In momentum dominated low-density jets, the instability is speculated to originate in the potential core. However, experiments have not succeeded in identifying the direct physical cause of the instability. For example, the theory predicts an oscillating mode for S<0.62 in the limit of zero momentum thickness, which contradicts with the experimental findings of Kyle and Sreenivasan. The analyses of momentum-dominated jets neglect buoyancy effects because of the small Richardson number. Although this assumption is appropriate in the potential core, the gravitational effects are important in the annular region surrounding the jet, where the density and velocity gradients are large. This reasoning provides basis for the hypothesis that the instability in low Richardosn number jets studied by Kyle and Sreenivasan and Monkewitz et al. is caused by buoyancy. The striking similarity in characteristics of the instability and virtually the identical conclusions reached by Subbarao and Cantwell in buoyant (Ri>0.5) helium jets on one hand and by Kyle and Sreenivasan in momentum-dominated (Ri<1x10(exp -3)) helium jets on the other support this hypothesis. However, quantitative experiments in normal and microgravity are necessary to obtain direct physical evidence of buoyancy effects on the flow instability and structure of momentum-dominated low-density jets. The primary objective of this new research project is to quantify how buoyancy affects the flow instability and structure in the near field of low-density jets. The flow will be described by the spatial and temporal evolutions of the instability, length and time scales of the oscillating mode, and the mean and fluctuating concentration fields. To meet this objective, concentration measurements will be obtained across the whole field using quantitative Rainbow Schlieren Deflectometry, providing spatial resolution of 0.1mm and temporal resolution of 0.017s to 1ms. The experimental effort will be supplemented with linear stability analysis of low-density jets by considering buoyancy. The first objective of this research is to investigate the effects of gravity on the flow instability and structure of low-density jets. The flow instability in these jets has been attributed to buoyancy. By removing buoyancy in our experiments, we seek to obtain the direct physical evidence of the instability mechanism. In the absence of the instability, the flow structure will undergo a significant change. We seek to quantify these changes by mapping the flow field (in terms of the concentration profiles) of these jets at non-buoyant conditions. Such information is presently lacking in the existing literature. The second objective of this research is to determine if the instability in momentum-driven, low-density jets is caused by buoyancy. At these conditions, the buoyancy effects are commonly ignored because of the small Richardson based on global parameters. By eliminating buoyancy in our experiments, globally as well as locally, we seek to examine the possibility that the instability mechanism in self-excited, buoyant or momentum-driven jets is the same. To meet this objective, we would quantify the jet flow in normal and microgravity, while systematically decreasing the Richardson number from buoyancy-driven to momentum driven flow regime. The third objective of this research is to perform a linear stability analysis of low-density gas jets by including the gravitational effects. The flow oscillations in these jets are attributed to an absolute instability, whereby the disturbance grows exponentially at the site to ultimately contaminate the entire flow field. We seek to study the characteristics of both convective and absolute instabilities and demarcate the boundary between them.

Agrawal A. K.↗

A finite element formulation for supersonic flows around complex configurations

The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element are assumed to be constant and equal to its value at the centroid of the element. This yields a set of linear algebraic equations whose coefficients are given by source and doublet integrals over the surface elements. Closed form evaluations of the integrals are presented.

Morino, L.↗

A finite-element analysis for steady and oscillatory supersonic flows around complex configurations

The problem of small perturbation potential supersonic flow around complex configurations is considered. This problem requires the solution of an integral equation relating the values of the potential on the surface of the body to the values of the normal derivative, which is known from the small perturbation boundary conditions. The surface of the body is divided into small (hyperboloidal quadrilateral) surface elements, sigma sub i, which are described in terms of the Cartesian components of the four corner points. The values of the potential (and its normal derivative) within each element is assumed to be constant and equal to its value at the centroid of the element, and this yields a set of linear algebraic equations. The coefficients of the equation are given by source and doublet integrals over the surface elements, sigma sub i. The results obtained using the above formulation are compared with existing analytical and experimental results.

Morino, L.↗

Predicted aerodynamics for a proposed personnel launch vehicle

The Langley Aerothermodynamic Upwind Relaxation Algorithm has been modified to compute equilibrium flows over shapes at hypersonic velocities. As an initial test of the modified code, solutions were obtained for a paraboloid at Mach 6 wind-tunnel test conditions in CF4 gas and at Mach 16.85 flight conditions in equilibrium air that approximated the normal shock density ratio for the experimental test. Solutions for a lifting-body concept that is proposed for a Personnel Launch System were compared with test results at Mach 6 in CF4 over an angle-of-attack range of 15 to 35 deg and with results for nominal entry flight conditions at Mach 15.

Greene, Francis A.↗

Pattern Formation in Diffusion Flames Embedded in von Karman Swirling Flows

Pattern formation is observed in nature in many so-called excitable systems that can support wave propagation. It is well-known in the field of combustion that premixed flames can exhibit patterns through differential diffusion mechanism between heat and mass. However, in the case of diffusion flames where fuel and oxidizer are separated initially there have been only a few observations of pattern formation. It is generally perceived that since diffusion flames do not possess an inherent propagation speed they are static and do not form patterns. But in diffusion flames close to their extinction local quenching can occur and produce flame edges which can propagate along stoichiometric surfaces. Recently, we reported experimental observations of rotating spiral flame edges during near-limit combustion of a downward-facing polymethylmethacrylate disk spinning in quiescent air. These spiral flames, though short-lived, exhibited many similarities to patterns commonly found in quiescent excitable media including compound tip meandering motion. Flame disks that grow or shrink with time depending on the rotational speed and in-depth heat loss history of the fuel disk have also been reported. One of the limitations of studying flame patterns with solid fuels is that steady-state conditions cannot be achieved in air at normal atmospheric pressure for experimentally reasonable fuel thickness. As a means to reproduce the flame patterns observed earlier with solid fuels, but under steady-state conditions, we have designed and built a rotating, porous-disk burner through which gaseous fuels can be injected and burned as diffusion flames. The rotating porous disk generates a flow of air toward the disk by a viscous pumping action, generating what is called the von K rm n boundary layer which is of constant thickness over the entire burner disk. In this note we present a map of the various dynamic flame patterns observed during the combustion of methane in air as a function of fuel flow rate and the burner rotational speed.

Nayagam, Vedha↗

Effect of Pressure on Piloted Ignition Delay of PMMA

In order to reduce the risk of decompression sickness associated with spacewalks, NASA is considering designing the next generation of exploration vehicles and habitats with a different cabin environment than used previously. The proposed environment uses a total cabin pressure of 52.7 to 58.6 kPa with an oxygen concentration of 30 to 34% by volume and was chosen with material flammability in mind. Because materials may burn differently under these conditions and there is little information on how this new environment affects the flammability of the materials onboard, it is important to conduct material flammability experiments at the intended exploration atmosphere. One method to evaluate material flammability is by its ease of ignition. To this end, piloted ignition delay tests were conducted in the Forced Ignition and Spread Test (FIST) apparatus subject to this new environment. In these tests, polymethylmethacylate (PMMA) was exposed to a range of oxidizer flow velocities and externally applied heat fluxes. The ultimate goal is to determine the individual effect of pressure and the combined effect of pressure and oxygen concentration on the ignition delay. Tests were conducted for a baseline case of normal pressure and oxygen concentration, low pressure (58.6 kPa) with normal oxygen (21%). Future work will focus on low pressure with 32% oxygen concentration (space exploration atmosphere - SEA) conditions. It was found that reducing the pressure while keeping the oxygen concentration at 21% reduced the ignition time by 17% on average. It was also noted that the critical heat flux for ignition decreases in low-pressure conditions. Because tests conducted in standard atmospheric conditions will underpredict the flammability of materials intended for use on spacecraft, fire safety onboard at exploration atmospheres may be compromised.

McAllister, Sara↗

Supersonic jet mixing enhancement by vortex generators

Experiments were conducted to explore the effects of vortex generators, in the form of tabs projecting normally into the jet, on the mixing and the far-field noise characteristics of a jet. A converging-diverging nozzle with a design Mach number of 1.36 was used in the experiments. The flow regimes from subsonic to highly underexpanded supersonic conditions were studied. One, two, and four tabs were used and some of the findings of previous investigators were examined and confirmed. The tabs eliminated screech noise from moderately overexpanded cases to highly underexpanded cases. Detailed flow visualizations and measurements showed that two tabs bifurcated the jet at all Mach numbers. While the effect of two tabs was persistent and the jet remained bifurcated, the distortions produced by one and four tabs disappeared by a streamwise distance of approximately 16 jet diameters. Two and four tabs significantly increased the entrainment of ambient air into the jet.

Samimy, M.↗

Pulmonary blood flow distribution in sheep: effects of anesthesia, mechanical ventilation, and change in posture

BACKGROUND: Recent studies providing high-resolution images of pulmonary perfusion have questioned the classical zone model of pulmonary perfusion. Hence the present work was undertaken to provide detailed maps of regional pulmonary perfusion to examine the influence of anesthesia, mechanical ventilation, and posture. METHODS: Pulmonary perfusion was analyzed with intravenous fluorescent microspheres (15 microm) in six sheep studied in four conditions: prone and awake, prone with pentobarbital-anesthesia and breathing spontaneously, prone with anesthesia and mechanical ventilation, and supine with anesthesia and mechanical ventilation. Lungs were air dried at total lung capacity and sectioned into approximately 1,100 pieces (about 2 cm3) per animal. The pieces were weighed and assigned spatial coordinates. Fluorescence was read on a spectrophotometer, and signals were corrected for piece weight and normalized to mean flow. Pulmonary blood flow heterogeneity was assessed using the coefficient of variation of flow data. RESULTS: Pentobarbital anesthesia and mechanical ventilation did not influence perfusion heterogeneity, but heterogeneity increased when the animals were in the supine posture (P < 0.01). Gravitational flow gradients were absent in the prone position but present in the supine (P < 0.001 compared with zero). Pulmonary perfusion was distributed with a hilar-to-peripheral gradient in animals breathing spontaneously (P < 0.05). CONCLUSIONS: The influence of pentobarbital anesthesia and mechanical ventilation on pulmonary perfusion heterogeneity is small compared with the effect of changes in posture. Analysis of flow gradients indicate that gravity plays a small role in determining pulmonary blood flow distribution.

NASA Discipline Cardiopulmonary↗

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗

Nonlinear Görtler Vortices and Their Secondary Instability in a Hypersonic Boundary Layer

Nonlinear development of the Görtler instability over a concave surface gives rise to a highly distorted inflectional flow field in the boundary layer that leads to both wall-normal and spanwise gradients in the flow. Such nonlinear structures are susceptible to strong, high-frequency secondary instabilities that may lead to the onset of laminar-turbulent transition. The present numerical study uses direct numerical simulations and linear secondary instability theory to investigate finite amplitude Görtler vortices and their secondary instability characteristics, respectively, in the hypersonic flow over an axisymmetric cone with a concave aft body. To complement previous studies in the literature wherein the Görtler instability was usually studied for a flat plate and initiated at some upstream location by imposing an eigenfunction as the inflow condition or by blowing and suction at the wall, the present investigation is focused on fully realizable Görtler instability that is excited by an azimuthally periodic array of surface protuberances. Furthermore, while the previous work had mostly focused on the secondary instability of Görtler vortices with cross-plane velocity contours that resembled bell-shaped structures, the present results confirm that fully developed mushroom structures also exist in the hypersonic regime when the Görtler vortex amplitude is sufficiently large. Computations further reveal that the dominant modes of secondary instability in these mushroom-shaped structures correspond to an antisymmetic (i.e., sinuous) “stem” mode that concentrates within the strong, nearly wall-normal internal shear layers surrounding the stem regions underneath the caps of the mushroom structures. Additionally, there exist a multitude of other significantly unstable secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes of both symmetric and antisymmetric types. Analogous to the secondary instability of crossflow vortices in hypersonic flows, secondary instability modes originating from the Mack mode instability play an important role during the nonlinear breakdown process.

Laminar turbulent transition↗

Low-speed wind-tunnel tests of an advanced eight-bladed propeller

As part of a research program on advanced turboprop aircraft aerodynamics, a low-speed wind-tunnel investigation was conducted to document the basic performance and force and moment characteristics of an advanced eight-bladed propeller. The results show that in addition to the normal force and pitching moment produced by the propeller/nacelle combination at angle of attack, a significant side force and yawing moment are also produced. Furthermore, it is shown that for test conditions wherein compressibility effects can be ignored, accurate simulation of propeller performance and flow fields can be achieved by matching the nondimensional power loading of the model propeller to that of the full-scale propeller.

Coe, P. L., Jr.↗

Time dependent inflow-outflow boundary conditions for 2D acoustic systems

An analysis of the number and form of the required inflow-outflow boundary conditions for the full two-dimensional time-dependent nonlinear acoustic system in subsonic mean flow is performed. The explicit predictor-corrector method of MacCormack (1969) is used. The methodology is tested on both uniform and sheared mean flows with plane and nonplanar sources. Results show that the acoustic system requires three physical boundary conditions on the inflow and one on the outflow boundary. The most natural choice for the inflow boundary conditions is judged to be a specification of the vorticity, the normal acoustic impedance, and a pressure gradient-density gradient relationship normal to the boundary. Specification of the acoustic pressure at the outflow boundary along with these inflow boundary conditions is found to give consistent reliable results. A set of boundary conditions developed earlier, which were intended to be nonreflecting is tested using the current method and is shown to yield unstable results for nonplanar acoustic waves.

Watson, Willie R.↗

Axisymmetric Two-Dimensional Modeling of No Vent Filling of a Cryogenic Tank using Generalized Fluid System Simulation Program

Filling a tank with cryogenic fluid is more challenging than filling a tank with water or any other fluid that is in liquid state at atmospheric condition. Filling a tank with cryogenic fluid is a two-step process. First the tank and the transfer line mustbe chilled. Liquid cryogens start flowing into the tank only after the tank and transfer lines are chilled to the fluid saturation temperature. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor, caused by heat transfer from the warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. Filling a cryogenic tank in the absence of gravity is more challenging because in a non-stratified environment liquid propellantmay not settle at the tank bottom as it does on earth. There is a strong possibility that liquid propellant may exit through the vent valve, which is typically located at the top of the tank to vent propellant vapor. There are several methodsof filling atank in space: Charge-Hold-Vent method [1], Vented-Chill / No-Vent-Fill (VCNVF)[2],and No Vent Fill with TVS assisted injector [3].This paper addressesnumerical modelingofthe No Vent Fill (NVF) processwhere a TVS (Thermodynamic Vent System) augmented injector was used with the vent valve closed during the entire filling process. NVF tests were conducted in the CRYOgenic Orbital Testbed (CRYOTE) tankat NASA/Marshall Space Flight Centerin 2018.The Generalized Fluid System Simulation Program (GFSSP), a general-purpose flow network code [4] was used to simulate the chilland fill process.

Cryogenics↗

A finite-element method for lifting surfaces in steady incompressible subsonic flow

The problem of potential steady subsonic flow for lifting surfaces is considered. This problem requires the solution of an integral equation relating the value of the potential discontinuity on the lifting surface and its wake to the values of the normal derivative of the potential which are known from the boundary conditions. The lifting surface is divided into small (quadrilateral hyperboloidal) surface elements, which are described in terms of the Cartesian components of the four corner points. The values of the potential discontinuity and the normal derivative of the potential are assumed to be constant within each element and equal to their values at the centroids of the elements. This yields a set of linear algebraic equations. Numerical results are in good agreement with existing ones.

Morino, L.↗

Effect of slotted casing treatment on performance of a multistage compressor

A J85-13 engine was equipped with a compressor case that allowed changes to the case wall over the rotor tips of six of its eight stages. The engine was tested with four inlet configurations: uniform inlet flow, 180 degree circumferential distortion, hub radial distortion, and tip radial distortion. Slotted inserts were installed in the first three stages, and the compressor was mapped under similar conditions. Overall compressor performance obtained with tip treatment was inferior to the performance for the compressor's normal operating range. Pumping capacity with the slotted inserts was reduced. Overall compressor efficiency was reduced 1 to 2 percent with the slotted rings installed for 90 and 100 percent corrected engine speeds.

Moss, J. E., Jr.↗

Comparison of predicted and experimental external heat transfer around a film cooled cylinder in crossflow

Calculations were made of the film cooling provided by rows of holes around the circumference of a cylinder in crossflow and the results were compared to experimental data. The calculations and experimental data were for conditions that simulate most of those that are typical of air cooled turbine vane leading edges. Injection was from a single and multiple rows of holes located at different angular locations from the stagnation line. The holes in the rows were angled normal to the flow direction and at a 25 degree angle to the cylinder wall. The calculations and experimental data were for several constant values of blowing ratios for all rows and for different blowing ratios for each row, representing a simulation of a common coolant plenum supply to multiple rows of holes. The calculations were made using a finite difference boundary layer code, STAN5. Contrary to initial expectations that injection would trip the boundary layer flow into the turbulent regime, the results indicated that the high free stream acceleration apparently kept the flow laminar for holes in the first 45 degrees past stagnation. The trend in Stanton number reduction due to coolant injection was predicted with generally good agreement at the lower blowing rates, but for multile rows of holes, agreement was poor beyond the first row.

Stepka, F. S.↗

Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver

A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-e two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user specified region. The latter approach is applied for the solutions obtained using other one- and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.

Pandya, Mohagna J.↗

Implementation of Flow Tripping Capability in the USM3D Unstructured Flow Solver

A flow tripping capability is added to an established NASA tetrahedral unstructured parallel Navier-Stokes flow solver, USM3D. The capability is based on prescribing an appropriate profile of turbulence model variables to energize the boundary layer in a plane normal to a specified trip region on the body surface. We demonstrate this approach using the k-epsilon two-equation turbulence model of USM3D. Modification to the solution procedure primarily consists of developing a data structure to identify all unstructured tetrahedral grid cells located in the plane normal to a specified surface trip region and computing a function based on the mean flow solution to specify the modified profile of the turbulence model variables. We leverage this data structure and also show an adjunct approach that is based on enforcing a laminar flow condition on the otherwise fully turbulent flow solution in user-specified region. The latter approach is applied for the solutions obtained using other one-and two-equation turbulence models of USM3D. A key ingredient of the present capability is the use of a graphical user-interface tool PREDISC to define a trip region on the body surface in an existing grid. Verification of the present modifications is demonstrated on three cases, namely, a flat plate, the RAE2822 airfoil, and the DLR F6 wing-fuselage configuration.

Pandya, Mohagna J.↗