Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable 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 1,153 records · Page 64

Experimental and analytical study of cryogenic propellant boiloff to develop and verify alternate pressurization concepts for Space Shuttle external tank using a scaled down tank

Self pressurization by propellant boiloff is experimentally studied as an alternate pressurization concept for the Space Shuttle external tank (ET). The experimental setup used in the study is an open flow system which is composed of a variable area test tank and a recovery tank. The vacuum jacketed test tank is geometrically similar to the external LOx tank for the Space Shuttle. It is equipped with instrumentation to measure the temperature and pressure histories within the liquid and vapor, and viewports to accommodate visual observations and Laser-Doppler Anemometry measurements of fluid velocities. A set of experiments were conducted using liquid Nitrogen to determine the temperature stratification in the liquid and vapor, and pressure histories of the vapor during sudden and continuous depressurization for various different boundary and initial conditions. The study also includes the development and calibration of a computer model to simulate the experiments. This model is a one-dimensional, multi-node type which assumes the liquid and the vapor to be under non-equilibrium conditions during the depressurization. It has been tested for a limited number of cases. The preliminary results indicate that the accuracy of the simulations is determined by the accuracy of the heat transfer coefficients for the vapor and the liquid at the interface which are taken to be the calibration parameters in the present model.

Akyuzlu, K. M.↗

Robust Airfoil Optimization in High Resolution Design Space

The robust airfoil shape optimization is a direct method for drag reduction over a given range of operating conditions and has three advantages: (1) it prevents severe degradation in the off-design performance by using a smart descent direction in each optimization iteration, (2) it uses a large number of B-spline control points as design variables yet the resulting airfoil shape is fairly smooth, and (3) it allows the user to make a trade-off between the level of optimization and the amount of computing time consumed. The robust optimization method is demonstrated by solving a lift-constrained drag minimization problem for a two-dimensional airfoil in viscous flow with a large number of geometric design variables. Our experience with robust optimization indicates that our strategy produces reasonable airfoil shapes that are similar to the original airfoils, but these new shapes provide drag reduction over the specified range of Mach numbers. We have tested this strategy on a number of advanced airfoil models produced by knowledgeable aerodynamic design team members and found that our strategy produces airfoils better or equal to any designs produced by traditional design methods.

Li, Wu↗

Extinction Criteria for Opposed-Flow Flame Spread in a Microgravity Environment

A simplified analysis is presented to extend a previous work on flame extinction in a quiescent microgravity environment to a more likely situation of a mild opposing flow. The energy balance equation, that includes surface re-radiation, is solved to yield a closed form spread rate expression in terms of its thermal limit, and a radiation number that can be evaluated from the known parameters of the problem. Based on this spread rate expression, extinction criterions for a flame over solid fuels, both thin and thick, have been developed that are qualitatively verified with experiments conducted at the MGLAB in Japan. Flammability maps with oxygen level, opposing flow velocity and fuel thickness as independent variables are extracted from the theory that explains the well-established trends in the existing experimental data.

Bhattacharjee, Subrata↗

Adjoint Formulation for an Embedded-Boundary Cartesian Method

Many problems in aerodynamic design can be characterized by smooth and convex objective functions. This motivates the use of gradient-based algorithms, particularly for problems with a large number of design variables, to efficiently determine optimal shapes and configurations that maximize aerodynamic performance. Accurate and efficient computation of the gradient, however, remains a challenging task. In optimization problems where the number of design variables dominates the number of objectives and flow- dependent constraints, the cost of gradient computations can be significantly reduced by the use of the adjoint method. The problem of aerodynamic optimization using the adjoint method has been analyzed and validated for both structured and unstructured grids. The method has been applied to design problems governed by the potential, Euler, and Navier-Stokes equations and can be subdivided into the continuous and discrete formulations. Giles and Pierce provide a detailed review of both approaches. Most implementations rely on grid-perturbation or mapping procedures during the gradient computation that explicitly couple changes in the surface shape to the volume grid. The solution of the adjoint equation is usually accomplished using the same scheme that solves the governing flow equations. Examples of such code reuse include multistage Runge-Kutta schemes coupled with multigrid, approximate-factorization, line-implicit Gauss-Seidel, and also preconditioned GMRES. The development of the adjoint method for aerodynamic optimization problems on Cartesian grids has been limited. In contrast to implementations on structured and unstructured grids, Cartesian grid methods decouple the surface discretization from the volume grid. This feature makes Cartesian methods well suited for the automated analysis of complex geometry problems, and consequently a promising approach to aerodynamic optimization. Melvin e t al. developed an adjoint formulation for the TRANAIR code, which is based on the full-potential equation with viscous corrections. More recently, Dadone and Grossman presented an adjoint formulation for the Euler equations. In both approaches, a boundary condition is introduced to approximate the effects of the evolving surface shape that results in accurate gradient computation.

Nemec, Marian↗

Base Heating Sensitivity Study for a 4-Cluster Rocket Motor Configuration in Supersonic Freestream

In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.

Mehta, Manish↗

Radiation-Spray Coupling for Realistic Flow Configurations

Three Large Eddy Simulations (LES) for a lean-direct injection (LDI) combustor are performed and compared. In addition to the cold flow simulation, the effect of radiation coupling with the multi-physics reactive flow is analyzed. The flame let progress variable approach is used as a subgrid combustion model combined with a stochastic subgrid model for spray atomization and an optically thin radiation model. For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is utilized. To achieve realistic inflow, a simple recycling technique is performed at the inflow section upstream of the swirler. Good comparison is shown with the experimental data mean and root mean square profiles. The effect of combustion is found to change the shape and size of the central recirculation zone. Radiation is found to change the spray dynamics and atomization by changing the heat release distribution and the local temperature values impacting the evaporation process. The simulation with radiation modeling shows wider range of droplet size distribution by altering the evaporation rate. The current study proves the importance of radiation modeling for accurate prediction in realistic spray combustion configurations, even for low pressure systems.

El-Asrag, Hossam↗

Transient Liquefaction on the Lunar or Martian Surface Operational Demonstration

The Cryogenic Fluid In-situ Liquefaction for Landers (CryoFILL) activity has been investigating concepts for the liquefaction of cryogenic fluids produced on the surface of the Moon and Mars. Liquefaction testing of oxygen in a scaled lander tank using an integrated, industrial cryocooler system was completed. The testing covered the determination of the nominal performance operation of the system, constant liquefaction performance, and transient liquefaction performance. The goal of the testing was to demonstrate cryogenic liquefaction operational capabilities on the Lunar and Martian surfaces for landers and In Situ Resource Utilization. The testing met all key performance parameter threshold values including liquefaction rate (demonstrated 1.6 kg/hr compared to threshold of 1.1 kg/hr), number of variables as transient (goal of 3: liquefaction flow, environmental temperature, and cryocooler input power – all demonstrated), and number of fill levels (goal of 3: < 5%, 50%, and 90% - all demonstrated). Key testing results from oxygen liquefaction testing will be discussed including sensitivities provided by the analysis team.

Liquefaction↗

Advanced supersonic propulsion study, phase 2

A continuation of the NASA/P and WA study to evaluate various types of propulsion systems for advanced commercial supersonic transports has resulted in the identification of two very promising engine concepts. They are the Variable Stream Control Engine which provides independent temperature and velocity control for two coannular exhaust streams, and a derivative of this engine, a Variable Cycle Engine that employs a rear flow-inverter valve to vary the bypass ratio of the cycle. Both concepts are based on advanced engine technology and have the potential for significant improvements in jet noise, exhaust emissions and economic characteristics relative to current technology supersonic engines. Extensive research and technology programs are required in several critical areas that are unique to these supersonic Variable Cycle Engines to realize these potential improvements. Parametric cycle and integration studies of conventional and Variable Cycle Engines are reviewed, features of the two most promising engine concepts are described, and critical technology requirements and required programs are summarized.

Howlett, R. A.↗

High frequency sound attenuation in short flow ducts

A geometrical acoustics approach is proposed as a practical design tool for absorbent liners in such short flow ducts as may be found in turbofan engine nacelles. As an example, a detailed methodology is presented for three different types of sources in a parallel plate duct containing uniform ambient flow. A plane wave whose wavefronts are not normal to the duct walls, an arbitrarily located point source, and a spatially harmonic line source are each considered. Optimal wall admittance distributions are found, and it is shown how to estimate the insertion loss for any admittance distribution. The extension of the methodology to realistic source distributions in variable area cylindrical or annular ducts containing arbitrary flow is shown to be conceptually straightforward and computationally practical on a vector-hardware digital computer.

Posey, J. W.↗

Validation of the geostrophic method for estimating zonal currents at the equator from Geosat altimeter data

The applicability of satellite altimeter data for estimating zonal current variability at the equator is assessed using the meriodionally differenced form of the geostrophic balance. Estimates of geostrophic zonal flow anomalies in the equatorial Pacific have been deduced from 17-day collinear altimeter data during the first year of the Geosat Exact Repeat Mission. Altimeter-derived geostrophic estimates agree well with in situ zonal current variability. Comparison of flow-frequency near-surface zonal current observed from equatorial moorings at 165 deg E, 140 deg W, and 110 deg W yield correlations of 0.83, 0.85, and 0.51, respectively, with a mean rms difference of 23 cm/sec. The inclusion of up to 11 ascending and descending Geosat tracks within the 9-deg band for every 17-day repeat effectively reduced the temporal sampling interval to 1.5 days at 165 deg E and 140 deg W. The 6.8-km along track spacing of the altimeter measurements provides sufficient resolution for the effective filtering of small-scale meridional noise, both instrumental and oceanic.

Picaut, Joel↗

Time-derivative preconditioning for viscous flows

A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.

Choi, Yunho↗

Time-derivative preconditioning for viscous flows

A time-derivative preconditioning algorithm that is effective over a wide range of flow conditions from inviscid to very diffusive flows and from low speed to supersonic flows was developed. This algorithm uses a viscous set of primary dependent variables to introduce well-conditioned eigenvalues and to avoid having a nonphysical time reversal for viscous flow. The resulting algorithm also provides a mechanism for controlling the inviscid and viscous time step parameters to be of order one for very diffusive flows, thereby ensuring rapid convergence at very viscous flows as well as for inviscid flows. Convergence capabilities are demonstrated through computation of a wide variety of problems.

Choi, Yunho↗

Choked flow effects in the NSI driven pin puller

This paper presents an analysis for pyrotechnic combustion and pin motion in the NASA Standard Initiator (NSI) actuated pin puller. The conservation principles and constitutive relations for a multi-phase system are posed and reduced to a set of eight ordinary differential equations which are solved to predict the system performance. The model tracks the interactions of the unreacted, incompressible solid pyrotechnic, incompressible condensed phase combustion products, and gas phase combustion products. The model accounts for multiple pyrotechnic grains, variable burn surface area, and combustion product mass flow rates through an orifice located within the device. Pressure-time predictions compare favorably with experimental data. Results showing model sensitivity to changes in the cross-sectional area of the orifice are presented.

Gonthier, Keith A.↗

The Sensitivity of Hyporheic Exchange to Fractal Properties of Riverbeds

Hyporheic exchange in riverbeds is driven by current-bed topography interactions. Because riverbeds exhibit topographic roughness across scales, from individual grains to bedforms and bars, they can exhibit fractal patterns. This study analyzed the influence of fractal properties of riverbed topography on hyporheic exchange. A set of synthetic fractal riverbeds with different scaling statistics was used as inputs to sequentially coupled numerical simulations of turbulent channel flow and hyporheic flow. In the analysis, the maximum power spectrum (dune size) and the fractal dimension (topographic complexity) were considered as independent variables and we then investigated how interfacial fluxes and hyporheic travel times are functionally related to these variables. As the maximum power spectrum increases (i.e., dune height to flow depth ratio), the average interfacial flux increases logarithmically whereas it increases exponentially with an increase in fractal dimension. Hyporheic exchange is more sensitive to additional roughness (larger fractal dimensions) than to bedform size (larger maximum power). Our results imply that fractal properties of riverbeds are crucial to predicting hyporheic exchange. The predictive relationships we propose could be integrated with reduced complexity, large-scale models. Further, they can also be used to design artificial topographies that target hyporheic ecosystem services.

54 ENVIRONMENTAL SCIENCES↗

Modal element method for potential flow in non-uniform ducts: Combining closed form analysis with CFD

An analytical procedure is presented, called the modal element method, that combines numerical grid based algorithms with eigenfunction expansions developed by separation of variables. A modal element method is presented for solving potential flow in a channel with two-dimensional cylindrical like obstacles. The infinite computational region is divided into three subdomains; the bounded finite element domain, which is characterized by the cylindrical obstacle and the surrounding unbounded uniform channel entrance and exit domains. The velocity potential is represented approximately in the grid based domain by a finite element solution and is represented analytically by an eigenfunction expansion in the uniform semi-infinite entrance and exit domains. The calculated flow fields are in excellent agreement with exact analytical solutions. By eliminating the grid surrounding the obstacle, the modal element method reduces the numerical grid size, employs a more precise far field boundary condition, as well as giving theoretical insight to the interaction of the obstacle with the mean flow. Although the analysis focuses on a specific geometry, the formulation is general and can be applied to a variety of problems as seen by a comparison to companion theories in aeroacoustics and electromagnetics.

Baumeister, Kenneth J.↗

Lagrangian and Dirac constraints for the ideal incompressible fluid and magnetohydrodynamics

The incompressibility constraint for fluid flow was imposed by Lagrange in the so-called Lagrangian variable description using his method of multipliers in the Lagrangian (variational) formulation. An alternative is the imposition of incompressibility in the Eulerian variable description by a generalization of Dirac’s constraint method using noncanonical Poisson brackets. Here it is shown how to impose the incompressibility constraint using Dirac’s method in terms of both the canonical Poisson brackets in the Lagrangian variable description and the noncanonical Poisson brackets in the Eulerian description, allowing for the advection of density. Both cases give the dynamics of infinite-dimensional geodesic flow on the group of volume preserving diffeomorphisms and explicit expressions for this dynamics in terms of the constraints and original variables is given. Since Lagrangian and Eulerian conservation laws are not identical, comparison of the various methods is made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Boundary conditions and the simulation of low Mach number flows

The problem of accurately computing low Mach number flows, with the specific intent of studying the interaction of sound waves with incompressible flow structures, such as concentrations of vorticity is considered. This is a multiple time (and/or space) scales problem, leading to various difficulties in the design of numerical methods. Concentration is on one of these difficulties - the development of boundary conditions at artificial boundaries which allow sound waves and vortices to radiate to the far field. Nonlinear model equations are derived based on assumptions about the scaling of the variables. Then these are linearized about a uniform flow and exact boundary conditions are systematically derived using transform methods. Finally, useful approximations to the exact conditions which are valid for small Mach number and small viscosity are computed.

Hagstrom, Thomas↗

3-D Geologic Controls of Hydrothermal Fluid Flow at Brady Geothermal Field, Nevada using PCA

In many hydrothermal systems, fracture permeability along faults provides pathways for groundwater to transport heat from depth. Faulting generates a range of deformation styles that cross-cut heterogeneous geology, resulting in complex patterns of permeability, porosity, and hydraulic conductivity. Vertical connectivity (a through going network of permeable areas that allows advection of heat from depth to the shallow subsurface) is rare and is confined to relatively small volumes that have highly variable spatial distribution. This local compartmentalization of connectivity represents a significant challenge to understanding hydrothermal circulation and for exploring, developing, and managing hydrothermal resources. Here, we present an evaluation of the geologic characteristics that control this compartmentalization in hydrothermal systems through 3-D analysis of the Brady geothermal field in western Nevada. A published 3-D geologic map of the Brady area is used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity. The 3-D distribution of these variables is compared to the distribution of productive and non-productive fluid flow intervals along production wells and non-productive wells via principal component analysis (PCA). This comparison elucidates which geologic and structural variables are most closely associated with productive fluid flow intervals. Results indicate that production intervals at Brady are located: (1) within or near to known and stress-loaded macro-scale faults, and (2) in areas of high fault and fracture density. This submission includes the published journal article detailing this work, the published 3-D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3-D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using existing R programs.

15 GEOTHERMAL ENERGY↗