Computational fluid dynamics for improved bioreactor design and 3D culture

The complex relationship between the hydrodynamic environment and surrounding tissues directly impacts on the design and production of clinically useful grafts and implants. Tissue engineers have generally seen bioreactors as 'black boxes' within which tissue engineering constructs (TECs) are cultured. It is accepted that a more detailed description of fluid mechanics and nutrient transport within process equipment can be achieved by using computational fluid dynamics (CFD) technology. This review discusses applications of CFD for tissue engineering-related bioreactors -- fluid flow processes have direct implications on cellular responses such as attachment, migration and proliferation. We conclude that CFD should be seen as an invaluable tool for analyzing and visualizing the impact of fluidic forces and stresses on cells and TECs.

Computational fluid dynamics analysis of a flat plate photocatalytic reactor for storm and wastewater reuse

A computational fluid dynamics (CFD) analysis has been performed for a flat plate photocatalytic reactor using CFD code FLUENT. Under the simulated conditions (Reynolds number, Re around 2650), a detailed time accurate computation shows the different stages of flow evolution and the effects of finite length of the reactor in creating flow instability, which is important to improve the performance of the reactor for storm and wastewater reuse. The efficiency of a photocatalytic reactor for pollutant decontamination depends on reactor hydrodynamics and configurations. This study aims to investigate the role of different parameters on the optimization of the reactor design for its improved performance. In this regard, more modelling and experimental efforts are ongoing to better understand the interplay of the parameters that influence the performance of the flat plate photocatalytic reactor.

Computational fluid dynamics simulation and turbomachinery code validation of a high pressure ratio radial-inflow turbine

The present study explores reproducing the closest geometry of a high pressure ratio single stage radial-inflow turbine applied in the Sundstrans Power Systems T-100 Multipurpose Small Power Unit. The commercial software ANSYS-Vista RTD along with a built in module, BladeGen, is used to conduct a meanline design and create 3D geometry of one flow passage. Carefully examining the proposed design against the geometrical and experimental data, ANSYS-TurboGrid is applied to generate computational mesh. CFD simulations are performed with ANSYS-CFX in which three-dimensional Reynolds-Averaged Navier-Stokes equations are solved subject to appropriate boundary conditions. Results are compared with numerical and experimental data published in the literature in order to generate the exact geometry of the existing turbine and validate the numerical results against the experimental ones.

A nonconservative Lagrangian framework for statistical fluid registration - SAFIRA

In this paper, we used a nonconservative Lagrangian mechanics approach to formulate a new statistical algorithm for fluid registration of 3-D brain images. This algorithm is named SAFIRA, acronym for statistically-assisted fluid image registration algorithm. A nonstatistical version of this algorithm was implemented, where the deformation was regularized by penalizing deviations from a zero rate of strain. In, the terms regularizing the deformation included the covariance of the deformation matrices Σ and the vector fields (q). Here, we used a Lagrangian framework to reformulate this algorithm, showing that the regularizing terms essentially allow nonconservative work to occur during the flow. Given 3-D brain images from a group of subjects, vector fields and their corresponding deformation matrices are computed in a first round of registrations using the nonstatistical implementation. Covariance matrices for both the deformation matrices and the vector fields are then obtained and incorporated (separately or jointly) in the nonconservative terms, creating four versions of SAFIRA. We evaluated and compared our algorithms' performance on 92 3-D brain scans from healthy monozygotic and dizygotic twins; 2-D validations are also shown for corpus callosum shapes delineated at midline in the same subjects. After preliminary tests to demonstrate each method, we compared their detection power using tensor-based morphometry (TBM), a technique to analyze local volumetric differences in brain structure. We compared the accuracy of each algorithm variant using various statistical metrics derived from the images and deformation fields. All these tests were also run with a traditional fluid method, which has been quite widely used in TBM studies. The versions incorporating vector-based empirical statistics on brain variation were consistently more accurate than their counterparts, when used for automated volumetric quantification in new brain images. This suggests the advantages of this approach for large-scale neuroimaging studies.

Collaborative learning approach to introduce computational fluid dynamics

This paper presents the design, implementation and evaluation of a collaborative learning activity designed to replace traditional face-to-face lectures in a large classroom. This activity aims to better engage the students with their learning and improve the students’ experience and outcomes. This project is implemented in the Fluid Mechanics unit of the Mechanical Engineering degree at the Queensland University of Technology to introduce students with the concept, terminology and process of Computational Fluid Dynamics (CFD). The approach integrates a constructive collaborative assignment which is a key element in the overall quality of teaching and learning, and an integral component of the students’ experience. A detailed survey, given to the students, showed an overall high level of satisfaction. However, the results also highlighted the gap between students’ expectations both for contents and assignment and teacher expectations. Discussions to address this issue are presented in the paper based on a critical reflection.

Fluid dynamics of the monsoon

The monsoonal regions of the world are characterized by a seasonal reversal in the direction of winds associated with the excursion of the equatorial trough (or the ITCZ) in response to the variation in the latitude of maximum insolation. This monsoonal circulation is a planetary scale phenomenon. However, the associated precipitation is critically dependent on the organization of the cumulus clouds (typically a few kilometers in horizontal extent) over the scale of synoptic vortices (typically a few hundred kilometers in horizontal extent). Thus modelling of the seasonal transitions and intraseasonal fluctuations requires an understanding of the fluid mechanics of these three scales of organizations and their interactions. The present paper is an attempt to outline the current state of understanding of these phenomena.

A note on the application of a radiation condition for a source in a rotating stratified fluid

The motion due to an oscillatory point source in a rotating stratified fluid has been studied by Sarma & Naidu (1972) by using threefold Fourier transforms. The solution obtained by them in the hyperbolic case is wrong since they did not make use of any radiation condition, which is always necessary to get the correct solution. Whenever the motion is created by a source, the condition of radiation is that the sources must remain sources, not sinks of energy and no energy may be radiated from infinity into the prescribed singularities of the field. The purpose of the present note is to explain how Lighthill's (1960) radiation condition can be applied in the hyperbolic case to pick the correct solution. Further, the solution thus obtained is reiterated by an alternative procedure using Sommerfeld's (1964) radiation condition.

Computational Fluid Dynamics in Combustion

Computational fluid dynamics has reached a stage where flow field in practical situation can be predicted to aid the design and to probe into the fundamental flow physics to understand and resolve the issues in fundamental fluid mechanics The study examines the computation of reacting flows After exploring the conservation equations for species and energy, the methods of closing the reaction rate terms in turbulent flow have been examined briefly Two cases of computation where combustion-flow interaction plays important role, have been discussed to illustrate the computational aspects and the physical insight that can be gained by the reacting flow computation

Breakage of a drop of inviscid fluid due to a pressure fluctuation at its surface

In this work, an attempt is made to gain a better understanding of the breakage of low-viscosity drops in turbulent flows by determining the dynamics of deformation of an inviscid drop in response to a pressure variation acting on the drop surface. Known scaling relationships between wavenumbers and frequencies, and between pressure fluctuations and velocity fluctuations in the inertial subrange are used in characterizing the pressure fluctuation. The existence of a maximum stable drop diameter d(max) follows once scaling laws of turbulent flow are used to correlate the magnitude of the disruptive forces with the duration for which they act. Two undetermined dimensionless quantities, both of order unity, appear in the equations of continuity, motion, and the boundary conditions in terms of pressure fluctuations applied on the surface. One is a constant of proportionality relating root-mean-square values of pressure and velocity differences between two points separated by a distance l. The other is a Weber number based on turbulent stresses acting on the drop and the resisting stresses in the drop due to interfacial tension. The former is set equal to 1, and the latter is determined by studying the interaction of a drop of diameter equal to d(max) with a pressure fluctuation of length scale equal to the drop diameter. The model is then used to study the breakage of drops of diameter greater than d(max) and those with densities different from that of the suspending fluid. It is found that, at least during breakage of a drop of diameter greater than d(max) by interaction with a fluctuation of equal length scale, a satellite drop is always formed between two larger drops. When very large drops are broken by smaller-length-scale fluctuations, highly deformed shapes are produced suggesting the possibility of further fragmentation due to instabilities. The model predicts that as the dispersed-phase density increases, d(max) decreases.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Stability of fluid flow past a membrane

The stability of fluid flow past a membrane of infinitesimal thickness is analysed in the limit of zero Reynolds number using linear and weakly nonlinear analyses. The system consists of two Newtonian fluids of thickness R* and H R*, separated by an infinitesimally thick membrane, which is flat in the unperturbed state. The dynamics of the membrane is described by its normal displacement from the flat state, as well as a surface displacement field which provides the displacement of material points from their steady-state positions due to the tangential stress exerted by the fluid flow. The surface stress in the membrane (force per unit length) contains an elastic component proportional to the strain along the surface of the membrane, and a viscous component proportional to the strain rate. The linear analysis reveals that the fluctuations become unstable in the long-wave (alpha --> 0) limit when the non-dimensional strain rate in the fluid exceeds a critical value Lambda(t), and this critical value increases proportional to alpha(2) in this limit. Here, alpha is the dimensionless wavenumber of the perturbations scaled by the inverse of the fluid thickness R*(-1), and the dimensionless strain rate is given by Lambda(t) = ((gamma) over dot* R*eta*/Gamma*), where eta* is the fluid viscosity, Gamma* is the tension of the membrane and (gamma) over dot* is the strain rate in the fluid. The weakly nonlinear stability analysis shows that perturbations are supercritically stable in the alpha --> 0 limit.

Flow of a stratified fluid in a wavy channel

The flow of a stratified fluid in a channel with small and large deformations is investigated. The analogy of this flow with swirling flow in tubes with non-uniform cross-sections is studied. The flow near the wall is blocked when the Froude number takes certain critical values. The possibility of preventing the stagnation zones in the flow field is also discussed

On an oscillatory point force in a rotating stratified fluid

The forced oscillations due to a point forcing effect in an infinite or contained, inviscid, incompressible, rotating, stratified fluid are investigated taking into account the density variation in the inertia terms in the linearized equations of motion. The solutions are obtained in closed form using generalized Fourier transforms. Solutions are presented for a medium bounded by a finite cylinder when the oscillatory forcing effect is acting at a point on the axis of the cylinder. In both the unbounded and bounded case, there exist characteristic cones emanating from the point of application of the force on which either the pressure or its derivatives are discontinuous. The perfect resonance existing at certain frequencies in an unbounded or bounded homogeneous fluid is avoided in the case of a confined stratified fluid.

A note on the damping and oscillations of a fluid drop moving in another fluid

The oscillations of a drop moving in another fluid medium have been studied at low values of Reynolds number and Weber number by taking into consideration the shape of the drop and the viscosities of the two phases in addition to the interfacial tension. The deformation of the drop modifies the Lamb's expression for frequency by including a correction term while the viscous effects split the frequency into a pair of frequencies—one lower and the other higher than Lamb's. The lower frequency mode has ample experimental support while the higher frequency mode has also been observed. The two modes almost merge with Lamb's frequency for the asymptotic cases of a drop in free space or a bubble in a dense viscous fluid but the splitting becomes large when the two fluids have similar properties. Instead of oscillations, aperiodic damping modes are found to occur in drops with sizes smaller than a critical size ($\sim\hat{\rho}\hat{\nu}^2/T $). With the help of these calculations, many of the available experimental results are analyzed and discussed.