Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Sidewall boundary-layer measurements with upstream suction in the Langley 0.3-meter transonic cryogenic tunnel

The Langley 0.3 Meter Transonic Cryogenic Tunnel has provision for boundary removal from the sidewalls to reduce sidewall interference effects on the test data. The tests carried out to determine the change in the empty test section sidewall boundary layer thickness at the model station with upstream boundary layer mass removal are described. The boundary layer measurements showed that the upstream removal region is effective in reducing the boundary layer thickness at the model station. The boundary layer displacement thickness reduced from about 1.2 percent to about .4 percent of the test section width. The boundary layer velocity profiles followed a power law variation in the outer region and showed good correlation when plotted in terms of boundary layer momentum thickness.

Murthy, A. V.↗
NASA NTRS1988

Mineralogical and geochemical anomalous data of the K-T boundary samples

Cretaceous-Tertiary boundary problem has been discussed previously from the geological research, mainly by fossil changes. Although geochemical bulk data of Ir anomaly suggest the extraterrestrial origin of the K-T boundary, the exact formation process discussed mainly by mineralogical and geochemical study has been started recently, together with noble gas contents. The K-T boundary sample at Kawaruppu River, Hokkaido was collected, in order to compare with the typical K-T boundary samples of Bubbio, Italy, Stevns Klint, Denmark, and El Kef, Tunisia. The experimental data of the silicas and calcites in these K-T boundary samples were obtained from the X-ray unit-cell dimension (i.e., density), ESR signal and total linear absorption coefficient, as well as He and Ne contents. The K-T boundary samples are usually complex mixture of the terrestrial activities after the K-T boundary event. The mineralogical and geochemical anomalous data indicate special terrestrial atmosphere at the K-T boundary formation probably induced by asteroid impact, followed the many various terrestrial activities (especially the strong role of sea-water mixture, compared with terrestrial highland impact and impact craters in the other earth-type planetary bodies).

Miura, Y.↗
NASA NTRS1989

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

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

Watson, Willie R.↗
NASA NTRS1991

Unified aeroacoustics analysis for high speed turboprop aerodynamics and noise. Volume 5: Propagation of propeller tone noise through a fuselage boundary layer

An analysis of tone noise propagation through a boundary layer and fuselage scattering effects was derived. This analysis is a three dimensional and the complete wave field is solved by matching analytical expressions for the incident and scattered waves in the outer flow to a numerical solution in the boundary layer flow. The outer wave field is constructed analytically from an incident wave appropriate to the source and a scattered wave in the standard Hankel function form. For the incident wave, an existing function - domain propeller noise radiation theory is used. In the boundary layer region, the wave equation is solved by numerical methods. The theoretical analysis is embodied in a computer program which allows the calculation of correction factors for the fuselage scattering and boundary layer refraction effects. The effects are dependent on boundary layer profile, flight speed, and frequency. Corrections can be derived for any point on the fuselage, including those on the opposite side from the source. The theory was verified using limited cases and by comparing calculations with available measurements from JetStar tests of model prop-fans. For the JetStar model scale, the boundary layer refraction effects produce moderate fuselage pressure reinforcements aft of and near the plane of rotation and significant attenuation forward of the plane of rotation at high flight speeds. At lower flight speeds, the calculated boundary layer effects result in moderate amplification over the fuselage area of interest. Apparent amplification forward of the plane of rotation is a result of effective changes in the source directivity due to boundary layer refraction effects. Full scale effects are calculated to be moderate, providing fuselage pressure amplification of about 5 dB at the peak noise location. Evaluation using available noise measurements was made under high-speed, high-altitude flight conditions. Comparisons of calculations made of free field noise, using a current frequency-domain propeller noise prediction method, and fuselage effects using this new procedure show good agreement with fuselage measurements over a wide range of flight speeds and frequencies. Correction factors for the JetStar measurements made on the fuselage are provided in an Appendix.

Magliozzi, B.↗
NASA NTRS1991

A study of methods to investigate nozzle boundary layer transition

To further investigate nozzle flow, numerical computations are employed. The computations produce complete flow velocity and temperature fields within the nozzle. As a check, these results can be compared with experimental data at the wall. Once an accurate numerical scheme has been validated, it can be used as a design tool to predict the performance of other nozzle designs without the cost of experimental testing. Typically, the numerical analysis assumes either a laminar boundary layer or a fully turbulent boundary layer which is steady and two-dimensional. Boundary layer transition is not considered. Computing both the completely laminar boundary layer and the completely turbulent boundary layer conditions gives the minimum and maximum wall heat flux possible for a specified geometry. When the experimental heat flux measurements lie between these two values, the nature of the boundary layer is unknown. The boundary layer may have transitioned from laminar to turbulent; three-dimensional structures may be present in the boundary layer, or the inlet flow conditions may not be correctly specified in the computation.

Pauley, Laura L.↗
NASA NTRS1992

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗
NASA NTRS1992

A fast, uncoupled, compressible, two-dimensional, unsteady boundary layer algorithm with separation for engine inlets

A finite difference boundary layer algorithm was developed to model viscous effects when an inviscid core flow solution is given. This algorithm solved each boundary layer equation separately, then iterated to find a solution. Solving the boundary layer equations sequentially was 2.4 to 4.0 times faster than solving the boundary layer equations simultaneously. This algorithm used a modified Baldwin-Lomax turbulence model, a weighted average of forward and backward differencing of the pressure gradient, and a backward sweep of the pressure. With these modifications, the boundary layer algorithm was able to model flows with and without separation. The number of grid points used in the boundary layer algorithm affected the stability of the algorithm affected the stability of the algorithm as well as the accuracy of the predictions of friction coefficients and momentum thicknesses. Results of this boundary layer algorithm compared well with experimental observations of friction coefficients and momentum thicknesses. In addition, when used interactively with an inviscid flow algorithm, this boundary layer algorithm corrected for viscous effects to give a good match with experimental observations for pressures in a supersonic inlet.

Roach, Robert L.↗
NASA NTRS1994

A non-local computational boundary condition for duct acoustics

A non-local boundary condition is formulated for acoustic waves in ducts without flow. The ducts are two dimensional with constant area, but with variable impedance wall lining. Extension of the formulation to three dimensional and variable area ducts is straightforward in principle, but requires significantly more computation. The boundary condition simulates a nonreflecting wave field in an infinite duct. It is implemented by a constant matrix operator which is applied at the boundary of the computational domain. An efficient computational solution scheme is developed which allows calculations for high frequencies and long duct lengths. This computational solution utilizes the boundary condition to limit the computational space while preserving the radiation boundary condition. The boundary condition is tested for several sources. It is demonstrated that the boundary condition can be applied close to the sound sources, rendering the computational domain small. Computational solutions with the new non-local boundary condition are shown to be consistent with the known solutions for nonreflecting wavefields in an infinite uniform duct.

Zorumski, William E.↗
NASA NTRS1994

Modelling the low-latitude boundary layer with reconnection entry

We develop a one-dimensional Low Latitude Boundary Layer (LLBL) model for northward interplanetary magnetic field (IMF). The boundary layer in this model is uniform in the direction normal to the magnetopause, a 'plateau-type' boundary layer. The boundary layer motion is decoupled from the magnetosheath motion and driven by the plasma pressure associated with the incoming solar wind plasma near local noon, which has become entrained on closed field lines as a result of reconnection in the cusp region. Dissipation in the ionosphere at the feet of the boundary layer field lines opposes this motion. There are two physical solutions for the model. In one, the boundary layer reaches a terminal velocity in the tail as the boundary layer plasma effectively joins the solar wind flow. In the other solution, the flow is nearly stopped in the far tail. In combination with other mechanisms, this latter solution may correspond to the case in which the boundary layer plasma participates in magnetospheric convection and returns sunward. The density, velocity, and thickness as functions of distance from local noon are studied, assuming that the magnetopause hasa elliptical shape and the magnetospheric field is dipolar.

Song, P.↗
NASA NTRS1994

Low-latitude boundary layer near noon: An open field line model

We propose that many features of the cusp and low-latitude boundary layer (LLBL) observed near noon MLT can be explained by interpreting the LLBL as being on open lines with an inner boundary at the separatrix between open and closed magnetic field lines. This interpretation places the poleward boundary of the LLBL and equatorward boundary of the cusp along the field line that bifurcates at the cusp neutral point. The interpretation accounts for the abrupt boundary of magnetosheath particles at the inner edge of the LLBL, a feature that is inconsistent with LLBL formation by diffusion onto closed field lines, and for the distribution of magnetosheath particles appearing more as one continuous region than as two distinct regions across the noon cusp/LLBL boundary. Furthermore, we can explain the existence of energetic radiation belt electrons and protons with differing pitch angle distributions within the LLBL and their abrupt cutoff at the poleward boundary of the LLBL. By modeling the LLBL and cusp region quantitatively, we can account for a hemispherical difference in the location of the equatorial boundary of the cusp that is observed to be dependent on the dipole tilt angle but not on the interplanetary magnetic field (IMF) x component. We also find important variations and hemispherical differences in that the size of the LLBL that should depend strongly upon the x component of the IMF. This prediction is observationally testable. Finally, we find that when the IMF is strongly northward, the LLBL may include a narrow region adjacent to the magnetopause where field lines are detached (i.e., have both ends connected to the IMF).

Lyons, L. R.↗
NASA NTRS1996

Effective Boundary Treatment for the Biharmonic Dirichlet Problem

The biharmonic equation can be rewritten as a system of two Poisson equations. Multigrid solution of this system is expected to converge with the same amount of work as solving two Poisson equations, requiring less than 70 floating point operations (scalar multiply or addition) per fine grid point to reach a solution using an FMG algorithm. For periodic boundary conditions, this goal is attained by simple, straightforward application of multigrid. For Dirichlet boundary conditions, however, convergence is impeded by poor interaction with the boundaries. Attempts to overcome the slowness without specifically addressing the boundaries have resulted in multigrid algorithms not attaining the Poisson convergence rate. We present three methods of boundary treatment with which full multigrid efficiency can be obtained. All implement an approach described by Brandt, concentrating some additional effort near the boundary. The first approach simply adds a number of relaxation sweeps over points close to the boundary. The second uses joint relaxation on near-boundary points. The third method takes something from each of the first two methods, resulting in a solver more suitable for highly parallel applications.

Brandt, A.↗
NASA NTRS1997

External Boundary Conditions for Three-Dimensional Problems of Computational Aerodynamics

We consider an unbounded steady-state flow of viscous fluid over a three-dimensional finite body or configuration of bodies. For the purpose of solving this flow problem numerically, we discretize the governing equations (Navier-Stokes) on a finite-difference grid. The grid obviously cannot stretch from the body up to infinity, because the number of the discrete variables in that case would not be finite. Therefore, prior to the discretization we truncate the original unbounded flow domain by introducing some artificial computational boundary at a finite distance of the body. Typically, the artificial boundary is introduced in a natural way as the external boundary of the domain covered by the grid. The flow problem formulated only on the finite computational domain rather than on the original infinite domain is clearly subdefinite unless some artificial boundary conditions (ABC's) are specified at the external computational boundary. Similarly, the discretized flow problem is subdefinite (i.e., lacks equations with respect to unknowns) unless a special closing procedure is implemented at this artificial boundary. The closing procedure in the discrete case is called the ABC's as well. In this paper, we present an innovative approach to constructing highly accurate ABC's for three-dimensional flow computations. The approach extends our previous technique developed for the two-dimensional case; it employs the finite-difference counterparts to Calderon's pseudodifferential boundary projections calculated in the framework of the difference potentials method (DPM) by Ryaben'kii. The resulting ABC's appear spatially nonlocal but particularly easy to implement along with the existing solvers. The new boundary conditions have been successfully combined with the NASA-developed production code TLNS3D and used for the analysis of wing-shaped configurations in subsonic (including incompressible limit) and transonic flow regimes. As demonstrated by the computational experiments and comparisons with the standard (local) methods, the DPM-based ABC's allow one to greatly reduce the size of the computational domain while still maintaining high accuracy of the numerical solution. Moreover, they may provide for a noticeable increase of the convergence rate of multigrid iterations.

Tsynkov, Semyon V.↗
NASA NTRS1998

Auroral Boundaries: Comparison Between UV Images, In Situ Precipitation, and Groundbased Optical Observations

The location of poleward and equatorward boundaries are used as diagnostics of substorm phase and energy storage in the magnetosphere. Boundaries are estimated from groundbased observations (all sky imagers, meridian scanning photometers, and magnetometers), from in situ particle observations, and from spacebased global auroral observations. Each observational technique has relative advantages. However, recent studies suggest caution when comparing boundaries derived from different methods. in the study, FUV auroral images from the POLAR Ultraviolet Imager are used to estimate auroral boundaries. These boundaries are then compared with similar boundaries derived from DMSP in situ particle precipitation observations and from groundbased meridian scanning photometers. The goal is to see how well the boundaries from these three dissimilar observations correlate, to estimate dependencies on instrumental capabilities, i.e. imager resolution, and to serve as a potential first step to derive an algorithm for finding oval boundaries from groundbased optical data.

Germany, G. A.↗
NASA NTRS1997

A Smoothed Boundary Condition for Reducing Nonphysical Field Effects

In this paper, we examine the problem associated with abruptly mixing boundary conditions in the context of a two-dimensional semiconductor device simulator. Explicitly, this paper addresses the transition between an ohmic-type Dirichlet condition and a passivated Neumann boundary. In the traditional setting, the details or the transition between the two boundary types are not addressed and an abrupt transition is assumed. Subsequently, the calculated observables (most notably the potential) exhibit discontinuous derivatives near the surface at the point where the boundary type switches. This paper proposes an alternative condition which models the progression between the two boundary types through the use of a finite length, smoothed boundary whereby the numerical discontinuities are eliminated. The physical and mathematical basis for this smoothed boundary condition is discussed and examples of the technique's implementation given. It is found that the proposed boundary condition is numerically efficient and can be implemented in pre-existing device simulators with relative ease.

Smith, Arlynn W.↗
NASA NTRS2003

The Disposition of Pt, Pd, Ir, Os, and Ru in Marine Sediments and the K/T Boundary

The marine record of platinum group elements (PGEs) and Os isotopic compositions provides information on different inputs of PGEs into the oceans. Some studies based on a smaller subset of the PGEs suggest that the PGEs may suffer post-depositional mobility during diagenesis. In some K/T boundary clays, Kyte and others showed that the relative abundances of Pt, Pd, Ir, and Os can differ significantly from chondritic, which is the signature expected from fallout of the meteorite impact. In some K/T boundary sections, elevated Ir concentrations are observed as far as 1 meter from the cm-thick boundary clay containing the meteoritic ejecta. The purpose of this study was to characterize Pt, Pd, Ir, Os, and Ru abundances in zones including the K/T boundary. We determined PGE abundances of boundary clays at two hemipelagic sites (Stevns Klint, Denmark and Caravaca, Spain) in which previous studies by Kyte and others showed that the Ir anomaly is confined to within a few cm. We also analyzed two pelagic Pacific sites: a boundary clay from the north Pacific (Hole 465A) characterized by a 0.5 m thick Ir anomaly and a transect across the K/T boundary from the south Pacific (Hole 596) where the Ir anomaly spans 2 m. The Stevns Klint, Caravaca, and north Pacific sites are characterized by abundant marls and limestones in the section, whereas the south Pacific site is dominated by clays. Samples were spiked with isotopic tracers, mixed with a flux, S and Ni, and equilibrated by fusion. PGEs were extracted from the Ni and analyzed on a Finnigan Element ICP-MS. We find that the narrow Caravaca and Stevns Klint boundary clays have relative PGE abundance patterns indistinguishable from chondritic values. The two Pacific sites were found to have nearly identical PGE patterns but have ratios at the peak, which differ from chondritic values as found earlier by Evans et al. The Pacific sites were found to have nearly identical PGE patterns but are extremely depleted in OS (Os/Ir = 0.07-0.15) and slightly enriched in Pd and Pt relative to Ir.

Lee, Cin-Ty↗
NASA NTRS2001

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗
NASA NTRS2006

Shuttle Damage/Repair from the Perspective of Hypersonic Boundary Layer Transition - Experimental Results

An overview is provided of the experimental wind tunnel program conducted at the NASA Langley Research Center Aerothermodynamics Laboratory in support of an agency-wide effort to prepare the Shuttle Orbiter for Return-to-Flight. The effect of an isolated protuberance and an isolated rectangular cavity on hypersonic boundary layer transition onset on the windward surface of the Shuttle Orbiter has been experimentally characterized. These experimental studies were initiated to provide a protuberance and cavity effects database for developing hypersonic transition criteria to support on-orbit disposition of thermal protection system damage or repair. In addition, a synergistic experimental investigation was undertaken to assess the impact of an isolated mass-flow entrainment source (simulating pyrolysis/outgassing from a proposed tile repair material) on boundary layer transition. A brief review of the relevant literature regarding hypersonic boundary layer transition induced from cavities and localized mass addition from ablation is presented. Boundary layer transition results were obtained using 0.0075-scale Orbiter models with simulated tile damage (rectangular cavities) of varying length, width, and depth and simulated tile damage or repair (protuberances) of varying height. Cavity and mass addition effects were assessed at a fixed location (x/L = 0.3) along the model centerline in a region of near zero pressure gradient. Cavity length-to-depth ratio was systematically varied from 2.5 to 17.7 and length-to-width ratio of 1 to 8.5. Cavity depth-to-local boundary layer thickness ranged from 0.5 to 4.8. Protuberances were located at several sites along the centerline and port/starboard attachment lines along the chine and wing leading edge. Protuberance height-to-boundary layer thickness was varied from approximately 0.2 to 1.1. Global heat transfer images and heating distributions of the Orbiter windward surface using phosphor thermography were used to infer the state of the boundary layer (laminar, transitional, or turbulent). Test parametrics include angles-of-attack of 30 deg and 40 deg, sideslip angle of 0 deg, freestream Reynolds numbers from 0.02x106 to 7.3x106 per foot, edge-to-wall temperature ratio from 0.4 to 0.8, and normal shock density ratios of approximately 5.3, 6.0, and 12 in Mach 6 air, Mach 10 air, and Mach 6 CF4, respectively. Testing to simulate the effects of ablation from a proposed tile repair concept indicated that transition was not a concern. The experimental protuberance and cavity databases highlighted in this report were used to formulate boundary layer transition correlations that were an integral part of an analytical process to disposition observed Orbiter TPS damage during STS- 114.

Horvath, Thomas J.↗
NASA NTRS2009

Numerical Study of Outlet Boundary Conditions for Unsteady Turbulent Internal Flows Using the NCC

This paper presents the results of studies on the outlet boundary conditions for turbulent internal flow simulations. Several outlet boundary conditions have been investigated by applying the National Combustion Code (NCC) to the configuration of a LM6000 single injector flame tube. First of all, very large eddy simulations (VLES) have been performed using the partially resolved numerical simulation (PRNS) approach, in which both the nonlinear and linear dynamic subscale models were employed. Secondly, unsteady Reynolds averaged Navier- Stokes (URANS) simulations have also been performed for the same configuration to investigate the effects of different outlet boundary conditions in the context of URANS. Thirdly, the possible role of the initial condition is inspected by using three different initial flow fields for both the PRNS/VLES simulation and the URANS simulation. The same grid is used for all the simulations and the number of mesh element is about 0.5 million. The main purpose of this study is to examine the long-time behavior of the solution as determined by the imposed outlet boundary conditions. For a particular simulation to be considered as successful under the given initial and boundary conditions, the solution must be sustainable in a physically meaningful manner over a sufficiently long period of time. The commonly used outlet boundary condition for steady Reynolds averaged Navier-Stokes (RANS) simulation is a fixed pressure at the outlet with all the other dependent variables being extrapolated from the interior. The results of the present study suggest that this is also workable for the URANS simulation of the LM6000 injector flame tube. However, it does not work for the PRNS/VLES simulation due to the unphysical reflections of the pressure disturbances at the outlet boundary. This undesirable situation can be practically alleviated by applying a simple unsteady convection equation for the pressure disturbances at the outlet boundary. The numerical results presented in this paper suggest that this unsteady convection of pressure disturbances at the outlet works very well for all the unsteady simulations (both PRNS/VLES and URANS) of the LM6000 single injector flame tube.

Liu, Nan-Suey↗