Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA one dimensional modal approach for flows controlled by contact line motion
64
0
0.0
(
0
)
Added by
Mauro Sbragaglia Dr
Publication date
2017
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The hydrodynamics of a liquid-vapour interface in contact with an heterogeneous surface is largely impacted by the presence of defects at the smaller scales. Such defects introduce morphological disturbances on the contact line and ultimately determine the force exerted on the wedge of liquid in contact with the surface. From the mathematical point of view, defects introduce perturbation modes, whose space-time evolution is governed by the interfacial hydrodynamic equations of the contact line. In this paper we derive the response function of the contact line to such generic perturbations. The contact line response may be used to design simplified 1+1 dimensional models accounting for the complexity of interfacial flows coupled to nanoscale defects, yet offering a more tractable mathematical framework to include thermal fluctuations and explore thermally activated contact line motion through a disordered energy landscape.
rate research
Read More
Common modal decomposition techniques for flowfield analysis, data-driven modeling and flow control, such as proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD) are usually performed in an Eulerian (fixed) frame of reference with snapshots from measurements or evolution equations. The Eulerian description poses some difficulties, however, when the domain or the mesh deforms with time as, for example, in fluid-structure interactions. For such cases, we first formulate a Lagrangian modal analysis (LMA) ansatz by a posteriori transforming the Eulerian flow fields into Lagrangian flow maps through an orientation and measure-preserving domain diffeomorphism. The development is then verified for Lagrangian variants of POD and DMD using direct numerical simulations (DNS) of two canonical flow configurations at Mach 0.5, the lid-driven cavity and flow past a cylinder, representing internal and external flows, respectively, at pre- and post-bifurcation Reynolds numbers. The LMA is demonstrated for several situations encompassing unsteady flow without and with boundary and mesh deformation as well as non-uniform base flows that are steady in Eulerian but not in Lagrangian frames. We show that LMA application to steady nonuniform base flow yields insights into flow stability and post-bifurcation dynamics. LMA naturally leads to Lagrangian coherent flow structures and connections with finite-time Lyapunov exponents (FTLE). We examine the mathematical link between FTLE and LMA by considering a double-gyre flow pattern. Dynamically important flow features in the Lagrangian sense are recovered by performing LMA with forward and backward (adjoint) time procedures.
A hybrid parallel approach for fully resolved simulations of particle-laden flows in sediment transport is proposed. To overcome the challenges of load imbalance in the traditional domain decomposition method when encountering highly uneven distributions of particles in space, we develop a hybrid parallel approach adopting the domain decomposition method for the carrier phase and a mirror domain technique for the disperse phase. We modify the mirror domain technique originally developed for point particles to fully resolved particle simulations, which are more challenging since a finite-sized particle may be split into different subdomains; thus, more complex treatments of particle-fluid interactions are needed. By utilizing the mirror domain technique, in which each processor handles nearly the same number of particles regardless of the particle spatial distribution, excellent load balance is achieved. The present hybrid parallel approach also shows strong scalability and high parallel efficiency in a test of a fully resolved simulation case of sediment transport. Furthermore, a novel memory optimization method is proposed for spherical particles of equal size, which can substantially reduce the memory cost and enable the simulation of millions of fully resolved particles on a common highly parallel computing platform. Our code is validated by several benchmark cases, and the results show good agreement with experimental and computational data in the literature.
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.
We model the advective Marangoni spreading of insoluble surfactant at the free surface of a viscous fluid that is confined within a two-dimensional rectangular cavity. Interfacial deflections are assumed small, with contact lines pinned to the walls of the cavity, and inertia is neglected. Linearizing the surfactant transport equation about the equilibrium state allows a modal decomposition of the dynamics, with eigenvalues corresponding to decay rates of perturbations. Computation of the family of mutually orthogonal two-dimensional eigenfunctions reveals singular flow structures near each contact line, resulting in spatially oscillatory patterns of wall shear stress and a pressure field that diverges logarithmically. These singularities at a stationary contact line are associated with dynamic compression of the surfactant monolayer; we show how they can be regularized by weak surface diffusion. Their existence highlights the need for careful treatment in computations of unsteady advection-dominated surfactant transport in confined domains.
We present numerical simulations of three-dimensional thermal convective flows in a cubic cell at high Rayleigh number using thermal lattice Boltzmann (LB) method. The thermal LB model is based on double distribution function approach, which consists of a D3Q19 model for the Navier-Stokes equations to simulate fluid flows and a D3Q7 model for the convection-diffusion equation to simulate heat transfer. Relaxation parameters are adjusted to achieve the isotropy of the fourth-order error term in the thermal LB model. Two types of thermal convective flows are considered: one is laminar thermal convection in side-heated convection cell, which is heated from one vertical side and cooled from the other vertical side; while the other is turbulent thermal convection in Rayleigh-Benard convection cell, which is heated from the bottom and cooled from the top. In side-heated convection cell, steady results of hydrodynamic quantities and Nusselt numbers are presented at Rayleigh numbers of $10^6$ and $10^7$, and Prandtl number of 0.71, where the mesh sizes are up to $257^3$; in Rayleigh-Benard convection cell, statistical averaged results of Reynolds and Nusselt numbers, as well as kinetic and thermal energy dissipation rates are presented at Rayleigh numbers of $10^6$, $3times 10^6$, and $10^7$, and Prandtl numbers of 0.7 and 7, where the nodes within thermal boundary layer are around 8. Compared with existing benchmark data obtained by other methods, the present LB model can give consistent results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
