Researcher profile

Fathollah Varnik

Fathollah Varnik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

Multi-phase-field modeling of microstructure evolution in metallic foams

This paper represents a model for microstructure formation in metallic foams based on the multi-phase-field (MPF) approach. By the use of a no-coalescence boundary condition within this MPF-framework, it is possible to completely prevent coalescence of bubbles and thus focus on the formation of a closed porous microstructure. A modification of this non-wetting criterion allows for the controlled initiation of coalescence and the evolution of open structures. The method is validated and used to simulate foam structure formation both in two and three dimensions.

preprint2020arXiv

Self similarity of the expanding universe as understood by quantum-phase-fields

An idea of the universe as a self-contained system of interacting fields as a closed doublon network is developed. The characteristic scale of this system is considered emergent from general principles. Self similarity of patterns in the universe can be understood by introducing a characteristic scale with respect towhich patterns are evaluated. The developed physical-mathematical model is in good agreement with cosmological length scales and gives a rationale of the empirical fact for an expanding universe. Within the framework of the de~Broglie-Bohm double solution program, the model can be applied to rationalize the existence of a ground potential and a cosmological constant.

preprint2020arXiv

Simulation of Capillary-Driven Kinetics with Multi-Phase-Field and Lattice-Boltzmann Method

We propose a combined computational approach based on the multi-phase-field and the lattice Boltzmann method for the motion of solid particles under the action of capillary forces. The accuracy of the method is analyzed by comparison with the analytic solutions for the motion of two parallel plates of finite extension connected by a capillary bridge. The method is then used to investigate the dynamics of multiple spherical solid bodies connected via capillary bridges. The amount of liquid connecting the spheres is varied, and the influence of the resulting liquid-morphology on their dynamics is investigated. It is shown that the method is suitable for a study of liquid-phase sintering which includes both phase transformation and capillary driven rigid body motion.

preprint2019arXiv

A crossover in spatio-temporal correlations of strain fluctuations in glass forming liquids

Via molecular dynamics simulations of a generic glass former in the supercooled and normal liquid states, it is shown that spatial correlations of strain fluctuations exhibit a crossover from the well-established power-law $\sim 1/r^3$-decay at long wavelengths to an exponential behavior, $\sim \exp(-r/l_{\text {c}})$ at intermediate distances. The characteristic length of the exponential decay grows both with temperature and time via, $l_{\text {c}}^2 \propto D(T)\, t$, with $D(T)$ being the temperature-dependent diffusion coefficient. This suggests that the crossover between the power-law and exponential decays is governed by a diffusion process.

preprint2019arXiv

Controlling bubble coalescence in metallic foams: A simple phase field-based approach

The phase-field method is used as a basis to develop a strictly mass conserving, yet simple, model for simulation of two-phase flow. The model is aimed to be applied for the study of structure evolution in metallic foams. In this regard, the critical issue is to control the rate of bubble coalescence compared to concurrent processes such as their rearrangement due to fluid motion. In the present model, this is achieved by tuning the interface energy as a free parameter. The model is validated by a number of benchmark tests. First, stability of a two dimensional bubble is investigated by the Young-Laplace law for different values of the interface energy. Then, the coalescence of two bubbles is simulated until the system reaches equilibrium with a circular shape. To address the major capability of the present model for the formation of foam structure, the bubble coalescence is simulated for various values of interface energy in order to slow down the merging process. These simulations are repeated in the presence of a rotational flow to highlight the fact that the model allows to suppress the coalescence process compared to the motion of bubbles relative to each other.

preprint2019arXiv

Diffuse interface models of solidification with convection: The choice of a finite interface thickness

The thin interface limit aims at minimizing the effects arising from a numerical interface thickness, inherent in diffuse interface models of solidification and microstructure evolution such as the phase field model. While the original formulation of this problem is restricted to transport by diffusion, we consider here the case of melt convection. Using an analysis of the coupled phase field-fluid dynamic equations, we show here that such a thin interface limit does also exist if transport contains both diffusion and convection. This prediction is tested by comparing simulation studies, which make use of the thin-interface condition, with an analytic sharp-interface theory for dendritic tip growth under convection.

preprint2019arXiv

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules [The Journal of Chemical Physics 150 (2019) 024903]. In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere (HS) mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided.

preprint2013arXiv

Interfacial roughening in non-ideal fluids: Dynamic scaling in the weak- and strong-damping regime

Interfacial roughening denotes the nonequilibrium process by which an initially flat interface reaches its equilibrium state, characterized by the presence of thermally excited capillary waves. Roughening of fluid interfaces has been first analyzed by Flekkoy and Rothman [Phys. Rev. Lett. 75, 260 (1995)], where the dynamic scaling exponents in the weakly damped case in two dimensions were found to agree with the Kardar-Parisi-Zhang universality class. We extend this work by taking into account also the strong-damping regime and perform extensive fluctuating hydrodynamics simulations in two dimensions using the Lattice Boltzmann method. We show that the dynamic scaling behavior is different in the weakly and strongly damped case.

preprint2011arXiv

Dynamics of cylindrical droplets on flat substrate: Lattice Boltzmann modeling versus simple analytic models

The steady state motion of cylindrical droplets under the action of external body force is investigated both theoretically and via lattice Boltzmann simulation. As long as the shape-invariance of droplet is maintained, the droplet's center-of-mass velocity linearly scales with both the force density and the square of droplet radius. However, a non-linear behavior appears as the droplet deformation becomes significant. This deformation is associated with the drop elongation occurring at sufficiently high external forcing. Yet, independent of either the force density or the droplet size, the center-of-mass velocity is found to be linear in terms of the inverse of dynamic viscosity. In addition, it is shown that the energy is mainly dissipated in a region near the substrate particularly close to the three phase contact line. The total viscous dissipation is found to be proportional to both the square of force density and the inverse of dynamic viscosity. Moreover, the dependence of the center-of-mass velocity on the equilibrium contact angle is investigated. A simple analytic model is provided reproducing the observed behavior.

preprint2010arXiv

Efficient and accurate simulations of deformable particles immersed in a fluid using a combined immersed boundary lattice Boltzmann finite element method

The deformation of an initially spherical capsule, freely suspended in simple shear flow, can be computed analytically in the limit of small deformations [D. Barthes-Biesel, J. M. Rallison, The Time-Dependent Deformation of a Capsule Freely Suspended in a Linear Shear Flow, J. Fluid Mech. 113 (1981) 251-267]. Those analytic approximations are used to study the influence of the mesh tessellation method, the spatial resolution, and the discrete delta function of the immersed boundary method on the numerical results obtained by a coupled immersed boundary lattice Boltzmann finite element method. For the description of the capsule membrane, a finite element method and the Skalak constitutive model [R. Skalak et al., Strain Energy Function of Red Blood Cell Membranes, Biophys. J. 13 (1973) 245-264] have been employed. Our primary goal is the investigation of the presented model for small resolutions to provide a sound basis for efficient but accurate simulations of multiple deformable particles immersed in a fluid. We come to the conclusion that details of the membrane mesh, as tessellation method and resolution, play only a minor role. The hydrodynamic resolution, i.e., the width of the discrete delta function, can significantly influence the accuracy of the simulations. The discretization of the delta function introduces an artificial length scale, which effectively changes the radius and the deformability of the capsule. We discuss possibilities of reducing the computing time of simulations of deformable objects immersed in a fluid while maintaining high accuracy.

preprint2010arXiv

Fall and rise of small droplets on rough hydrophobic substrates

Liquid droplets on patterned hydrophobic substrates are typically observed either in the Wenzel or the Cassie state. Here we show that for droplets of comparable size to the roughness scale an additional local equilibrium state exists, where the droplet is immersed in the texture, but not yet contacts the bottom grooves. Upon evaporation, a droplet in this state enters the Cassie state, leading to a qualitatively new self-cleaning mechanism. The effect is of generic character and is expected to occur in any hydrophobic capillary wetting situation where a spherical liquid reservoir is involved.

preprint2010arXiv

Inertia effects and stress accumulation in a constricted duct: A combined experimental and lattice Boltzmann study

We experimentally and numerically investigate the flow of a Newtonian fluid through a constricted geometry for Reynolds numbers in the range $0.1 - 100$. The major aim is to study non-linear inertia effects at larger Reynolds numbers (>10) on the shear stress evolution in the fluid. This is of particular importance for blood flow as some biophysical processes in blood are sensitive to shear stresses, e.g., the initialization of blood clotting. We employ the lattice Boltzmann method for the simulations. The conclusion of the predictions is that the peak value of shear stress in the constriction grows disproportionally fast with the Reynolds number which leads to a non-linear shear stress accumulation. As a consequence, the combination of constricted blood vessel geometries and large Reynolds numbers may increase the risk of undesired blood clotting.

preprint2010arXiv

Shear stress in lattice Boltzmann simulations

A thorough study of shear stress within the lattice Boltzmann method is provided. Via standard multiscale Chapman-Enskog expansion we investigate the dependence of the error in shear stress on grid resolution showing that the shear stress obtained by the lattice Boltzmann method is second order accurate. This convergence, however, is usually spoiled by the boundary conditions. It is also investigated which value of the relaxation parameter minimizes the error. Furthermore, for simulations using velocity boundary conditions, an artificial mass increase is often observed. This is a consequence of the compressibility of the lattice Boltzmann fluid. We investigate this issue and derive an analytic expression for the time-dependence of the fluid density in terms of the Reynolds number, Mach number and a geometric factor for the case of a Poiseuille flow through a rectangular channel in three dimensions. Comparison of the analytic expression with results of lattice Boltzmann simulations shows excellent agreement.