Euler Equations Lab AA Computer Project 2
|
|
- Emmeline McKenzie
- 5 years ago
- Views:
Transcription
1 Euler Equations Lab AA Computer Project 2 Mishaal Aleem February 26, 2015 Contents 1 Introduction Algorithms Lax-Wendroff MacCormack Roe (1st-order TVD scheme) Harten-Yee (2nd-order TVD scheme) Shock Tube Theoretical Numerical Supersonic Diverging Nozzle Theoretical Numerical Conclusion References Appendix A-1 1 Introduction The Euler Equations Lab is a MATLAB computational fluid dynamics (CFD) program that allows the user to study the behavior of several algorithms and compare the results to those that are physically expected for the pseudo-one-dimensional Euler equations as applied to a shock tube and a nozzle. 2 Algorithms The pseudo-one-dimensional Euler equations are written as where Q t + F x = S ρa ρua Q = ρua F = (ρu 2 + p)a S = ea u(e + p)a Each of the algorithms described briefly in the following sections were used to solve these equations for the case of a Sod shock tube and a nozzle. 2.1 Lax-Wendroff The Lax-Wendroff algorithm is a centered differencing scheme with artificial viscosity for stability that is applicable to the one-dimensional nonlinear Euler equations. However, this algorithm does not have enough dissipation to solve the psuedo-one-dimensional equations with even moderate expansions. 0 p da dx 0 1
2 2.2 MacCormack The MacCormack algorithm is a predictor-correct method. It is easily applicable to the pseudo-onedimensional Euler equations because it simply modifies the standard one-dimensional Euler MacCormack solver with a source term for S. 2.3 Roe (1st-order TVD scheme) Unlike the Lax-Wendroff and MacCormack algorithms which are centered-differencing methods with diffusion applied for stability, the Roe method is a first order total variation diminishing (TVD) scheme that applies the principle of flux-splitting into positive and negative components. 2.4 Harten-Yee (2nd-order TVD scheme) The Harten-Yee algorithm is also a TVD flux splitting method and is of second order. 3 Shock Tube For the Sod shock tube, the area is set to A = 1 throughout the nozzle making da dx = 0 and reducing the psuedo-one-dimensional Euler equations to the standard unsteady one-dimensional Euler equations. In the shock tube problem, a tube is filled with a gas and has a diaphragm in the middle. The gas is at different pressures and densities on either side of the diaphragm at t < 0. At t = 0 the diaphragm is removed. A diagram for the resulting motion in the tube is presented in Fig. 1. Figure 1: Motion in Sod shock tube [2] As can be seen from Fig 1, there is a left-moving shock wave within the tube resulting from the pressure difference between regions 1 and 4. For this shock tube, the diaphragm pressure ratio is given as p 4 = 4. The gas is air with a heat capacity ratio of γ = 1.4. The initial pressures are set to p 4 = 4/γ and p 4 = 1/γ and the initial densities are set to ρ 4 = 4 and ρ 1 = Theoretical The speed of the shock wave was found analytically using the shock strength and shock tube relationships. First, the speed of sound in regions 1 and 4 was calculated using Eq. 1. γp a = (1) ρ Note that because of the initial conditions, for region 1, a 1 = 1. Next, the pressure ratio between regions 1 and 2, p 2, was calculated using Eq. 2 and the known diaphragm pressure ratio. p 4 = p 2 [ ] 2γ 1 (γ 1)(a γ 1 1/a 4 )(p2/p1 1) 2γ 2γ + (γ + 1)(p2/p1 1) (2) 2
3 Next, Eq. 3 was used to find the Mach number of the shock wave, M s. γ 1 M s = 2γ + γ + 1 p 2 (3) 2γ Finally, the shock Mach was converted to speed using Eq. 4. u s = M s a 1 (4) Per this analysis, the shock speed was found to be u s = Note that because a 1 = 1, the numerical value of M s is equivalent to the numerical value of u s, and thus in the discussion, the physical units of the problem may be ignored. 3.2 Numerical The Euler Equations Lab was used to study each algorithms ability to simulate the conditions in the shock tube over time. To find the numerical shock speed, the program was run for 100 iterations and 200 iterations for each algorithm, and the time was recorded at the end of each run. The shock location was measured for each run using the pressure plots (for consistency). Finally, the shock speed was calculated using the definition of speed, Eq. 5. u s = x 2 x 1 (5) t 2 t 1 As an example, the plots used to calculate the numerical shock speed for the MacCormack algorithm is presented in Fig. 2. Figure 2: Shock speed calculation for MacCormack algorithm Presented in Tab. 1 are the numerical shock speed calculations from each algorithm and their errors. Note that all are within 3.13% of the theoretical value, supporting the ability of these algorithms to simulate the Sod shock tube problem well. Table 1: Shock Speed Calculations Algorithm MacCormack Lax-Wendroff Roe Harten-Yee Shock Speed Error 0.24% 1.42% 0.34% 3.13% Each algorithm was able to simulate the physics inside the Sod shock tube, though each did exhibit certain non-physical numerical computation phenomena. Figure 3 shows the plot of the pressure in the shock tube after 180 iterations for CFL=0.5 for all the algorithms. 3
4 Figure 3: Pressure in shock tube The MacCormack algorithm was clearly able to capture the shock, however it was visibly dispersive, as seen Fig 3. As the CFL value was decreased, the dispersion got worse, represented in Fig. 4. Figure 4: Density in shock tube, MacCormack algorithm, 300 grid points The Lax-Wendroff was very similar in behavior to MacCormack; the algorithm was able to capture the shock, however it was at the expense of dispersion which got worse as the CFL number was decreased from CFL=1 to CFL=0. For both the MacCormack and the Lax-Wendroff algorithm, because of the dispersion there are magnitude oscillations in the solutions for pressure, density, and energy. Dispersion is useful for shock capturing as it keeps the shock jump discontinuity relatively vertical and the expansion fan relatively linear. However, this is at the expense of diminishing true local information with the oscillations. The Roe algorithm, however, was more visibly diffusive. Shown in Fig 5 is the density solution using the Roe algorithm at t = 0.8 for various CFL numbers. Clearly, the diffusion did increase as the CFL number was varied, though not significantly. Diffusion is detrimental to shock capturing because shocks are locations of sharp discontinuity, but diffusion causes the discontinuity to spread over grid points rather than occurring at the single true location. As seen in Fig. 5, this causes the shock and contact discontinuity to be significantly non-vertical and the expansion fan to be non-linear (particularly near the top and bottom of the line). However, this diffusion keeps the algorithm very stable. 4
5 Figure 5: Density in shock tube, Roe algorithm, 300 grid points The Harten-Yee algorithm also had slight diffusion for the shock tube problem. There was also visible dispersion for higher CFL numbers. Interestingly, the dispersion became more pronounced as the grid points were increased ( 500+) for larger CFL numbers (.8+). This is represented in Fig. 6. Figure 6: Density in shock tube, Harten-Yee algorithm While all the algorithms were able to capture the shock, each had a trade-off with either diffusion or dispersion to varying degrees in the solution. 4 Supersonic Diverging Nozzle For the supersonic diverging nozzle problem, the area change of a nozzle is given by Eq. 6 for x=[0,10]. A(x) = tanh(0.8x4) (6) The inflow Mach number is 1.26 and the exit pressure is fixed such that the ratio gas is air with a heat capacity ratio of γ = Theoretical p exit p entrance = The The shock location for this nozzle was be found analytically by employing isentropic relations, shock jump relations, and principles of supersonic diverging nozzle. 5
6 First, Eq. 6 was used to calculate the area at the inlet, x = 0. Then, Eq. 7 was used to calculate A. ( ) A 2 A = 1 [ 2 M 2 γ + 1 (1 + γ 1 M 2 ) 2 ] γ+1 γ 1 (7) Next, a shock location, x s, was assumed. From this assumed location the area at the shock, A s, was calculated using Eq. 6 and from Eq. 7 an upstream Mach number, M 1, was calculated. Then, using the shock jump relations given by Eq. 8, the downstream Mach number, M 2, was calculated, with the subsonic solution chosen. M2 2 = 1 + γ 1 2 M 1 2 γm1 2 γ 1 (8) 2 Across the shock, the pressure jump was calculated using the Eq. 9 normal shock pressure jump relation. p 2 = 2γM12 (γ 1) (9) γ + 1 Using the downstream M 2 at x s and the nozzle area A s, A was once again calculated. Then, using the nozzle area at the exit calculated by Eq. 6 at x = 10, and the second calculated A, the exit Mach number was calculated. The total pressure to static pressure was calculated at the entrance, upstream of the shock, downstream of the shock, and at the exit using Eq. 10 and the Mach number at each of these locations. ( p t p = 1 + γ 1 ) γ M 2 2 p Finally, the pressure ratio exit p entrance was calculated using Eq. 11 and knowing that p t,exit p p t,2 = t,1 p t,entrance = 1. If the calculated value was higher than 1.931, the process was repeated using a shock location assumption closer to the entrance, and vice-versa. γ 1 (10) p exit p entrance = p exit p t,exit p t,exit p t,2 p t,2 p 2 p 2 p t,1 p t,1 p t,entrance p t,entrance p entrance (11) Per this analytical method, the shock location was found to be x s = The nozzle and shock location are represented in Fig 7. Figure 7: Supersonic Diverging Nozzle 4.2 Numerical The Euler Equations Lab was used to study the ability of the various algorithms to capture the supersonic diverging nozzle behavior and compare the numerical results to the theoretical shock location. 6
7 Expectedly, the Lax-Wendroff algorithm was unacceptable for the nozzle problem. This algorithm cannot handle even moderate area expansions which, as shown in the Fig 7 diagram of the nozzle, do occur for this nozzle. After just 50 iterations, the instability became apparent as the solution simply grew unchecked and, of course, did not capture the shock location. The Lax-Wendroff numerical solution after 1000 iterations is presented in Fig. 8. Figure 8: Lax-Wendroff Algorithm for Nozzle The other algorithms, however, were appropriate for the nozzle. Each of them were able to successfully capture the behavior in the nozzle. Reaching the steady solution, however, did take thousands of iterations. For stability, CFL numbers of greater than 1 were not used, and 1 was avoided due to potential non-linearity-induced stability issues. The MacCormack algorithm successfully simulated a solution, however with visible dispersion. With this method, after about iterations for a 600 grid point solution with CFL number = 0.99, the shock location was found to be 4.908, presented in Fig. 9. As the CFL number was incrementally decreased from CFL=1 to CFL=0, the solution became more and more dispersive. While varying the CFL number did change the amount of dispersion, the variations were not as significant as for the CFL variations seen for some algorithms in the Linear Advection Lab. Figure 9: MacCormack Algorithm for Nozzle The Roe algorithm neatly captured the behavior of the nozzle problem, and with much less numerical computation noise than the MacCormack algorithm. With this method, after about iterations for 7
8 a 600 grid point solution with CFL number = 0.99, the shock location was again found to be 4.908, presented in Fig. 10. Similar to the MacCormack algorithm, as the CFL number was incrementally decreased from CFL=1 to CFL=0, the solution became more and more dispersive. Figure 10: Roe Algorithm for Nozzle The Harten-Yee method was more unstable than the previous two methods. For large CFL numbers, the solution had a lot of noise. For this method, after about iterations for a 600 grid point solution with CFL number = 0.5, the shock location was found to be Figure 11: Harten-Yee Algorithm for Nozzle The shock location as computed by each algorithm is presented along with error value in Tab. 2. Note again the low error values, all less that 0.13%, pointing to the success of these algorithms in capturing the physical behavior in the nozzle. Table 2: Shock Location Calculations Algorithm MacCormack Roe Harten-Yee Shock Location Error 0.09% 0.09% 0.13% 8
9 5 Conclusion The Euler Equations Lab was used to study behavior in a Sod shock tube and in a supersonic diverging nozzle using the pseudo-one-dimensional Euler equations. Theoretical solution values were compared to numerical results and the ability of various algorithms to capture the physical behavior was explored. For the Sod shock tube, the Lax-Wendroff, MacCormack, Roe, and Harten-Yee algorithms were found to be able to capture the shock. The Lax-Wendroff and MacCormack algorithms tended to have dispersive behavior in the solutions. The Roe algorithm had diffusion in the solution, and the Harten-Yee algorithm had significant dispersion for high grid points coupled with high CFL numbers. However, each algorithm was able to capture the shock speed within 3.13% of the theoretical value pointing to their overall success. For the supersonic diverging nozzle, Lax-Wendroff was found to be unstable because there was not enough damping to compensate for the increasing area of the nozzle. However, the other three algorithms were able to simulate the physical behavior in the nozzle. The MacCormack algorithm showed dispersion and the Harten-Yee was again very noisy for large CFL numbers. The Roe algorithm was arguably the best with minimal diffusion, even as the CFL number was decreased. Each of the three algorithms, however, were able to numerically simulate the shock location within 0.13% of the theoretical solution. Through working with the Euler Equations Lab,the fundamental principle of using CFD to model model fluid dynamics equations was reiterated: there is an inherent numerical trade-off when capturing a physical solution. Diffusion caused loss of information but ensured stability of the algorithms. Dispersion caused oscillations in the solution but allowed for simulation sharp discontinuities. Ultimately, the appropriate algorithm is dependent on the application, the desired information to be retained, and the desired accuracy of the solution globally and locally. References [1] Eberhardt, D.S. and Shumlak, U. AA543 Computational Fluid Dynamics I, University of Washington, Seattle, WA, [2] Liepmann, H.W. and Roshko, A. Elements of Gasdynamics, Dover Publications, Inc., Mineola, NY,
10 Appendix MacCormack Algorithm The MATLAB code used to implement the MacCormack algorithm in the Euler Equations Lab is presented below. if strcmp(algorithm,'maccormack')==1 % Start MacCormack Algorithm here. [flux]=calfx(q,flux,1,imax,gamma1); for i=2:imax for nx=1:3 qbar(i,nx) = (q(i,nx)) - (dt/dx)*(flux(i,nx)-flux(i-1,nx)); end end [flux]=calfx(qbar,flux,1,imax-1,gamma1); for i=2:imax-1 for nx = 1:3 dq(i,nx) = -q(i,nx)/2 +qbar(i,nx)/2 -(dt/dx)/2 * (flux(i+1,nx) - flux(i,nx)); end end for i = 2:imax rhoa = q(i,1) + dq(i,1); rhoua = q(i,2) + dq(i,2); ea = q(i,3) + dq(i,3); u = rhoua / rhoa; p = gamma1 * (ea - (.5 * rhoua * u)) / area(i,1); dq(i,2) = dq(i,2) + (dt * p * s(i,2)); end Shock Tube Shock Speed The MATLAB code used to calculate the shock speed for the Sod shock tube problem is presented below. %gamma g=1.4; %initial conditions p1=1/g; p4=4/g; rho1=1; rho4=4; %calculate speed of sound a1=sqrt(g*p1/rho1); a4=sqrt(g*p4/rho4); %solve for p2/p1 p4 p1=4; syms p2 p1 alpha=((g-1)*(a1/a4)*(p2 p1-1))/(sqrt(2*g)*sqrt(2*g+(g+1)*(p2 p1-1))); p2 p1=solve(p4 p1==p2 p1*(1-alpha)ˆ(-2*g/(g-1)),p2 p1); p2 p1=double(p2 p1(2)); %solve for shock Mach Ms=((g-1)/(2*g)+(g+1)/(2*g)*p2 p1)ˆ(1/2); %solve for shock speed us=ms*a1 A-1
11 Supersonic Nozzle Shock Location The MATLAB code used to calculate the shock location for the supersonic diverging nozzle is presented below. % Calculate A* from the LHS M0=1.26; A0= *tanh(0.8*0-4); syms Astar1 Astar1=solve((A0/Astar1)ˆ2==1/M0ˆ2*(2/(g+1)*(1+(g-1)/2*M0ˆ2))ˆ((g+1)/(g-1)),Astar1); Astar1=double(Astar1(1)); % calculate M1 based on a guess x x guess=4.9036; A xs= *tanh(0.8.*x guess-4); syms M1 M1=solve((A xs/astar1)ˆ2==1/m1ˆ2*(2/(g+1)*(1+(g-1)/2*m1ˆ2))ˆ((g+1)/(g-1)),m1); M1=double(M1(2)); % calculate M2 M2=sqrt((1+(g-1)/2*M1ˆ2)/(g*M1ˆ2-(g-1)/2)); % calculate A* from RHS syms Astar2 Astar2=solve((A xs/astar2)ˆ2==1/m2ˆ2*(2/(g+1)*(1+(g-1)/2*m2ˆ2))ˆ((g+1)/(g-1)),astar2); Astar2=double(Astar2(1)); % calculate Mexit syms Mexit Aend= *tanh(0.8*10-4); Mexit=solve((Aend/Astar2)ˆ2==1/Mexitˆ2*(2/(g+1)*(1+(g-1)/2*Mexitˆ2))ˆ((g+1)/(g-1)),Mexit); Mexit=double(Mexit(1)); % pressure ratios pe pte=((1+(g-1)/2*mexitˆ2)ˆ(g/(g-1)))ˆ(-1); pte pt2=1; pt2 p2=((1+(g-1)/2*m2ˆ2)ˆ(g/(g-1))); p2 p1=(2*g*m1ˆ2-(g-1))/(g+1); p1 pt1=((1+(g-1)/2*m1ˆ2)ˆ(g/(g-1)))ˆ(-1); pt1 pint=1; pint pin=((1+(g-1)/2*m0ˆ2)ˆ(g/(g-1))); calc pratio=pe pte*pte pt2*pt2 p2*p2 p1*p1 pt1*pt1 pint*pint pin A-2
Computer Project 3. AA Computational Fluid Dyanmics University of Washington. Mishaal Aleem March 17, 2015
Computer Project 3 AA 543 - Computational Fluid Dyanmics University of Washington Mishaal Aleem March 17, 2015 Contents Introduction........................................... 1 3.1 Grid Generator.......................................
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 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 informationEVALUATE SHOCK CAPTURING CAPABILITY WITH THE NUMERICAL METHODS IN OpenFOAM
THERMAL SCIENCE: Year 2013, Vol. 17, No. 4, pp. 1255-1260 1255 Open forum EVALUATE SHOCK CAPTURING CAPABILITY WITH THE NUMERICAL METHODS IN OpenFOAM by Reza KHODADADI AZADBONI a*, Mohammad Rahim MALEKBALA
More informationA Review on the Numerical Solution of the 1D Euler Equations. Hudson, Justin. MIMS EPrint:
A Review on the Numerical Solution of the D Euler Equations Hudson, Justin 6 MIMS EPrint: 6.9 Manchester Institute for Mathematical Sciences School of Mathematics The University of Manchester Reports available
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 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 informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationModeling Supersonic Jet Screech Noise Using Direct Computational Aeroacoustics (CAA) 14.5 Release
Modeling Supersonic Jet Screech Noise Using Direct Computational Aeroacoustics (CAA) 14.5 Release 2011 ANSYS, Inc. November 7, 2012 1 Workshop Advanced ANSYS FLUENT Acoustics Introduction This tutorial
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 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 informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationMathematical modeling of fluid flow using the numerical scheme with artificial viscosity
Mathematical modeling of fluid flow using the numerical scheme with artificial viscosity Ing. Tomáš Sommer, Ing. Martin Helmich Thesis supervised by: Doc. Svatomír Slavík, CSc., Doc. Luboš Janko, CSc.
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 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 information4A2 MODULE: WRITING AN EULER SOLVER FOR DUCT FLOW
1 4A2 MODULE: WRITING AN EULER SOLVER FOR DUCT FLOW INTRODUCTION The objective of this exercise is to write a program for solving the Euler equations for two dimensional internal flows using a time marching
More informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationSIMULATION OF A COMPARTMENTED AIRBAG DEPLOYMENT USING AN EXPLICIT, COUPLED EULER/LAGRANGE METHOD WITH ADAPTIVE EULER DOMAINS
SIMULATION OF A COMPARTMENTED AIRBAG DEPLOYMENT USING AN EXPLICIT, COUPLED EULER/LAGRANGE METHOD WITH ADAPTIVE EULER DOMAINS W.A. van der Veen MSC.Software Corporation Groningenweg 6 Gouda, The Netherlands
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 information2DNS Documentation. Release Sayop Kim
2DNS Documentation Release 0.0.1 Sayop Kim September 25, 2014 Contents 1 Project description 3 1.1 Given task................................................ 3 1.2 Governing Equations...........................................
More informationLax-Wendroff and McCormack Schemes for Numerical Simulation of Unsteady Gradually and Rapidly Varied Open Channel Flow
Archives of Hydro-Engineering and Environmental Mechanics Vol. 60 (2013), No. 1 4, pp. 51 62 DOI: 10.2478/heem-2013-0008 IBW PAN, ISSN 1231 3726 Lax-Wendroff and McCormack Schemes for Numerical Simulation
More informationA New Simulation Method for Airbags With Multiple Compartments (Adaptive Euler Domains)
4 th European LS -DYNA Users Conference LS-Dyna Environment II A New Simulation Method for Airbags With Multiple Compartments (Adaptive Euler Domains) Authors: R.S. Sadeghi and W.A. van der Veen* MSC.Software
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 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 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 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 informationNUMERICAL VISCOSITY. Convergent Science White Paper. COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved.
Convergent Science White Paper COPYRIGHT 2017 CONVERGENT SCIENCE. All rights reserved. This document contains information that is proprietary to Convergent Science. Public dissemination of this document
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 informationSupersonic Flow Over a Wedge
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 2 Supersonic Flow Over a Wedge Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification A uniform supersonic stream encounters a wedge
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 informationValidation of Internal Flow Prediction Codes for Solid Propellant Rocket Motors
Validation of Internal Flow Prediction Codes for Solid Propellant Rocket Motors H. Tuğrul TINAZTEPE, Mine E.YUMUŞAK ROKETSAN Missiles Industries Inc. Ankara-Samsun Karayolu 40.km PO. Box 30 Elmadağ 06780
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 informationRotorcraft Noise Prediction with Multi-disciplinary Coupling Methods. Yi Liu NIA CFD Seminar, April 10, 2012
Rotorcraft Noise Prediction with Multi-disciplinary Coupling Methods Yi Liu NIA CFD Seminar, April 10, 2012 Outline Introduction and Background Multi-disciplinary Analysis Approaches Computational Fluid
More informationDesign Optimization of a Subsonic Diffuser. for a Supersonic Aircraft
Chapter 5 Design Optimization of a Subsonic Diffuser for a Supersonic Aircraft 5. Introduction The subsonic diffuser is part of the engine nacelle leading the subsonic flow from the intake to the turbo-fan
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
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 informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
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 informationCarl Leake Embry-Riddle Aeronautical University,
College of Engineering Honors Program 8-11-2016 Testing the Accuracy of Commercial Computational Fluid Dynamics Codes: Vortex Transport by Uniform Flow and Transonic Ringleb Flow Problems Using ANSYS Fluent
More informationStudy on the Numerical Accuracy for the CFD
Study on the Numerical Accuracy for the CFD T.Yamanashi 1, H.Uchida, and M.Morita 1 Department of Mathematics, Master s Research Course of Faculty of Science, Tokyo University of Science,1-3 Kagurazaka,
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 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 informationThe WENO Method in the Context of Earlier Methods To approximate, in a physically correct way, [3] the solution to a conservation law of the form u t
An implicit WENO scheme for steady-state computation of scalar hyperbolic equations Sigal Gottlieb Mathematics Department University of Massachusetts at Dartmouth 85 Old Westport Road North Dartmouth,
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 informationDevelopment of the Compliant Mooring Line Model for FLOW-3D
Flow Science Report 08-15 Development of the Compliant Mooring Line Model for FLOW-3D Gengsheng Wei Flow Science, Inc. October 2015 1. Introduction Mooring systems are common in offshore structures, ship
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationChapter 6. Petrov-Galerkin Formulations for Advection Diffusion Equation
Chapter 6 Petrov-Galerkin Formulations for Advection Diffusion Equation In this chapter we ll demonstrate the difficulties that arise when GFEM is used for advection (convection) dominated problems. Several
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 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 informationCOMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY. Abstract. I. Introduction
COMPUTATIONAL AND EXPERIMENTAL INTERFEROMETRIC ANALYSIS OF A CONE-CYLINDER-FLARE BODY John R. Cipolla 709 West Homeway Loop, Citrus Springs FL 34434 Abstract A series of computational fluid dynamic (CFD)
More informationINTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012
INTERNATIONAL JOURNAL OF CIVIL AND STRUCTURAL ENGINEERING Volume 2, No 3, 2012 Copyright 2010 All rights reserved Integrated Publishing services Research article ISSN 0976 4399 Efficiency and performances
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
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 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 informationThe viscous forces on the cylinder are proportional to the gradient of the velocity field at the
Fluid Dynamics Models : Flow Past a Cylinder Flow Past a Cylinder Introduction The flow of fluid behind a blunt body such as an automobile is difficult to compute due to the unsteady flows. The wake behind
More informationEdge Detection Free Postprocessing for Pseudospectral Approximations
Edge Detection Free Postprocessing for Pseudospectral Approximations Scott A. Sarra March 4, 29 Abstract Pseudospectral Methods based on global polynomial approximation yield exponential accuracy when
More informationModeling Unsteady Compressible Flow
Tutorial 4. Modeling Unsteady Compressible Flow Introduction In this tutorial, FLUENT s density-based implicit solver is used to predict the timedependent flow through a two-dimensional nozzle. As an initial
More informationNumerical Analysis of the Bleed Slot Design of the Purdue Mach 6 Wind Tunnel
43rd AIAA Aerospace Sciences Meeting and Exhibit, Jan 10 13, Reno, Nevada Numerical Analysis of the Bleed Slot Design of the Purdue Mach 6 Wind Tunnel Ezgi S. Taskinoglu and Doyle D. Knight Dept of Mechanical
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 informationPrediction of Helicopter Blade- Vortex Interaction Noise using Motion Data from Experiment
Prediction of Helicopter Blade- Vortex Interaction Noise using Motion Data from Experiment Yoshinobu Inada, Choongmo Yang, and Takashi Aoyama Computational Science Research Group (CSRG) Institute of Aerospace
More informationStore Separation Simulation using Oct-tree Grid Based Solver
SAROD 2009 142 Symposium on Applied Aerodynamics and Design of Aerospace Vehicles (SAROD 2009) December 10-12, 2009, Bengaluru, India Store Separation Simulation using Oct-tree Grid Based Solver Saurabh
More informationSimulating Sinkage & Trim for Planing Boat Hulls. A Fluent Dynamic Mesh 6DOF Tutorial
Simulating Sinkage & Trim for Planing Boat Hulls A Fluent Dynamic Mesh 6DOF Tutorial 1 Introduction Workshop Description This workshop describes how to perform a transient 2DOF simulation of a planing
More informationDiscontinuous Fluctuation Distribution
Discontinuous Fluctuation Distribution Matthew Hubbard School of Computing, University of Leeds, Leeds, LS2 9JT, UK. Abstract This paper describes a new numerical scheme for the approximation of steady
More informationDetached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads
Detached Eddy Simulation Analysis of a Transonic Rocket Booster for Steady & Unsteady Buffet Loads Matt Knapp Chief Aerodynamicist TLG Aerospace, LLC Presentation Overview Introduction to TLG Aerospace
More informationOn the order of accuracy and numerical performance of two classes of finite volume WENO schemes
On the order of accuracy and numerical performance of two classes of finite volume WENO schemes Rui Zhang, Mengping Zhang and Chi-Wang Shu November 29, 29 Abstract In this paper we consider two commonly
More informationBASICS OF FLUID MECHANICS AND INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS
BASICS OF FLUID MECHANICS AND INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS Numerical Methods and Algorithms Volume 3 Series Editor: Claude Brezinski Université des Sciences et Technologies de Lille, France
More informationComparison of Classic and Finned Piston Reciprocating Linear Air Compressor Using COMSOL Multiphysics
Comparison of Classic and Finned Piston Reciprocating Linear Air Compressor Using COMSOL Multiphysics M. Heidari*, P. Barrade, and A. Rufer LEI, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
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 informationSTAR-CCM+: Wind loading on buildings SPRING 2018
STAR-CCM+: Wind loading on buildings SPRING 2018 1. Notes on the software 2. Assigned exercise (submission via Blackboard; deadline: Thursday Week 3, 11 pm) 1. NOTES ON THE SOFTWARE STAR-CCM+ generates
More informationScuola Politecnica DIME
Scuola Politecnica DIME Ingegneria Meccanica - Energia e Aeronautica Anno scolastico 2017-2018 Fluidodinamica Avanzata Aircraft S-shaped duct geometry optimization Professor Jan Pralits Supervisor Joel
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationPARALLEL DNS USING A COMPRESSIBLE TURBULENT CHANNEL FLOW BENCHMARK
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS Computational Fluid Dynamics Conference 2001 Swansea, Wales, UK, 4-7 September 2001 ECCOMAS PARALLEL DNS USING A COMPRESSIBLE
More informationCFD design tool for industrial applications
Sixth LACCEI International Latin American and Caribbean Conference for Engineering and Technology (LACCEI 2008) Partnering to Success: Engineering, Education, Research and Development June 4 June 6 2008,
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 informationSolver Settings. Introductory FLUENT Training ANSYS, Inc. All rights reserved. ANSYS, Inc. Proprietary
Solver Settings Introductory FLUENT Training 2006 ANSYS, Inc. All rights reserved. 2006 ANSYS, Inc. All rights reserved. 5-2 Outline Using the Solver Setting Solver Parameters Convergence Definition Monitoring
More informationNASA Rotor 67 Validation Studies
NASA Rotor 67 Validation Studies ADS CFD is used to predict and analyze the performance of the first stage rotor (NASA Rotor 67) of a two stage transonic fan designed and tested at the NASA Glenn center
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 informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationCoupling of STAR-CCM+ to Other Theoretical or Numerical Solutions. Milovan Perić
Coupling of STAR-CCM+ to Other Theoretical or Numerical Solutions Milovan Perić Contents The need to couple STAR-CCM+ with other theoretical or numerical solutions Coupling approaches: surface and volume
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationHigh Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study
Commun. Comput. Phys. doi:.48/cicp.375.6a Vol., No. 4, pp. 835-869 October 6 High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study Liang
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
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 informationBackward facing step Homework. Department of Fluid Mechanics. For Personal Use. Budapest University of Technology and Economics. Budapest, 2010 autumn
Backward facing step Homework Department of Fluid Mechanics Budapest University of Technology and Economics Budapest, 2010 autumn Updated: October 26, 2010 CONTENTS i Contents 1 Introduction 1 2 The problem
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 informationOptimal design of supersonic nozzle contour for altitude test facility
Journal of Mechanical Science and Technology 26 (8) (2012) 2589~2594 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-012-0634-x Optimal design of supersonic nozzle contour for altitude test facility
More informationLocal Fourier-spectral filtering for non-linear stabilization of high-order Flux Reconstruction schemes
Local Fourier-spectral filtering for non-linear stabilization of high-order Flux Reconstruction schemes Manuel R. López Morales, Kartikey Asthana Aerospace Computing Laboratory Stanford University November
More informationEVALUATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE URBAN FLOW
EVALATION OF A GENERAL CFD-SOLVER FOR A MICRO-SCALE RBAN FLOW Jarkko Saloranta and Antti Hellsten Helsinki niversity of Technology, Laboratory of Aerodynamics, Finland INTRODCTION In this work we study
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 informationVector Visualization
Vector Visualization Vector Visulization Divergence and Vorticity Vector Glyphs Vector Color Coding Displacement Plots Stream Objects Texture-Based Vector Visualization Simplified Representation of Vector
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 informationA Parallel Explicit/Implicit Time Stepping Scheme on Block-Adaptive Grids
A Parallel Explicit/Implicit Time Stepping Scheme on Block-Adaptive Grids Gábor Tóth, Darren L. De Zeeuw, Tamas I. Gombosi, Kenneth G. Powell Center for Space Environment Modeling, University of Michigan,
More informationANSYS AIM Tutorial Compressible Flow in a Nozzle
ANSYS AIM Tutorial Compressible Flow in a Nozzle Author(s): Sebastian Vecchi Created using ANSYS AIM 18.1 Problem Specification Pre-Analysis & Start Up Pre-Analysis Start-Up Geometry Import Geometry Mesh
More informationMaurice Holt University of California at Berkeley, Berkeley, California
NASA/CR-1998-208958 ICASE Report No. 98-54 3D Characteristics Maurice Holt University of California at Berkeley, Berkeley, California Institute for Computer Applications in Science and Engineering NASA
More informationComputational simulation of an altitude adaptive nozzle concept K. Schomberg 1, a, G. Doig 1, b and J. Olsen 1, c
Applied Mechanics and Materials Vol. 553 (2014) pp 223-228 (2014) Trans Tech Publications, Switzerland doi:10.4028/www.scientific.net/amm.553.223 Computational simulation of an altitude adaptive nozzle
More informationUNCERTAINTIES FOR A TRANSONIC DIFFUSER PROBLEM. William H. Mason, and Layne T. Watson
QUANTITATIVE RELATIVE COMPARISON OF CFD SIMULATION UNCERTAINTIES FOR A TRANSONIC DIFFUSER PROBLEM Serhat Hosder, Bernard Grossman, Raphael T. Haftka, William H. Mason, and Layne T. Watson Multidisciplinary
More information