Influence of Mg2+ on CaCO3 precipitation during subsurface reactive transport in a homogeneous silicon-etched pore network
Explore the source record for details and available documents.
nuclear (including radiation effects), carbon sequ↗
SEARCH · Engineering Papers
Results for “Reactive Flow”
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.
Modeling of fault activation and seismicity by injection directly into a fault zone associated with hydraulic fracturing of shale-gas reservoirs
We conducted three-dimensional coupled fluid-flow and geomechanical modeling of fault activation and seismicity associated with hydraulic fracturing stimulation of a shale-gas reservoir. We simulated a case in which a horizontal injection well intersects a steeply dip- ping fault, with hydraulic fracturing channeled within the fault, during a 3-hour hydraulic fracturing stage. Consistent with field observations, the simulation results show that shale-gas hydraulic fracturing along faults does not likely induce seismic events that could be felt on the ground surface, but rather results in numerous small microseismic events, as well as aseismic deformations along with the fracture propagation. The calculated seismic moment magnitudes ranged from about -2.0 to 0.5, except for one case assuming a very brittle fault with low residual shear strength, for which the magnitude was 2.3, an event that would likely go unnoticed or might be barely felt by humans at its epicenter. The calculated moment magnitudes showed a dependency on injection depth and fault dip. We attribute such dependency to variation in shear stress on the fault plane and associated variation in stress drop upon reactivation. Our simulations showed that at the end of the 3-hour injection, the rupture zone associated with tensile and shear failure extended to a maximum radius of about 200 m from the injection well. The results of this modeling study for steeply dipping faults at 1000 to 2500 m depth is in agreement with earlier studies and field observations showing that it is very unlikely that activation of a fault by shale-gas hydraulic fracturing at great depth (thousands of meters) could cause felt seismicity or create a new flow path (through fault rupture) that could reach shallow groundwater resources.
54 ENVIRONMENTAL SCIENCES↗
Numerical investigation of spontaneous flame propagation under RCCI conditions
This paper presents results from one and two-dimensional direct numerical simulations under Reactivity Controlled Compression Ignition (RCCI) conditions of a primary reference fuel (PRF) mixture consisting of n-heptane and iso-octane. RCCI uses in-cylinder blending of two fuels with different autoignition characteristics to control combustion phasing and the rate of heat release. These simulations employ an improved model of compression heating through mass source/sink terms developed in a previous work by Bhagatwala et al. (2014), which incorporates feedback from the flow to follow a predetermined experimental pressure trace. Two-dimensional simulations explored parametric variations with respect to temperature stratification, pressure profiles and n-heptane concentration. Furthermore, statistics derived from analysis of diffusion/reaction balances locally normal to the flame surface were used to elucidate combustion characteristics for the different cases. Both deflagration and spontaneous ignition fronts were observed to co-exist, however it was found that higher n-heptane concentration provided a greater degree of flame propagation, whereas lower n-heptane concentration (higher fraction of iso-octane) resulted in more spontaneous ignition fronts. A significant finding was that simulations initialized with a uniform initial temperature and a stratified n-heptane concentration field, resulted in a large fraction of combustion occurring through flame propagation. The proportion of spontaneous ignition fronts increased at higher pressures due to shorter ignition delay when other factors were held constant. For the same pressure and fuel concentration, the contribution of flame propagation to the overall combustion was found to depend on the level of thermal stratification, with higher initial temperature gradients resulting in more deflagration and lower gradients generating more ignition fronts. Statistics of ignition delay are computed to assess the Zel’dovich (1980) theory for the mode of combustion propagation based on ignition delay gradients.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Plutonium reactive transport in fractured granite: Multi-species experiments and simulations
Plutonium (Pu) in the subsurface environment can transport in different oxidation states as an aqueous solute or as colloidal particles. The transport behavior of Pu is affected by the relative abundances of these species and can be difficult to predict when they simultaneously exist. This study investigates the concurrent transport of Pu intrinsic colloids, Pu(IV) (aq) and Pu(V-VI) (aq) through a combination of controlled experiments and semi-analytical dual-porosity transport modeling. Pu transport experiments were conducted in a fractured granite at high and low flow rates to elucidate sorption processes and their scaling behavior. In the experiments, Pu(IV) (aq) was the least mobile of the Pu species, Pu(V-VI) (aq) had intermediate mobility, and the colloidal Pu, which consisted mainly of precipitated and/or hydrolyzed Pu(IV), was the most mobile. The semi-analytical modeling revealed that the sorption of each Pu species was rate-limited, as the sorption could not be described by assuming local equilibrium in the experiments. The model was able to describe the sorption of the different Pu species that occurring either on fracture surfaces, in the pores of the rock matrix, or simultaneously in both locations. While equally good fits to the data could be achieved using any of these assumptions, a fracture-dominated process was considered to be the most plausible because it provided the most reasonable estimates of sorption rate constants. Importantly, a key result of this work is that the sorption rate constant of all Pu species tends to decrease with increasing time scales, which implies that Pu will tend to be more mobile at longer time scales than observations at shorter time scales suggest. This result has important implications for predicting the environmental impacts of Pu in the safety assessments of geologic repositories for radioactive waste disposal, and we explore potential mechanistic bases for upscaling the sorption rate constants to time and distance scales that cannot be practically evaluated in experiments.
58 GEOSCIENCES↗
On velocity and reactive scalar spectra in turbulent premixed flames
Abstract Kinetic energy and reactive scalar spectra in turbulent premixed flames are studied from compressible three-dimensional direct numerical simulations (DNS) of a temporally evolving rectangular slot-jet premixed flame, a statistically one-dimensional configuration. The flames correlate to a lean premixed hydrogen–air mixture at an equivalence ratio of 0.7, preheated to 700 K and at 1 atm, and three DNS are considered with a fixed jet Reynolds number of 10 000 and a jet Damköhler number varying between 0.13 and 0.54. For the study of spectra, motivated by the need to account for density change, which can be locally strong in premixed flames, a new density-weighted definition for two-point velocity/scalar correlations is proposed. The density-weighted two-point correlation tensor retains the essential properties of its constant-density (incompressible) counterpart and recovers the density-weighted Reynolds stress tensor in the limit of zero separation. The density weighting also allows the derivation of balance equations for velocity and scalar spectrum functions in the wavenumber space that illuminate physics unique to combusting flows. Pressure–dilatation correlation is a source of kinetic energy at high wavenumbers and, analogously, reaction rate–scalar fluctuation correlation is a high-wavenumber source of scalar energy. These results are verified by the spectra constructed from the DNS data. The kinetic energy spectra show a distinct inertial range with a $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}-5/3$ scaling followed by a ‘diffusive–reactive’ range at higher wavenumbers. The exponential drop-off in this range shows a distinct inflection in the vicinity of the wavenumber corresponding to a laminar flame thickness, $\delta _L$ , and this is attributed to the contribution from the pressure–dilatation term in the energy balance in wavenumber space. Similarly, a clear spike in spectra of major reactant species (hydrogen) arising from the reaction-rate term is observed at wavenumbers close to $\delta _L$ . It appears that in the inertial range classical scaling laws for the spectra involving the Kolmogorov scale are applicable, but in the high-wavenumber range where chemical reactions have a strong signature the laminar flame thickness produces a better collapse. It is indicated that a full scaling should perhaps involve the Kolmogorov scale, laminar flame thickness, Damköhler number and Karlovitz number.
42 ENGINEERING↗
Large eddy simulation of particle dispersion analysis in a pressurized reactive fluid–particle system
Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new combustion technologies. Because of high pressure, the volumetric fraction of particles in a pressurized oxy-combustor could be higher than that in the conventional combustor, consequently, the particle dynamics in the pressurized vessel need to be evaluated accordingly. Specifically, the present computational work analyzes the pattern of the particle dispersion by extracting the particle mass flow rate and the particle size distribution in the pressurized combustor. A large eddy simulation (LES) with the environment pressure 15-bar and 100 kWth thermal input of fuel particle has been performed. The particle mass flow rates sampled in several cross-section planes strongly fluctuate, indicating the impacts of the flame as well as the gas-phase flow on the particle movement. Animations of the particle-gas flow show that the flow recirculation plays a major role in producing the pulsing phenomenon of the particle mass flow rate.
Li, Lei↗
Experimental and modeling analysis of detonation in circular arcs of the conventional high explosive PBX 9501
We examine the diffraction dynamics of a two-dimensional (2D) detonation in a circular arc of the conventional HMX-based, high performance, solid explosive PBX 9501, for which the detonation reaction zone length scale is estimated to be of the order of 100–150 µm. In this configuration, a steady propagating detonation will develop, sweeping around the arc with constant angular speed. We report on results from three PBX 9501 arc experiments, exploring the variation in linear speed on the inner and outer arc surfaces for the steady wave along with the structure of the curved detonation front, as a function of varying inner surface radius and arc thickness. Comparisons of the properties of the motion of the steady wave for each arc configuration are then made with a spatially-distributed PBX 9501 reactive burn model, calibrated to detonation performance properties in a 2D planar slab geometry. We show that geometry-induced curvature of the detonation near the inner arc surface has a significant effect on the detonation motion even for conventional high explosives. We also examine the detonation driving zone structure for each arc case, and thus the subsonic regions of the flow that determine the influence of the arc geometry on the detonation propagation. In addition, streamline paths and reaction progress isolines are calculated. Overall, we conclude that a common approximation for modeling conventional high explosive detonation, wherein the shock-normal detonation speed is assumed equal to the Chapman–Jouguet speed, can lead to significant errors in describing the speed at which the detonation propagates.
42 ENGINEERING↗
Macroporous Carbon Supported Zerovalent Iron for Remediation of Trichloroethylene
Groundwater contamination with chlorinated hydrocarbons has become a widespread problem that threatens water quality and human health. Permeable reactive barriers (PRBs), which employ zerovalent iron, are effective for remediation; however, a need exists to reduce the economic and environmental costs associated with constructing PRBs. Here, we present a method to produce zerovalent iron supported on macroporous carbon using only lignin and magnetite. Biochar- ZVI (BC-ZVI) produced by this method exhibits a broad pore size distribution with micrometer sized ZVI phases dispersed throughout a carbon matrix. X-ray diffraction revealed that pyrolysis at 900 °C of a 50/50 lignin$-$magnetite mixture resulted in almost complete reduction of magnetite to ZVI and that compression molding promotes iron reduction in pyrolysis due to mixing of starting materials. High temperature pyrolysis of lignin yields some graphite in BC-ZVI due to reduction of carbonaceous gases on iron oxides. TCE was removed from water as it passed through a column packed with BC-ZVI at flow rates representative of average and high groundwater flow. Lastly, one-dimensional convection$-$dispersion modeling revealed that adsorption by biochar influences TCE transport and that BC-ZVI facilitated removal of TCE from contaminated water by both adsorption and degradation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The near-infrared spectrum of ethynyl radical
Transient diode laser absorption spectroscopy has been used to measure three strong vibronic bands in the near infrared spectrum of the C{sub 2}H, ethynyl, radical not previously observed in the gas phase. The radical was produced by ultraviolet excimer laser photolysis of either acetylene or (1,1,1)-trifluoropropyne in a slowly flowing sample of the precursor diluted in inert gas, and the spectral resolution was Doppler-limited. The character of the upper states was determined from the rotational and fine structure in the observed spectra and assigned by measurement of ground state rotational combination differences. The upper states include a {sup 2}Σ{sup +} state at 6696 cm{sup −1}, a second {sup 2}Σ{sup +} state at 7088 cm{sup −1}, and a {sup 2}Π state at 7110 cm{sup −1}. By comparison with published calculations [R. Tarroni and S. Carter, J. Chem. Phys 119, 12878 (2003); Mol. Phys. 102, 2167 (2004)], the vibronic character of these levels was also assigned. The observed states contain both X{sup 2}Σ{sup +} and A{sup 2}Π electronic characters. Several local rotational level perturbations were observed in the excited states. Kinetic measurements of the time-evolution of the ground state populations following collisional relaxation and reactive loss of the radicals formed in a hot, non-thermal, population distribution were made using some of the strong rotational lines observed. The case of C{sub 2}H may be a good place to investigate the behavior at intermediate pressures of inert colliders, where the competition between relaxation and reaction can be tuned and observed to compare with master equation models, rather than deliberately suppressed to measure thermal rate constants.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Modeling the Transient Reactor Test Loop in TRACE for PWR Critical Heat Flux Benchmarking
Understanding transient critical heat flux (CHF) resulting from a Reactivity Insertion Accident (RIA) in nuclear systems is an extremely undeveloped field. Though many experiments have been done for CHF, remarkably few have been done for rods at PWR operating temperature and pressure (300 deg. C, 2200 psia) with a short time constant ramp in power. This leads to many concerns regarding understanding how RIA scenarios develop and evolve in commercial nuclear reactors. With a limited understanding of the phenomenology governing transient CHF, there is little hope in being able to accurately model these scenarios using today's thermal hydraulic codes. The Accident Tolerant Fuels (ATF) Program established by the Department of Energy (DOE) is seeking out new fuel and cladding systems that can be used to reduce the probability and consequences of severe accidents. To better understand this phenomenon, the DOE is restarting the Transient Reactor Test (TREAT) Facility at the Idaho National Laboratory (INL). TREAT will be able to provide representative neutronic boundary conditions for testing coupled thermal hydraulics on the Transient Water Environment Recirculating Loop (TWERL) to study transient CHF with the best representation possible. To support this task, Oregon State University (OSU) is designing and building a representative flow loop that will be able to mimic the conditions of TWERL (sans the neutronic boundary conditions) to create a benchmark study for existing lumped parameter thermal-hydraulic codes such as RELAP5 and TRACE. This flow loop, the Transient Reactor Test Facility (TRTL), will support the ATF program by providing much needed NQA-1 data for transient CHF at operating temperature and pressure.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Enabling BWR fuel rod analysis in the BISON fuel performance code
Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach
In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effects on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Biogeochemical Transport and Reaction Model (BeTR) v1
The Biogeochemical Transport and Reaction Model (BeTR) is a F90 code that enables reactive transport modeling in land modules of earth system models (e.g. CESM, ACME). The code adopts the Objective-Oriented-Design, and allows users to plug in their own biogeochemical (BGC) formulations/codes, and compare them to other existing BGC codes in those ESMs. The code takes information of soil physics variables, such as variables, such as temperature, moisture, soil density profile; water flow, etc., from a land model to track the movement of different chemicals in presence of biogeochemical reactions.
Tang, Jinyun↗
Drying model of a high salt content cementitious waste form: Effect of capillary forces and salt solution
Highlights: • Drying model for a high salt content cementitious waste form is developed. • Water vapor diffusion and capillary liquid flow are distinguished. • Capillary and salt solution effects are considered in vapor-liquid equilibrium. A water transport model coupling capillary liquid flow with vapor diffusion is developed to describe the drying process for a cementitious waste form with high salinity porewater. Vapor-liquid equilibrium is formulated as the driving force for vapor diffusion and the model accounts for pore capillary and high salinity effects on water thermodynamic activity. Pore filling and porewater surface tension as a function of pore size distribution and water saturation have been quantified for the material. Geochemical speciation modeling is used to simulate porewater activity as a function of composition over the range of saturation. The theoretical relationship between relative humidity and water saturation generally agrees with experimental measurement, and the developed model is capable of predicting drying rates under various external relative humidity conditions. The model was developed to be incorporated into reactive transport models considering the effects of drying such as salt redistribution and efflorescence.
36 MATERIALS SCIENCE↗
Predicting Water Cycle Characteristics from Percolation Theory and Observational Data
The fate of water and water-soluble toxic wastes in the subsurface is of high importance for many scientific and practical applications. Although solute transport is proportional to water flow rates, theoretical and experimental studies show that heavy-tailed (power-law) solute transport distribution can cause chemical transport retardation, prolonging clean-up time-scales greatly. However, no consensus exists as to the physical basis of such transport laws. In percolation theory, the scaling behavior of such transport rarely relates to specific medium characteristics, but strongly to the dimensionality of the connectivity of the flow paths (for example, two- or three-dimensional, as in fractured-porous media or heterogeneous sediments), as well as to the saturation characteristics (i.e., wetting, drying, and entrapped air). In accordance with the proposed relevance of percolation models of solute transport to environmental clean-up, these predictions also prove relevant to transport-limited chemical weathering and soil formation, where the heavy-tailed distributions slow chemical weathering over time. The predictions of percolation theory have been tested in laboratory and field experiments on reactive solute transport, chemical weathering, and soil formation and found accurate. Recently, this theoretical framework has also been applied to the water partitioning at the Earth’s surface between evapotranspiration, ET, and run-off, Q, known as the water balance. A well-known phenomenological model by Budyko addressed the relationship between the ratio of the actual evapotranspiration (ET) and precipitation, ET/P, versus the aridity index, ET 0 /P, with P being the precipitation and ET 0 being the potential evapotranspiration. Existing work was able to predict the global fractions of P represented by Q and ET through an optimization of plant productivity, in which downward water fluxes affect soil depth, and upward fluxes plant growth. In the present work, based likewise on the concepts of percolation theory, we extend Budyko’s model, and address the partitioning of run-off Q into its surface and subsurface components, as well as the contribution of interception to ET. Using various published data sources on the magnitudes of interception and information regarding the partitioning of Q, we address the variability in ET resulting from these processes. The global success of this prediction demonstrated here provides additional support for the universal applicability of percolation theory for solute transport as well as guidance in predicting the component of subsurface run-off, important for predicting natural flow rates through contaminated aquifers.
54 ENVIRONMENTAL SCIENCES↗
Fenofibrate down-regulates the expressions of androgen receptor (AR) and AR target genes and induces oxidative stress in the prostate cancer cell line LNCaP
Highlights: ► Fenofibrate induces cell cycle arrest in G1 phase and apoptosis in LNCaP cells. ► Fenofibrate reduces the expressions of androgen receptor in LNCaP cells. ► Fenofibrate induces oxidative stress in the prostate cancer cell line LNCaP. -- Abstract: Fenofibrate, a peroxisome proliferator-androgen receptor-alpha agonist, is widely used in treating different forms of hyperlipidemia and hypercholesterolemia. Recent reports have indicated that fenofibrate exerts anti-proliferative and pro-apoptotic properties. This study aims to investigate the effects of fenofibrate on the prostate cancer (PCa) cell line LNCaP. The effects of fenofibrate on LNCaP cells were evaluated by flow cytometry, reverse transcription-polymerase chain reaction, enzyme-linked immunosorbent assays, Western blot analysis, and dual-luciferase reporter assay. Fenofibrate induces cell cycle arrest in G1 phase and apoptosis in LNCaP cells, reduces the expressions of androgen receptor (AR) and AR target genes (prostate-specific antigen and TMPRSS2), and inhibits Akt phosphorylation. Fenofibrate can induce the accumulation of intracellular reactive oxygen species and malondialdehyde, and decrease the activities of total anti-oxidant and superoxide dismutase in LNCaP cells. Fenofibrate exerts an anti-proliferative property by inhibiting the expression of AR and induces apoptosis by causing oxidative stress. Therefore, our data suggest fenofibrate may have beneficial effects in fenofibrate users by preventing prostate cancer growth through inhibition of androgen activation and expression.
60 APPLIED LIFE SCIENCES↗
Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis
Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).
42 - ENGINEERING↗
Asymmetric Wettability Directs Leidenfrost Droplets
Leidenfrost phenomena on nano- and microstructured surfaces are of great importance for increasing control over heat transfer in high power density systems utilizing boiling phenomena. They also provide an elegant means to direct droplet motion in a variety of recently emerging fluidic systems. Here, we report the fabrication and characterization of tilted nanopillar arrays (TNPAs) that exhibit directional Leidenfrost water droplets under dynamic conditions, namely on impact with Weber numbers 40 at T 325 C. The batch fabrication of the TNPAs was achieved by glancing-angle anisotropic reactive ion etching of a thermally dewet platinum mask, with mean pillar diameters of 100 nm and heights of 200-500 nm. In contrast to previously implemented macro- and microscopic Leidenfrost ratchets, our TNPAs induce no preferential directional movement of Leidenfrost droplets under conditions approaching steady-state film boiling, suggesting that the observed droplet directionality is not a result of asymmetric vapor flow. Using high-speed imaging, phase diagrams were constructed for the boiling behavior upon impact for droplets falling onto TNPAs, straight nanopillar arrays, and smooth silicon surfaces. The asymmetric impact and directional trajectory of droplets was exclusive to the TNPAs for impacts corresponding to the transition boiling regime, revealing that asymmetric wettability upon impact is the mechanism for the droplet directionality.
asymmetric rebound↗