High-Fidelity Multidisciplinary Design Using an Integrated Design Environment
|
|
- Gervase Small
- 5 years ago
- Views:
Transcription
1 High-Fidelity Multidisciplinary Design Using an Integrated Design Environment Antony Jameson and Juan J. Alonso Department of Aeronautics and Astronautics Stanford University, Stanford CA AFOSR Joint Contractors Meeting for Applied Analysis and Computational Mathematics Programs June
2 Outline Aero-Structural Wing Shape and Planform Optimization Fast Time Integration Methods for Unsteady Problems : application to shape optimization of a pitching airfoil Filtering the Navier Stokes Equations with an Invertible Filter: implication for LES subgrid modeling 2
3 Aero-Structural Wing Shape and Planform Optimization Shape Optimization Baseline profile (red) Redesigned profile (blue) Weakened Shock Boeing 747 -Planform Optimization Baseline planform (green) Redesigned planform (blue) Unstructured Mesh Optimization 3
4 Background: Levels of CFD HIGHEST Automatic Design Integrate the predictive capability into an automatic design method that incorporates computer optimization. Interactive Calculation Attainable when flow calculation can be performed fast enough But does NOT provide any guidance on how to change the shape if performance is unsatisfactory. LOWEST Flow Prediction Predict the flow past an airplane or its important components in different flight regimes such as take-off or cruise and off-design conditions such as flutter. Substantial progress has been made during the last decade. 4
5 Optimization and Design using Sensitivities Calculated by the Finite Difference Method The simplest approach is to define the geometry as f (x) = Âa b i i (x) where a i = weight, b i (x) = set of shape functions Then using the finite difference method, a cost function has sensitivities If the shape changes is I = I(w,a) (such as C D at constant C L ) I ª I(a + da ) - I(a ) i i i a i da i f(x) a n +1 = a n - l I a i (with small positive l) The resulting improvements is I + di = I - IT IT I da = I - l a a a < I More sophisticated search may be used, such as quasi - Newton. 5
6 Disadvantage of the Finite Difference Method The need for a number of flow calculations proportional to the number of design variables Using 4224 mesh points on the wing as design variables 4225 flow calculations ~ 30 minutes each (RANS) Boeing 747 Too Expensive 6
7 Application of Control Theory GOAL : Drastic Reduction of the Computational Costs Drag Minimization Optimal Control of Flow Equations subject to Shape(wing) Variations Define the cost function I = I(w,F) (for example C D at fixed C L ) and a change in F results in a change È di = I Î Í w T È dw + I Î Í F Suppose that the governing equation R which expresses the dependencd of w and F as R(w,F) = 0 and È dr = R È R Î Í w dw + Î Í F df = 0 T df 7
8 4224 design variables Application of Control Theory Since the variation dr is zero, it can be multiplied by a Lagrange Multiplier y and subtracted from the variation di without changing the result. di = IT w Ï dw + IT F = Ì IT w -yt Ó Ê È R È R df -yt Á Î Í w dw + Î Í F Ë df ˆ È R Ï IT Î Í w dw + Ì Ó Choosing y to satisfy the adjoint equation È R Î Í w T y = I w the first term is eliminated, and we find that where Ï G T = Ì I Ó F di = G T df T F -yt È R -y T Î Í F È R Î Í F df (Adjoint Equation) (Gradient) One Flow Solution + One Adjoint Solution 8
9 Advantages of the Adjoint Method: Gradient for N design variables with cost equivalent to two flow solutions Minimal memory requirement in comparison with automatic differentiation Enables shapes to be designed as free surface No need for user defined shape function No restriction on the design space 4224 design variables 9
10 Outline of the Design Process Flow solution Adjoint solution Gradient calculation Sobolev gradient Shape & Grid Modification Repeated until Convergence to Optimum Shape 10
11 Discrete versus Continuous Adjoint Methods The discrete adjoint method evaluates the adjoint and gradient equations algebraically from the discretized flow equations. The continuous adjoint method evaluates the costate solution from the partial differential adjoint equation. The continuous adjoint method leads to no inconsistency as long as it is combined with a compatible search method In the limit of grid convergence the two approaches yield identical gradients. Numerical tests of a model problem verify slightly superior accuracy with the continuous formulation ( Jameson and Vassberg 2000) 11
12 Summary of the Continuous Flow and Adjoint Equations With computational coordinates x i Euler equations for the flow : (1) x i S ij f j (w) = 0 where S ij are metrices, f j (w) the fluxes. Adjoint equation y f (2) C i = 0, C i = S j ij x i w Boundary condition for the Inverse problem (3) I = 1 2 Ú (p - p t ) 2 ds y 2 n x +y 3 n y +y 3 n z = p - p t Gradient y T (4) di = - Ú ds ij f j dd - Ú Ú ( ds D 21 y 2 + ds 22 y 3 + ds 23 y 4 )pdx 1 dx x b 3 w i 12
13 Point-wise Gradient and Shape Parameterizations If the shape is parameterized as then di = f (x) = Ú Â = da k k   k a k b k (x) g(x)df (x)dx Ú = G k da k where G k = k Ú g(x)b k (x)dx g(x)b k (x)dx 13
14 Sobolev Gradient Key issue for successful implementation of the Continuous adjoint method. Define the gradient with respect to the Sobolev inner product di = < g,df > = Ú ( gdf + eg'df ')dx Set df = - lg, di = - l < g,g > This approximates a continuous descent process df dt = -g The Sobolev gradient g is obtained from the simple gradient g by the smoothing equation g - x e g x = g. Continuous descent path 14
15 Computational Costs with N Design Variables Cost of Search Algorithm Steepest Descent Quasi-Newton Sobolev Gradient (Note: K is independent of N) Total Computational Cost of Design Finite Difference Gradients + Steepest Descent Finite Difference Gradients + Quasi-Newton Search or Response surface Adjoint Gradient + Quasi-Newton Search Adjoint Gradient + Sobolev Gradient (Note: K is independent of N) O(N 2 ) O(N ) O(K ) O(N 3 ) O(N 2 ) O(N ) O(K ) - N~ Big Savings - Enables Calculations on a Laptop 15
16 Trailing Edge Crossover Leading term in gradient for drag reduction may lead to trailing edge crossover crossover This corresponds to a sink in a free stream and hence negative drag. This is prevented by Sobolev gradient or by shape parameterizations which don t allow crossover. 16
17 Example: Viscous RAE Drag Minimization Initial Shape and Non-Smoothed Gradient Initial Shape and Smoothed Gradient 17
18 Example: Viscous RAE Drag Minimization Final Shape and Non-Smoothed Gradient Final Shape and Smoothed Gradient 18
19 Good Shape Parameterization of an Airfoil via Conformal Mapping C Define the mapping by On C z log ds dq + i(a -q - p 2 ) = n  C ke -ikq k= 0 n log dz ds = C  k s k =  ( a k cos(kq) + b k sin(kq) ) + iâ b k cos(kq) + a k sin(kq) k= 0 a(s) n k= 0 ( ) Fourier coefficients define the mapping Spectral convergence with n Trialing edge gap is 2piC 1, so prevent crossover by setting C 1 = 0. s 19
20 Sobolev Gradient for Partial Redesign One may wish to freeze some fo the profile, e.g. the structure box Fixed (Structure box) In this case shape changes needed to be smoothly blended into the frozen geometry. This is accomplished by adding more derivatives to the Sobolev gradient ( ) Then the corresponding smoothed gradient satisfies g - Ê x e g ˆ Á Ê 2 g ˆ e Ë x x 2 Á 4 Ë x 2 = g, Hence a patch can be smoothly blended. Ú < u,v >= uv + e 2 u v + e 4 u v dx g = g x = 0 at both end points. 20
21 Potential Modification of F8 Bearcat Racer to win the RENO Air Racer 21
22 Bearcat Shock Free Result Initial C p Final C p 22
23 RANS Calculations Redesign of the Boeing 747 Wing at its Cruise Mach Number Constraints : Fixed C L = 0.42 : Fixed span-load distribution : Fixed thickness 10% wing drag saving (3 hrs cpu time - 16proc.) ~5% aircraft drag saving baseline redesign 23
24 Redesign of the Boeing 747 Wing at Mach 0.9 Sonic Cruiser Constraints : Fixed C L = 0.42 : Fixed span-load distribution RANS Calculations : Fixed thickness Same C We can fly faster at the same drag. 24
25 Redesign of the Boeing 747: Drag Rise ( Three-Point Design ) Constraints : Fixed C L = 0.42 : Fixed span-load distribution : Fixed thickness RANS Calculations Improved Wing L/D Improved M DD benefit benefit Lower drag at the same Mach Number Fly faster with the same drag 25
26 Planform and Aero-Structural Optimization Design tradeoffs suggest an multidisciplinary design and optimization Range = VL D 1 sfc logw o + W f W o Maximize Minimize Planform variations can further maximize VL/D but affects W O 26
27 Aerodynamic Design Tradeoffs The drag coefficient can be split into C D = C DO + C 2 L pear L is maximized if the two terms are equal. D Induced drag is half of the total drag. If we want to have large drag reduction, we should target the induced drag. Design dilemma Increase b D i = 2L 2 perv 2 b 2 D i decreases W O increases Change span by changing planform 27
28 Break Down of Drag Boeing 747 at C L ~.47 (including fuselage lift ~ 15%) Item Wing Pressure Wing friction Fuselage Tail Nacelles Other Total C D 120 counts (15 shock, 105 induced) Cumulative C D 120 counts Induced Drag is the largest component 28
29 Wing Planform Optimization Simplified Planform Model Wing planform modification can yield larger improvements BUT affects structural weight. 1 I = a 1 C D + a 2 2 where C W = Ú (p - p d ) 2 ds Structural Weight q S ref + a 3 C W Can be thought of as constraints 29
30 Additional Features Needed Structural Weight Estimation Large scale gradient : span, sweep, etc Adjoint gradient formulation for dc w /dx Choice of a 1, a 2, and a 3 Use fully-stressed wing box to estimate the structural weight. Large scale gradient Use summation of mapped gradients to be large scale gradient 30
31 Choice of Weighting Constants Breguet range equation R = VL 1 D sfc logw O + W f W O With fixed V, L, sfc, and (W O + W f W TO ), the variation of R can be stated as Ê dr Á R = - Á dc D + C Á D Ë 1 log W TO W O dw O W O ˆ Ê Á = -Á dc D + Á C D Á Ë 1 log C W TO C WO dc WO C WO ˆ Maximizing Range Minimizing I = C D + a 3 a 1 C W using a 3 a 1 = C D C WO log C W TO C W0 31
32 Planform Optimization of Boeing 747 Constraints : Fixed CL=0.42 Redesign Baseline 1) Longer span reduces the induced drag 2) Less sweep and thicker wing sections reduces structure weight 3) Section modification keeps shock drag minimum Overall: Drag and Weight Savings C D C W Baseline Optimize Section at Fixed planform Optimize both section and planform
33 Planform Optimization of MD11 Constraints : Fixed CL=0.45 Baseline Redesign C D C W Baseline Optimize Section at Fixed planform Optimize both section and planform
34 Pareto Front: Expanding the Range of Designs Use multiple a 3 /a 1 ==> Multiple Optimal Shapes Boundary of realizable designs Pareto front of Boeing
35 Super 747 Design a new wing for the Boeing 747 Strategy Use the same fuselage of Boeing 747 Use a new planform (from the Planform Optimization result) Use new airfoil sections (AJ airfoils) Optimized for fixed lift coefficient 35
36 Baseline Boeing
37 Super
38 Comparison C L C D counts C W counts Boeing (112 pressure, 34.3 viscous) (82,550 lbs) Super (83.8 pressure, 37.2 viscous) (70,620 lbs) Save ~10% of airplane drag Save ~ 7% of airplane structural weight 38
39 Automatic design for the Complete Aircraft Geometry on an Unstructured Mesh (SYNPLANE) Key step: reduce the gradient to a surface integral independent of the mesh perturbation (Jameson, A., and Kim, S., "Reduction of the Adjoint Gradient Formula in the Continuous Limit", 41 st AIAA Aerospace Sciences Meeting & Exhibit, AIAA Paper , Reno, NV, January 6-9, ) di = Ú Ú y T ( ds 2 j f j + C 2 dw * )dx 1 dx 3 - Ú Ú ( ds 21 y 2 + ds 22 y 3 + ds 23 y 4 )pdx 1 dx 3 b w Compared to the previous formulation y T di = - Ú ds ij f j dd - Ú Ú ( ds D 21 y 2 + ds 22 y 3 + ds 23 y 4 )pdx 1 dx x b 3 w i b w This field integral is converted to boundary integral 39
40 Question : If the entire profile or wing is translated as a rigid body, the flow is unchanged. Therefore if one calculates the point-wise gradients g xi for movement in the x i direction, they should satisfy the relation Ú g xi ds = 0 Is this true for the computed gradients? 40
41 Actually it can be verified for some implementations of the field integral formula for the gradient. But the boundary integral formula presents a difficulty in the neighborhood of the front and rear stagnation points where w is not well defined. 41
42 Validation of the Reduced Adjoint RAE 2822 airfoil,grid 128x32, Euler calculation Translate the mesh as a rigid body (Far-field boundary is not changed.) Expect summed gradient around the airfoil to be zero. Initial Grid x Translated Grid (1 chord length in x direction) 42
43 Validation of the Reduced Adjoint Gradient Comparison of Original Adjoint, Reduced Adjoint, and Complex-Step methods Comparison of Original Adjoint and Reduced Adjoint methods 43
44 Gradient Comparison Point-wise gradient Integral of the point-wise gradient over the airfoil surface FD 1E-6 Grad 1E-7 Theoretically there is no change in the flow, the integral of the gradient should become zero. 44
45 Deformation of Unstructured Meshes Movement of an individual surface point influences entire mesh Method 1) Spring method: treat edges as springs. Computationally inexpensive but doesn t guarantee to prevent crossovers. 2) Treat the mesh as a pseudo elastic body. Solve x j s ij = 0, s ij = u i x j + u j x i s = stress, u = displacement for a given boundary displacement. 3) Traction method : treat the mesh as an elastic body, but subject to force inputs proportional to the gradient, instead of displacements - alternative to the Sobolev-gradient. 45
46 Redesign of Falcon Complete aircraft calculation on Unstructured Mesh Shock C D = 234 counts 46
47 Redesign of Falcon Using SYNPLANE Drag reduction 18 counts at fixed C L = 0.4 Weakened Shock C D = 216 counts 47
48 Fast Time Integration Methods for Unsteady Problems 48
49 Potential Applications Flutter Analysis, Flow past Helicopter blades, Rotor-Stator Combinations in Turbomachinery, Zero-Mass Synthetic Jets for Flow Control 49
50 Dual Time Stepping BDF The kth-order accurate backward difference formula (BDF) is of the form D t = 1 Dt The non-linear BDF is solved by inner iterations which advance in pseudo-time t* The second-order BDF solves Implementation via Dt * k 1 Â (D - ) q where D - w n +1 = w n +1 - w n q q=1 dw dt + È 3w - 4wn + w n-1 + R(w) * Í = 0 Î 2Dt RK dual time stepping scheme with variable local (RK-BDF) Nonlinear SGS dual time stepping scheme (SGS-BDF) with Multigrid 50
51 Test Case: NACA64A010 pitching airfoil (CT6 Case) Mach Number Pitching amplitude Reduced Freq. Reynolds Number /- 1.01deg million Cycling to limit cycle Pressure Contours at Various Time Instances (AGARD 702) Results of SGS-BDF Scheme (36 time steps per pitching cycle, 3 iterations per time step ) 51
52 Payoff of Dual-time Stepping BDF Schemes Accurate simulations with an order of magnitude reduction in time steps. order For the pitching airfoil: from ~ 1000 to 36 time steps per pitching cycle with three sub-iterations in each step. 52
53 Frequency Domain and Global Space-Time Multigrid Spectral Methods Application : Time-periodic flows Using a Fourier representation in time, the time period T is divided into N steps. ˆ w k = 1 N N-1 w n e -ikndt  n= 0 Then, the discretization operator is given by D t w n = N 2-1  ikw ˆ k e ikndt k= -N 2 53
54 Method 1 (McMullen et.al.) : Transform the equations into frequency domain and solve them in pseudo-time t* d ˆ w k dt * + ik w ˆ k + R ˆ k = 0 Method 2 (Gopinath et.al.) : Solve the equations in the timedomain. The space-time spectral discretization operator is N 2-1 D t w n = Â d m w n +m, d m = 1 2 (-1)m +1 cot( pm N ),m 0 m=- N 2 +1 This is a central difference operator connecting all time levels, yielding an integrated space-time formulation which requires simultaneous solution of the equations at all time levels. 54
55 Comparison with Experimental Data - C L vs. a RANS Time-Spectral Solution with 4 and 8 intervals per pitching cycle Computed Results - Time spectral 4 time intervals - Time spectral 8 time intervals * AGARD-702:Davis Experimental Data 55
56 Payoff of Time Spectral Schemes Engineering accuracy with very small number of time intervals and same rate of convergence as the BDF. Spectral accuracy for sufficiently smooth solutions. Periodic solutions directly without the need to evolve through 5-10 cycles, yielding an order of magnitude reduction in computing cost beyond the reduction already achieved with the BDF, for a total of two orders of magnitude. 56
57 Example of Shape Optimization in Periodic Unsteady Flow (Nadarajah, McMullen, and Jameson AIAA ) Find the (fixed) shape which minimizes the averaged drag coefficient of a pitching airfoil with the constraint that the average lift coefficient is maintained. 57
58 Optimization Results Non-Linear Frequency Domain result Matched Time Accurate result Initial Shape Final Shape Comparison of the Final Airfoil Geometry between the Time Accurate and NLFD method with Lift Constraints. Comparison of the Initial and Final Airfoil Geometry using the NLFD method with Lift Constraints. 58
59 Pressure Distributions Initial Shape Final Shape Comparison of the Final Pressure Distribution between the Time Accurate and NLFD method with Lift Constraints. Comparison of the Initial and Final Pressure Distribution using the NLFD method with Lift Constraints. 59
60 Convergence History Time Accurate Vs. Non-Linear Frequency Domain method 60
61 Flow Past Helicopter Blades 61
62 Challenges Blade runs into its own wake. Resolve blade-vortex interaction Simulate fully articulated hub that is capable of lead-lag, flapping motions coupled with aeroelastic solver Reduce simulation time through convergence acceleration techniques such as space-time multigrid spectral methods 62
63 Test Case: Non-Lifting Hover 2-bladed NACA0012 Mesh size: Tip Mach number 0º Collective pitch Only one section of the blade needs to be calculated because of periodicity Euler Calculation using BDF scheme 63
64 Preliminary Results using the Backward Different Formula Plots of -Cp at different span stations Results look reasonable but needs refinement Subtle modification to the mesh should yield better agreement with experiments 64
65 Future Work Apply space-time multigrid spectral methods This will enable solving Navier-Stokes calculations in reasonable time Couple the flow solution with an aeroelastic model Simulation of forward flight that will allow blade flapping, lagging and change of cyclic pitch Eventually, introduce aero-structural shape optimization to improve rotor performance 65
66 Filtering the Navier Stokes Equations with an Invertible Filter 66
67 Consider the incompressible Navier--Stokes equations where r u i t + ru j u i + r p = m 2 u i (1) x j x i x i x j u i x i = 0 In large eddy simulation (LES) the solution is filtered to remove the small scales. Typically one sets u i (x) = Ú G(x - x )u( x )dx (2) where the kernel G is concentrated in a band defined by the filter width. Then the filtered equations contain the extra virtual stress t ij = u i u j - ui u j (3) because the filtered value of a product is not equal to the product of the filtered values. This stress has to be modeled. 67
68 A filter which completely cuts off the small scales or the high frequency components is not invertible. The use, on the other hand, of an invertible filter would allow equation (1) to be directly expressed in terms of the filtered quantities. Thus one can identify desirable properties of a filter as 1. Attenuation of small scales 2. Commutativity with the differential operator 3. Invertibility Suppose the filter has the form ui = Pu i (4) which can be inverted as Qui = u i (5) where Q = P -1. Moreover Q should be coercieve, so that Qu > c u (6) for some positive constant c. 68
69 Note that if Q commutes with x i then so does Q -1, since for any quantity f which is sufficiently differentiable Also since Q commutes with (Q -1 f ) = Q -1 Q (Q -1 f ) x i x i x i, = Q -1 x i (QQ -1 f ) = Q -1 x i ( f ) u x i = 0 (7) As an example P can be the inverse Helmholtz operator, so that one can write Ê Qui = 1-a 2 2 ˆ Á ui = u i (8) Ë x k x k where a is a length scale proportional to the largest scales to be retained. One may also introduce a filtered pressure p, satisfying the equation Ê Qp = 1-a 2 2 ˆ Á p = p (9) Ë x k x k 69
70 Now one can substitute equation (8) and(9) for u i and p in equation (1) to get r Ê 2 ˆ Ê Á t 1-a2 ui + r 1-a 2 2 ˆ Á u j Ë x k x k Ë x k x k = m 2 Ê 1-a 2 2 ˆ Á x j x j Ë x k x k Ê 1-a 2 2 ˆ Á ui + Ê 1-a 2 2 ˆ Á p x j Ë x l x l x i Ë x k x k Because the order of the differentiations can be interchanged and the Helmholtz operator satisfies condition(6), it can be removed. The product term can be written as r Ï Ê 1-a 2 2 ˆ Ê Á ui 1-a 2 2 ˆ Ì Á u j x j Ó Ë x k x k Ë x l x l = r x j Ï Ì uiu j -a 2 ui Ó = r Ï 2 Ì uiu j -a 2 x j Ó x k x k = rq x j 2 u j x k x k -a 2 u j ui 2 ui x k x k + a 4 2 ui x k x k 2 u j x l x l ( ui u j ) + 2a 2 ui u j 2 + a 4 ui x k x k x k x k Ï Ê uiu j + a 2 Q -1 2 ui u j 2 + a 2 ui 2 u ˆ j Ì Á Ó Ë x k x k x k x k x l x l 2 u j x l x l According to condition (6), if Qf = 0 for any sufficiently differentiable quantity f, then f = 0. 70
71 Thus the filtered equation finally reduces to with the virtual stress r u i t + r x j ( ui u j ) + p = m 2 ui - r t ij (10) x i x k x k x j Ê t ij = a 2 Q -1 2 ui u j + a 2 2 ui 2 u ˆ j Á (11) Ë x k x k x k x k x l x l The virtual stress may be calculated by solving Ê 1-a 2 2 ˆ Ê Á t ij = a 2 2 ui u j + a 2 2 ui Á Ë x k x k Ë x k x k x k x k 2 u ˆ j (12) x l x l Taking the divergence of equation (10), it also follows that p satisfies the Poisson equation 2 p + r x i x i x i x j ui u j ( ) + r 2 x i x j t ij = 0 (13) 71
72 In a discrete solution scales smaller than the mesh width would not be resolved, amounting to an implicit cut off. There is the possibility of introducing an explicit cut off off in t ij. Also one could use equation (8) to restore an estimate of the unfiltered velocity. In order to avoid solving the Helmholtz equation (12), the inverse Helmholtz operator could be expanded formally as ( 1- a 2 D) -1 = 1+a 2 D +a 4 D where D denotes the Laplacian the approximate virtual stress tensor assumes the form 2 x k x k.now retaining terms up to the fourth power of a, u t ij = 2a 2 i u È Ê j +a 4 2D u i u j Í Á x k x k Î Í Ë x k x k ˆ + Du idu j One ay regard the forms (11) or (14) as prototypes for subgrid scale (SGS) models. (14) 72
73 The inverse Helmholtz operator cuts off the smaller scales quite gradually. One could design filters with a sharper cut off by shaping their frequency response. Denote the Fourier transform of f as where (in one space dimension) ˆ f (k) = f (k) = ˆ f = Ff 1 2p 1 2p Then the general form of an invertible filter is F Pf Ÿ f (x)e -ikx dx Ú - f ˆ (k)e -ikx dk Ú - = S(k) ˆ f (k) Ÿ 1 F Qf = S(k) ˆ f (k) where S(k) should decrease rapidly beyond a cut off wave number inversely proportional to a length scale a. 73
74 In the case of a general filter with inverse Q, the virtual stress follows from the relation Then Qu i u j = u i u j = Qu i Qu j t ij = u i u j - u i u j = Q -1 (Qu i Qu j - Q(u i u j )) This formula provides the form for a family of subgrid-scale models. 74
75 Quotation From John Vassberg You are better off with a design method which optimizes 5000 design variables with 10 function evaluations than one which requires 5000 function evaluations to optimize 10 design variables. 75
Solution Algorithms for Viscous Flow
Solution Algorithms for Viscous Flow Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA Indian Institute of Science Bangalore, India September 20, 2004 Flo3xx Computational
More informationChallenges and Complexity of Aerodynamic Wing Design
Chapter 1 Challenges and Complexity of Aerodynamic Wing Design Kasidit Leoviriyakit and Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA kasidit@stanford.edu and
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics and Astronautics AIAA Aerospace Sciences Meeting, Reno, NV AIAA Paper 2007-0037
More informationSolution Algorithms for High-Speed Viscous Flows Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA
Solution Algorithms for High-Speed Viscous Flows Antony Jameson Department of Aeronautics and Astronautics Stanford University, Stanford CA 13 th Conference on Finite Element for Flow Problems University
More informationJun 22-25, 2009/San Antonio, TX
19th AIAA Computational Fluid Dynamics 22-25 June 2009, San Antonio, Texas AIAA 2009-4273 AIAA 2009 4273 An Assessment of Dual-Time Stepping, Time Spectral and Artificial Compressibility based Numerical
More informationApplication of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations
Application of Dual Time Stepping to Fully Implicit Runge Kutta Schemes for Unsteady Flow Calculations Antony Jameson Department of Aeronautics and Astronautics, Stanford University, Stanford, CA, 94305
More informationApplications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow
Applications of adjoint based shape optimization to the design of low drag airplane wings, including wings to support natural laminar flow Antony Jameson and Kui Ou Aeronautics & Astronautics Department,
More informationA Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization
A Crash-Course on the Adjoint Method for Aerodynamic Shape Optimization Juan J. Alonso Department of Aeronautics & Astronautics Stanford University jjalonso@stanford.edu Lecture 19 AA200b Applied Aerodynamics
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Matthew McMullen and Antony Jameson and Juan J. Alonso Dept. of Aeronautics & Astronautics Stanford University
More informationNonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations p.1/50
Nonlinear Frequency Domain Methods Applied to the Euler and Navier-Stokes Equations Matthew McMullen Advisor: Antony Jameson Co-advisor: Juan Alonso Sponsor: Accelerated Strategic Computing Initiative
More informationAn Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond
45 th Aerospace Sciences Meeting and Exhibit, January 8 11, 2007, Reno, Nevada An Investigation of the Attainable Efficiency of Flight at Mach One or Just Beyond Antony Jameson Department of Aeronautics
More informationCOMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD
COMPLETE CONFIGURATION AERO-STRUCTURAL OPTIMIZATION USING A COUPLED SENSITIVITY ANALYSIS METHOD Joaquim R. R. A. Martins Juan J. Alonso James J. Reuther Department of Aeronautics and Astronautics Stanford
More informationAERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS
AERO-STRUCTURAL WING DESIGN OPTIMIZATION USING HIGH-FIDELITY SENSITIVITY ANALYSIS Joaquim R. R. A. Martins and Juan J. Alonso Department of Aeronautics and Astronautics Stanford University, Stanford, CA
More informationApplication of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations
Application of a Non-Linear Frequency Domain Solver to the Euler and Navier-Stokes Equations Contours of Pressure.8.6.4...4.6.8 Matthew McMullen Antony Jameson Juan J. Alonso Aero Astro Dept. Stanford
More informationA COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION
A COUPLED-ADJOINT METHOD FOR HIGH-FIDELITY AERO-STRUCTURAL OPTIMIZATION Joaquim Rafael Rost Ávila Martins Department of Aeronautics and Astronautics Stanford University Ph.D. Oral Examination, Stanford
More informationActive Flutter Control using an Adjoint Method
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-844 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Active Flutter Control using an
More informationTurbomachinery Applications with the Time Spectral Method
35th AIAA Fluid Dynamics Conference and Exhibit, June 6-9,005, Toronto,Ontario Turbomachinery Applications with the Time Spectral Method Edwin van der Weide, Arathi K Gopinath and Antony Jameson Stanford
More informationHigh Speed Aerodynamics. Copyright 2009 Narayanan Komerath
Welcome to High Speed Aerodynamics 1 Lift, drag and pitching moment? Linearized Potential Flow Transformations Compressible Boundary Layer WHAT IS HIGH SPEED AERODYNAMICS? Airfoil section? Thin airfoil
More informationIndustrial Applications of Aerodynamic Shape Optimization
Industrial Applications of Aerodynamic Shape Optimization John C. Vassberg Boeing Technical Fellow Advanced Concepts Design Center Boeing Commercial Airplanes Long Beach, CA 90846, USA Antony Jameson T.
More informationSimulation of Aeroelastic System with Aerodynamic Nonlinearity
Simulation of Aeroelastic System with Aerodynamic Nonlinearity Muhamad Khairil Hafizi Mohd Zorkipli School of Aerospace Engineering, Universiti Sains Malaysia, Penang, MALAYSIA Norizham Abdul Razak School
More informationDrag Computation (1)
Drag Computation (1) Why drag so concerned Its effects on aircraft performances On the Concorde, one count drag increase ( C D =.0001) requires two passengers, out of the 90 ~ 100 passenger capacity, be
More informationFeedback Control of Aerodynamic Flows
44th AIAA Aerospace Sciences Meeting and Exhibit 9-12 January 26, Reno, Nevada AIAA 26-843 44th AIAA Aerospace Sciences Meeting and Exhibit, Reno, Nevada, 9 12 Jan, 26. Feedback Control of Aerodynamic
More informationAIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION
AIRFRAME NOISE MODELING APPROPRIATE FOR MULTIDISCIPLINARY DESIGN AND OPTIMIZATION AIAA-2004-0689 Serhat Hosder, Joseph A. Schetz, Bernard Grossman and William H. Mason Virginia Tech Work sponsored by NASA
More informationOptimum aerodynamic design using the Navier-Stokes equations
Copyright 1997, American Institute of Aeronautics and Astronautics, Inc. AIAA Meeting Papers on Disc, January 1997 A9715189, AF-AFOSR-91-0391, AIAA Paper 97-0101 Optimum aerodynamic design using the Navier-Stokes
More informationTHE current generation of civilian transport aircraft are typically
JOURNAL OF AIRCRAFT Vol. 47, No. 2, March April 2010 Aerodynamic Structural Design Studies of Low-Sweep Transonic Wings Antony Jameson Stanford University, Stanford, California 94305-3030 John C. Vassberg
More informationAeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior
Aeroelastic Analysis of Engine Nacelle Strake Considering Geometric Nonlinear Behavior N. Manoj Abstract The aeroelastic behavior of engine nacelle strake when subjected to unsteady aerodynamic flows is
More informationModeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique
Modeling Unsteady Flow in Turbomachinery Using a Harmonic Balance Technique Torsten Palenschat 317220 30.04.2014 Master seminar CES WS 2013/2014 Center for Computational Engineering Sciences RWTH Aachen
More informationComparative Analysis of Computational Methods for Limit-Cycle Oscillations
47th AIAA Structures, Structural Dynamics and Materials Conference, May 1-4, 26 Newport, Rhode Island Comparative Analysis of Computational Methods for Limit-Cycle Oscillations Arathi K. Gopinath, Stanford
More informationThe Truth About Elliptic Spanloads or Optimum Spanloads Incorporating Wing Structural Weight
The Truth About Elliptic Spanloads or Optimum Spanloads Incorporating Wing Structural Weight Sergio Iglesias and William H. Mason AIAA Paper 2001-5234 as presented at the 1st AIAA Aircraft Technology,
More informationCOMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL
COMPUTATIONAL STUDY OF SEPARATION CONTROL MECHANISM WITH THE IMAGINARY BODY FORCE ADDED TO THE FLOWS OVER AN AIRFOIL Kengo Asada 1 and Kozo Fujii 2 ABSTRACT The effects of body force distribution on the
More informationA HARMONIC BALANCE APPROACH FOR MODELING THREE-DIMENSIONAL NONLINEAR UNSTEADY AERODYNAMICS AND AEROELASTICITY
' - ' Proceedings of ASME International Mechanical Engineering Conference and Exposition November 17-22, 22, New Orleans, Louisiana, USA IMECE-22-3232 A HARMONIC ALANCE APPROACH FOR MODELING THREE-DIMENSIONAL
More information41st AIAA Aerospace Sciences Meeting and Exhibit January 6 9, 2003/Reno, NV
AIAA 3 267 Calculation of Unsteady Transonic Flow by an Euler Method with Small Perturbation Boundary Conditions C. Gao, S. Luo, F. Liu Department of Mechanical and Aerospace Engineering University of
More informationAirfoils and Wings. Eugene M. Cliff
Airfoils and Wings Eugene M. Cliff 1 Introduction The primary purpose of these notes is to supplement the text material related to aerodynamic forces. We are mainly interested in the forces on wings and
More informationMULTIGRID CALCULATIONS FOB. CASCADES. Antony Jameson and Feng Liu Princeton University, Princeton, NJ 08544
MULTIGRID CALCULATIONS FOB. CASCADES Antony Jameson and Feng Liu Princeton University, Princeton, NJ 0544 1. Introduction Development of numerical methods for internal flows such as the flow in gas turbines
More informationIs My CFD Mesh Adequate? A Quantitative Answer
Is My CFD Mesh Adequate? A Quantitative Answer Krzysztof J. Fidkowski Gas Dynamics Research Colloqium Aerospace Engineering Department University of Michigan January 26, 2011 K.J. Fidkowski (UM) GDRC 2011
More informationAerodynamic Inverse Design and Shape Optimization via Control Theory
Aerodynamic Inverse Design and Shape Optimization via Control Theory 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics & Astronautics Stanford University SciTech 2015 January 5, 2015
More informationAeroelastic Gust Response
Aeroelastic Gust Response Civil Transport Aircraft - xxx Presented By: Fausto Gill Di Vincenzo 04-06-2012 What is Aeroelasticity? Aeroelasticity studies the effect of aerodynamic loads on flexible structures,
More informationAProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy
AProofoftheStabilityoftheSpectral Difference Method For All Orders of Accuracy Antony Jameson 1 1 Thomas V. Jones Professor of Engineering Department of Aeronautics and Astronautics Stanford University
More informationInvestigation potential flow about swept back wing using panel method
INTERNATIONAL JOURNAL OF ENERGY AND ENVIRONMENT Volume 7, Issue 4, 2016 pp.317-326 Journal homepage: www.ijee.ieefoundation.org Investigation potential flow about swept back wing using panel method Wakkas
More informationTransonic Flutter Prediction of Supersonic Jet Trainer with Various External Store Configurations
Transonic Flutter Prediction of Supersonic Jet Trainer with Various External Store Configurations In Lee * Korea Advanced Institute of Science and Technology, Daejeon, 305-701, Korea Hyuk-Jun Kwon Agency
More informationCOMPUTATION OF CASCADE FLUTTER WITH A COUPLED AERODYNAMIC AND STRUCTURAL MODEL
COMPUTATION OF CASCADE FLUTTER WITH A COUPLED AERODYNAMIC AND STRUCTURAL MODEL M. SADEGHI AND F. LIU Department of Mechanical and Aerospace Engineering University of California, Irvine, CA 92697-3975,
More informationPaul Garabedian s Contributions to Transonic Airfoil and Wing Design
Paul Garabedian s Contributions to Transonic Airfoil and Wing Design Antony Jameson October 13, 010 Abstract This note on Paul Garabedian s work on transonic airfoil and wing design is written from the
More informationOptimization of Flapping Airfoils for Maximum Thrust and Propulsive Efficiency I. H. Tuncer, M. Kay
Czech Technical University in Prague Acta Polytechnica Vol. 44 No. 1/2004 Optimization of Flapping Airfoils for Maximum Thrust and Propulsive Efficiency I. H. Tuncer, M. Kay A numerical optimization algorithm
More informationNumerical Optimization Algorithms
Numerical Optimization Algorithms 1. Overview. Calculus of Variations 3. Linearized Supersonic Flow 4. Steepest Descent 5. Smoothed Steepest Descent Overview 1 Two Main Categories of Optimization Algorithms
More informationDerivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint
Derivatives for Time-Spectral Computational Fluid Dynamics using an Automatic Differentiation Adjoint Charles A. Mader University of Toronto Institute for Aerospace Studies Toronto, Ontario, Canada Joaquim
More informationAirfoil shape optimization using adjoint method and automatic differentiation
Airfoil shape optimization using adjoint method and automatic differentiation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 5665
More informationAdjoint Formulations for Topology, Shape and Discrete Optimization
45 th Aerospace Sciences Meeting and Exhibit, January 8 11, 27, Reno, Nevada Adjoint Formulations for Topology, Shape and iscrete Optimization Sriram and Antony Jameson epartment of Aeronautics and Astronautics
More informationImplementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible Flapping Wings within Overture
10 th Symposimum on Overset Composite Grids and Solution Technology, NASA Ames Research Center Moffett Field, California, USA 1 Implementing a Partitioned Algorithm for Fluid-Structure Interaction of Flexible
More informationTHE calculation of unsteady flows continues to present a severe
AIAA JOURNAL Vol. 44, No. 7, July 2006 Demonstration of Nonlinear Frequency Domain Methods Matthew McMullen, Antony Jameson, and Juan Alonso Stanford University, Stanford, California 94305 USA DOI: 10.2514/1.15127
More informationNumerical Solution of Partial Differential Equations governing compressible flows
Numerical Solution of Partial Differential Equations governing compressible flows Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore
More informationToward Practical Aerodynamic Design Through Numerical Optimization
Toward Practical Aerodynamic Design Through Numerical Optimization David W. Zingg, and Laura Billing, Institute for Aerospace Studies, University of Toronto 4925 Dufferin St., Toronto, Ontario M3H 5T6,
More informationDrag (2) Induced Drag Friction Drag Form Drag Wave Drag
Drag () Induced Drag Friction Drag Form Drag Wave Drag Outline Nomenclature and Concepts Farfield Drag Analysis Induced Drag Multiple Lifting Surfaces Zero Lift Drag :Friction and Form Drag Supersonic
More informationDrag Characteristics of a Low-Drag Low-Boom Supersonic Formation Flying Concept
Drag Characteristics of a Low-Drag Low-Boom Supersonic Formation Flying Concept Yuichiro Goto, Shigeru Obayashi and Yasuaki Kohama Tohoku University, Sendai, Japan In this paper, a new concept for low-drag,
More information39th AIAA Aerospace Sciences Meeting and Exhibit January 8 11, 2001/Reno, NV
AIAA 2 529 Unsteady Flow Investigations in an Axial Turbine Using the Massively Parallel Flow Solver TFLO Jixian Yao, Roger L. Davis, Juan J. Alonso, and Antony Jameson Stanford University, Stanford, CA
More informationThree-Dimensional Unsteady Multi-stage Turbomachinery Simulations using the Harmonic Balance Technique
Three-Dimensional Unsteady Multi-stage Turbomachinery Simulations using the Harmonic Balance Technique Arathi K. Gopinath, Edwin van der Weide, Juan J. Alonso, Antony Jameson Stanford University, Stanford,
More informationCHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION
CHAPTER 7 SEVERAL FORMS OF THE EQUATIONS OF MOTION 7.1 THE NAVIER-STOKES EQUATIONS Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationDEVELOPMENT OF A THREE-DIMENSIONAL TIGHTLY COUPLED EULER/POTENTIAL FLOW SOLVER FOR TRANSONIC FLOW
DEVELOPMENT OF A THREE-DIMENSIONAL TIGHTLY COUPLED EULER/POTENTIAL FLOW SOLVER FOR TRANSONIC FLOW Yeongmin Jo*, Se Hwan Park*, Duck-Joo Lee*, and Seongim Choi *Korea Advanced Institute of Science and Technology,
More informationAerodynamic Resonance in Transonic Airfoil Flow. J. Nitzsche, R. H. M. Giepman. Institute of Aeroelasticity, German Aerospace Center (DLR), Göttingen
Aerodynamic Resonance in Transonic Airfoil Flow J. Nitzsche, R. H. M. Giepman Institute of Aeroelasticity, German Aerospace Center (DLR), Göttingen Source: A. Šoda, PhD thesis, 2006 Slide 2/39 Introduction
More informationComplete Configuration Aero-Structural Optimization Using a Coupled Sensitivity Analysis Method
AIAA 2002 5402 Complete Configuration Aero-Structural Optimization Using a Coupled Sensitivity Analysis Method Joaquim R. R. A. Martins, Juan J. Alonso Stanford University, Stanford, CA 94305 James J.
More informationFurther Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow
29th AIAA Applied Aerodynamics Conference 27-30 June 2011, Honolulu, Hawaii AIAA 2011-3509 Further Studies of Airfoils Supporting Non-unique Solutions in Transonic Flow Antony Jameson, John C. Vassberg,
More informationHelicopter Rotor Unsteady Aerodynamics
Helicopter Rotor Unsteady Aerodynamics Charles O Neill April 26, 2002 Helicopters can t fly; they re just so ugly the earth repels them. 1 Introduction Unsteady helicopter rotor aerodynamics are complicated.
More informationA method for identification of non-linear multi-degree-of-freedom systems
87 A method for identification of non-linear multi-degree-of-freedom systems G Dimitriadis and J E Cooper School of Engineering, Aerospace Division, The University of Manchester Abstract: System identification
More informationFar Field Noise Minimization Using an Adjoint Approach
Far Field Noise Minimization Using an Adjoint Approach Markus P. Rumpfkeil and David W. Zingg University of Toronto Institute for Aerospace Studies 4925 Dufferin Street, Toronto, Ontario, M3H 5T6, Canada
More information( ) Notes. Fluid mechanics. Inviscid Euler model. Lagrangian viewpoint. " = " x,t,#, #
Notes Assignment 4 due today (when I check email tomorrow morning) Don t be afraid to make assumptions, approximate quantities, In particular, method for computing time step bound (look at max eigenvalue
More informationHigh-Fidelity Aero-Structural Design Optimization of a Supersonic Business Jet
AIAA 2002 1483 High-Fidelity Aero-Structural Design Optimization of a Supersonic Business Jet Joaquim R. R. A. Martins, Juan J. Alonso Stanford University, Stanford, CA 94305 James J. Reuther NASA Ames
More informationA Multi-Dimensional Limiter for Hybrid Grid
APCOM & ISCM 11-14 th December, 2013, Singapore A Multi-Dimensional Limiter for Hybrid Grid * H. W. Zheng ¹ 1 State Key Laboratory of High Temperature Gas Dynamics, Institute of Mechanics, Chinese Academy
More informationIntroduction to Aerospace Engineering
4. Basic Fluid (Aero) Dynamics Introduction to Aerospace Engineering Here, we will try and look at a few basic ideas from the complicated field of fluid dynamics. The general area includes studies of incompressible,
More informationAero-Propulsive-Elastic Modeling Using OpenVSP
Aero-Propulsive-Elastic Modeling Using OpenVSP August 8, 213 Kevin W. Reynolds Intelligent Systems Division, Code TI NASA Ames Research Center Our Introduction To OpenVSP Overview! Motivation and Background!
More informationNumerical Simulation of Flow Separation Control using Multiple DBD Plasma Actuators
Journal of Applied Fluid Mechanics, Vol. 9, No. 4, pp. 1865-1875, 2016. Available online at www.jafmonline.net, ISSN 1735-3572, EISSN 1735-3645. DOI: 10.18869/acadpub.jafm.68.235.25325 Numerical Simulation
More informationSENSITIVITY ANALYSIS OF THE FACTORS AFFECTING FORCE GENERATION BY WING FLAPPING MOTION
Proceedings of the ASME 2013 International Mechanical Engineering Congress and Exposition IMECE2013 November 15-21, 2013, San Diego, California, USA IMECE2013-65472 SENSITIVITY ANALYSIS OF THE FACTORS
More informationComparative Analysis of Computational Methods for Limit-Cycle Oscillations
47th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, May 1-4, 26 Newport, Rhode Island Comparative Analysis of Computational Methods for Limit-Cycle Oscillations Arathi
More informationThe Importance of drag
Drag Computation The Importance of drag Its effects on aircraft performances On the Concorde, one count drag increase (ΔC D =.0001) requires two passengers, out of the 90 ~ 100 passenger capacity, be taken
More informationREYNOLDS NUMBER EFFECTS ON LEADING EDGE RADIUS VARIATIONS OF A SUPERSONIC TRANSPORT AT TRANSONIC CONDITIONS
REYNOLDS NUMBER EFFECTS ON LEADING EDGE RADIUS VARIATIONS OF A SUPERSONIC TRANSPORT AT TRANSONIC CONDITIONS S. M. B. Rivers *, R. A. Wahls, and L. R. Owens Aerodynamics, Aerothermodynamics, and Acoustics
More informationAN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION. Abstract
nd Workshop on CFD Uncertainty Analysis - Lisbon, 19th and 0th October 006 AN UNCERTAINTY ESTIMATION EXAMPLE FOR BACKWARD FACING STEP CFD SIMULATION Alfredo Iranzo 1, Jesús Valle, Ignacio Trejo 3, Jerónimo
More informationAN IMPLICIT-EXPLICIT FLOW SOLVER FOR COMPLEX UNSTEADY FLOWS
AN IMPLICIT-EXPLICIT FLOW SOLVER FOR COMPLEX UNSTEADY FLOWS A DISSERTATION SUBMITTED TO THE DEPARTMENT OF AERONAUTICAL ASTRONAUTICAL ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY
More informationLARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE
20 th Annual CFD Symposium, August 09-10, 2018, Bangalore LARGE EDDY SIMULATION OF FLOW OVER NOZZLE GUIDE VANE OF A TRANSONIC HIGH PRESSURE TURBINE Bharathan R D, Manigandan P, Vishal Tandon, Sharad Kapil,
More informationApplications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin
Applications of Harmonic Balance Method in Periodic Flows Gregor Cvijetić and Inno Gatin Faculty of Mechanical Engineering and Naval Architecture, Zagreb 10000, Croatia, gregor.cvijetic@gmail.com The Harmonic
More informationAdjoint-based aerodynamic shape optimization
IT Licentiate theses 2003-012 Adjoint-based aerodynamic shape optimization OLIVIER AMOIGNON UPPSALA UNIVERSITY Department of Information Technology Adjoint-based aerodynamic shape optimization BY OLIVIER
More informationTime Spectral Method for Periodic and Quasi- Periodic Unsteady Computations on Unstructured Meshes
University of Wyoming Wyoming Scholars Repository Mechanical Engineering Faculty Publications Mechanical Engineering 1-211 Time Spectral Method for Periodic and Quasi- Periodic Unsteady Computations on
More informationSIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 2011 GT2011 June 6-10, 2011, ancouver, Canada GT2011-46367 SIMULATION OF UNSTEADY TURBOMACHINERY FLOWS USING AN IMPLICITLY COUPLED NONLINEAR HARMONIC BALANCE METHOD Jonathan
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationEvaluation of the Drag Reduction Potential and Static Stability Changes of C-130 Aft Body Strakes
U.S. Air Force T&E Days 2009 10-12 February 2009, Albuquerque, New Mexico AIAA 2009-1721 Evaluation of the Drag Reduction Potential and Static Stability Changes of C-130 Aft Body Strakes Heather G. Pinsky
More informationCFD COMPUTATION OF THE GROUND EFFECT ON AIRPLANE WITH HIGH ASPECT RATIO WING
28 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES CFD COMPUTATION OF THE GROUND EFFECT ON AIRPLANE WITH HIGH ASPECT RATIO WING Sun Tae Kim*, Youngtae Kim**, Tae Kyu Reu* *Agency for Defense Development,
More informationRoyal Aeronautical Society 2016 Applied Aerodynamics Conference Tuesday 19 th Thursday 21 st July Science Centre, Bristol, UK
Assessment and validation of aerodynamic performance results for a Natural Laminar Flow transonic wing tested in cryogenic conditions via simulation of turbulent wedges in CFD. Royal Aeronautical Society
More informationStability-Constrained Aerodynamic Shape Optimization of a Flying Wing Configuration
13th AIAA/ISSMO Multidisciplinary Analysis Optimization Conference 13-15 September 2010, Fort Worth, Texas AIAA 2010-9199 13th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference September
More informationResolving the dependence on free-stream values for the k-omega turbulence model
Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok Resolving the dependence on free-stream values for the k-omega turbulence model J.C. Kok This report is based on
More informationA Simple Explanation of the Sobolev Gradient Method
A Simple Explanation o the Sobolev Gradient Method R. J. Renka July 3, 2006 Abstract We have observed that the term Sobolev gradient is used more oten than it is understood. Also, the term is oten used
More informationMixing-Plane Method for Flutter Computation in Multi-stage Turbomachines
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 2009, Orlando, Florida AIAA 2009-862 Mixing-Plane Method for Flutter Computation in Multi-stage
More informationHigh Order Discontinuous Galerkin Methods for Aerodynamics
High Order Discontinuous Galerkin Methods for Aerodynamics Per-Olof Persson Massachusetts Institute of Technology Collaborators: J. Peraire, J. Bonet, N. C. Nguyen, A. Mohnot Symposium on Recent Developments
More informationFlow field in the compressor blade cascade NACA
Flow field in the compressor blade cascade NACA 65-100 Tomáš Turek Thesis Supervisor: Ing. Tomáš Hyhlík, Ph.D. Abstract An investigation is made into the effects of a flow field in the compressor blade
More informationReduced order model of three-dimensional Euler equations using proper orthogonal decomposition basis
Journal of Mechanical Science and Technology 24 (2) (2010) 601~608 www.springerlink.com/content/1738-494x DOI 10.1007/s12206-010-0106-0 Reduced order model of three-dimensional Euler equations using proper
More informationTowards Reduced Order Modeling (ROM) for Gust Simulations
Towards Reduced Order Modeling (ROM) for Gust Simulations S. Görtz, M. Ripepi DLR, Institute of Aerodynamics and Flow Technology, Braunschweig, Germany Deutscher Luft und Raumfahrtkongress 2017 5. 7. September
More informationThe Study on Re Effect Correction for Laminar Wing with High Lift
The Study on Re Effect Correction for Laminar Wing with High Lift Jieke Yao, Wenliang Feng, Lingying Lv and Bin Chen Chengdu Aircraft Industrial (group) CO.LTD, 692, Chengdu, China Abstract. In the past
More informationGerman Aerospace Center (DLR)
German Aerospace Center (DLR) AEROGUST M30 Progress Meeting 23-24 November 2017, Bordeaux Presented by P. Bekemeryer / J. Nitzsche With contributions of C. Kaiser 1, S. Görtz 2, R. Heinrich 2, J. Nitzsche
More informationA hybrid adjoint approach applied to RANS equations
Center for urbulence Research Annual Research Briefs 2013 303 A hybrid adjoint approach applied to RANS equations By. W. R. aylor, F. Palacios, K. Duraisamy AND J. J. Alonso 1. Motivation and objectives
More informationTotal Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation
Total Pressure Losses Minimization in Turbomachinery Cascades, using a New Continuous Adjoint Formulation D.I. Papadimitriou, K.C. Giannakoglou Laboratory of Thermal Turbomachines, National Technical University
More informationMestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 2012/13
Mestrado Integrado em Engenharia Mecânica Aerodynamics 1 st Semester 212/13 Exam 2ª época, 2 February 213 Name : Time : 8: Number: Duration : 3 hours 1 st Part : No textbooks/notes allowed 2 nd Part :
More informationSeveral forms of the equations of motion
Chapter 6 Several forms of the equations of motion 6.1 The Navier-Stokes equations Under the assumption of a Newtonian stress-rate-of-strain constitutive equation and a linear, thermally conductive medium,
More informationComputational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2
Computational Fluid Dynamics Study Of Fluid Flow And Aerodynamic Forces On An Airfoil S.Kandwal 1, Dr. S. Singh 2 1 M. Tech Scholar, 2 Associate Professor Department of Mechanical Engineering, Bipin Tripathi
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More information