Search

It is essential that mesoscopic simulations of reactive systems reproduce the correct statistical distributions at thermodynamic equilibrium. By considering a compressible fluctuating hydrodynamics (FHD) simulation method of ideal gas mixtures undergoing reversible reactions described by the chemical Langevin equations, we show that thermodynamic consistency in reaction rates and the use of instantaneous temperatures for the evaluation of reaction rates is required for fluctuations for the overall system to be correct. We then formulate the required properties of a thermodynamically-consistent reaction (TCR) model. As noted in the literature, while reactions are often discussed in terms of forward and reverse rates, these rates should not be modeled independently because they must be compatible with thermodynamic equilibrium for the system. Using a simple TCR model where each chemical species has constant heat capacity, we derive the explicit condition that the forward and reverse reaction rate constants must satisfy in order for the system to be thermodynamically consistent. We perform equilibrium and non-equilibrium simulations of ideal gas mixtures undergoing a reversible dimerization reaction to measure the fluctuational behavior of the system numerically. We confirm that FHD simulations with the TCR model give the correct static structure factor of equilibrium fluctuations. For the statistically steady simulation of a gas mixture between two isothermal walls with different temperatures, we show using the TCR model that the temperature variance agrees with the corresponding thermodynamic-equilibrium temperature variance in the interior of the system, whereas noticeable deviations are present in regions near walls, where chemistry is far from equilibrium.

We simulate the mesoscopic dynamics of droplets formed by phase separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn-Hillard formulation with incompressible fluctuating hydrodynamics; for sessile droplets the fluid-solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (greater than 90 degrees) on both slip and no-slip surfaces.

These notes are an introduction to fluctuating hydrodynamics (FHD) and the formulation of numerical schemes for the resulting stochastic partial differential equations (PDEs). Fluctuating hydrodynamics was originally introduced by Landau and Lifshitz as a way to put thermal fluctuations into a continuum framework by including a stochastic forcing to each dissipative transport process (e.g., heat flux). While FHD has been useful in modeling transport and fluid dynamics at the mesoscopic scale, theoretical calculations have been feasible only with simplifying assumptions. As such there is great interest in numerical schemes for Computational Fluctuating Hydrodynamics (CFHD). There are a variety of algorithms (e.g., spectral, finite element, lattice Boltzmann) but in this introduction we focus on finite volume schemes. Accompanying these notes is a demonstration program in Python available on GitHub.

State redistribution (SRD) is a recently developed technique for stabilizing cut cells that result from finite-volume embedded boundary methods. SRD has been successfully applied to a variety of compressible and incompressible flow problems. When used in conjunction with adaptive mesh refinement (AMR), additional steps are needed to preserve the accuracy and conservation properties of the solution if the embedded boundary is not restricted to a single level of the mesh hierarchy. In this work, we extend the weighted state redistribution algorithm to cases where cut cells live at or near a coarse-fine interface within the domain. The resulting algorithm maintains conservation and is demonstrated on several two- and three-dimensional example problems., Comment: 29 pages, 10 figures, journal article

We study the role of steric effects on the induced-charge electro-osmosis (ICEO) phenomenon using a recently developed mesoscale fluid model. A hybrid Eulerian-Lagrangian method is used to simulate the dynamics of discrete immersed ions in a thermally fluctuating solvent near a metallic plate embedded in the dielectric interface. We observe that the characteristic velocity scales almost linearly with electric field when the generated $\zeta$-potentials exceed the order of the thermal voltage, as opposed to a quadratic scaling predicted by Helmholtz-Smoluchowski equation, although qualitative agreement with experiments and theories is obtained at low electric fields. Our simulations reveal that the steric effects play a crucial role at strong electric fields, which is observed from the aggregation of ions towards the center of the metal plate instead of at the edges, and the overcharging of co-ions to the surface charge near the electric double layer. A comparison to a continuum electrolyte model also highlights significant differences in charge distribution and flow field that are attributed to the steric repulsion between ions., Comment: 34 pages, 11 figures

The Rayleigh-Plateau instability occurs when surface tension makes a fluid column become unstable to small perturbations. At nanometer scales, thermal fluctuations are comparable to surface energy densities. Consequently, at these scales, thermal fluctuations play a significant role in the dynamics of the instability. These microscopic effects have previously been investigated numerically using particle-based simulations, such as molecular dynamics, and stochastic partial differential equation based hydrodynamic models, such as stochastic lubrication theory. In this paper we present an incompressible fluctuating hydrodynamics model with a diffuse-interface formulation for binary fluid mixtures designed for the study of stochastic interfacial phenomena. An efficient numerical algorithm is outlined and validated in numerical simulations of stable equilibrium interfaces. We present results from simulations of the Rayleigh-Plateau instability for long cylinders pinching into droplets for Ohnesorge numbers of Oh = 0.5 and 5.0. Both stochastic and perturbed deterministic simulations are analyzed and ensemble results show significant differences in the temporal evolution of the minimum radius near pinching. Short cylinders, with lengths less than their circumference, were also investigated. As previously observed in molecular dynamics simulations, we find that thermal fluctuations cause these to pinch in cases where a perturbed cylinder would be stable deterministically. Finally we show that the fluctuating hydrodynamics model can be applied to study a broader range of surface-tension driven phenomena.

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, $N$, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable $N$'s, the finite size effect is significantly larger when $N$ is even than when $N$ is odd. Moreover, as $N$ increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

We present a numerical formulation for the solution of non-isothermal, compressible, Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are significantly simplified and compact discretization of the diffusive and stochastic momentum fluxes, and an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously in Balakrishnan \emph{et al.} [Phys. Rev. E 89, 013017 (2014)]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients, and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo (DSMC) method., Comment: accepted version; published in Phys. Rev. E

In this article we develop an algorithm for the efficient simulation of electrolytes in the presence of physical boundaries. In previous work the Discrete Ion Stochastic Continuum Overdamped Solvent (DISCOS) algorithm was derived for triply periodic domains, and was validated through ion-ion pair correlation functions and Debye-H{\"u}ckel-Onsager theory for conductivity, including the Wien effect for strong electric fields. In extending this approach to include an accurate treatment of physical boundaries we must address several important issues. First, the modifications to the spreading and interpolation operators necessary to incorporate interactions of the ions with the boundary are described. Next we discuss the modifications to the electrostatic solver to handle the influence of charges near either a fixed potential or dielectric boundary. An additional short-ranged potential is also introduced to represent interaction of the ions with a solid wall. Finally, the dry diffusion term is modified to account for the reduced mobility of ions near a boundary, which introduces an additional stochastic drift correction. Several validation tests are presented confirming the correct equilibrium distribution of ions in a channel. Additionally, the methodology is demonstrated using electro-osmosis and induced charge electro-osmosis, with comparison made to theory and other numerical methods. Notably, the DISCOS approach achieves greater accuracy than a continuum electrostatic simulation method. We also examine the effect of under-resolving hydrodynamic effects using a `dry diffusion' approach, and find that considerable computational speedup can be achieved with a negligible impact on accuracy., Comment: 27 pages, 15 figures

State redistribution is an algorithm that stabilizes cut cells for embedded boundary grid methods. This work extends the earlier algorithm in several important ways. First, state redistribution is extended to three spatial dimensions. Second, we discuss several algorithmic changes and improvements motivated by the more complicated cut cell geometries that can occur in higher dimensions. In particular, we introduce a weighted version with less dissipation in an easily generalizable framework. Third, we demonstrate that state redistribution can also stabilize a solution update that includes both advective and diffusive contributions. The stabilization algorithm is shown to be effective for incompressible as well as compressible reacting flows. Finally, we discuss the implementation of the algorithm for several exascale-ready simulation codes based on AMReX, demonstrating ease of use in combination with domain decomposition, hybrid parallelism and complex physics.

Using fluctuating hydrodynamics we investigate the effect of thermal fluctuations in the dissipation range of homogeneous, isotropic turbulence. Simulations confirm theoretical predictions that the energy spectrum is dominated by these fluctuations at length scales comparable to the Kolmogorov length. We also find that the extreme intermittency in the far-dissipation range predicted by Kraichnan is replaced by Gaussian thermal equipartition.

We present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an $hp$ adaptive mesh refinement ($hp$-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The $hp$-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The $hp$-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

In this paper we develop a methodology for the mesoscale simulation of strong electrolytes. The methodology is an extension of the Fluctuating Immersed Boundary (FIB) approach that treats a solute as discrete Lagrangian particles that interact with Eulerian hydrodynamic and electrostatic fields. In both cases the Immersed Boundary (IB) method of Peskin is used for particle-field coupling. Hydrodynamic interactions are taken to be overdamped, with thermal noise incorporated using the fluctuating Stokes equation, including a "dry diffusion" Brownian motion to account for scales not resolved by the coarse-grained model of the solvent. Long range electrostatic interactions are computed by solving the Poisson equation, with short range corrections included using a novel immersed-boundary variant of the classical Particle-Particle Particle-Mesh (P3M) technique. Also included is a short range repulsive force based on the Weeks-Chandler-Andersen (WCA) potential. The new methodology is validated by comparison to Debye-H{\"u}ckel theory for ion-ion pair correlation functions, and Debye-H{\"u}ckel-Onsager theory for conductivity, including the Wein effect for strong electric fields. In each case good agreement is observed, provided that hydrodynamic interactions at the typical ion-ion separation are resolved by the fluid grid., Comment: 30 pages, 12 figures, 2 tables

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.

A new algorithm is developed to tackle the issue of sampling non-Gaussian model parameter posterior probability distributions that arise from solutions to Bayesian inverse problems. The algorithm aims to mitigate some of the hurdles faced by traditional Markov Chain Monte Carlo (MCMC) samplers, through constructing proposal probability densities that are both, easy to sample and that provide a better approximation to the target density than a simple Gaussian proposal distribution would. To achieve that, a Gaussian proposal distribution is augmented with a Gaussian Process (GP) surface that helps capture non-linearities in the log-likelihood function. In order to train the GP surface, an iterative approach is adopted for the optimal selection of points in parameter space. Optimality is sought by maximizing the information gain of the GP surface using a minimum number of forward model simulation runs. The accuracy of the GP-augmented surface approximation is assessed in two ways. The first consists of comparing predictions obtained from the approximate surface with those obtained through running the actual simulation model at hold-out points in parameter space. The second consists of a measure based on the relative variance of sample weights obtained from sampling the approximate posterior probability distribution of the model parameters. The efficacy of this new algorithm is tested on inferring reaction rate parameters in a 3-node and 6-node network toy problems, which imitate idealized reaction networks in combustion applications.

We propose an algorithm for low Mach number reacting flows subjected to electric field that includes the chemical production and transport of charged species. This work is an extension of a multi-implicit spectral deferred correction (MISDC) algorithm designed to advance the conservation equations in time at scales associated with advective transport. The fast and nontrivial interactions of electrons with the electric field are treated implicitly using a Jacobian-Free Newton Krylov approach for which a preconditioning strategy is developed. Within the MISDC framework, this enables a close and stable coupling of diffusion, reactions and dielectric relaxation terms with advective transport and is shown to exhibit second-order convergence in space and time. The algorithm is then applied to a series of steady and unsteady problems to demonstrate its capability and stability. Although developed in a one-dimensional case, the algorithmic ingredients are carefully designed to be amenable to multidimensional applications.

A fluid in a non-equilibrium state exhibits long-ranged correlations of its hydrodynamic fluctuations. In this article, we examine the effect of a transpiration interface on these correlations -- specifically, we consider a dilute gas in a domain bisected by the interface. The system is held in a non-equilibrium steady state by using isothermal walls to impose a temperature gradient. The gas is simulated using both direct simulation Monte Carlo (DSMC) and fluctuating hydrodynamics (FHD). For the FHD simulations two models are developed for the interface based on master equation and Langevin approaches. For appropriate simulation parameters, good agreement is observed between DSMC and FHD results with the latter showing a significant advantage in computational speed. For each approach we quantify the effects of transpiration on long-ranged correlations in the hydrodynamic variables.

At mesoscopic scales electrolyte solutions are modeled by the fluctuating generalized Poisson-Nernst-Planck (PNP) equations [J.-P. P\'eraud et al., Phys. Rev. F, 1(7):074103, 2016]. However, at length and time scales larger than the Debye scales, electrolytes are effectively electroneutral, and the charged-fluid PNP equations become too stiff to solve numerically. Here we formulate the isothermal incompressible equations of fluctuating hydrodynamics for reactive multispecies mixtures involving charged species in the electroneutral limit, and design a numerical algorithm to solve these equations. Our model does not assume a dilute electrolyte solution but rather treats all species on an equal footing, accounting for cross-diffusion and non-ideality using Maxwell-Stefan theory. By enforcing local electroneutrality as a constraint, we obtain an elliptic equation for the electric potential that replaces the Poisson equation in the fluctuating PNP equations. We develop a second-order midpoint predictor-corrector algorithm to solve either the charged-fluid or electroneutral equations with only a change of the elliptic solver. We use the electroneutral algorithm to study a gravitational fingering instability, triggered by thermal fluctuations, at an interface where an acid and base react to neutralize each other. Our results demonstrate that, because the four ions diffuse with very different coefficients, one must treat each ion as an individual species, and cannot treat the acid, base, and salt as neutral species. This emphasizes the differences between electrodiffusion and classical Fickian diffusion, even at electroneutral scales., Comment: To appear in Phys. Rev. Fluids

We apply fluctuating hydrodynamics to strong electrolyte mixtures to compute the concentration corrections for chemical potential, diffusivity, and conductivity. We show these corrections to be in agreement with the limiting laws of Debye, H\"uckel, and Onsager. We compute explicit corrections for a symmetric ternary mixture and find that the co-ion Maxwell-Stefan diffusion coefficients can be negative, in agreement with experimental findings., Comment: Submitted to Current Opinion in Electrochemistry

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.

In this paper we present a fourth-order in space and time block-structured adaptive mesh refinement algorithm for the compressible multicomponent reacting Navier-Stokes equations. The algorithm uses a finite volume approach that incorporates a fourth-order discretization of the convective terms. The time stepping algorithm is based on a multi-level spectral deferred corrections method that enables explicit treatment of advection and diffusion coupled with an implicit treatment of reactions. The temporal scheme is embedded in a block-structured adaptive mesh refinement algorithm that includes subcycling in time with spectral deferred correction sweeps applied on levels. Here we present the details of the multi-level scheme paying particular attention to the treatment of coarse-fine boundaries required to maintain fourth-order accuracy in time. We then demonstrate the convergence properties of the algorithm on several test cases including both nonreacting and reacting flows. Finally we present simulations of a vitiated dimethyl ether jet in 2D and a turbulent hydrogen jet in 3D, both with detailed kinetics and transport.

Fluctuating hydrodynamics (FHD) provides a framework for modeling microscopic fluctuations in a manner consistent with statistical mechanics and nonequilibrium thermodynamics. This paper presents an FHD formulation for isothermal reactive incompressible liquid mixtures with stochastic chemistry. Fluctuating multispecies mass diffusion is formulated using a Maxwell-Stefan description without assuming a dilute solution, and momentum dynamics is described by a stochastic Navier-Stokes equation for the fluid velocity. We consider a thermodynamically consistent generalization for the law of mass action for non-dilute mixtures and use it in the chemical master equation (CME) to model reactions as a Poisson process. The FHD approach provides remarkable computational efficiency over traditional reaction-diffusion master equation methods when the number of reactive molecules is large, while also retaining accuracy even when there are as few as ten reactive molecules per hydrodynamic cell. We present a numerical algorithm to solve the coupled FHD and CME equations and validate it on both equilibrium and nonequilibrium problems. We simulate a diffusively-driven gravitational instability in the presence of an acid-base neutralization reaction, starting from a perfectly flat interface. We demonstrate that the coupling between velocity and concentration fluctuations dominate the initial growth of the instability.

In this letter we analyze the effects of an externally applied electric field on thermal fluctuations for a fluid containing charged species. We show in particular that the fluctuating Poisson-Nernst-Planck equations for charged multispecies diffusion coupled with the fluctuating fluid momentum equation, result in enhanced charge transport. Although this transport is advective in nature, it can macroscopically be represented as electrodiffusion with renormalized electric conductivity. We calculate the renormalized electric conductivity by deriving and integrating the structure factor coefficients of the fluctuating quantities and show that the renormalized electric conductivity and diffusion coefficients are consistent although they originate from different noise terms. In addition, the fluctuating hydrodynamics approach recovers the electrophoretic and relaxation corrections obtained by Debye-Huckel-Onsager theory, and provides a quantitative theory that predicts a non-zero cross-diffusion Maxwell-Stefan coefficient that agrees well with experimental measurements. Finally, we show that strong applied electric fields result in anisotropically enhanced velocity fluctuations and reduced fluctuations of salt concentrations., Comment: 12 pages, 1 figure

We develop numerical methods for reaction-diffusion systems based on the equations of fluctuating hydrodynamics (FHD). While the FHD formulation is formally described by stochastic partial differential equations (SPDEs), it becomes similar to the reaction-diffusion master equation (RDME) description when those SPDEs are spatially discretized and reactions are modeled as a source term having Poisson fluctuations. However, unlike the RDME, the FHD description naturally extends from the regime where fluctuations are strong, i.e., each hydrodynamic cell has few (reactive) molecules, to regimes with moderate or weak fluctuations, and ultimately to the deterministic limit. By treating diffusion implicitly, we avoid the severe restriction on time step size that limits all methods based on explicit treatments of diffusion, and construct numerical methods that are more efficient than RDME methods, without compromising accuracy. We find that an implicit midpoint tau leaping scheme attains second-order weak accuracy for the linearized FHD, and gives an accurate and stable structure factor for a time step size an order of magnitude larger than the hopping time scale of diffusing molecules. We study the numerical accuracy of our methods for the Schl\"ogl reaction-diffusion model both in and out of thermodynamic equilibrium. We demonstrate and quantify the importance of thermodynamic fluctuations to the formation of a two-dimensional Turing-like pattern, and examine the effect of fluctuations on three-dimensional chemical front propagation. By comparing FHD simulations to deterministic reaction-diffusion simulations, we show that fluctuations accelerate pattern formation in spatially homogeneous systems, and lead to a qualitatively-different disordered pattern behind a traveling wave., Comment: To appear in J. Chem. Phys., 2017

In parameter estimation problems one computes a posterior distribution over uncertain parameters defined jointly by a prior distribution, a model, and noisy data. Markov Chain Monte Carlo (MCMC) is often used for the numerical solution of such problems. An alternative to MCMC is importance sampling, which can exhibit near perfect scaling with the number of cores on high performance computing systems because samples are drawn independently. However, finding a suitable proposal distribution is a challenging task. Several sampling algorithms have been proposed over the past years that take an iterative approach to constructing a proposal distribution. We investigate the applicability of such algorithms by applying them to two realistic and challenging test problems, one in subsurface flow, and one in combustion modeling. More specifically, we implement importance sampling algorithms that iterate over the mean and covariance matrix of Gaussian or multivariate t-proposal distributions. Our implementation leverages massively parallel computers, and we present strategies to initialize the iterations using "coarse" MCMC runs or Gaussian mixture models.

We formulate and study computationally the low Mach number fluctuating hydrodynamic equations for electrolyte solutions. We are interested in studying transport in mixtures of charged species at the mesoscale, down to scales below the Debye length, where thermal fluctuations have a significant impact on the dynamics. Continuing our previous work on fluctuating hydrodynamics of multicomponent mixtures of incompressible isothermal miscible liquids (A. Donev, et al., Physics of Fluids, 27, 3, 2015), we now include the effect of charged species using a quasielectrostatic approximation. Localized charges create an electric field, which in turn provides additional forcing in the mass and momentum equations. Our low Mach number formulation eliminates sound waves from the fully compressible formulation and leads to a more computationally efficient quasi-incompressible formulation. We demonstrate our ability to model saltwater (NaCl) solutions in both equilibrium and nonequilibrium settings. We show that our algorithm is second-order in the deterministic setting, and for length scales much greater than the Debye length gives results consistent with an electroneutral/ambipolar approximation. In the stochastic setting, our model captures the predicted dynamics of equilibrium and nonequilibrium fluctuations. We also identify and model an instability that appears when diffusive mixing occurs in the presence of an applied electric field., Comment: 37 pages, 5 figures

The dynamics of helium shell convection driven by nuclear burning establish the conditions for runaway in the sub-Chandrasekhar mass, double detonation model for Type Ia supernovae, as well as for a variety of other explosive phenomena. We explore these convection dynamics for a range of white dwarf core and helium shell masses in three dimensions using the low Mach number hydrodynamics code Maestro. We present calculations of the bulk properties of this evolution, including time-series evolution of global diagnostics, lateral averages of the 3D state, and the global 3D state. We find a variety of outcomes including quasi-equilibrium, localized runaway, and nova-like runaway. Our results suggest the double detonation progenitor model is promising, that 3D, dynamic convection plays a key role, and that these systems warrant further study., Comment: 23 pages, 18 figures, 2 tables, accepted in ApJ

Fluctuating hydrodynamics provides a model for fluids at mesoscopic scales where thermal fluctuations can have a significant impact on the behavior of the system. Here we investigate a model for fluctuating hydrodynamics of a single component, multiphase flow in the neighborhood of the critical point. The system is modeled using a compressible flow formulation with a van der Waals equation of state, incorporating a Korteweg stress term to treat interfacial tension. We present a numerical algorithm for modeling this system based on an extension of algorithms developed for fluctuating hydrodynamics for ideal fluids. The scheme is validated by comparison of measured structure factors and capillary wave spectra with equilibrium theory. We also present several non-equilibrium examples to illustrate the capability of the algorithm to model multi-phase fluid phenomena in a neighborhood of the critical point. These examples include a study of the impact of fluctuations on the spinodal decomposition following a rapid quench, as well as the piston effect in a cavity with supercooled walls. The conclusion in both cases is that thermal fluctuations affect the size and growth of the domains in off-critical quenches.