scholarly journals Level-set simulations of a 2D topological rearrangement in a bubble assembly: effects of surfactant properties

2018 ◽  
Vol 838 ◽  
pp. 222-247 ◽  
Author(s):  
A. Titta ◽  
M. Le Merrer ◽  
F. Detcheverry ◽  
P. D. M. Spelt ◽  
A.-L. Biance

A liquid foam is a dispersion of gas bubbles in a liquid matrix containing surface-active agents. Its flow involves the relative motion of bubbles, which switch neighbours during a so-called topological rearrangement of type 1 (T1). The dynamics of T1 events, as well as foam rheology, have been extensively studied, and experimental results point to the key role played by surfactants in these processes. However, the complex and multiscale nature of the system has so far impeded a complete understanding of the mechanisms involved. In this work, we investigate numerically the effect of surfactants on the rheological response of a 2D sheared bubble cluster. To do so, a level-set method previously employed for simulation of two-phase flow has been extended to include the effects of surfactants. The dynamical processes of the surfactants – diffusion in the liquid and along the interface, adsorption/desorption at the interface – and their coupling with the flow – surfactant advection and Laplace and Marangoni stresses at the interface – are all taken into account explicitly. Through a systematic study of the Biot, capillary and Péclet numbers that characterise the surfactant properties in the simulation, we find that the presence of surfactants can affect the liquid/gas hydrodynamic boundary condition (from a rigid-like situation to a mobile one), which modifies the nature of the flow in the volume from a purely extensional situation to a shear. Furthermore, the work done by surface tension (the 2D analogue of the work by pressure forces), resulting from surfactant and interface dynamics, can be interpreted as an effective dissipation, which reaches a maximum for a Péclet number of order unity. Our results, obtained at high liquid fraction, should provide a reference point, with which experiments and models of T1 dynamics and foam rheology can be compared.

Author(s):  
Mamta Raju Jotkar ◽  
Daniel Rodriguez ◽  
Bruno Marins Soares

2017 ◽  
Vol 34 (3) ◽  
pp. 709-724 ◽  
Author(s):  
Amirmahdi Ghasemi ◽  
R. Nikbakhti ◽  
Amirreza Ghasemi ◽  
Faraz Hedayati ◽  
Amir Malvandi

Purpose A numerical method is developed to capture the interaction of solid object with two-phase flow with high density ratios. The current computational tool would be the first step of accurate modeling of wave energy converters in which the immense energy of the ocean can be extracted at low cost. Design/methodology/approach The full two-dimensional Navier–Stokes equations are discretized on a regular structured grid, and the two-step projection method along with multi-processing (OpenMP) is used to efficiently solve the flow equations. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. The full two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the flow field. Level set and immersed boundary methods are used to capture the free surface of liquid and solid object, respectively. A proper contact angle between the solid object and the fluid is used to enhance the accuracy of the advection of the mass and momentum of the fluids in three-phase cells. Findings The computational tool is verified based on numerical and experimental data with two scenarios: a cylinder falling into a rectangular domain due to gravity and a dam breaking in the presence of a fixed obstacle. In the former validation simulation, the accuracy of the immersed boundary method is verified. However, the accuracy of the level set method while the computational tool can model the high-density ratio is confirmed in the dam-breaking simulation. The results obtained from the current method are in good agreement with experimental data and other numerical studies. Practical/implications The computational tool is capable of being parallelized to reduce the computational cost; therefore, an OpenMP is used to solve the flow equations. Its application is seen in the following: wind energy conversion, interaction of solid object such as wind turbine with water waves, etc. Originality/value A high efficient CFD approach method is introduced to capture the interaction of solid object with a two-phase flow where they have high-density ratio. The current method has the ability to efficiently be parallelized.


2008 ◽  
Vol 11 (4-6) ◽  
pp. 221-235 ◽  
Author(s):  
S. P. van der Pijl ◽  
A. Segal ◽  
C. Vuik ◽  
P. Wesseling

2014 ◽  
Vol 100 ◽  
pp. 138-154 ◽  
Author(s):  
Lanhao Zhao ◽  
Jia Mao ◽  
Xin Bai ◽  
Xiaoqing Liu ◽  
Tongchun Li ◽  
...  

2017 ◽  
Vol 14 (06) ◽  
pp. 1750063 ◽  
Author(s):  
A. M. Hegab ◽  
S. A. Gutub ◽  
A. Balabel

This paper presents the development of an accurate and robust numerical modeling of instability of an interface separating two-phase system, such as liquid–gas and/or solid–gas systems. The instability of the interface can be refereed to the buoyancy and capillary effects in liquid–gas system. The governing unsteady Navier–Stokes along with the stress balance and kinematic conditions at the interface are solved separately in each fluid using the finite-volume approach for the liquid–gas system and the Hamilton–Jacobi equation for the solid–gas phase. The developed numerical model represents the surface and the body forces as boundary value conditions on the interface. The adapted approaches enable accurate modeling of fluid flows driven by either body or surface forces. The moving interface is tracked and captured using the level set function that initially defined for both fluids in the computational domain. To asses the developed numerical model and its versatility, a selection of different unsteady test cases including oscillation of a capillary wave, sloshing in a rectangular tank, the broken-dam problem involving different density fluids, simulation of air/water flow, and finally the moving interface between the solid and gas phases of solid rocket propellant combustion were examined. The latter case model allowed for the complete coupling between the gas-phase physics, the condensed-phase physics, and the unsteady nonuniform regression of either liquid or the propellant solid surfaces. The propagation of the unsteady nonplanar regression surface is described, using the Essentially-Non-Oscillatory (ENO) scheme with the aid of the level set strategy. The computational results demonstrate a remarkable capability of the developed numerical model to predict the dynamical characteristics of the liquid–gas and solid–gas flows, which is of great importance in many civilian and military industrial and engineering applications.


1996 ◽  
Vol 270 (1) ◽  
pp. C313-C320 ◽  
Author(s):  
C. N. Mayne ◽  
H. Sutherland ◽  
J. C. Jarvis ◽  
S. J. Gilroy ◽  
A. J. Craven ◽  
...  

Chronic electrical stimulation of skeletal muscle at 10 Hz induces fast-to-slow fiber type transformation. Does a lower aggregate amount of activity lead to a less complete transformation, or does it produce the same transformation over a longer time course? We examined this question by subjecting adult rabbit tibialis anterior and extensor digitorum longus muscles to continuous stimulation at 2.5 Hz for 2-12 wk. Most of the fibers acquired the histochemical and immunocytochemical characteristics of type 2A, not type 1, fibers. There was a corresponding rise in oxidative activity, but this was accompanied by a marked decline in anaerobic glycolysis. The activities of hexokinase and 3-oxoacid CoA-transferase stopped increasing after 2 wk, glutamate oxaloacetate transaminase after 4 wk, and beta-hydroxyacyl-CoA dehydrogenase after 6 wk of stimulation. Succinate dehydrogenase, citrate synthase, lactate dehydrogenase, and creatine phosphokinase continued to change up to 12 wk of stimulation. Changes in enzyme activity were not as rapid or as marked as those observed for stimulation at 10 Hz, and none showed the typical two-phase response of oxidative enzyme activities to stimulation at 10 Hz. The latter may therefore be dependent on induction of type 1 myosin isoforms.


Sign in / Sign up

Export Citation Format

Share Document