Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA Low Mach Number Fluctuating Hydrodynamics Model For Ionic Liquids
89
0
0.0
(
0
)
Authors
Katherine Klymko
Ask ChatGPT about the research
No Arabic abstract
We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.
rate research
Read More
We develop numerical schemes for solving the isothermal compressible and incompressible equations of fluctuating hydrodynamics on a grid with staggered momenta. We develop a second-order accurate spatial discretization of the diffusive, advective and stochastic fluxes that satisfies a discrete fluctuation-dissipation balance, and construct temporal discretizations that are at least second-order accurate in time deterministically and in a weak sense. Specifically, the methods reproduce the correct equilibrium covariances of the fluctuating fields to third (compressible) and second (incompressible) order in the time step, as we verify numerically. We apply our techniques to model recent experimental measurements of giant fluctuations in diffusively mixing fluids in a micro-gravity environment [A. Vailati et. al., Nature Communications 2:290, 2011]. Numerical results for the static spectrum of non-equilibrium concentration fluctuations are in excellent agreement between the compressible and incompressible simulations, and in good agreement with experimental results for all measured wavenumbers.
We develop a one-dimensional model for the unsteady fluid--structure interaction (FSI) between a soft-walled microchannel and viscous fluid flow within it. A beam equation, which accounts for both transverse bending rigidity and nonlinear axial tension, is coupled to a one-dimensional fluid model obtained from depth-averaging the two-dimensional incompressible Navier--Stokes equations across the channel height. Specifically, the Navier--Stokes equations are scaled in the viscous lubrication limit relevant to microfluidics. The resulting set of coupled nonlinear partial differential equations is solved numerically through a segregated approach employing fully-implicit time stepping. We explore both the static and dynamic FSI behavior of this example microchannel system by varying a reduced Reynolds number $Re$, which necessarily changes the Strouhal number $St$, while we keep the geometry and a modified dimensionless Youngs modulus $Sigma$ fixed. At steady state, an order-of-magnitude analysis (balancing argument) shows that the axially-averaged pressure in the flow, $langle Prangle$, exhibits two different scaling regimes, while the maximum deformation of the top wall of the channel, $H_{mathrm{max}}$, can fall into four different regimes, depending on the magnitudes of $Re$ and $Sigma$. These regimes are physically explained as resulting from the competition between the inertial and viscous forces in the fluid flow as well as the bending resistance and tension in the elastic wall. Finally, the linear stability of the steady inflated microchannel shape is assessed via a modal analysis, showing the existence of many highly oscillatory but stable modes, which further highlights the computational challenge of simulating unsteady FSIs.
Flows in which the primary features of interest do not rely on high-frequency acoustic effects, but in which long-wavelength acoustics play a nontrivial role, present a computational challenge. Integrating the entire domain with low-Mach-number methods would remove all acoustic wave propagation, while integrating the entire domain with the fully compressible equations can in some cases be prohibitively expensive due to the CFL time step constraint. For example, simulation of thermoacoustic instabilities might require fine resolution of the fluid/chemistry interaction but not require fine resolution of acoustic effects, yet one does not want to neglect the long-wavelength wave propagation and its interaction with the larger domain. The present paper introduces a new multi-level hybrid algorithm to address these types of phenomena. In this new approach, the fully compressible Euler equations are solved on the entire domain, potentially with local refinement, while their low-Mach-number counterparts are solved on subregions of the domain with higher spatial resolution. The finest of the compressible levels communicates inhomogeneous divergence constraints to the coarsest of the low-Mach-number levels, allowing the low-Mach-number levels to retain the long-wavelength acoustics. The performance of the hybrid method is shown for a series of test cases, including results from a simulation of the aeroacoustic propagation generated from a Kelvin-Helmholtz instability in low-Mach-number mixing layers. It is demonstrated that compared to a purely compressible approach, the hybrid method allows time-steps two orders of magnitude larger at the finest level, leading to an overall reduction of the computational time by a factor of 8.
We report the statistical properties of temperature and thermal energy dissipation rate in low-Prandtl number turbulent Rayleigh-Benard convection. High resolution two-dimensional direct numerical simulations were carried out for the Rayleigh number ($Ra$) of $10^{6} le Ra le 10^{7}$ and the Prandtl number ($Pr$) of 0.025. Our results show that the global heat transport and momentum scaling in terms of Nusselt number ($Nu$) and Reynolds number ($Re$) are $Nu=0.21Ra^{0.25}$ and $Re=6.11Ra^{0.50}$, respectively, indicating that the scaling exponents are smaller than those for moderate-Prandtl number fluids (such as water or air) in the same convection cell. In the central region of the cell, probability density functions (PDFs) of temperature profiles show stretched exponential peak and the Gaussian tail; in the sidewall region, PDFs of temperature profiles show a multimodal distribution at relative lower $Ra$, while they approach the Gaussian profile at relative higher $Ra$. We split the energy dissipation rate into contributions from bulk and boundary layers and found the locally averaged thermal energy dissipation rate from the boundary layer region is an order of magnitude larger than that from the bulk region. Even if the much smaller volume occupied by the boundary layer region is considered, the globally averaged thermal energy dissipation rate from the boundary layer region is still larger than that from the bulk region. We further numerically determined the scaling exponents of globally averaged thermal energy dissipation rates as functions of $Ra$ and $Re$.
We propose a new model of turbulence for use in large-eddy simulations (LES). The turbulent force, represented here by the turbulent Lamb vector, is divided in two contributions. The contribution including only subfilter fields is deterministically modeled through a classical eddy-viscosity. The other contribution including both filtered and subfilter scales is dynamically computed as solution of a generalized (stochastic) Langevin equation. This equation is derived using Rapid Distortion Theory (RDT) applied to the subfilter scales. The general friction operator therefore includes both advection and stretching by the resolved scale. The stochastic noise is derived as the sum of a contribution from the energy cascade and a contribution from the pressure. The LES model is thus made of an equation for the resolved scale, including the turbulent force, and a generalized Langevin equation integrated on a twice-finer grid. The model is validated by comparison to DNS and is tested against classical LES models for isotropic homogeneous turbulence, based on eddy viscosity. We show that even in this situation, where no walls are present, our inclusion of backscatter through the Langevin equation results in a better description of the flow.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
