Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dimensionless Numbers”

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 181 records · Page 10

Effects of a temperature-dependent rheology on large scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗

Effects of a temperature-dependent rheology on large-scale continental extension

The effects of a temperature-dependent rheology on large-scale continental extension are investigated using a thin viscous sheet model. A vertically-averaged rheology is used that is consistent with laboratory experiments on power-law creep of olivine and that depends exponentially on temperature. Results of the calculations depend principally on two parameters: the Peclet number, which describes the relative rates of advection and diffusion of heat, and a dimensionless activation energy, which controls the temperature dependence of the rheology. At short times following the beginning of extension, deformation occurs with negligible change in temperature, so that only small changes in lithospheric strength occur due to attenuation of the lithosphere. However, after a certain critical time interval, thermal diffusion lowers temperatures in the lithosphere, strongly increasing lithospheric strength and slowing the rate of extension. This critical time depends principally on the Peclet number and is short compared with the thermal time constant of the lithosphere. The strength changes cause the locus of high extensional strain rates to shift with time from regions of high strain to regions of low strain. Results of the calculations are compared with observations from the Aegean, where maximum extensional strains are found in the south, near Crete, but maximum present-day strain rates are largest about 300 km further north.

Sonder, Leslie J.↗

Dynamic Formation of a Transient Jet from Arched Magnetized Laboratory Plasma

Abstract A laboratory plasma experiment was built to explore the eruptive behavior of arched magnetized plasmas with dimensionless parameters relevant to the Sun’s photosphere ( β ≈ 10 −3 , Lundquist number ≈10 4 , plasma radius/ion gyroradius ≈20, ion–neutral collision frequency ≫ion cyclotron frequency). Dynamic formation of a transient plasma jet was observed in the presence of the strapping magnetic field. The eruption leading to the jet is unintuitive because the arched plasma is both kink- and torus-stable. The jet structure erupts within a few Alfvén transit times from the formation of the arched plasma. Extensive measurements of plasma temperature, density, magnetic field, and flows are presented. In its early stages, the jet plasma flows away from the arch with supersonic speeds (Mach 1.5). This high-speed flow persists up to the resistive diffusion time in the arched plasma and is driven by large gradients in the magnetic and thermal pressures near the birthplace of jets. There are two distinct electric current channels within the jet, one consisting of outgoing electrons and another composed of electrons returning to the anode footpoint. Significant current density around the jet is a consequence of the diamagnetic current produced by a large thermal pressure gradient in the jet. Ion–neutral charge-exchange collisions provide an efficient mechanism to produce the cross-field current and control the dynamics of the complex current channels of the jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An algebraic model for the turbulent flux of a passive scalar

Direct numerical simulations of the unsteady incompressible Navier-Stokes equations have been performed in order to investigate the behavior of passive-scalar fields resulting from mean scalar gradients in each of three orthogonal directions in homogeneous turbulent shear flow. For all orientations of the mean scalar gradient, the sum of the pressure-scalar gradient and velocity gradient-scalar gradient terms in the turbulent scalar flux balance equation are found to be approximately aligned with the scalar flux vector itself. The simulation results are used to obtain dimensionless model coefficients as a function of the turbulence Reynolds and Peclet numbers.

Rogers, Michael M.↗

Nearly sonic and transsonic convective motions in the solar atmosphere related to the solar wind origin

MHD equations are considered for the solar atmosphere. 15 different simplest MHD regimes are indicated for the momentum transport equation depending on the mutual binary interplay between 6 possible and locally dominant terms: non-stationarity and inhomogeneity of the flow, gas pressure, magnetic tensions, viscous and gravity forces. These regimes are delimited by five physically independent dimensionless parameters, for example, Strouhal, sonic Mach, alfvenic Mach, Reynolds and Froude numbers or their combinations. Another partially overlapping classification of the simplest regimes may be introduced based on the energy conservation equation. There are also 15 independent binary combinations between nonstationary and inhomogeneous convective terms, dissipative energy sinks and sources (viscous, heat-conductive, Joule and radiative ones) in the energy equation. More complicated regimes are considered with multiple dominated terms. All these MHD regimes play their important role somewhere in the solar atmosphere complicated by the tensor transport coefficients in the magnetically dominated regions of the upper atmosphere. Nearly sonic and transsonic nonstationary convective motions with ascending and descending flows are observed in the solar chromosphere. the transition region and the lower pans of the solar corona together with related horizontal velocity components. This convection represents a kind of the 'cocoonery' manufacturing nonstationary vortices generated here and partially connected to the photosphere and to the solar wind. The solar wind originates from this powerful transsonic muddle in the solar atmosphere as a tiny fraction of the streamlines which are temporarily getting detached from the 'cocoons; and going to the infinity. The topologically complicated instantaneous 'runaway surface' around the Sun, i.e., the surface which separates outgoing to the infinity streams from other finite flows in the solar atmosphere was not described in the literature and needs additional investigation. We conclude that a simple one-connected smooth and quasistationary 'critical surface' which is often supposed to be placed somewhere in the solar corona (say, at 4 solar radii) in reality does not exist. Observational and theoretical arguments favor instead of this a highly structured, disordered, patched and time-dependent sonic transition surface in the solar atmosphere permanently fluctuating at positional dependent heights starting sometimes from photospheric (and maybe subphotospheric) levels and more frequently from the chromosphenc network up to several solar radii in the solar corona.

Veselovsky, I. S.↗

Correlating isothermal compressibility to nucleon fluctuations in the inner crust of neutron stars

The question of how and which physical observables or thermodynamic parameters can best predict the onset of a possible phase transition in the inner crust of neutron stars remains largely unresolved. Here, using semiclassical Monte Carlo simulations, we investigate the isothermal compressibility and density fluctuations in a region of relevance to the dynamics of the inner crust. We show that the isothermal compressibility serves as a robust observable to characterize the transition from the nonuniform crust to the uniform core for proton fractions over 0.2. Moreover, we show explicitly how the two-component isothermal compressibility, computed using the Kirkwood-Buff theory, is directly connected to the fluctuations in the number density, recorded in the grand canonical ensemble by monitoring the number of particles in a small volume located at the center of the simulation box. That is, we compute mean-square particle fluctuations and compare them against the isothermal compressibility for different proton fractions. Although our results show that the mean-square particle fluctuations are proportional to the isothermal compressibility, the lack of a perfect correlation is attributed to the relatively small number of particles included in the simulations. The nonunity slope observed in the dimensionless isothermal compressibility—total nucleon fluctuation variance relationship suggests that the inner crust of neutron stars is composed of anisotropic and inhomogeneous matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Flow boiling with enhancement devices for cold plate coolant channel design

Future space exploration and commercialization will require more efficient heat rejection systems. For the required heat transfer rates, such systems must use advanced heat transfer techniques. Forced two phase flow boiling heat transfer with enhancements falls in this category. However, moderate to high quality two phase systems tend to require higher pressure losses. This report is divided into two major parts: (1) Multidimensional wall temperature measurement and heat transfer enhancement for top heated horizontal channels with flow boiling; and (2) Improved analytical heat transfer data reduction for a single side heated coolant channel. Part 1 summarizes over forty experiments which involve both single phase convection and flow boiling in a horizontal channel heated externally from the top side. Part 2 contains parametric dimensionless curves with parameters such as the coolant channel radius ratio, the Biot number, and the circumferential coordinate.

Boyd, Ronald D.↗

Full-scale Force and Pressure-distribution Tests on a Tapered U.S.A. 45 Airfoil

This report presents the results of force and pressure-distribution tests on a 2:1 tapered USA 45 airfoil as determined in the full-scale wind tunnel. The airfoil has a constant-chord center section and rounded tips and is tapered in thickness from 18 percent at the root to 9 percent at the tip. Force tests were made throughout a Reynolds Number range of approximately 2,000,000 to 8,000,000 providing data on the scale effect in addition to the conventional characteristics. Pressure-distribution data were obtained from tests at a Reynolds Number of approximately 4,000,000. The aerodynamic characteristics given by the usual dimensionless coefficients are presented graphically.

Parsons, John F↗

Acoustics Research of Propulsion Systems

The liftoff phase induces some of the highest acoustic loading over a broad frequency for a launch vehicle. These external acoustic environments are used in the prediction of the internal vibration responses of the vehicle and components. Thus, predicting these liftoff acoustic environments is critical to the design requirements of any launch vehicle but there are challenges. Present liftoff vehicle acoustic environment prediction methods utilize stationary data from previously conducted hold-down tests; i.e. static firings conducted in the 1960's, to generate 1/3 octave band Sound Pressure Level (SPL) spectra. These data sets are used to predict the liftoff acoustic environments for launch vehicles. To facilitate the accuracy and quality of acoustic loading, predictions at liftoff for future launch vehicles such as the Space Launch System (SLS), non-stationary flight data from the Ares I-X were processed in PC-Signal in two forms which included a simulated hold-down phase and the entire launch phase. In conjunction, the Prediction of Acoustic Vehicle Environments (PAVE) program was developed in MATLAB to allow for efficient predictions of sound pressure levels (SPLs) as a function of station number along the vehicle using semiempirical methods. This consisted, initially, of generating the Dimensionless Spectrum Function (DSF) and Dimensionless Source Location (DSL) curves from the Ares I-X flight data. These are then used in the MATLAB program to generate the 1/3 octave band SPL spectra. Concluding results show major differences in SPLs between the hold-down test data and the processed Ares IX flight data making the Ares I-X flight data more practical for future vehicle acoustic environment predictions.

Gao, Ximing↗

Acoustics Research of Propulsion Systems

The liftoff phase induces high acoustic loading over a broad frequency range for a launch vehicle. These external acoustic environments are used in the prediction of the internal vibration responses of the vehicle and components. Present liftoff vehicle acoustic environment prediction methods utilize stationary data from previously conducted hold-down tests to generate 1/3 octave band Sound Pressure Level (SPL) spectra. In an effort to update the accuracy and quality of liftoff acoustic loading predictions, non-stationary flight data from the Ares I-X were processed in PC-Signal in two flight phases: simulated hold-down and liftoff. In conjunction, the Prediction of Acoustic Vehicle Environments (PAVE) program was developed in MATLAB to allow for efficient predictions of sound pressure levels (SPLs) as a function of station number along the vehicle using semi-empirical methods. This consisted of generating the Dimensionless Spectrum Function (DSF) and Dimensionless Source Location (DSL) curves from the Ares I-X flight data. These are then used in the MATLAB program to generate the 1/3 octave band SPL spectra. Concluding results show major differences in SPLs between the hold-down test data and the processed Ares I-X flight data making the Ares I-X flight data more practical for future vehicle acoustic environment predictions.

Gao, Ximing↗

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗

On hydrologic similarity: A dimensionless flood frequency model using a generalized geomorphologic unit hydrograph and partial area runoff generation

One of the shortcomings of the original theory of the geomorphologic unit hydrograph (GUH) is that it assumes that runoff is generated uniformly from the entire catchment area. It is now recognized that in many catchments much of the runoff during storm events is produced on partial areas which usually form on narrow bands along the stream network. A storm response model that includes runoff generation on partial areas by both Hortonian and Dunne mechanisms was recently developed by the authors. In this paper a methodology for integrating this partial area runoff generation model with the GUH-based runoff routing model is presented; this leads to a generalized GUH. The generalized GUH and the storm response model are then used to estimate physically based flood frequency distributions. In most previous work the initial moisture state of the catchment had been assumed to be constant for all the storms. In this paper we relax this assumption and allow the initial moisture conditions to vary between storms. The resulting flood frequency distributions are cast in a scaled dimensionless framework where issues such as catchment scale and similarity can be conveniently addressed. A number of experiments are performed to study the sensitivity of the flood frequency response to some of the 'similarity' parameters identified in this formulation. The results indicate that one of the most important components of the derived flood frequency model relates to the specification of processes within the runoff generation model; specifically the inclusion of both saturation excess and Horton infiltration excess runoff production mechanisms. The dominance of these mechanisms over different return periods of the flood frequency distribution can significantly affect the distributional shape and confidence limits about the distribution. Comparisons with observed flood distributions seem to indicate that such mixed runoff production mechanisms influence flood distribution shape. The sensitivity analysis also indicated that the incorporation of basin and rainfall storm scale also greatly influences the distributional shape of the flood frequency curve.

Sivapalan, Murugesu↗

Microgravity Geyser and Flow Field Prediction

Modeling and prediction of flow fields and geyser formation in microgravity cryogenic propellant tanks was investigated. A computational simulation was used to reproduce the test matrix of experimental results performed by other investigators, as well as to model the flows in a larger tank. An underprediction of geyser height by the model led to a sensitivity study to determine if variations in surface tension coefficient, contact angle, or jet pipe turbulence significantly influence the simulations. It was determined that computational geyser height is not sensitive to slight variations in any of these items. An existing empirical correlation based on dimensionless parameters was re-examined in an effort to improve the accuracy of geyser prediction. This resulted in the proposal for a re-formulation of two dimensionless parameters used in the correlation; the non-dimensional geyser height and the Bond number. It was concluded that the new non-dimensional geyser height shows little promise. Although further data will be required to make a definite judgement, the reformulation of the Bond number provided correlations that are more accurate and appear to be more general than the previously established correlation.

Hochstein, J. I.↗

Heat-Transfer and Pressure Measurements on a Flat-Face Cylinder at a Mach Number Range of 2.49 to 4.44

Heat-transfer coefficients and pressure distributions were obtained on a 4-inch-diameter flat-face cylinder in the Langley Unitary Plan wind tunnel. The measured stagnation heat-transfer coefficient agrees well with 55 percent of the theoretical value predicted by the modified Sibulkin method for a hemisphere. Pressure measurements indicated the dimensionless velocity gradient parameter r du\ a(sub t) dx, where x=0 at the stagnation point was approximately 0.3 and invariant throughout the Mach number range from 2.49 to 4.44 and the Reynolds number range from 0.77 x 10(exp 6) to 1.46 x 10(exp 6). The heat-transfer coefficients on the cylindrical afterbody could be predicted with reasonable accuracy by flat-plate theory at an angle of attack of 0 deg. At angles of attack the cylindrical afterbody stagnation-line heat transfer could be computed from swept-cylinder theory for large distances back of the nose when the Reynolds number is based on the distance from the flow reattachment points.

Burbank, Paige B.↗

The study of coronal plasma structures and fluctuations with Faraday rotation measurements

We report dual-frequency, polarimetric measurements of Faraday rotation of extragalactic radio sources viewed through the solar corona. The observations were made at the Very Large Array in 1990 during solar maximum. Of the nine observed, an excess rotation measure of -12.6 rad/sq m was detected for one source (0010+005), which was observed at an elongation of about 9 solar radii. This measurement is in fair agreement with an a priori model rotation measure of -8.6 rad/sq m estimated from coronal potential field models and the electron density model of Paetzold et al. (1992). Our measurement provides a value for the coronal magnetic field strength at 9 solar radii given a knowledge of the magnetic field sector structure, of 12.5 +/- 2.3 mG. Rotation measurements of 0010+005 were made approximately once per hour over an 11 hr period. During this interval, a slow change of about 1 rad/sq m/hr in rotation measure was detected. Although we are not absolutely certain that this drift is not unremoved ionospheric Faraday rotation, extensive analysis of data from the other sources suggests that this is not the case (Sakurai & Spangler 1994). The very long timescale for this variation argues against the agency of magnetohydrodynamics (MHD) waves, and we suggest occultation of 0010+005 by relatively static plasma structures in the corona. We filtered our rotation measure time series to search for variations on an hourly timescale, such as those reported by Hollweg et al. (1992), which could be attributed to coronal MHD waves. We were unable to detect such fluctuations and can report only an upper limit to the rms variation of 1.6 rad/sq m. This upper limit is of the same order, but slightly larger than the values typically reported by Hollweg et al. (1982). This upper limit to the rotation measure fluctuations limits the dimensionless wave amplitude (delta B)/B in the corona to be less than 0.7. Using the number, we estimate the MHD wave flux at the coronal base to be less than 1.6 x 10(exp 5) ergs/sq cm/s. This is less than the amount of wave energy flux required by wave-driven models of the solar wind. Finally, we discuss a number of ways in which such observations could be improved in the future.

Sakurai, Takayuki↗

Measurements in Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Detailed velocity measurements were made along a flat plate subject to the same dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil. Reynolds numbers based on wetted plate length and nominal exit velocity were varied from 50,000 to 300,000, covering cruise to takeoff conditions. Low and high inlet free-stream turbulence intensities (0.2% and 7%) were set using passive grids. The location of boundary-layer separation does not depend strongly on the free-stream turbulence level or Reynolds number, as long as the boundary layer remains non-turbulent prior to separation. Strong acceleration prevents transition on the upstream part of the plate in all cases. Both free-stream turbulence and Reynolds number have strong effects on transition in the adverse pressure gradient region. Under low free-stream turbulence conditions transition is induced by instability waves in the shear layer of the separation bubble. Reattachment generally occurs at the transition start. At Re = 50,000 the separation bubble does not close before the trailing edge of the modeled airfoil. At higher Re, transition moves upstream, and the boundary layer reattaches. With high free-stream turbulence levels, transition appears to occur in a bypass mode, similar to that in attached boundary layers. Transition moves upstream, resulting in shorter separation regions. At Re above 200,000, transition begins before separation. Mean velocity, turbulence and intermittency profiles are presented.

Volino, Ralph J.↗

Extensional Flow of a Polystyrene Boger Fluid Through a 4:1:4 Axisymmetric Contraction/Expansion

The creeping flow of a dilute (0.025 wt%) monodisperse polystyrene/polystyrene Boger fluid through a 4:1:4 axisymmetric contraction/expansion is experimentally observed for a wide range of Deborah numbers. Pressure drop measurements across the orifice plate show a large extra pressure drop that increases monotonically with Deborah number above the value observed for a similar Newtonian fluid at the same flow rate. This enhancement in the dimensionless pressure drop is not associated with the onset of a flow instability, yet it is not predicted by existing steady-state or transient numerical computations with simple dumbbell models. It is conjectured that this extra pressure drop is the result of an additional dissipative contribution to the polymeric stress arising from a stress-conformation hysteresis in the strong non-homogeneous extensional flow near the contraction plane. Such a hysteresis has been independently measured and computed in recent studies of homogeneous transient uniaxial stretching of PS/PS Boger fluids. Flow visualization and velocity field measurements using digital particle image velocimetry (DPIV) show large upstream growth of the corner vortex with increasing Deborah number. At large Deborah numbers, the onset of an elastic instability is observed, first locally as small amplitude fluctuations in the pressure measurements, and then globally as an azimuthal precessing of the upstream corner vortex accompanied by periodic oscillations in the pressure drop across the orifice.

Rothstein, Jonathan P.↗

Transition Experiments on Large Bluntness Cones with Distributed Roughness in Hypersonic Flight

Large bluntness cones with smooth nosetips and roughened frusta were flown in the NASA Ames hypersonic ballistic range at a Mach number of 10 through quiescent air environments. Global surface intensity (temperature) distributions were optically measured and analyzed to determine transition onset and progression over the roughened surface. Real-gas Navier-Stokes calculations of model flowfields, including laminar boundary layer development in these flowfields, were conducted to predict values of key dimensionless parameters used to correlate transition on such configurations in hypersonic flow. For these large bluntness cases, predicted axial distributions of the roughness Reynolds number showed (for each specified freestream pressure) that this parameter was a maximum at the physical beginning of the roughened zone and decreased with increasing run length along the roughened surface. Roughness-induced transition occurred downstream of this maximum roughness Reynolds number location, and progressed upstream towards the beginning of the roughened zone as freestream pressure was systematically increased. Roughness elements encountered at the upstream edge of the roughened frusta thus acted like a finite-extent trip array, consistent with published results concerning the tripping effectiveness of roughness bands placed on otherwise smooth surfaces.

Reda, Daniel. C.↗