ترغب بنشر مسار تعليمي؟ اضغط هنا

Multicomponent Flow on Curved Surfaces: A Vielbein Lattice Boltzmann Approach

95   0   0.0 ( 0 )
 نشر من قبل Halim Kusumaatmaja
 تاريخ النشر 2019
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We develop and implement a novel lattice Boltzmann scheme to study multicomponent flows on curved surfaces, coupling the continuity and Navier-Stokes equations with the Cahn-Hilliard equation to track the evolution of the binary fluid interfaces. Standard lattice Boltzmann method relies on regular Cartesian grids, which makes it generally unsuitable to study flow problems on curved surfaces. To alleviate this limitation, we use a vielbein formalism to write down the Boltzmann equation on an arbitrary geometry, and solve the evolution of the fluid distribution functions using a finite difference method. Focussing on the torus geometry as an example of a curved surface, we demonstrate drift motions of fluid droplets and stripes embedded on the surface of a torus. Interestingly, they migrate in opposite directions: fluid droplets to the outer side while fluid stripes to the inner side of the torus. For the latter we demonstrate that the global minimum configuration is unique for small stripe widths, but it becomes bistable for large stripe widths. Our simulations are also in agreement with analytical predictions for the Laplace pressure of the fluid stripes, and their damped oscillatory motion as they approach equilibrium configurations, capturing the corresponding decay timescale and oscillation frequency. Finally, we simulate the coarsening dynamics of phase separating binary fluids in the hydrodynamics and diffusive regimes for tori of various shapes, and compare the results against those for a flat two-dimensional surface. Our lattice Boltzmann scheme can be extended to other surfaces and coupled to other dynamical equations, opening up a vast range of applications involving complex flows on curved geometries.



قيم البحث

اقرأ أيضاً

We propose a novel approach to the numerical simulation of thin film flows, based on the lattice Boltzmann method. We outline the basic features of the method, show in which limits the expected thin film equations are recovered and perform validation tests. The numerical scheme is applied to the viscous Rayleigh-Taylor instability of a thin film and to the spreading of a sessile drop towards its equilibrium contact angle configuration. We show that the Cox-Voinov law is satisfied, and that the effect of a tunable slip length on the substrate is correctly captured. We address, then, the problem of a droplet sliding on an inclined plane, finding that the Capillary number scales linearly with the Bond number, in agreement with experimental results. At last, we demonstrate the ability of the method to handle heterogenous and complex systems by showcasing the controlled dewetting of a thin film on a chemically structured substrate.
281 - M. Sbragaglia , H. Chen , X. Shan 2009
It is shown that the Shan-Chen (SC) model for non-ideal lattice fluids can be made compliant with a pseudo free-energy principle by simple addition of a gradient force, whose expression is uniquely specified in terms of the fluid density. This additi onal term is numerically shown to provide fairly negligible effects on the system evolution during phase-separation. To the best of our knowledge, these important properties of the SC model were not noted before. The approach developed in the present work is based on a continuum analysis: further extensions, more in line with a discrete lattice theory (X. Shan, {it Phys Rev E}, {bf 77} 066702 (2008)) can be envisaged for the future.
The capability to simulate a two-way coupled interaction between a rarefied gas and an arbitrary-shaped colloidal particle is important for many practical applications, such as aerospace engineering, lung drug deliver and semiconductor manufacturing. By means of numerical simulations based on the Direct Simulation Monte Carlo (DSMC) method, we investigate the influence of the orientation of the particle and rarefaction on the drag and lift coefficients, in the case of prolate and oblate ellipsoidal particles immersed in a uniform ambient flow. This is done by modelling the solid particles using a cut-cell algorithm embedded within our DSMC solver. In this approach, the surface of the particle is described by its analytical expression and the microscopic gas-solid interactions are computed exactly using a ray-tracing technique. The measured drag and lift coefficients are used to extend the correlations available in the continuum regime to the rarefied regime, focusing on the transitional and free-molecular regimes. The functional forms for the correlations for the ellipsoidal particles are chosen as a generalisation from the spherical case. We show that the fits over the data from numerical simulations can be extended to regimes outside the simulated range of $Kn$ by testing the obtained predictive model on values of $Kn$ that where not included in the fitting process, allowing to achieve an higher precision when compared with existing predictive models from literature. Finally, we underline the importance of this work in providing new correlations for non-spherical particles that can be used for point-particle Euler-Lagrangian simulations to address the problem of contamination from finite-size particles in high-tech mechanical systems.
We present a novel method for fluid structure interaction (FSI) simulations where an original 2nd-order curved space lattice Boltzmann fluid solver (LBM) is coupled to a finite element method (FEM) for thin shells. The LBM can work independently on a standard lattice in curved coordinates without the need for interpolation, re-meshing or an immersed boundary. The LBM distribution functions are transformed dynamically under coordinate change. In addition, force and momentum can be calculated on the nodes exactly in any geometry. Furthermore, the FEM shell is a complete numerical tool with implementations such as growth, self-contact and strong external forces. We show resolution convergent error for standard tests under metric deformation. Mass and volume conservation, momentum transfer, boundary-slip and pressure maintenance are verified through specific examples. Additionally, a brief deformation stability analysis is carried out. Next, we study the interaction of a square fluid flow channel to a deformable shell. Finally, we simulate a flag at moderate Reynolds number, air flow channel. The scheme is limited to small deformations of O(10%) relative to domain size, by improving its stability the method can be naturally extended to multiple applications without further implementations.
72 - Y. Yu , Q. Li , 2019
The lattice Boltzmann (LB) method has gained much success in a variety of fields involving fluid flow and/or heat transfer. In this method, the bounce-back scheme is a popular boundary scheme for treating nonslip boundaries. However, this scheme lead s to staircase-shaped boundaries for curved walls. Therefore many curved boundary schemes have been proposed, but mostly suffer from mass leakage at the curved boundaries. Several correction schemes have been suggested for simulating single-phase flows, but very few discussions or studies have been made for two-phase LB simulations with curved boundaries. In this paper, the performance of three well-known types of curved boundary schemes in two-phase LB simulations is investigated through modeling a droplet resting on a circular cylinder. For all of the investigated schemes, the results show that the simulated droplet rapidly evaporates under the nonslip and isothermal conditions, owing to the imbalance between the mass streamed out of the system by the outgoing distribution functions and the mass streamed into the system by the incoming distribution functions at each boundary node. Based on the numerical investigation, we formulate two modified mass-conservative curved boundary schemes for two-phase LB simulations. The accuracy of the modified curved boundary schemes and their capability of conserving mass in two-phase LB simulations are numerically demonstrated.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا