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 127 records · Page 7

User's Manual for Aerofcn: a FORTRAN Program to Compute Aerodynamic Parameters

The computer program AeroFcn is discussed. AeroFcn is a utility program that computes the following aerodynamic parameters: geopotential altitude, Mach number, true velocity, dynamic pressure, calibrated airspeed, equivalent airspeed, impact pressure, total pressure, total temperature, Reynolds number, speed of sound, static density, static pressure, static temperature, coefficient of dynamic viscosity, kinematic viscosity, geometric altitude, and specific energy for a standard- or a modified standard-day atmosphere using compressible flow and normal shock relations. Any two parameters that define a unique flight condition are selected, and their values are entered interactively. The remaining parameters are computed, and the solutions are stored in an output file. Multiple cases can be run, and the multiple case solutions can be stored in another output file for plotting. Parameter units, the output format, and primary constants in the atmospheric and aerodynamic equations can also be changed.

Conley, Joseph L.↗

Laminar Dust Flames: A Program of Microgravity and Ground Based Studies at McGill

Fundamental knowledge of heterogeneous combustion mechanisms is required to improve utilization of solid fuels (e.g. coal), safe handling of combustible dusts in industry, and solid propulsion systems. The objective of the McGill University research program on dust combustion is to obtain a reliable set of data on basic combustion parameters for dust suspensions (i.e. laminar burning velocity, flame structure, quenching distance, flammability limits, etc.) over a range of particle sizes, dust concentrations, and types of fuel. This set of data then permits theoretical models to be validated and, when necessary, new models to be developed to describe the detailed reaction mechanisms and transport processes. Microgravity is essential to the generation of a uniform dust suspension of arbitrary particle size and concentration. When particles with a characteristic size on the order of tens of microns are suspended, they rapidly settle in a gravitational field. To maintain a particulate in suspension for time duration adequate to carry out combustion experiments invariably requires continuous convective flow in excess of the gravitational settling velocity (which is comparable with and can even exceed the dust laminar burning velocity). This makes the experiments turbulent in nature and thus renders it impossible to study laminar dust flames. Even for small particle sizes on the order of microns, a stable laminar dust flow can be maintained only for relatively low dust concentrations at normal gravity conditions. High dust loading leads to gravitational instability of the dust cloud and to the formation of recirculation cells in the dust suspension in a confined volume, or to the rapid sedimentation of the dense dust cloud, as a whole, in an unconfined volume. Many important solid fuels such as carbon and boron also have low laminar flame speeds (of the order of several centimeters per second). Convection that occurs in combustion products due to buoyancy disrupts the low speed dust flames and makes observation of such flames at normal gravity difficult.

Goroshin, Sam↗

Deformation, Stress Relaxation, and Crystallization of Lithium Silicate Glass Fibers Below the Glass Transition Temperature

The deformation and crystallization of Li(sub 2)O (center dot) 2SiO2 and Li(sub 2)O (center dot) 1.6SiO2 glass fibers subjected to a bending stress were measured as a function of time over the temperature range -50 to -150 C below the glass transition temperature (Tg). The glass fibers can be permanently deformed at temperatures about 100 C below T (sub)g, and they crystallize significantly at temperatures close to, but below T,, about 150 C lower than the onset temperature for crystallization for these glasses in the no-stress condition. The crystallization was found to occur only on the surface of the glass fibers with no detectable difference in the extent of crystallization in tensile and compressive stress regions. The relaxation mechanism for fiber deformation can be best described by a stretched exponential (Kohlrausch-Williams-Watt (KWW) approximation), rather than a single exponential model.The activation energy for stress relaxation, Es, for the glass fibers ranges between 175 and 195 kJ/mol, which is considerably smaller than the activation energy for viscous flow, E, (about 400 kJ/mol) near T, for these glasses at normal, stress-free condition. It is suspected that a viscosity relaxation mechanism could be responsible for permanent deformation and crystallization of the glass fibers below T,

Ray, Chandra S.↗

Real gas flow parameters for NASA Langley 22-inch Mach 20 helium tunnel

A computational procedure was developed which can be used to determine the flow properties in hypersonic helium wind tunnels in which real gas behavior is significant. In this procedure, a three-coefficient virial equation of state and the assumption of isentropic nozzle flow are employed to determine the tunnel reservoir, nozzle, throat, freestream, and post-normal shock conditions. This method was applied to a range of conditions which encompasses the operational capabilities of the LaRC 22-Inch Mach 20 Helium Tunnel. Results are presented graphically in the form of real gas correction factors which can be applied to perfect gas calculations. Important thermodynamic properties of helium are also plotted versus pressure and temperature. The computational scheme used to determine the real-helium flow parameters was incorporated into a FORTRAN code which is discussed.

Hollis, Brian R.↗

Test Outline for Flutter Analysis of Rectangular Panels in Rarefied Flow Conditions

Jet plume impingement forces acting on large flexible space structures may precipitate dynamically unstable behavior during space flights. Typical operating conditions in space involve rarefied gas flow regimes which are intrinsically distinct from continuum gas flow and are normally modeled using the kinetic theory of gas flow. Docking and undocking operations of the Space Shuttle with the Russian Mir space laboratory represent a scenario in which the stability boundaries of solar panels may be of interest. Extensive literature review of research work on the dynamic stability of rectangular panels in rarefied gas flow conditions indicated the lack of published reports dealing with this phenomenon. A recently completed preliminary study for NASA JSC dealing with the mathematical analysis of the stability of two-degree-of-freedom elastically supported rigid panels under the effect of rarefied gas flow was reviewed. A test plan outline is prepared for the purpose of conducting a series of experiments on four rectangular rigid test articles in a vacuum chamber under the effect of continuous and pulsating Nitrogen jet plumes. The purpose of the test plan is to gather enough data related to a number of key parameters to allow the validation of the two-degree-of-freedom mathematical model. The hardware required careful design to select a very lightweight material while satisfying rigidity and frequency requirements within the constraints of the test environment. The data to be obtained from the vacuum chamber tests can be compared with the predicted behavior of the theoretical two-degree-of-freedom model. Using the data obtained in this study, further research can identify the limitations of the mathematical model. In addition modifications to the mathematical model can be made, if warranted, to accurately predict the behavior of rigid panels under rarefied gas flow regimes.

Akl, Fred A.↗

Performance of Thermal Mass Flow Meters in a Variable Gravitational Environment

The performance of five thermal mass flow meters, MKS Instruments 179A and 258C, Unit Instruments UFM-8100, Sierra Instruments 830L, and Hastings Instruments HFM-200, were tested on the KC-135 Reduced Gravity Aircraft in orthogonal, coparallel, and counterparallel orientations relative to gravity. Data was taken throughout the parabolic trajectory where the g-level varied from 0.01 to 1.8 times normal gravity. Each meter was calibrated in normal gravity in the orthogonal position prior to flight followed by ground testing at seven different flow conditions to establish a baseline operation. During the tests, the actual flow rate was measured independently using choked-flow orifices. Gravitational acceleration and attitude had a unique effect on the performance of each meter. All meters operated within acceptable limits at all gravity levels in the calibrated orthogonal position. However, when operated in other orientations, the deviations from the reference flow became substantial for several of the flow meters. Data analysis indicated that the greatest source of error was the effect of orientation, followed by the gravity level. This work emphasized that when operating thermal flow meters in a variable gravity environment, it is critical to orient the meter in the same direction relative to gravity in which it was calibrated. Unfortunately, there was no test in normal gravity that could predict the performance of a meter in reduced gravity. When operating in reduced gravity, all meters indicated within 5 percent of the full scale reading at all flow conditions and orientations.

Brooker, John E.↗

Prediction of Shock Wave Structure in Weakly Ionized Gas Flow by Solving MGD Equation

This paper reports the recent research results of shockwave structure predictions using a new developed code. The modified Rankine-Hugoniot relations across a standing normal shock wave are discussed and adopted to obtain jump conditions. Coupling a electrostatic body force to the Burnett equations, the weakly ionized flow field across the shock wave was solved. Results indicated that the Modified Rankine-Hugoniot equations for shock wave are valid for a wide range of ionization fraction. However, this model breaks down with small free stream Mach number and with large ionization fraction. The jump conditions also depend on the value of free stream pressure, temperature and density. The computed shock wave structure with ionization provides results, which indicated that shock wave strength may be reduced by existence of weakly ionized gas.

Deng, Z. T.↗

Left ventricular outflow tract mean systolic acceleration as a surrogate for the slope of the left ventricular end-systolic pressure-volume relationship

OBJECTIVE: The goal of this study was to analyze left ventricular outflow tract systolic acceleration (LVOT(Acc)) during alterations in left ventricular (LV) contractility and LV filling. BACKGROUND: Most indexes described to quantify LV systolic function, such as LV ejection fraction and cardiac output, are dependent on loading conditions. METHODS: In 18 sheep (4 normal, 6 with aortic regurgitation, and 8 with old myocardial infarction), blood flow velocities through the LVOT were recorded using conventional pulsed Doppler. The LVOT(Acc) was calculated as the aortic peak velocity divided by the time to peak flow; LVOT(Acc) was compared with LV maximal elastance (E(m)) acquired by conductance catheter under different loading conditions, including volume and pressure overload during an acute coronary occlusion (n = 10). In addition, a clinically validated lumped-parameter numerical model of the cardiovascular system was used to support our findings. RESULTS: Left ventricular E(m) and LVOT(Acc) decreased during ischemia (1.67 +/- 0.67 mm Hg.ml(-1) before vs. 0.93 +/- 0.41 mm Hg.ml(-1) during acute coronary occlusion [p < 0.05] and 7.9 +/- 3.1 m.s(-2) before vs. 4.4 +/- 1.0 m.s(-2) during coronary occlusion [p < 0.05], respectively). Left ventricular outflow tract systolic acceleration showed a strong linear correlation with LV E(m) (y = 3.84x + 1.87, r = 0.85, p < 0.001). Similar findings were obtained with the numerical modeling, which demonstrated a strong correlation between predicted and actual LV E(m) (predicted = 0.98 [actual] -0.01, r = 0.86). By analysis of variance, there was no statistically significant difference in LVOT(Acc) under different loading conditions. CONCLUSIONS: For a variety of hemodynamic conditions, LVOT(Acc) was linearly related to the LV contractility index LV E(m) and was independent of loading conditions. These findings were consistent with numerical modeling. Thus, this Doppler index may serve as a good noninvasive index of LV contractility.

Non-NASA Center↗

Stability aspects of diverging subsonic flow

The paper describes the computation of two-dimensional, subsonic, diverging internal flows and how they differ from the corresponding converging flows. Such diverging or decelerating flows occur in such obvious places as subsonic diffusers and inlets; however, such flows also occur in supersonic nozzles in the presence of a normal shock. The flow instability and its relation to the numerical method used, boundary conditions, and viscous effects are assessed both analytically and numerically. The inviscid flow is shown to be physically unstable and a poor representation of the true viscous flow.

Cline, M. C.↗

Mixing layer simulation by an improved three-dimensional vortex-in-cell algorithm

An extension is presented of the three-dimensional vortex-in-cell (VIC) method that allows for the simulation of incompressible turbulent flows such as plane mixing layers and wakes with exact boundary conditions in the normal direction. The vorticity is assumed to be confined between two parallel planes and the turbulence is assumed to be homogeneous in the two directions parallel to the planes. As an application, an analysis is presented of the evolution of the plane mixing layer subjected to random two- and three-dimensional initial disturbances.

Couet, B.↗

Sound radiation from annular ducts/nozzles using modal decomposition of in-duct acoustic power

Using a refined acoustic impulse technique, an experimental program was carried out to determine the proper procedure for measuring the acoustic power for both incident and reflected sound fields for annular ducts with and without nozzles at various flow conditions. Several radial and azimuthal pressure measurements were made, and using these complex pressure data, the modal content of the pressure field was determined. A mathematical approach for the modal decomposition of annular duct acoustic power was developed. The acoustic power due to each individual mode and also their total sum were derived. The acoustic characteristics of annular duct-nozzle systems were studied at various flow conditions. The results derived from the experiments include the termination reflection coefficient, the normalized transmission coefficients, the power transfer functions and the power imbalance.

Salikuddin, M.↗

An integral analysis of transonic normal shock wave/turbulent boundary layer interactions in internal flow

An approximate integral viscous-inviscid interaction method is presented for calculating the development of a turbulent boundary layer subjected to a normal shock wave induced adverse pressure gradient in an internal axisymmetric flow. The inflow conditions and the downstream pressure are provided for the computation. In the supersonic region of shock pressure rise, the Prandtl-Meyer function is used to couple the viscous and inviscid flows. An analytical model for the coupling process is postulated and appropriate equations are defined. Downstream of the sonic point, one-dimensional inviscid flow is assumed for coupling with the viscous flow. The turbulent boundary layer is calculated using Green's integral lag-entrainment method. Comparisons of the solutions with the experimental data are made for interactions which are unseparated, near separation and separated. For comparison purposes, solutions to the time-dependent, mass-averaged, Navier-Stokes equations incorporating a two-equation, Wilcox-Rubesin turbulence model are also shown. The computed results from the integral method show good agreement with experimental data for unseparated interactions and reasonable agreement with the trend of the viscous effects when the interaction becomes increasingly separated.

Om, D.↗

Surface-Tension Driven Convection

The present experiments are directed to examine surface-tension gradients. The surface-tension induced convection under reduced gravity it's comparison with normal gravity convection with all other identical conditions are studied. The intensities of such flows and their penetration depth below the free surface are also determined. In this way, the dominant role of a free surface in a reduced-gravity environment is indicated. The interface between two fluid phases can influence the motion of fluids when either the interface has finite curvature or when the interfacial tension varies from point to point. In both cases, forces appear in the interfacial region that can affect or generate fluid motions.

Ostrach, S.↗

Calculation of particle dynamics effects on laser velocimeter data

In order to analytically assess the effects of particle dynamics on laser velocimeter measurements, the particle response to the fluid in which it is entrained must be modeled. The equation of particle motion derived at Arnold Engineering and Development Center (AEDC), is used to model the motion of the particle. The equation is derived from the Navier-Stokes equation under the assumption of unsteady Stokes flow. It should be noted that the Stokes flow eliminates the need for a far field boundary condition, and hence eliminates any normal (lift) forces. To extend the equation to higher particle slip velocities, the Stokes drag is replaced with the particle drag law, modified to include heat transfer effects. The drag coefficient accounts for inertial, rarefaction, compressibility, shape, and heat transfer effects.

Nichols, R. H.↗

A New Formulation of Time Domain Boundary Integral Equation for Acoustic Wave Scattering in the Presence of a Uniform Mean Flow

It has been well-known that under the assumption of a constant uniform mean flow, the acoustic wave propagation equation can be formulated as a boundary integral equation, in both the time domain and the frequency domain. Compared with solving partial differential equations, numerical methods based on the boundary integral equation have the advantage of a reduced spatial dimension and, hence, requiring only a surface mesh. However, the constant uniform mean flow assumption, while convenient for formulating the integral equation, does not satisfy the solid wall boundary condition wherever the body surface is not aligned with the uniform mean flow. In this paper, we argue that the proper boundary condition for the acoustic wave should not have its normal velocity be zero everywhere on the solid surfaces, as has been applied in the literature. A careful study of the acoustic energy conservation equation is presented that shows such a boundary condition in fact leads to erroneous source or sink points on solid surfaces not aligned with the mean flow. A new solid wall boundary condition is proposed that conserves the acoustic energy and a new time domain boundary integral equation is derived. In addition to conserving the acoustic energy, another significant advantage of the new equation is that it is considerably simpler than previous formulations. In particular, tangential derivatives of the solution on the solid surfaces are no longer needed in the new formulation, which greatly simplifies numerical implementation. Furthermore, stabilization of the new integral equation by Burton-Miller type reformulation is presented. The stability of the new formulation is studied theoretically as well as numerically by an eigenvalue analysis. Numerical solutions are also presented that demonstrate the stability of the new formulation.

Hu, Fang↗

Experimental Investigation of Jet-Induced Mixing of a Large Liquid Hydrogen Storage Tank

Experiments have been conducted to investigate the effect of fluid mixing on the depressurization of a large liquid hydrogen storage tank. The test tank is approximately ellipsoidal, having a volume of 4.89 m(exp 3) and an average wall heat flux of 4.2 W/m(exp 2) due to external heat input. A mixer unit was installed near the bottom of the tank to generate an upward directed axial jet flow normal to the liquid-vapor interface. Mixing tests were initiated after achieving thermally stratified conditions in the tank either by the introduction of hydrogen gas into the tank or by self-pressurization due to ambient heat leak through the tank wall. The subcooled liquid jet directed towards the liquid-vapor interface by the mixer induced vapor condensation and caused a reduction in tank pressure. Tests were conducted at two jet submergence depths for jet Reynolds numbers from 80,000 to 495,000 and Richardson numbers from 0.014 to 0.52. Results show that the rate of tank pressure change is controlled by the competing effects of subcooled jet flow and the free convection boundary layer flow due to external tank wall heating. It is shown that existing correlations for mixing time and vapor condensation rate based on small scale tanks may not be applicable to large scale liquid hydrogen systems.

Lin, C. S.↗

Measurement and Analysis of Terminal Shock Oscillation and Buffet Forcing Functions on a Launch Vehicle Payload Fairing

The buffet loads on a launch vehicle payload shroud can be impacted by the unsteadiness associated with a terminal shock at high subsonic speeds. At these conditions, flow accelerates to supersonic speeds on the nose of the payload fairing and is terminated by a normal shock on the cylindrical section downstream of the nose cone/cylinder shoulder. The location of the terminal shock and associated separated boundary layer is affected by the freestream Mach number, Reynolds number, and the pitch/yaw of the launch vehicle. Furthermore, even when the freestream conditions and vehicle attitude are constant, this terminal shock oscillates on the surface of the vehicle. The time-varying surface pressure associated with the terminal shock results in unsteady aerodynamic loads that may interact with vehicle structural dynamic modes and the guidance and control of the vehicle. Buffet testing of a 3-percent scale rigid buffet model of a launch vehicle cargo configuration with a tangent-ogive payload shroud was conducted in 2012 and in 2016. Initial buffet forcing functions (BFFs) utilized a coarse pressure sensor distribution on the vehicle surface in which a single longitudinal station with eight sensors observed the terminal shock environment at Mach 0.90. An examination of these circumferential pressures reveal large impulse-like pressure fluctuations and an asymmetry in pressure when the vehicle is at a nonzeroangle of attack that result in high BFFs. Revisions to the shock integration region were made based on computational fluid dynamics and shadowgraph video of shock motion to better represent the BFFs and reduce the high loads resulting from this environment. To more clearly understand this terminal shock environment, a second wind tunnel test was conducted with a dense distribution of 256 sensors at the terminal shock location. These sensor arrays presents a unique opportunity to observe the unsteady terminal shock environment and to characterize the impact of various integration schemes on the BFFs. This paper presents a summary of the development of BFFs for this terminal shock and a detailed analyses of shock region pressure coefficients, coherence, BFFs, shock location time histories, and power spectral density to help guide development of BFFs for other launch vehicle test and analysis programs.

Piatak, David J.↗