Criteria Of Turbulent Transition In Parallel Flows

Based on the energy gradient method, criteria for turbulent transition are proposed for pressure driven flow and shear driven flow, respectively. For pressure driven flow, the necessary and sufficient condition for turbulent transition is the presence of the velocity inflection point in the averaged flow. For shear driven flow, the necessary and sufficient condition for turbulent transition is the existence of zero velocity gradient in the averaged flow profile. It is shown that turbulent transition can be effected via a singularity of the energy gradient function which may be associated with the chaotic attractor in dynamic system. The role of disturbance in the transition is also clarified in causing the energy gradient function to approach the singularity. Finally, it is interesting that turbulence can be controlled by modulating the distribution of the energy gradient function.Comment: The Third International Symposium on Physics of Fluids (ISPF3) 15-18 June, 2009, Jiuzhaigou, China. 6 pages; 1 figur

A Conservative Front Tracking Algorithm

The discontinuities in the solutions of systems of conservation laws are widely considered as one of the difficulties in numerical simulation. A numerical method is proposed for solving these partial differential equations with discontinuities in the solution. The method is able to track these sharp discontinuities or interfaces while still fully maintain the conservation property. The motion of the front is obtained by solving a Riemann problem based on the state values at its both sides which are reconstructed by using weighted essentially non oscillatory (WENO) scheme. The propagation of the front is coupled with the evaluation of "dynamic" numerical fluxes. Some numerical tests in 1D and preliminary results in 2D are presented.Singapore-MIT Alliance (SMA

Numerical comparison between Maxwell stress method and equivalent multipole approach for calculation of the dielectrophoretic force in octupolar cell traps

This work presents detailed numerical calculations of the dielectrophoretic force in octupolar traps designed for single-cell trapping. A trap with eight planar electrodes is studied for spherical and ellipsoidal particles using an indirect implementation of the boundary element method (BEM). Multipolar approximations of orders one to three are compared with the full Maxwell stress tensor (MST) calculation of the electrical force on spherical particles. Ellipsoidal particles are also studied, but in their case only the dipolar approximation is available for comparison with the MST solution. The results show that the full MST calculation is only required in the study of non-spherical particles.Singapore-MIT Alliance (SMA

An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains

We present an immersed interface method for the incompressible Navier Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular forces are applied on the fluid at the immersed boundary. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines. The strength of singular forces is determined by solving a small system of equations at each time step. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method for pressure and velocity.Singapore-MIT Alliance (SMA

An Immersed Interface Method for the Incompressible Navier-Stokes Equations

We present an immersed interface algorithm for the incompressible Navier Stokes equations. The interface is represented by cubic splines which are interpolated through a set of Lagrangian control points. The position of the control points is implicitly updated using the fluid velocity. The forces that the interface exerts on the fluid are computed from the constitutive relation of the interface and are applied to the fluid through jumps in the pressure and jumps in the derivatives of pressure and velocity. A projection method is used to time advance the Navier-Stokes equations on a uniform cartesian mesh. The Poisson-like equations required for the implicit solution of the diffusive and pressure terms are solved using a fast Fourier transform algorithm. The position of the interface is updated implicitly using a quasi-Newton method (BFGS) within each timestep. Several examples are presented to illustrate the flexibility of the presented approach.Singapore-MIT Alliance (SMA

Bubble Simulation Using Level Set-Boundary Element Method

In bubble dynamics, an underwater bubble may evolve from being singly-connected to being toroidal. Furthermore, two or more individual bubbles may merge to form a single large bubble. These dynamics involve significant topological changes such as merging and breaking, which may not be handled well by front-tracking boundary element methods. In the level set method, topological changes are handled naturally through a higher-dimensional level set function. This makes it an attractive method for bubble simulation. In this paper, we present a method that combines the level set method and the boundary element method for the simulation of bubble dynamics. We propose a formulation for the update of a potential function in the level set context. This potential function is non-physical off the bubble surface but consistent with the physics on the bubble surface. We consider only axisymmetric cavitation bubbles in this paper. Included in the paper are some preliminary results and findings.Singapore-MIT Alliance (SMA

Acceleration of Gas Bubble-Free Surface Interaction Computation Using Basis Preconditioners

The computation of gas bubble-free surface interaction entails a time-stepping algorithm whereby a linear system is solved at each time-iteration. In our investigation, the linear systems are derived from a desingularized boundary integral formulation and are poorly conditioned. This leads to poor convergence rates when Krylov subspace methods are used to solve these systems. The convergence rates may however be improved with proper preconditioning. We limit our investigation to gas bubbles initiated at depths sufficiently small such that a spike forms on the free surface during the later stages of evolution. Bubble dynamics dictate that for gas bubbles initiated at such depths, the stages through which the gas bubble and free surface evolve are similar. Based on this fact, we propose to perform one computation run for a gas bubble initiated at one particular depth, obtain a judicious set of a priori basis preconditioners from this run and thereafter, use this set of preconditioners on computation runs for gas bubble initiated at different depths. The computation time taken by the proposed method is, in general, 50% and 20% of the time taken by the present method (without preconditioning) with terminating criteria of 1.0e-5 and 1.0e-7 in the infinity-norm respectively using the Bi-conjugate Gradient Stabilized solver. The present method further enables computation to an infinity-norm terminating criterion of 1.0e-10 in a shorter time compared to the present method with a criterion of 1.0e-5.Singapore-MIT Alliance (SMA

Model Order Reduction for Determining Bubble Parameters to Attain a Desired Fluid Surface Shape

In this paper, a new methodology for predicting fluid free surface shape using Model Order Reduction (MOR) is presented. Proper Orthogonal Decomposition combined with a linear interpolation procedure for its coefficient is applied to a problem involving bubble dynamics near to a free surface. A model is developed to accurately and efficiently capture the variation of the free surface shape with different bubble parameters. In addition, a systematic approach is developed within the MOR framework to find the best initial locations and pressures for a set of bubbles beneath the quiescent free surface such that the resultant free surface attained is close to a desired shape. Predictions of the free surface in two-dimensions and three-dimensions are presented.Singapore-MIT Alliance (SMA