Search NASA⌕ Search

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.

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

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Dynamics of hard sphere colloidal dispersions

Our objective is to perform on homogeneous, fully equilibrated dispersions the full set of experiments characterizing the transition from fluid to solid and the properties of the crystalline and glassy solid. These include measurements quantifying the nucleation and growth of crystallites, the structure of the initial fluid and the fully crystalline solid, and Brownian motion of particles within the crystal, and the elasticity of the crystal and the glass. Experiments are being built and tested for ideal microgravity environment. Here we describe the ground based effort, which exploits a fluidized bed to create a homogeneous, steady dispersion for the studies. The differences between the microgravity environment and the fluidized bed is gauged by the Peclet number Pe, which measures the rate of convection/sedimentation relative to Brownian motion. We have designed our experiment to accomplish three types of measurements on hard sphere suspensions in a fluidized bed: the static scattering intensity as a function of angle to determine the structure factor, the temporal autocorrelation function at all scattering angles to probe the dynamics, and the amplitude of the response to an oscillatory forcing to deduce the low frequency viscoelasticity. Thus the scattering instrument and the colloidal dispersion were chosen such as that the important features of each physical property lie within the detectable range for each measurement.

Zhu, J. X.↗
NASA NTRS1994

Penetration below a convective zone

Two-dimensional numerical simulations are used to investigate how fully compressible nonlinear convection penetrates into a stably stratified zone beneath a stellar convection zone. Estimates are obtained of the extent of penetration as the relative stability S of the stable to the unstable zone is varied over a broad range. The model deals with a perfect gas possessing a constant dynamic viscosity. The dynamics is dominated by downward-directed plumes which can extend far into the stable material and which can lead to the excitation of a broad spectrum of internal gravity waves in the lower stable zone. The convection is highly time dependent, with the close coupling between the lateral swaying of the plumes and the internal gravity waves they generate serving to modulate the strength of the convection. The depth of penetration delta, determined by the position where the time-averaged kinetic flux has its first zero in the stable layer, is controlled by a balance between the kinetic energy carried into the stable layer by the plumes and the buoyancy braking they experience there. A passive scalar is introduced into the unstable layer to evaluate the transport of chemical species downward. Such a tracer is effectively mixed within a few convective overturning times down to a depth of delta within the stable layer. Analytical estimates based on simple scaling laws are used to interpret the variation of delta with S, showing that it first involves an interval of adiabatic penetration if the local Peclet number of the convection exceeds unity, followed by a further thermal adjustment layer, the depths of each interval scaling in turn as S(exp -1) and S(exp -1/4). These estimates are in accord with the penetration results from the simulations.

Hurlburt, Neal E.↗
NASA NTRS1990

Direct Numerical Simulations of Turbulent Convection: Variable Gravity and Uniform Rotation - Part 1

Turbulent convection may have played a major role in determining the structure and evolution of the primordial solar nebula, but current, incomplete models of convection and turbulence give very different results and remain largely untested in the absence of detailed astronomical observations. Numerical simulations provide an 'experimental' database for comparison with these models, and, to this end, direct numerical simulations of turbulent convection were performed with modifications intended to mimic some of the unique physical features of thin accretion disks, such as the primordial solar nebula: internal heating, a gravitational acceleration that is linearly proportional to the distance from midplane of the nebula, and rapid rotation. Peclet numbers in the simulations are comparable to those in solar nebula models; Rossby numbers in the simulations are an order of magnitude larger than those in solar nebula models because of the unrealistically high Prandtl and low, Reynolds numbers required to resolve all scales of the convective flow. We find that, despite the loss of buoyancy at midplane, turbulent motions easily penetrate the midplane region with little loss of intensity, providing efficient transport of heat and turbulent kinetic energy throughout the interior. A simple mixing length model modified to include rotation is found to give convective heat fluxes for the interior flow in rough agreement with the numerical simulations. We discuss these preliminary results with regard to assumptions about heating distributions and convective heat fluxes made in standard solar nebula models. More definitive comparisons with solar nebula modelling will become possible when more realistic effects of shear, density, stratification, and compression are included.

Cabot, W.↗
NASA NTRS1996

Nonlinear Response of Layer Growth Dynamics in the Mixed Kinetics-Bulk-Transport Regime

In situ high-resolution interferometry on horizontal facets of the protein lysozyme reveal that the local growth rate R, vicinal slope p, and tangential (step) velocity v fluctuate by up to 80% of their average values. The time scale of these fluctuations, which occur under steady bulk transport conditions through the formation and decay of step bunches (macrosteps), is of the order of 10 min. The fluctuation amplitude of R increases with growth rate (supersaturation) and crystal size, while the amplitude of the v and p fluctuations changes relatively little. Based on a stability analysis for equidistant step trains in the mixed transport-interface-kinetics regime, we argue that the fluctuations originate from the coupling of bulk transport with nonlinear interface kinetics. Furthermore, step bunches moving across the interface in the direction of or opposite to the buoyancy-driven convective flow increase or decrease in height, respectively. This is in agreement with analytical treatments of the interaction of moving steps with solution flow. Major excursions in growth rate are associated with the formation of lattice defects (striations). We show that, in general, the system-dependent kinetic Peclet number, Pe(sub k) , i.e., the relative weight of bulk transport and interface kinetics in the control of the growth process, governs the step bunching dynamics. Since Pe(sub k) can be modified by either forced solution flow or suppression of buoyancy-driven convection under reduced gravity, this model provides a rationale for the choice of specific transport conditions to minimize the formation of compositional inhomogeneities under steady bulk nutrient crystallization conditions.

Vekilov, Peter G.↗
NASA NTRS2000

The Dynamics of Miscible Interfaces: A Space Flight Experiment

Experiments as well as accompanying simulations are described that serve in preparation of a space flight experiment to study the dynamics of miscible interfaces. The investigation specifically addresses the importance of both nonsolenoidal effects as well as nonconventional Korteweg stresses in flows that give rise to steep but finite concentration gradients. The investigation focuses on the flow in which a less viscous fluid displaces one of higher viscosity and different density within a narrow capillary tube. The fluids are miscible in all proportions. An intruding finger forms that occupies a fraction of the total tube diameter. Depending on the flow conditions, as expressed by the Peclet number, a dimensionless viscosity ratio, and a gravity parameter, this fraction can vary between approximately 0.9 and 0.2. For large Pe values, a quasi-steady finger forms, which persists for a time of O(Pe) before it starts to decay, and Poiseuille flow and Taylor dispersion are approached asymptotically. Depending on the specific flow conditions, we observe a variety of topologically different streamline patterns, among them some that leak fluid from the finger tip. For small Pe values, the flow decays from the start and asymptotically reaches Taylor dispersion after a time of O(Pe). Comparisons between experiments and numerical simulations based on the 'conventional' assumption of solenoidal velocity fields and without Korteweg stresses yield poor agreement as far as the Pe value is concerned that distinguishes these two regimes. As one possibility, we attribute this lack of agreement to the disregard of these terms. An attempt is made to use scaling arguments in order to evaluate the importance of the Korteweg stresses and of the assumption of solenoidality. While these effects should be strongest in absolute terms when steep concentration fronts exist, i.e., at large Pe, they may be relatively most important at lower values of Pe. We subsequently compare these conventional simulations to more complete simulations that account for nonvanishing divergence as well as Korteweg stresses. While the exact value of the relevant stress coefficients are not known, ballpark numbers do exist, and their use in the simulations indicates that these stresses may indeed be important. We plan to evaluate these issues in detail by means of comparing a space experiment with corresponding simulations, in order to extract more accurate Korteweg stress coefficients, and to confirm or deny the importance of such stresses.

Maxworthy, Tony↗
NASA NTRS2001

The Dynamics of Miscible Fluids: A Space Flight Experiment (MIDAS)

We propose a space flight experiment to study the dynamics of miscible interfaces. A less viscous fluid displaces one of higher viscosity within a tube. The two fluids are miscible in all proportions. An intruding "finger" forms that occupies a fraction of the tube. As time progresses diffusion at the interface combined with flow induced straining between the two fluids modifies the concentration and velocity distributions within the whole tube. Also, under such circumstances it has been proposed that the interfacial stresses could depend on the local concentration gradients (Korteweg stresses) and that the divergence of the velocity need not be zero, even though the flow is incompressible. We have obtained reasonable agreement for the tip velocity between numerical simulations (that ignored the Korteweg stress and divergence effects) and physical experiments only at high Peclet Numbers. However at moderate to low Pe agreement was poor. As one possibility we attributed this lack of agreement to the disregard of these effects. We propose a space experiment to measure the finger shape, tip velocity, and the velocity and concentration fields. From intercomparisons between the experiment and the calculations we can then extract values for the coefficients of the Korteweg stress terms and confirm or deny the importance of these stresses.

Maxworthy, T.↗
NASA NTRS2003

Free Dendritic Growth of Succinonitrile-Acetone Alloys with Thermosolutal Melt Convection

A stagnant film model of the effects of thermosolutal convection on free dendritic growth of alloys is developed, and its predictions are compared to available earth-based experimental data for succinonitrileacetone alloys. It is found that the convection model gives excellent agreement with the measured dendrite tip velocities and radii for low solute concentrations. However, at higher solute concentrations the present predictions show some deviations from the measured data, and the measured (thermal) Peclet numbers tend to fall even below the predictions from diffusion theory. Furthermore, the measured selection parameter (sigma*) is significantly above the expected value of 0.02 and exhibits strong scatter. It is shown that convection is not responsible for these discrepancies. Some of the deviations between the predicted and measured data at higher supercoolings could be caused by measurement difficulties. The systematic disagreement in the selection parameter for higher solute concentrations and all supercoolings examined, indicates that the theory for the selection of the dendrite tip operating state in alloys may need to be reexamined.

Beckermann, Christoph↗
NASA NTRS2002

Solutal Convection Around Growing Protein Crystal and Diffusional Purification in Space

This work theoretically addressed two subjects: 1) onset of convection, 2) distribution of impurities. Onset of convection was considered analytically and numerically. Crystal growth was characterized by slow surface incorporation kinetics, i.e. growth kinetic coefficient beta (cm/s) small as compared to the typical bulk diffusion rate, D(sub 1)/h, where D(sub 1) is diffusivity of major crystallizing protein and h is the crystal size. Scaling type analysis predicted two laws on how the convection rate, v, essentially the Peclet number, Pe exactly equal to vh/D(sub 1), depends on dimensionless kinetic coefficient a exactly equal to beta h/D(sub 1). Namely: Pe = C(sub 2/5)(aRa(sup 2/5)) and Pe = C(sub 1) aRa. Here, Reynolds number Ra = rho(sub 1)(sup 0)gh(sup 3)(rho(sub p) - rho(sub w))/rho(sup p)rho(sub 1)vD(sub 1), v being solution viscosity. The constants C(sub 2/5), exactly equal to 0.28 and C(sub 1) exactly equal to 10(exp -2) found from the full scale computer simulation for a cylindrical crystal inside big cylindrical vessel. The linear boundary conditions connecting protein and impurity concentration at the interface with the flux to/from the interface was applied. No-slip condition for Navier-Shocker equations was employed. With these conditions, flow and concentration distributions were calculated. Validity of the Pe(Ra) dependencies follows for wide range of parameters for which numerical calculations have been accomplished and presented by various points.

Chernov, A. A.↗
NASA NTRS2003

Concentration and Velocity Gradients in Fluidized Beds

In this work we focus on the height dependence of particle concentration, average velocity components, fluctuations in these velocities and, with the flow turned off, the sedimentation velocity. The latter quantities are measured using Particle Imaging Velocimetry (PIV). The PIV technique uses a 1-megapixel camera to capture two time-displaced images of particles in the bed. The depth of field of the imaging system is approximately 0.5 cm. The camera images a region with characteristic length of 2.6 cm for the small particles and 4.7 cm. for the large particles. The local direction of particle flow is determined by calculating the correlation function for sub-regions of 32 x 32 pixels. The velocity vector map is created from this correlation function using the time between images (we use 15 to 30 ms). The software is sensitive variations of 1/64th of a pixel. We produce velocity maps at various heights, each consisting of 3844 velocities. We break this map into three vertical zones for increased height information. The concentration profile is measured using an expanded (1 cm diameter) linearly polarized HeNe Laser incident on the fluidized bed. A COHU camera (gamma=1, AGC off) with a lens and a polarizer images the transmitted linearly polarized light to minimize the effects of multiply scattered light. The intensity profile (640 X 480 pixels) is well described by a Gaussian fit and the height of the Gaussian is used to characterize the concentration. This value is compared to the heights found for known concentrations. The sedimentation velocity is estimated using by imaging a region near the bottom of the bed and using PIV to measure the velocity as a function of time. With a nearly uniform concentration profile, the time can be converted to height information. The stable fluidized beds are made from large pseudo-monodisperse particles (silica spheres with radii (250-300) microns and (425-500) microns) dispersed in a glycerin/water mix. The Peclet number is sufficiently large that Brownian motion of the particles can be ignored and the Reynolds number sufficiently small that particle inertia is negligible. A packed particle bed is used to randomize and disperse the flowing fluid introduced by a peristaltic pump. The bed itself is a rectangular glass cell 8 cm wide (x), 0.8 cm deep and a height of 30.5 cm (z). The depth of field of the camera is approximately 0.5 cm so depth information is averaged. Over flow fluid is returned to the reservoir making a closed loop system. In these experiments the particles form a sediment approximately 5.7 cm high with the pump off and expand to 22 cm with the pump on. For the smaller particles the pump velocity is .5 millimeters per second and 1.1 millimeters per second for the large particles. At this concentration the bed has a very well defined top where particle concentration rapidly drops to zero.

McClymer, James P.↗
NASA NTRS2000

Thermal Edge-Effects Model for Automated Tape Placement of Thermoplastic Composites

Two-dimensional thermal models for automated tape placement (ATP) of thermoplastic composites neglect the diffusive heat transport that occurs between the newly placed tape and the cool substrate beside it. Such lateral transport can cool the tape edges prematurely and weaken the bond. The three-dimensional, steady state, thermal transport equation is solved by the Green's function method for a tape of finite width being placed on an infinitely wide substrate. The isotherm for the glass transition temperature on the weld interface is used to determine the distance inward from the tape edge that is prematurely cooled, called the cooling incursion Delta a. For the Langley ATP robot, Delta a = 0.4 mm for a unidirectional lay-up of PEEK/carbon fiber composite, and Delta a = 1.2 mm for an isotropic lay-up. A formula for Delta a is developed and applied to a wide range of operating conditions. A surprise finding is that Delta a need not decrease as the Peclet number Pe becomes very large, where Pe is the dimensionless ratio of inertial to diffusive heat transport. Conformable rollers that increase the consolidation length would also increase Delta a, unless other changes are made, such as proportionally increasing the material speed. To compensate for premature edge cooling, the thermal input could be extended past the tape edges by the amount Delta a. This method should help achieve uniform weld strength and crystallinity across the width of the tape.

Costen, Robert C.↗
NASA NTRS1998

Tape-Drop Transient Model for In-Situ Automated Tape Placement of Thermoplastic Composites

Composite parts of nonuniform thickness can be fabricated by in-situ automated tape placement (ATP) if the tape can be started and stopped at interior points of the part instead of always at its edges. This technique is termed start/stop-on-the-part, or, alternatively, tape-add/tape-drop. The resulting thermal transients need to be managed in order to achieve net shape and maintain uniform interlaminar weld strength and crystallinity. Starting-on-the-part has been treated previously. This paper continues the study with a thermal analysis of stopping-on-the-part. The thermal source is switched off when the trailing end of the tape enters the nip region of the laydown/consolidation head. The thermal transient is determined by a Fourier-Laplace transform solution of the two-dimensional, time-dependent thermal transport equation. This solution requires that the Peclet number Pe (the dimensionless ratio of inertial to diffusive heat transport) be independent of time and much greater than 1. Plotted isotherms show that the trailing tape-end cools more rapidly than the downstream portions of tape. This cooling can weaken the bond near the tape end; however the length of the affected region is found to be less than 2 mm. To achieve net shape, the consolidation head must continue to move after cut-off until the temperature on the weld interface decreases to the glass transition temperature. The time and elapsed distance for this condition to occur are computed for the Langley ATP robot applying PEEK/carbon fiber composite tape and for two upgrades in robot performance. The elapsed distance after cut-off ranges from about 1 mm for the present robot to about 1 cm for the second upgrade.

Costen, Robert C.↗
NASA NTRS2004

Structure of Propagating and Attached Hydrocarbon Flames

Direct numerical simulations with C3-chemistry and radiative heat-loss models have been performed to reveal the internal structure of propagating and attached flames in an axisymmetric fuel jet of methane, ethane, ethylene, acetylene, or propane in air under normal and zero gravity. Observations of the flames were also made at the NASA Glenn 2.2-Second Drop Tower. In computations, the fuel issued into quasi-quiescent air for a fixed mixing time before it was ignited along the centerline at stoichiometry. The edge of the flame propagated through a flammable layer at the laminar flame speed of the stoichiometric fuel-air mixture independent of gravity. For all cases, a peak reactivity spot, i.e., reaction kernel, was formed in the flame base, thereby holding a trailing diffusion flame. The location of the reaction kernel in the attached flames depended inversely on the reactivity. The reaction-kernel correlations between the reactivity and the velocity were developed further using variables related to local Damkahler and Peclet numbers.

Takahashi, Fumiaki↗
NASA NTRS2007

Catalytic Ignition and Upstream Reaction Propagation in a Platinum Tube

A challenge for catalytic combustion in monolithic reactors at elevated temperatures is the start-up or "light-off" from a cold initial condition. In this work, we demonstrate a concept called "back-end catalytic ignition that potentially can be utilized in the light-off of catalytic monoliths. An external downstream flame or Joule heating raises the temperature of a small portion of the catalyst near the outlet initiating a localized catalytic reaction that propagates upstream heating the entire channel. This work uses a transient numerical model to demonstrate "back-end" ignition within a single channel which can characterize the overall performance of a monolith. The paper presents comparisons to an experiment using a single non-adiabatic channel but the concept can be extended to the adiabatic monolith case. In the model, the time scales associated with solid heat-up are typically several orders of magnitude larger than the gas-phase and chemical kinetic time-scales. Therefore, the model assumes a quasi-steady gas-phase with respect to a transient solid. The gas phase is one-dimensional. Appropriate correlations, however, account for heat and mass transfer in a direction perpendicular to the flow. The thermally-thin solid includes axial conduction. The gas phase, however, does not include axial conduction due to the high Peclet number flows. The model includes both detailed gas-phase and catalytic surface reactions. The experiment utilizes a pure platinum circular channel oriented horizontally though which a CO/O2 mixture (equivalence ratios ranging from 0.6 to 0.9) flows at 2 m/s.

Struk, P. M.↗
NASA NTRS1999

Shape Parameter for a Non-Axisymmetric Isothermal Dendrite

In previous work, we found approximate solutions for paraboloids having perturbations with four-fold axial symmetry in order to model dendritic growth in cubic materials. These solutions provide self-consistent corrections through second order in a shape parameter e to the Peclet number-supercooling relation of the Ivantsov solution. The parameter e is proportional to the amplitude of the four-fold correction to the dendrite shape, as measured from the Ivantsov paraboloid of revolution. We calculate e by comparing the dendrite tip shape to the portion of the equilibrium shape near the growth direction, (001), for anisotropic surface free energy, where the ni are components of the unit normal of the crystal surface. This comparison results in epsilon = -2(epsilon 4), independent of the Peclet number. From the experimental value of epsilon 4, we find epsilon approximately 0.011, in good agreement with the measured value epsilon approximately 0.008 of LaCombe et al.

McFadden, G. B.↗
NASA NTRS1991

Fractal interfaces and product generation in the two-dimensional mixing layer

The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two-dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This reaction is verified using product generation measurements and fractal dimensions derived from the box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion of the flow field during the first pairing interaction. It is suggested that the structures reponsible for this transition are nondeterministic, nonrandom, inhomogeneous fractals. In the range of Schmidt numbers investigated (0.25-4), only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Jimenez, Javier↗
NASA NTRS1990

A fractal transition in the two dimensional shear layer

The dependence of product generation with the Peclet and Reynolds number in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is used to compute the fractal dimension of passive scalar interfaces. A transition from a low dimension of 4/3 to a higher one of 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transition is found for Reynolds numbers up to 20,000.

Jimenez, Javier↗
NASA NTRS1991

Fractal interfaces and product generation in the two dimensional mixing layer

The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This relation is verified using product generation measurements and dimensions derived from a standard box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Jimenez, Javier↗
NASA NTRS1986

Numerical modeling of physical vapor transport in a vertical cylindrical ampoule, with and without gravity

Numerical modeling has been performed of the fluid dynamics in a prototypical physical vapor transport crystal growing situation. Cases with and without gravity have been computed. Dependence of the flows upon the dimensionless parameters aspect ratio and Peclet, Rayleigh, and Schmidt numbers is demonstrated to a greater extent than in previous works. Most notably, it is shown that the effects of thermally-induced buoyant convection upon the mass flux on the growth interface crucially depend upon the temperature boundary conditions on the sidewall (e.g., whether adiabatic or of a fixed profile, and in the latter case the results depend upon the shape of the profile assumed).

Miller, T. L.↗