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 775 records · Page 43

An approximate theoretical method for modeling the static thrust performance of non-axisymmetric two-dimensional convergent-divergent nozzles

An analytical/numerical method has been developed to predict the static thrust performance of non-axisymmetric, two-dimensional convergent-divergent exhaust nozzles. Thermodynamic nozzle performance effects due to over- and underexpansion are modeled using one-dimensional compressible flow theory. Boundary layer development and skin friction losses are calculated using an approximate integral momentum method based on the classic karman-Polhausen solution. Angularity effects are included with these two models in a computational Nozzle Performance Analysis Code, NPAC. In four different case studies, results from NPAC are compared to experimental data obtained from subscale nozzle testing to demonstrate the capabilities and limitations of the NPAC method. In several cases, the NPAC prediction matched experimental gross thrust efficiency data to within 0.1 percent at a design NPR, and to within 0.5 percent at off-design conditions.

Hunter, Craig A.↗

Nonlinear evolution of the first mode supersonic oblique waves in compressible boundary layers. Part 1: Heated/cooled walls

The nonlinear stability of an oblique mode propagating in a two-dimensional compressible boundary layer is considered under the long wave-length approximation. The growth rate of the wave is assumed to be small so that the concept of unsteady nonlinear critical layers can be used. It is shown that the spatial/temporal evolution of the mode is governed by a pair of coupled unsteady nonlinear equations for the disturbance vorticity and density. Expressions for the linear growth rate show clearly the effects of wall heating and cooling and in particular how heating destabilizes the boundary layer for these long wavelength inviscid modes at O(1) Mach numbers. A generalized expression for the linear growth rate is obtained and is shown to compare very well for a range of frequencies and wave-angles at moderate Mach numbers with full numerical solutions of the linear stability problem. The numerical solution of the nonlinear unsteady critical layer problem using a novel method based on Fourier decomposition and Chebychev collocation is discussed and some results are presented.

Gajjar, J. S. B.↗

Cross correlation and length scales in turbulent flows near surfaces

Two kinds of length scales are used in turbulent flows; 'functional length scales' such as mixing length, dissipation length L(sub epsilon), etc., and 'flow-field length scales' derived from cross correlations of velocity, pressure, etc. in the flow. Some connection between these scales are derived here. We first consider the cross correlation R(sub vv)(y,y(sub 1)) of the normal components u at two heights y, y(sub 1) above a rigid surface, normalized by the velocity y(sub 1) (greater than y). For shear-free boundary layers it is found theoretically, and in field and numerical experiments that R(sub vv) approximately equals y/y(sub 1). For shear layers it is also found that R(sub vv) approximately equals f(y/y(sub 1)) less than or equal to y,y(sub 1). This function f differs slightly between low Reynolds number numerical simulations and field experiments. The lateral structure defined by R(sub vv)(y,r(sub 3); y(sub 1),0) is also self similar and shows that the eddies centered at about y(sub 1) appear to have constant lateral width a(sub 3) above and below y(sub 1), where a(sub 3, sup +) approximately equals 7+1/(1.4dU(sup +)/dy(sup +)), when normalized on u(sub *) and v, where U is the mean velocity. Results for L(sub epsilon, sup -1) from direct numerical simulation are found to compare well with the formula L(sub epsilon, sup -1) = A(sub B)/y + A(sub S)dU/dy/v, for unidirectional and reversing turbulent boundary layers and channel flow, except near where dU/dy approximately equals 0. The conclusion is that the large-scale eddy structure and length scales in these flows are determined by a combination of shear and blocking, and that the vertical component of turbulence has a self-similar structure in both kinds of boundary layer.

Hunt, J. C. R.↗

The structure of the atmospheric boundary layer in the central equatorial Pacific during January and February of FGGE

The structure of the atmospheric boundary layer for a region between Hawaii and the equator is studied using dropwindsonde data (about 900 soundings) collected in January and February 1979 as part of the First GARP Global Experiment. These data were used to establish the longitudinal and latitudinal variations in the thermodynamic and wind structure of the boundary layer and describe differences in the structure for convectively active and inactive conditions. Low-level inversions (at a pressure of approximately 850 mb) were found to be an important feature of the structure of the boundary layer for undisturbed conditions. Although the frequency of low-level inversions decreased equatorward from Hawaii, the average height of the inversions observed in these regions did not vary significantly in the north-south direction. For convectively undisturbed conditions, low-level inversions are on the average of sufficient strength to suppress deep convection. However, it is shown they could easily be modified by upward vertical motion to give a boundary layer structure similar to that observed in areas where deep convection was observed.

Firestone, James K.↗

Nonlinear Excitation of Inviscid Stationary Vortex in a Boundary-Layer Flow

We examine the excitation of inviscid stationary crossflow instabilities near an isolated surface hump (or indentation) underneath a three-dimensional boundary layer. As the hump height (or indentation depth) is increased from zero, the receptivity process becomes nonlinear even before the stability characteristics of the boundary layer are modified to a significant extent. This behavior contrasts sharply with earlier findings on the excitation of the lower branch Tollmien-Schlichting modes and is attributed to the inviscid nature of the crossflow modes, which leads to a decoupling between the regions of receptivity and stability. As a result of this decoupling, similarity transformations exist that allow the nonlinear receptivity of a general three-dimensional boundary layer to be studied with a set of canonical solutions to the viscous sublayer equations. The parametric study suggests that the receptivity is likely to become nonlinear even before the hump height becomes large enough for flow reversal to occur in the canonical solution. We also find that the receptivity to surface humps increases more rapidly as the hump height increases than is predicted by linear theory. On the other hand, receptivity near surface indentations is generally smaller in comparison with the linear approximation. Extension of the work to crossflow receptivity in compressible boundary layers and to Gortler vortex excitation is also discussed.

Choudhari, Meelan↗

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗

Breakdown of the $v=1$ integer quantum Hall effect in the single particle and collective localization regimes

We report large signal current–voltage characteristics measurements near the $v=1$ integer quantum Hall plateau in a high mobility GaAs/AlGaAs sample that exhibits both single particle and collective localization. According to prior results, the central part of the $v=1$ plateau corresponds to single particle localization and hence an Anderson insulator. In contrast, the two regions in the flanks of this Anderson insulator were associated with collective localization and were identified with the integer quantum Hall Wigner solid. We find that, at low temperatures and under large bias currents, both the Anderson insulator and integer quantum Hall Wigner solid have a sharp onset of conduction above a threshold. Furthermore, this threshold shows an approximately monotonic dependence on the filling factor across the boundary between the two localized phases. Finally, to our surprise, we find that these large signal characteristics do not clearly distinguish between the Anderson insulator and the integer quantum Hall Wigner solid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Complexion dictated thermal resistance with interface density in reactive metal multilayers

Multilayers composed of aluminum (A1) and platinum (Pt) exhibit a non-monotonic trend in thermal resistance with bilayer thickness as measured by time domain thermoreflectance. The thermal resistance initially increases with reduced bilayer thickness only to reach a maximum and then decrease with further shrinking of the multilayer period. These observations are attributed to the evolving impact of an intermixed amorphous complexion approximately 10 nm in thickness, which forms at each boundary between Al- and Pt-rich layers. Scanning transmission electron microscopy combined with energy dispersive X-ray spectroscopy find that the elemental composition of the complexion varies based on bilayer periodicity as does the fraction of the multilayer composed of this interlayer. These variations in complexion mitigate boundary scattering within the multilayers as shown by electronic transport calculations employing density functional theory and non-equilibrium Green's functions on amorphous structures obtained via finite temperature molecular dynamics. The lessening of boundary scattering reduces the total resistance to thermal transport leading to the observed non-monotonic trend thereby highlighting the central role of complexion on thermal transport within reactive metal multilayers.

36 MATERIALS SCIENCE↗

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗

High latitude proton precipitation and light-ion density profiles during the magnetic storm initial phase

Measurements of precipitating protons and light ion densities by experiments on OGO-4 indicate that widespread proton precipitation occurs in predawn hours during the magnetic storm initial phase from the latitude of the high-latitude ion trough, or plasmapause , up to Lambda 75 deg. A softening of the proton spectrum is apparent as the plasmapause is approached. The separation of the low-latitude precipitation boundaries for 7.3 kev and 23.8 kev protons is approximately 1 deg, compared with a 3.6 deg separation which has been computed using the formulas of Gendrin and Eather and Carovillano. Consideration of probable proton drift morphology leads to the conclusion that protons ase injected in predawn hours, with widespread precipitation occurring in the region outside the plasmapause. Protons less energetic than approximately 7 kev drift eastward, while the more energetic protons drift westward, producing the observed dawn-dusk asymmetry for the lower-energy protons.

Burch, J. L.↗

Influence of solar magnetic sector structure on terrestrial atmosphere vorticity

The solar magnetic sector structure has a sizable and reproducible influence on tropospheric and lower stratospheric vorticity. The average vorticity during winter in the Northhern Hemisphere north of 20 deg N latitude reaches a minimum approximately one day after the passing of a sector boundary, and then increases during the following two or three days. The effect is found at all heights within the troposphere, but is not prominent in the stratosphere, except at the lower levels. No single longitudinal interval appears to dominate the effect.

Wilcox, J. M.↗

Influence of solar magnetic sector structure on terrestrial atmospheric vorticity

The solar magnetic sector structure has a sizable and reproducible influence on tropospheric and lower stratospheric vorticity. The average vorticity during winter in the Northern Hemisphere north of 20N latitude reaches a minimum approximately one day after the passing of a sector boundary, and then increases during the following two or three days. The effect is found at all heights within the troposphere, but is not prominent in the stratosphere, except at the lower levels. No single longitudinal interval appears to dominate the effect.

Wilcox, J. M.↗

Airloads near the open port of a one-meter airborne telescope

An infrared telescope is mounted in a cavity in the fuselage of a C-141 airplane. The cavity is sealed off from the rest of the fuselage and can be opened to the outside for viewing. A retractable porous spoiler is installed ahead of the opening to inhibit cavity resonance. The paper presents the mean pressures, pressure fluctuations, spectra, phase, and coherence of the pressures in the cavity and on the fuselage near the cavity, along with pertinent acceleration and torque information for the aircraft structure and telescope. It is shown that a moderate spoiler extension controls resonance in normal flight attitudes and produces a cavity environment with approximately the same turbulence intensity as an attached turbulent boundary layer.

Buell, D. A.↗

The determination of orbits using Picard iteration

The determination of orbits by using Picard iteration is reported. This is a direct extension of the classical method of Picard that has been used in finding approximate solutions of nonlinear differential equations for a variety of problems. The application of the Picard method of successive approximations to the initial value and the two point boundary value problems is given.

Mikkilineni, R. P.↗

On a spherical star system with a collapsed core

Exact solutions of Einstein's equations for a perfect fluid that describe the gravitational collapse of matter around a collapsed core are presented. The collapsed region grows monotonically until an absolute Schwarzschild horizon is formed. It is found that the minimum redshift of light emitted radially from the boundary of the system and received at infinity is approximately 4.62. These systems are expected to describe the very last stages in the gravitational collapse of a globular star cluster with a central collapsed core.

Glass, E. N.↗

Rotating stratified flow over finite isolated topography

Inviscid, steady, stratified rotating flow over a finite, isolated topographic feature is critically analyzed. The formulation is based on approximating the horizontal momentum by the geostrophic momentum. A boundary value problem governs the perturbation pressure field. The solution is an anticyclonic, topographically-bound vortex whose characteristics are independent of the upstream velocity but do depend on stratification, rotation, and the nature of the topography. The vortex is baroclinic in the vicinity of the mountain but barotropic in the far field. The velocity field is a combination of the bound vortex and an upstream velocity interacting with the perturbation pressure field. The effects of stratification, upstream velocity and the nature of the topography are investigated, and fluid trajectories are plotted.

Merkine, L.-O.↗

Crack border stress and displacement equations revisited

It is more or less accepted in fracture mechanics that the elastic stress and displacements very near to the tip of a plane line crack can be approximated with sufficient accuracy, for all geometries and outer boundary loading conditions, by a one-parameter representation, i.e., strictly in terms of the stress intensity factors KI and/or KII. It is shown here that this presumption which appears to be reasonable on face value, quantitatively speaking, is nevertheless unacceptable as a general proposition. The reason lies with the quite arbitrary practice of omitting the second term of the series representation for the stresses, a contribution which is independent of distance from the crack tip. It is not difficult to show by way of specific examples how such omission can lead to error of serious qualitative nature in the prediction of stress and displacement related quantities of interest.

Eftis, J.↗

Noncompressive density enhancements in the solar wind

When the bulk flow speed is nearly constant or falling, high densities are sometimes observed in the solar wind. These densities do not appear to be generated in interplanetary space. It is noted that the magnetic field is not enhanced within these events, and that the proton and/or electron temperatures are low, varying in opposition to the density. About 1/3 of these density events contains interplanetary magnetic field reversals, some of which are noisy and do not qualify as sector boundaries. It is estimated that the average event contains approximately 10 to the 16th g of material and 2.6 x 10 to the 31st ergs, so that aggregated events, when they are common, make a negligible contribution to the total mass and energy budget of the solar wind at 1 AU. It is suggested that there may be an association between density enhancements and solar coronal mass ejection events.

Gosling, J. T.↗