Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximate boundary”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Turbulent mixing layers in the interstellar medium of galaxies

We propose that turbulent mixing layers are common in the interstellar medium (ISM). Injection of kinetic energy into the ISM by supernovae and stellar winds, in combination with density and temperature inhomogeneities, results in shear flows. Such flows will become turbulent due to the high Reynolds number (low viscosity) of the ISM plasma. These turbulent boundary layers will be particularly interesting where the shear flow occurs at boundaries of hot (approximately 10(exp 6) K) and cold or warm (10(exp 2) - 10(exp 4) K) gas. Mixing will occur in such layers producing intermediate-temperature gas at T is approximately equal to 10(exp 5.0) - 10(exp 5.5) that radiates strongly in the optical, ultraviolet, and EUV. We have modeled these layers under the assumptions of rapid mixing down to the atomic level and steady flow. By including the effects of non-equilibrium ionization and self-photoionization of the gas as it cools after mixing, we predict the intensities of numerous optical, infrared, and ultraviolet emission lines, as well as absorption column densities of C 4, N 5, Si 4, and O 6.

Slavin, J. D.↗
NASA NTRS1993

Aerothermodynamics of a MESUR Martian entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley Aerothermodynamics Upwind Relaxation Algorithm (LAURA) for Martian Atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Martian-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Martian-entry conditions characteristic of the Mars Environmental Survey (MESUR) mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/sq cm. Radiative heating is negligible.

Mitcheltree, R. A.↗
NASA NTRS1993

Sun-perturbed Earth-to-moon transfers with ballistic capture

A method is described for constructing a new type of low energy transfer trajectory from the Earth to the moon leading to ballistic capture. This is accomplished by utilizing the nonlinear Earth-moon-sun perturbations on a point mass in three dimensions. The interaction of the gravitational fields of the bodies defines transition regions in the position-velocity space where the dynamic effects on the point mass tend to balance. These are termed weak stability boundaries. The transfer is obtained by the use of trajectories connecting the weak stability boundaries. It uses approximately 18 percent less Delta-V than the Hohmann transfer to insert a spacecraft into a circular orbit about the moon. The use of this transfer has recently been demonstrated by Japan's Hiten spacecraft, which arrived at the moon on October 2, 1991. Application of the transfer method is also made to the Lunar Observer Mission.

Belbruno, Edward A.↗
NASA NTRS1948

The effects of aerodynamic heating and heat transfer on the surface temperature of a body of revolution in steady supersonic flight

An approximate method for determining the convective cooling requirement in the laminar boundary-layer region of a body of revolution in high-speed flight was developed and applied to an example body. The cooling requirement for the example body was determined as a function of Mach number, altitude, size, and a surface-temperature parameter. The maximum value of Mach number considered was 3.0 and the altitudes considered were those within the lower constant-temperature region of the atmosphere (40,000 to 120,000 ft.). The extent of the laminar boundary layer was determined approximately at each condition as a function of the variables considered.

Scherrer, Richard↗
NASA NTRS1994

Aerothermodynamic methods for a Mars environmental survey Mars entry

Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley aerothermodynamics upwind relaxation algorithm (LAURA) for Mars atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Mars-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Mars-entry conditions characteristics of the Mars environmental survey mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/cm(exp 2). Radiative heating is negligible.

AEROTHERMODYNAMIC METHODS↗
NASA NTRS1995

Turbulence measurements in hypersonic boundary layers using constant-temperature anemometry and Reynolds stress measurements in hypersonic boundary layers

The primary objective in the two research investigations performed under NASA Langley sponsorship (Turbulence measurements in hypersonic boundary layers using constant temperature anemometry and Reynolds stress measurements in hypersonic boundary layers) has been to increase the understanding of the physics of hypersonic turbulent boundary layers. The study began with an extension of constant-temperature thermal anemometry techniques to a Mach 11 helium flow, including careful examinations of hot-wire construction techniques, system response, and system calibration. This was followed by the application of these techniques to the exploration of a Mach 11 helium turbulent boundary layer (To approximately 290 K). The data that was acquired over the course of more than two years consists of instantaneous streamwise mass flux measurements at a frequency response of about 500 kHz. The data are of exceptional quality in both the time and frequency domain and possess a high degree of repeatability. The data analysis that has been performed to date has added significantly to the body of knowledge on hypersonic turbulence, and the data reduction is continuing. An attempt was then made to extend these thermal anemometry techniques to higher enthalpy flows, starting with a Mach 6 air flow with a stagnation temperature just above that needed to prevent liquefaction (To approximately 475 F). Conventional hot-wire anemometry proved to be inadequate for the selected high-temperature, high dynamic pressure flow, with frequent wire breakage and poor system frequency response. The use of hot-film anemometry has since been investigated for these higher-enthalpy, severe environment flows. The difficulty with using hot-film probes for dynamic (turbulence) measurements is associated with construction limitations and conduction of heat into the film substrate. Work continues under a NASA GSRP grant on the development of a hot film probe that overcomes these shortcomings for hypersonic flows. Each of the research tasks performed during the NASA Langley research grants is discussed separately below.

Spina, Eric F.↗
NASA NTRS1997

Multilevel methods for elliptic problems on unstructured grids

Multilevel methods on unstructured grids for elliptic problems are reviewed. The advantages of these techniques are the flexible approximation of the boundaries of complicated physical domains and the ability to adapt the grid to the resolution of fine scaled structures. Multilevel methods, which include multigrid methods and domain decomposition methods, depend on the correct splitting of appropriate finite element spaces. The standard splittings used in the structured grid case cannot be directly extended to unstructured grids due to their requirement for a hierarchical grid structure. Issues related to the application of multilevel methods to unstructured grids are discussed, including how the coarse spaces and transfer operators are defined and how different types of boundary conditions are treated. An obvious way to generate a coarse mesh is to regrid the physical domain several times. Several alternatives are proposed and discussed: node nested coarse spaces, agglomerated coarse spaces and algebraically generated coarse spaces.

Chan, Tony F.↗
NASA NTRS1959

On Fully Developed Channel Flows: Some Solutions and Limitations, and Effects of Compressibility, Variable Properties, and Body Forces

An examination of the effects of compressibility, variable properties, and body forces on fully developed laminar flow has indicated several limitations on such streams. In the absence of a pressure gradient, but presence of a body force (e.g., gravity), an exact fully developed gas flow results. For a liquid this follows also for the case of a constant streamwise pressure gradient. These motions are exact in the sense of a Couette flow. In the liquid case two solutions (not a new result) can occur for the same boundary conditions. An approximate analytic solution was found which agrees closely with machine calculations.In the case of approximately exact flows, it turns out that for large temperature variations across the channel the effects of convection (due to, say, a wall temperature gradient) and frictional heating must be negligible. In such a case the energy and momentum equations are separated, and the solutions are readily obtained. If the temperature variations are small, then both convection effects and frictional heating can consistently be considered. This case becomes the constant-property incompressible case (or quasi-incompressible case for free-convection flows) considered by many authors. Finally there is a brief discussion of cases wherein streamwise variations of all quantities are allowed but only a such form that independent variables are separable. For the case where the streamwise velocity varies inversely as the square root distance along the channel a solution is given.

Maslen, Stephen H.↗
NASA NTRS1999

Cryogenic Insulation System for Soft Vacuum

The development of a cryogenic insulation system for operation under soft vacuum is presented in this paper. Conventional insulation materials for cryogenic applications can be divided into three levels of thermal performance, in terms of apparent thermal conductivity [k-value in milliwatt per meter-kelvin (mW/m-K)]. System k-values below 0.1 can be achieved for multilayer insulation operating at a vacuum level below 1 x 10(exp -4) torr. For fiberglass or powder operating below 1 x 10(exp -3) torr, k-values of about 2 are obtained. For foam and other materials at ambient pressure, k-values around 30 are typical. New industry and aerospace applications require a versatile, robust, low-cost thermal insulation with performance in the intermediate range. The target for the new composite insulation system is a k-value below 4.8 mW/m-K (R-30) at a soft vacuum level (from 1 to 10 torr) and boundary temperatures of approximately 77 and 293 kelvin (K). Many combinations of radiation shields, spacers, and composite materials were tested from high vacuum to ambient pressure using cryostat boiloff methods. Significant improvement over conventional systems in the soft vacuum range was demonstrated. The new layered composite insulation system was also shown to provide key benefits for high vacuum applications as well.

Augustynowicz, S. D.↗
NASA NTRS2001

Aerogel Beads as Cryogenic Thermal Insulation System

An investigation of the use of aerogel beads as thermal insulation for cryogenic applications was conducted at the Cryogenics Test Laboratory of NASA Kennedy Space Center. Steady-state liquid nitrogen boiloff methods were used to characterize the thermal performance of aerogel beads in comparison with conventional insulation products such as perlite powder and multilayer insulation (MLI). Aerogel beads produced by Cabot Corporation have a bulk density below 100 kilograms per cubic meter (kg/cubic m) and a mean particle diameter of 1 millimeter (mm). The apparent thermal conductivity values of the bulk material have been determined under steady-state conditions at boundary temperatures of approximately 293 and 77 kelvin (K) and at various cold vacuum pressures (CVP). Vacuum levels ranged from 10(exp -5) torr to 760 torr. All test articles were made in a cylindrical configuration with a typical insulation thickness of 25 mm. Temperature profiles through the thickness of the test specimens were also measured. The results showed the performance of the aerogel beads was significantly better than the conventional materials in both soft-vacuum (1 to 10 torr) and no-vacuum (760 torr) ranges. Opacified aerogel beads performed better than perlite powder under high-vacuum conditions. Further studies for material optimization and system application are in progress.

Fesmire, J. E.↗
NASA NTRS2004

Experimental, Theoretical, and Computational Investigation of Separated Nozzle Flows

A detailed experimental, theoretical, and computational study of separated nozzle flows has been conducted. Experimental testing was performed at the NASA Langley 16-Foot Transonic Tunnel Complex. As part of a comprehensive static performance investigation, force, moment, and pressure measurements were made and schlieren flow visualization was obtained for a sub-scale, non-axisymmetric, two-dimensional, convergent- divergent nozzle. In addition, two-dimensional numerical simulations were run using the computational fluid dynamics code PAB3D with two-equation turbulence closure and algebraic Reynolds stress modeling. For reference, experimental and computational results were compared with theoretical predictions based on one-dimensional gas dynamics and an approximate integral momentum boundary layer method. Experimental results from this study indicate that off-design overexpanded nozzle flow was dominated by shock induced boundary layer separation, which was divided into two distinct flow regimes; three- dimensional separation with partial reattachment, and fully detached two-dimensional separation. The test nozzle was observed to go through a marked transition in passing from one regime to the other. In all cases, separation provided a significant increase in static thrust efficiency compared to the ideal prediction. Results indicate that with controlled separation, the entire overexpanded range of nozzle performance would be within 10% of the peak thrust efficiency. By offering savings in weight and complexity over a conventional mechanical exhaust system, this may allow a fixed geometry nozzle to cover an entire flight envelope. The computational simulation was in excellent agreement with experimental data over most of the test range, and did a good job of modeling internal flow and thrust performance. An exception occurred at low nozzle pressure ratios, where the two-dimensional computational model was inconsistent with the three-dimensional separation observed in the experiment. In general, the computation captured the physics of the shock boundary layer interaction and shock induced boundary layer separation in the nozzle, though there were some differences in shock structure compared to experiment. Though minor, these differences could be important for studies involving flow control or thrust vectoring of separated nozzles. Combined with other observations, this indicates that more detailed, three-dimensional computational modeling needs to be conducted to more realistically simulate shock-separated nozzle flows.

Hunter, Craig A.↗
NASA NTRS2008

Unsteady Analysis of Blade and Tip Heat Transfer as Influenced by the Upstream Momentum and Thermal Wakes

The effect of the upstream wake on the blade heat transfer has been numerically examined. The geometry and the flow conditions of the first stage turbine blade of GE s E3 engine with a tip clearance equal to 2 percent of the span was utilized. Based on numerical calculations of the vane, a set of wake boundary conditions were approximated, which were subsequently imposed upon the downstream blade. This set consisted of the momentum and thermal wakes as well as the variation in modeled turbulence quantities of turbulence intensity and the length scale. Using a one-blade periodic domain, the distributions of unsteady heat transfer rate on the turbine blade and its tip, as affected by the wake, were determined. Such heat transfer coefficient distribution was computed using the wall heat flux and the adiabatic wall temperature to desensitize the heat transfer coefficient to the wall temperature. For the determination of the wall heat flux and the adiabatic wall temperatures, two sets of computations were required. The results were used in a phase-locked manner to compute the unsteady or steady heat transfer coefficients. It has been found that the unsteady wake has some effect on the distribution of the time averaged heat transfer coefficient on the blade and that this distribution is different from the distribution that is obtainable from a steady computation. This difference was found to be as large as 20 percent of the average heat transfer on the blade surface. On the tip surface, this difference is comparatively smaller and can be as large as four percent of the average.

Ameri, Ali A.↗
NASA NTRS2008

Saturn's Auroral Response to the Solar Wind: Centrifugal Instability Model

We describe a model initially presented by Sittler et al. [2006] which attempts to explain the global response of Saturn's magnetosphere and its corresponding auroral behavior to variations in the solar wind. The model was derived from published simultaneous Hubble Space Telescope (HST) auroral images and Cassini upstream measurements taken during the month of January 2004. These observations show a direct correlation between solar wind dynamic pressure and (1) auroral brightening toward dawn local time, (2) an increase of rotational movement of auroral features to as much as 75% of the corotation speed, (3) the movement of the auroral oval to higher latitudes and (4) an increase in the intensity of Saturn Kilometric Radiation (SKR). This model is an alternative to the reconnection model of Cowley et al. [2004a,b; 2005] which is more Earth-like while ours stresses rotation. If angular momentum is conserved in a global sense, then when compressed the magnetosphere will tend to spin up and when it expands will tend to spin down. With the plasma sheet outer boundary at L approximates 15 we argue this region to be the dominant source region for the precipitating particles. If radial transport is dominated by centrifugal driven flux tube interchange motions, then when the magnetosphere spins up, outward transport will increase, the precipitating particles will move radially outward and cause the auroral oval to move to higher latitudes as observed. The Kelvin-Helmholtz instability may contribute to the enhanced emission along the dawn meridian as observed by HST. We present this model in the context of presently published observations by Cassini.

Sittler, Edward C.↗
NASA NTRS2011

Inflow/Outflow Boundary Conditions with Application to FUN3D

Several boundary conditions that allow subsonic and supersonic flow into and out of the computational domain are discussed. These boundary conditions are demonstrated in the FUN3D computational fluid dynamics (CFD) code which solves the three-dimensional Navier-Stokes equations on unstructured computational meshes. The boundary conditions are enforced through determination of the flux contribution at the boundary to the solution residual. The boundary conditions are implemented in an implicit form where the Jacobian contribution of the boundary condition is included and is exact. All of the flows are governed by the calorically perfect gas thermodynamic equations. Three problems are used to assess these boundary conditions. Solution residual convergence to machine zero precision occurred for all cases. The converged solution boundary state is compared with the requested boundary state for several levels of mesh densities. The boundary values converged to the requested boundary condition with approximately second-order accuracy for all of the cases.

Carlson, Jan-Renee↗
NASA NTRSDate not supplied

Thermal Control in Hypersonic Leading Edges Using Liquid Metal High Temperature Oscillating Heat Pipes

The desire to increase the range and speed of hypersonic vehicles requires sharp, shape-stable Wing Leading Edges (WLE) with performance capabilities well beyond the current state of the art. The concentrated heat flux on sharp leading edges can lead to runaway thermal and mechanical failures. Heat pipes provide a passive solution to minimize the peak temperature, thermal gradients, and resulting thermal stresses near the WLE. Conventional heat pipes have been shown to address these issues at conditions as high as Mach 8. Oscillating heat pipes (OHPs) could be used to extend this benefit to missions with higher Mach numbers since their heat transport capacity typically exceeds that of conventional heat pipes. In the current investigation, several identical liquid metal high-temperature oscillating heat pipes (LMHOHPs) were constructed and tested to demonstrate the concept. The LMHOHPs had six turns, internal diameters of 2 mm, and total lengths of 100 mm. They were additively manufactured from C103 niobium alloy and charged with lithium to reach a fill ratio of approximately 50%. The boundary conditions of the experiment were designed to resemble those of various hypersonic flight profiles. Testing was performed with a system consisting of flow-controlled acetylene torch heating and radiation dominated heat rejection. The LMHOHP’s performance was measured using temperature data from an infrared camera and pyrometer. The experimental results were used to evaluate the capability of C103/lithium LMHOHPs under various flight regimes. A heat transfer model of the LMHOHP was proposed and validated against test data. Insights from the model were then extrapolated to evaluate a Tungsten-RHC/Gallium LMHOHP with anticipated capabilities beyond that of the C103/lithium system.

Max Pawlick↗
NASA NTRSDate not supplied

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The planned microgravity experiment Growth and Extinction Limit (GEL) is designed to study the material flammability and extinction of thermally thick solid with an emphasis to quantify the effect of conductive heat loss into the solid interior by measuring the in-depth temperature profile. To obtain a more detailed understanding of the extinction phenomena, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is expanded to account for the excessive heat loss into the solid beyond the steady burning state. This is important because for thick solids, this is important because the burning period can have a long and slow solid transient with a quasi-steady gas -phase. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper and extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The classical U-shaped flammability boundary in the oxygen vs. stretch rate plot now becomes a series of U-shaped boundaries as a function of Ψ. For a given oxygen ambient, extinction boundary is a function of Ψ. The boundary consists of a blowoff branch and a radiative quenching branch. This is similar to the steady burning flammability boundary (Ψ =0) in the oxygen vs stretch rate plot. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, and reaction rate are presented along the extinction boundaries. In particular, we note that the non-dimensional gas reaction rate is the ratio of diffusion reaction time to the reaction diffusion time, so it has the classical Damhkoler number character. By examining the maximum non-dimensional reaction rates, we observe that along the extinction boundaries, they are approximate constant across the entire range of stretch rates (from quenching to blow-off). These critical values are of the order unity but may vary in different oxygen ambient. The implication will be discussed in the presentation together with the near-limit flame structures.

Chengyao Li↗
NASA NTRS1977

On the stability of Galerkin methods for initial-boundary value problems for hyperbolic systems

The stability of approximating the solution of mixed initial-boundary value problems for hyperbolic systems by semidiscrete Galerkin methods is studied. It is shown that a particular straightforward Galerkin method yields an unstable approximation, and that this numerical instability is caused by an improper treatment of the boundary. Stable schemes are then presented, one of which differs from the unstable scheme only insofar as the treatment of the boundary is concerned. These stable schemes make use of a particular matrix which symmetrizes the differential system. It is therefore shown that the use of this matrix is crucial to the stability of the computations as well as for obtaining a priori bounds on the energy of the continuous system. This symmetrizing matrix is also related to the diagonalizing matrix for the system of hyperbolic equations and to the Liapunov matrix for the system of ordinary differential equations resulting from the application of Galerkin's method.

Gunzburger, M. D.↗
NASA NTRS2018

Using Adaptive Mesh Refinement to Study Grid Resolution Effects for Shock-Boundary Layer Interactions

Adaptive Mesh Refinement (AMR) promises a much more computationally efficient means to obtain a discrete approximation to a continuous boundary value problem of a specified accuracy than classic isotropic grid refinement. The AMR capability of OVERFLOW (a computational fluid dynamics (CFD) code) is utilized to provide estimates of the exact analytical solutions to problems of interest to turbulence modeling. Predictions of surface pressure and skin friction, essentially the state of stress at the surface, shows little difference with grids believed to be "grid resolved." Velocity profiles, on the other hand, show marked differences in flows with shocks. The AMR method, as implemented in OVERFLOW 2.2k, appears to provide the ability to produce arbitrarily accurate solutions at a predictable cost much smaller than classic uniform mesh refinement.

Turbulence Modeling↗