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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Supersonic Wind Tunnel Tests of a Half-axisymmetric 12 Deg-spike Inlet to a Rocket-based Combined-cycle Propulsion System

Results of an isolated inlet test for NASA's GTX air-breathing launch vehicle concept are presented. The GTX is a Vertical Take-off/ Horizontal Landing reusable single-stage-to-orbit system powered by a rocket-based combined-cycle propulsion system. Tests were conducted in the NASA Glenn 1- by 1-Foot Supersonic Wind Tunnel during two entries in October 1998 and February 1999. Tests were run from Mach 2.8 to 6. Integrated performance parameters and static pressure distributions are reported. The maximum contraction ratios achieved in the tests were lower than predicted by axisymmetric Reynolds-averaged Navier-Stokes computational fluid dynamics (CFD). At Mach 6, the maximum contraction ratio was roughly one-half of the CFD value of 16. The addition of either boundary-layer trip strips or vortex generators had a negligible effect on the maximum contraction ratio. A shock boundary-layer interaction was also evident on the end-walls that terminate the annular flowpath cross section. Cut-back end-walls, designed to reduce the boundary-layer growth upstream of the shock and minimize the interaction, also had negligible effect on the maximum contraction ratio. Both the excessive turning of low-momentum comer flows and local over-contraction due to asymmetric end-walls were identified as possible reasons for the discrepancy between the CFD predictions and the experiment. It is recommended that the centerbody spike and throat angles be reduced in order to lessen the induced pressure rise. The addition of a step on the cowl surface, and planar end-walls more closely approximating a plane of symmetry are also recommended. Provisions for end-wall boundary-layer bleed should be incorporated.

DeBonis, J. R.↗
NASA NTRS2016

Simulation of Sweep-Jet Flow Control, Single Jet and Full Vertical Tail

This work is a simulation technology demonstrator, of sweep jets used to suppress boundary layer separation and increase maximum achievable load coefficients. A sweep jet is a discrete Coanda jet that oscillates in the plane parallel to an aerodynamic surface. It injects mass and momentum in the approximate stream wise direction. It also generate turbulent eddies at the oscillation frequency, which are typically large relative to boundary layer turbulence, and which augmenting mixing across the boundary layer to attack flow separation. Simulations of a fluidic oscillator, the sweep jet emerging from the oscillator, and the suppression of boundary layer separation by an array of sweep jets are performed. Simulation results are compared to data from a dedicated CFD validation experiment of a single oscillator and its sweep jet, and from a study of a full-scale Boeing 757 vertical tail augmented with an array of sweep jets.2, 20 A critical step in the work is the development of realistic time-dependent sweep-jet in flow boundary conditions, derived from the results of the single-oscillator simulations, which create the sweep jets in the full-tail simulations. Simulations were performed using the Over flow CFD solver, with high-order spatial discretization and a range of turbulence modeling. Good results were obtained for all flows simulated, when suitable turbulence modeling was used.

flow control↗
NASA NTRS1994

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Siegel, R.↗
NASA NTRS1975

Theoretical face pressure and drag characteristics of forward-facing steps in supersonic turbulent boundary layers

A theoretical investigation of the pressure distributions and drag characteristics was made for forward facing steps in turbulent flow at supersonic speeds. An approximate solution technique proposed by Uebelhack has been modified and extended to obtain a more consistent numerical procedure. A comparison of theoretical calculations with experimental data generally indicated good agreement over the experimentally available range of ratios of step height to boundary layer thickness from 7 to 0.05.

Patel, D. K.↗
NASA NTRS2013

Atmospheric Turbulence Estimates from a Pulsed Lidar

Estimates of the eddy dissipation rate (EDR) were obtained from measurements made by a coherent pulsed lidar and compared with estimates from mesoscale model simulations and measurements from an in situ sonic anemometer at the Denver International Airport and with EDR estimates from the last observation time of the trailing vortex pair. The estimates of EDR from the lidar were obtained using two different methodologies. The two methodologies show consistent estimates of the vertical profiles. Comparison of EDR derived from the Weather Research and Forecast (WRF) mesoscale model with the in situ lidar estimates show good agreement during the daytime convective boundary layer, but the WRF simulations tend to overestimate EDR during the nighttime. The EDR estimates from a sonic anemometer located at 7.3 meters above ground level are approximately one order of magnitude greater than both the WRF and lidar estimates - which are from greater heights - during the daytime convective boundary layer and substantially greater during the nighttime stable boundary layer. The consistency of the EDR estimates from different methods suggests a reasonable ability to predict the temporal evolution of a spatially averaged vertical profile of EDR in an airport terminal area using a mesoscale model during the daytime convective boundary layer. In the stable nighttime boundary layer, there may be added value to EDR estimates provided by in situ lidar measurements.

Pruis, Matthew J.↗
NASA NTRS2000

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗
NASA NTRS1992

Analytical shape sensitivities and approximations of modal response of generally laminated tapered skew plates

In the present determination of the derivatives of natural frequencies and mode shapes of a generally laminated tapered skew plate, with respect to various shape parameters, springs are used to simulate the essential boundary conditions. The independent shape parameters are plate surface area, aspect ratio, taper ratio, and sweep angle. Eigenvalues and eigenvectors are approximated over the range of the variable using linear, exponential, and pseudoexponential approximation schemes, and are compared with reanalysis-obtained values. Numerical results are presented for symmetrically and unsymmetrically laminated plates.

Singhvi, Sarvesh↗
NASA NTRS1983

Experimental investigation of a two-dimensional shock-turbulent boundary layer interaction with bleed

The two-dimensional interaction of an oblique shock wave with a turbulent boundary layer that included the effect of bleed was examined experimentally using a shock generator mounted across a supersonic wind tunnel The studies were performed at Mach numbers 2.5 and 2.0 and unit Reynolds number of approximately 2.0 x 10 to the 7th/meter. The study includes surface oil flow visualization, wall static pressure distributions and boundary layer pitot pressure profiles. In addition, the variation of the local bleed rates were measured. The results show the effect of the bleed on the boundary layer as well as the effect of the flow conditions on the local bleed rate.

Hingst, W. R.↗
NASA NTRS1983

Experimental investigation of a two-dimensional shock-turbulent boundary layer interaction with bleed

The two-dimensional interaction of an oblique shock wave with a turbulent boundary layer that included the effect of bleed was examined experimentally using a shock generator mounted across a supersonic wind tunnel. The studies were performed at Mach numbers 2.5 and 2.0 and unit Reynolds number of approximately 2.0 x 10 to the 7th/meter. The study includes surface oil flow visualization, wall static pressure distributions and boundary layer pitot pressure profiles. In addition, the variation of the local bleed rates were measured. The results show the effect of the bleed on the boundary layer as well as the effect of the flow conditions on the local bleed rate.

Hingst, W. R.↗
NASA NTRS1999

Static Performance of a Fixed-Geometry Exhaust Nozzle Incorporating Porous Cavities for Shock-Boundary Layer Interaction Control

An investigation was conducted in the model preparation area of the Langley 16-Foot Transonic Tunnel to determine the internal performance of a fixed-geometry exhaust nozzle incorporating porous cavities for shock-boundary layer interaction control. Testing was conducted at static conditions using a sub-scale nozzle model with one baseline and 27 porous configurations. For the porous configurations, the effects of percent open porosity, hole diameter, and cavity depth were determined. All tests were conducted with no external flow at nozzle pressure ratios from 1.25 to approximately 9.50. Results indicate that baseline nozzle performance was dominated by unstable, shock-induced, boundary-layer separation at over-expanded conditions. Porous configurations were capable of controlling off-design separation in the nozzle by either alleviating separation or encouraging stable separation of the exhaust flow. The ability of the porous nozzle concept to alternately alleviate separation or encourage stable separation of exhaust flow through shock-boundary layer interaction control offers tremendous off-design performance benefits for fixed-geometry nozzle installations. In addition, the ability to encourage separation on one divergent flap while alleviating it on the other makes it possible to generate thrust vectoring using a fixed-geometry nozzle.

Asbury, Scott C.↗
NASA NTRS2008

Role of Hf on Phase Formation in Ti45Zr(38-x)Hf(x)Ni17 Liquids and Solids

Hafnium and zirconium are very similar, with almost identical sizes and chemical bonding characteristics. However, they behave differently when alloyed with Ti and Ni. A sharp phase formation boundary near 18-21 at.% Hf is observed in rapidly-quenched and as-cast Ti45Zr38-xHfxNi17 alloys. Rapidly-quenched samples that contain less than 18 at.% Hf form the icosahedral quasicrystal phase, whiles samples containing more than 21 at.% form the 3/2 rational approximant phase. In cast alloys, a C14 structure is observed for alloys with Hf lower than the boundary concentration, while a large-cell (11.93 ) FCC Ti2Ni-type structure is found in alloys with Hf concentrations above the boundary. To better understand the role of Hf on phase formation, the structural evolution with supercooling and the solidification behavior of liquid Ti45Zr38-xHfxNi17 alloys (x=0, 12, 18, 21, 38) were studied using the Beamline Electrostatic Levitation (BESL) technique using 125keV x-rays on the 6ID-D beamline at the Advanced Photon Source, Argonne National Laboratory. For all liquids primary crystallization was to a BCC solid solution phase; interestly, an increase in Hf concentration leads to a decrease in the BCC lattice parameter in spite of the chemical similarity between Zr and Hf. A Reitveld analysis confirmed that as in the cast alloys, the secondary phase that formed was the C14 below the phase formation boundary and a Ti2Ni-type structure at higher Hf concentrations. Both the liquidus temperature and the reduced undercooling change sharply on traversing the phase formation boundary concentration, suggesting a change in the liquid structure. Structural information from a Honeycutt-Anderson index analysis of reverse Monte Carlo fits to the S(q) liquid data will be presented to address this issue.

Wessels, V.↗
NASA NTRS2012

MESSENGER and Mariner 10 Flyby Observations of Magnetotail Structure and Dynamics at Mercury

The first (M1), second (M2), and third (M3) MESSENGER flybys of Mercury traversed the planet's magnetotail from 1.25 to 3.25 RM downstream of the planet, where R(sub M) is Mercury's radius (2440 km). The encounters took place under northward, southward, and variable-polarity interplanetary magnetic field (IMF), respectively. The magnetic field strength B in Mercury's magnetotail follows a power law decrease with increasing antisunward distance ∣X∣, B approximately ∣X∣(sup G), with G varying from -5.4 for northward to -1.6 for southward IMF. Low-latitude boundary layers (LLBLs) containing strong northward magnetic field were detected at the tail flanks during two of the flybys. The observed thickness of the LLBL was 33% and 16% of the radius of the tail during M1 and M3, respectively, but the boundary layer was completely absent during M2. Clear signatures of tail reconnection are evident in the M2 and M3 magnetic field measurements. Plasmoids and traveling compression regions were observed during M2 and M3 with typical durations of approximately 1-3 s, suggesting diameters of approximately 500-1500 km. Overall, the response of Mercury's magnetotail to the steady southward IMF during M2 appeared very similar to steady magnetospheric convection events at Earth, which are believed to be driven by quasi-continuous reconnection. In contrast, the M3 measurements are dominated by tail loading and unloading events that resemble the large-scale magnetic field reconfigurations observed during magnetospheric substorms at Earth.

Magnetotail↗
NASA NTRS1986

Algebraic grid generation about wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated. Each grid has six exterior surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived from the intersection of planes with the Coons' patch definition. Single-valued functions relating approximate arc lengths along the grid curves to a computational coordinate define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grid.

Smith, R. E.↗
NASA NTRS2011

Two Azimuthally Separated Regions of Cusp Ion Injection Observed via Energetic Neutral Atoms

The low-energy neutral atom (LENA) imager on the IMAGE spacecraft can detect energetic neutral atoms produced by ion injection into the cusp through a charge exchange with the Earth's hydrogen exosphere. We examined the occurrence of the LENA cusp signal during positive IMF B(sub z) in terms of the arrival direction and the IMF clock angle theta(sub CA). Results of statistical analyses show that the occurrence frequency is high on the postnoon side when theta(sub CA) is between approximately 20 degrees and approximately 50 degrees. This is ascribed to ion injection caused by cusp reconnection typical of positive IMF B(sub z). Our results also show that there is another situation of high occurrence frequency, which can be identified with theta(sub CA) of approximately 30 degrees to approximately 80 degrees. When theta(sub CA) is relatively large (60 degrees - 80 degrees), occurrence frequencies are high at relatively low latitudes over a wide extent spanning both prenoon and postnoon sectors. This feature suggests that the ion injection is caused by reconnection at the dayside magnetopause. Its postnoon side boundary shifts toward the prenoon as theta(sub CA) decreases. When theta(sub CA) is less than approximately 50 degrees, the high occurrence frequency exists well inside the prenoon sector, which is azimuthally separated from the postnoon region ascribed to cusp reconnection. The prenoon region, which is thought due to ion injection caused by dayside reconnection, may explain the recent report that proton aurora brightening occurs in the unanticipated prenoon sector of the northern high-latitude ionosphere for IMF B(sub y) greater than 0 and B(sub z) greater than 0.

Abe, M.↗
NASA NTRS1999

Experimental Investigation of Pool Boiling Heat Transfer Enhancement in Microgravity in the Presence of Electric Fields

In boiling high heat fluxes are possible driven by relatively small temperature differences, which make its use increasingly attractive in aerospace applications. The objective of the research is to develop ways to overcome specific problems associated with boiling in the low gravity environment by substituting the buoyancy force with the electric force to enhance bubble removal from the heated surface. Previous studies indicate that in terrestrial applications nucleate boiling heat transfer can be increased by a factor of 50, as compared to values obtained for the same system without electric fields. The goal of our research is to experimentally explore the mechanisms responsible for EHD heat transfer enhancement in boiling in low gravity conditions, by visualizing the temperature distributions in the vicinity of the heated surface and around the bubble during boiling using real-time holographic interferometry (HI) combined with high-speed cinematography. In the first phase of the project the influence of the electric field on a single bubble is investigated. Pool boiling is simulated by injecting a single bubble through a nozzle into the subcooled liquid or into the thermal boundary layer developed along the flat heater surface. Since the exact location of bubble formation is known, the optical equipment can be aligned and focused accurately, which is an essential requirement for precision measurements of bubble shape, size and deformation, as well as the visualization of temperature fields by HI. The size of the bubble and the frequency of bubble departure can be controlled by suitable selection of nozzle diameter and mass flow rate of vapor. In this approach effects due to the presence of the electric field can be separated from effects caused by the temperature gradients in the thermal boundary layer. The influence of the thermal boundary layer can be investigated after activating the heater at a later stage of the research. For the visualization experiments a test cell was developed. All four vertical walls of the test cell are transparent, and they allow transillumination with laser light for visualization experiments by HI. The bottom electrode is a copper cylinder, which is electrically grounded. The copper block is heated with a resistive heater and it is equipped with 6 thermocouples that provide reference temperatures for the measurements with HI. The top electrode is a mesh electrode. Bubbles are injected with a syringe into the test cell through the bottom electrode. The working fluids presently used in the interferometric visualization experiments, water and PF 5052, satisfy requirements regarding thermophysical, optical and electrical properties. A 30kV power supply equipped with a voltmeter allows to apply the electric field to the electrodes during the experiments. The magnitude of the applied voltage can be adjusted either manually or through the LabVIEW data acquisition and control system connected to a PC. Temperatures of the heated block are recorded using type-T thermocouples, whose output is read by a data acquisition system. Images of the bubbles are recorded with 35mm photographic and 16mm high-speed cameras, scanned and analyzed using various software packages. Visualized temperature fields HI allows the visualization of temperature fields in the vicinity of bubbles during boiling in the form of fringes. Typical visualized temperature distributions around the air bubbles injected into the thermal boundary layer in PF5052 are shown. The temperature of the heated surface is 35 C. The temperature difference for a pair of fringes is approximately 0.05 C. The heat flux applied to the bottom surface is moderate, and the fringe patterns are regular. In the image a bubble penetrating the thermal boundary layer is visible. Because of the axial symmetry of the problem, simplified reconstruction techniques can be applied to recover the temperature field. The thermal plume developing above the heated surface for more intensive heating is shown. The temperature distribution in the liquid is clearly 3D, and tomographic techniques have to be applied to recover the temperature distribution in such a physical situation. A sequence of interferometric images showing the temperature distribution around the rising bubble, recorded with a high-speed camera is shown. Again, the temperature distribution is 3D, and a more complex approach to the evaluation, the tomographic reconstruction has to be taken. Measurement of the temperature distribution from the fringe pattern temperature distributions that yield important information regarding heat transfer are determined. Two algorithms that allow the quantitative evaluation of interferometric fringe patterns and the reconstruction of temperature fields during boiling have been developed at the Heat Transfer Laboratory of the Johns Hopkins University. In the first algorithm the bubble is assumed to be axially symmetrical, which significantly reduces the computational effort for quantifying the temperature distribution around the bubble. For this purpose the thermal boundary layer around the bubble is divided into equidistant concentric shells, and the refractive index is assumed to be constant in each of the shells. Since large temperature gradients are expected in the vicinity of the bubble during boiling, the deflection of the light beam cannot be neglected in boiling experiments. Since the exit angle of the light beam is known, this allows to account for the deflections and phase shifts outside the boundary layer (in the bulk fluid and in the windows of the test cell). Three dimensional temperature distributions in the vicinity of the bubble are reconstructed using tomographic techniques. In tomography, the measurement volume is sliced into 2D planes. In the present study these planes are parallel to the heated surface. The objective is to determine the values of the field parameter of interest in form of the field function in these 2D planes. The field parameter is the change of the refractive index of the liquid in the measurement volume caused by temperature changes. By superimposing data for many 2D planes recorded at the same time instant, the 3D temperature distribution in the measurement volume is recovered.

Herman, Cila↗
NASA NTRS2015

Investigation of Gas Seeding for Planar Laser-Induced Fluorescence in Hypersonic Boundary Layers

Numerical simulations of the gas-seeding strategies required for planar laser-induced fluorescence in a Mach 10 (approximately Mach 8.2 postshock) airflow were performed. The work was performed to understand and quantify the adverse effects associated with gas seeding and to assess various types of seed gas that could potentially be used in future experiments. In prior experiments, NO and NO2 were injected through a slot near the leading edge of a flatplate wedge model used in NASA Langley Research Center's 31 in. Mach 10 air tunnel facility. In this paper, nitric oxide, krypton, and iodine gases were simulated at various injection rates. Simulations showing the deflection of the velocity boundary layer for each of the cases are presented. Streamwise distributions of velocity and concentration boundary-layer thicknesses, as well as vertical distributions of velocity, temperature, and mass distributions, are presented for each of the cases. A comparison between simulated streamwise velocity profiles and experimentally obtained molecular tagging velocimetry profiles using a nitric oxide seeding strategy is performed to verify the influence of such a strategy on the boundary layer. The relative merits of the different seeding strategies are discussed. The results from a custom solver based on OpenFOAM version 2.2.1 are compared against results obtained from ANSYS® Fluent version 6.3.

Arisman, C. J.↗
NASA NTRS1981

Boundary layer thermal stresses in angle-ply composite laminates, part 1

Thermal boundary-layer stresses (near free edges) and displacements were determined by a an eigenfunction expansion technique and the establishment of an appropriate particular solution. Current solutions in the region away from the singular domain (free edge) are found to be excellent agreement with existing approximate numerical results. As the edge is approached, the singular term controls the near field behavior of the boundary layer. Results are presented for cases of various angle-ply graphite/epoxy laminates with (theta/-theta/theta/theta) configurations. These results show high interlaminar (through-the-thickness) stresses. Thermal boundary-layer thicknesses of different composite systems are determined by examining the strain energy density distribution in composites. It is shown that the boundary-layer thickness depends on the degree of anisotropy of each individual lamina, thermomechanical properties of each ply, and the relative thickness of adjacent layers. The interlaminar thermal stresses are compressive with increasing temperature. The corresponding residual stresses are tensile and may enhance interply delaminations.

Wang, S. S.↗
NASA NTRS1986

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.

Gottlieb, D.↗