Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PRANDTL NUMBER”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Statistical simulation of the magnetorotational dynamo

We analyze turbulence and dynamo induced by the magnetorotational instability (MRI) using quasi-linear statistical simulation methods. We find that homogenous turbulence is unstable to a large scale dynamo instability, which saturates to an inhomogenous equilibrium with a very strong dependence on the magnetic Prandtl number (Pm). Despite its enormously reduced nonlinearity, the quasi-linear model exhibits the same qualitative scaling of angular momentum transport with Pm as fully nonlinear turbulence. This demonstrates the relationship of recent convergence problems to the large scale dynamo and suggests possible methods for studying astrophysically relevant regimes at very low or high Pm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

GRAND MINIMA AND EQUATORWARD PROPAGATION IN A CYCLING STELLAR CONVECTIVE DYNAMO

The 3D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo-generated magnetic fields possesses many timescales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulation’s relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the α-Ω type. The timescales that appear to be relevant to the magnetic polarity reversal are also identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Heat-Exchange in One-Phase Liquid-Metal Flows: Databank (Temperature Distribution in Circular Pipes)

In order to develop accurate methods of calculating heat-exchange in forced channel flows, it is necessary to have knowledge about the temperature distribution over the cross section in differ layers of the flow. For ordinary liquids (gas, water, where the Prandtl number Pr ~ 1), such measurements are difficult to perform because of the small thickness of the boundary layer. The problem simplifies for liquid metals, where Pr << 1. This article contains the format of the data as well as an analysis of the primary data from 37 home and foreign works (mercury, sodium-potassium alloy, sodium, lithium, lead-bismuth alloy). The data are presented in electronic form.

42 ENGINEERING↗
DOE OSTI2018

Morphological Simulation of Phase Separation Coupled Oscillation Shear and Varying Temperature Fields

This paper explores the effect of the shear frequency and Prandtl number (Pr) on the procedure and pattern formation of phase separation in symmetric and asymmetric systems. For the symmetric system, the periodic shear significantly prolongs the spinodal decomposition stage and enlarges the separated domain in domain growth stage. By adjusting the Pr and shear frequency, the number and orientation of separated steady layer structures can be controlled during domain stretch stage. The numerical results indicate that the increase in Pr and decrease in the shear frequency can significantly increase in the layer number of the lamellar structure, which relates to the decrease in domain size. Furthermore, the lamellar orientation parallel to the shear direction is altered into that perpendicular to the shear direction by further increasing the shear frequency, and also similar results for larger systems. For asymmetric system, the quantitative analysis shows that the decrease in the shear frequency enlarges the size of separated minority phases. These numerical results provide guidance for setting the optimum condition for the phase separation under periodic shear and slow cooling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Equatorial symmetry of Boussinesq convective solutions in a rotating spherical shell allowing rotation of the inner and outer spheres

We investigate properties of convective solutions of the Boussinesq thermal convection in a moderately rotating spherical shell allowing the respective rotation of the inner and outer spheres due to the viscous torque of the fluid. The ratio of the inner and outer radii of the spheres, the Prandtl number, and the Taylor number are fixed to 0.4, 1, and 500{sup 2}, respectively. The Rayleigh number is varied from 2.6 × 10{sup 4} to 3.4 × 10{sup 4}. In this parameter range, the behaviours of obtained asymptotic convective solutions are almost similar to those in the system whose inner and outer spheres are restricted to rotate with the same constant angular velocity, although the difference is found in the transition process to chaotic solutions. The convective solution changes from an equatorially symmetric quasi-periodic one to an equatorially symmetric chaotic one, and further to an equatorially asymmetric chaotic one, as the Rayleigh number is increased. This is in contrast to the transition in the system whose inner and outer spheres are assumed to rotate with the same constant angular velocity, where the convective solution changes from an equatorially symmetric quasi-periodic one, to an equatorially asymmetric quasi-periodic one, and to equatorially asymmetric chaotic one. The inner sphere rotates in the retrograde direction on average in the parameter range; however, it sometimes undergoes the prograde rotation when the convective solution becomes chaotic.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Power-law wrinkling turbulence-flame interaction model for astrophysical flames

We extend a model for turbulence-flame interactions (TFI) to consider astrophysical flames with a particular focus on combustion in Type Ia supernovae. The inertial range of the turbulent cascade is nearly always under-resolved in simulations of astrophysical flows, requiring the use of a model in order to quantify the effects of subgrid-scale wrinkling of the flame surface. We provide implementation details to extend a well-tested TFI model to low-Prandtl number flames for use in the compressible hydrodynamics code FLASH. A local, instantaneous measure of the turbulent velocity is calibrated for FLASH and verification tests are performed. Particular care is taken to consider the relation between the subgrid rms turbulent velocity and the turbulent flame speed, especially for high-intensity turbulence where the turbulent flame speed is not expected to scale with the turbulent velocity. Finally, we explore the impact of different TFI models in full-star, three-dimensional simulations of Type Ia supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

MHD effects and heat transfer for the UCM fluid along with Joule heating and thermal radiation using Cattaneo-Christov heat flux model

Present study examines the numerical analysis of MHD flow of Maxwell fluid with thermal radiation and Joule heating by considering the recently developed Cattaneo-Christov heat flux model which explains the time relaxation characteristics for the heat flux. The objective is to analyze the governing parameters such as viscoelastic fluid parameter, Magnetic parameter, Eckert and Prandtl number’s impact on the velocity and temperature profiles through graphs and tables. Suitable similarity transformations have been used to reduce the formulated PDEs into a system of coupled non-linear ODEs. Shooting technique has been invoked for finding the numerical solutions of the dimensionless velocity and temperature profiles. Additionally, the MATLAB built-in routine bvp4c has also been used to verify and strengthen the results obtained by shooting method. From some special cases of the present work, a comparison with the previously published results has been presented.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Dissipation range turbulent cascades in plasmas

Dissipation range cascades in plasma turbulence are described and spectra are formulated from the scaled attenuation in wavenumber space of the spectral energy transfer rate. This yields spectra characterized by the product of a power law and exponential fall-off, applicable to all scales. Spectral indices of the power law and exponential fall-off depend on the scaling of the dissipation, the strength of the nonlinearity, and nonlocal effects when dissipation rates of multiple fluctuation fields are different. The theory is used to derive spectra for MHD turbulence with magnetic Prandtl number greater than unity, extending previous work. The theory is also applied to generic plasma turbulence by considering the spectrum from damping with arbitrary wavenumber scaling. The latter is relevant to ion temperature gradient turbulence modeled by gyrokinetics. The spectrum in this case has an exponential component that becomes weaker at small scale, giving a power law asymptotically. Results from the theory are compared to three very different types of turbulence. These include the magnetic plasma turbulence of the Madison Symmetric Torus, the MHD turbulence of liquid metal in the Madison Dynamo Experiment, and gyrokinetic simulation of ion temperature gradient turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Episodic jet power extracted from a spinning black hole surrounded by a neutrino-dominated accretion flow in gamma-ray bursts

It was suggested that the relativistic jets in gamma-ray bursts (GRBs) are powered via the Blandford-Znajek (BZ) mechanism or the annihilation of neutrinos and anti-neutrinos from a neutrino cooling-dominated accretion flow (NDAF). The advection and diffusion of the large-scale magnetic field of an NDAF is calculated, and the external magnetic field is found to be dragged inward efficiently by the accretion flow for a typical magnetic Prandtl number P{sub m}=η/ν∼1. The maximal BZ jet power can be ∼10{sup 53}-10{sup 54} erg s{sup –1} for an extreme Kerr black hole, if an external magnetic field with 10{sup 14} Gauss is advected by the NDAF. This is roughly consistent with the field strength of the disk formed after a tidal disrupted magnetar. The accretion flow near the black hole horizon is arrested by the magnetic field if the accretion rate is below than a critical value for a given external field. The arrested accretion flow fails to drag the field inward and the field strength decays, and then the accretion re-starts, which leads to oscillating accretion. The typical timescale of such episodic accretion is of an order of one second. This can qualitatively explain the observed oscillation in the soft extended emission of short-type GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE LARGE-SCALE MAGNETIC FIELDS OF THIN ACCRETION DISKS

Large-scale magnetic field threading an accretion disk is a key ingredient in the jet formation model. The most attractive scenario for the origin of such a large-scale field is the advection of the field by the gas in the accretion disk from the interstellar medium or a companion star. However, it is realized that outward diffusion of the accreted field is fast compared with the inward accretion velocity in a geometrically thin accretion disk if the value of the Prandtl number P{sub m} is around unity. In this work, we revisit this problem considering the angular momentum of the disk to be removed predominantly by the magnetically driven outflows. The radial velocity of the disk is significantly increased due to the presence of the outflows. Using a simplified model for the vertical disk structure, we find that even moderately weak fields can cause sufficient angular momentum loss via a magnetic wind to balance outward diffusion. There are two equilibrium points, one at low field strengths corresponding to a plasma-beta at the midplane of order several hundred, and one for strong accreted fields, {beta} {approx} 1. We surmise that the first is relevant for the accretion of weak, possibly external, fields through the outer parts of the disk, while the latter one could explain the tendency, observed in full three-dimensional numerical simulations, of strong flux bundles at the centers of disk to stay confined in spite of strong magnetororational instability turbulence surrounding them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Aspects of the theory of incompressible MHD turbulence with cross-helicity and applications to the solar wind

Solar wind observations have shown that the normalized cross-helicity {sigma}{sub c}, the ratio of the cross-helicity spectrum to the energy spectrum, is approximately constant, independent of wavenumber, throughout the inertial range. This means that the correlation between velocity and magnetic field fluctuations is the same at every scale, that the ratio of the two Elsasser energies (w{sup +}/w{sup -}){sup 2} is the same at every scale, and that the ratio of the energy cascade times of the two Elsasser energies {tau}{sup +}/{tau}{sup -} is the same at every scale. In the case when the magnetic Prandtl number is unity, it can be shown from the equations of incompressible MHD that if {sigma}{sub c} is a constant, then the cascade times of the two Elsasser energies are equal so that {tau}{sup +}/{tau}{sup -} = 1. This is an important constraint for turbulence theories. Using this result, the Goldreich and Sridhar theory and the Boldyrev theory are generalized to MHD turbulence with nonvanishing cross-helicity in such a way that the scaling laws of the original two theories are unchanged. The derivation and some of the important properties of these more general theories shall be presented. Solar wind measurements in support of these theoretical models will also be discussed. For example, new solar wind measurements of the total energy spectrum (kinetic plus magnetic) show that the power-law exponent is closer to 3/2 than 5/3, consistent with simulations of 3D incompressible MHD turbulence with a strong mean meagnetic field that show a 3/2 scaling. For highly Alfvenic, high cross-helicity solar wind turbulence, new measurements presented her show that the average spectral index is 1.540 {+-} 0.033.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

ANALYTICAL MODELS OF EXOPLANETARY ATMOSPHERES. I. ATMOSPHERIC DYNAMICS VIA THE SHALLOW WATER SYSTEM

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical, and spherical), rotation, magnetic tension, and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag, and magnetic drag), and magnetic tension are included. The global atmospheric structure is largely controlled by a single key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations, and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to develop decisively physical intuition and as a reference point for three-dimensional magnetohydrodynamic simulations of atmospheric circulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Turbulent dynamo in a conducting fluid and a partially ionized gas

By following the Kazantsev theory and taking into account both microscopic and turbulent diffusion of magnetic fields, we develop a unified treatment of the kinematic and nonlinear stages of a turbulent dynamo process, and we study the dynamo process for a full range of magnetic Prandtl number P {sub m} and ionization fractions. We find a striking similarity between the dependence of dynamo behavior on P {sub m} in a conducting fluid and R (a function of ionization fraction) in a partially ionized gas. In a weakly ionized medium, the kinematic stage is largely extended, including not only exponential growth but a new regime of dynamo characterized by a linear-in-time growth of magnetic field strength, and the resulting magnetic energy is much higher than the kinetic energy carried by viscous-scale eddies. Unlike the kinematic stage, the subsequent nonlinear stage is unaffected by microscopic diffusion processes and has a universal linear-in-time growth of magnetic energy with the growth rate as a constant fraction 3/38 of the turbulent energy transfer rate, showing good agreement with earlier numerical results. Applying the analysis to the first stars and galaxies, we find that the kinematic stage is able to generate a field strength only an order of magnitude smaller than the final saturation value. But the generation of large-scale magnetic fields can only be accounted for by the relatively inefficient nonlinear stage and requires longer time than the free-fall time. It suggests that magnetic fields may not have played a dynamically important role during the formation of the first stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

Intermittency and geometrical statistics of three-dimensional homogeneous magnetohydrodynamic turbulence: A wavelet viewpoint

Scale-dependent and geometrical statistics of three-dimensional incompressible homogeneous magnetohydrodynamic turbulence without mean magnetic field are examined by means of the orthogonal wavelet decomposition. The flow is computed by direct numerical simulation with a Fourier spectral method at resolution 512{sup 3} and a unit magnetic Prandtl number. Scale-dependent second and higher order statistics of the velocity and magnetic fields allow to quantify their intermittency in terms of spatial fluctuations of the energy spectra, the flatness, and the probability distribution functions at different scales. Different scale-dependent relative helicities, e.g., kinetic, cross, and magnetic relative helicities, yield geometrical information on alignment between the different scale-dependent fields. At each scale, the alignment between the velocity and magnetic field is found to be more pronounced than the other alignments considered here, i.e., the scale-dependent alignment between the velocity and vorticity, the scale-dependent alignment between the magnetic field and its vector potential, and the scale-dependent alignment between the magnetic field and the current density. Finally, statistical scale-dependent analyses of both Eulerian and Lagrangian accelerations and the corresponding time-derivatives of the magnetic field are performed. It is found that the Lagrangian acceleration does not exhibit substantially stronger intermittency compared to the Eulerian acceleration, in contrast to hydrodynamic turbulence where the Lagrangian acceleration shows much stronger intermittency than the Eulerian acceleration. The Eulerian time-derivative of the magnetic field is more intermittent than the Lagrangian time-derivative of the magnetic field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Laminar free convection boundary layer heat transfer over non-isothermal surface

Accurate heat transfer results are provided for the case of nonisothermal objects. A steady, laminar, free convection boundary layer flow over two-dimensional or rotationally symmetrical bodies of nonuniform surface temperature situated in an ambient fluid of undisturbed temperature is considered analytically. The surface heat flux is given in terms of the Nusselt number and wall derivatives of universal functions for Prandtl numbers of 0.72 and 100 are provided. The method is shown to be valid up to a temperature/radius ratio of 130 deg.

Lin, F. N.↗
NASA NTRS1982

Finite analytic numerical solution of axisymmetric Navier-Stokes and energy equations

Convective heat transfer for steady state laminar flow in axisymmetric coordinates is considered. Numerical solutions for flow pattern and temperature distribution are obtained by the finite analytic numerical method applied to the Navier-Stokes equations expressed in terms of vorticity and stream function, and the energy equation. The finite analytic numerical method differs from other numerical methods in that it utilizes a local analytic solution in an element of the problem to construct the total numerical solution. Finite analytic solutions of vorticity, stream function, temperature and heat transfer coefficients for flow with Reynolds number of 5.0, 100.0, 1000.0, and 2000.0, and Prandtl number of 0.1, 1.0, and 10.0 with uniform grid sizes are reported for axisymmetric pipe with sudden expansion and contraction. The wall temperature is considered to be isothermal and differs from the inlet temperature. It is shown that the finite analytic solution is stable, converges rapidly, and simulates the convection of fluid flow accurately since the local analytic solution is capable of simulating automatically the influence of skewed convection through the element boundary on the interior nodal values thereby minimizing the false numerical diffusion.

Chen, C.-J.↗
NASA NTRS1994

Multilayer fluid dynamics of immiscible liquids

A comprehensive analytical, numerical and experimental study of convective flow in multilayer immiscible liquids has been performed. Studies include transparent high Prandtl number of liquid as well as opaque low Prandtl metallic melts. A new radioscopic flow visualization has been developed for the latter studies.

Koster, Jean N.↗
NASA NTRS1959

Analysis of Turbulent Flow and Heat Transfer in Noncircular Passages

Previous work on turbulent heat transfer and flow in tubes was generalized and applied to flow in non-circular passages of equilateral triangular and square cross section. Expressions for eddy diffusivity that had been verified for flow and heat transfer in tubes were assumed to apply in general along lines normal to a wall. Velocity distributions, wall shear-stress distributions, and friction factors, as well as wall heat-transfer distributions, wall temperature distributions, and average heat-transfer coefficients were calculated. In addition, results from a previous analysis for axial flow between rods were compared with new experimental data. For calculating wall temperature distributions, uniform heat generation in the passage wall and uniform heat transfer at the outer surface were assumed. The application of the results is restricted to moderately small peripheral wall temperature variations. Calculations were made for Reynolds numbers from 20,000 to 900,000 and Prandtl numbers from 0.73 to 300. Results show that velocities, shear stresses, and heat transfer in the region near the corner were lower than average values and were zero at the corner. Friction factors and average Nusselt numbers were lower than in a tube.

Taylor, Maynard F.↗