Implementation of TOPOPT in a Microfluidic Heat Transfer Model
|
|
- Marian Neal
- 5 years ago
- Views:
Transcription
1 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 MIC Department of Micro and Nanotechnology Technical University of Denmark 22 June 2005
2 ii
3 Abstract The material distribution method in topology optimization can be applied in order to solve numerous physical problems numerically. we have used a modified version of a recently developed high-level programming implementation of topology optimization in order to model and simulate heat transfer in cross-flow in two dimensions. Inspired by the article Constructal multi-scale cylinders in cross-flow by T. Bello-Ochende and A.Bejan [1], where the finite elements method (FEM) is used to calculate the effect of cooling cylinders in a liquid flow. We have performed a vast amount of simulations using FEMLAB and MAT- LAB in order to determine the optimum placing of heat dissipating material in the cooling liquid. Our work has shown that it is possible to implement the model and simulate it in MATLAB with FEMLAB but fast processors and a huge (minimum 1 GB) RAM facilities are needed to carry out the iterations in a reasonable time. We have also stated that the optimum shape of the material to cool is far from cylindrical but has large surface areas which is expected in order to get the biggest contact area where the heat transport takes place. iii
4 iv ABSTRACT
5 Contents 1 Introduction The system Inspiration Distributing material in channel Theoretical Aspects Fluid Dynamics The continuum hypothesis The Navier-Stokes equation HeatEquation The γ-field Dimensionless equation Implementation of TOPOPT Setting up Femlab Adapting the code Heat dissipation The Objective function Global effects of local changes Permeable material Adjusting the material parameters The conservative heat equation Usolve Results Conservative Vs. Non-conservative Making Q independent of the γ-area Lowering Darcy s number Varying C p Varying k Varying ρ Varying the initial conditions Further work 29 v
6 vi CONTENTS 6 Conclusion 31
7 Chapter 1 Introduction In development of Lab On a Chip (LOC) systems cooling can be of high importance. Detecting bacteria on chickens requires amplification of DNA which is done with PCR (Polymerase Chain Reaction) on chip, basically copying of DNA. Cooling is needed in order to control the fluid temperature to a high precision. The denaturation of the DNA (splitting the double stranded DNA) involve one temperature and the enzymatic reactions taking care of copying DNA involves another. This has lead to the idea of optimizing the cooling process, basically the heat transfer between a cooling liquid and a heat dissipating material (e.g.iron or cobber). Looking into the problem it is clear that the shape and thermodynamically properties plays a superior role. In this report it is discovered how the transfer of heat into a cooling liquid can be optimized by choosing a best-shape of the heaters. 1.1 The system The idea is to setup a closed channel with inlet and outlet and stalactites-like tubes going from top to bottom. Inducing a flow of a cooling liquid will transfer the heat from the stalactites, out of the system by convection. The questions are how should the shape of the stalactites be and what should be their individual placement Inspiration A group of scientists have worked with heat transfer in cross-flow as well [1]. In 2 dimensions with the channel seen from above, they have placed material in center of the fluid stream with a circular shape, illustrating a tube. The characteristic of the heat transfer from tube to fluid was numerically solved and the optimal distance from tube to channel wall was found. Further on smaller tubes were placed in between the channel wall and the initial tube in order to find an optimal number of tubes for transferring the heat to the liquid, taking into account that the channel must not be blocked by tubes. Again the optimal placement was found by numerical calculations. In the article of Bejan A. the shape of the tubes were cylindrical. This is exactly the point where it is believed that further optimization can be done. 1
8 2 CHAPTER 1. INTRODUCTION Figure 1.1: Picture of dripstones hanging from ceiling to floor in a cave. This is how the channel is supposed to look a like. When a flow of liquid is induced a cross flow is established and a temperature difference between solid material and fluid makes heat transfer possible
9 1.1. THE SYSTEM Distributing material in channel A special type of topology optimization, abbreviated the material distribution method (MDM), has been utilized in diverse areas such as structural mechanics, acoustics, optics and fluidics for optimizing performance. Basically the method is used for optimizing a Figure 1.2: Schematic of the setup used in the simulations. Due to horizontal symmetry only half of the model is evaluated, easing the calculation process. A pressure difference between inlet and outlet, p, drives the cooling liquid through the main design region. Heat dissipating material is modeled by both spatially varying the permeability in the design region and varying the material properties towards metallic values. Furthermore the objective function to be minimized is activated in this central region. given objective function, Φ, by mapping (i.e. distributing material) a designated design region. Recently the method has been applied to a full Navier-Stokes flow through the implementation in a high-level programming language based on FEMLAB and MATLAB [2]. The development of new and innovative fluidic and micro fluidic devices will benefit from this progress. In the present report we have applied the MDM in order to study the model of a 2D heat exchanger for use in LOC. The basic setup is simple. Fluid with a fixed inlet temperature is pumped through the design region by a pressure difference p. By iteratively varying the material distribution in the design region, the placement of heat dissipating matter is varied in order to minimize the average material temperature. Albeit the setup is simple, refer to Fig. 1.2, the task of finding the optimal material distribution is a complex Navier- Stokes conduction convection problem linked to the minimization of the temperature, which can only be solved numerically. The main objective has been to implement the TOPOPT-algorithm into a heat transfer problem. Explained in Chapter 3.
10 4 CHAPTER 1. INTRODUCTION
11 Chapter 2 Theoretical Aspects 2.1 Fluid Dynamics In this section the theory concerning the motions of fluids will be discussed. First with a brief description of the continuum hypothesis, then the different parts of the Navier-Stokes equation are explained The continuum hypothesis This work is based on the continuum hypothesis, in which a fluid is not treated as individual molecules. Instead large numbers of molecules are blocked together to form fluid particles. These fluid particles are usually cubic and have a side length, λ of about 10 nm, such a fluid particle contains about molecules. The reason why the fluid is not treated as individual molecules is the fact that these have very large fluctuation. Repeated measurements would therefore yield different result and hence not be reproducible. Statistics then gives a lower boundary to the number of molecules in a fluid particle, the upper boundary is needed to keep the fluid particle small enough to avoid large variations in external forces. When using the continuum hypothesis it is important to work on length scales much larger than λ 10 nm otherwise theory breaks down. However, this is not a big concern as most Lab-on-a-Chip systems are made with length scales in the order of several µm The Navier-Stokes equation The Navier-Stokes equation can be derived from Newton s second law m a = F. First the acceleration, a is substituted with time-derivative of the velocity, d t v. Then the force is expressed as the sum of several different external forses, F = j F j. Now Newton s second equation look as follows m d t v = F j. (2.1) j However, since this equation is going to be used on the fluid particles and not on the individual molecules it is necessary to divide with the volume of such a fluid particle. This 5
12 6 CHAPTER 2. THEORETICAL ASPECTS leaves the density, ρ instead of the mass, m, the force densities f j instead of the forces F j And the material time-derivative, D t instead of d t. The equation now takes the following form; ρ D t v = f j. (2.2) j Since Newton s second law follows the Eulerian description, where the Navier-Stokes equation follows the Lagrangian description, it is necessary to use the material time-derivative, D t = t + (v ). Now the Navier-Stokes equation can be expressed as; ( ) ρ t v + (v )v = f j, (2.3) j All that now remains is to specify the external forces. Normally four forces must be taken into account, these forces are caused by; gravity, pressure-gradient, viscosity and electricity. All forces must be used as force densities. The gravitational force density has the form f grav = ρg. The force density caused by the pressure gradient is f pres = p. When dealing with incompressible fluids, as in this entire work, the viscous force density is of the form f visc = η 2 v. The fourth and final force density is that of electricity, this force density has the form f el = ρ el E where ρ el is the charge density. If all of these force densities are substituted into Eq. (2.3) the Navier-Stokes equation looks like this: ( ) ρ t v + (v )v = p + η 2 v + ρ g + ρ el E, (2.4) It is often possible to neglect the gravitational force, as it is balanced out by an opposite normal force.
13 2.2. HEATEQUATION HeatEquation The total time dependent heat equation including conduction and convection reads δ ts ρc p T t + ( k T ) = Q ρc pu T (2.5) For the steady state problem the equation is time independent and the first term disappears The properties of the elements are ( k T ) = Q ρc p u T (2.6) C p k Q u ρ heat capacity thermal conductivity tensor heat source velocity field vector the density The equation Eq. (2.6) written above is in its non-conservative form. When the nonconservative mode is chosen it is assumed that the continuity equation for the velocity field, u = 0, is fulfilled. From lectures in Theoretical microfluidics we know that this is true for incompressible fluids. The non-conservative mode also ensures that no unphysical source term arises from a flow field where the incompressibility is not fulfilled. The conservative form of the time independent heat equation allows us to alter the physical parameters; density, heat capacity, and conductivity since they are placed inside the integral as seen ( k T +ρc p ut ) = Q (2.7) The practical consequences are further explained in Section
14 8 CHAPTER 2. THEORETICAL ASPECTS 2.3 The γ-field The γ-field describes the ratio between the amount of cross-flowing liquid and heat dissipating material in any area. If γ = 1 no material is present in this area of the channel and the liquid can flow freely. If γ = 0 then this part of the channel is blocked by the material. To simulate that the material is obstructing the flow, the Darcy s force has been implemented. The Darcy s force looks like this Where F x = α v (2.8) α 1 γ, γ = (2.9) Da and Da is the Darcy s number and has a value about This results in a strong force working in the opposite direction than that of v.
15 2.4. DIMENSIONLESS EQUATION Dimensionless equation In the article Constructal multi-scale cylinders in cross-flow of T. Bello-Ochende, Bejan [1] a dimensionless heat transfer density rate is derived, in order to be able to optimize the spacing between the cylinders placed in the fluid stream (see page 1376). To set up the scenario for the optimization it has therefore been necessary to rearrange the conservation equation for mass, momentum and, energy, into dimensionless equations. To get a feeling for the method we will rearrange the momentum equation into a dimensionless momentum equation. The equation for conservation of momentum reads u v x + v v y = 1 P ρ y + ν 2 v (2.10) The dimensionless variables needed for the rearrangement are ( x, ỹ) = (x, y) D 0, (ũ, ṽ) = (u, v) P D 0 /µ, P P = P (2.11) We also need the Prandtl number defined P r = ν α (2.12) and the dimensionless pressure drop number, named the Bejan number Be = P D 0 αµ (2.13) Starting with the right-hand side of the momentum conservation equation, ρ is from its SI-units seen to be equal ν /µ. It gives u v x + v v y = ν P µ y + ν 2 v (2.14) substituting the variables leads to u v x + v v y = ν µ ( P P ) ( (ỹ D 0 ) + ν 2 ṽ P D ) 0 µ Introducing the Prandtl number and the dimensionless Bejan number gives u v x + v v ( Be α 2 y = P r ( P P ) P D0 2 (ỹ D 0 ) + αµ ) 2 ṽ P D 0 u v x + v v ( Be α 2 y = P r ( P P ) Be ) α2 P D ṽ (ỹ D 0 ) D 0 Putting the Bejan number and D 0 outside the parenthesis u v x + v v y r Be ( α 2 = P P ) D 0 D0 2 ỹ + α2 2 ṽ (2.15) (2.16) (2.17) (2.18)
16 10 CHAPTER 2. THEORETICAL ASPECTS Extracting D 2 0 and α2 u v x + v v y P r Be ( α2 P ) = D0 3 ỹ + D2 0 2 ṽ (2.19) Keeping in mind that 2 ṽ is equal to ṽ reads u v x + v v y x, ṽ y, we get D 2 0 with Eq. (2.11) the right side P r Be ( α2 P ) = D0 3 ỹ + 2 ṽ Substituting the variables with dimensionless variables on the left hand side ) ) ( P D0 ) ṽ( P D0 ( µ P D0 ) ṽ( P D0 µ P r Be ( α2 P ) ũ + ṽ = µ x D 0 µ ỹd 0 D0 3 ỹ + 2 ṽ (2.20) (2.21) Collecting the tilde variables P 2 D0 2 1 ( µ 2 ũ ṽ D 0 x + ṽ ṽ ) ỹ () Collecting terms outside the on the left side = P r Be ( α2 P ) D0 3 ỹ + 2 ṽ P 2 D 0 D 3 ( 0 µ 2 P r Be α 2 ũ ṽ x + ṽ ṽ ) ( P ) = ỹ ỹ + 2 ṽ (2.22) (2.23) and the total dimen- With the definition in Eq. (2.13) the prefactor can be written as Be P r sionless equation for conservation of momentum is Be ( ũ ṽ P r x + ṽ ṽ ) ỹ ( P ) = ỹ + 2 ṽ (2.24) This agrees with the dimensionless equation for conservation of momentum seen in the article of T. Bello-Ochende and A. Bejan
17 Chapter 3 Implementation of TOPOPT Our main objective have been to modify the source code of the high-level programming implementation of topology optimization in order to simulate heat transfer in a well defined design region. 3.1 Setting up Femlab The MDM is used to model the heat dissipating material density. The problem has a complex multi physical nature composed of a Poiseuille flow in which both convection and conduction of the heat is present. Furthermore the temperature is evaluated in the main design region, in which the heat dissipating material density is iteratively altered by the algorithm in order to achieve the minimal integral value of the total temperature. In FEMLAB we have modeled the heat transfer by choosing appropriate subdomain and boundary settings. For the incompressible Navier-Stokes part, which is activated in all regions, we have defined the pressure at the inlet and outlet and imposed the no slip boundary conditions on all outer walls, the slip/symmetry condition was chosen for the all symmetry boundaries. In the case of the convection and conduction of the heat we have defined T 0 = 0 at the inlet and furthermore imposed a convective flux boundary condition at the outlet. The heat dissipation is solely activated in the main design region and as a function of γ in order to simulate that the material are not allowed to move in or out of the system. In order to include the design variable γ, describing the spatial permeability (i.e. the material density), we finally define a PDE, General Form with linear lagrange elements in the design region. By modifying the source code as describe in notes handed out by Fridolin Okkels and adapting the functionality (as described below) to heat transfer, rather than metabolism, the program was able to calculate heat transfer. In Fig. 3.1 we have included an example of the graphs produced by running the simulation. By running numerous simulations we can thus investigate the impact of both geometrical and physical parameters on the outlet concentration. 11
18 12 CHAPTER 3. IMPLEMENTATION OF TOPOPT Figure 3.1: Example of the output from running a simulation. The output displays the material distribution, the flow velocity, the temperature. The mesh size is in the design region and on the boundaries between the inlet/design region and design region/outlet. 3.2 Adapting the code Heat dissipation The first problem that was encountered during the implementation was how to simulate heat dissipation. It seemed that Femlab offered two different possibilities, either one should set the surface temperature, T of the material to a fixed value and then let Femlab work out how this would affect the surrounding liquid, or one could set the amount of energy dissipated,q from the material into the liquid. The problem with fixing the surface temperature, is that the system no longer works as a cooling devise, instead it works as a water heater. In the case of fixing Q the same amount of energy is released into the liquid no matter how fast it is moving past the material or how hot it already is. It was decided to use the fixed Q value, but since Q is a energy flux it is dependent on the area of the material it is dissipating from. This means that as the MDM algorithm places more or less materiel in the design region the amount of dissipated energy will increase or decrease, respectively. Unfortunately this leads to the obvious conclusion that if no material is placed at all, then no energy will be dissipated and hence the overall temperature will be the lowest possible. thereby, tricking the algorithm to believe that the ultimate cooling structure has been achieved. This was fixed by adding two modifications; the first was to set a minimum value of material volume, the second was to divide the Q by the area of the materia, thereby fixing the amount of energy dissipated no matter the amount of
19 3.2. ADAPTING THE CODE 13 material The Objective function The objective function is the function the MDM algorithm tries to optimize by distributing material. In the code this function is called Phi. In our case of a cooling system, the obvious choice would be to maximize the amount of energy carried out of the system by the liquid. However, this is not possible since it was chosen to use fixed Q instead of fixed T. When specifying Q and working with steady state systems, it is given that the same energy which comes into the system must come out. Hence the amount of energy carried away by the liquid is independent of material distribution. To avoid this problem the objective function was chosen to depend on the temperature in the design region. This also caused some troubles since it would be better to place the material as late in the design region as possible. This was not the desired solution, so instead of looking at the temperature in the entire design region, only the temperature of the material should be evaluated. This was done by changing the objective function to: Where γ = 1 is pure liquid and γ = 0 is solid material Global effects of local changes Φ = T (1 γ) (3.1) During the simulations it was discovered that a important feature had been over looked. Namely the fact that when dividing the amount of dissipated heat, Q with the area of the γ-field, then local changes causes Global effects. In other words, when the γ-field is changed in one place, the temperature is changed everywhere. This might make it very difficult for MDM-algorithm to optimize the structure. However, in the relative simple cases studied here the MDM-algorithm seems to handle the job just fine. If this code is further expanded, one should keep an eye on this effect Permeable material Many of the initial simulations were made with a Darcy s number of Da = This made the γ-field settle at some value between 0 and 1. Thereby meaning that the optimal structure would be a permeable heat dissipating material. This problem could be avoided in two ways; one idea was simply to decrease the darcy s number to Da = making the permeability of the material much less. Another idea would simply be to refine the mesh size, but due to the lack of computational powers and limited time, the mesh was not further decreased. With these changes it becomes much easier for the MDM-algorithm to decide on whether γ should be one or zero Adjusting the material parameters In order to make the simulation as physical as possible it is necessary to adjust the material parameters as γ is varied between zero and one. When γ = 1 the material is the liquid
20 14 CHAPTER 3. IMPLEMENTATION OF TOPOPT which in this case is water. All material parameters such as the density ρ, the specific heat capacity C p and the thermal conductivity k should then take values corresponding to those of water. On the other hand when γ = 0 the parameters should match those of the material (copper in our case). Hence all material parameters should be functions of γ. Below is shown how this has been implemented in the case of the thermal conductivity, the procedure is exactly the same for the other parameters. k water = 0.6 k copper = 385 k material = k vand γ+k copper (1 γ), γ = (3.2) All these parameters was never implemented at the same time, due to several crashes in the algorithm The conservative heat equation To prevent the algorithm from crashing due to changes in the material parameters, it was necessary to change the non-conservative heat equation to a conservative one, see Section 2.2. Unfortunately the conservative equation give much rougher result than the nonconservative, this means that in order to obtain the same resolution as earlier the mesh must be further refined. This was however not possible due to the lack of computational powers and limited time. The switch from the default non-conservative equation to the conservative is made by the follwing line. fem.appl{- x -}.prop.equform = cons ; where - x - is the application number corresponding to convection and conduction. With this line implemented it is possible change the material parameters as a function of γ Usolve FEMLAB has a built in scaling function resulting in division by a high number ( ) when running our main-file. This forces the pc to calculate with numbers close to its uncertainty value, making it crash. To avoid the problem we have disabled the scaling by the following script, when calling the femlin function fem.sol=femnlin(fem, init,fem.sol, solcomp,{ u, v, p },... outcomp,{ u, v, p, gamma },... nonlin, on, uscale, none );
21 Chapter 4 Results In this chapter the process is described and the results are discussed. The first attempt to run the MatLab-code modeling our heat exchanger failed. Every parameter of the material had to be disabled or smeared out giving the heat dissipating material the same physical properties as the cooling fluid (water), in order to ease the calculations and making it possible for the PC to find a solution. Restarting with the simplest setup, and adding one parameter at a time gave us the ability to learn how Femlab responded to different codes. Parameters: for both the fluid and the material distributed the parameters are; Da = k = 0.6 C p = 4186 ρ = Conservative Vs. Non-conservative The first test was conducted to check that two simulations, one using the conservative heat equation and one using the non-conservative, would yield the same result. In all three simulations the exact same parameters were used and the simulations lasted for 100 iterations. The first test was made with a grid size of and using the conservative equation. The result is shown in Fig The second simulation was made with the non-conservative heat equation, The result is shown in Fig When comparing Fig. 4.1 and Fig. 4.2 it can be seen that the overall structure in the γ-field is almost the same. However, when comparing the temperature figures it is seen that the non-conservative seems to be of much higher resolution. This was later found out to an effect of using the conservative equation. To compensate for this a third simulation was conducted with a smaller grid. The result is shown in Fig When comparing Fig. 4.3 and Fig. 4.2 it seems that the smaller grid made the two solution more alike, but since neither of the simulations did converge within the 100 iterations it could not be concluded that the conservative equation could be implemented for later use. 15
22 16 CHAPTER 4. RESULTS Figure 4.1: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. Figure 4.2: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the nonconservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
23 4.2. MAKING Q INDEPENDENT OF THE γ-area 17 Figure 4.3: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. 4.2 Making Q independent of the γ-area As described in Section it was believed that dividing Q with the area of the γ-field could confuse the algorithm, it was therefore decided to make a pair of simulations without this factor. The results are shown in Fig. 4.4 and Fig Strangely enough the removal of the area dependence, did not help at all. If the two results are compared it can be seen that they are much more different than those of Section 4.1. The area dependence was added again.
24 18 CHAPTER 4. RESULTS Figure 4.4: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. Q is independent the γ-area. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. Figure 4.5: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the nonconservative heat equation. Q is independent the γ-area. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
25 4.3. LOWERING DARCY S NUMBER Lowering Darcy s number In order to make the solutions converge faster further simulations were made with Da = In the hope that this would yield a much nicer outcome. The results are shown in Fig. 4.6 and Fig. 4.7 As expected the results were much better, if the results are Figure 4.6: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. compared, it can be seen that they now are very alike. Hence, it was concluded that the conservative heat equation could be used in further simulations.
26 20 CHAPTER 4. RESULTS Figure 4.7: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the nonconservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. 4.4 Varying C p Once it had been verified that the conservative heat equation could be used, it was possible to change the material parameters. The first test was made with a slightly varying heat capacity C p. The difference between the value for the liquid and for the material was only about 5%. This test was conducted both with Q dependent and independent of the γ-area. the results are shown in Fig and Fig In both cases the algorithm seems to work fine.
27 4.4. VARYING C P 21 Figure 4.8: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p 4186, ρ = 1000, size grid = and the conservative heat equation. Q is independent the γ-area. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. Figure 4.9: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p 4186, ρ = 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
28 22 CHAPTER 4. RESULTS 4.5 Varying k A pair of simulations were also conducted with the thermal conductivity being varied about 5% between the liquid and the material. Again the simulation was tested both with and without Q depending on the γ-area. The results are shown in Fig.?? and Fig.??. Figure 4.10: This figure shows the result after 100 iteration using the following parameters: Da = , k 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. Q is independent the γ-area. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. These results also show that the algorithm actually works better when Q is divided with the area of γ. Fig.?? also shows that the algorithm can work when k is changing.
29 4.6. VARYING ρ 23 Figure 4.11: This figure shows the result after 100 iteration using the following parameters: Da = , k 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. 4.6 Varying ρ For the last of the three parameters two test was also made, again both with and without Q depending on the area of γ. ρ is varied about 5%. The results are shown in Fig.?? and Fig.??. Again the results show that the algorithm works better when Q is dependent of the area of γ. Fig.?? also shows that the algorithm can work when ρ is changing.
30 24 CHAPTER 4. RESULTS Figure 4.12: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ 1000, size grid = and the conservative heat equation. Q is independent the γ-area. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. Figure 4.13: This figure shows the result after 100 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ 1000, size grid = and the conservative heat equation. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
31 4.7. VARYING THE INITIAL CONDITIONS Varying the initial conditions In order to see how results depend on the initial conditions of γ three simulations were made. All parameters was independent of γ. All the simulation were given the time to converge. The first simulation is shown in Fig This simulation was started with γ as a sine function of x. The second simulation is shown in Fig This simulation was started with γ as a sine function of y. The third simulation is shown in Fig This simulation was started with γ as a constant. Figure 4.14: This figure shows the result after 986 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The initial conditions for the γ-field is a sin(x)-wave. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature. If these three results are compared it is easy to see that the final result is dependent on the initial conditions for γ. Hence, one cannot be certain that the optimal structure is achieved when using the MDM-algorithm.
32 26 CHAPTER 4. RESULTS Figure 4.15: This figure shows the result after 145 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The initial conditions for the γ-field is a sin(y)-wave. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
33 4.7. VARYING THE INITIAL CONDITIONS 27 Figure 4.16: This figure shows the result after 900 iteration using the following parameters: Da = , k = 0.6, C p = 4186, ρ = 1000, size grid = and the conservative heat equation. The initial conditions for the γ-field is a constant γ 0. The top window shows the γ-field, in the middle the flow velocity is shown and in the bottom the temperature.
34 28 CHAPTER 4. RESULTS
35 Chapter 5 Further work During this work we did not have sufficient time to do all the worked we would have liked. We only managed to simulate a small change in each of the three parameters C p, kandρ, but not on the same time. Had there been more time, it would be of great benefit to add one parameter at a time, and watching how the algorithm handles it. Finally it should hopefully be possible to change all the material parameters, thereby giving the simulations a real physical relevance. It would also be interesting to see how the structures would change if the simulations were made with a much smaller mesh, this task could not be accomplished as the computational time increases greatly as the mesh is decreased. It would also be interesting to see how efficient self-designed-structures would be, and how the MDM-algorithm would try to improve those. 29
36 30 CHAPTER 5. FURTHER WORK
37 Chapter 6 Conclusion We have successfully implemented the physics of heat transfer in cross-flow in FEMLAB with MATLAB. We have learned the importance of starting from scratch adding more and more parameters making the solutions physical. We have gathered knowledge about the way FEMLAB iterates and how to take this into account when implementing the physics of the model. Setting γ 0 to the value zero resulted in a zero denominator which were avoided by adding a sine function to describe the initial value of gamma in the design field. Unfortunately the initial guess seemed to influence the solution. This suggests that the FEMLAB solution is not singular. The built in scaling function described in Chapter 3 needed to be enabled in order to hinder the pc to crash this was done by adding usolve in the script calling the femlin function. To prevent the algorithm from crashing due to changes in the material parameters, we changed the heat equation to the conservative version, allowing a less refined resolution. Acknowledgement We thank Fridolin Okkels, Laurits Olesen, Niels Asger Mortensen, and Henrik Bruus for interesting discussions on the functionality of the developed high-level program used for the implementation of the algorithm and the parameters making our model physical realistic for simulating the 2D Heat transfer model. 31
38 32 CHAPTER 6. CONCLUSION
39 Bibliography [1] Constructal multi-scale cylinders in cross-flow, T. Bello-Ochende and A. Bejan, International Journal of Heat and Mass Transfer. 48, (2005). [2] A high-level programming implementation of topology optimization applied to steadystate Navier-Stokes flow, Laurits Højgaard, Fridolin Okkels and Henrik Bruus, Phys. Rev. E. (11 Oct. 2004). 33
The 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 informationTOPOLOGY OPTIMIZATION OF NAVIER-STOKES FLOW IN MICROFLUIDICS
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ä, 24 28 July 2004
More informationThis section develops numerically and analytically the geometric optimisation of
7 CHAPTER 7: MATHEMATICAL OPTIMISATION OF LAMINAR-FORCED CONVECTION HEAT TRANSFER THROUGH A VASCULARISED SOLID WITH COOLING CHANNELS 5 7.1. INTRODUCTION This section develops numerically and analytically
More informationThis chapter focuses on the study of the numerical approximation of threedimensional
6 CHAPTER 6: NUMERICAL OPTIMISATION OF CONJUGATE HEAT TRANSFER IN COOLING CHANNELS WITH DIFFERENT CROSS-SECTIONAL SHAPES 3, 4 6.1. INTRODUCTION This chapter focuses on the study of the numerical approximation
More informationDESIGN 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 informationProject #1 Internal flow with thermal convection
Project #1 Internal flow with thermal convection MAE 494/598, Fall 2017, Project 1 (20 points) Hard copy of report is due at the start of class on the due date. The rules on collaboration will be released
More informationThe effect of Entry Region on Thermal Field
The effect of Entry Region on Thermal Field Flow Fractionation Nick Cox Supervised by Professor Bruce Finlayson University of Washington Department of Chemical Engineering June 6, 2007 Abstract Achieving
More informationMAE 598 Project #1 Jeremiah Dwight
MAE 598 Project #1 Jeremiah Dwight OVERVIEW A simple hot water tank, illustrated in Figures 1 through 3 below, consists of a main cylindrical tank and two small side pipes for the inlet and outlet. All
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 informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
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 informationLecture 2: Hydrodynamics at milli micrometer scale
1 at milli micrometer scale Introduction Flows at milli and micro meter scales are found in various fields, used for several processes and open up possibilities for new applications: Injection Engineering
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 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 informationNUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A CAPILLARY TUBE
Proceedings of the Asian Conference on Thermal Sciences 2017, 1st ACTS March 26-30, 2017, Jeju Island, Korea ACTS-P00786 NUMERICAL INVESTIGATION OF THERMOCAPILLARY INDUCED MOTION OF A LIQUID SLUG IN A
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationMaximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection
Maximum Heat Transfer Density From Finned Tubes Cooled By Natural Convection Ahmed Waheed Mustafa 1 Mays Munir Ismael 2 AL-Nahrain University College of Engineering Mechanical Engineering Department ahmedwah@eng.nahrainuniv.edu.iq
More informationFluid Mechanics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S. K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 15 Conservation Equations in Fluid Flow Part III Good afternoon. I welcome you all
More informationFINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT OF A HEATED SQUARE HOLLOW CYLINDER IN A LID-DRIVEN RECTANGULAR ENCLOSURE
Proceedings of the International Conference on Mechanical Engineering 2011 (ICME2011) 18-20 December 2011, Dhaka, Bangladesh ICME11-TH-014 FINITE ELEMENT ANALYSIS OF MIXED CONVECTION HEAT TRANSFER ENHANCEMENT
More informationUNIT II CONVECTION HEAT TRANSFER
UNIT II CONVECTION HEAT TRANSFER Convection is the mode of heat transfer between a surface and a fluid moving over it. The energy transfer in convection is predominately due to the bulk motion of the fluid
More informationHeat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay
Heat and Mass Transfer Prof. S.P. Sukhatme Department of Mechanical Engineering Indian Institute of Technology, Bombay Lecture No. 18 Forced Convection-1 Welcome. We now begin our study of forced convection
More informationTransport by convection. Coupling convection-diffusion
Transport by convection. Coupling convection-diffusion 24 mars 2017 1 When can we neglect diffusion? When the Peclet number is not very small we cannot ignore the convection term in the transport equation.
More informationSimulation and improvement of the ventilation of a welding workshop using a Finite volume scheme code
1 st. Annual (National) Conference on Industrial Ventilation-IVC2010 Feb 24-25, 2010, Sharif University of Technology, Tehran, Iran IVC2010 Simulation and improvement of the ventilation of a welding workshop
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 informationBasic Fluid Mechanics
Basic Fluid Mechanics Chapter 3B: Conservation of Mass C3B: Conservation of Mass 1 3.2 Governing Equations There are two basic types of governing equations that we will encounter in this course Differential
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 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 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 informationCooling of a multi-chip power module
Cooling of a multi-chip power module G. CAMMARAA, G. PERONE Department of Industrial Engineering University of Catania Viale A. Doria 6, 953 Catania IALY gcamma@dii.unict.it, gpetrone@dii.unict.it Abstract:
More informationCOURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour. Basic Equations in fluid Dynamics
COURSE NUMBER: ME 321 Fluid Mechanics I 3 credit hour Basic Equations in fluid Dynamics Course teacher Dr. M. Mahbubur Razzaque Professor Department of Mechanical Engineering BUET 1 Description of Fluid
More informationMixing in Flow Devices:
Mixing in Flow Devices: Spiral microchannels in two and three dimensions Prepared by Ha Dinh Mentor: Professor Emeritus Bruce A. Finlayson Department of Chemical Engineering University of Washington June
More informationNumerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel
Numerical analysis of fluid flow and heat transfer in 2D sinusoidal wavy channel Arunanshu Chakravarty 1* 1 CTU in Prague, Faculty of Mechanical Engineering, Department of Process Engineering,Technická
More informationApplied CFD Project 1. Christopher Light MAE 598
Applied CFD Project 1 Christopher Light MAE 598 October 5, 2017 Task 1 The hot water tank shown in Fig 1 is used for analysis of cool water flow with the heat from a hot plate at the bottom. For all tasks,
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 information4.1 Derivation and Boundary Conditions for Non-Nipped Interfaces
Chapter 4 Roller-Web Interface Finite Difference Model The end goal of this project is to allow the correct specification of a roller-heater system given a general set of customer requirements. Often the
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 informationNicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA.
Transport Limitations in Thermal Diffusion Nicholas Cox, Pawel Drapala, and Bruce F. Finlayson Department of Chemical Engineering, University of Washington, Seattle, WA, USA Abstract Numerical simulations
More informationSimplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis
1 Portál pre odborné publikovanie ISSN 1338-0087 Simplified Model of WWER-440 Fuel Assembly for ThermoHydraulic Analysis Jakubec Jakub Elektrotechnika 13.02.2013 This work deals with thermo-hydraulic processes
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 9 (214 ) 599 64 1th International Conference on Mechanical Engineering, ICME 213 Validation criteria for DNS of turbulent heat
More information(Refer Slide Time 1:25)
Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 2 Lecture - 24 Transient Response of Pressure Transducers
More informationContents. Microfluidics - Jens Ducrée Physics: Laminar and Turbulent Flow 1
Contents 1. Introduction 2. Fluids 3. Physics of Microfluidic Systems 4. Microfabrication Technologies 5. Flow Control 6. Micropumps 7. Sensors 8. Ink-Jet Technology 9. Liquid Handling 10.Microarrays 11.Microreactors
More informationTHE EFFECTS OF LONGITUDINAL RIBS ON ENTROPY GENERATION FOR LAMINAR FORCED CONVECTION IN A MICROCHANNEL
THE EFFECTS OF LONGITUDINAL RIBS ON ENTROPY GENERATION FOR LAMINAR FORCED CONVECTION IN A MICROCHANNEL Nader POURMAHMOUD, Hosseinali SOLTANIPOUR *1,, Iraj MIRZAEE Department of Mechanical Engineering,
More informationLectures on Applied Reactor Technology and Nuclear Power Safety. Lecture No 6
Lectures on Nuclear Power Safety Lecture No 6 Title: Introduction to Thermal-Hydraulic Analysis of Nuclear Reactor Cores Department of Energy Technology KTH Spring 2005 Slide No 1 Outline of the Lecture
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 informationAcoustic Streaming Driven Mixing
Acoustic Streaming Driven Mixing Nitesh Nama, Po-Hsun Huang, Francesco Costanzo, and Tony Jun Huang Department of Engineering Science and Mechanics The Pennsylvania State University, State College, PA,
More informationCONVECTIVE HEAT TRANSFER
CONVECTIVE HEAT TRANSFER Mohammad Goharkhah Department of Mechanical Engineering, Sahand Unversity of Technology, Tabriz, Iran CHAPTER 3 LAMINAR BOUNDARY LAYER FLOW LAMINAR BOUNDARY LAYER FLOW Boundary
More informationCHAPTER 4. Basics of Fluid Dynamics
CHAPTER 4 Basics of Fluid Dynamics What is a fluid? A fluid is a substance that can flow, has no fixed shape, and offers little resistance to an external stress In a fluid the constituent particles (atoms,
More informationFluid Mechanics II Viscosity and shear stresses
Fluid Mechanics II Viscosity and shear stresses Shear stresses in a Newtonian fluid A fluid at rest can not resist shearing forces. Under the action of such forces it deforms continuously, however small
More informationOE4625 Dredge Pumps and Slurry Transport. Vaclav Matousek October 13, 2004
OE465 Vaclav Matousek October 13, 004 1 Dredge Vermelding Pumps onderdeel and Slurry organisatie Transport OE465 Vaclav Matousek October 13, 004 Dredge Vermelding Pumps onderdeel and Slurry organisatie
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 49 Introduction to Turbulent Flow part -II Good morning I welcome you all to this
More informationMOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow
TRANSPORT PHENOMENA MOMENTUM TRANSPORT Velocity Distributions in Turbulent Flow Introduction to Turbulent Flow 1. Comparisons of laminar and turbulent flows 2. Time-smoothed equations of change for incompressible
More informationA NUMERICAL APPROACH FOR ESTIMATING THE ENTROPY GENERATION IN FLAT HEAT PIPES
A NUMERICAL APPROACH FOR ESTIMATING THE ENTROPY GENERATION IN FLAT HEAT PIPES Dr. Mahesh Kumar. P Department of Mechanical Engineering Govt College of Engineering, Kannur Parassinikkadavu (P.O), Kannur,
More informationLesson 6 Review of fundamentals: Fluid flow
Lesson 6 Review of fundamentals: Fluid flow The specific objective of this lesson is to conduct a brief review of the fundamentals of fluid flow and present: A general equation for conservation of mass
More informationTurbulent Natural Convection in an Enclosure with Colliding Boundary Layers
Turbulent Natural Convection in an Enclosure with Colliding Boundary Layers Abstract Mutuguta John Wanau 1* 1. School of Pure and Applied Sciences, Murang a University of Technology, P.O box 75-10200,
More informationNatural Convection in Vertical Channels with Porous Media and Adiabatic Extensions
Natural Convection in Vertical Channels with Porous Media and Adiabatic Extensions Assunta Andreozzi 1,a, Bernardo Buonomo 2,b, Oronzio Manca 2,c and Sergio Nardini 2,d 1 DETEC, Università degli Studi
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 informationPhone: , For Educational Use. SOFTbank E-Book Center, Tehran. Fundamentals of Heat Transfer. René Reyes Mazzoco
8 Fundamentals of Heat Transfer René Reyes Mazzoco Universidad de las Américas Puebla, Cholula, Mexico 1 HEAT TRANSFER MECHANISMS 1.1 Conduction Conduction heat transfer is explained through the molecular
More information3.8 The First Law of Thermodynamics and the Energy Equation
CEE 3310 Control Volume Analysis, Sep 30, 2011 65 Review Conservation of angular momentum 1-D form ( r F )ext = [ˆ ] ( r v)d + ( r v) out ṁ out ( r v) in ṁ in t CV 3.8 The First Law of Thermodynamics and
More informationFluid Dynamics and Balance Equations for Reacting Flows
Fluid Dynamics and Balance Equations for Reacting Flows Combustion Summer School 2018 Prof. Dr.-Ing. Heinz Pitsch Balance Equations Basics: equations of continuum mechanics balance equations for mass and
More informationFluid Dynamics Exercises and questions for the course
Fluid Dynamics Exercises and questions for the course January 15, 2014 A two dimensional flow field characterised by the following velocity components in polar coordinates is called a free vortex: u r
More informationChapter 5. The Differential Forms of the Fundamental Laws
Chapter 5 The Differential Forms of the Fundamental Laws 1 5.1 Introduction Two primary methods in deriving the differential forms of fundamental laws: Gauss s Theorem: Allows area integrals of the equations
More informationCHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE
CHAPTER 7 NUMERICAL MODELLING OF A SPIRAL HEAT EXCHANGER USING CFD TECHNIQUE In this chapter, the governing equations for the proposed numerical model with discretisation methods are presented. Spiral
More informationFLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study
FLOW MALDISTRIBUTION IN A SIMPLIFIED PLATE HEAT EXCHANGER MODEL - A Numerical Study Nityanand Pawar Mechanical Engineering, Sardar Patel College of Engineering, Mumbai, Maharashtra, India nitya.pawar@gmail.com
More informationMath background. Physics. Simulation. Related phenomena. Frontiers in graphics. Rigid fluids
Fluid dynamics Math background Physics Simulation Related phenomena Frontiers in graphics Rigid fluids Fields Domain Ω R2 Scalar field f :Ω R Vector field f : Ω R2 Types of derivatives Derivatives measure
More informationShell/Integral Balances (SIB)
Shell/Integral Balances (SIB) Shell/Integral Balances Shell or integral (macroscopic) balances are often relatively simple to solve, both conceptually and mechanically, as only limited data is necessary.
More informationTheoretical microfluidics
Lecture notes third edition, fall 2006 Theoretical microfluidics Henrik Bruus + + + + MIC Department of Micro and Nanotechnology Technical University of Denmark ii Preface In the fall 2003 MIC launched
More informationFluid Animation. Christopher Batty November 17, 2011
Fluid Animation Christopher Batty November 17, 2011 What distinguishes fluids? What distinguishes fluids? No preferred shape Always flows when force is applied Deforms to fit its container Internal forces
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 informationOverall Heat Transfer Coefficient
Overall Heat Transfer Coefficient A heat exchanger typically involves two flowing fluids separated by a solid wall. Heat is first transferred from the hot fluid to the wall by convection, through the wall
More informationMass Transfer in a Stirred Batch Reactor
Mass Transfer in a Stirred Batch Reactor In many processes, efficient reactor usage goes hand in hand with efficient mixing. The ability to accurately examine the effects of impeller placement, speed,
More informationLecture 30 Review of Fluid Flow and Heat Transfer
Objectives In this lecture you will learn the following We shall summarise the principles used in fluid mechanics and heat transfer. It is assumed that the student has already been exposed to courses in
More informationNumerical Heat and Mass Transfer
Master Degree in Mechanical Engineering Numerical Heat and Mass Transfer 15-Convective Heat Transfer Fausto Arpino f.arpino@unicas.it Introduction In conduction problems the convection entered the analysis
More informationNumerical Study of PCM Melting in Evacuated Solar Collector Storage System
Numerical Study of PCM Melting in Evacuated Collector Storage System MOHD KHAIRUL ANUAR SHARIF, SOHIF MAT, MOHD AFZANIZAM MOHD ROSLI, KAMARUZZAMAN SOPIAN, MOHD YUSOF SULAIMAN, A. A. Al-abidi. Energy Research
More informationEFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION
EFFECT OF DISTRIBUTION OF VOLUMETRIC HEAT GENERATION ON MODERATOR TEMPERATURE DISTRIBUTION A. K. Kansal, P. Suryanarayana, N. K. Maheshwari Reactor Engineering Division, Bhabha Atomic Research Centre,
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 informationA Case Study of Hooke s Law and the Variation of Spring Constant
A Case Study of Hooke s Law and the Variation of Spring Constant James Freeman PHYS 113; Professor Khalifa Freeman 2 Introduction Bungee jumping is more than an extreme sport when viewed through the lenses
More informationHomework #4 Solution. μ 1. μ 2
Homework #4 Solution 4.20 in Middleman We have two viscous liquids that are immiscible (e.g. water and oil), layered between two solid surfaces, where the top boundary is translating: y = B y = kb y =
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 information2.3 The Turbulent Flat Plate Boundary Layer
Canonical Turbulent Flows 19 2.3 The Turbulent Flat Plate Boundary Layer The turbulent flat plate boundary layer (BL) is a particular case of the general class of flows known as boundary layer flows. The
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 informationTheoretical microfluidics
Lecture notes second edition, fall 2005 Theoretical microfluidics Henrik Bruus + + + + MIC Department of Micro and Nanotechnology Technical University of Denmark ii Preface In the fall 2003 MIC launched
More informationSupplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer
Supplementary Information for Engineering and Analysis of Surface Interactions in a Microfluidic Herringbone Micromixer Thomas P. Forbes and Jason G. Kralj National Institute of Standards and Technology,
More informationComputational Astrophysics
Computational Astrophysics Lecture 1: Introduction to numerical methods Lecture 2:The SPH formulation Lecture 3: Construction of SPH smoothing functions Lecture 4: SPH for general dynamic flow Lecture
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 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 informationMOMENTUM PRINCIPLE. Review: Last time, we derived the Reynolds Transport Theorem: Chapter 6. where B is any extensive property (proportional to mass),
Chapter 6 MOMENTUM PRINCIPLE Review: Last time, we derived the Reynolds Transport Theorem: where B is any extensive property (proportional to mass), and b is the corresponding intensive property (B / m
More informationB.1 NAVIER STOKES EQUATION AND REYNOLDS NUMBER. = UL ν. Re = U ρ f L μ
APPENDIX B FLUID DYNAMICS This section is a brief introduction to fluid dynamics. Historically, a simplified concept of the boundary layer, the unstirred water layer, has been operationally used in the
More informationME 144: Heat Transfer Introduction to Convection. J. M. Meyers
ME 144: Heat Transfer Introduction to Convection Introductory Remarks Convection heat transfer differs from diffusion heat transfer in that a bulk fluid motion is present which augments the overall heat
More informationEffect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process
, pp. 844 850 Effect of Static Magnetic Field Application on the Mass Transfer in Sequence Slab Continuous Casting Process Baokuan LI and Fumitaka TSUKIHASHI 1) Department of Thermal Engineering, The School
More informationV (r,t) = i ˆ u( x, y,z,t) + ˆ j v( x, y,z,t) + k ˆ w( x, y, z,t)
IV. DIFFERENTIAL RELATIONS FOR A FLUID PARTICLE This chapter presents the development and application of the basic differential equations of fluid motion. Simplifications in the general equations and common
More informationUNIVERSITY OF OSLO. Faculty of Mathematics and Natural Sciences. MEK4300/9300 Viscous flow og turbulence
UNIVERSITY OF OSLO Faculty of Mathematics and Natural Sciences Examination in: Day of examination: Friday 15. June 212 Examination hours: 9. 13. This problem set consists of 5 pages. Appendices: Permitted
More informationMATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall
MATLAB Solution of Flow and Heat Transfer through a Porous Cooling Channel and the Conjugate Heat Transfer in the Surrounding Wall James Cherry, Mehmet Sözen Grand Valley State University, cherryj1@gmail.com,
More information( ) II. Why statistics? SIO 221B, Rudnick adapted from Davis!1
SIO 221B, Rudnick adapted from Davis 1 II. Why statistics? This course addresses the analysis of oceanographic observations and, as they become more complex, ocean models. Much of the material involves
More information6. Laminar and turbulent boundary layers
6. Laminar and turbulent boundary layers John Richard Thome 8 avril 2008 John Richard Thome (LTCM - SGM - EPFL) Heat transfer - Convection 8 avril 2008 1 / 34 6.1 Some introductory ideas Figure 6.1 A boundary
More informationDEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M.
DEVELOPED LAMINAR FLOW IN PIPE USING COMPUTATIONAL FLUID DYNAMICS M. Sahu 1, Kishanjit Kumar Khatua and Kanhu Charan Patra 3, T. Naik 4 1, &3 Department of Civil Engineering, National Institute of technology,
More informationISSN Article
Entropy 23, 5, 28-299; doi:.339/e5628 OPEN ACCESS entropy ISSN 99-43 www.mdpi.com/journal/entropy Article Analysis of Entropy Generation Rate in an Unsteady Porous Channel Flow with Navier Slip and Convective
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 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 informationINTRODUCTION TO FLUID MECHANICS June 27, 2013
INTRODUCTION TO FLUID MECHANICS June 27, 2013 PROBLEM 3 (1 hour) A perfect liquid of constant density ρ and constant viscosity µ fills the space between two infinite parallel walls separated by a distance
More information