The gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last
|
|
- Lewis Fowler
- 5 years ago
- Views:
Transcription
1 Parallel Implementation of Gas-Kinetic BGK Scheme on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer Middle East Technical University, Ankara, 653, Turkey The gas-kinetic BGK scheme on unstructured hybrid grids is presented. In order to accurately resolve the boundary layers in wall bounded viscous flow solutions, quadrilateral grid cells are employed in the boundary layer regions normal to solid surfaces while the rest of the domain is discretized by triangular cells. The computation time, which is a significant deficiency of the gas-kinetic schemes, is improved by performing computations in parallel. The parallel algorithm for hybrid grids is based on the domain decomposition using METIS, a graph partitioning software. Several numerical test cases are presented to show the accuracy and robustness of the proposed approach. I. Introduction The gas-kinetic methods have become popular for the simulation of compressible fluid flows in the last decade. The most promising ones are the Equilibrium Flux Method, the Kinetic Flux Vector Splitting 2, 3 and the Gas-Kinetic BGK method. 4 Due to the inclusion of intermolecular collisions with BGK simplification, the gas-kinetic BGK scheme gives a more complete and realistic description of the flow and has been well studied. 5 2 Although the first studies about the gas-kinetic BGK method were based on finite volume formulation on structured grids, the obvious advantages of unstructured grids over structured counterparts for complex configurations have led researcher to adopt the scheme on unstructured meshes. 3 5 Recently, in the previous work of the present authors, the gas-kinetic flux vector splitting and the gas-kinetic BGK schemes on unstructured grids were given. 6 The solutions were obtained in parallel in order to reduce the computational inefficiency common in gas-kinetic schemes. Several numerical test cases were presented and it was shown that the gas-kinetic schemes are efficient in parallel computations and do not experience numerical instabilities as their classical counterparts. However, in general, unstructured grids are not so suitable for the simulation of viscous flows especially for the well resolution of the boundary layers. In the present work, a high-order gas-kinetic scheme for the Navier-Stokes equations on unstructured hybrid grids is given. The finite volume formulations are combined with explicit Runge-Kutta time-stepping scheme. The parallel algorithm and the domain decomposition for unstructured hybrid meshes are explained and several numerical test cases are presented to show the accuracy and robustness of the present approach. II. Gas-Kinetic Theory In gas-kinetic theory, gases are comprised of small particles and each particle has a mass and velocity. Due to the existence of large number of particles in a small volume at standard conditions (e.g. at STP, 9 air molecules in cm 3 volume), it is very difficult to trace individual particle motion. Instead, a particle distribution function is defined to describe the probability of particles to be located in a certain velocity interval Senior Research Engineer, Aerodynamics Division, AIAA Student Member. Professor, Department of Aerospace Engineering, AIAA Member. of 3
2 f(x i, t, u i ) () Here x i = (x, y, z) is the position, t is the time and u i = (u, v, w) are the particle velocities. The macroscopic properties of the gas can be obtained as the moments of the distribution function. For example, the gas density can be written as ρ = m n i (2) i where m is the particle mass, n i is the number density. Since, by definition, distribution function is the particle density in phase space, it is concluded that m n i = f(x i, t, u i ), (3) ρ = f dudvdw. (4) The time evolution of the distribution function is governed by the Boltzmann equation 7 f t + u i f xi + a i f ui = Q(f, f). (5) Here a i shows the external force on the particle in the ith direction and Q(f, f) is the collision operator. When the collision operator is equal to zero, collisionless Boltzmann equation is obtained and the solution of this equation gives the Maxwellian (equilibrium) distribution function g = ρ ( λ π ) N+3 2 exp { λ [(u i U i ) 2 + ξ i 2 ]} (6) where ξ i = (ξ, ξ 2,..., ξ N ) are the particle internal velocities, N is the internal degrees of freedom, U i = (U, V, W ) are the macroscopic velocities of the gas and λ is a function of temperature given by R being the gas constant. λ = 2 R T (7) III. Numerical Implementation The proposed numerical method is based on the cell-centered finite volume formulation. The numerical fluxes are calculated at the edges of the cells and the conservative variables at the cell centers are updated using a third-order explicit Runge-Kutta time-stepping scheme. In the following, the numerical implementation of the proposed approach is presented in terms of the initial reconstruction of the conservative flow variables and gas evolution stages. A. Initial Reconstruction The initial reconstruction of conservative variables is based on the method proposed by Frink. 8 In this method, a universal expression for the Taylor series expansion within a triangular cell is derived, which requires the knowledge of the state at the nodes. The Taylor series expansion of the conservative variables in the neighborhood of each cell center (x cc, y cc ) is Q(x, y) = Q(x cc, y cc ) + Q cc r + O( r 2 ) (8) where Q = ρ ρu ρv ρe. (9) 2 of 3
3 Figure. Arbitrary triangular cell. This formulation requires the gradients of the flow variables at the cell centers. The evaluation of the gradient vectors is based on the Green s theorem: Q da = Q ˆn dl () Ω The average gradient vector at the cell center, Q cc, is then given by the closed boundary integral over the cell Q cc = Q ˆn dl () A Ω where the unit normal ˆn is assumed to point outward from the domain. Finally, the flow variables at the middle of the cell edges, Q ci, (Fig. ) is given by Frink 8 Ω Ω Q ci = Q cc + 3 [ 2 (Q 2 + Q 3 ) Q ] (2) where the nodal values are calculated from the surrounding cell center values. Mitchell 9 stated that this expression is an average of one central difference term and two upwind terms and recovers a family of schemes common with structured grid schemes. Rewriting Eq. (2) Q ci = 3 [ 2 (Q 2 + Q 3 ) + 2( 3 2 Q cc 2 Q )], (3) Q ci = w[c.d.] + ( w)[upwind] (4) where w = is fully central differencing, w = is fully upwinding. Note that w = /3 recovers Eq. (2) which is incidently suggested by Mitchell as the optimum value for a variety of test functions and mesh distributions. 9 The reconstruction algorithm given for triangular cells may similarly be applied to quadrilateral cells. For an arbitrary quadrilateral cell shown in Fig. 2, the triangular cell in red is used for the reconstruction of the flow variables on cell edge 2-3: Q ci = Q cc + 3 [ 2 (Q 2 + Q 3 ) Q cc ] (5) For the remaining edges, the corresponding triangular cells within the quadrilateral cell are used for reconstruction. It should be noted that since the reconstruction procedure involves the simple averaging of the flow variables, the formal accuracy of the method may be less than second-order. B. Gas Evolution The gas evolution is based on the gas-kinetic BGK scheme. The gas-kinetic BGK scheme is based on the Boltzmann BGK equation where the collision operator (see Eq. (5)) is replaced by the Bhatnagar-Gross- Krook model. 2 The Boltzmann BGK equation in two-dimensions can be written as f t + u f x + v f y = g f τ (6) 3 of 3
4 Figure 2. Arbitrary quadrilateral cell. where f is the particle distribution function, g is the equilibrium state approached by f over particle collision time τ, u and v are the particle velocities in x- and y-directions, respectively. The equilibrium state is usually assumed to be a Maxwellian g = ρ ( λ π ) K+2 2 exp { λ [(u U) 2 + (v V ) 2 + ξ 2 ]} (7) where ρ is the density, U and V are the macroscopic velocities in x- and y-directions and K is the dimension of the internal velocities. The general solution of the particle distribution function f at the cell edge ci and time t is f(s ci, t, u, v, ξ) = τ t g(s, t, u, v, ξ) exp [ (t t )/τ] dt + exp ( t/τ) f (s ci ut vt). (8) Here s = s ci u (t t ) v (t t ) is the particle trajectory of a particle and f is the initial gas distribution function at the beginning of each time step. The relation between mass ρ, momentum ρu, ρv and total energy ρe densities with the distribution function is defined based on Eq. (4) ρ ρu ρv ρe where ψ represents the vector of moments for the distribution function u ψ = v 2 (u2 + v 2 + ξ 2 ) = f ψ dξ (9) and dξ is the volume element in phase space. Since mass, momentum and energy are conserved during particle collisions, f and g must satisfy the conservation constraint of (f g) ψ dξ =, (2) at all (x, y) and t. The high-order gas-kinetic BGK scheme on unstructured hybrid meshes is implemented as follows: Consider the control volumes (triangular or quadrilateral), their neighbors and local coordinates given in Fig. 3 where dots represent cell centers, the triangles or quadrilaterals in red show the control volumes, L and R the left and right states, x n and x t the local coordinate system normal and tangent to the cell edge ci, respectively, x and y the global coordinate system. In cell-centered, finite volume gas-kinetic BGK scheme on unstructured hybrid meshes, the velocity components U n and U t in local (normal and tangential) coordinate system are found from global coordinate system counterparts U and V for both mesh types. (2) 4 of 3
5 Figure 3. Control volumes, their neighbors and local coordinate system. In the present work, the initial gas distribution function f and the equilibrium state g are assumed to be { g f = L [ + a L (x n x nci )], if x n x nci (22) g R [ + a R (x n x nci )], if x n x nci and g = { g [ + a L (x n x nci ) + A t], if x n x nci (23) g [ + a R (x n x nci ) + A t], if x n x nci where g L, g R and g are the equilibrium distribution functions to the left, right and in the middle of the cell edge ci, respectively, of the form g = ρ ( λ π ) K+2 2 exp { λ [(u n U n) 2 + (u t U t ) 2 + ξ 2 ]}, (24) and a L, a R, a L, and a R are the spatial slopes and A is the time slope, given by a L = a L + a L 2 u n + a L 3 u t + 2 al 4 (u 2 n + u 2 t + ξ 2 ) a R = a R + a R 2 u n + a R 3 u t + 2 ar 4 (u 2 n + u 2 t + ξ 2 ) a L = a L + a L 2 u n + a L 3 u t + 2 al 4 (u 2 n + u 2 t + ξ 2 ) a R = a R + a R 2 u n + a R 3 u t + 2 ar 4 (u 2 n + u 2 t + ξ 2 ) A = A + A 2 u n + A 3 u t + 2 A 4 (u 2 n + u 2 t + ξ 2 ) (25) Here u n and u t are the particle velocities normal and tangent to the cell edge, respectively. With the reconstructed values obtained in III.A, the equilibrium distribution functions to the left and right of the cell edge as well as their slopes can be determined g L ψ ci dξ = ρ L ci ρu L n ci ρu L t ci ρe L ci, g R ψ ci dξ = ρ R ci ρu R n ci ρu R t ci ρe R ci (26) 5 of 3
6 g L a L ψ ci dξ = (ρ L ci ρl cc)/ s (ρu L n ci ρu L n cc )/ s (ρu L t ci ρu L t cc )/ s (ρe L ci ρel cc)/ s, g R a R ψ ci dξ = (ρ R cc ρ R ci )/ s (ρu R n cc ρu R n ci )/ s (ρu R t cc ρu R t ci )/ s (ρe R cc ρe R ci )/ s (27) where the subscript ci corresponds to the reconstructed values at the cell edge, the subscript cc refers to the cell center values, s is the distance from the cell edge to the cell center and ψ ci stands for the vector of moments at the cell edge u ψ ci = n (28) u t 2 (u2 n + u 2 t + ξ 2 ) Considering Eq. (24), all the parameters in g L and g R can be uniquely determined in Eq. (26). Once g L and g R are obtained, the slopes a L and a R can be computed from Eq. (27). After determining f, the equilibrium state g can be found using the compatibility condition g ψ ci dξ = g L ψ ci dξ + g R ψ ci dξ (29) u> from which the values of ρ, U n, U t and λ in g are determined. Then the slopes a L and a R can be obtained through the relations (ρ ρ L cc)/ s (ρ R cc ρ g a L (ρu ψ ci dξ = n ρun L )/ s cc )/ s (ρu t ρut L cc )/ s, g a R (ρu ψ ci dξ = n cc ρu n )/ s (ρu R t cc ρu t )/ s (3) (ρe ρecc)/ s L (ρecc R ρe )/ s Up to this point, the equilibrium states to the left, to the right and in the middle of the cell edge as well as the corresponding spatial slopes are determined. The only unknown term is the time slope term A. Since both f and g contain the time slope A, the conservation constraint (Eq. (2)) at the cell edge ci can be applied and integrated over the time step t t (f g) ψ ci dξ dt =, (3) from which A can be calculated. Substituting Eqs. (22) and (23) into Eq. (8), the final gas distribution function at the cell edge ci is expressed as u< f(x nci, t, u, v, ξ) = [ exp ( t/τ)] g +{τ [ + exp ( t/τ)] + t exp ( t/τ)} [a L H(u) + a R ( H(u))] u g +τ [t/τ + exp ( t/τ)] A g +exp ( t/τ) [( u t a L ) H(u) g L + ( u t a R ) ( H(u)) g L ] (32) where H(u) is the Heaviside function H(u) = {, if u <, if u >. (33) The local numerical fluxes for the mass, momentum and total energy across the cell edge ci can then be computed as F ρ F ρun F ρut F ρe = t u n f xnci ψ ci dξ dt (34) The fluxes at the cell faces are first computed and the total flux for both triangular and quadrilateral cells is then obtained by adding the flux contributions. 6 of 3
7 IV. Parallel Processing Parallel processing is based on domain decomposition. The unstructured hybrid grid is partitioned using METIS software package. METIS needs the graph file for the unstructured hybrid mesh, which is actually the neighbor connectivity of the cells. A sample unstructured hybrid mesh and the corresponding graph is given in Fig. 4. The first line in the graph represents the number of cells, number of nodes, weighting option and number of weights, respectively. The remaining lines show the number of faces and the indices of the neighbors of the corresponding cell. Partitioning of the graph is performed using kmetis program. During the partitioning, each cell is weighted by its number of edges so that each partition has about the same number of total edges to improve the load balancing in parallel computations. Figure 4. Sample computational mesh and the corresponding graph. Parallel Virtual Machine (PVM) message-passing library routines are employed in a master-worker algorithm. The master process performs all the input-output, starts up pvm, spawns worker processes and sends the initial data to the workers. The worker processes first receive the initial data, apply the interface and the flow boundary conditions, and solve the flow field within the partition. The flow variables at the interface boundaries are exchanged among the neighboring partitions at each time step for the implementation of inter-partition boundary conditions. V. Results and Discussions In order to show the accuracy and robustness of the present approach, the gas-kinetic BGK scheme is applied to a laminar flow over a flat plate and a transonic viscous flow over an airfoil. The results are compared to the analytical solutions and available experimental data. A. Case : Laminar Flow over a Flat Plate A flat plate at zero angle of attack is selected as a first test case for which the analytical Blasius solution is available. The freestream conditions and the computational mesh for this case are given in Table and Fig. 5, respectively. Table. Freestream conditions for Case. Mach Number Angle of Attack, deg Pressure, psia Temperature, R The above conditions correspond to a Reynolds number of 2. The flat plate is placed between and and the boundary layer region is meshed with quadrilateral cells while triangular cells are used outside the boundary layer. Fig. 6 shows the velocity vectors inside the boundary layer on the flat plate and the 7 of 3
8 .8.6 y x Figure 5. Computational mesh for Case (zoomed). y.6 Mach x Figure 6. Mach number contours and velocity vectors (zoomed). 8 of 3
9 8 7 6 Exact x =.25 x =.4 x =.55 5 η u/u Figure 7. u-velocity profile in the boundary layer. comparison of u-velocity profile at three different locations with the Blasius solution is given in Fig. 7. The results compare quite well with the exact Blasius solution. It should be noted that the gas-kinetic BGK method produces Navier-Stokes solution in smooth flow regions as expected. Moreover, the use of quadrilateral cells near the wall region leads to the better resolution of the boundary layer profile while the use of triangular cells outside reduce the overall number of cells. The laminar flow solution given above verifies the accuracy and robustness of the present approach. B. Case 2: Transonic Flow over an Airfoil RAE 2822 transonic airfoil, which has available experimental data, is selected as a second test case. The freestream conditions for this case are given in Table 2. The computational mesh and the partitions generated using the METIS software is presented in Fig. 8. Table 2. Freestream condition for Case 2. Mach Number Angle of Attack, deg Pressure, psia Temperature, R Fig. 9 shows the Mach number contours obtained from the gas-kinetic BGK scheme on the unstructured hybrid grid. As shown, the flow solution is smooth and the shock is captured well. The comparison of the surface pressure distribution with the experimental data is given in Fig.. As seen, the shock location and the strength compare quite well and the overall pressure distribution is in excellent agreement. In Fig., the present viscous solution is compared to the inviscid Euler solution. The viscous effects are clearly seen in the prediction of shock location and the pressure distribution downstream of the shock. The parallel computations are performed on an Itanium cluster running on Linux. Dual Itanium2 processors operate at.3ghz with 3MB L2 cache and 2GB of memory for each. The parallel efficiency of 9 of 3
10 Figure 8. Computational mesh and partitions for Case 2 (zoomed)..2 y Mach x Figure 9. Mach number contours. of 3
11 - Pressure Coefficient Experiment Gas-Kinetic BGK x Figure. Pressure distribution on the airfoil surface (every other three data point is displayed). - Pressure Coefficient Inviscid Solution Viscous Solution x Figure. Viscous versus inviscid solution (every other three data point is displayed). of 3
12 Table 3. Parallel efficiency of computations. Number of Nodes Computational Efficiency, sec/iter the computations is given in Table 3 and Fig. 2. It is observed that the high parallel efficiency of the gaskinetic BGK scheme is maintained as the number of processors is increased, which is attributed to the high computing to communication ratio. Although the gas-kinetic BGK schemes are computationally expensive and require longer computational time than classical schemes in serial computations, 6 it may not be an issue in parallel computations Ideal Gas-Kinetic BGK Speed Up Number of Processors Figure 2. Computational speed up. VI. Conclusion In this paper, the gas-kinetic BGK scheme on unstructured hybrid grids are presented. The high-order finite volume formulations are given. The solutions are obtained in parallel and the results are compared to the available analytical/experimental data. The unstructured hybrid grids used in the present study removes the difficulties faced in viscous flow computations with triangular grids. The viscous flow solutions with the gas-kinetic BGK scheme on unstructured hybrid grids agree well with analytical/experimental data. In addition, the parallel computations improve the computational efficiency of the gas-kinetic BGK schemes significantly. 2 of 3
13 Acknowledgments The authors acknowledge the support of the Defense Industries Research and Development Institute (TUBITAK-SAGE) under project SAM. References Pullin, D. I., Direct Simulation Methods for Compressible Inviscid Ideal Gas Flow, J. Comp. Phys., Vol. 34, 98, pp Mandal, J. C., and Deshpande, S. M., Kinetic Flux Vector Splitting for Euler Equations, Comp. Fluids, No. 23-2, 994, p Chou, S. Y., and Baganoff, D., Kinetic Flux Vector Splitting for the Navier-Stokes Equations, J. Comp. Phys., Vol. 3, 997, pp Prendergast, K. H., and Xu, K., Numerical Hydrodynamics from Gas-Kinetic Theory, J. Comp. Phys., Vol. 9, 993, pp Xu, K., Martinelli, L., and Jameson, A., Gas-Kinetic Finite Volume Methods, Flux Vector Splitting and Artificial Diffusion, J. Comp. Phys., Vol. 2, 995, pp Xu, K., and Jameson, A., Gas-Kinetic Relaxation (BGK-Type) Schemes for the Compressible Euler Equations, AIAA Paper , Xu, K., BGK-Based Scheme for Multicomponent Flow Calculations, J. Comp. Phys., Vol. 34, 997, pp Xu, K., A Gas-Kinetic Scheme for the Euler Equations with Heat Transfer, SIAM J. Sci. Comp., Vol. 2-4, 997, pp Xu, K., and Hu, J., Projection Dynamics in Godunov-Type Schemes, J. Comp. Phys., Vol. 42, 998, pp Lian, Y. S., Xu, K., A Gas-Kinetic Scheme for Multimaterial Flows and Its Application in Chemical Reactions, J. Comp. Phys., Vol. 63, 2, pp Chae D., Kim C., and Rho O., Development of an Improved Gas-Kinetic BGK Scheme for Inviscid and Viscous Flows, J. Comp. Phys., Vol. 58, 2, pp Xu, K., A Gas-Kinetic BGK Scheme for the Navier-Stokes Equations and Its Connection with Artificial Dissipation and Godunov Method, J. Comp. Phys., Vol. 7, 2, pp Kim, C., and Jameson, A., A Robust and Accurate LED-BGK Solver on Unstructured Adaptive Meshes, J. Comp. Phys., Vol. 43, 998, pp May, G., and Jameson, A., Improved Gaskinetic Multigrid Method for Three-Dimensional Computation of Viscous Flows, AIAA Paper 25-56, May, G., Srinivasan, B., and Jameson, A., Three Dimensional Flow Calculations on Arbitrary Meshes Using a Gas- Kinetic BGK Finite-Volume Method, AIAA Paper , Ilgaz, M., and Tuncer, I., H., Parallel Implementation of Gas-Kinetic Schemes for 2-D Flows on Unstructured Grids, 3rd Ankara International Aerospace Conference [CD-ROM], AIAC-25-8, Ankara, Turkey, Cercignani, C., The Boltzmann Equation and Its Applications, Springer-Verlag, Frink, N., T., Three-Dimensional Upwind Scheme for Solving the Euler Equations on Unstructured Tetrahedral Grids, Ph.D. Thesis, Virginia Polytechnic Institute and State University, Mitchell, C. R., Improved Reconstruction Scheme for the Navier-Stokes Equations on Unstructured Meshes, AIAA Paper , Bhatnagar, P. L., Gross, E. P., and Krook, M., A Model for Collision Processes in Gases I: Small Amplitude Processes in Charged and Neutral One-Component Systems, Phys. Rev., Vol. 94, pp. 5, of 3
The use of gas-kinetic schemes for the simulation of compressible flows become widespread in the two last
A Gas-Kinetic BGK Scheme for Parallel Solution of 3-D Viscous Flows on Unstructured Hybrid Grids Murat Ilgaz Defense Industries Research and Development Institute, Ankara, 626, Turkey and Ismail H. Tuncer
More informationDebojyoti Ghosh. Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering
Debojyoti Ghosh Adviser: Dr. James Baeder Alfred Gessow Rotorcraft Center Department of Aerospace Engineering To study the Dynamic Stalling of rotor blade cross-sections Unsteady Aerodynamics: Time varying
More informationModeling External Compressible Flow
Tutorial 3. Modeling External Compressible Flow Introduction The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a nonzero angle of attack. You will use the Spalart-Allmaras
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationThree dimensional meshless point generation technique for complex geometry
Three dimensional meshless point generation technique for complex geometry *Jae-Sang Rhee 1), Jinyoung Huh 2), Kyu Hong Kim 3), Suk Young Jung 4) 1),2) Department of Mechanical & Aerospace Engineering,
More informationAn Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid
An Investigation of Directional-Coarsening And Line-Implicit Smoothing Applied to Agglomeration Multigrid J. V. Lassaline Ryerson University 35 Victoria St, Toronto, ON, M5B 2K3, Canada D. W. Zingg University
More informationOn the high order FV schemes for compressible flows
Applied and Computational Mechanics 1 (2007) 453-460 On the high order FV schemes for compressible flows J. Fürst a, a Faculty of Mechanical Engineering, CTU in Prague, Karlovo nám. 13, 121 35 Praha, Czech
More informationAirfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition
Airfoil Design Optimization Using Reduced Order Models Based on Proper Orthogonal Decomposition.5.5.5.5.5.5.5..5.95.9.85.8.75.7 Patrick A. LeGresley and Juan J. Alonso Dept. of Aeronautics & Astronautics
More informationThe Development of a Navier-Stokes Flow Solver with Preconditioning Method on Unstructured Grids
Proceedings of the International MultiConference of Engineers and Computer Scientists 213 Vol II, IMECS 213, March 13-15, 213, Hong Kong The Development of a Navier-Stokes Flow Solver with Preconditioning
More informationAxisymmetric Viscous Flow Modeling for Meridional Flow Calculation in Aerodynamic Design of Half-Ducted Blade Rows
Memoirs of the Faculty of Engineering, Kyushu University, Vol.67, No.4, December 2007 Axisymmetric Viscous Flow Modeling for Meridional Flow alculation in Aerodynamic Design of Half-Ducted Blade Rows by
More informationDiscontinuous Galerkin Spectral Element Approximations for CFD
Discontinuous Galerkin Spectral Element Approimations for CFD D.A. Kopriva Florida State Universit Tallahassee, FL 3236 G.B. Jacobs San Diego State Universit San Diego, CA 92182 September 3, 211 1 Code
More informationLATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS
14 th European Conference on Mixing Warszawa, 10-13 September 2012 LATTICE-BOLTZMANN METHOD FOR THE SIMULATION OF LAMINAR MIXERS Felix Muggli a, Laurent Chatagny a, Jonas Lätt b a Sulzer Markets & Technology
More informationAil implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations
Ail implicit finite volume nodal point scheme for the solution of two-dimensional compressible Navier-Stokes equations Vimala Dutta Computational and Theoretical Fluid Dynamics Division National Aerospace
More informationGrid. Apr 09, 1998 FLUENT 5.0 (2d, segregated, lam) Grid. Jul 31, 1998 FLUENT 5.0 (2d, segregated, lam)
Tutorial 2. Around an Airfoil Transonic Turbulent Flow Introduction: The purpose of this tutorial is to compute the turbulent flow past a transonic airfoil at a non-zero angle of attack. You will use the
More informationNumerical Methods for PDEs. SSC Workgroup Meetings Juan J. Alonso October 8, SSC Working Group Meetings, JJA 1
Numerical Methods for PDEs SSC Workgroup Meetings Juan J. Alonso October 8, 2001 SSC Working Group Meetings, JJA 1 Overview These notes are meant to be an overview of the various memory access patterns
More informationComputational Fluid Dynamics for Engineers
Tuncer Cebeci Jian P. Shao Fassi Kafyeke Eric Laurendeau Computational Fluid Dynamics for Engineers From Panel to Navier-Stokes Methods with Computer Programs With 152 Figures, 19 Tables, 84 Problems and
More informationA Fourth-Order Gas-Kinetic CPR Method for the Navier-Stokes Equations on Unstructured Meshes
Tenth International Conference on Computational Fluid Dnamics (ICCFD), Barcelona,Spain, Jul 9-, 8 ICCFD-8- A Fourth-Order Gas-Kinetic CPR Method for the Navier-Stokes Equations on Unstructured Meshes Chao
More informationA higher-order finite volume method with collocated grid arrangement for incompressible flows
Computational Methods and Experimental Measurements XVII 109 A higher-order finite volume method with collocated grid arrangement for incompressible flows L. Ramirez 1, X. Nogueira 1, S. Khelladi 2, J.
More informationHigh-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method
High-Order Navier-Stokes Simulations using a Sparse Line-Based Discontinuous Galerkin Method Per-Olof Persson University of California, Berkeley, Berkeley, CA 9472-384, U.S.A. We study some of the properties
More informationComparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils
Comparisons of Compressible and Incompressible Solvers: Flat Plate Boundary Layer and NACA airfoils Moritz Kompenhans 1, Esteban Ferrer 2, Gonzalo Rubio, Eusebio Valero E.T.S.I.A. (School of Aeronautics)
More informationEFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES
EFFICIENT SOLUTION ALGORITHMS FOR HIGH-ACCURACY CENTRAL DIFFERENCE CFD SCHEMES B. Treidler, J.A. Ekaterineris and R.E. Childs Nielsen Engineering & Research, Inc. Mountain View, CA, 94043 Abstract Preliminary
More informationSHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY
SHOCK WAVES IN A CHANNEL WITH A CENTRAL BODY A. N. Ryabinin Department of Hydroaeromechanics, Faculty of Mathematics and Mechanics, Saint-Petersburg State University, St. Petersburg, Russia E-Mail: a.ryabinin@spbu.ru
More information(LSS Erlangen, Simon Bogner, Ulrich Rüde, Thomas Pohl, Nils Thürey in collaboration with many more
Parallel Free-Surface Extension of the Lattice-Boltzmann Method A Lattice-Boltzmann Approach for Simulation of Two-Phase Flows Stefan Donath (LSS Erlangen, stefan.donath@informatik.uni-erlangen.de) Simon
More informationComputation of Fictitious Gas Flow with Euler Equations
1 Computation of Fictitious Gas Flow with Euler Equations Pei Li and Helmut Sobieczky DLR Göttingen, Germany Abstract The Fictitious Gas Concept supports some computational design methods to construct
More informationAlgorithmic Developments in TAU
Algorithmic Developments in TAU Ralf Heinrich, Richard Dwight, Markus Widhalm, and Axel Raichle DLR Institute of Aerodynamics and Flow Technology, Lilienthalplatz 7, 38108, Germany ralf.heinrich@dlr.de,
More informationMid-Year Report. Discontinuous Galerkin Euler Equation Solver. Friday, December 14, Andrey Andreyev. Advisor: Dr.
Mid-Year Report Discontinuous Galerkin Euler Equation Solver Friday, December 14, 2012 Andrey Andreyev Advisor: Dr. James Baeder Abstract: The focus of this effort is to produce a two dimensional inviscid,
More informationNIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011
NIA CFD Seminar, October 4, 2011 Hyperbolic Seminar, NASA Langley, October 17, 2011 First-Order Hyperbolic System Method If you have a CFD book for hyperbolic problems, you have a CFD book for all problems.
More informationImplementing third order compressible flow solver for hexahedral meshes in OpenFoam
Tutorial/Report in OpenFoam Course 8 Implementing third order compressible flow solver for hexahedral meshes in OpenFoam Martin Olausson, Chalmers University of Technology, SE-1 9 Gothenburg, Sweden Abstract
More informationAn Object-Oriented Serial and Parallel DSMC Simulation Package
An Object-Oriented Serial and Parallel DSMC Simulation Package Hongli Liu and Chunpei Cai Department of Mechanical and Aerospace Engineering, New Mexico State University, Las Cruces, New Mexico, 88, USA
More informationCFD-1. Introduction: What is CFD? T. J. Craft. Msc CFD-1. CFD: Computational Fluid Dynamics
School of Mechanical Aerospace and Civil Engineering CFD-1 T. J. Craft George Begg Building, C41 Msc CFD-1 Reading: J. Ferziger, M. Peric, Computational Methods for Fluid Dynamics H.K. Versteeg, W. Malalasekara,
More informationNumerical Analysis of Shock Tube Problem by using TVD and ACM Schemes
Numerical Analysis of Shock Tube Problem by using TVD and Schemes Dr. Mukkarum Husain, Dr. M. Nauman Qureshi, Syed Zaid Hasany IST Karachi, Email: mrmukkarum@yahoo.com Abstract Computational Fluid Dynamics
More informationOn the Resolution Necessary to Capture Dynamics of Unsteady Detonation
On the Resolution Necessary to Capture Dynamics of Unsteady Detonation Christopher M. Romick, University of Notre Dame, Notre Dame, IN Tariq D. Aslam, Los Alamos National Laboratory, Los Alamos, NM and
More informationNonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws
Nonoscillatory Central Schemes on Unstructured Triangulations for Hyperbolic Systems of Conservation Laws Ivan Christov Bojan Popov Department of Mathematics, Texas A&M University, College Station, Texas
More informationEXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS
EXPLICIT AND IMPLICIT TVD AND ENO HIGH RESOLUTION ALGORITHMS APPLIED TO THE EULER AND NAVIER-STOKES EQUATIONS IN THREE-DIMENSIONS RESULTS Edisson Sávio de Góes Maciel, edissonsavio@yahoo.com.br Mechanical
More informationThis is an author-deposited version published in: Eprints ID: 4362
This is an author-deposited version published in: http://oatao.univ-toulouse.fr/ Eprints ID: 4362 To cite this document: CHIKHAOUI Oussama, GRESSIER Jérémie, GRONDIN Gilles. Assessment of the Spectral
More informationExpress Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil
Express Introductory Training in ANSYS Fluent Workshop 04 Fluid Flow Around the NACA0012 Airfoil Dimitrios Sofialidis Technical Manager, SimTec Ltd. Mechanical Engineer, PhD PRACE Autumn School 2013 -
More informationFAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD
TASK QUARTERLY 12 No 3, 273 287 FAST ALGORITHMS FOR CALCULATIONS OF VISCOUS INCOMPRESSIBLE FLOWS USING THE ARTIFICIAL COMPRESSIBILITY METHOD ZBIGNIEW KOSMA Institute of Applied Mechanics, Technical University
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationNumerical and theoretical analysis of shock waves interaction and reflection
Fluid Structure Interaction and Moving Boundary Problems IV 299 Numerical and theoretical analysis of shock waves interaction and reflection K. Alhussan Space Research Institute, King Abdulaziz City for
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationComputation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables
Computation of Sensitivity Derivatives of Navier-Stokes Equations using Complex Variables By Veer N. Vatsa NASA Langley Research Center, Hampton, VA 23681 Mail Stop 128, email: v.n.vatsa@larc.nasa.gov
More informationPossibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics
Possibility of Implicit LES for Two-Dimensional Incompressible Lid-Driven Cavity Flow Based on COMSOL Multiphysics Masanori Hashiguchi 1 1 Keisoku Engineering System Co., Ltd. 1-9-5 Uchikanda, Chiyoda-ku,
More informationTutorial 2. Modeling Periodic Flow and Heat Transfer
Tutorial 2. Modeling Periodic Flow and Heat Transfer Introduction: Many industrial applications, such as steam generation in a boiler or air cooling in the coil of an air conditioner, can be modeled as
More informationVerification of Moving Mesh Discretizations
Verification of Moving Mesh Discretizations Krzysztof J. Fidkowski High Order CFD Workshop Kissimmee, Florida January 6, 2018 How can we verify moving mesh results? Goal: Demonstrate accuracy of flow solutions
More informationIntroduction to CFX. Workshop 2. Transonic Flow Over a NACA 0012 Airfoil. WS2-1. ANSYS, Inc. Proprietary 2009 ANSYS, Inc. All rights reserved.
Workshop 2 Transonic Flow Over a NACA 0012 Airfoil. Introduction to CFX WS2-1 Goals The purpose of this tutorial is to introduce the user to modelling flow in high speed external aerodynamic applications.
More informationA Hybrid Cartesian Grid and Gridless Method for Compressible Flows
rd AIAA Aerospace Sciences Meeting and Exhibit,, January 5, Reno, Nevada A Hybrid Cartesian Grid and Gridless Method for Compressible Flows Hong Luo and Joseph D. Baum Science Applications International
More informationTHE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS
March 18-20, 2013 THE EFFECTS OF THE PLANFORM SHAPE ON DRAG POLAR CURVES OF WINGS: FLUID-STRUCTURE INTERACTION ANALYSES RESULTS Authors: M.R. Chiarelli, M. Ciabattari, M. Cagnoni, G. Lombardi Speaker:
More informationFaculty of Mechanical and Manufacturing Engineering, University Tun Hussein Onn Malaysia (UTHM), Parit Raja, Batu Pahat, Johor, Malaysia
Applied Mechanics and Materials Vol. 393 (2013) pp 305-310 (2013) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.393.305 The Implementation of Cell-Centred Finite Volume Method
More informationNumerical Analysis of a Blast Wave Using CFD-CAA Hybrid Method
Numerical Analysis of a Blast Wave Using CFD-CAA Hybrid Method In Cheol Lee * and Duck-Joo Lee. Korea Advanced Institute of Science and Technology, Daejeon, 305-701, Republic of Korea Sung Ho Ko and Dong
More informationA STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS
HEFAT2012 9 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 16 18 July 2012 Malta A STUDY ON THE UNSTEADY AERODYNAMICS OF PROJECTILES IN OVERTAKING BLAST FLOWFIELDS Muthukumaran.C.K.
More informationParallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes
Parallelization study of a VOF/Navier-Stokes model for 3D unstructured staggered meshes L. Jofre, O. Lehmkuhl, R. Borrell, J. Castro and A. Oliva Corresponding author: cttc@cttc.upc.edu Centre Tecnològic
More informationA BGK-based Two-Equation Turbulence Model Algorithm for Solving Compressible Navier-Stokes Equations
Ong et al. CFD Letters Vol. 6(2) 2013 www.cfdl.issres.net Vol. 6 (2) June 2013 A BGK-based Two-Equation Turbulence Model Algorithm for Solving Compressible Navier-Stokes Equations Jiunn Chit Ong 1, Ashraf
More informationA COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS
A COUPLED FINITE VOLUME SOLVER FOR THE SOLUTION OF LAMINAR TURBULENT INCOMPRESSIBLE AND COMPRESSIBLE FLOWS L. Mangani Maschinentechnik CC Fluidmechanik und Hydromaschinen Hochschule Luzern Technik& Architektur
More informationDevelopment of a Maxwell Equation Solver for Application to Two Fluid Plasma Models. C. Aberle, A. Hakim, and U. Shumlak
Development of a Maxwell Equation Solver for Application to Two Fluid Plasma Models C. Aberle, A. Hakim, and U. Shumlak Aerospace and Astronautics University of Washington, Seattle American Physical Society
More informationMesh deformation and shock capturing techniques for high-order simulation of unsteady compressible flows on dynamic meshes
AIAA SciTech 5-9 January 2015, Kissimmee, Florida 53rd AIAA Aerospace Sciences Meeting AIAA 2015-1741 Mesh deformation and shock capturing techniques for high-order simulation of unsteady compressible
More informationSPH: Towards the simulation of wave-body interactions in extreme seas
SPH: Towards the simulation of wave-body interactions in extreme seas Guillaume Oger, Mathieu Doring, Bertrand Alessandrini, and Pierre Ferrant Fluid Mechanics Laboratory (CNRS UMR6598) Ecole Centrale
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationNumerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods
Numerical Simulation of Coupled Fluid-Solid Systems by Fictitious Boundary and Grid Deformation Methods Decheng Wan 1 and Stefan Turek 2 Institute of Applied Mathematics LS III, University of Dortmund,
More informationContinuum-Microscopic Models
Scientific Computing and Numerical Analysis Seminar October 1, 2010 Outline Heterogeneous Multiscale Method Adaptive Mesh ad Algorithm Refinement Equation-Free Method Incorporates two scales (length, time
More informationApplication of the MCMC Method for the Calibration of DSMC Parameters
Application of the MCMC Method for the Calibration of DSMC Parameters James S. Strand and David B. Goldstein Aerospace Engineering Dept., 1 University Station, C0600, The University of Texas at Austin,
More informationInvestigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM)
Computational Methods and Experimental Measurements XVII 235 Investigation of cross flow over a circular cylinder at low Re using the Immersed Boundary Method (IBM) K. Rehman Department of Mechanical Engineering,
More informationComputation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow
Excerpt from the Proceedings of the COMSOL Conference 8 Boston Computation of Velocity, Pressure and Temperature Distributions near a Stagnation Point in Planar Laminar Viscous Incompressible Flow E. Kaufman
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationPreliminary Spray Cooling Simulations Using a Full-Cone Water Spray
39th Dayton-Cincinnati Aerospace Sciences Symposium Preliminary Spray Cooling Simulations Using a Full-Cone Water Spray Murat Dinc Prof. Donald D. Gray (advisor), Prof. John M. Kuhlman, Nicholas L. Hillen,
More informationParallel Direct Simulation Monte Carlo Computation Using CUDA on GPUs
Parallel Direct Simulation Monte Carlo Computation Using CUDA on GPUs C.-C. Su a, C.-W. Hsieh b, M. R. Smith b, M. C. Jermy c and J.-S. Wu a a Department of Mechanical Engineering, National Chiao Tung
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationVerification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard
Verification and Validation in CFD and Heat Transfer: ANSYS Practice and the New ASME Standard Dimitri P. Tselepidakis & Lewis Collins ASME 2012 Verification and Validation Symposium May 3 rd, 2012 1 Outline
More informationA Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids
A Multigrid Block LU-SGS Algorithm for Euler Equations on Unstructured Grids Ruo Li, Xin Wang and Weibo Zhao January 26, 27 Abstract We proposed an efficient and robust algorithm to solve the steady Euler
More informationHigh-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering
High-order solutions of transitional flow over the SD7003 airfoil using compact finite-differencing and filtering Daniel J. Garmann and Miguel R. Visbal Air Force Research Laboratory, Wright-Patterson
More informationA fully implicit Navier-Stokes algorithm for unstructured grids incorporating a two-equation turbulence model
Copyright 1996, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1996 A9618376, AIAA Paper 96-0414 A fully implicit Navier-Stokes algorithm for unstructured
More informationThermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA
Thermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA Jingxiao Xu 1, Jason Wang 2 1 LSTC 2 LSTC 1 Abstract Smooth particles hydrodynamics is a meshfree, Lagrangian particle method
More informationCFD Post-Processing of Rampressor Rotor Compressor
Gas Turbine Industrial Fellowship Program 2006 CFD Post-Processing of Rampressor Rotor Compressor Curtis Memory, Brigham Young niversity Ramgen Power Systems Mentor: Rob Steele I. Introduction Recent movements
More informationNon-Newtonian Transitional Flow in an Eccentric Annulus
Tutorial 8. Non-Newtonian Transitional Flow in an Eccentric Annulus Introduction The purpose of this tutorial is to illustrate the setup and solution of a 3D, turbulent flow of a non-newtonian fluid. Turbulent
More informationTwo-dimensional laminar shock wave / boundary layer interaction
Two-dimensional laminar shock wave / boundary layer interaction J.-Ch. Robinet (), V. Daru (,) and Ch. Tenaud () () SINUMEF Laboratory, ENSAM-PARIS () LIMSI-CNRS 5, Bd. de l Hôpital, PARIS 753, France
More informationAvailable online at ScienceDirect. Procedia Engineering 99 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 99 (2015 ) 575 580 APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology, APISAT2014 A 3D Anisotropic
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
1 DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Direct Simulation-Based Study of Radiance in a Dynamic Ocean LONG-TERM GOALS Dick K.P. Yue Center for Ocean Engineering
More informationA High-Order Accurate Unstructured GMRES Solver for Poisson s Equation
A High-Order Accurate Unstructured GMRES Solver for Poisson s Equation Amir Nejat * and Carl Ollivier-Gooch Department of Mechanical Engineering, The University of British Columbia, BC V6T 1Z4, Canada
More informationLecture 1: Finite Volume WENO Schemes Chi-Wang Shu
Lecture 1: Finite Volume WENO Schemes Chi-Wang Shu Division of Applied Mathematics Brown University Outline of the First Lecture General description of finite volume schemes for conservation laws The WENO
More informationMESHLESS METHOD FOR SIMULATION OF COMPRESSIBLE REACTING FLOW
MESHLESS METHOD FOR SIMULATION OF COMPRESSIBLE REACTING FLOW Jin Young Huh*, Kyu Hong Kim**, Suk Young Jung*** *Department of Mechanical & Aerospace Engineering, Seoul National University, **Department
More informationFinite Volume Discretization on Irregular Voronoi Grids
Finite Volume Discretization on Irregular Voronoi Grids C.Huettig 1, W. Moore 1 1 Hampton University / National Institute of Aerospace Folie 1 The earth and its terrestrial neighbors NASA Colin Rose, Dorling
More informationEfficient Unstructured Mesh Flow Solver Using Hamiltonian Paths and Strand Grids
Efficient Unstructured Mesh Flow Solver Using Hamiltonian Paths and Strand Grids *Yong Su Jung 1), Bharath Govindarajan 2) and James Baeder 3) 1), 2), 3) Department of Aerospace Engineering, University
More informationCARTESIAN grids have been gaining popularity in
42nd AIAA Aerospace Sciences Meeting and Exhibit 5-8 January 24, Reno, Nevada AIAA 24-581 An Embedded Boundary Cartesian Grid Scheme for Viscous Flos using a Ne Viscous Wall Boundary Condition Treatment
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More informationA Scalable GPU-Based Compressible Fluid Flow Solver for Unstructured Grids
A Scalable GPU-Based Compressible Fluid Flow Solver for Unstructured Grids Patrice Castonguay and Antony Jameson Aerospace Computing Lab, Stanford University GTC Asia, Beijing, China December 15 th, 2011
More informationModule 1: Introduction to Finite Difference Method and Fundamentals of CFD Lecture 13: The Lecture deals with:
The Lecture deals with: Some more Suggestions for Improvement of Discretization Schemes Some Non-Trivial Problems with Discretized Equations file:///d /chitra/nptel_phase2/mechanical/cfd/lecture13/13_1.htm[6/20/2012
More informationRecent developments for the multigrid scheme of the DLR TAU-Code
www.dlr.de Chart 1 > 21st NIA CFD Seminar > Axel Schwöppe Recent development s for the multigrid scheme of the DLR TAU-Code > Apr 11, 2013 Recent developments for the multigrid scheme of the DLR TAU-Code
More informationSensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method
2th AIAA Computational Fluid Dynamics Conference 27-3 June 2, Honolulu, Hawaii AIAA 2-348 Sensitivity Analysis for the Compressible Navier-Stokes Equations Using a Discontinuous Galerkin Method Li Wang
More informationMulti-Element High-Lift Configuration Design Optimization Using Viscous Continuous Adjoint Method
JOURNAL OF AIRCRAFT Vol. 41, No. 5, September October 2004 Multi-Element High-Lift Configuration Design Optimization Using Viscous Continuous Adjoint Method Sangho Kim, Juan J. Alonso, and Antony Jameson
More informationFEMLAB Exercise 1 for ChE366
FEMLAB Exercise 1 for ChE366 Problem statement Consider a spherical particle of radius r s moving with constant velocity U in an infinitely long cylinder of radius R that contains a Newtonian fluid. Let
More informationThermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA
Thermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA INTRODUCTION: Jingxiao Xu, Jason Wang LSTC Heat transfer is very important in many industrial and geophysical problems. Many
More informationA MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES
A MULTI-DOMAIN ALE ALGORITHM FOR SIMULATING FLOWS INSIDE FREE-PISTON DRIVEN HYPERSONIC TEST FACILITIES Khalil Bensassi, and Herman Deconinck Von Karman Institute for Fluid Dynamics Aeronautics & Aerospace
More informationAn Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm
An Optimization Method Based On B-spline Shape Functions & the Knot Insertion Algorithm P.A. Sherar, C.P. Thompson, B. Xu, B. Zhong Abstract A new method is presented to deal with shape optimization problems.
More informationTVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II
TVD Flux Vector Splitting Algorithms Applied to the Solution of the Euler and Navier-Stokes Equations in Three-Dimensions Part II EDISSON SÁVIO DE GÓES MACIEL IEA- Aeronautical Engineering Division ITA
More informationNumerical Methods for Hyperbolic and Kinetic Equations
Numerical Methods for Hyperbolic and Kinetic Equations Organizer: G. Puppo Phenomena characterized by conservation (or balance laws) of physical quantities are modelled by hyperbolic and kinetic equations.
More informationAn Implicit Gradient Method for Cell-Centered Finite-Volume Solver on Unstructured Grids
AIAA SciTech Forum 7-11 January 2019, San Diego, California AIAA Scitech 2019 Forum 10.2514/6.2019-1155 An Implicit Gradient Method for Cell-Centered Finite-Volume Solver on Unstructured Grids Hiroaki
More informationNUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING
Review of the Air Force Academy No.3 (35)/2017 NUMERICAL 3D TRANSONIC FLOW SIMULATION OVER A WING Cvetelina VELKOVA Department of Technical Mechanics, Naval Academy Nikola Vaptsarov,Varna, Bulgaria (cvetelina.velkova1985@gmail.com)
More informationWe present a high-order accurate space-time discontinuous Galerkin method for solving two-dimensional
A High-Order Discontinuous Galerkin Method with Unstructured Space-Time Meshes for Two-Dimensional Compressible Flows on Domains with Large Deformations Luming Wang a, Per-Olof Persson a, a Department
More informationEstimating Vertical Drag on Helicopter Fuselage during Hovering
Estimating Vertical Drag on Helicopter Fuselage during Hovering A. A. Wahab * and M.Hafiz Ismail ** Aeronautical & Automotive Dept., Faculty of Mechanical Engineering, Universiti Teknologi Malaysia, 81310
More informationModeling & Simulation of Supersonic Flow Using McCormack s Technique
Modeling & Simulation of Supersonic Flow Using McCormack s Technique M. Saif Ullah Khalid*, Afzaal M. Malik** Abstract In this work, two-dimensional inviscid supersonic flow around a wedge has been investigated
More information