Fluid Instabilities Research Articles

Graphene Microflower by Photothermal Marangoni-Induced Fluid Instability for Omnidirectional Broadband Photothermal Conversion.

2-D carbon-based materials are well-known for their broadband absorption properties for efficient solar energy conversion. However, their high reflectivity poses a challenge for achieving efficient omnidirectional light absorption. Inspired by the multilevel structures of the flower, a Graphene Microflower (GM) material with gradient refractive index surface was fabricated on polymer substrates using the UV-intense laser-induced phase explosion technique under the synergistic design of the photothermal Marangoni effect and the fluid instability principle. The refractive index gradient reduces light reflection and absorbs at least 96% of light at incident angles of 0-60° across the entire solar wavelength range (200-2500 nm). Over 90% absorption even at 75° angle of incidence. The light absorption is enhanced by the multiple interferometric phase cancelation and localized surface plasmon resonance, resulting in a steady-state temperature 60 °C higher than ambient conditions under one solar irradiation. The max rate of temperature rise can reach up to 62 °C s-1. The device is then integrated at the hot end of the temperature difference generator at high altitude to ensure continuous and efficient power generation, producing a steady-state power of 196 mW.

Dyakonov-Shur instability of the electronic fluid: Spectral effect of weak magnetic field

Dyakonov-Shur instability of the electronic fluid: Spectral effect of weak magnetic field

Double Layer SiO2-Coated Water-Stable Halide Perovskite as a Promising Antimicrobial Photocatalyst under Visible Light.

Halide perovskite nanocrystals (HPNCs) have emerged as promising materials for various light harvesting applications due to their exceptional optical and electronic properties. However, their inherent instability in water and biological fluids has limited their use as photocatalysts in the aqueous phase. In this study, we present highly water-stable SiO2-coated HPNCs as efficient photocatalysts for antimicrobial applications. The double SiO2 layer coating method confers long-term structural and optical stability to HPNCs in water, while the in situ synthesis of lead- and bismuth-based perovskite NCs into the SiO2 shell enhances their versatility and tunability. We demonstrate that the substantial generation of singlet oxygen via energy transfer from HPNCs enables efficient photoinduced antibacterial efficacy under aqueous conditions. More than 90% of Escherichia coli was inactivated under mild visible light irradiation for 6 h. The excellent photocatalytic antibacterial performance suggests that SiO2-coated HPNCs hold great potential for various aqueous phase photocatalytic applications.

Crystal Patterning from Aqueous Solutions via Solutal Instabilities.

Fluid instabilities can be harnessed for facile self-assembly of patterned structures on the nano- and microscale. Evaporative self-assembly from drops is one simple technique that enables a range of patterning behaviors due to the multitude of fluid instabilities that arise due to the simultaneous existence of temperature and solutal gradients. However, the method suffers from limited controllability over patterns that can arise and their morphology. Here, we demonstrate that a range of distinct crystalline patterns including hexagonal arrays, branches, and sawtooth structures emerge from evaporation of water drops containing calcium sulfate on hydrophilic and superhydrophilic substrates. Different pattern regimes emerge as a function of contact line dynamics and evaporation rates, which dictate which fluid instabilities are most likely to emerge. The underlying physical mechanisms behind instability for controlled self-assembly involve Marangoni flows and forced wetting/dewetting. We also demonstrate that these patterns composed of water-soluble inorganic crystals can serve as sustainable and easily removable masks for applications in microscale fabrication.

Giant Decrease in Interfacial Energy of Liquid Metals by Native Oxides.

Native oxides form on the surface of many metals. Here, using gallium-based liquid metal alloys, Johnson-Kendall-Roberts (JKR) measurements are employed to show that native oxide dramatically lower the tension of the metal interface from 724 to 10 mNm-1. Like conventional surfactants, the oxide has asymmetry between the composition of its internal and external interfaces. Yet, in comparison to conventional surfactants, oxides are an order of magnitude more effective at lowering tension and do not need to be added externally to the liquid (i.e., oxides form naturally on metals). This surfactant-like asymmetry explains the adhesion of oxide-coated metals to surfaces. The resulting low interfacial energy between the metal and the interior of the oxide helps stabilize non-spherical liquid metal structures. In addition, at small enough macroscopic contact angles, the finite tension of the liquid within the oxide can drive fluid instabilities that are useful for separating the oxide from the metal to form oxide-encased bubbles or deposit thin oxide films (1-5nm) on surfaces. Since oxides form on many metals, this work can have implications for a wide range of metals and metal oxides in addition to explaining the physical behavior of liquid metal.

Lateral migration of a deformable fluid particle in a square channel flow of viscoelastic fluid

Cross-stream migration of a deformable fluid particle is investigated computationally in a pressure-driven channel flow of a viscoelastic fluid via interface-resolved simulations. Flow equations are solved fully coupled with the Giesekus model equations using an Eulerian–Lagrangian method and extensive simulations are performed for a wide range of flow parameters to reveal the effects of particle deformability, fluid elasticity, shear thinning and fluid inertia on the particle migration dynamics. Migration rate of a deformable particle is found to be much higher than that of a solid particle under similar flow conditions mainly due to the free-slip condition on its surface. It is observed that the direction of particle migration can be altered by varying shear thinning of the ambient fluid. With a strong shear thinning, the particle migrates towards the wall while it migrates towards the channel centre in a purely elastic fluid without shear thinning. An onset of elastic flow instability is observed beyond a critical Weissenberg number, which in turn causes a path instability even for a nearly spherical particle. An inertial path instability is also observed once particle deformation exceeds a critical value. Shear thinning is found to be suppressing the path instability in a viscoelastic fluid with a high polymer concentration whereas it reverses its role and promotes path instability in a dilute polymer solution. It is found that migration of a deformable particle towards the wall induces a secondary flow with a velocity that is approximately an order of magnitude higher than the one induced by a solid particle under similar flow conditions.

Electro-osmotic flow instability of viscoelastic fluids in a nanochannel

The study of the complex rheological properties of viscoelastic fluids in nanochannels will facilitate the application of nanofluidics in biomedical and other fields. However, the flow of viscoelastic fluids in nanochannels has significant instabilities, and numerical simulation failures are prone to occur at high Weissenberg numbers (Wi). In this study, the simplified Phan-Thien-Tanner viscoelastic fluid model is solved using the log-conformation tensor approach, and the effects of rheological parameters of the viscoelastic fluid, such as the Weissenberg number (Wi), extensibility parameter (ε), and viscosity ratio (β), on the flow characteristics and flow instability within the nanochannel are investigated. The results indicate that the variation of rheological parameters of viscoelastic fluids has a significant effect on the flow state and flow instability of fluids in nanochannels. When the rheological parameters are in a specific range, the flow velocity and outlet current in the nanochannel exhibit relatively regular periodic fluctuations. As the flow transitions from an up-and-down moving single-vortex state to a symmetric double-vortex state, the average velocity of the central axis in the nanochannel is increased by about 15%. Furthermore, when Wi increases from 150 to 400, the length and height of the vortex increase by 50% and 100%, respectively.

Marginal stability analyses for thermochemical convection and its implications for the dynamics of continental lithosphere and core–mantle boundary regions

SUMMARY Significant compositional differences may exist in the lithospheric mantle and above the core–mantle boundary (CMB) relative to the ambient mantle. The intrinsic density differences may affect the development of thermal boundary layer (TBL) instabilities associated with lithospheric delamination and formation of thermochemical plumes. In this study, we explored the instability of two-layer thermochemical fluid using two different techniques: marginal stability analysis with a propagator-matrix method and finite element modelling. We investigated both the instabilities in lithospheric mantle (i.e. lithospheric instability) and the mantle above the CMB (i.e. plume-forming instability) using a background temperature Tbg(z) with the TBL. For lithospheric instability, we found that two-layer fluid with free-slip boundary conditions mainly undergoes the same three different convective modes (i.e. two oscillatory convection modes and one layered convection regime) as that with no-slip boundary conditions reported in Jaupart et al. However, with free-slip boundary conditions, the transitions between these convection modes occur at larger values of buoyancy number B. Free-slip boundary conditions lead to smaller critical Rayleigh number Rac, but larger convective wavelength and oscillation frequency ωc, compared with those with no-slip boundary conditions. Our numerical modelling results demonstrate that Rac and ωc predicted from the classical marginal stability analyses using Tbg(z) with TBL temperature may have significant errors when the oscillatory period is comparable with or larger than the timescale of lithospheric thermal diffusion that causes Tbg(z) to vary with time significantly. In this case, using a more gently sloped background temperature profile ignoring the TBL temperature, the stability analysis predicts more accurate stability conditions, thus presenting an effective remedy to the stability analysis. For plume-forming instability, because of the reduced viscosity in the hot and compositionally dense bottom layer, the transition to the layered convection occurs at significantly smaller B values, and in the oscillatory convection regime, Rac is larger but ωc is smaller, compared with those for lithospheric instability. Finally, our study provides a successful benchmark of numerical models of thermochemical convection by comparing Rac and ωc from numerical models with those from the marginal stability analysis.

Kinetic and hydrodynamic post-Newtonian equations in the Brans–Dicke theory

Abstract A kinetic theory for the post-Newtonian Brans–Dicke theory is developed. The Boltzmann equation and the equilibrium Maxwell-Jüttner distribution function are determined from the knowledge of the components of the metric tensor and Christoffel symbols in the post-Newtonian Brans–Dicke theory. The hydrodynamic equations for the mass density, momentum density and mass-energy density are obtained from a transfer equation derived from the Boltzmann equation. The problem of self-gravitating fluid instabilities in the post-Newtonian Brans–Dicke theory is investigated.

Oscillation-free implicit pressure explicit concentration discontinuous Galerkin methods for compressible miscible displacements with applications in viscous fingering

The system of compressible miscible displacements is widely adopted to model surfactant flooding in enhanced oil recovery techniques, where a low-viscosity fluid is injected underground to replace the high-viscosity oil. When the mobility ratio of the injected fluid to oil is high, the waterflood front tends to be unstable and exhibits a finger-like growth pattern, known as viscous fingering. Due to its unstable nature, the viscous fingering phenomenon is sensitive to mesh orientation and numerical discretization. Therefore, high-order numerical methods are preferable to reduce numerical artifacts and mesh dependence. In this paper, we propose a high-order discontinuous Galerkin method for the coupled nonlinear system of compressible miscible displacements to simulate the viscous fingering fluid instability in porous media. We adopt the IMplicit Pressure Explicit Concentration time marching approach based on implicit-explicit Runge-Kutta methods to achieve high-order temporal accuracy. Additionally, we introduce an oscillation-free damping term to control the spurious oscillations encountered in the waterflood front due to the large gradient of saturation. We have conducted ample numerical tests in two space dimensions to demonstrate the effectiveness and robustness of the proposed schemes in recovering viscous fingering.

Nonlinear, single-mode, two-dimensional Rayleigh-Taylor instability in ideal media.

A model for the single mode, two-dimensional Rayleigh-Taylor instability in ideal, incompressible, immiscible, and inviscid fluids is developed as an extension of a previous linear model based on the Newton's second law [A. R. Piriz et al., Am. J. Phys. 74, 1095 (2006)0002-950510.1119/1.2358158]. It describes the transition from linear to nonlinear regimes and takes into account the mass of fluids participating in the motion during the instability evolution, including the laterally displaced mass. This latter feature naturally leads to the bubble and spike velocity saturation without requiring the usual drag term necessary in the well-known buoyancy-drag model (BDM). In addition, it also provides an explanation to the latter phase of bubble reacceleration without appealing to the vorticity generation due to the Kelvin-Helmholtz instability. The model is in perfect agreement with the BDM buoyancy-drag model, but, apart from extending its range of application, it solves many of its issues of concern and provides a more consistent physical picture.

Ambient printing of native oxides for ultrathin transparent flexible circuit boards.

Metal oxide films are essential in most electronic devices, yet they are typically deposited at elevated temperatures by using slow, vacuum-based processes. We printed native oxide films over large areas at ambient conditions by moving a molten metal meniscus across a target substrate. The oxide gently separates from the metal through fluid instabilities that occur in the meniscus, leading to uniform films free of liquid residue. The printed oxide has a metallic interlayer that renders the films highly conductive. The metallic character of the printed films promotes wetting of trace amounts of evaporated gold that would otherwise form disconnected islands on conventional oxide surfaces. The resulting ultrathin (<10 nanometers) conductors can be patterned into flexible circuits that are transparent, mechanically robust, and electrically stable, even at elevated temperatures.

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

Viscoelastic flow instabilities for enhanced heat transfer in battery pack cooling

Effective thermal management plays a critical role in ensuring peak performance for heavy-duty electric vehicle battery packs and high-performance electronics. To address the need for effective heat dissipation, liquid coolants are employed, which exhibit higher heat transfer performance compared to traditional air-cooled systems. However, the associated increased viscosity of such liquid coolants and constricted flow path lead to reduced mixing, thereby deteriorating the thermal performance of the system. Here, we explore the potential of utilizing viscoelastic liquid coolants for enhancing heat transfer from the heated walls of a constricted flow path with obstacles, modeling features of battery coolant channels. At low-Reynolds number regime, viscoelastic fluids exhibit instabilities which can disrupt the thermal boundary layer and enhance fluid mixing in the domain, thereby facilitating wall heat transfer. We employed computational fluid dynamics simulations, validated with experimental observations, to explore the effect of the flow geometry, Weissenberg number, solvent viscosity ratio, polymer extensibility, and Reynolds number on thermal performance. Our study indicates that the heat transfer is enhanced with flow instabilities in viscoelastic fluids, which can be generated exclusively in the flow channel with obstacles by increasing the Weissenberg number, reducing the solvent viscosity ratio, increasing the polymer extensibility, and increasing the Reynolds number. The enhancement in the heat transfer can exclusively be attributed to the elastic properties of the fluid in contrast to the inelastic shear-thinning fluids exhibiting no instability. The insights from the study can play an important role in the development of advanced, high-efficiency battery cooling systems, combining surface engineering and synthesis of optimized heat transfer liquids.

Gravitational instability in magnetized viscoelastic fluid with radiation pressure effect

This study examines the gravitational instability of a finitely conducting magnetized viscoelastic fluid, under the influence of radiation pressure within the framework of a generalized hydrodynamic fluid model. The general dispersion relation, valid for strongly and weakly coupled fluids under kinetic and hydrodynamic limits respectively, is derived using normal mode analysis for a transverse wave propagation mode. The obtained Jeans instability criteria, in terms of critical Jeans wavenumbers under both limits, are modified by the presence of radiative effects, viscoelastic compressional speed, and Alfvén wave velocity. The findings highlight that viscoelastic compressional speed and radiative pressure slow the Jeans instability's growth rate, exerting a stabilizing effect on the initiation of gravitational collapse. The present investigation provides valuable insights into the role of radiative mechanisms in the gravitational collapse within the strongly coupled region of molecular cloud clumps.

The Electro-Elastic Instability of Viscoelastic Fluid in a Microchannel with Obstacles Under Heterogeneous Surface Potential

The Electro-Elastic Instability of Viscoelastic Fluid in a Microchannel with Obstacles Under Heterogeneous Surface Potential

Water-Based Continuous Fabrication of Highly Elastic Electromagnetic Fibers.

Elastic electromagnetic fibers are promising building blocks for next-generation flexible electronics. However, fabrication of elastic fibers is still difficult and usually requires organic solvents or high temperatures, restricting their widespread applications. Furthermore, the continuous production of electromagnetic fibers has not been realized previously. In this study, we propose an ionic chelation strategy to continuously produce electromagnetic fibers with a magnetic liquid metal (MLM) as the core and elastic polyurethane as the sheath in water at room temperature. Sodium alginate (SA) has been introduced to rapidly chelate with calcium ions (Ca2+) in a coagulation bath to support the continuous spinning of waterborne polyurethane (WPU) as a sheath. Meanwhile, WPU-encapsulated MLM microparticles efficiently suppress the fluid instability of MLM for continuous extrusion as the core. The resultant fiber exhibits excellent mechanical performances (tensile strength and toughness up to 32 MPa and 124 MJ/m3, respectively), high conductive stability in large deformations (high conductivity of 7.6 × 104 S/m at 580% strain), and magnetoactive properties. The applications of this electromagnetic fiber have been demonstrated by conductance-stable wires, sensors, actuation, and electromagnetic interference shielding. This work offers a water-based molecular principle for efficient and green fabrication of multifunctional fibers and will inspire a series of applications.

Numerical investigation of mixing and heat transfer in a 7.9 m JP-8 pool fire

The response of objects engulfed in, and adjacent to, large-scale pool fires is of interest in accident and safety assessments. In this study, Fuego, a low-Mach turbulent reacting flow code was used to study conjugate heat transfer in a 7.9 m diameter JP-8 pool fire. Simulations were designed to replicate past experimental measurements (Blanchat et al., 2006) of incident heat flux to three cylindrical calorimeters in and around the pool fire. Two turbulent combustion models were compared directly - the eddy dissipation concept and a more recently developed unsteady flamelet model. First and second order spatial and temporal discretization schemes were also compared to assess the performance of low-dissipation numerical operators. Heat flux predictions to the transportation size calorimeter outside the fire were within experimental uncertainties. Inside the fire, experimental measurements were higher than predicted values and may have been a consequence of soot deposition and augmented participating media radiation from soot and fuel vapor. Simulation predictions improved in cases where turbulent kinetic energy and mixing were more resolved. This work, and others referenced herein, suggest that spatial resolution on the order of 0.5–1.0 cm may be required to fully resolve fluid instabilities, vortex production between fire plumes and crosswind, soot production, and fuel-air mixing. This presents a substantial computational challenge for safety assessments of engulfed objects in fully turbulent pool fires.

Effects of gravity and anisotropy on the convective instability in a nanofluid porous medium

Controlling thermal loads in heat transfer systems that produce high heat flux is a real mechanical challenge. Various cooling/heating techniques are available to achieve the reliability requirements of electronic devices. Among them, the use of nanofluids is one of the important techniques that provide long term stability, higher cooling rates, decreased pumping power needs, and cost-effectiveness. Instability in fluids typically occurs when a disturbance affects the equilibrium of forces of external, inertial, or viscous origin acting on a system. The disturbances may arise from the surroundings due to small changes in the boundary conditions, surface roughness, free stream turbulence, noise, etc. In several circumstances, these disturbances are very small and hence cannot be measured with the available instruments but can be observed after the instability sets in. In this study, we focus on examining the effect of time-periodic vertical gravity modulation on the initiation of natural convection in a porous medium saturated with a nanofluid. The porous medium is characterized by being heated at the bottom and possessing anisotropic properties in terms of permeability and thermal diffusivity. Specifically, we investigate the onset of convection in water-based nanofluids that contain metallic (Cu) and metal oxide ( A l 2 O 3 ) particles. The Khanafer-Vafai-Lightstone (KVL) nanofluid model is employed to describe the behavior of the nanofluids, and the thresholds for convective instability are determined using the Mathieu functions. In addition, we study the instabilities of synchronous and subharmonic modes and explore their transitions along with practical applications.

Generation of monoclonal antibody against 6-Keto PGF1α and development of ELISA for its quantification in culture medium

Prostacyclin or prostaglandin I2 (PGI2), a metabolite of arachidonic cyclooxygenase pathway, has been demonstrated as an effector of adipocyte differentiation. However, due to its instability in biological fluid, it is difficult to evaluate the role of PGI2 in regulating adipocyte differentiation in different stages in culture. Therefore, this study aimed to establish a simple and rapid method for the production of monoclonal antibody against 6-Keto PGF1α, a stable PGI2 metabolite, and its quantification to determine the role of PGI2 in culture medium. Eight-week-old female BALB/c mice were immunized with the hapten of 6-Keto PGF1α and BSA for several weeks until a higher antibody titer (absorbance value > 0.9 at 1000-times dilution) against 6-Keto PGF1α was found. Then, fusion of antibody-producing spleen lymphocytes with SP-2 myeloma cells and thymocytes was performed and cultured in HAT-medium supplemented with hypoxanthine, aminopterin, and thymine. Specific antibody-producing cells (M2-A4-B8-D10) against 6-Keto PGF1α were identified and separated. A standard ELISA calibration curve was developed with 100% reactivity for 6-Keto-PGF 1 α ranging from 0.26 pg to 6.44 ng corresponding to 90% and 10% of the maximum binding capacity for the immobilized antigen respectively. This method can easily be applied to monitor PGI2 regulation in different stages of cultured adipocytes to reveal the regulatory roles of PGI2 in maintaining homeostasis and adipocyte differentiation.