TOPOLOGY OPTIMIZATION OF NAVIER-STOKES FLOW IN MICROFLUIDICS
|
|
- Phyllis Andrews
- 5 years ago
- Views:
Transcription
1 European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, July 2004 TOPOLOGY OPTIMIZATION OF NAVIER-STOKES FLOW IN MICROFLUIDICS Laurits Højgaard Olesen, Fridolin Okkels, and Henrik Bruus MIC - Department of Micro and Nanotechnology, Technical University of Denmark DK-2800 Kongens Lyngby, Denmark Key words: Topology optimization, microfluidics, inertial effects, Femlab Abstract. We study how the inertial term in the steady state Navier-Stokes flow can drive transitions in the optimal topology of a microfluidic network. Our method is a generalization of the work by Borrvall and Petersson on topology optimization of fluids in Stokes flow [Int. J. Num. Meth. Fluids 41, 77 (2003)]. We state the Navier-Stokes flow problem in a generic divergence form that allows the problem to be solved using of the commercial finite element software package Femlab. Moreover we show how Femlab can be used as a high-level programming language that allows a wide range of optimization objectives to be easily dealt with. 1
2 u max γ 2 l γ 1 2l L 2l Figure 1: Schematic illustration of the flow-problem. Two inlets and two outlets (white areas) of width l and length 2l lead to the design region (gray area) of width 5l and length L. The maximal velocity of the fluid at the inlet is denoted u max, while its density and viscosity are ρ and η, respectively. The flow is characterized by the Reynolds number Re = ρ l u max /η. 1 Introduction Topology optimization has been used for many physical systems, such as structural mechanics, optics, and acoustics [1, 2, 3, 4], and recently Borrvall and Petersson introduced the method for fluids in Stokes flow [5]. However, it is desirable to extend the method to fluids described in a full Navier-Stokes flow; a direction pioneered by the work of Sigmund and Gersborg-Hansen [6, 7]. In Section 2 we formulate the full incompressible Navier-Stokes equation for flows through a porous medium. The topology optimization is introduced with the goal of finding the microfluidic network with the least power dissipation for the given boundary conditions. In Section 3 we cast the entire problem into a generic divergence form that encompasses a wide range of coupled non-linear second order partial differential equation problems. This form is the one used in Femlab. In Section 4 we work out the sensitivity analysis in such a way that the built-in symbolic differentiation tools of Femlab can be exploited. In Section 5 we introduce our model system, a four-terminal microfluidic device, and in Section 6 we present the results. We find that there is a transition in the optimal topology of the microfluidic device when the Reynolds number Re is increased. Finally in Section 7 we conclude. 2 Topology optimization for Navier-Stokes flow in steady state We consider a given computational domain Ω with boundary conditions for the flow given on the domain boundary Ω. The goal of the optimization is to distribute a certain amount of solid material inside Ω such that the material defines a microfluidic device or network that is optimal in the sense of minimizing the dissipated energy inside the system for a prescribed flow rate. 2
3 Following Ref. [5] the discrete transition from solid wall to open channel is made continuous by assuming the design domain Ω to be filled with a porous medium of spatially varying pore size. Then solid wall and open channels correspond to the limits of very small or very large pore sizes, respectively. The fluid in the porous medium is subject to a friction force f, c.f. Darcy s law, which is proportional to the fluid velocity f = αv, where v is the fluid velocity, and α(r) is the inverse of the local permeability of the medium at position r. Note that permeability will depend strongly on the local pore size. The flow problem is described in terms of the fluid velocity field v(r) and pressure p(r). The governing equations are the steady state Navier-Stokes equation and the incompressibility constraint ρ(v )v = σ αv, (1) v = 0, (2) where ρ is the mass density of the fluid. For an incompressible Newtonian fluid the Cauchy stress tensor σ is given by ( vi σ ij = p δ ij + η + v ) j, (3) x j x i where η is the dynamic viscosity. For simplicity we consider here only two-dimensional problems, i.e., we assume translational invariance in the third dimension and set r = (x 1, x 2 ) and v = ( v 1 (r), v 2 (r) ). In the optimization process, the material distribution is controlled by introducing a design variable field γ(r) describing the local porosity of the medium. We allow γ to vary between zero in the regions of vanishing pore size and unity in the open channels. Following Ref. [5] we relate the inverse permeability α to the porosity γ as α(r) α max + (α min α max ) γ(r) q + 1 q + γ(r), (4) where α min and α max denote the minimum and maximum values of α obtained when γ = 1 and γ = 0, respectively. We use α min = 0 and α max = 10 3, while q is a real and positive parameter. As discussed in [5], choosing a nonzero α min the problem Eqs. (1) can be thought of as modelling a plane flow between parallel plates with a slowly varying separation h(r). In this case α(r) = 8η/h(r) 2 and α min = 8η/h 2 max. The total power dissipation inside the fluidic system is given by [8] Φ(v, p, γ) = Ω ( i,j σ ij v i x j + i α(γ)v 2 i ) dr. Notice that in steady state this must equal the power delivered to the system, and it is simply ( Φ(v, p, γ) = n σ v n v [ 1 ρv2]) ds, (5b) 2 Ω 3 (5a)
4 where n is the unit outward normal vector, such that n σ v is the external work done on the system, while n v( 1 2 ρv2 ) is the kinetic energy convected into the system. In the common case where the boundary conditions are chosen such that on the external solid walls v = 0, while on the inlet and outlet boundaries v is parallel to n and (n ) v = 0, Eq. (5b) reduces to Φ(v, p, γ) = ( n) v [ p + 1ρv2] ds. (5c) 2 Ω This makes it clear that minimizing the dissipated power at fixed flow rate is equivalent to minimizing the pressure drop across the system. Conversely, if the boundary conditions fix the pressure drop, then the natural design objective will be to maximize the flow rate or equivalently to maximize the dissipated power. 3 Formulation of the flow problem in divergence form In this section we cast the governing equations Eqs. (1) and (2) into a generic divergence form that encompasses a wide range of coupled second order non-linear partial differential equation problems. In particular, this form is the one used to specify a partial differential equation problem in strong form in Femlab. We first introduce the velocity-pressure vector u = [v 1, v 2, p ] and define for i = 1, 2, 3 the quantities Γ i and F i as and Γ 1 [ σ11 σ 21 ], Γ 2 [ σ12 σ 22 ], Γ 3 [ ] 0 0 F 1 ρ(v )v 1 + αv 1, F 2 ρ(v )v 2 + αv 2, F 3 v. (7) Using this, Eqs. (1) and (2) can be written as Γ i = F i in Ω, (8a) R i = 0 on Ω, (8b) R j n Γ i = G i + µ j on Ω, (8c) u i where Γ i and F i are functions of the solution u, its gradient u, and of the design variable γ. The quantities R i (u, γ) in Eq. (8b) describe Dirichlet type boundary conditions. For example, fluid no-slip boundary conditions are obtained by defining R 1 v 1 and R 2 v 2 on the external solid walls. The quantities G i (u, γ) in Eq. (8c) describe Neumann type boundary conditions, and µ i denote the Lagrange multiplier necessary to enforce the constraint R i = 0, e.g., the force with which the solid wall has to act upon the fluid to enforce the no-slip boundary condition. (6) 4
5 The optimal design problem can now be stated as a constrained non-linear optimization problem: min γ imize Φ(u, γ) (9a) subject to Γ i = F i The three governing equations, (9b) Ω R i = 0 n Γ i = G i + R j u i µ j 0 γ 1 Design variable bounds, (9c) γ(r) dr βv Ω 0 Volume constraint. (9d) where V Ω is the volume of the design domain and β is the largest fraction of the volume that may be occupied by open channels. This optimization problem can be solved using several different approaches, and we have chosen the popular method of moving asymptotes (MMA) [9]. The MMA is designed for problems with a large number of variables such as ours. It is a gradient based algorithm and thus needs information about the derivative with respect to γ of both the objective function Φ and the constraints. Notice that for any γ the governing equations allow us to solve for u, so in fact they define u(γ) as an implicit function. Therefore using the chain rule we obtain [ ] Φ(u(γ), γ) = Φ + Ω Φ u i dr. (10) u i The complication here is that changes in γ(r) generally results in non-local changes in the implicit function u(γ). In other words the matrix u(r )/(r) is not diagonal but dense. This problem is dealt with using the standard adjoint method [10] to eliminate u/ from Eq. (10) as discussed below. The optimization process is iterative. The kth iteration consists of three steps: (I) For a given guess for the optimal material distribution γ (k) we solve Eqs. (8a) to (8c) for u(γ) as a finite element problem using Femlab. (II) Next we perform a sensitivity analysis by solving the adjoint problem to eliminate u/ from Eq. (10). This is also done using Femlab. (III) Finally we use MMA to obtain a new guess for the optimal design γ (k+1) based on the gradient information and the past iteration history. Of those, step (I) is the most expensive since it involves the solution of a non-linear partial differential equation. 5
6 4 Discretization and sensitivity analysis We solve the incompressible Navier-Stokes flow problem with the finite element method using the commercial software package Femlab. In this section we discuss how the generic problem Eqs. (8a) to (8c) is discretized and solved in Femlab. We also work out the sensitivity analysis for such generic problems. The starting point is to approximate the solution component u i using a set of finite element basis functions {ϕ i,n (r)}, u i (r) = n u i,n ϕ i,n (r), (11) where u i,n are the expansion coefficients. We have used the standard Taylor-Hood element pair with piecewise quadratic velocity approximation and piecewise linear pressure. Similiarly, the design variable field γ(r) is expressed as γ(r) = n γ n ϕ 4,n (r), (12) where we have chosen the basis functions to be continuous and piecewise linear. The problem of Eqs. (8a) to (8c) is discretized by the standard Galerkin method and takes the form L i (U, γ) N T jiλ j = 0, (13a) M i (U, γ) = 0, (13b) where U i, Λ i, and γ are column vectors holding all the expansion coefficients for the solution u i,n, the Lagrange multipliers µ i,n, and the design variable field γ n, respectively. The components of the column vector L i contain the projection of Eq. (8a) onto ϕ i,n and are given by L i,n = ϕ i,n F i + ϕ i,n Γ i dr + ϕ i,n G i ds, (14) Ω whereas the column vector M i holds the pointwise enforcement of the Dirichlet constraint Eq. (8b) given by M i,n = R i (r i,n ). (15) Finally, the matrix N ij = M i / U j describe the coupling to the Lagrange multipliers on Ω. The solution of the non-linear system Eqs. (13a) and (13b) above correspond to step (I) in the kth iteration. Then in step (II) the task is to compute [ ] Φ(γ, U(γ)) = Φ + 6 Ω Φ. (16)
7 Now, since by construction U(γ) and Λ(γ) are such that L i (U(γ), γ) 3 NT jiλ j (γ) = 0 and M i (U(γ), γ) = 0 for any γ, the derivative of those quantities with respect to γ is also zero and therefore adding any multiple, e.g. Ũ and Λ, of those to the equation above does not change the result [ ] Φ(γ, U(γ)) = Φ + = Φ + + Φ + [ŨT i [ Φ + L i + Λ T i Ũ T j Ũ T i M i [ L i ] ] N T jiλ j + Λ T i L j + Λ T j N ji ] Ui + i, Ũ T j N T Λ i ij ] [M i We see that the derivatives / and Λ i / of the implicit functions can be eliminated by choosing Ũi and Λ i as follows [10]: and Φ + (17) [ Lj ] T Ũj + N T ji Λ j = 0, i = 1, 2, 3 (18) N ij Ũ j = 0, i = 1, 2, 3. (19) In Femlab a very convenient feature is that all the matrices L i / U j, N i,j, L i /, and M i / are immediately available from the solution of the problem Eqs. (13a) and (13b) above. And the vectors Φ/ and Φ/ can also be obtained with only little extra effort. This makes the Femlab programming of the sensitivity analysis almost trivial. 5 Model system Once the topology optimization is implemented, we can move on and study how, as the Reynolds number increases from zero, the non-linear inertial term in the Navier-Stokes equation can drive transitions in the topology of networks of fluidic channels, given that the dissipation is to be minimized. To demonstrate such a transition we choose the model system presented in Fig. 1. It consists of a design region, shaped as a rectangular domain in 2D (shown in gray), to which two rectangular inlets and two outlets (shown in white) are attached symmetrically. The latter are introduced to allow for well-defined inlet and two outlet flow profiles without these imposing spurious effects in the design region. Defining the Reynolds number Re = u max lρ/η, where u max is the maximal velocity of the inlet, the two main parameters of this setup are Re and the normalized length L/l of the design region. 7
8 A C B D Figure 2: Topology optimization of the Navier-Stokes flow for two inlets and two outlets as a function of Reynolds number. Panel (A) and (C) shows the pressure-fields with flow-arrows for Re = and Re = 187.5, respectively. Panel (B) and (D) show the resulting channels (white) corresponding to (A) and (C) carved into the solid medium (black). The boundary conditions of the problem are: (1) Parabolic flow profiles at the inlets and outlets. (2) Zero pressure at the outlets. (3) No-slip boundary conditions on all other external boundaries. (4) A fixed volume-ratio β = 0.4 of the channels, see Eq. (9d). To get an idea of how the solutions of this problem will change topology as Re increases, we imagine the length L being fixed at a moderate value, say L/l = 3, and let Re increase from zero. For Re 1 the fluid is viscous, making the viscous drag from the side-walls dominate the choice of optimal solution. In order to minimize this drag, the topology optimization comes up with the configuration with the shortest possible side-walls, which is a pair of simple U-turns connecting the inlet and outlet on the same side of the design region, see panel (A) and (B) in Fig. 2. As Re increases, an increasing amount of energy is dissipated due to the change of flow direction enforced by the U-turns, and at some point this contribution dominates that of the viscous drag. At this point, this induces a transition in the topology from the two U-turns to two parallel channels connecting the inlet and outlet on opposite sides of the design region, see panel (C) and (D) in Fig. 2. For high Re, other more complicated solutions exist, such as the the top right inset in Fig. 3. 8
9 4 L/l Re Figure 3: Contour plot of the geometry parameter g = γ 1 + γ 2 as function of Re and L/l, with typical geometries inserted. Three contours at g = 0.45, 0.55, 0.65 are shown. High g (gray) corresponds to the state of straight parallel channels and low g (white) corresponds to U-turns. Throughout this work we have chosen the values α min = 0 and α max = 10 3 in Eq. (4). For the choice of q we use the following convergence-procedure adapted from Ref. [7]: initially, we set q = q 0 = , where this low value is chosen in order to find the global optimal solution. Once a steady solution is found, q is increased by a factor 1.15, and a new steady solution is found and so forth. The increase of q results in more well-defined boundaries between γ = 0 (solid) and γ = 1 (channel), and the final solution is found when q exceeds unity. 6 Results As argued in the last section and presented in Fig. 2, the parameter space (Re, L/l) will roughly be divided between the two classes of extreme solutions, the U-turns and the straight parallel channels. In addition some intermediate cases appear for high Re. The remaining part of this section will first check if the solutions are indeed optimal, and further characterize the solutions and their parameter-space. Given a parameter space (Re, L/l), it is desirable to quantify the topological changes of the solutions by the variation of a single quantity, which would then characterize the different domains in the parameter space. We have chosen to study the averages γ 1 9
10 and γ 2 of the design variable γ along the two symmetry-axes of the domain in the x 1 x 2 plane, as illustrated in Fig. 1. This choice is natural as only one of these averages are non-zero for the two extreme solutions; γ 2 0 for the U-turn, while γ 1 0 for the straight parallel channels. As these two non-zero values are different, we use their sum for characterizing the topology of the solutions: g = γ 1 + γ 2. The intermediate solutions are also resolved, since they give rise to other values of g than the two extreme cases. Fig. 3 shows a contour plot of g as a function of Re and L/l, with three contours of value g = 0.45, 0.55, 0.65 where high γ (gray) corresponds to the state of straight parallel channels and low g (white) corresponds to U-turns. The geometries of characteristic solutions of the different regions have been inserted. As expected from the discussion in Section 5, the straight parallel channels dominate for large Re and small L/l, while the U-turns dominate most of the remaining parameter-space, with exception of mixed solutions for both large Re and L/l. 7 Conclusion We have presented a new implementation of topology optimization for the full Navier- Stokes flow problem. In particular we have chosen to cast the problem into a generic divergence form suitable for the commercial finite element software package Femlab. There are several advantages of this approach. Firstly, we need not write the code for the basic meshing and matrix manipulations. Secondly, we have utilized the symbolic differentiation to ease the sensitivity analysis. Thirdly, we have used the simple and convenient way to work with any number of fields (here velocity, pressure, and design variable). Finally, our formulation is in fact general, and it can quite easily be extended to take more physical variables into account such as electrical fields and concentration gradients. As our main example we have studied how the inertial term in the steady state Navier- Stokes flow can drive transitions in the optimal topology of a microfluidic network. By treating full Navier-Stokes flows, our work is a natural extension of the initial work by Borrvall and Petersson on topology optimization of fluids in Stokes flow. We are in the process of using the topology optimization to study microfluidic systems, where the non-linearity of the flow field plays an important role even though the Reynolds numbers are small. Fluidic diodes is an example of such a system, where our method may be used to find optimal designs for real devices. For updates of this work see our web-site 8 Acknowledgement We are grateful to Ole Sigmund and Allan Gersborg-Hansen for illuminating discussions on the topology optimization method, and to Krister Svanberg for providing us with the MATLAB code for the MMA routine. This work is partly supported by The Danish Technical Research Council, Grant No
11 REFERENCES [1] Topology Optimization: Theory, Methods and Applications by M. P. Bendsøe and O. Sigmund, Springer Verlag (Berlin, 2003). [2] H. A. Eschenauer and N. Olhoff, Appl. Mech. Rev. 54, 331 (2001). [3] J. S. Jensen, J. Sound and Vibration 266, 1053 (2003). [4] J. S. Jensen and O. Sigmund, Appl. Phys. Lett. 84, 2022 (2004). [5] T. Borrvall and J. Petersson, Int. J. Num. Meth. Fluids 41, 77 (2003). [6] O. Sigmund, A. Gersborg-Hansen, and R. B. Haber, proc. Nordic Matlab Conference 2003, L. Gregersen (Ed.), Comsol A/S, Søborg, Denmark, pp [7] Topology optimization of incompressible Newtonian flows at moderate Reynolds numbers, A. Gersborg-Hansen, M.Sc.-thesis, DTU (2003). [8] Course of Theortical Physics, Vol. 6, Fluid Mechanics by L. D. Landau and E. M. Lifshitz, Butterworth and Heinemann, 2nd edition, (Oxford 2000). [9] K. Svanberg, SIAM J. Optim. 12, 555 (2002). [10] P. Michaleris, D. A. Tortorelli, and C. A. Vidal, Int. J. Numer. Methods Eng. 37, 2471 (1994). 11
DESIGN OF MICRO-FLUIDIC BIO-REACTORS USING TOPOLOGY OPTIMIZATION
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 DESIGN OF MICRO-FLUIDIC BIO-REACTORS USING TOPOLOGY OPTIMIZATION
More informationA high-level programming-language implementation of topology optimization applied to steady-state Navier Stokes flow
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2006; 65:975 1001 Published online 1 September 2005 in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/nme.1468
More informationarxiv:physics/ v2 1 Feb 2006
A high-level programming-language implementation of topology optimization applied to steady-state Navier Stokes flow Laurits Højgaard Olesen, Fridolin Okkels, and Henrik Bruus MIC Department of Micro and
More informationTopology Optimization of Micro Tesla Valve in low and moderate Reynolds number
Sen Lin (CIOMP.CAS) E-mail: SenLin41@gmail.com of Micro Tesla Valve 1 / 22 Topology of Micro Tesla Valve in low and moderate Reynolds number (Summary of my works during the past year) Sen Lin Advised by
More informationA Methodology for Microchannel Heat Sink Design Based on Topology Optimization
2 nd International Conference on Engineering Optimization September 6-9, 2010, Lisbon, Portugal A Methodology for Microchannel Heat Sink Design Based on opology Optimization Cícero R. de Lima 1, Adriano
More informationLevel Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level Set-based Topology Optimization Method for Viscous Flow Using Lattice Boltzmann Method
More informationImplementation of TOPOPT in a Microfluidic Heat Transfer Model
Microfluidics Theory and Simulation, c33443 Implementation of TOPOPT in a Microfluidic Heat Transfer Model Kristian Krüger Grøndahl, (s011426) and Thomas Glasdam Jensen, (s011427) Supervisor: Henrik Bruus
More informationDiscrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction
Discrete Projection Methods for Incompressible Fluid Flow Problems and Application to a Fluid-Structure Interaction Problem Jörg-M. Sautter Mathematisches Institut, Universität Düsseldorf, Germany, sautter@am.uni-duesseldorf.de
More informationSimulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions
Simulating Drag Crisis for a Sphere Using Skin Friction Boundary Conditions Johan Hoffman May 14, 2006 Abstract In this paper we use a General Galerkin (G2) method to simulate drag crisis for a sphere,
More informationNumerical methods for the Navier- Stokes equations
Numerical methods for the Navier- Stokes equations Hans Petter Langtangen 1,2 1 Center for Biomedical Computing, Simula Research Laboratory 2 Department of Informatics, University of Oslo Dec 6, 2012 Note:
More informationThe Influence of Channel Aspect Ratio on Performance of Optimized Thermal-Fluid Structures
Excerpt from the Proceedings of the COMSOL Conference 2010 Boston The Influence of Channel Aspect Ratio on Performance of Optimized Thermal-Fluid Structures Ercan M. Dede 1* 1 Technical Research Department,
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationSound Attenuation by Hearing Aid Earmold Tubing
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Sound Attenuation by Hearing Aid Earmold Tubing Mads J. Herring Jensen Widex A/S Ny Vestergaardsvej 19, DK-3500 Vaerloese, Denmark, mjj@widex.dk
More informationThis report shows the capabilities of Tdyn for modelling the fluid flow through porous media.
Flow through porous media 1 Flow through porous media Introduction This report shows the capabilities of Tdyn for modelling the fluid flow through porous media. Modelling approach Modeling of flows through
More informationOptimization with COMSOL Multiphysics COMSOL Tokyo Conference Walter Frei, PhD Applications Engineer
Optimization with COMSOL Multiphysics COMSOL Tokyo Conference 2014 Walter Frei, PhD Applications Engineer Product Suite COMSOL 5.0 Agenda An introduction to optimization A lot of concepts, and a little
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationNonlinear elasticity and gels
Nonlinear elasticity and gels M. Carme Calderer School of Mathematics University of Minnesota New Mexico Analysis Seminar New Mexico State University April 4-6, 2008 1 / 23 Outline Balance laws for gels
More information2 D.D. Joseph To make things simple, consider flow in two dimensions over a body obeying the equations ρ ρ v = 0;
Accepted for publication in J. Fluid Mech. 1 Viscous Potential Flow By D.D. Joseph Department of Aerospace Engineering and Mechanics, University of Minnesota, MN 55455 USA Email: joseph@aem.umn.edu (Received
More informationChapter 2: Fluid Dynamics Review
7 Chapter 2: Fluid Dynamics Review This chapter serves as a short review of basic fluid mechanics. We derive the relevant transport equations (or conservation equations), state Newton s viscosity law leading
More information7 The Navier-Stokes Equations
18.354/12.27 Spring 214 7 The Navier-Stokes Equations In the previous section, we have seen how one can deduce the general structure of hydrodynamic equations from purely macroscopic considerations and
More informationHomework 4 in 5C1212; Part A: Incompressible Navier- Stokes, Finite Volume Methods
Homework 4 in 5C11; Part A: Incompressible Navier- Stokes, Finite Volume Methods Consider the incompressible Navier Stokes in two dimensions u x + v y = 0 u t + (u ) x + (uv) y + p x = 1 Re u + f (1) v
More informationConvection-Diffusion in Microchannels
33443 Microfluidics Theory and Simulation Convection-Diffusion in Microchannels Catharina Hartmann (s97176) Luigi Sasso (s06094) Institute of Micro and Nanotechnology, DTU, 800 Kongens Lyngby, DK June,
More informationFluid flow in a channel partially filled with porous material
Fluid flow in a channel partially filled with porous material M. Hriberkk', R. Jec12& L. Skerget' 'Faculty of Mechanical Engineering, University of Maribor, Smetanova 17, SI-2000 Maribor, Slovenia 2Facultyof
More informationA Robust Preconditioned Iterative Method for the Navier-Stokes Equations with High Reynolds Numbers
Applied and Computational Mathematics 2017; 6(4): 202-207 http://www.sciencepublishinggroup.com/j/acm doi: 10.11648/j.acm.20170604.18 ISSN: 2328-5605 (Print); ISSN: 2328-5613 (Online) A Robust Preconditioned
More informationEntropy generation and transport
Chapter 7 Entropy generation and transport 7.1 Convective form of the Gibbs equation In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationFREE BOUNDARY PROBLEMS IN FLUID MECHANICS
FREE BOUNDARY PROBLEMS IN FLUID MECHANICS ANA MARIA SOANE AND ROUBEN ROSTAMIAN We consider a class of free boundary problems governed by the incompressible Navier-Stokes equations. Our objective is to
More informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationViscous Fluids. Amanda Meier. December 14th, 2011
Viscous Fluids Amanda Meier December 14th, 2011 Abstract Fluids are represented by continuous media described by mass density, velocity and pressure. An Eulerian description of uids focuses on the transport
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 6A: Internal Incompressible Viscous Flow 4/16/2018 C6A: Internal Incompressible Viscous Flow 1 6.1 Introduction For the present chapter we will limit our study to incompressible
More informationNumerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall
Numerical Analysis of MHD Flow of Fluid with One Porous Bounding Wall Ramesh Yadav 1 & Vivek Joseph 2 1Assistant Professor, Department of Mathematics BBDNITM Lucknow U P 2Professor, Department of Mathematics
More informationPressure Losses for Fluid Flow Through Abrupt Area. Contraction in Compact Heat Exchangers
Pressure Losses for Fluid Flow Through Abrupt Area Contraction in Compact Heat Exchangers Undergraduate Research Spring 004 By Bryan J. Johnson Under Direction of Rehnberg Professor of Ch.E. Bruce A. Finlayson
More informationCranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom. Cranfield University, Cranfield, Bedfordshire, MK43 0AL, United Kingdom
MultiScience - XXX. microcad International Multidisciplinary Scientific Conference University of Miskolc, Hungary, 21-22 April 2016, ISBN 978-963-358-113-1 NUMERICAL INVESTIGATION OF AN INCOMPRESSIBLE
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationWall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes
Excerpt from the Proceedings of the COMSOL Conference 9 Boston Wall Effects in Convective Heat Transfer from a Sphere to Power Law Fluids in Tubes Daoyun Song *1, Rakesh K. Gupta 1 and Rajendra P. Chhabra
More informationShell Balances in Fluid Mechanics
Shell Balances in Fluid Mechanics R. Shankar Subramanian Department of Chemical and Biomolecular Engineering Clarkson University When fluid flow occurs in a single direction everywhere in a system, shell
More informationAnalysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications
Analysis of Mixing Chambers for the Processing of Two-Component Adhesives for Transport Applications P. Steinert* 1, I. Schaarschmidt 1, R. Paul 1, M. Zinecker 1, M. Hackert-Oschätzchen 1, Th. Muschalek
More informationFundamentals of Fluid Dynamics: Elementary Viscous Flow
Fundamentals of Fluid Dynamics: Elementary Viscous Flow Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research
More informationNUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL
Proceedings of ALGORITMY 212 pp. 29 218 NUMERICAL SIMULATION OF INTERACTION BETWEEN INCOMPRESSIBLE FLOW AND AN ELASTIC WALL MARTIN HADRAVA, MILOSLAV FEISTAUER, AND PETR SVÁČEK Abstract. The present paper
More informationParabolic Flow in Parallel Plate Channel ME 412 Project 4
Parabolic Flow in Parallel Plate Channel ME 412 Project 4 Jingwei Zhu April 12, 2014 Instructor: Surya Pratap Vanka 1 Project Description The objective of this project is to develop and apply a computer
More informationThe dependence of the cross-sectional shape on the hydraulic resistance of microchannels
3-weeks course report, s0973 The dependence of the cross-sectional shape on the hydraulic resistance of microchannels Hatim Azzouz a Supervisor: Niels Asger Mortensen and Henrik Bruus MIC Department of
More information6.2 Governing Equations for Natural Convection
6. Governing Equations for Natural Convection 6..1 Generalized Governing Equations The governing equations for natural convection are special cases of the generalized governing equations that were discussed
More informationCHAPTER 8 ENTROPY GENERATION AND TRANSPORT
CHAPTER 8 ENTROPY GENERATION AND TRANSPORT 8.1 CONVECTIVE FORM OF THE GIBBS EQUATION In this chapter we will address two questions. 1) How is Gibbs equation related to the energy conservation equation?
More informationA Study on Numerical Solution to the Incompressible Navier-Stokes Equation
A Study on Numerical Solution to the Incompressible Navier-Stokes Equation Zipeng Zhao May 2014 1 Introduction 1.1 Motivation One of the most important applications of finite differences lies in the field
More informationApplication of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder with Regressing Walls
Mechanics and Mechanical Engineering Vol. 21, No. 2 (2017) 379 387 c Lodz University of Technology Application of Reconstruction of Variational Iteration Method on the Laminar Flow in a Porous Cylinder
More informationGetting started: CFD notation
PDE of p-th order Getting started: CFD notation f ( u,x, t, u x 1,..., u x n, u, 2 u x 1 x 2,..., p u p ) = 0 scalar unknowns u = u(x, t), x R n, t R, n = 1,2,3 vector unknowns v = v(x, t), v R m, m =
More informationLaminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid
Laminar Forced Convection Heat Transfer from Two Heated Square Cylinders in a Bingham Plastic Fluid E. Tejaswini 1*, B. Sreenivasulu 2, B. Srinivas 3 1,2,3 Gayatri Vidya Parishad College of Engineering
More informationAA210A Fundamentals of Compressible Flow. Chapter 1 - Introduction to fluid flow
AA210A Fundamentals of Compressible Flow Chapter 1 - Introduction to fluid flow 1 1.2 Conservation of mass Mass flux in the x-direction [ ρu ] = M L 3 L T = M L 2 T Momentum per unit volume Mass per unit
More informationTurbulent drag reduction by streamwise traveling waves
51st IEEE Conference on Decision and Control December 10-13, 2012. Maui, Hawaii, USA Turbulent drag reduction by streamwise traveling waves Armin Zare, Binh K. Lieu, and Mihailo R. Jovanović Abstract For
More informationA Permeability Function for Brinkman s Equation
Proceedings of the th WSEAS Int. Conf. on MATHEMATICAL METHOS, COMPUTATIONAL TECHNIQUES AN INTELLIGENT SYSTEMS A Permeability Function for Brinman s Equation M.H. HAMAN, M.T. KAMEL, and H.I. SIYYAM 3,
More informationFlow Focusing Droplet Generation Using Linear Vibration
Flow Focusing Droplet Generation Using Linear Vibration A. Salari, C. Dalton Department of Electrical & Computer Engineering, University of Calgary, Calgary, AB, Canada Abstract: Flow focusing microchannels
More informationMicro Cooling of SQUID Sensor
Excerpt from the Proceedings of the COMSOL Conference 2008 Hannover Micro Cooling of SQUID Sensor B.Ottosson *,1, Y. Jouahri 2, C. Rusu 1 and P. Enoksson 3 1 Imego AB, SE-400 14 Gothenburg, Sweden, 2 Mechanical
More informationUniversität des Saarlandes. Fachrichtung 6.1 Mathematik
Universität des Saarlandes U N I V E R S I T A S S A R A V I E N I S S Fachrichtung 6. Mathematik Preprint Nr. 228 A Variational Multiscale Method for Turbulent Flow Simulation with Adaptive Large Scale
More informationChapter 9: Differential Analysis
9-1 Introduction 9-2 Conservation of Mass 9-3 The Stream Function 9-4 Conservation of Linear Momentum 9-5 Navier Stokes Equation 9-6 Differential Analysis Problems Recall 9-1 Introduction (1) Chap 5: Control
More information2. FLUID-FLOW EQUATIONS SPRING 2019
2. FLUID-FLOW EQUATIONS SPRING 2019 2.1 Introduction 2.2 Conservative differential equations 2.3 Non-conservative differential equations 2.4 Non-dimensionalisation Summary Examples 2.1 Introduction Fluid
More informationStudy of Forced and Free convection in Lid driven cavity problem
MIT Study of Forced and Free convection in Lid driven cavity problem 18.086 Project report Divya Panchanathan 5-11-2014 Aim To solve the Navier-stokes momentum equations for a lid driven cavity problem
More informationCFD ANALYSIS OF TRIANGULAR ABSORBER TUBE OF A SOLAR FLAT PLATE COLLECTOR
Int. J. Mech. Eng. & Rob. Res. 2013 Basavanna S and K S Shashishekar, 2013 Research Paper ISSN 2278 0149 www.imerr.com Vol. 2, No. 1, January 2013 2013 IJMERR. All Rights Reserved CFD ANALYSIS OF TRIANGULAR
More informationOptimization in Continuum Flow Problems
Linköping Studies in Science and Technology. Dissertations. No. 1196 Optimization in Continuum Flow Problems Niclas Wiker Division of Mechanics Department of Management and Engineering Linköping University,
More informationThe effective slip length and vortex formation in laminar flow over a rough surface
The effective slip length and vortex formation in laminar flow over a rough surface Anoosheh Niavarani and Nikolai V. Priezjev Movies and preprints @ http://www.egr.msu.edu/~niavaran A. Niavarani and N.V.
More informationThe effect of geometric parameters on the head loss factor in headers
Fluid Structure Interaction V 355 The effect of geometric parameters on the head loss factor in headers A. Mansourpour & S. Shayamehr Mechanical Engineering Department, Azad University of Karaj, Iran Abstract
More informationMicrofluidics 1 Basics, Laminar flow, shear and flow profiles
MT-0.6081 Microfluidics and BioMEMS Microfluidics 1 Basics, Laminar flow, shear and flow profiles 11.1.2017 Ville Jokinen Outline of the next 3 weeks: Today: Microfluidics 1: Laminar flow, flow profiles,
More informationChapter 9: Differential Analysis of Fluid Flow
of Fluid Flow Objectives 1. Understand how the differential equations of mass and momentum conservation are derived. 2. Calculate the stream function and pressure field, and plot streamlines for a known
More informationEntropic Evaluation of Dean Flow Micromixers
COMSOL Conference, Boston, 2013 Brian Vyhnalek, Petru S. Fodor and Miron Kaufman Physics Department Cleveland State University Entropic Evaluation of Dean Flow Micromixers ABSTRACT We study the use of
More informationStudies on flow through and around a porous permeable sphere: II. Heat Transfer
Studies on flow through and around a porous permeable sphere: II. Heat Transfer A. K. Jain and S. Basu 1 Department of Chemical Engineering Indian Institute of Technology Delhi New Delhi 110016, India
More informationEffects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
Transport in Porous Media (2006) 64: 1 14 Springer 2006 DOI 10.1007/s11242-005-1126-6 Effects of Viscous Dissipation on Unsteady Free Convection in a Fluid past a Vertical Plate Immersed in a Porous Medium
More informationImproved Two-Phase Projection Topology Optimization
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Improved Two-Phase Projection Topology Optimization Josephine V. Carstensen and James K. Guest
More informationCFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel
CFD Analysis of Forced Convection Flow and Heat Transfer in Semi-Circular Cross-Sectioned Micro-Channel *1 Hüseyin Kaya, 2 Kamil Arslan 1 Bartın University, Mechanical Engineering Department, Bartın, Turkey
More informationTopology optimisation of passive coolers for light-emitting diode lamps
Downloaded from orbit.dtu.dk on: Aug 15, 2018 Topology optimisation of passive coolers for light-emitting diode lamps Alexandersen, Joe Published in: Proceedings of WCSMO-11 Publication date: 2015 Document
More informationASTR 320: Solutions to Problem Set 2
ASTR 320: Solutions to Problem Set 2 Problem 1: Streamlines A streamline is defined as a curve that is instantaneously tangent to the velocity vector of a flow. Streamlines show the direction a massless
More informationA microfluidic-based hydrodynamic trap: Design and implementation
SUPPLEMENTARY MATERIAL A microfluidic-based hydrodynamic trap: Design and implementation Melikhan Tanyeri, a Mikhil Ranka, a Natawan Sittipolkul a and Charles M. Schroeder* a,b a Department of Chemical
More informationCESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP
8 th GRACM International Congress on Computational Mechanics Volos, 2 July 5 July 205 CESSATION OF VISCOPLASTIC POISEUILLE FLOW IN A RECTANGULAR DUCT WITH WALL SLIP Yiolanda Damianou, George Kaoullas,
More informationSYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL
ISSN 1726-4529 Int j simul model 7 (2008) 4, 186-197 Original scientific paper SYMMETRY BREAKING PHENOMENA OF PURELY VISCOUS SHEAR-THINNING FLUID FLOW IN A LOCALLY CONSTRICTED CHANNEL Ternik, P. University
More informationEffect of Couple Stresses on the MHD of a Non-Newtonian Unsteady Flow between Two Parallel Porous Plates
Effect of Couple Stresses on the MHD of a Non-Newtonian Unsteady Flow between Two Parallel Porous Plates N. T. M. Eldabe A. A. Hassan and Mona A. A. Mohamed Department of Mathematics Faculty of Education
More informationMHD flow and heat transfer due to a linearly stretching sheet. with induced magnetic field: Exact solution. Tarek M. A.
MHD flow and heat transfer due to a linearly stretching sheet with induced magnetic field: Exact solution Tarek M. A. El-Mistikawy Dept. Eng. Math. & Phys., Faculty of Engineering, Cairo University, Giza
More informationENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD
HEFAT 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 6 8 July Malta ENTROPY GENERATION IN HEAT AND MASS TRANSFER IN POROUS CAVITY SUBJECTED TO A MAGNETIC FIELD Nawaf
More informationDifferential relations for fluid flow
Differential relations for fluid flow In this approach, we apply basic conservation laws to an infinitesimally small control volume. The differential approach provides point by point details of a flow
More informationp + µ 2 u =0= σ (1) and
Master SdI mention M2FA, Fluid Mechanics program Hydrodynamics. P.-Y. Lagrée and S. Zaleski. Test December 4, 2009 All documentation is authorized except for reference [1]. Internet access is not allowed.
More informationBoundary-Layer Theory
Hermann Schlichting Klaus Gersten Boundary-Layer Theory With contributions from Egon Krause and Herbert Oertel Jr. Translated by Katherine Mayes 8th Revised and Enlarged Edition With 287 Figures and 22
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationHelical Coil Flow: a Case Study
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Helical Coil Flow: a Case Study Marco Cozzini Renewable Energies and Environmental Technologies (REET) Research Unit, Fondazione Bruno Kessler
More informationProject 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions
Project 4: Navier-Stokes Solution to Driven Cavity and Channel Flow Conditions R. S. Sellers MAE 5440, Computational Fluid Dynamics Utah State University, Department of Mechanical and Aerospace Engineering
More informationBasic concepts in viscous flow
Élisabeth Guazzelli and Jeffrey F. Morris with illustrations by Sylvie Pic Adapted from Chapter 1 of Cambridge Texts in Applied Mathematics 1 The fluid dynamic equations Navier-Stokes equations Dimensionless
More informationInterpreting Differential Equations of Transport Phenomena
Interpreting Differential Equations of Transport Phenomena There are a number of techniques generally useful in interpreting and simplifying the mathematical description of physical problems. Here we introduce
More informationLattice Boltzmann Method for Fluid Simulations
1 / 16 Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas Simon Fraser University April 7, 2011 2 / 16 Ludwig Boltzmann and His Kinetic Theory of Gases The Boltzmann Transport
More informationTHE EFFECT OF SLIP CONDITION ON UNSTEADY MHD OSCILLATORY FLOW OF A VISCOUS FLUID IN A PLANER CHANNEL
THE EFFECT OF SLIP CONDITION ON UNSTEADY MHD OSCILLATORY FLOW OF A VISCOUS FLUID IN A PLANER CHANNEL A. MEHMOOD, A. ALI Department of Mathematics Quaid-i-Azam University 4530, Islamabad 44000 Pakistan
More informationMHD Free convection flow of couple stress fluid in a vertical porous layer
Available online at www.pelagiaresearchlibrary.com Advances in Applied Science Research,, (6:5- ISSN: 976-86 CODEN (USA: AASRFC MHD Free convection flow of couple stress fluid in a vertical porous layer
More informationLattice Boltzmann Method for Fluid Simulations
Lattice Boltzmann Method for Fluid Simulations Yuanxun Bill Bao & Justin Meskas April 14, 2011 1 Introduction In the last two decades, the Lattice Boltzmann method (LBM) has emerged as a promising tool
More informationPublication 97/2. An Introduction to Turbulence Models. Lars Davidson, lada
ublication 97/ An ntroduction to Turbulence Models Lars Davidson http://www.tfd.chalmers.se/ lada Department of Thermo and Fluid Dynamics CHALMERS UNVERSTY OF TECHNOLOGY Göteborg Sweden November 3 Nomenclature
More informationChapter 1. Continuum mechanics review. 1.1 Definitions and nomenclature
Chapter 1 Continuum mechanics review We will assume some familiarity with continuum mechanics as discussed in the context of an introductory geodynamics course; a good reference for such problems is Turcotte
More informationDouble penalisation method for porous-fluid problems and numerical simulation of passive and active control strategies for bluff-body flows
Double penalisation method for porous-fluid problems and numerical simulation of passive and active control strategies for bluff-body flows Iraj Mortazavi Institut Polytechnique de Bordeaux INRIA Bordeaux
More informationSTEADY AND UNSTEADY 2D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS FLOW. Radka Keslerová, Karel Kozel
Conference Applications of Mathematics 1 in honor of the th birthday of Michal Křížek. Institute of Mathematics AS CR, Prague 1 STEADY AND UNSTEADY D NUMERICAL SOLUTION OF GENERALIZED NEWTONIAN FLUIDS
More informationExternal and Internal Incompressible Viscous Flows Computation using Taylor Series Expansion and Least Square based Lattice Boltzmann Method
Available online at http://ijim.srbiau.ac.ir/ Int. J. Industrial Mathematics (ISSN 2008-5621) Vol. 10, No. 2, 2018 Article ID IJIM-00726, 8 pages Research Article External and Internal Incompressible Viscous
More informationComputation of turbulent Prandtl number for mixed convection around a heated cylinder
Center for Turbulence Research Annual Research Briefs 2010 295 Computation of turbulent Prandtl number for mixed convection around a heated cylinder By S. Kang AND G. Iaccarino 1. Motivation and objectives
More informationThe Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011
The Shape of a Rain Drop as determined from the Navier-Stokes equation John Caleb Speirs Classical Mechanics PHGN 505 December 12th, 2011 Derivation of Navier-Stokes Equation 1 The total stress tensor
More informationForchheimer law derived by homogenization of gas flow in turbomachines
Journal of Computational and Applied Mathematics 215 (2008) 467 476 www.elsevier.com/locate/cam Forchheimer law derived by homogenization of gas flow in turbomachines Gottfried Laschet ACCESS e.v., IntzestraYe
More information