Parallel domain decomposition meshless modeling of dilute convection-diffusion of species
|
|
- Leon Mosley
- 5 years ago
- Views:
Transcription
1 Boundary Elements XXVII 79 Parallel domain decomposition meshless modeling of dilute convection-diffusion of species E. Divo 1,. Kassab 2 & Z. El Zahab 2 1 Department of Engineering Technology, University of Central Florida, Orlando, FLLLLLLLlorida L , U.S.. 2 Department of Mechanical, Materials, and erospace Engineering, University of Central Florida, Orlando, Florida , U.S.. bstract meshless formulation and computational modeling of non-reacting dilute species transport is presented in this paper. Two species are considered for the problems in the example section: methane (CH 4 ) and air. parallel domain decomposition meshless method based on Multiquadric RBF collocation is developed herein. This approach reduces conditioning numbers of the resulting coefficient matrices while accelerating the solution process and reducing the memory requirements. Numerical examples are presented to validate the approach by comparing the meshless solutions to CFD-FVM solutions provided by the commercial package, Fluent Introduction The underlying algorithms of classical numerical methods such as finite element methods (FEM) and finite volume methods (FVM) have matured and afford the analyst the ability to routinely and accurately analyze a wide range of complex phenomena. However, the most commonly used methods require domain meshing due to the fact that rely on locally-supported polynomial expansions for the unknown field variables. Mesh generation is far from automated and in particular in computational fluid dynamics (CFD) applications require a large amount of effort and human interaction to generate meshes of sufficient quality to produce acceptable results.
2 80 Boundary Elements XXVII In recent years, meshless or mesh-free methods, reminiscent of spectral tech- niques [8] in their algorithmic structure, have been under intensive development in an attempt to overcome the difficulty of meshing the solution domain. The term mesh-free method refers to any numerical technique that relies on a set of unstructured nodes, that is, without pre-defined connectivity or relationship among the nodes, to solve a particular field problem. The key advantage of mesh-free methods as compared to traditional mesh-based methods is the expected large savings in labor dedicated to mesh generation. Unlike spectral methods that rely on global orthogonal expansion using typically, Legendre or Chebyshev polynomials, the meshless method of interest in this paper belongs to a class of that techniques that adopts radial basis function (RBFs) [9] expansions for the dependent variable(s). Domain decomposition [5, 10] has been shown to effectively reduce storage and RM access requirements as well as significantly reduce matrix conditioning in applications to scalar problems. In this paper, we develop a parallel domain decomposition collocation meshless method for dilute convection-diffusion of species. The species densities is approximated using the inverse Multiquadrics expansion RBF. Two cases are considered: (1) a parallel plates channel problem with and without domain decomposition and (2) a T-shape channel decomposed into five sub-domains. 2 Meshless formulation of dilute convection-diffusion of species Transient incompressible fluid flow is governed by the coupled set of non-linear partial differential equations, and these are the continuity, Navier-Stokes, and energy transport equations which are collectively often referred to as the Navier- Stokes equations: V =0 ρ V t + ρ( V ) V = p + µ 2 V + ρ f ρc T t + ρc( V )T = k 2 T +Φ where V is the flow velocity vector, ρ is the bulk density of the flow, p is the field pressure, µ is the absolute fluid viscosity, f is a specific body force, c is the fluid specific heat, T is the field temperature, k is the fluid thermal conductivity,and Φ is the viscous heat dissipation. ll field variables are functions of space (x) and time (t), in a fixed domain Ω surrounded by a closed boundary Γ. If a new species () is introduced into the bulk flow in dilute form (small concentration), then the above equations can be solved independently under the assumption that the overall flow field remains approximately incompressible. The transport of such species is governed by the following differential equation, coupled one-way to the Navier-Stokes equations (Eqn. (1)): ρ +( V t )ρ = D B 2 ρ (2) (1)
3 Boundary Elements XXVII 81 Where ρ is the partial density of species, andd B is the binary diffusivity of species into the bulk flow. Equation (2) can be solved independently if the bulk flow velocity field, V, is predetermined. In this paper we use a pressurecorrection based meshless method as described in [12]. Expressing the species density transport Eqn. (2) in two-dimensional Cartesian form, see [13], yields: ρ t with generalized boundary conditions: + u ρ x + v ρ y = D B 2 ρ (3) β ρ n + γρ = σ (4) where the space-time dependency (x, y, t) has been omitted for simplicity. Here u and v are the horizontal and vertical velocity components of the bulk flow respectively. The parameters β, γ,andσ are specified and determine the type of boundary condition. The meshless approximation begins by globally expanding the species density over a number NB + NI of independent points in the domain Ω(x, y) and on the boundary Γ(x,y)as: ρ (x, y, t) = NB+NI j=1 α j (t)χ j (x, y) (5) where NB: boundary expansion points, NI: internal expansion points, α j (t):transient expansion coefficient, and χ j (x, y): global expansion functions. mong the family of global expansion functions the radial-basis functions (RBF) have shown to effectively model scattered data, see [6, 9]. The Multiquadric RBF, see [5, 7], stands out in comparison to other RBF s for its ability to accurately interpolate field variables over generalized domains and point distributions. The Multiquadric RBF s are defined as: χ j (x, y) = [ r 2 j (x, y)+c 2] n 3 2 (6) where n: exponent parameter (positive integer), c: shape parameter, and r j : Euclidean distance from (x, y) to (x j,y j ). In Eqn. (5), the profile of the expansion flattens as the values of n and c increase which will give a smooth and accurate global representation of the field variable. Nevertheless, the drawback of having a flat representation of the field variable lies in the linear dependence that arises among the system of linear equations resulting in a high matrix conditioning number [7]. Thus, it is necessary to choose n and c values to yield an admissible matrix conditioning number. For the current formulation, n is taken as 1 yielding the socalled inverse Multiquadric. The species density transport Eqn. (3) can be numerically evaluated between time steps p and p +1by replacing the time derivative with a first order backward differencing approximation as:
4 82 Boundary Elements XXVII ρ p+1 ρp + u t p+1 ρp+1 ρp+1 + vp+1 = D B 2 ρ p+1 (7) x y or, placing the current time step values, p +1, on the left-hand side yields: D B 2 ρ p+1 ρp+1 up+1 ρp+1 vp+1 x y ρp+1 t = ρp t t this point, the global expansion of the species density in Eqn. (5) can be introduced. Notice that the spatial differential operators and the summations are interchanged and the differential operators are applied directly to the expansion functions χ j (x, y) as: (8) NB+NI j=1 α p+1 j (D B 2 χ j u p+1 χ j x χ vp+1 j y χ j t )= ρp t (9) where ρ p is known from the previous time step starting at ρ0. The same procedure can be applied to the boundary condition relation in Eqn. (4) as: NB+NI j=1 α p+1 j (β χ j n + γχ j)=σ (10) The coefficients α p+1 j corresponding to the time step p +1can be found by collocating the boundary condition expansion in Eqn. (10) at points i =1..NB and the governing equation expansion in Eqn. (9) at points i = NB+1..NB+NI. This process yields the following set of linear simultaneous equations: [ ] ψj { } { } (x i,y i ) α p+1 σ(xi,y η j (x i,y i ) j = i ) 1 t ρp (x i,y i ) where: NB+NI NB+NI ψ j (x i,y i )=β χ j(x i,y i ) n η j (x i,y i )=D B 2 χ j (x i,y i ) u p+1 (x i,y i ) χ j(x i,y i ) x NB+NI 1 + γχ j (x i,y i ) v p+1 (x i,y i ) χ j(x i,y i ) y NB+NI 1 (11) χ j(x i,y i ) t (12) the left-hand side coefficient matrix is fully-populated and normally found to have very large conditioning numbers. Both of these features are major draw-backs for the implementation of this technique to large-scale multi-dimensional problems. Evidently, pre-conditioning this matrix is a mandatory step in the solution of realistic systems, see [5, 10] for more detail. This issue will be addressed next and the salient features of a typical coefficient matrix resulting from mesh-free expansions will be illustrated.
5 Boundary Elements XXVII 83 Boundary Point Internal Point Figure 1: Problem domain and typical boundary and internal point collocation. 3 Parallel domain decomposition implementation In order to pre-condition a set of algebraic equations, one may pre-multiply the system of equations by a matrix that closely approximates the inverse of the coefficient matrix and that is also cheaply evaluated from the computational standpoint. Incomplete LU pre-conditioning is such an approach. n alternative to mitigate the effect of large conditioning numbers is to reduce the number of degrees of freedom in the algebraic system. This can be accomplished by decomposing the problem domain Ω into multiple sub-domains Ω k separated by artificially created interfaces Γ I, see [10]. The problem is started from initial conditions which are used as boundary conditions at the artificially created interfaces at time level p =0.Once each sub-domain is a well-posed boundary value problem, the standard collocation approach described in the previous section may be followed to render independent expansions at each sub-domain. simple iteration implementation can follow to force continuity conditions of the field variable and normal derivatives at the interfaces. If a single-region point collocation approach as the one shown in Fig. 1 is followed, a system of simultaneous linear equations of size N,whereN = NB+NI, will be formed as: Ω [] N N {α} N 1 = {b} N 1 (13) where {α} represents the vector of coefficients for the global expansion. The number of floating point operations required to arrive at the algebraic system above is proportional to N 2 as well as direct memory allocation also proportional to N 2. The solution to the algebraic system for the coefficients {α} can be performed using a direct solution method such as LU decomposition requiring floating point operations proportional to N 3 or an indirect method such as Bi-conjugate Gradient or Generalized Minimum Residual (GMRes) which, in general, require floating point operations proportional to N 2 to achieve convergence. On the other hand, if a the multi-region iteration approach were to be adopted, the domain is divided into K sub-domains and each one is independently filled with collocation points. See Fig. 2 for a typical domain decomposition with K =4.
6 84 Boundary Elements XXVII 3(1) 3(2) 3(3) 3(4) 4(1) I 2(1) + I 4(2) II 2(2) + II 4(3) III 2(3) + III 4(4) 2(4) (1) 1(2) 1(3) 1(4) Figure 2: Typical domain decomposition and collocation. The information between neighboring sub-domains separated by an interface will be passed through using their respective radial-basis expansions. Therefore, the collocation points at the artificially created interfaces do not need to coincide. The boundary value problem will now be solved independently over each subdomain starting by imposing the initial condition at the artificially created interface. For instance, the boundary value problem of sub-domain Ω 1 is transformed into the algebraic analogue of corresponding influence coefficient matrices and nodal boundary values as: Ω 1 [] n n {α} n 1 = {b} n 1 (14) The composition of this algebraic system requires floating point operations proportional to the square of the number of boundary and internal collocation points (n 2 ) in the sub-domain Ω 1. Memory allocation is also proportional to the same number (n 2 ). This new proportionality number n is roughly equivalent to: n N K (15) as long as the collocation density along the boundaries is equivalent to the collocation density in the interior and the artificially created interfaces are selected to perpendicularly sub-section a longitudinal direction. Direct memory allocation requirement for later algebraic manipulation is now reduced to a proportion of n 2 as the system matrices can easily be stored on disk for later use after the boundary value problems on remaining sub-domains have been effectively solved. For the example shown here where the number of subdomains is K = 4 the new proportionality value n is approximately equal to n N/4. This simple multi-region example reduces the memory requirements to about n 2 /N 2 =(1/16) = 6.25% of the single-region approach. The solution of all of the four resulting algebraic systems for all four sub-domains require a number of floating point operations proportional to K (n 3 /N 3 )=4 (1/64) = 6.25% of the single-region approach if a direct algebraic
7 Boundary Elements XXVII 85 solution method is employed, or a number of floating point operations proportional to K (n 2 /N 2 )=4 (1/16) = 25% of the single-region approach if an indirect algebraic solution method is employed. For both, floating point operation count and direct memory requirement the reduction is dramatic. However, as the first set of solutions for the sub-domains were obtained using guessed boundary conditions along the interfaces, the global solution needs to follow an iteration process and convergence criteria. The more significant reduction is revealed in the RM requirements as only the memory needs for one of the sub-domains must be allocated at a time. The rest of the system matrices for the remaining sub-domains are temporarily stored on disk until access and manipulation is required. With respect to the algebraic solution of the system of equations for each subdomain of the domain decomposition, if a direct approach such as the LU decomposition method is employed for all sub-domains, the LU factors of the coefficient matrices for all sub-domains can be computed only once at the first iteration step and stored on disc memory for later use during the iteration process, unless the sub-domains are kept small enough to store the LU factors in RM on each processor when the algorithm is implemented in parallel. Thereon, only a forward and a backward substitution will be required at each iteration to arrive at a solution for the system in hand. This feature allows a significant reduction in the operational count through the iteration process until convergence is achieved, as only a number of floating point operations proportional to n as opposed to n 3 is required at each iteration step. To this computation time is added the access to disk (if the option of storing the factors in RM is not feasible) at each iteration step which is usually larger than the access to RM. 4 Time progression iteration process Once each boundary value problem is independently solved in each sub-domain, a time progression iteration process must follow to continuously ensure continuity of the field variable and its normal derivative across the interfaces. This is accomplished by setting up a process in which each interface is treated differently on both sides: on one side a first kind boundary condition is applied while at the other side a second kind boundary condition is applied. The resulting boundary values are then averaged with their corresponding values across the interfaces. This process is illustrated below in Fig. 3 for a single point on the interface between sub-domain Ω 1 and Ω 2. In Fig. 3 the imposed boundary conditions are ρ p on the left-hand side of the interface and ρ p / n on the right-hand side. fter the solution of the boundary value problem on each sub-domain the resulting (provisional) values are ρ p+1 / n on the left-hand side of the interface and ρp+1 on the right-hand side. These resulting values are averaged to generate the boundary conditions for the
8 86 Boundary Elements XXVII I 2(1) I 4(2) I 2(1) I 4(2) I 2(1) I 4(2) p ˆ ˆ n p p n 1 p 1 ˆ 1 p p ˆ n 1 Figure 3: Iteration process averaging across an interface. next time step iteration as: ρ p+1 ρ p+1 n = ρp + ρp+1 2 = ρp n 1 ( ρp 2 n + ρp+1 n ) (16) ensuring continuity of the field variable and its normal derivative at each time step iteration. Notice that in order to achieve truly time-accurate analysis, sublevel iterations must be performed at each time step. The iteration process may be stopped once steady-state is achieved or at any time level. Steady-state may be verified by defining a residual as: r p+1 = 1 NB + NI NB+NI i=1 (D B 2 ρ p+1 i u p+1 i ρ p+1 i x vp+1 i ρ p+1 i y )2 (17) and comparing it to that one of the previous time step, p. Once the residual difference reaches a preset convergence criterion or an oscillative behavior, the iteration process may be stopped and the resulting field values approximate the steady-state solution. 5 Numerical results Two examples are now presented to validate the formulation by comparing the meshless results to Finite Volume Method (FVM) solutions provided by the commercial CFD package Fluent 6.1. The first example consists of a 10 2cm channel through which air enters at 0.25m/s. leak of Methane enters through a section of 0.2cm at the center of the air inlet. The FVM mesh, 1-region, and 5-region meshless point distribution is presented in Fig. 4 along with the corresponding velocity magnitudes. The residual evolution for the 1-region meshless solution is shown in Fig. 5 along with the resulting Methane partial density distribution. The second example consists of a T-valve with 3mm inlets and outlets. ir enters at 0.25m/s through the two inlets while a leak of Methane enters through a section of 0.34mm at the center of the left wall. The FVM mesh and 5-region meshless point distribution is presented in Fig. along with the corresponding
9 Boundary Elements XXVII 87 Figure 4: FVM, 1-region meshless, 5-region meshless, point distribution and velocity magnitude contours Mass X-Momentum Y-Momentum Energy Species Step Figure 5: Residual evolution of the 1-region meshless solution and FVM, 1-region meshless, 5-region meshless Methane partial density distribution ρ. velocity magnitudes. The FVM and 5-region meshless pressure distribution is presented in Fig. 6 along with the resulting Methane partial density distribution. The residual evolution for the 1-region meshless solution is shown in Fig Conclusion meshless formulation and computational modeling of non-reacting dilute species transport is presented in this paper. Two species are considered for the problems in the example section: Methane (CH4 ) and air. parallel domain decomposition meshless method based on Multiquadric RBF collocation is developed herein. This approach reduces conditioning numbers of the resulting coefficient matrices while accelerating the solution process and reducing the memory requirements. Numerical examples are presented to validate the approach by comparing the meshless solutions to CFD-FVM solutions provided by the commercial package, Fluent 6.1.
10 88 Boundary Elements XXVII Figure 6: FVM and 5-region meshless, point distribution, velocity magnitude contours, and Methane partial density distribution ρ. 102 Mass X-Momentum Y-Momentum Energy Species Step Figure 7: Residual evolution of the 5-region meshless solution. References [1] Belytscho, T., Lu, Y.Y., and Gu, L., Element-free Galerkin methods, Int. J. Num. Methods, Vol. 37, 1994, pp [2] tluri, S.N. and Zhu, T., new meshless local Petrov- Galerkin (MLPG) approach in computational mechanics, Computational Mechanics, Vol. 22, 1998, pp
11 Boundary Elements XXVII 89 [3] Melenk, J.M. and Babuska, I., The partition of unity finite element method: basic theory and application, Comp. Meth. ppl. Mechanics and Eng., Vol. 139, 1996, pp [4] Kansa, E.J., Multiquadrics- a scattered data approximation scheme with applications to computational fluid dynamics I, Comp. Math. ppl., Vol. 19, 1990, pp [5] tluri,s.n. and Zhu, T., new meshless local Petrov-Galerkin (MLPG) approach in computational mechanics, Computational Mechanics, Vol. 22, 1998, pp [6] Franke, R., Scattered data interpolation: Test of some methods, Math Comput., Vol. 38, 1982, pp [7] Cheng,.H.-D., Golberg, M.., Kansa, E.J., Zammito, G., Exponential Convergence and H-c Multiquadric Collocation Method for Partial Differential Equations, Numerical Methods in P.D.E., Vol. 19, No. 5, 2003, pp [8] Gottlieb, D. and Orzag, S.., Numerical nalysis of Spectral Methods: theory and applications, Society for Industrial and pplied Mathematics, Bristol, England, [9] Powell, M.J.D., The Theory of Radial Basis Function pproximation, in dvances in Numerical nalysis, Vol. II, Light, W., ed., Oxford Science Publications, Oxford, pp [10] Divo, E., Kassab,.J., Mitteff, E., and Quintana, L. Parallel Domain Decomposition Technique for Meshless Methods pplications to Large- Scale Heat Transfer Problems, SME Paper: HT-FED [11] Harlow, F.H. and Welch, J.E., Physics of Fluids, Vol. 8, pp , [12] Divo, E. and Kassab,., Effective Domain Decomposition Meshless Formulation of Fully-Viscous Incompressible Fluid Flows, Proc. of BEM 27, WIT Press, Southampton, 2005, (in press.) [13] Turns, S., R., n Introduction to Combustion, 2nd Ed., McGraw-Hill, New York, 2000.
CHAPTER 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 informationA Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer
Copyright c 2008 ICCES ICCES, vol.6, no.1, pp.13-18 A Meshless Radial Basis Function Method for Fluid Flow with Heat Transfer K agamani Devi 1,D.W.Pepper 2 Summary Over the past few years, efforts have
More informationEfficient Large Scale Transient Heat Conduction Analysis Using A Parallelized Boundary Element Method
University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Efficient Large Scale Transient Heat Conduction Analysis Using A Parallelized Boundary Element Method 2006
More informationSolution Methods. Steady State Diffusion Equation. Lecture 04
Solution Methods Steady State Diffusion Equation Lecture 04 1 Solution methods Focus on finite volume method. Background of finite volume method. Discretization example. General solution method. Convergence.
More information5. FVM discretization and Solution Procedure
5. FVM discretization and Solution Procedure 1. The fluid domain is divided into a finite number of control volumes (cells of a computational grid). 2. Integral form of the conservation equations are discretized
More informationPoint interpolation method based on local residual formulation using radial basis functions
Structural Engineering and Mechanics, Vol. 14, No. 6 (2002) 713-732 713 Point interpolation method based on local residual formulation using radial basis functions G.R. Liu, L. Yan, J.G. Wang and Y.T.
More informationNumerical solution of Maxwell equations using local weak form meshless techniques
Journal of mathematics and computer science 13 2014), 168-185 Numerical solution of Maxwell equations using local weak form meshless techniques S. Sarabadan 1, M. Shahrezaee 1, J.A. Rad 2, K. Parand 2,*
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationPDE Solvers for Fluid Flow
PDE Solvers for Fluid Flow issues and algorithms for the Streaming Supercomputer Eran Guendelman February 5, 2002 Topics Equations for incompressible fluid flow 3 model PDEs: Hyperbolic, Elliptic, Parabolic
More informationSolving the One Dimensional Advection Diffusion Equation Using Mixed Discrete Least Squares Meshless Method
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 2015 Paper No. 292 Solving the One Dimensional Advection Diffusion Equation
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 informationComputation Time Assessment of a Galerkin Finite Volume Method (GFVM) for Solving Time Solid Mechanics Problems under Dynamic Loads
Proceedings of the International Conference on Civil, Structural and Transportation Engineering Ottawa, Ontario, Canada, May 4 5, 215 Paper o. 31 Computation Time Assessment of a Galerkin Finite Volume
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 informationIn Proc. of the V European Conf. on Computational Fluid Dynamics (ECFD), Preprint
V European Conference on Computational Fluid Dynamics ECCOMAS CFD 2010 J. C. F. Pereira and A. Sequeira (Eds) Lisbon, Portugal, 14 17 June 2010 THE HIGH ORDER FINITE ELEMENT METHOD FOR STEADY CONVECTION-DIFFUSION-REACTION
More informationCFD Analysis of High Temperature and High Velocity System: Plasma Torch
CFD Analysis of High Temperature and High Velocity System: Plasma Torch Abhishek Pratap Singh Bhadauria 1 1 Department of Mechanical Engineering, K. J. Somaiya College of Engineering, Mumbai, Maharashtra,
More informationComputation of Unsteady Flows With Moving Grids
Computation of Unsteady Flows With Moving Grids Milovan Perić CoMeT Continuum Mechanics Technologies GmbH milovan@continuummechanicstechnologies.de Unsteady Flows With Moving Boundaries, I Unsteady flows
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 informationInternational Journal of Scientific & Engineering Research, Volume 6, Issue 5, May ISSN
International Journal of Scientific & Engineering Research, Volume 6, Issue 5, May-2015 28 CFD BASED HEAT TRANSFER ANALYSIS OF SOLAR AIR HEATER DUCT PROVIDED WITH ARTIFICIAL ROUGHNESS Vivek Rao, Dr. Ajay
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 informationLocal Radial Basis Partition of Unity for Modeling Quantum Dot Semiconductor
Local Radial Basis Partition of Unity for Modeling Quantum Dot Semiconductor Abstract Solving the Schrodinger equation in quantum systems presents a challenging task for numerical methods due to the high
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 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 informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationIntroduction to Heat and Mass Transfer. Week 7
Introduction to Heat and Mass Transfer Week 7 Example Solution Technique Using either finite difference method or finite volume method, we end up with a set of simultaneous algebraic equations in terms
More informationIMPLEMENTATION OF A PARALLEL AMG SOLVER
IMPLEMENTATION OF A PARALLEL AMG SOLVER Tony Saad May 2005 http://tsaad.utsi.edu - tsaad@utsi.edu PLAN INTRODUCTION 2 min. MULTIGRID METHODS.. 3 min. PARALLEL IMPLEMENTATION PARTITIONING. 1 min. RENUMBERING...
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 informationALGEBRAIC FLUX CORRECTION FOR FINITE ELEMENT DISCRETIZATIONS OF COUPLED SYSTEMS
Int. Conf. on Computational Methods for Coupled Problems in Science and Engineering COUPLED PROBLEMS 2007 M. Papadrakakis, E. Oñate and B. Schrefler (Eds) c CIMNE, Barcelona, 2007 ALGEBRAIC FLUX CORRECTION
More informationCFD with OpenSource software
CFD with OpenSource software The Harmonic Balance Method in foam-extend Gregor Cvijetić Supervisor: Hrvoje Jasak Faculty of Mechanical Engineering and Naval Architecture University of Zagreb gregor.cvijetic@fsb.hr
More informationAn example of the Rvachev function method
arxiv:1603.00320v1 [physics.flu-dyn] 1 Mar 2016 An example of the Rvachev function method Alexander V. Proskurin Altai State University, Altai State Technical University, k210@list.ru Anatoly M. Sagalakov
More informationSOE3213/4: CFD Lecture 3
CFD { SOE323/4: CFD Lecture 3 @u x @t @u y @t @u z @t r:u = 0 () + r:(uu x ) = + r:(uu y ) = + r:(uu z ) = @x @y @z + r 2 u x (2) + r 2 u y (3) + r 2 u z (4) Transport equation form, with source @x Two
More informationUniversity of Illinois at Urbana-Champaign College of Engineering
University of Illinois at Urbana-Champaign College of Engineering CEE 570 Finite Element Methods (in Solid and Structural Mechanics) Spring Semester 03 Quiz # April 8, 03 Name: SOUTION ID#: PS.: A the
More informationAnalysis of Steady State Heat Conduction Problem Using EFGM
International Journal of Engineering and Management Research, Vol.-2, Issue-6, December 2012 ISSN No.: 2250-0758 Pages: 40-47 www.ijemr.net Analysis of Steady State Heat Conduction Problem Using EFGM Manpreet
More informationPerformance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds Number
Copyright c 2006 Tech Science Press CMES, vol.15, no.3, pp.137-146, 2006 Performance of Multiquadric Collocation Method in Solving Lid-driven Cavity Flow Problem with Low Reynolds umber S. Chantasiriwan
More informationThe Method of Fundamental Solutions with Eigenfunction Expansion Method for Nonhomogeneous Diffusion Equation
The Method of Fundamental Solutions with Eigenfunction Expansion Method for Nonhomogeneous Diffusion Equation D. L. Young, 1 C. W. Chen, 1 C. M. Fan, 1 C. C. Tsai 2 1 Department of Civil Engineering and
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 informationCFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) 31 (2011) 33-41 CFD STUDY OF MASS TRANSFER IN SPACER FILLED MEMBRANE MODULE Sharmina Hussain Department of Mathematics and Natural Science BRAC University,
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
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 informationFLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.
FLUID FLOW AND HEAT TRANSFER INVESTIGATION OF PERFORATED HEAT SINK UNDER MIXED CONVECTION 1 Mr. Shardul R Kulkarni, 2 Prof.S.Y.Bhosale 1 Research scholar, 2 Head of department & Asst professor Department
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
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 informationExcerpt from the Proceedings of the COMSOL Users Conference 2006 Boston
Using Comsol Multiphysics to Model Viscoelastic Fluid Flow Bruce A. Finlayson, Professor Emeritus Department of Chemical Engineering University of Washington, Seattle, WA 98195-1750 finlayson@cheme.washington.edu
More informationPressure-velocity correction method Finite Volume solution of Navier-Stokes equations Exercise: Finish solving the Navier Stokes equations
Today's Lecture 2D grid colocated arrangement staggered arrangement Exercise: Make a Fortran program which solves a system of linear equations using an iterative method SIMPLE algorithm Pressure-velocity
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationDraft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer
Draft Notes ME 608 Numerical Methods in Heat, Mass, and Momentum Transfer Instructor: Jayathi Y. Murthy School of Mechanical Engineering Purdue University Spring 00 c 1998 J.Y. Murthy and S.R. Mathur.
More informationA note on the numerical treatment of the k-epsilon turbulence model Λ
A note on the numerical treatment of the k-epsilon turbulence model Λ Adrián J. Lew y, Gustavo C. Buscaglia z and Pablo M. Carrica Centro Atómico Bariloche and Instituto Balseiro, 8400 Bariloche, Argentina.
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 informationINTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING---
ILASS-Europe 2002 Zaragoza 9 11 September 2002 INTERNAL FLOW IN A Y-JET ATOMISER ---NUMERICAL MODELLING--- Z. Tapia, A. Chávez e-mail: ztapia@imp.mx Instituto Mexicano del Petróleo Blvd. Adolfo Ruiz Cortines
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 informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationCFD in COMSOL Multiphysics
CFD in COMSOL Multiphysics Mats Nigam Copyright 2016 COMSOL. Any of the images, text, and equations here may be copied and modified for your own internal use. All trademarks are the property of their respective
More informationFinite volume method for CFD
Finite volume method for CFD Indo-German Winter Academy-2007 Ankit Khandelwal B-tech III year, Civil Engineering IIT Roorkee Course #2 (Numerical methods and simulation of engineering Problems) Mentor:
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationCompact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems
Copyright 2011 Tech Science Press SL, vol.6, no.2, pp.93-107, 2011 Compact Local Stencils Employed With Integrated RBFs For Fourth-Order Differential Problems T.-T. Hoang-Trieu 1, N. Mai-Duy 1 and T. Tran-Cong
More informationApplication of meshless EFG method in fluid flow problems
Sādhan ā Vol. 29, Part 3, June 2004, pp. 285 296. Printed in India Application of meshless EFG method in fluid flow problems I V SINGH Mechanical Engineering Group, Birla Institute of Technology and Science,
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationPartial Differential Equations
Partial Differential Equations Introduction Deng Li Discretization Methods Chunfang Chen, Danny Thorne, Adam Zornes CS521 Feb.,7, 2006 What do You Stand For? A PDE is a Partial Differential Equation This
More informationEVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS
EVALUATION OF FOUR TURBULENCE MODELS IN THE INTERACTION OF MULTI BURNERS SWIRLING FLOWS A Aroussi, S Kucukgokoglan, S.J.Pickering, M.Menacer School of Mechanical, Materials, Manufacturing Engineering and
More informationTermination criteria for inexact fixed point methods
Termination criteria for inexact fixed point methods Philipp Birken 1 October 1, 2013 1 Institute of Mathematics, University of Kassel, Heinrich-Plett-Str. 40, D-34132 Kassel, Germany Department of Mathematics/Computer
More informationdifferential problems
A numerical study of 2D integrated RBFNs incorporating Cartesian grids for solving 2D elliptic differential problems N. Mai-Duy and T. Tran-Cong Computational Engineering and Science Research Centre Faculty
More informationThe Effect of Undulated Wall on Natural Convection Heat Transfer in Vertical Channels
Proceedings of the 3 rd International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT 16) Ottawa, Canada May 2 3, 2016 Paper No. 123 The Effect of Undulated Wall on Natural Convection Heat Transfer
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationNumerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led Application by using CFD
GRD Journals- Global Research and Development Journal for Engineering Volume 1 Issue 8 July 2016 ISSN: 2455-5703 Numerical Investigation of Convective Heat Transfer in Pin Fin Type Heat Sink used for Led
More informationConstruction of a New Domain Decomposition Method for the Stokes Equations
Construction of a New Domain Decomposition Method for the Stokes Equations Frédéric Nataf 1 and Gerd Rapin 2 1 CMAP, CNRS; UMR7641, Ecole Polytechnique, 91128 Palaiseau Cedex, France 2 Math. Dep., NAM,
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 informationNumerical Modeling of Inclined Negatively Buoyant Jets
Numerical Modeling of Inclined Negatively Buoyant Jets Presentation by: Hossein Kheirkhah Graduate Student in Civil Engineering Dep. of Civil Engineering University of Ottawa CANADA ICDEMOS April2014 Outline
More informationTowards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below
Copyright 2008 Tech Science Press FDMP, vol.4, no.4, pp.263-269, 2008 Towards a Numerical Benchmark for 3D Low Mach Number Mixed Flows in a Rectangular Channel Heated from Below G. Accary 1, S. Meradji
More informationMasters in Mechanical Engineering. Problems of incompressible viscous flow. 2µ dx y(y h)+ U h y 0 < y < h,
Masters in Mechanical Engineering Problems of incompressible viscous flow 1. Consider the laminar Couette flow between two infinite flat plates (lower plate (y = 0) with no velocity and top plate (y =
More informationCOMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES OF DIFFERENT ASPECT RATIOS
HEFAT214 1 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 214 Orlando, Florida COMPUTATIONAL ANALYSIS OF LAMINAR FORCED CONVECTION IN RECTANGULAR ENCLOSURES
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More informationLeast-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations
Least-Squares Spectral Collocation with the Overlapping Schwarz Method for the Incompressible Navier Stokes Equations by Wilhelm Heinrichs Universität Duisburg Essen, Ingenieurmathematik Universitätsstr.
More informationSimulation of sympathetic detonation by a CIP Eulerian code
Computational Ballistics II 107 Simulation of sympathetic detonation by a CIP Eulerian code S. Kubota 1, Z. Liu 2, T. Saburi 1, Y. Ogata 1 & M. Yoshida 1 1 Research Center for Explosion Safety, National
More informationApplications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin
Applications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin Faculty of Mechanical Engineering and Naval Architecture, Zagreb 10000, Croatia, gregor.cvijetic@gmail.com The Harmonic
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationNumerical investigation of swirl flow inside a supersonic nozzle
Advances in Fluid Mechanics IX 131 Numerical investigation of swirl flow inside a supersonic nozzle E. Eslamian, H. Shirvani & A. Shirvani Faculty of Science and Technology, Anglia Ruskin University, UK
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 informationIterative Solvers in the Finite Element Solution of Transient Heat Conduction
Iterative Solvers in the Finite Element Solution of Transient Heat Conduction Mile R. Vuji~i} PhD student Steve G.R. Brown Senior Lecturer Materials Research Centre School of Engineering University of
More informationMulti-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media
Multi-Point Special Boundary-Value Problems and Applications to Fluid Flow Through Porous Media Mohamed A. Hajji Abstract In this paper we propose a numerical scheme based on finite differences for the
More informationCalculating equation coefficients
Fluid flow Calculating equation coefficients Construction Conservation Equation Surface Conservation Equation Fluid Conservation Equation needs flow estimation needs radiation and convection estimation
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 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 informationLinear Solvers. Andrew Hazel
Linear Solvers Andrew Hazel Introduction Thus far we have talked about the formulation and discretisation of physical problems...... and stopped when we got to a discrete linear system of equations. Introduction
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationA numerical study of heat transfer and fluid flow over an in-line tube bank
Fluid Structure Interaction VI 295 A numerical study of heat transfer and fluid flow over an in-line tube bank Z. S. Abdel-Rehim Mechanical Engineering Department, National Research Center, Egypt Abstract
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 informationNumerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method To cite this article: R
More informationCFD study for cross flow heat exchanger with integral finned tube
International Journal of Scientific and Research Publications, Volume 6, Issue 6, June 2016 668 CFD study for cross flow heat exchanger with integral finned tube Zena K. Kadhim *, Muna S. Kassim **, Adel
More informationAnalysis of the Cooling Design in Electrical Transformer
Analysis of the Cooling Design in Electrical Transformer Joel de Almeida Mendes E-mail: joeldealmeidamendes@hotmail.com Abstract This work presents the application of a CFD code Fluent to simulate the
More informationNumerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions
Advanced Computational Methods in Heat Transfer X 25 Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions F. Selimovic & B. Sundén
More informationAnalysis of Heat Transfer in Pipe with Twisted Tape Inserts
Proceedings of the 2 nd International Conference on Fluid Flow, Heat and Mass Transfer Ottawa, Ontario, Canada, April 30 May 1, 2015 Paper No. 143 Analysis of Heat Transfer in Pipe with Twisted Tape Inserts
More information20. A Dual-Primal FETI Method for solving Stokes/Navier-Stokes Equations
Fourteenth International Conference on Domain Decomposition Methods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 23 DDM.org 2. A Dual-Primal FEI Method for solving Stokes/Navier-Stokes
More informationENERGY PERFORMANCE IMPROVEMENT, FLOW BEHAVIOR AND HEAT TRANSFER INVESTIGATION IN A CIRCULAR TUBE WITH V-DOWNSTREAM DISCRETE BAFFLES
Journal of Mathematics and Statistics 9 (4): 339-348, 2013 ISSN: 1549-3644 2013 doi:10.3844/jmssp.2013.339.348 Published Online 9 (4) 2013 (http://www.thescipub.com/jmss.toc) ENERGY PERFORMANCE IMPROVEMENT,
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationIncompressible Flow Solver by Means of Pseudo-Compressibility Method
Incompressible Flow Solver by Means of Pseudo-Compressibility Method Árpád Veress Department of Aircraft and Ships Budapest University of Technology and Economics H Budapest, Sztoczek u. 6. J. ép., Hungary
More informationPerformance characteristics of turbo blower in a refuse collecting system according to operation conditions
Journal of Mechanical Science and Technology 22 (2008) 1896~1901 Journal of Mechanical Science and Technology www.springerlink.com/content/1738-494x DOI 10.1007/s12206-008-0729-6 Performance characteristics
More information3D Numerical Simulation of Supercritical Flow in Bends of Channel
3D Numerical Simulation of Supercritical Flow in Bends of Channel Masoud. Montazeri-Namin, Reyhaneh-Sadat. Ghazanfari-Hashemi, and Mahnaz. Ghaeini- Hessaroeyeh Abstract An attempt has been made to simulate
More informationA PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM
A PRACTICALLY UNCONDITIONALLY GRADIENT STABLE SCHEME FOR THE N-COMPONENT CAHN HILLIARD SYSTEM Hyun Geun LEE 1, Jeong-Whan CHOI 1 and Junseok KIM 1 1) Department of Mathematics, Korea University, Seoul
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 information