Search NASA⌕ Search

SEARCH · Search NASA

Results for “Hydrodynamic stability”

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
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Electroconvection near an ion-selective surface with Butler–Volmer kinetics

We study the effects of interfacial kinetics on the electro-hydrodynamics of ion transport near an ion-selective surface using a combination of linear stability analysis and numerical simulation. The finite kinetics of the electrolyte–electrode interface affects the ion transfer and electroconvection in many ways. On a surface of fixed topography, such as a metal surface of slow and stable ion deposition or covered by a polymer membrane, the finite kinetics reduces the current in one-dimensional ion diffusion/migration, increases the critical voltage for the onset of the electroconvective instability, changes the dynamics of the electroconvection and the overlimiting current, and enhances the lateral ion diffusion within the interfacial layer. The first three effects are indirectly caused by the reaction kinetics and can be characterized by an effective voltage difference across the liquid electrolyte. In comparison, the last effect is controlled by a direct interplay between kinetics and nonlinear electroconvection. Scaling laws for ion transport and features of electroconvection are proposed. We also analyse the linear stability of a surface which evolves under ion deposition and find that the finite kinetics decreases the growth rate of both electroconvective and morphological instabilities and therefore modifies the wavenumber of the most unstable mode.

Mechanics↗
DOE OSTI2020

Correlating inter-particle forces and particle shape to shear-induced aggregation/fragmentation and rheology for dilute anisotropic particle suspensions: A complementary study via capillary rheometry and in-situ small and ultra-small angle X-ray scattering

Predicting the stability of suspensions composed of anisotropic particles is difficult due to the complexity and interplay of hydrodynamic and colloidal forces. For example, it is not well understood how non-spherical particle shapes affect suspension-rheological responses. This study presents evidence showing how the stacking of boehmite platelets dispersed in aqueous solution leads to moderate shear thinning, even at low boehmite mass fractions. Wide angle X-ray, small angle X-ray, and ultra-small angle X-ray scattering were used to simultaneously monitor changes to the size and fractal dimensions of boehmite aggregates from 6 – 10,000 Å as the sample was recirculated through an in situ capillary rheometer. The latter also provided simultaneous measurements of the viscosity of the suspensions. Computational fluid dynamics modeling of the apparatus was used to provide a more rigorous analysis of the fluid flow. In addition to shear-induced aggregation/fragmentation correlated to a complicated balance between hydrodynamic and colloidal forces, the results indicate that orientationally-dependent interactions, arising from non-spherical boehmite aggregate sub-units, may explain the relatively large experimental viscosity when the hydrodynamic force is small compared to colloidal forces. Ultimately, the results from this study can provide a foundation for future work investigating particle shape effects on colloidal and hydrodynamics forces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

High yields in direct-drive inertial confinement fusion using thin-ice DT liner targets

In this work, we report simulation results that suggest record neutron yields above 3 x 10 14 can be achieved in laser-driven fusion on OMEGA by imploding large-outer-diameter targets (~1010μm) containing ice layers as thin as 34μm. Historically, the performance of thin and fast capsules is thought to be highly degraded by hydrodynamic instabilities. We present pulse designs that push implosion velocities above 650 km/s while maintaining enough ablative stabilization and low enough in-flight aspect ratios to survive the acceleration phase. We show these targets also possess improved stability during deceleration stemming from augmented core temperatures and low fuel convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Flow-induced oscillation patterns for two tandem cylinders with turbulence stimulation and variable stiffness and damping

Herein, the oscillation patterns of two 1-DOF cylinder oscillators, undergoing VIV and galloping, are investigated in a free surface water channel for 3.2 × 104 ≤ Re ≤ 1.2 × 10 5 . The cylinders are arranged in tandem and supported by springs for a range of different spring stiffness and damping parameters. The efficiency of a current energy converter (CEC), based on flow-induced oscillation (FIO) of multiple cylinders in tandem, is critically related to the cylinder oscillation patterns. Due to the limited number of studies on the FIO for multiple cylindrical or prismatic bodies, their oscillation patterns have neither been identified nor classified. The surfaces of the cylinders are modified by turbulence stimulation to enhance FIO. Three different center-to-center spacing, five stiffness, and six damping ratios, for a total of 90 sets of experiments were conducted. The current velocity range is from 0.34 m/s to 1.32 m/s. From more than 2000 tests, five major patterns, nine sub-patterns are identified and classified. The patterns are defined based on the amplitudes, frequencies and the phase angle differences between the two cylinders. The characteristics and mechanics of each oscillation pattern are explored and explained from the perspective of fluid-structure interaction (FSI). By systematically varying the parameters, the underlying hydrodynamic mechanisms, including the coupling level between vortices and cylinders, the shielding effect, the wake effect, and the stability states are revealed. Few preliminary observations on the connection between oscillation pattern and harnessed power by the tandem cylinders are reported.

42 ENGINEERING↗
DOE OSTI2021

MARS-Q modeling of kink-peeling instabilities in DIII-D QH-mode plasma

In quiescent H-mode (QH-mode) regime, edge harmonic oscillations (EHOs) are believed to provide necessary radial transport to prevent occurrence of large edge localized modes. A systematic modeling study is performed here on the low- n EHOs in a DIII-D QH-mode plasma (Chen et al 2016 Nucl. Fusion 56 076011), by utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). Both the n = 1 and n = 2 instabilities are found to be strongly localized near the plasma edge, exhibiting the edge-peeling characteristics. The DIII-D resistive wall is found to have minor effects on these instabilities. The plasma resistivity is found to strongly modify the mode growth rate. Assuming the Spitzer model for the plasma resistivity, the computed mode growth rate scales as S sup>-1/3 with S being the Lundquist number. Toroidal flow of the plasma slightly stabilizes these edge localized kink-peeling modes. Drift kinetic effects all have a destabilization effect on these modes. Non-perturbative magneto-hydrodynamic-kinetic hybrid computations find that the drift kinetic effects associated with thermal particle species push the peak location of the eigenmode radially inward but still in the pedestal region. The modeled plasma temperature and density fluctuations in the plasma edge region, as well as the poloidal magnetic field perturbations along both the low and high field sides of the plasma surface, are in good agreement with experimental measurements. Finally, the quasi-linear initial value simulations find a strong non-linear interplay between the kink-peeling instability and the toroidal flow near the plasma edge. The combined effect of the damping of the flow amplitude and change of the edge flow shear is found to be the stabilizing factor for the kink-peeling mode, leading to the mode saturation and thus EHOs.

Physics↗
DOE OSTI2020

The High-energy Radiation Environment around a 10 Gyr M Dwarf: Habitable at Last?

Recent work has demonstrated that high levels of X-ray and UV activity on young M dwarfs may drive rapid atmospheric escape on temperate, terrestrial planets orbiting within the habitable zone. However, secondary atmospheres on planets orbiting older, less active M dwarfs may be stable and present more promising candidates for biomarker searches. In order to evaluate the potential habitability of Earth-like planets around old, inactive M dwarfs, we present new Hubble Space Telescope and Chandra X-ray Observatory observations of Barnard’s Star (GJ 699), a 10 Gyr old M3.5 dwarf, acquired as part of the Mega-MUSCLES program. Despite the old age and long rotation period of Barnard’s Star, we observe two FUV (δ 130 ≈ 5000 s; E 130 ≈ 10 29.5 erg each) and one X-ray (E X ≈ 10 29.2 erg) flares, and we estimate a high-energy flare duty cycle (defined here as the fraction of the time the star is in a flare state) of ~25%. A publicly available 5 Å to 10 μm spectral energy distribution of GJ 699 is created and used to evaluate the atmospheric stability of a hypothetical, unmagnetized terrestrial planet in the habitable zone (r HZ ~ 0.1 au). Both thermal and nonthermal escape modeling indicate (1) the quiescent stellar XUV flux does not lead to strong atmospheric escape: atmospheric heating rates are comparable to periods of high solar activity on modern Earth, and (2) the flare environment could drive the atmosphere into a hydrodynamic loss regime at the observed flare duty cycle: sustained exposure to the flare environment of GJ 699 results in the loss of ≈87 Earth atmospheres Gyr -1 through thermal processes and ≈3 Earth atmospheres Gyr -1 through ion loss processes. These results suggest that if rocky planet atmospheres can survive the initial ~5 Gyr of high stellar activity, or if a second-generation atmosphere can be formed or acquired, the flare duty cycle may be the controlling stellar parameter for the stability of Earth-like atmospheres around old M stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Self-Assembly Interactions in Magnetite-Coated Cellulose Nanocrystals: Implications for Magnetic Hyperthermia Applications

Magnetic cellulose nanocrystal (MCNC) nanocomposites are promising sustainable and biocompatible platforms for magnetic hyperthermia; however, the molecular mechanisms governing Fe 3 O 4 adsorption and deposition onto CNCs remain poorly understood. Here, sulfated (S-CNC) and TEMPO-oxidized CNCs (T-CNC) were used to prepare nanocomposites at 1:2 and 1:4 CNC:Fe 3 O 4 mass ratios, enabling a systematic evaluation of how surface chemistry and nanoparticle loading dictate interfacial interactions and magneto-colloidal behavior. Bare magnetite nanoparticles were 21 ± 5 nm by TEM but grew to 144 ± 18 in the DLS measurement at pH 7. The S-CNC nanocomposites had hydrodynamic sizes between 144 and 210 nm, not much larger than the 140 nm long CNC rods, suggesting an enhanced dispersion stability compared to Fe3O4 alone. X-ray photoelectron spectroscopy combined with density functional theory revealed that −OH and −COOH groups drive electrostatic adsorption with charge transfer from Fe 3 O 4 to the CNC surface, while T-CNCs showed more favorable adsorption energies and evidence of covalent Fe−O bonding. Vibrating sample magnetometry demonstrated superparamagnetic behavior for all samples, with S-CNC/Fe 3 O 4 1:4 and 1:2 displaying saturation magnetizations of 78 and 77 emu/g-Fe 3 O 4 , close to the 83 emu/g of bare magnetite. The T-CNC composites showed lower (60 and 66 emu/g-Fe3O4) saturation magnetizations. Zero-field-cooled/field-cooled measurements resulted in a blocking temperature of 112 K for all samples, except T-CNC/Fe 3 O 4 1:2 (100 K). Magnetic hyperthermia studies revealed that specific absorption rate (SAR) increased with field strength and Fe 3 O 4 content; however, S-CNC/Fe 3 O 4 (1:2) achieved the highest intrinsic SAR per gram of Fe 3 O 4 (649 W/g-Fe 3 O 4 ) likely due to its anisotropy and fast magnetic relaxation. Cytotoxicity assays confirmed that all nanocomposites were nontoxic toward mammalian cells. These results establish quantitative structure−property relationships between CNC surface chemistry, interfacial bonding mechanisms, and magnetic heating performance, providing a foundation for rational design of biocompatible magnetic nanocomposites for hyperthermia and related applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Measurement of mix at the fuel–ablator interface in indirectly driven capsule implosions on the National Ignition Facility

The interface between the capsule ablator and fuel ice layer is susceptible to hydrodynamic instabilities. The subsequent mixing of hot ablator material into the ice reduces fuel compression at stagnation and is a candidate for reduced capsule performance. The ability to diagnose ice–ablator mix is critical to understanding and improving stability at this interface. Combining the crystal backlighter imager with the single line of sight camera on the National Ignition Facility (NIF) allows direct measurement of ice–ablator mix by providing multiple quasi-monochromatic radiographs of layered capsule implosions per experiment with high spatial (∼12 μm) and temporal (∼35 ps) resolution. The narrow bandwidth of this diagnostic platform allows radiography of the inner edge of the capsule limb close to stagnation without capsule self-emission contaminating the data and removes opacity uncertainties typically associated with the spectral content of the radiograph. Analysis of radiographic data via a parameterized forward-fitting Abel inversion technique provides measurements of the distribution of mix mass inwards from the ice–ablator interface. The sensitivity of this mix measurement technique was demonstrated by applying it to layered experiments in which the stability of the ice–ablator interface was expected to vary significantly. Additional experiments suggest that high-density carbon capsules that employ a buried-layer dopant profile suffer from mixing at the innermost doped–undoped interface. Data from these experiments suggest that opacity models used in hydrodynamic simulations of NIF experiments can potentially over-predict the opacity of doped capsules. LLNL-JRNL-850535-DRAFT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Optical N-plasmon: topological hydrodynamic excitations in graphene from repulsive Hall viscosity

Abstract Edge states occurring in Chern and quantum spin-Hall phases are signatures of the topological electronic band structure in two-dimensional (2D) materials. Recently, a new topological electromagnetic phase of graphene characterized by the optical N-invariant was proposed. Optical N-invariant arises from repulsive Hall viscosity in hydrodynamic many-body electron systems, distinct from the Chern and Z 2 invariants. In this paper, we introduce the topologically protected edge excitation—optical N-plasmon of interacting many-body electron systems in the topological optical N-phase. These optical N-plasmons are signatures of the topological plasmonic band structure in 2D materials. We demonstrate that optical N-plasmons exhibit unique dispersion relations, stability against various boundary conditions, and edge profiles when compared with the topologically trivial edge magneto plasmons. Based on the optical N-plasmon, we design an ultra sub-wavelength broadband topological hydrodynamic circulator, which is a chiral quantum radio-frequency circuit component crucial for information routing and interfacing quantum–classical computing systems. Furthermore, we reveal that optical N-plasmons can be effectively tuned by the neighboring dielectric environment without breaking the topological properties. Our work provides a smoking gun signature of topological electromagnetic phases occurring in 2D materials arising from repulsive Hall viscosity.

Physics↗
DOE OSTI2024

Advancing our understanding of instabilities in high-energy-density systems through the study of benchmark datasets with advanced modeling tools

Numerous types of pulsed power driven inertial confinement fusion (ICF) and high energy density (HED) systems rely on implosion stability to achieve desired temperatures, pressures, and densities. Sandia National Laboratories Pulsed Power Sciences Center’s main ICF platform, Magnetized Liner Inertial Fusion (MagLIF), suffers from implosion instabilities which limit attainable fuel conditions and can compromise fuel confinement. This Truman Fellowship research primarily focused on computationally exploring (a) methods for improving our understanding of hydrodynamic and magnetohydrodynamic instabilities that form during cylindrical liner implosions, (b) methods for mitigating implosion instabilities, particularly those that degrade performance of MagLIF targets, and (c) novel MagLIF target designs intended to improve target performance primarily via enhanced implosion stability. Several multi-dimensional computational tools were used, including the magnetohydrodynamics code ALEGRA, the radiation-magnetohydrodynamics code HYDRA, and the magnetohydrodynamics code KRAKEN. This research succeeded in executing and analyzing simulations of automagnetizing liner implosions, shockless MagLIF implosions, dynamic screw pinch driven cylindrical liner implosions, and cylindrically convergent HED instability studies. The methods and tools explored and developed in this Truman Fellowship research have been published in several peer-reviewed journal articles and will serve as useful contributions to the fields of pulsed power science and engineering, particularly pertaining to pulsed power ICF and HED science.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mitigation of deceleration-phase Rayleigh–Taylor instability growth in inertial confinement fusion implosions

Rayleigh–Taylor growth during shell deceleration is one of the main limiting factors for target performance in inertial confinement fusion implosions. Here, using analytical scaling laws and hydrodynamic simulations, we show that such amplification can be mitigated by reducing the initial mass density in the central target region. The perturbation growth reduction is caused by a smaller hot-spot convergence ratio during deceleration, increased density scale length, and enhanced ablation stabilization. The required central density reduction can be achieved using the dynamic shell formation concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Predicting the Slowing of Stellar Differential Rotation by Instability-driven Turbulence

Abstract Differentially rotating stars and planets transport angular momentum (AM) internally due to turbulence at rates that have long been a challenge to predict reliably. We develop a self-consistent saturation theory, using a statistical closure approximation, for hydrodynamic turbulence driven by the axisymmetric Goldreich–Schubert–Fricke instability at the stellar equator with radial differential rotation. This instability arises when fast thermal diffusion eliminates the stabilizing effects of buoyancy forces in a system where a stabilizing entropy gradient dominates over the destabilizing AM gradient. Our turbulence closure invokes a dominant three-wave coupling between pairs of linearly unstable eigenmodes and a near-zero frequency, viscously damped eigenmode that features latitudinal jets. We derive turbulent transport rates of momentum and heat and provide them in analytic forms. Such formulae, free of tunable model parameters, are tested against direct numerical simulations; the comparison shows good agreement. They improve upon prior quasi-linear or “parasitic saturation” models containing a free parameter. Given model correspondences, we also extend this theory to heat and compositional transport for axisymmetric thermohaline-instability-driven turbulence in certain regimes.

Astronomy & Astrophysics↗
DOE OSTI2022

Comprehensive Analysis of the Neutrino Process in Core-collapsing Supernovae

We investigate the neutrino flavor change effects due to neutrino self-interaction and shock wave propagation, as well as the matter effects on the neutrino process in core-collapsing supernovae (CCSNe). For the hydrodynamics, we use two models: a simple thermal bomb model and a specified hydrodynamics model for SN1987A. For the presupernova model, we take an updated model, adjusted to explain SN1987A, which employs recent developments in the (n, γ) reaction rates for nuclei near the stability line (A ~ 100). As for the neutrino luminosity, we adopt two different models: equivalent neutrino luminosity and nonequivalent luminosity models. The latter is taken from a synthetic analysis of CCSN simulation data, which quantitatively presented the results obtained by various neutrino transport models. Relevant neutrino-induced reaction rates are calculated using a shell model for light nuclei and a quasiparticle random phase approximation model for heavy nuclei. For each model, we present abundances of the light nuclei ( 7 Li, 7 Be, 11 B, and 11 C) and the heavy nuclei ( 92 Nb, 98 Tc, 138 La, and 180 Ta) produced by the neutrino process. The light nuclei abundances turn out to be sensitive to the Mikheyev–Smirnov–Wolfenstein (MSW) region around O-Ne-Mg layer while the heavy nuclei are mainly produced prior to the MSW region. Through detailed analyses, we find that neutrino self-interaction becomes a key ingredient, in addition to the MSW effect, for understanding the neutrino process and the relevant nuclear abundances. The normal mass hierarchy is shown to be more compatible with the meteorite data. The main nuclear reactions for each nucleus are also investigated in detail.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Closed-loop control of active nematic flows

Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware/software/wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.

Nishiyama, Katsu [Brandeis Univ., Waltham, MA (Uni↗
DOE OSTI2025

Closed-Loop Control of Active Nematic Flows

Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware-software-wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.

Active nematics↗
DOE OSTI2022

Experimental and analytical study of the hydrodynamic and single and two-phase convective heat transfer performance of flexible PDMS microchannels with micropillar arrays

Various copper and silicon based thermal management systems are used in the cooling of electronics. However, the rigid nature of these materials along with their high thermal and electrical conductivity pose a difficulty in developing direct contact embedded flexible cooling systems that can offer robust cooling performance. The low density, thermal stability, chemical inertness, and electrical insulation of Polydimethylsiloxane (PDMS) make it an ideal material to develop lightweight direct contact thermal management systems for electronics. Its ease of fabrication with tunable flexibility provides the opportunity to go beyond traditional electronics and develop advanced active and passive thermal management systems for a wide–range of applications in foldable and wearable electronics, liquid cooling garments, microgravity, and electric motors. In this study, a flexible PDMS based microchannel with micropillar arrays, which enhance the thermal performance of the device through capillary-assisted flow, has been developed. The hydrodynamic and convection heat transfer performance of three PDMS wick pillar geometries, ranging from a porosity of 0.8–0.91, are investigated and compared under single-phase and two-phase conditions. Dielectric coolant FC-3283 is employed and permeability measurements are made for mass fluxes ranging from 53 kg/m 2 s to 369 kg/m 2 s. Given its conformability, the device demonstrates a deviation from Darcy’s Law, within the laminar regime, with an increasing permeability with mass flux at the rate of ~0.5–0.8 Darcy/(kg/m 2 s). A semi-analytical model has been developed and reported to quantify the conformability of the device. The heat transfer performance is experimentally evaluated using the same dielectric fluid for mass flux ranging from 105 kg/m 2 s to 420 kg/m 2 s with heat fluxes ranging from 1.5 W/cm 2 to 16 W/cm 2 . Heat transfer coefficients of up to 7000 W/m 2 K are observed, which are comparable to copper and silicon microchannels. The effect of porosity on the single phase thermal performance has been evaluated against the pumping power to provide a basis for thermal management system design. Finally, high-speed imaging is performed to study the two-phase flow characteristics to provide insight into the vapor formation and removal.

42 ENGINEERING↗