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.

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulation of flow and mixing in fracture intersections

Fluid transport through fractured geological formations is strongly influenced by the redistribution of solutes at fracture intersections. In this study, we perform detailed numerical simulations of flow and scalar transport within the intersection of two smooth, planar fractures. The analysis focuses on the mixing ratio, the proportion of solute flux exiting along the outlet branch aligned with the primary inlet flow direction, relative to the total solute flux at the outlets. We systematically investigate how the mixing ratio varies with four key parameters: Peclet number, Reynolds number, flow rate ratio between outlet branches, and fracture intersection angle. Results show that the mixing ratio decreases with increasing Peclet number and outlet flow rate ratio, consistent with reduced diffusive spreading and enhanced streamline routing. While low Reynolds numbers have minimal impact, inertial effects at higher Reynolds numbers significantly increase the mixing ratio. Additionally, acute and obtuse intersection angles alter flow partitioning and modify the solute distribution at the outlets. These findings provide a quantitative basis for incorporating physically realistic mixing behavior—intermediate between complete mixing and streamline-following assumptions—into network-scale transport models. The results have direct relevance to subsurface energy systems, including geothermal energy production, carbon sequestration, and contaminant remediation.

58 GEOSCIENCES↗
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↗
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↗
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↗
DOE OSTI2022

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↗
DOE OSTI2023

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↗
DOE OSTI2022

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↗
DOE OSTI2020

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

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↗
DOE OSTI2025

Unraveling the Determinant Mechanisms in Flow-Mediated Crystal Growth and Phase Behaviors

To uncover the critical mechanisms responsible for mesoscopic level development during flow-mediated crystal growth, we develop a semi-two-way hydrodynamic coupled structural phase-field crystal formalism (HXPFC-s2). The new formalism, inspired by previous attempts at coupling hydrodynamic and phase-field crystal (PFC) equations, allows for studying mesoscopic flow-mediated crystallization at diffusive timescales pertinent to industrial applications. Unlike previous efforts, the devised coupling to the structural PFC (XPFC) equations allows generalization to more complex crystal structures through explicit parameterization of the direct correlation function (DCF). Utilizing the HXPFC-s2 formalism, we seek to uncover the determinant physical mechanisms in crystallization under simple shear flows by comparing temperature-driven crystallization to flow-mediated crystallization under varying flow-strengths. Parallels and deviations of under-cooling and flow-strength effects on crystal growth are drawn using the crystal cluster-size and system ordering time evolutions. In doing so, we identify scaling behaviors with a Peclet-like number, Pe∼, a critical Peclet-like number, Pe∼*, and flow-field-crystal plane-dependent interactions. Our findings may be relevant for controlling crystal growth and phase behaviors in flow applications.

Willis, L. Connor (ORCID:0009000961321848)↗
DOE OSTI2023

Mass transfer in catalytic depolymerization: External effectiveness factors and serendipitous processivity in stagnant and stirred melts

Several heterogeneous catalysts are being developed to recycle plastics. Most operate in viscous polymer melts, where external mass transfer effects could limit the supply of co-reactants to active sites. External mass transfer can also impede the diffusion of long chain products away from the catalyst after each cut. Product egress limitations could potentially confer unintentional processivity to catalyst operation, i.e. a tendency for the catalyst to repeatedly cut the same chain after an initial encounter. We formulate reaction–diffusion equations to quantify mass transfer effects on the co-reactant transport to the catalyst and the degree of serendipitous processivity. Results are developed for catalysts in stagnant or stirred melts, with simple expressions involving Damkohler, Peclet, and Sherwood numbers, i.e. dimensionless combinations of rate constants, catalyst particle size, polymer diffusivities, and shear rates (where applicable). In conclusion, we estimate the impact of these effects for a spherical core–shell catalyst.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗