SEARCH · Search NASA

Results for “Peclet 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A note on thermal history kernel for unsteady heat transfer of a spherical particle

When a particle is subjected to an unsteady ambient flow, in terms of either time-dependent relative velocity or time-dependent temperature difference, the net heat transfer from the particle cannot be calculated based on the quasi-steady heat transfer correlation alone. Due to unsteady evolution of the thermal boundary layer, there is also a history contribution to heat transfer. The history contribution to heat transfer is expressed as a convolution integral of past evolution of temperature difference between the particle and the surrounding. While Basset history force and its finite Reynolds number extension have been well studied, similar understanding of unsteady heat transfer and thermal history kernel is lacking. Here, we use existing particle-resolved simulation results to develop a finite Peclet number thermal history kernel, which when used with the convolution integral is demonstrated to accurately predict unsteady heat transfer over a range of Peclet numbers and particle-to-fluid heat capacity ratio.

42 ENGINEERING↗

NASA NTRS1980

Velocity and concentration distribution in a Stefan diffusion tube

Numerical solutions for the diffusion- and gravity-driven flow in a cylindrical Stefan tube were obtained from the coupled diffusion and Navier-Stokes equations of Peclet numbers 0.3, 1 and 5. Distributions of binary component concentrations and velocities were calculated. The mass average velocity is parabolic in nature, except at high Peclet numbers. The solvent is not stagnant but recirculates, even in the absence of gravity. Radial concentration gradients develop which act convectively destabilizing. Consequences for the deduction of diffusion coefficients from Stefan tube experiments are discussed.

Markham, B. L.↗

NASA NTRS2003

Convective Effects During Diffusivity Measurements in Liquids with An Applied Magnetic Field

Convective contamination of self-diffusion experiments with an applied magnetic field is considered using a two-dimensional axisymmetric model. Constant, uniform, and an additional non-uniform heat fluxes are imposed along the sidewall of the cylinder while constant heat loss occurs through the top and bottom. In this model, due to a very small thermal Peclet number, convective heat transfer is neglected, and the flow is steady and inertialess. Time-dependent concentration is solved for various values of the mass Peclet number, Pe(sub m), (the ratio between the convective transport rate and the diffusive transport rate) and different magnetic field strengths represented by the Hartmann number Ha. Normalized values of these diffusivities vs. effective Pe(sub m) are presented for different imposed temperature profiles. In all cases, the diffusivity value obtained through the simulated measurement increases as the effective Pe(sub m) increases. The numerical results suggest that an additional periodic flux, or hot and cold spots, can significantly decrease the convective contamination in our geometry.

Khine, Yu Yu↗

DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

NASA NTRS1975

Zero-derivative boundary condition for pulsed distributed systems

To permit use of experimentally determined Peclet numbers in numerical simulations of pulsed distributed flow systems such as chromatograph columns, substitution of the zero-derivative boundary condition for the infinite boundary condition used in treating data is examined. Moment analysis shows that application of the zero-derivative condition external to the column will yield equivalent numerical results for the two boundary conditions. Criteria for locating this position are provided as a function of the Peclet number.

Lashmet, P. K.↗

NASA NTRS1988

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.↗

NASA NTRS1989

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

Sonder, Leslie J.↗

NASA NTRS1980

On the skewness of the temperature derivative in turbulent flows

An attempt is made to explain the fact that, contrary to the requirements of local isotropy, the skewness of the streamwise temperature derivative has been observed to be a nonzero constant of magnitude of about unity in high-Reynolds-number and high-Peclet-number turbulent shear flows. Measurements in slightly heated homogeneous shear flows and in unsheared grid turbulence suggest that the skewness is nonzero only when the mean shear and the mean temperature gradient are both nonzero.

Sreenivasan, K. R.↗

NASA NTRS1956

Review of experimental investigations of liquid-metal heat transfer

Experimental data of various investigators of liquid-metal heat-transfer characteristics were reevaluated using as consistent assumptions and methods as possible and then compared with each other and with theoretical results. The reevaluated data for both local fully developed and average Nusselt numbers in the turbulent flow region were found still to have considerable spread, with the bulk of the data being lower than predicted by existing analysis. An equation based on empirical grounds which represents most of the fully developed heat-transfer data is nu = 0.625 pe(0.4) where nu represents the Nusselt number and pe the Peclet number. The theoretical prediction of the heat transfer in the entrance region was found to give lower values, in most cases, than those found in the experimental work.

Lubarsky, Bernard↗

NASA NTRS1982

Finite analytic numerical solution of heat transfer and flow past a square channel cavity

A numerical solution of flow and heat transfer characteristics is obtained by the finite analytic method for a two dimensional laminar channel flow over a two-dimensional square cavity. The finite analytic method utilizes the local analytic solution in a small element of the problem region to form the algebraic equation relating an interior nodal value with its surrounding nodal values. Stable and rapidly converged solutions were obtained for Reynolds numbers ranging to 1000 and Prandtl number to 10. Streamfunction, vorticity and temperature profiles are solved. Local and mean Nusselt number are given. It is found that the separation streamlines between the cavity and channel flow are concave into the cavity at low Reynolds number and convex at high Reynolds number (Re greater than 100) and for square cavity the mean Nusselt number may be approximately correlated with Peclet number as Nu(m) = 0.365 Pe exp 0.2.

Chen, C.-J.↗

NASA NTRS1995

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

NASA NTRS1990

Thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity

Computer simulation of steady-state axisymmetrical heat transfer and fluid flow was conducted to study thermocapillary flow and melt/solid interfaces in floating-zone crystal growth under microgravity. The effects of key variables on the extent of thermocapillary flow in the melt zone, the shapes of melt/solid interfaces and the length of the melt zone were discussed. These variables are: (1) the temperature coefficient of surface tension (or the Marangoni number), (2) the pulling speed (or the Peclet number), (3) the feed rod radius, (4) the ambient temperature distribution, (5) the heat transfer coefficient (or the Biot number), and (6) the thermal diffusivity of the material (or the Prandtl number).

Lan, C. W.↗

NASA NTRS2004

The formation of spikes in the displacement of miscible fluids

We report on experiments in which a more viscous fluid displaces a less viscous one in a vertical cylindrical tube. These experiments were performed using silicone oils in a vertical pipette of small diameter. The more viscous fluid also had a slightly larger density than the less viscous fluid. In the initial configuration, the fluids were at rest, and the interface was nominally flat. A dye was added to the more viscous fluid for ease of observation of the interface between the fluids. The flow was initiated by pumping the more viscous fluid into the less viscous one. The displacement velocity was such that the Reynolds number was smaller than unity and the Peclet number for mass transfer between the fluids was large compared to unity. For upward displacement of the more viscous fluid from an initially stable configuration, an axisymmetric finger was observed under all conditions. However, a needle-shaped spike was seen to propagate from the main finger in many cases, similar to that observed by Petitjeans and Maxworthy for the displacement of a more viscous fluid by a less viscous one.

Rheology/methods↗

NASA NTRS2005

The Dynamics of Miscible Interfaces: A Space Flight Experiment

This report presents experiments of a less viscous liquid displaced by a more viscous, miscible liquid in a vertical tube of small diameter. The liquids used are silicone oils of various viscosities. The more viscous fluid has a higher density than the less viscous fluid. The experiments showed that the interface between the liquids can be unstable and exhibits an asymmetric, sinuous shape for downward displacement. The Reynolds number is negligibly small in these experiments and the Peclet number for mass transfer is O(50)-0(l0(exp 5)). The influence of gravity is characterized by the dimensionless parameter F=g(Delta)(rho)d(sup 2)/(mu)U where g is the acceleration due to gravity, Delta(rho) is the density difference between the fluids, d is the tube diameter, mu is the viscosity of the more viscous fluid, and U is the displacement speed. Visualization and PIV measurement of the unstable linger are reported. Experiments were also conducted to study immiscible displacements in tubes in order to compare with previously obtained results for miscible fluids.

Maxworthy, Tony↗

NASA NTRS1982

Finite element analysis of nonisothermal polymer processing operations

A finite element formulation for the analysis of polymer processing is presented and its use in some typical situation including entry flow, transient Couette flow, and the Graetz (forced convection) problem is illustrated. The element formulations are constructed on the premise that momentum convection can be neglected (polymer melt flows typically have very low Reynolds' numbers), but that convective heat transfer may be significant (high Peclet numbers). Nonisothermal effects are considered important in polymer processing, due in part to the significant heating which may occur due to viscous dissipation, and also to the very strong influence of temperature on fluid viscosity. The flow is treated as Newtonian with the flow field being coupled to the heat transfer equation only through the viscous heat generation.

Douglas, C.↗

NASA NTRS2011

Stellar Mixing: I. Formalism

In this paper we use the Reynolds stress models (RSM) to derive algebraic expressions for the following variables: a) heat fluxes; b) J.l fluxes; and c) momentum fluxes. These relations, which are fully 3D, include: 1) stable and unstable stratification, represented by the Brunt-Vaislila frequency, N(exp 2) =-g/H(sub p_(del - del(sub ad))(1 - RI(sub mu)); 2) double diffusion, salt-fingers, and semi-convection, represented by the density ratio R(sub mu) = del(sub mu)/(del - del(sub ad)); 3) shear (differential rotation), represented by the mean squared shear Sigma(exp 2) or by the Richardson number, Ri =N(exp 2)Sigma(exp -2); 4) radiative losses represented by a Peclet number, Pe; 5) a complete analytical solution of the ID version of the model. In general, the model requires the solution of two differential equations for the eddy kinetic energy K and its rate of dissipation, epsilon. In the local and stationary cases, when production equals dissipation, the model equations are all algebraic.

Canuto, V .M.↗

DOE OSTI2022

On the hot surface ignition of a wall-stagnating spray flame

The ignition of fuel sprays due to interaction with hot surfaces is an important phenomenon in the safety analysis of many engineering systems. We perform a parametric study of the hot surface ignition (HSI) of a fuel spray approaching a heated surface caused by the accidental leakage of a fuel line. Here, to this end, we employ a one-dimensional Eulerian-Eulerian formulation with a non-equilibrium evaporation model and a realistic chemical mechanism to describe $n$-dodecane fuel chemistry. We first describe and analyze the phenomenology of the unsteady processes leading to ignition using non-dimensionalized quantities. Through consideration of the temporal development of the most reactive mixture, we demonstrate that ignition occurs at a fuel-lean composition in a premixed region near the hot surface. Using non-dimensional parameters identified from the governing equations, we perform a parametric study of the time, location and local mixture composition at ignition and determine the ignition limits. We then identify the most important parametric sensitivities for physical analysis using a data-driven classification method. Our analysis demonstrates a contraction of the ignition limits with increased Stokes number and a regime of parametric insensitivity of igniting mixture composition. We also show that at high Damköhler numbers, the ignition location conforms to the parametric behavior of the thermal boundary layer, whereas at low Damköhler numbers approaching the ignition limit it reaches a near-unity value of the quenching Peclet number. We then compare the demonstrated parametric dependencies to the results of the quasi-steady asymptotic ignition literature, showing that our results are consistent with those obtained analytically within the limitations imposed by the simplified formulation of the latter.

42 ENGINEERING↗

NASA NTRS1988

Scalar entrainment in the mixing layer

New definitions of entrainment and mixing based on the passive scalar field in the plane mixing layer are proposed. The definitions distinguish clearly between three fluid states: (1) unmixed fluid, (2) fluid engulfed in the mixing layer, trapped between two scalar contours, and (3) mixed fluid. The difference betwen (2) and (3) is the amount of fluid which has been engulfed during the pairing process, but has not yet mixed. Trends are identified from direct numerical simulations and extensions to high Reynolds number mixing layers are made in terms of the Broadwell-Breidenthal mixing model. In the limit of high Peclet number (Pe = ReSc) it is speculated that engulfed fluid rises in steps associated with pairings, introducing unmixed fluid into the large scale structures, where it is eventually mixed at the Kolmogorov scale. From this viewpoint, pairing is a prerequisite for mixing in the turbulent plane mixing layer.

Sandham, N. D.↗