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 433 records · Page 24

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part II: Entropy Stability

In the hyperbolic research community, there exists the strong belief that a continuous Galerkin scheme is notoriously unstable and additional stabilization terms have to be added to guarantee stability. In the first part of the series, the application of simultaneous approximation terms for linear problems is investigated where the boundary conditions are imposed weakly. By applying this technique, the authors demonstrate that a pure continuous Galerkin scheme is indeed linearly stable if the boundary conditions are imposed in the correct way. In this work, we extend this investigation to the nonlinear case and focus on entropy conservation. Here, by switching to entropy variables, we provide an estimation of the boundary operators also for nonlinear problems, that guarantee conservation. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗

End-wall boundary layer measurements in a two-stage fan

Detailed flow measurements made in the casing boundary layer of a two-stage transonic fan are summarized. These measurements were taken at a station upstream of the fan, between all blade rows, and downstream of the last row. Conventional boundary layer parameters were calculated from the measured data. A classical two dimensional casing boundary layer was measured at the fan inlet and extended inward to approximately 15 percent of span. A highly three dimensional boundary layer was measured at the exit of each blade row and extended inward to approximately 10 percent of span. The steep radial gradient of axial velocity noted at the exit of the rotors was reduced substantially as the flow passed through the stators. This reduced gradient is attributed to flow mixing. The amount of flow mixing was reflected in the radial redistribution of total temperature as the flow passed through the stators. The blockage factors calculated from the measured data show an increase in blockage across the rotors and a decrease across the stators. For this fan the calculated blockages for the second stage were essentially the same as those for the first stage.

Ball, C. L.↗

End-wall boundary layer measurements in a two-stage fan

Detailed flow measurements made in the casing boundary layer of a two-stage transonic fan are summarized. These measurements were taken at stations upstream of the fan, between all blade rows, and downstream of the last blade row. At the design tip speed of 429 m/sec the fan achieved a peak efficiency of 0.846 at a pressure ratio of 2.471. The boundary layer data were obtained at three weight flows at the design speed: one near choke flow, one near peak efficiency, and one near stall. Conventional boundary layer parameters were calculated from the data measured at each measuring station for each of the three flows. A classical two dimensional casing boundary layer was measured at the fan inlet and extended inward to approximately 15 percent of span. A highly three dimensional boundary layer was measured at the exit of each blade row and extended inward to approximately 10 percent of span. The steep radial gradient of axial velocity noted at the exit of the rotors was reduced substantially as the flow passed through the stators. This reduced gradient is attributed to flow mixing. The amount of flow mixing was reflected in the radial redistribution of total temperature as the flow passed through the stators. The data also show overturning of the tip flow at the stator exits that is consistent with the expected effect of the secondary flow field. The blockage factors calculated from the measured data show an increase in blockage across the rotors and a decrease across the stators.

Ball, C. L.↗

Integral methods in the theory of the boundary layer

The application of the well-known basic principle of mechanics, the principle of Jourdain, to problems of the theory of the boundary layerleads tp an equation from which the equations of Von Karman, Leibenson, and Golubev are derived as special cases. The given equation may be employed in other integral methods. The present paper deals with the method of the variation of the thickness of the boundary layer. A number of new approximate formulas valuable in aerodynamic calculations for the fristion distribution are derived from this procedure. The method has been applied only to laminar boundary layers, but it seems probable that it may be generalized to include turbulent layers as well.

Loitsianskii, L G↗

Computational solution of the defect stream-function equation for nonequilibrium turbulent boundary layers

The derivation of the accurate, second-order, almost linear, approximate equation governing the defect stream function for nonequilibrium compressible turbulent boundary layers is reviewed. The similarity of this equation to the heat conduction equation is exploited in the development of an unconditionally stable, tridiagonal computational method which is second-order accurate in the marching direction and fourth-order accurate in the surface-normal direction. Results compare well with experimental data. Nonlinear effects are shown to be small. This two-dimensional method is simple and has been implemented on a programmable calculator.

Barnwell, Richard W.↗

Lear jet boundary layer/shear layer laser propagation experiments

Optical degradations of aircraft turbulent boundary layers with shear layers generated by aerodynamic fences are analyzed. A collimated 2.5 cm diameter helium-neon laser (0.63 microns) traversed the approximate 5 cm thick natural aircraft boundary layer in double pass via a reflective airfoil. In addition, several flights examined shear layer-induced optical degradation. Flight altitudes ranged from 1.5 to 12 km, while Mach numbers were varied from 0.3 to 0.8. Average line spread function (LSF) and Modulation Transfer Function (MTF) data were obtained by averaging a large number of tilt-removed curves. Fourier transforming the resulting average MTF yields an LSF, thus affording a direct comparison of the two optical measurements. Agreement was good for the aerodynamic fence arrangement, but only fair in the case of a turbulent boundary layer. Values of phase variance inferred from the LSF instrument for a single pass through the random flow and corrected for a large aperture ranged from 0.08 to 0.11 waves (lambda = .63 microns) for the boundary layer. Corresponding values for the fence vary from 0.08 to 0.16 waves. Extrapolation of these values to 10.6 microns suggests negligible degradation for a CO2 laser transmitted through a 5 cm thick, subsonic turbulent boundary layer.

Gilbert, K.↗

Free-air and Bouguer gravity anomalies and the Martian crustal dichotomy

Free-air and Bouguer gravity anomalies from a 50x50 field, derived from re-analysis of Viking Orbiter and Mariner 9 tracking data and using a 50x50 expansion of the current Mars topography and the GSFC degree 50 geoid as the equipotential reference surface, with the Martian crustal dichotomy are compared. The spherical harmonic topography used has zero mean elevation, and differs from the USGS maps by about 2 km. In this field the dichotomy boundary in eastern Mars lies mostly at -1 to -2 km elevation. Bouguer gravity anomalies are shown on a map of Noachian, Hesperian, and Amazonian age terrains, simplified from current geologic maps. The map is centered at 300 deg W to show the continuity of the dichotomy boundary. Contour interval is 100 mgals. Gravity and topography were compared along approximately 40 profiles oriented parallel to the dichotomy boundary topographic gradient, to determine how the geophysical character of the boundary changes along its length and what this implies for its origin and development.

Frey, Herbert↗

Inception length to a fully-developed fin-generated shock wave boundary-layer interaction

An experimental study of fin-generated shock wave turbulent boundary-layer interactions confirmed previous observations that, sufficiently far from the fin apex, such interactions become conical. The inception length to conical symmetry was found to increase weakly with Mach number for Mach numbers from 2.5 to 4 and fin angles from 4 to 22 deg. For the range of interactions examined, the inception length was found to depend primarily upon the inviscid shock angle, this angle ranging from 21 to 40 deg. The behavior of the inception length with shock angle can be broadly divided into two categories. For 'weak' interactions with shock angles less than about 35 deg, the inception length decreased as the shock angle increased. For 'strong' interactions with shock angles greater than about 35 deg, the inception region was small and was approximately constant at three boundary-layer thicknesses in length. In the latter, strong interaction case, the inception length was an order of magnitude smaller than that found in the weakest interactions examined, to the extent that strong interactions were practically fully-developed from the apex.

Lu, Frank K.↗

Ruthenium/Iridium Ratios in the Cretaceous-tertiary Boundary Clay: Implications for Global Dispersal and Fractionation Within the Ejecta Cloud

Ruthenium (Ru) and iridium (Ir) are the least mobile platinum group elements (PGE's) within the Cretaceous-Tertiary (K-T) boundary clay (BC). The Ru/Ir ratio is, therefore, the most useful PGE interelement ratio for distinguishing terrestrial and extraterrestrial contributions to the BC. The Ru/Ir ratio of marine K-T sections (1.77 +/- 0.53) is statistically different from that of the continental sections (0.93 +/- 0.28). The marine Ru/Ir ratios are chondritic (C1 = 1.48 +/- 0.09), but the continental ratios are not. We discovered an inverse correlation of shocked quartz size (or distance from the impact site) and Ru/Ir ratio. This correlation may arise from the difference in Ru and Ir vaporization temperature and/or fractionation during condensation from the ejecta cloud. Postsedimentary alteration, remobilization, or terrestrial PGE input may be responsible for the Ru/Ir ratio variations within the groups of marine and continental sites studied. The marine ratios could also be attained if approximately 15 percent of the boundary metals were contributed by Deccan Trap emissions. However, volcanic emissions could not have been the principal source of the PGE's in the BC because mantle PGE ratios and abundances are inconsistent with those measured in the clay. The Ru/Ir values for pristine Tertiary mantle xenoliths (2.6 +/- 0.48), picrites (4.1 +/- 1.8), and Deccan Trap basalt (3.42 +/- 1.96) are all statistically distinct from those measured in the K-T BC.

Evans, Noreen Joyce↗

Experimental Investigation of the Pressure Rise Required for the Incipient Separation of Turbulent Boundary Layers in Two-Dimensional Supersonic Flow

An experimental investigation has been made of turbulent boundary-layer separation associated with compression corners, curved surfaces of various radii, and incident shock waves. The purpose of the investigation was to provide design information, and to define significant physical trends, which would aid in the prediction of turbulent separation for various aerodynamic devices, such as compressor blades, flaps, spoilers, and diffusers. A characteristic change in the longitudinal static-pressure distribution (i.e., a change from a curve with one inflection point to a curve with three inflection points) was employed to detect the occurrence of separation. The effects of Reynolds number (10(exp 6) to 10(exp 7) per foot or l.5 x 10(exp 4) to 7.5 x 10(exp 4) based upon boundary-layer thickness) and Mach number (1.6 to 4.2) on the onset of turbulent boundary-layer separation were investigated. The pressure gradient of the boundary-layer flow ahead of the interaction region was essentially zero. The results show a considerable effect of Mach number on the pressure rise for incipient separation for all configurations. For a curved-surface model, the static pressure-rise ratio required to cause separation varied from about 2.5, at a Mach number of 2 to about 16, at a Mach number of 3.5. A substantial effect of Reynolds number on the pressure rise for incipient separation was observed in the upper Mach number range and in the lower Reynolds number range; namely, the pressure rise required for separation decreased with increasing Reynolds number. For low Mach numbers and high Reynolds numbers, there appeared to be no Reynolds number effect. The effects of Mach number and of Reynolds number were similar for all models. Model shape was also found to be an important variable affecting the onset of separation. Large gains were realized in the pressure-rise ratio with no separation when the radius of curvature of the model surface was increased. At a Mach number of 3.4, for instance, the pressure-rise ratio with no separation increased from about 5 to 15 as a result of an increase in the radius of curvature from approximately 0 to 30 boundary-layer thicknesses.

Kuehn, Donald M.↗

Subcloud and Cloud-Base Latent Heat Fluxes during Shallow Cumulus Convection

Doppler and Raman lidar observations of vertical velocity and water vapor mixing ratio are used to probe the physics and statistics of subcloud and cloud-base latent heat fluxes during cumulus convection at the ARM Southern Great Plains (SGP) site in Oklahoma, United States. The statistical results show that latent heat fluxes increase with height from the surface up to ~0.8 Z i (where Z i is the convective boundary layer depth) and then decrease to ~0 at Z i . Peak fluxes aloft exceeding 500 W m −2 are associated with periods of increased cumulus cloud cover and stronger jumps in the mean humidity profile. These entrainment fluxes are much larger than the surface fluxes, indicating substantial drying over the 0–0.8 Z i layer accompanied by moistening aloft as the CBL deepens over the diurnal cycle. We also show that the boundary layer humidity budget is approximately closed by computing the flux divergence across the 0–0.8 Z i layer. Composite subcloud velocity and water vapor anomalies show that clouds are linked to coherent updraft and moisture plumes. The moisture anomaly is Gaussian, most pronounced above 0.8 Z i and systematically wider than the velocity anomaly, which has a narrow central updraft flanked by downdrafts. This size and shape disparity results in downdrafts characterized by a high water vapor mixing ratio and thus a broad joint probability density function (JPDF) of velocity and mixing ratio in the upper CBL. We also show that cloud-base latent heat fluxes can be both positive and negative and that the instantaneous positive fluxes can be very large (~10 000 W m −2 ). However, since cloud fraction tends to be small, the net impact of these fluxes remains modest.

54 ENVIRONMENTAL SCIENCES↗

Fast permeability measurement for tight reservoir cores using only initial data of the one chamber pressure pulse decay test

Here, in this study, a mathematical model for fast determination of the permeabilities of tight rocks using measurements taken from the initial period of the One Chamber Pressure Pulse Decay (OC-PPD) test is presented. The model applies to measurements taken both before and after the pressure pulse front has reached the downstream end of the specimen. The analytical solutions for the pressure decay in the upstream chamber are derived based on a parabolic arc approximation of pore pressure distribution along the test specimen. This approximation allows converting the initial–boundary value problem of fluid diffusion in the specimen, governed by partial differential equations, to a system of ordinary differential equations that can be easily solved by explicit formulae. Thus, an explicit formula for the pressure decay rate is obtained, which enables inverse analysis of the initial experimental data to estimate the rock permeability. The proposed method expedites the pulse decay test as it does not require the system to reach equilibrium. The method is validated with three sets of experimental data of the OC-PPD test using helium as the diffusing fluid, for which the relative error of the permeability is found to be less than 6%. This method is particularly useful if the equilibrium time of the pulse decay test for rock specimens with permeabilities in the range of nano-Darcy takes hours or days.

early-time solution↗

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗

Numerical methods for stiff systems of two-point boundary value problems

Numerical procedures are developed for constructing asymptotic solutions of certain nonlinear singularly perturbed vector two-point boundary value problems having boundary layers at one or both endpoints. The asymptotic approximations are generated numerically and can either be used as is or to furnish a general purpose two-point boundary value code with an initial approximation and the nonuniform computational mesh needed for such problems. The procedures are applied to a model problem that has multiple solutions and to problems describing the deformation of thin nonlinear elastic beam that is resting on an elastic foundation.

Flaherty, J. E.↗

Cusp/cleft auroral activity in relation to solar wind dynamic pressure, interplanetary magnetic field B(sub z) and B(sub y)

Continuous optical observations of cusp/cleft auroral activities within approximately equal to 09-15 MLT and 70-76 deg magnetic latitude are studied in relation to changes in solar wind dynamic pressure and interplanetary magnetic field (IMF) variability. The observed latitudinal movements of the cusp/cleft aurora in response to IMF B(sub z) changes may be explained as an effect of a variable magnetic field intensity in the outer dayside magnetosphere associated with the changing intensity of region 1 field-aligned currents and associated closure currents. Ground magnetic signatures related to such currents were observed in the present case (January 10, 1993). Strong, isolated enhancements in solar wind dynamic pressure (Delta p/p is greater than or equal to 0.5) gave rise to equatorward shifts of the cusp/cleft aurora, characteristic auroral transients, and distinct ground magnetic signatures of enhanced convection at cleft latitudes. A sequence of auroral events of approximately equal to 5-10 min recurrence time, moving eastward along the poleward boundary of the persistent cusp/cleft aurora in the approximately equal to 10-14 MLT sector, during negative IMF B(sub z) and B(sub y) conditions, were found to be correlated with brief pulses in solar wind dynamic pressure (0.1 is less than Delta p/p is less than 0.5). Simultaneous photometer observations from Ny Alesund, Svalbard, and Danmarkshavn, Greenland, show that the events often appeared on the prenoon side (approximately equal to 10-12 MLT), before moving into the postnoon sector in the case we study here, when IMF B(sub y) is less than 0. In other cases, similar auroral event sequences have been observed to move westward in the prenoon sector, during intervals of positive B(sub y). Thus a strong prenoon/postnoon asymmetry of event occurence and motion pattern related to the IMF B(sub y) polarity is observed. We find that this category of auroral event sequence is stimulated bursts of electron precipitation that originate from magnetosheath plasma that has accessed that dayside magnetosphere in the noon or near-noon sector, possibly at high latitudes, partly governed by the IMF orientation as well as by solar wind dynamic pressure pulses.

Sandholt, P. E.↗

Factors Influencing Pitot Probe Centerline Displacement in a Turbulent Supersonic Boundary Layer

When a total pressure probe is used for measuring flows with transverse total pressure gradients, a displacement of the effective center of the probe is observed (designated Delta). While this phenomenon is well documented in incompressible flow and supersonic laminar flow, there is insufficient information concerning supersonic turbulent flow. In this study, three NASA Lewis Research Center Supersonic Wind Tunnels (SWT's) were used to investigate pitot probe centerline displacement in supersonic turbulent boundary layers. The relationship between test conditions and pitot probe centerline displacement error was to be determined. For this investigation, ten circular probes with diameter-to-boundary layer ratios (D/delta) ranging from 0.015 to 0.256 were tested in the 10 ft x 10 ft SWT, the 15 cm x 15 cm SWT, and the 1 ft x 1 ft SWT. Reynolds numbers of 4.27 x 10(exp 6)/m, 6.00 x 10(exp 6)/in, 10.33 x 10(exp 6)/in, and 16.9 x 10(exp 6)/m were tested at nominal Mach numbers of 2.0 and 2.5. Boundary layer thicknesses for the three tunnels were approximately 200 mm, 13 mm, and 30 mm, respectively. Initial results indicate that boundary layer thickness, delta, and probe diameter, D/delta play a minimal role in pitot probe centerline offset error, Delta/D. It appears that the Mach gradient, dM/dy, is an important factor, though the exact relationship has not yet been determined. More data is needed to fill the map before a conclusion can be drawn with any certainty. This research provides valuable supersonic, turbulent boundary layer data from three supersonic wind tunnels with three very different boundary layers. It will prove a valuable stepping stone for future research into the factors influencing pitot probe centerline offset error.

Grosser, Wendy I.↗

The Effect of Backward-Facing Step Height on Instability Growth and Breakdown in Swept Wing Boundary-Layer Transition

A low-speed experiment was performed on a swept at plate model with an imposed pressure gradient to determine the effect of a backward-facing step on transition in a stationary-cross flow dominated flow. Detailed hot-wire boundary-layer measurements were performed for three backward-facing step heights of approximately 36, 45, and 49% of the boundary-layer thickness at the step. These step heights correspond to a subcritical, nearly-critical, and critical case. Three leading-edge roughness configurations were tested to determine the effect of stationary-cross flow amplitude on transition. The step caused a local increase in amplitude of the stationary cross flow for the two larger step height cases, but farther downstream the amplitude decreased and remained below the baseline amplitude. The smallest step caused a slight local decrease in amplitude of the primary stationary cross flow mode, but the amplitude collapsed back to the baseline case far downstream of the step. The effect of the step on the amplitude of the primary cross flow mode increased with step height, however, the stationary cross flow amplitudes remained low and thus, stationary cross flow was not solely responsible for transition. Unsteady disturbances were present downstream of the step for all three step heights, and the amplitudes increased with increasing step height. The only exception is that the lower frequency (traveling crossflow-like) disturbance was not present in the lowest step height case. Positive and negative spikes in instantaneous velocity began to occur for the two larger step height cases and then grew in number and amplitude downstream of reattachment, eventually leading to transition. The number and amplitude of spikes varied depending on the step height and cross flow amplitude. Despite the low amplitude of the disturbances in the intermediate step height case, breakdown began to occur intermittently and the flow underwent a long transition region.

Eppink, Jenna L.↗