Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStaggered Schemes for Fluctuating Hydrodynamics
167
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
We develop a highly efficient numerical method to simulate small-amplitude flapping propulsion by a flexible wing in a nearly inviscid fluid. We allow the wings elastic modulus and mass density to vary arbitrarily, with an eye towards optimizing these distributions for propulsive performance. The method to determine the wing kinematics is based on Chebyshev collocation of the 1D beam equation as coupled to the surrounding 2D fluid flow. Through small-amplitude analysis of the Euler equations (with trailing-edge vortex shedding), the complete hydrodynamics can be represented by a nonlocal operator that acts on the 1D wing kinematics. A class of semi-analytical solutions permits fast evaluation of this operator with $O(N log N)$ operations, where $N$ is the number of collocation points on the wing. This is in contrast to the minimum $O(N^2)$ operations required by a direct 2D fluid solver. The coupled wing-fluid problem is thus recast as a PDE with nonlocal operator, which we solve using a preconditioned iterative method. These techniques yield a solver of near-optimal complexity, $O(N log N)$, allowing one to rapidly search the infinite-dimensional parameter space of all possible material distributions and even perform optimization over this space.
We present a hybrid spectral element-Fourier spectral method for solving the coupled system of Navier-Stokes and Cahn-Hilliard equations to simulate wall-bounded two-phase flows in a three-dimensional domain which is homogeneous in at least one direction. Fourier spectral expansions are employed along the homogeneous direction and $C^0$ high-order spectral element expansions are employed in the other directions. A critical component of the method is a strategy we developed in a previous work for dealing with the variable density/viscosity of the two-phase mixture, which makes the efficient use of Fourier expansions in the current work possible for two-phase flows with different densities and viscosities for the two fluids. The attractive feature of the presented method lies in that the two-phase computations in the three-dimensional space are transformed into a set of de-coupled two-dimensional computations in the planes of the non-homogeneous directions. The overall scheme consists of solving a set of de-coupled two-dimensional equations for the flow and phase-field variables in these planes. The linear algebraic systems for these two-dimensional equations have constant coefficient matrices that need to be computed only once and can be pre-computed. We present ample numerical simulations for different cases to demonstrate the accuracy and capability of the presented method in simulating the class of two-phase problems involving solid walls and moving contact lines.
We present a family of p-enrichment schemes. These schemes may be separated into two basic classes: the first, called emph{fixed tolerance schemes}, rely on setting global scalar tolerances on the local regularity of the solution, and the second, called emph{dioristic schemes}, rely on time-evolving bounds on the local variation in the solution. Each class of $p$-enrichment scheme is further divided into two basic types. The first type (the Type I schemes) enrich along lines of maximal variation, striving to enhance stable solutions in areas of highest interest. The second type (the Type II schemes) enrich along lines of maximal regularity in order to maximize the stability of the enrichment process. Each of these schemes are tested over a pair of model problems arising in coastal hydrology. The first is a contaminant transport model, which addresses a declinature problem for a contaminant plume with respect to a bay inlet setting. The second is a multicomponent chemically reactive flow model of estuary eutrophication arising in the Gulf of Mexico.
The uncertainty relations in hydrodynamics are numerically studied. We first give a review for the formulation of the generalized uncertainty relations in the stochastic variational method (SVM), following the paper by two of the present authors [Phys. Lett. A382, 1472 (2018)]. In this approach, the origin of the finite minimum value of uncertainty is attributed to the non-differentiable (virtual) trajectory of a quantum particle and then both of the Kennard and Robertson-Schr{o}dinger inequalities in quantum mechanics are reproduced. The same non-differentiable trajectory is applied to the motion of fluid elements in hydrodynamics. By introducing the standard deviations of position and momentum for fluid elements, the uncertainty relations in hydrodynamics are derived. These are applicable even to the Gross-Pitaevskii equation and then the field-theoretical uncertainty relation is reproduced. We further investigate numerically the derived relations and find that the behaviors of the uncertainty relations for liquid and gas are qualitatively different. This suggests that the uncertainty relations in hydrodynamics are used as a criterion to classify liquid and gas in fluid.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
