I. Multi-Phase Flows: Computational Multi-Fluid Dynamics (CMFD) Development, Application and Analysis
by D. Datta, Vinesh H. Gada, Saurish Das, Vishesh Aggarwal, Absar Lakdawala, Javed Shaikh and Nagesh Patil
Level Set Method (LSM) is a recent numerical method for capturing transient evolution of interfaces, encountered in various branches of science and engineering. As it was proposed by a mathematician, we were struggling with the physical interpretation of the function used in the method. In the LSM, there are two sets of equations: Navier-Stokes for flow-properties and level-set for interface. In our first work [1], we came up with a physical interpretation and used them in control volume based derivation (commonly used in fluid-mechanics) of the equations used in level set method.
In our second work [2], an improvement in the LSM was proposed to improve the mass-error; the biggest disadvantage of LSM. It was proposed to used finer grid for level-set and coarser grid for Navier-Stokes equations; called as Dual Resolution LSM (DR-LSM). Since the interface-equations are advection equations as compared to non-linear coupled flow-equations, the computational time for the level-set as compared to Navier-Stokes equations are much smaller. The slight increase in the computational cost - due to fine grid for LS equations - for DR-LSM as compared to LSM was envisaged [2] to greatly benefit in terms of substantial improvement the accuracy/mass-error of the result. Thus, a built-in grid globally-refined DR-LSM (one additional grid for level set function in-between two uniformly spaced grid points for velocity/pressure/temperature everywhere in the domain) was proposed for Cartesian [2] and Cylindrical [3] coordinate system; called as Dual-Grid Level-Set Method (DGLSM). It takes slightly more computational time as compared to coarse-grid LSM to achieve a computational accuracy slightly less than that obtained by fine-grid LSM. This was demonstrated for various two-phase flow problems in Cartesian Coordinate [2]; and for two fluid electrodynamic axi-symmetric flow [4]. (Refer ANIMATIONS).
DGLSM was also presented for 2D [2] and 3D [3] film boiling. LSM for casting solidification was developed [5] and a level set based Eulerian–Lagrangian technique was proposed for computation of feed-path - to determine hot-spots. Further modification of DGLSM was presented for two-phase electro-hydrodynamic flow [4] and for non-newtonian two-fluid flow [6]. For contact line motion due to impact dynamics of a droplet on a solid surface, DGLSM was presented recently [7]. Recently, a LSM based parallelization methodology and parallel performance study on five different 3D transient test-cases - corresponding to single/two phase flow, with/without phase change and Cartesian/Cylindrical coordinate system - was presented [8].
LSM based Computational Multi-Fluid Dynamics (CMFD) application and analysis was done for developing smooth-stratified flow in a horizontal channel and a square duct [9]. This was also done for developing oil–water smooth-stratified (SS) and wavy-stratified (WS) flow in a horizontal and inclined plane-channel at various inlet-velocity, inlet-interface-height, inclination-angle, reduced surface-tension and reduced gravity. For smooth-stratified flow in a pipe, a LSM based non-dimensional 3D study is done [11] at various viscosity-ratio, inlet area-fraction of the less-viscous fluid and Reynolds number. For two-phase flow with heat transfer, a DGLSM based study was done [12] for developing two-fluid stratified flow in a plane channel - subjected to various thermal boundary conditions - with and without phase change. DGLSM based CMFD application and analysis was also done for injection of a liquid jet injected upwards into another stationary or co-flowing immiscible liquid [13] at six different combination of the dispersed and continuous fluid and various injection velocity. The above sudy on upward injection of Newtonian jet was extended to downward injection of non-Newtonian jet recently [6]. DGLSM based study was also done for impact-dynamics of bouncing and non-bouncing droplet on superhydrophobic substrates was also studied recently [7]. DGLSM based 2D transient simulation of saturated pool film-boiling on horizontal surface was done for single as well as multi-mode film boiling [3, 14]. 3D transient conjugate heat-transfer simulations for film-boiling in horizontal heat generating rod was done [3] and a transition from synchronous to quasi-periodic bubble-release regime - at a certain critical value of Jacob number - was reported. 2D film-condensation under a horizontal plate was simulated and flow transition from dripping to jetting mode was presented, with increasing degree of sub-cooling of the plate [15]. (Refer ANIMATIONS).
II. Stationary and Moving Bluff-Body Flows: CFD Development, Application and Analysis
by Sachin B. Paramane, Kaushal Prasad and Mukul Srivastava
A Multi-block structured finite volume method based NS solver was developed and effect of various convection schemes was studied [1]. The code is used for unconfined [2-4] and channel-confined [5] flow and heat transfer across a rotating circular cylinder. For the cylinder maintained at constant-wall-temperature [2] and constant heat flux [3], the heat transfer is discussed with the help of heat-lines; proposed for the first time for complex geometry CFD problems. Under the influence of cross-stream buoyancy [4], origin of secondary frequency with the shedding of vortex-pair of two different sizes alternatively is reported at Re=40 and 100. Rotating and channel-confinement is found to have stabilizing effect on the flow [5]. Positive as compared to negative eccentricity - in the position of cylinder in a channel- is found to have more stabilization effect, for counter-clockwise rotating cylinder [6].
The physical principles underlying the extraordinary mobility of swimming and flying animals have been the subject of years of research and there is still much that is not understood. Fishes are smart swimmers which use jet-stream propulsion thereby achieving higher propulsive efficiency as compared to man made machines. The maneuverability and efficiency of fish is inspiring new styles of propulsion and maneuvering in Autonomous Underwater Vehicles (AUVs) for applications such as under-sea exploration and environmental monitoring. The unresolved issues pertain to the mechanisms of propulsion and the mechanical efficiency of self-propelling organisms. Experimental approaches are very difficult, since instrumenting these animals often disrupts the behavior of interest. An opportunity exists to bring modern numerical methods in fluid dynamics to bear on these problems. With this motivation, a level set based immersed boundary method was proposed and an in-house code was developed as well as tested for various moving boundary problems [7]. The code is used for CFD application and analysis of flish like locomotion. 2D simulations are done [8] for flow generated by translating and undulating fish-like body - modeled as NACA 0012 hydrofoil. Shape of hydrofoil needs to be transformed into fish-like body by adding time varying undulations at various X locations from the head of hydrofoil along the axis. A detailed parametric investigation - for various flow, geometric as well as locomotion governing flow parameters - is under-progress for the CFD simulation on hydrodynamics of single as well as multiple (tandem, side-by-side and school arrangement) fish like locomotion (Refer ANIMATIONS).