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 271 records · Page 15

Finite-temperature vibrational full configuration interaction

Thermodynamic functions of an ideal molecular gas due to its anharmonic vibrations are evaluated in a wide range of temperature (T) by the vibrational full-configuration-interaction (FCI) method using a quartic force field and a finite number (N) of harmonic-oscillator basis functions along each normal mode. The thermodynamic functions considered are the grand potential (Ω), internal energy (U), and entropy (S). They are compared with those obtained from the Bose–Einstein theory with or without truncation of the harmonic-oscillator basis functions after quantum number N–1. The comparison reveals that the finite-basis-set errors in Ω and U are, respectively, k B Tln(k B T/Nℏω) and k B T per mode in the high-T limit, obscuring anharmonic effects when k B T > ℏω, where ω is the lowest mode frequency. Here, the benchmark data for several low-order perturbation corrections to Ω, U, and S are also obtained as the numerical derivatives of their FCI values with respect to dimensionless perturbation strength, and the domain of T and N in which these data are reliable (for the N → ∞ limits) is discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗

Analysis of laminar flow between stationary and rotating disks with inflow

The laminar flow between a rotating and a stationary disk with inflow was analyzed. Solutions to the dimensionless governing equations are sought by expanding each of the velocity components in powers of inverse radius. The equations to leading order are those for the configuration with no inflow. The subsequent orders yield sets of linear ordinary differential equations. Solutions are obtained for the first two of these subsequent orders. The solutions indicate that inflow tends to increase the magnitude of the azimuthal velocity in the flow between the two disks and to decrease the torque on the rotating disk. For Prandtl number one, an energy integral is obtained which relates the temperature distribution to the velocity distribution for all Reynolds numbers and therefore eliminates the needs for separate solution of the energy equation.

Rohatgi, U.↗

Dimensionless parameters to categorize the failure modes of ductile plate perforation

In this work, two dimensionless parameters were derived organically from impact mechanics equations to remove the need for empirical constants. These two dimensionless parameters effectively collapse the ballistic performance of ductile metal plates along two distinct loci, depending on the projectile geometry and the corresponding failure mode. Projectiles that initiate ductile hole enlargement via cavity expansion fall along a linear curve, while shear-plugging failure modes tend to fall along a cubic curve. Results were demonstrated for various projectiles of different calibers impacting aluminum and Ti-6Al-4V alloys. The main advantage of the approach in this work is to drastically reduce the number of features for input into machine learning algorithms.

36 MATERIALS SCIENCE↗

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗

Dimensionless parameters to categorize the failure modes of ductile plate perforation

In this work, two dimensionless parameters were derived organically from impact mechanics equations to remove the need for empirical constants. These two dimensionless parameters effectively collapse the ballistic performance of ductile metal plates along two distinct loci, depending on the projectile geometry and the corresponding failure mode. Projectiles that initiate ductile hole enlargement via cavity expansion fall along a linear curve, while shear-plugging failure modes tend to fall along a cubic curve. Results were demonstrated for various projectiles of different calibers impacting aluminum and Ti-6Al-4V alloys, and further explored for more complex target combinations. The main advantage of the approach in this work is to drastically reduce the number of features for input into machine learning algorithms.

36 MATERIALS SCIENCE↗

Forced convection heat transfer to air/water vapor mixtures

Heat transfer coefficients were measured using both dry and humid air in the same forced convection cooling scheme and were compared using appropriate nondimensional parameters (Nusselt, Prandtl and Reynolds numbers). A forced convection scheme with a complex flow field, two dimensional arrays of circular jets with crossflow, was utilized with humidity ratios (mass ratio of water vapor to air) up to 0.23. The dynamic viscosity, thermal conductivity and specific heat of air, steam and air/steam mixtures are examined. Methods for determining gaseous mixture properties from the properties of their pure components are reviewed as well as methods for determining these properties with good confidence. The need for more experimentally determined property data for humid air is discussed. It is concluded that dimensionless forms of forced convection heat transfer data and empirical correlations based on measurements with dry air may be applied to conditions involving humid air with the same confidence as for the dry air case itself, provided that the thermophysical properties of the humid air mixtures are known with the same confidence as their dry air counterparts.

Richards, D. R.↗

A catalogue of potentially bright close binary gravitational wave sources

This is a current print-out of results of a survey, undertaken in the spring of 1985, to identify those known binary stars which might produce significant gravitational wave amplitudes at earth, either dimensionless strain amplitudes exceeding a threshold h = 10(exp -21), or energy fluxes exceeding F = 10(exp -12) erg cm(exp -2) s(exp -1). All real or putative binaries brighter than a certain limiting magnitude (calculated as a function of primary spectral type, orbital period, orbital eccentricity, and bandpass) are included. All double degenerate binaries and Wolf-Rayet binaries with known or suspected orbital periods have also been included. The catalog consists of two parts: a listing of objects in ascending order of Right Ascension (Equinox B1950), followed by an index, listing of objects by identification number according to all major stellar catalogs. The object listing is a print-out of the spreadsheets on which the catalog is currently maintained. It should be noted that the use of this spreadsheet program imposes some limitations on the display of entries. Text entries which exceed the cell size may appear in truncated form, or may run into adjacent columns. Greek characters are not available; they are represented here by the first two or three letters of their Roman names, the first letter appearing as a capital or lower-case letter according to whether the capital or lower-case Greek character is represented. Neither superscripts nor subscripts are available; they appear here in normal position and type-face. The index provides the Right Ascension and Declination of objects sorted by catalogue number.

Webbink, Ronald F.↗

Steady simple shear flow past a circular cylinder at moderate Reynolds numbers - A numerical solution

The two-dimensional steady flow of an incompressible viscous fluid past a circular cylinder, placed symmetrically in a simple shear field, has been studied for both the stationary and the freely rotating case by solving numerically the Navier-Stokes equations for values of the Reynolds number R in the range from 0.047 to 70. At R = 0.047, the results obtained are in substantial agreement with the analytic small-R perturbation solution given by Robertson and Acrivos (1970). Inertia effects were found, however, to play a significant role even at R = 1, and hence the calculated flow pattern for R greater than or equal to 1 differs significantly from that of the creeping-flow solution. Specifically, for the freely rotating case, the region of closed streamlines decreases rapidly in extent with increasing R, two symmetrically placed wakes are formed on either side of the cylinder, and the dimensionless rotational speed of the freely suspended cylinder decreases as the reciprocal of the square root of R.

Kossack, C. A.↗

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heat transfer in rotating serpentine passages with smooth walls

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, smooth-wall heat transfer model with both radially inward and outward flow. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages (coolant-to-wall temperature ratio, Rossby number, Reynolds number and radius-to-passage hydraulic diameter ratio). These four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. It was found that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs and that the effect of rotation on the heat transfer coefficients was markedly different depending on the flow direction. Local heat transfer coefficients were found to decrease by as much as 60 percent and increase by 250 percent from no rotation levels. Comparisons with a pioneering stationary vertical tube buoyancy experiment showed reasonably good agreement. Correlation of the data is achieved employing dimensionless parameters derived from the governing flow equations.

Wagner, J. H.↗

Low Reynolds number multiple-time-scale turbulence model and calculations of steady and pulsating shear layers

A low Reynolds number multiple-time-scale turbulence model (LMS) and its application to fully developed turbulent channel flows and pulsating pipe flows are presented. The LMS can describe the inequilibrium turbulence phenomena down to the viscous sublayer. The calculated fluid flow and turbulence fields for the channel flows are in better agreement with the direct numerical simulation (DNS) results than those obtained using a Reynolds stress turbulence model, and the calculated near-wall dissipation rates are in qualitatively correct agreement with the DNS results. The LMS also successfully predicts the rapidly varying phase-lead of the wall shearing stress that occurs in a narrow range of the dimensionless frequency omega (+) = (omega x nu)/(upsilon x tau(exp2)) for the pulsating pipe flows while various other turbulence models fail to predict this phenomenon, and the LMS yields significantly improved numerical results for a wide range of the dimensionless frequency compared with those obtained using a rapid distortion theory (RDT).

Kim, Sang-Wook↗

A hill-sliding strategy for initialization of Gaussian clusters in the multidimensional space

A hill sliding technique was devised to extract Gaussian clusters from the multivariate probability density estimate of sample data for the first step of iterative unsupervised classification. Each cluster was assumed to posses a unimodal normal distribution. A clustering function proposed distinguished elements of a cluster under formation from the rest in the feature space. Initial clusters were extracted one by one according to the hill sliding tactics. A dimensionless cluster compactness parameter was proposed as a universal measure of cluster goodness and used satisfactorily in test runs with LANDSAT multispectral scanner data. The normalized divergence, defined by the cluster divergence divided by the entropy of the entire sample data, was utilized as a general separability measure between clusters. An overall clustering objective function was set forth in terms of cluster covariance matrices, from which the cluster compactness measure could be deduced. Minimal improvement of initial data partitioning was evaluated by this objective function in eliminating scattered sparse data points. The hill sliding clustering technique developed herein has the potential applicability to decomposition any multivariate mixture distribution into a number of unimodal distributions when an appropriate distribution function to the data set is employed.

Park, J. K.↗

An energy model of droplet impingement on an inclined wall under isothermal and non-isothermal environments

The observation of spray-wall interaction is of great importance to understand the dynamics that occur during fuel impingement onto the chamber wall or piston surfaces in internal combustion engines. It is found that the maximum spreading length of an impinged droplet can provide a quantitative estimation of heat transfer and energy transformation for spray-wall interaction. Furthermore, it influences the air-fuel mixing and hydrocarbon and particle emissions at combusting conditions. In this paper, an analytical model of different droplet-wall impingement conditions is developed in terms of βm (dimensionless maximum spreading length, the ratio of maximum spreading length to initial droplet diameter) to understand the detailed impinging dynamic process. These conditions are grouped as: a single diesel droplet impinging on the wall with different inclined angles (α); cold wall - heated droplet and heated wall - cold droplet impingement when inclined angle of the wall is 0°, respectively. The analytical model is built up based on the energy conservation that considers kinetic energy, gravitation energy, and surface energy before impingement, as well as viscous dissipation, gravitation energy, adhesion energy, deformation energy, and heat energy after impingement. The experimental work of diesel droplet impinging on an inclined wall is performed at a certain range of the Weber number (We of 33 to 420) with various inclined angles (α of 0° to 45°), while for inclined angle is 0°, droplet and wall temperature are varied from 25°C to 150°C to study the effects of the inclined angle and temperature on the temporal evolution of the post-impingement characteristics (i.e. droplet spreading length, dynamic contact angle). The analytical model is validated and evaluated at the aforementioned experimental operating points. The validated model can be employed to predict maximum spreading length of the droplet impinged on the wall. It is further utilized to determine the transition from capillary regime to kinetic regime, then to viscous regime at different inclined angle of the wall.

42 ENGINEERING↗

Coherent-phase or random-phase acceleration of electron beams in solar flares

Time structures of electron beam signatures at radio wavelengths are investigated to probe correlated versus random behavior in solar flares. In particular we address the issue whether acceleration and injection of electron beams is coherently modulated by a single source, or whether the injection is driven by a stochastic (possibly spatially fragmented) process. We analyze a total of approximately = 6000 type III bursts observed by Ikarus (Zurich) in the frequency range of 100-500 MHz, during 359 solar flares with simultaneous greater than or = 25 keV hard X-ray emission, in the years 1890-1983. In 155 flares we find a total of 260 continuous type III groups, with an average number of 13 +/- 9 bursts per group, a mean duration of D = 12 +/- 14 s, a mean period of P = 2.0 +/- 1.2 s, with the highest burst rate at a frequency of nu = 310 +/- 120 MHz. Pulse periods have been measured between 0.5 and 10 s, and can be described by an exponential distribution, i.e., N(P) varies as e (exp -P/1.0s). The period shows a frequency dependence of P(nu)=46(exp-0.6)(sub MHz)s for different flares, but is invariant during a particular flare. We measure the mean period P and its standard deviation sigma (sub p) in each type III group, and quantify the degree of periodicity (or phase-coherence) by the dimensionless parameter sigma (sub p)P. The representative sample of 260 type III burst groups shows a mean periodicity of sigma (sub p/P) = 0.37 +/- 0.12, while Monte Carlo simulations of an equivalent set of truly random time series show a distinctly different value of sigma (sub p)P = 0.93 +/- 0.26. This result indicates that the injection of electron beams is coherently modulated by a particle acceleration source which is either compact or has a global organization on a timescale of seconds, in contrast to an incoherent acceleration source, which is stochastic either in time or space. We discuss the constraints on the size of the acceleration region resulting from electron beam propagation delays and from Alfvenic synchronization during a pulse period. We discuss two periodic preocesses in flares, which potentially control quasi-periodic particle acceleration: (1) MHD oscillations, and (2) current sheets with oscillatory dynamics.

Aschwanden, Markus J.↗

Mixing of Multiple Jets with a Confined Subsonic Crossflow in a Cylindrical Duct

This paper summarizes NASA-supported experimental and computational results on the mixing of a row of jets with a confined subsonic crossflow in a cylindrical duct. The studies from which these results were derived investigated flow and geometric variations typical of the complex 3-D flowfield in the combustion chambers in gas turbine engines. The principal observations were that the momentum-flux ratio and the number of orifices were significant variables. Jet penetration was critical, and jet penetration decreased as either the number of orifices increased or the momentum-flux ratio decreased. It also appeared that jet penetration remained similar with variations in orifice size, shape, spacing, and momentum-flux ratio when the number of orifices was proportional to the square-root of the momentum-flux ratio. In the cylindrical geometry, planar variances are very sensitive to events in the near wall region, so planar averages must be considered in context with the distributions. The mass-flow ratios and orifices investigated were often very large (mass-flow ratio greater than 1 and ratio of orifice area-to-mainstream cross-sectional area up to 0.5), and the axial planes of interest were sometimes near the orifice trailing edge. Three-dimensional flow was a key part of efficient mixing and was observed for all configurations. The results shown also seem to indicate that non-reacting dimensionless scalar profiles can emulate the reacting flow equivalence ratio distribution reasonably well. The results cited suggest that further study may not necessarily lead to a universal 'rule of thumb' for mixer design for lowest emissions, because optimization will likely require an assessment for a specific application.

Holdeman, James D.↗

Pre-Test Assessment of the Use Envelope of the Normal Force of a Wind Tunnel Strain-Gage Balance

The relationship between the aerodynamic lift force generated by a wind tunnel model, the model weight, and the measured normal force of a strain-gage balance is investigated to better understand the expected use envelope of the normal force during a wind tunnel test. First, the fundamental relationship between normal force, model weight, lift curve slope, model reference area, dynamic pressure, and angle of attack is derived. Then, based on this fundamental relationship, the use envelope of a balance is examined for four typical wind tunnel test cases. The first case looks at the use envelope of the normal force during the test of a light wind tunnel model at high subsonic Mach numbers. The second case examines the use envelope of the normal force during the test of a heavy wind tunnel model in an atmospheric low-speed facility. The third case reviews the use envelope of the normal force during the test of a floor-mounted semi-span model. The fourth case discusses the normal force characteristics during the test of a rotated full-span model. The wind tunnel model's lift-to-weight ratio is introduced as a new parameter that may be used for a quick pre-test assessment of the use envelope of the normal force of a balance. The parameter is derived as a function of the lift coefficient, the dimensionless dynamic pressure, and the dimensionless model weight. Lower and upper bounds of the use envelope of a balance are defined using the model's lift-to-weight ratio. Finally, data from a pressurized wind tunnel is used to illustrate both application and interpretation of the model's lift-to-weight ratio.

strain-gage balance↗

Currents and Flows in Distant Magnetospheres

Space scientists have explored, described, and explained the terrestrial magnetosphere for four decades. Rarely do they point out that the planetary and solar wind parameters controlling the size, shape, and activity of Earth's magnetosphere map out only a small portion of the space of dimensionless parameters that govern magnetospheric properties. With the discovery of Ganymede's magnetosphere, the range of parameters relevant to magnetospheric studies has grown by orders of magnitude. Consider the extremes of Ganymede's and Jupiter's magnetospheres. Jupiter's magnetosphere forms within a plasma flowing at super-Alfvenic speed, whereas Ganymede's forms in a sub-Alfvenic flow. The scale sizes of these magnetospheres, characterized by distances to the magnetopause of order 7x10(exp 6) km and 5x10(exp 3) km, respectively, differ by three orders of magnitude, ranging from 100 to 0.1 times the scale of Earth's magnetosphere. The current systems that control the structure and dynamics of a magnetosphere depend on specific plasma and field properties. Magnetopause currents at Ganymede differ greatly from the forms familiar for Earth and Jupiter, principally because the Mach number of the ambient plasma flow greatly influences the shape of the magnetosphere. A magnetodisk current, present at Jupiter because of its rapid rotation, is absent at Earth and Ganymede. The ring current, extensively investigated at Earth, is probably unimportant at Ganymede because the dynamical variations of the external flow are slow. The ring current is subsumed within the magnetodisk current at Jupiter. This paper describes and contrasts aspects of these and other current systems for the three bodies.

Kivelson, Margaret Galland↗

Drop Interactions with the Conical Shock Structure Generated by a Mach 4.5 Projectile

This work presents measurements of liquid drop deformation and breakup time behind approximately conical shock waves and evaluates the predictive capabilities of low-order models and correlations developed using planar shock experiments. A conical shock was approximated by firing a bullet at Mach 4.5 past a vertical column of water drops with a mean initial diameter of $192$ $\mu$m. The time-resolved drop position and maximum transverse dimension were characterized using backlit stereo images taken at 500 kHz. The gas density and velocity fields experienced by the drops were estimated using a Reynolds-averaged Navier–Stokes simulation of the bullet. Classical correlations predict drop breakup times and deformation in error by a factor of 3 or more. The Taylor analogy breakup (TAB) model predicts deformed drop diameters that agree within the confidence bounds of the ensemble-averaged experimental values using a dimensionless constant $C$ 2 = $2$ compared to the accepted value $C$ 2 = $2/3$. In conclusion, results demonstrate existing correlations are inadequate for predicting the drop response to the three-dimensional relaxation of the flowfield downstream of a conical-like shock and suggest the TAB model results represent a path toward improved predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗