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 55 records · Page 3

Characteristics-based methods applied to infinite Prandtl number thermal convection in the hard turbulent regime

Characteristics-based methods for the advection-diffusion equation are presented and directly applied to study thermal convection with extremely large Rayleigh number (Ra). It is shown that the operator-splitting method for advection-diffusion problems is very accurate for determining the advected field at extremely high Peclet number (Pe). The technique presented is considered to have great potential for solving advection-dominated problems, while the Langrangian method is more accurate for lower Pe. It is noted that the accuracy of these characteristics-based methods strongly depends on the quality of interpolation. The computational time for the operator-splitting method grows with the number of time steps employed. The Langrangian method was used for simulations of convection at very high Ra, up to 3 x 10 to the 9th, and time-dependent, thermal convection solutions were obtained for infinite Prandtl number.

Malevsky, A. V.↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Finite element analyses of plane thermal entry-length flows

It is pointed out that the design of actively cooled engine and airframe structures for hypersonic flight requires strong interaction between fluid, thermal, and structural analyses. Research programs are being conducted to develop an integrated thermal-structural analysis capability based upon the finite element method. In connection with these programs, a finite element engineering approach has been developed based upon a number of assumptions customarily used in practical heat transfer analysis. An evaluation is conducted of the finite element engineering approach for entry-length flows. The evaluation makes use of an analysis of plane thermal entry-length flows by a finite element approach utilizing mean fluid temperatures with a convection coefficient. Another analysis is based on a finite element continuum approach, taking into account the rigorous momentum, mass, and energy equations of viscous, incompressible flow. The engineering approach was found to give excellent agreement with the rigorous continuum approach for Peclet number values not less than 100.

Thornton, E. A.↗

Gastight rotating cylinder electrode: Toward decoupling mass transport and intrinsic kinetics in electrocatalysis

Abstract Decoupling and understanding the various mass, charge, and heat transport phenomena involved in the electrocatalytic transformation of small molecules (i.e., CO 2 , CO, H 2 , N 2 , NH 3 , O 2 , and CH 4 ) is challenging but it can be readily achieved using dimensionless quantities (i.e., Reynolds, Sherwood, Schmidt, Damköhler, Nusselt, Prandtl, and Peclet Numbers) to simplify the characterization of systems with multiple interacting physical phenomena. Herein we report the development of a gastight rotating cylinder electrode cell with well‐defined mass transport characteristics that can be applied to experimentally decouple mass transfer effects from intrinsic kinetics in electrocatalytic systems. The gastight rotating cylinder electrode cell enables the dimensionless analysis of electrocatalytic systems and should enable the rigorous research and development of electrocatalytic technologies.

Jang, Joonbaek↗

Nanoparticle Diffusion in Miscible Polymer Nanocomposite Melts

X-ray photon correlation spectroscopy measurements were used to quantify the dynamics of bare and bimodal grafted silica nanoparticles mixed with PEO melts of different molecular weights. In dilute polymer nanocomposite (PNC) samples, we find diffusive NP behavior as described by the Stokes–Einstein relationship so long as the adsorbed PEO polymer layer is taken into account in determining both the effective NP size and its role on composite viscosity. The size of this bound layer was found to be approximately 2R g , where R g is the chain radius of gyration. We also expanded our system to investigate how the dynamics of grafted NPs differ from bare NPs with an adsorbed layer. We showed that the dynamics again can be determined by an effective NP radius at a scale smaller than the effective interparticle spacing; however, at larger length scales, the morphology and grafting parameters play a major role in the system dynamics. These results allow us to quantify NP ordering driven by polymer crystallization. It has previously been speculated that behavior is controlled by the relative ratio of time scale of crystal growth and the diffusive time scale of the NPs, a Peclet number. When the former time scale is longer, then the NPs are expected to be segregated into the interlamellar amorphous zones, while the NPs are too slow to be reorganized in the opposite case. We show here that this conjecture is quantitatively correct and the demarcation in behavior occurs for Pe = 1. Thus, we provide a way to estimate a critical spherulite growth rate for any semicrystalline PNC, at which a given NP can be ordered. Together, the results of this study permit us to tunably design PNCs through directed dispersion of NPs in a polymer matrix.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Reactive Transport Modeling of Mineral Precipitation and Carbon Trapping in Discrete Fracture Networks

In this study we use numerical experiments to analyze reactive flow and transport behavior in discrete intersecting fracture networks, focusing on (a) how reaction-induced changes in physical and chemical properties affect flow connectivity and (b) how fracture networks developed in the Earth's critical zone contribute to carbon sequestration via mineral weathering reactions. In the first part of the study, we used two-dimensional reactive flow and transport simulations to analyze the impacts of mixing in a natural discrete fracture network. We concluded that reaction-induced changes can substantially alter the flow connectivity, especially at fracture intersections. The second set of simulations considered the problem of natural weathering of fractured mafic and ultramafic rocks in the partially saturated Earth's critical zone as a function of infiltration rates, fracture permeability, and partially saturated flow parameters. As a model system, we considered an incongruent reaction network with dissolution of forsterite and precipitation of magnesite. The behavior is complex in terms of the rate-controlling processes because of the multicomponent nature of the system as shown by the grid Peclet number: the CO 2 behavior is gas diffusion-controlled in the partially saturated zone, while the rate of water flow via the Damkӧhler number controls Mg 2+ transport through the fracture network. The amounts of carbon that can be trapped are modest, but the naturally fractured domain considered here provides a useful “base case” against which various engineered solutions can be compared.

58 GEOSCIENCES↗

Two dimensional topology optimization of heat exchangers with the volume fraction method

We perform a comparison study of two topology optimizations methods applied to the design of a two fluids heat exchanger modeled with a coupled thermal-flow problem. The flow follows an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation without internal generation and high Peclet number. To keep the two fluid phases separated, we solve two Stokes-Brinkman equations, where the Brinkman term models the other phase as a solid. These two velocity fields are then fed to the heat transfer equation. Our goal is to maximize the enthalpy at the cold outlet while constraining the pressure drop. We first solve the design modeling the solid and fluid phases with a volume fraction variable. A SIMP-like penalization in the Brinkman term drives the optimization to a discrete design. The cost and constraint function derivatives are calculated with the library pyadjoint and the optimization is performed by IPOPT. We present optimized designs in two dimensions and discuss the influence of the parameters.

Beck, VictorA.↗

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗

Analysis of impingement heat transfer for two parallel liquid-metal slot jets

An analytical method is developed for determining heat transfer by impinging liquid-metal slot jets. The method involves mapping the jet flow region, which is bounded by free streamlines, into a potential plane where it becomes a uniform flow in a channel of constant width. The energy equation is transformed into potential plane coordinates and is solved in the channel flow region. Conformal mapping is then used to transform the solution back into the physical plane and obtain the desired heat-transfer characteristics. The analysis given here determines the heat-transfer characteristics for two parallel liquid-metal slot jets impinging normally against a uniformly heated flat plate. The liquid-metal assumptions are made that the jets are inviscid and that molecular conduction is dominating heat diffusion. Wall temperature distributions along the heated plate are obtained as a function of spacing between the jets and the jet Peclet number.

Siegel, R.↗

Analysis of liquid-metal-jet impingement cooling in a corner region and for a row of jets

A conformal mapping method was used to analyze liquid-metal-jet impingement heat transfer. The jet flow region and energy equation are transformed to correspond to uniform flow in a parallel plate channel with nonuniform heat addition along a portion of one wall. The exact solution for the wall-temperature distribution was obtained in the transformed channel, and the results are mapped back into the physical plane. Two geometries are analyzed. One is for a single slot jet directed either into an interior corner formed by two flat plates, or over the external sides of the corner; the flat plates are uniformly heated, and the corner can have various included angles. The heat-transfer coefficient at the stagnation point at the apex of the plates is obtained as a function of the corner angle, and temperature distributions are calculated along the heated walls. The second geometry is an infinite row of uniformly spaced parallel slot jets impinging normally against a uniformly heated plate. The heat-transfer behavior is obtained as a function of the spacing between the jets. Results are given for several jet Peclet numbers from 5 to 50.

Siegel, R.↗

Further calculations based upon a theory of flame spread across solid fuels

A theory of flame spread above solid fuel beds of finite or infinite thickness is considered. Heat transfer ahead of the flame is allowed through both the solid and gas phases. The simplifying assumption of a surface reaction is made. Temperature profiles in both phases, concentration profiles in the gas phase, heat transfer at the surface, and flame spreading rates are determined. It is found that chemical kinetic parameters together with transport properties strongly influence the spreading rate. Heat transfer through the gas phase is the dominant mechanism for energy transfer ahead of the spreading flame over a wide range of Peclet numbers for the gas and solid phases.

Feng, C. C.↗

Tabulation and summary of thermodynamic effects data for developed cavitation on ogive-nosed bodies

Thermodynamic effects data for developed cavitation on zero and quarter caliber ogives in Freon 113 and water are tabulated and summarized. These data include temperature depression (delta T), flow coefficient (C sub Q), and various geometrical characteristics of the cavity. For the delta T tests, the free-stream temperature varied from 35 C to 95 C in Freon 113 and from 60 C to 125 C in water for a velocity range of 19.5 m/sec to 36.6 m/sec. Two correlations of the delta T data by the entrainment method are presented. These correlations involve different combinations of the Nusselt, Reynolds, Froude, Weber, and Peclet numbers and dimensionless cavity length.

Holl, J. W.↗

Stellar convection 2: A multi-mode numerical solution for convection in spheres

The convective flow of a self gravitating sphere of Boussinesq fluid for small Reynolds and Peclet numbers is numerically determined. The decomposition of the equations of motion into modes is reviewed and a relaxation method is developed and presented to compute the solutions to these equations. The stable equilibrium flow for a Rayleigh number of 10 to the 4th power and a Prandtl number of 10 is determined. The 2 and 3 dimensional spectra of the kinetic and thermal energies and the convective flux as a function of wavelengths are calculated in terms of modes. The anisotropy of the flow as a function of wavelength is defined.

Marcus, P. S.↗

Boundary integral equation method calculations of surface regression effects in flame spreading

A solid-phase conduction problem that is a modified version of one that has been treated previously in the literature and is applicable to flame spreading over a pyrolyzing fuel is solved using a boundary integral equation (BIE) method. Results are compared to surface temperature measurements that can be found in the literature. In addition, the heat conducted through the solid forward of the flame, the heat transfer responsible for sustaining the flame, is also computed in terms of the Peclet number based on a heated layer depth using the BIE method and approximate methods based on asymptotic expansions. Agreement between computed and experimental results is quite good as is agreement between the BIE and the approximate results.

Altenkirch, R. A.↗

The effect of fluid flow due to the crystal-melt density change on the growth of a parabolic isothermal dendrite

The Ivantsov (1947) analysis of an isolated isothermal dendrite (with zero surface tension) growing into a supercooled liquid is extended to include the effects of the fluid flow due to volume contraction or expansion upon solidification. For an axisymmetric paraboloidal dendrite, an analytic solution to the Navier-Stokes equations is obtained. The magnitude of the flow is proportional to the relative density change epsilon, and the flow becomes negligible far from the surface of the dendrite. The temperature field consistent with this flow can also be found explicitly. The well-known expression that relates the dimensionless supercooling to the Peclet number in the absence of fluid flow is modified for nonzero epsilon, but the effect is of order epsilon and hence is seen to be minor for most values of epsilon and dimensionless supercooling that occur in practice.

Mcfadden, G. B.↗

Vortex-scalar element calculations of a diffusion flame stabilized on a plane mixing layer

The vortex-scalar element method, a scheme which utilizes vortex elements to discretize the region of high vorticity and scalar elements to represent species or temperature fields, is utilized in the numerical simulations of a two-dimensional reacting mixing layer. Computations are performed for a diffusion flame at high Reynolds and Peclet numbers without resorting to turbulence models. In the nonreacting flow, the mean and fluctuation profiles of a conserved scalar show good agreement with experimental measurements. Results for the reacting flow indicate that for temperature independent kinetics, the chemical reaction begins immediately downstream of the splitter plate where mixing starts. Results for the reacting flow with Arrhenius kinetics show an ignition delay, which depends on reactant temperature, before significant chemical reaction occurs. Harmonic forcing changes the structure of the layer, and concomitantly the rates of mixing and reaction, in accordance with experimental results. Strong stretch within the braids in the nonequilibrium kinetics case causes local flame quenching due to the temperature drop associated with the large convective fluxes.

Ghoniem, Ahmed F.↗

CONDIF - A modified central-difference scheme for convective flows

The paper presents a method, called CONDIF, which modifies the CDS (central-difference scheme) by introducing a controlled amount of numerical diffusion based on the local gradients. The numerical diffusion can be adjusted to be negligibly low for most problems. CONDIF results are significantly more accurate than those obtained from the hybrid scheme when the Peclet number is very high and the flow is at large angles to the grid.

Runchal, Akshai K.↗

Fermi acceleration at shocks with arbitrary velocity profiles

Diffusive particle acceleration in modified shock fronts is considered. A method by means of which to obtain the resulting particle spectrum for an arbitrary flow profile is outlined; this method is based on a reformulation of the transport equation as a nonlinear first-order equation. Besides the shape of the flow profile, the spectral index of particles depends on the Peclet number of the shock. Results are given for some sample profiles.

Schneider, Peter↗