Stabilization of Explicit Flow Solvers Using a Proper Orthogonal Decomposition Technique
|
|
- Job Perry
- 5 years ago
- Views:
Transcription
1 50th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition January 2012, Nashville, Tennessee AIAA Stabilization of Explicit Flow Solvers Using a Proper Orthogonal Decomposition Technique Kivanc Ekici 1,, Kenneth C. Hall 2, and Huang Huang 1, 1 University of Tennessee, Knoxville, Tennessee Duke University, Durham, North Carolina A new numerical technique for the stabilization of explicit computational fluid dynamic (CFD) solvers is presented. When using CFD codes to solve practical problems, one sometimes finds the solution fails to converge due to the presence of a small number of eigenvalues of the flow solver outside the unit circle. In this paper, we use the proper orthogonal decomposition technique to estimate the unstable eigenvalues and eigenmodes of the flow solver as the solution diverges for linear problems, or exhibits limit cycle oscillations for nonlinear problems. We use the resulting eigen-information to construct a preconditioner that repositions the unstable eigenvalues from outside the unit circle to a point inside the unit close to the origin, resulting in a stable flow solver. In this work, the proposed method is applied to a nonlinear steady cascade solver. However, the technique can be easily applied to external flow solvers, and to unsteady frequency-domain (time-linearized and harmonic balance) flow solvers. I. Nomenclature A i, B i i th -stage coefficients of Runge-Kutta integration D artificial dissipation vector E total energy F, G flux vectors h total enthalpy p pressure P preconditioner P r l laminar Prandtl number P r t turbulent Prandtl number Q discretized convective and viscous fluxes R residual vector S source vector S t source term for turbulence model t time U vector of conservation variables u, v Cartesian velocities V control volume x, y Cartesian coordinates X matrix of right eigenvectors γ specific heat ratio µ l dynamic laminar viscosity µ t dynamic turbulent viscosity ν kinematic viscosity ν working variable for Spalart-Allmaras turbulence model ρ density Assistant Professor, Dept. of Mechanical, Aerospace and Biomedical Engineering. Senior Member AIAA. Julian Francis Abele Professor, Dept. of Mechanical Engineering and Materials Science. Fellow AIAA. Graduate Research Assistant, Dept. of Mechanical, Aerospace and Biomedical Engineering. Student Member AIAA. Copyright c 2012 by Kivanc Ekici, Kenneth C. Hall and Huang Huang. Published by the American Institute of Aeronautics and Astronautics, Inc. with permission. Copyright 2012 by Kivanc Ekici, Kenneth C. Hall and Huang Huang. Published by the, Inc., with permission. 1 of 15
2 II. Introduction Frequency-domain computational fluid dynamic (CFD) models of aerodynamic and aeroelastic phenomena have improved dramatically over the past several decades. In the turbomachinery community, investigators have used two main families of theories based on frequency-domain modeling: time-linearized methods 1 6 and nonlinear harmonic balance methods 7 14 to simulate unsteady flows. The time-linearized approach assumes that the unsteady disturbances are small compared to the mean flow, and are harmonic in time with frequency ω. These assumptions allow one to linearize the governing equations about a nonlinear steady or mean operating condition. The resulting linearized equations can be solved very efficiently using convergence acceleration techniques such as multigrid and local time stepping. In the nonlinear harmonic balance approach, the unsteady flow is also assumed to be temporally periodic with a fundamental excitation frequency ω. However, unlike the time-linearized approach, the unsteady disturbances are not assumed to be small. Using the harmonic balance technique, the dependent variables of the governing Navier-Stokes equations are stored at a number of sub-time levels spanning a single period of the unsteadiness. The various physical time level solutions are only coupled to one another through a pseudo-spectral operator that approximates the time derivatives in the Navier-Stokes equations. As in the time-linearized flow solvers, pseudo time-marching is used to drive the solution to convergence, allowing one to model dynamically nonlinear flows very efficiently. Often, nominally steady or frequency-domain unsteady flow solutions fail to converge. In the case of the timelinearized flow solvers, the solution may unexpectedly diverge exponentially. For steady or unsteady harmonic balance solvers, the solution may diverge, or more commonly, exhibit small limit cycle oscillations. The root cause of the lack of convergence may be physical, as in the case of small scale shedding behind a blunt trailing edge, or a spurious numerical artifact. In any event, the instability in the flow solution is often limited to a very small number of unstable modes. In most of the work reported in the literature, frequency-domain and steady flow solvers have not exhibited convergence problems due to numerical instabilities associated with explicit pseudo-time marching. A few exceptions are the work of Ning and He, 15 Campobasso and Giles, 16, 17 and Elliot. 18 As an example, Campobasso and Giles 16, 17 investigated a configuration in which a steady-state flow solver underwent a low-level limit cycle oscillation due to the existence of a suction side separation bubble in a turbine blade row. This limit cycle oscillation did not adversely affect the accuracy of the mean flow solution because of the small size of the oscillation. However, for the accompanying time-linearzed flow solver, the problem appeared as an exponentially diverging solution. In their work, Campobasso and Giles developed two techniques based on the GMRES algorithm and the Recursive Projection Method (RPM) to stabilize their time-linearized flow solver. In this paper, we propose a new stabilization technique based on proper orthogonal decomposition (POD). POD has been used by a number of investigators to analyze turbulence, 19 to construct reduced models of fluid dynamic systems, 20, 21 and for other purposes. Here we use the POD technique to analyze flow solver instabilities, and then use this information to precondition the flow solver to eliminate any instabilities. The current approach is similar in some ways to the RPM method developed by Campobasso and Giles in that we seek to identify the unstable eigenvalues and eigenmodes of the flow solver, and then modify the the eigenvalues to produce a stable solver. Our method is designed for use with steady nonlinear solvers, but may also be used with time-linearized or harmonic balance solvers. III. Flow Solver A. Governing Equations To demonstrate the POD stabilization method, we apply the method to the two-dimensional Reynolds-Averaged Navier-Stokes (RANS) equations, with the Spalart-Allmaras 22 turbulence model to model Reynolds stresses. In strong conservation form, these equations are given by U t + F x + G y = S (1) 2 of 15
3 where U is the vector of conservation variables, i.e., ρ ρu U = ρv ρe ρ ν The last entry ν, is the working variable in the Spalart-Allmaras turbulence model, from which the eddy viscosity is computed. The flux vectors F, G and the source term S are given by ρu ρv 0 ρuv τ yx F = ρu 2 + p τ xx ρuv τ xy ρuh τ xh ρu ν τ xν, G = ρv 2 + p τ yy, S = ρvh τ yh ρv ν τ yν The shear stresses, τ xx and τ xy, in the flux vectors are defined as τ xx = 2 ( 3 (µ l + µ t ) 2 u x v ) y ( u τ xy = (µ l + µ t ) y + v ) x In addition, the τ xh term in the energy equation and the τ xν term in the turbulence model equation are given by ( µl τ xh = uτ xx + vτ xy + + µ ) t h P r l P r t x τ xν = 3 2 (µ l + ρ ν) ν x The remaining components of the shear stresses and terms in the energy and turbulence model equations are defined similarly. For an ideal gas with a constant specific heat ratio, the pressure is related to the conserved variables through [ p = (γ 1)ρ E 1 ( u 2 + v 2)] 2 Finally, the total enthalpy is defined as h = ρe + p ρ It should be noted that the inviscid flux and source vectors depend on the conservation variables and the Cartesian coordinates. The viscous fluxes depend on the gradients of the flow velocities, temperature and the working turbulent variable. B. Numerical Integration Having presented the governing equations to be used in our analysis, we next integrate Eq. (1) on a computational grid of H-O-H topology using an explicit time-marching scheme introduced by Jameson. 23 The fact that the steady flow is periodic in the circumferential direction allows us to use a computational grid spanning a single blade, thereby reducing the computational cost. In addition, to simplify the treatment of boundary conditions, we use a vertex-based approach 24 in our discretization. Because the steady-state solution is of interest in this work, convergence acceleration techniques such as local time stepping are incorporated in the flow solver. The details of Jameson s scheme is well-documented in the literature However, its salient features are presented here for completeness. S t 3 of 15
4 In the vertex-based Jameson scheme, 24 a control volume is formed by the union of the cells that meet at the vertex. In the two-dimensional case the control volume is made of four cells surrounding the vertex. Written in a semi-discrete form, the Navier-Stokes are given by d (VU) + Q (U) D (U) = 0 dt where V is the volume of union cells, Q (U) denotes to the discretized convective and viscous flux terms, and D (U) denotes to the artificial dissipative term. The convective flux at the cell surface center is evaluated by averaging the conservation variables at the corresponding cell vertex. The flux balance at each vertex can be ensured by summing all the fluxes from each constituting cell. To eliminate odd-even decoupling and oscillations around shock waves, a blend of second and fourth-order dissipative terms, D (U), are included in Jameson s formulation. The system of differential equations are marched temporally toward a steady state solution with a modified multi-stage Runge-Kutta scheme. The conserved variables at the (i + 1)-st stage of each iteration are computed from U (i+1) = U (0) A (i+1) t V [ ( Q U (i)) B (i+1) D (U (i)) ( ) 1 B (i+1) D (U (0))] where U (0) are the conserved variables from the previous iteration. For the five-stage Runge-Kutta integration used in this work, coefficients are given by A 1 = 1 4, A 2 = 1 6, A 3 = 3 8, A 4 = 1 2, A 5 = 1 B 1 = 1, B 2 = 0, B 3 = 14 25, B 4 = 0, B 5 = To reduce the computational cost and maintain a good high frequency error damping property, the viscous fluxes are frozen at the first stage and the artificial dissipative terms are evaluated only at the first, third and fifth stages. IV. Stabilization In general, an explicit CFD scheme can be thought as an iterative process that drives a vector of residuals R to convergence, that is, U n+1 = U n + R (U n ) (2) where n is the iteration number and U is the vector of conserved variables. The residual operator, R, is nonlinear for a steady flow solver. However, to motivate the development of our method (and to simplify the discussion), assume for the moment that R can be written in the linear form given by, R (U n ) = AU n b (3) where the matrix A and the vector b are functions of the mean flow. Note that, in this work we use an explicit Jameson scheme, but the method described in this paper can be applied to any scheme having the generic form given by Eq. (2). Inserting Eq. (3) into Eq. (2), one obtains, U n+1 = (I + A) U n b (4) This scheme is unstable whenever the magnitude of any eigenvalue z of (I + A) is greater than or equal to unity. Therefore, if one can identify those unstable eigenmodes and modify the eigenvalues of the scheme so that all the eigenvalues lie inside the unit circle, then the numerical scheme will be stable ensuring convergence of the explicit scheme. Consider for the moment the homogeneous problem (b = 0). In this case, the residual operator given in Eq. (4) can be rewritten as, U n+1 = U n + AU n (5) Further suppose the ith eigenvalue z i is the only eigenvalue that lies outside the unit circle. To stabilize the homogeneous solver in this simple case, we propose a modification to the explicit operator of the form U n+1 = U n + ( I + β i x i y T ) i AU n (6) 4 of 15
5 where x i and y i are the right and left eigenvectors associated with the unstable eigenvalue and β i is a constant to be determined. Note that if the numerical solution of the steady flow field is stable, β i is taken to be zero and the original form of the operator is recovered. However, if the solution is unstable (including the scenario where the residual value goes into a limit cycle) the value of β i will be non-zero. In general, since the linear homogeneous problem is considered, the solution can be thought of as a linear combination of the eigenvectors of the system, so that U = Xα (7) where X is the matrix of right eigenvalues and α is a vector that represents the weighting coefficients of each eigenvector. Substituting Eq. (7) into Eq. (6), one gets Xα n+1 = Xα n + ( I + β i x i y T ) i AXα n (8) Next, pre-multiplying Eq. (8) by the transpose of the left eigenvector (y T i ) and making use of the orthogonality of the right and the left eigenvectors, i.e., { y T 0 i j i x j = 1 i = j and one gets y T i (I + A) x j = { 0 i j i = j z i α n+1 i = α n i + [z i + β i (z i 1)] α n i (9) where α i is the generalized coordinate associated with the unstable eigenvalue z i. The effect of the preconditioning, then, is to change the eigenvalue of the unstable mode from z i to z i + β i (z i 1), leaving the stable eigenvalues and eigenmodes unchanged. To completely eliminate the unstable error mode in one iteration, we set z i β i = z i 1 (10) so that the modified eigenvalue is zero. Recall that to simplify the discussion for the development of the technique, we have considered that only one eigenvalue was unstable and that the problem was homogeneous. For the more general case where more than one eigenvalue is unstable, and b 0, the preconditioned scheme can be written as [ ] U n+1 = U n + I + i β i x i y T i [AU n b] (11) where the sum is over the set of unstable eigenmodes. Similarly, for nonlinear CFD schemes, the preconditioning technique can be expressed as [ ] or U n+1 = U n + I + i β i x i y T i R (U n ) (12) U n+1 = U n + PR (U n ) (13) where the eigenvalues and the eigenvectors are those of the nonlinear problem linearized about some mean flow and P is the preconditioning matrix. Therefore, to stabilize a diverging or limit cycling CFD solver, we use eigen-information associated with the unstable or limit cycling modes to precondition the solver. In the next section, we present a proper orthogonal decomposition (POD) technique used to compute approximately the unstable eigenvalues and eigenmodes. 5 of 15
6 V. Determination of POD Vectors For the case where the steady solution residual is diverging or exhibiting limit cycle, the unsteadiness in the solution will be dominated by just a few modes of motion. One way to identify these modes is to use proper orthogonal decomposition (POD). The process is as follows. Every few iterations, the entire flow solution is stored. Each of these several stored solutions or snapshots is loaded as a column vector into a tall, narrow matrix (S), so that S = { U 1 Ū} { U 2 Ū} {... U K Ū} (14) where U 1, U 2, are the solution snapshots and Ū = 1 K K U k (15) is the average of the solution snapshots. We seek to construct optimal basis vectors that are linear combinations of the columns of the snapshot matrix, that is Φ = Sv (16) The so-called proper orthogonal decomposition vectors are found by first solving the small eigenvalue problem defined by k=1 S T Sv k = λ k v k, k = 1,..., K (17) where λ k and v k are the eigenvalues and eigenvectors of the matrix S T S. Those eigenmodes with the largest eigenvalues are dominant. Next, we use several of the dominant POD vectors to estimate the eigenvalues and eigenmodes of the flow solver, that is, we assume that, x i Φγ i (18) where γ i is a vector whose length is equal to the number of dominate POD vectors retained in Φ. Thus the eigenvalue problem of the flow solver, (I + A) x i = z i x i (19) can be rewritten as (I + A) Φγ i = z i Φγ i (20) Finally, because this system is over constrained, we project the error against the basis vectors to obtain the reduced order eigenvalue problem, Φ T (I + A) Φγ i = z i Φ T Φγ i ; i = 1,, I (21) The above equation approximates the eigenvalues of the explicit solution scheme using a reduced order model, which is small (on the order of 12 12) and very inexpensive to calculate. Note that, although A may not be easily determined, the product AΦ can be easily obtained using numerical differentiation eliminating the need to explicitly linearize the code. In this work, we use second-order finite differencing to calculate AΦ in a matrix-free approach. VI. Numerical Results To demonstrate the effectiveness of the novel stabilization technique proposed in this work, we consider the Eleventh Standard Configuration, 27 which is a two-dimensional viscous turbine cascade (see Fig. 1). The turbine row is comprised of airfoils with stagger angle of degrees and a chord-to-pitch ratio (solidity) of Experimental measurements for this configuration were reported by Fransson et al. 27 In addition, a number of investigators presented unsteady computations 16, 17, for this test case to validate their flow solvers. Experimental data is available for a number of different flow conditions for this turbine. However, we consider only the transonic off-design condition in this paper because this (nominally steady) case demonstrates a numerical 6 of 15
7 limit cycle oscillation caused by a separation bubble that forms on the suction side of the blades near the leading edge. Note that, the same test case was also used by Campobasso and Giles 16, 17 to validate their numerical stabilization techniques based on the GMRES and Recursive Projection Methods (RPM) applied to unsteady time-linearized flow solvers. Before presenting the results of the POD stabilization scheme, we present results of the unmodified CFD scheme to illustrate this limit cycle. A. Unmodified Flow Solver Solutions (No POD Stabilization) Generally speaking, because of the use of local time-stepping, the solution of nonlinear steady flows is not timeaccurate. However, the pseudo time-marching of the governing equations reflects some physical properties of the flow field causing a limit cycle behavior in convergence. 16 Figure 1 shows the H-O-H viscous computational grid we used to compute the flow field; the grid contains nodes in the O-grid block, and nodes in each of the two H-grid blocks. The flow enters the blade Figure 1. H-O-H grid used for flow computations. row with a Mach number of 0.38 and a flow angle of 34 deg with respect to the axial direction. The Reynolds number, based on inlet flow conditions and the chord length, is approximately 1,100,000. The isentropic Mach number at the exit of the blade row is 0.99 causing a weak shock formation on the suction side of the blades at approximately 80% of the chord. As mentioned earlier, a separation bubble forms between the leading edge and 30% chord of the suction side [see Fig. 2] that introduces a numerical instability associated with the bubble, and hence, the steady flow solver fails to converge. As can be seen from Fig. 3, the steady solution residual drops one order of magnitude within the first 750 iterations, but then grows, and eventually goes into a periodic limit cycle around 3000 iterations. A fully converged solution cannot be achieved for this case without the POD stabilization. Nevertheless, the (unconverged) contour plots from our nominal flow solver, shown in Fig. 4 at the end of 5000 iterations, 16, 17 agree very well with the computations of Campobasso and Giles. The pressure and absolute Mach number contours within the passage clearly show the shock formation as well as the separation bubble on the suction side. As reported by Campobasso and Giles, 16 a small-amplitude limit-cycle behavior about a mean flow does not necessarily prevent the steady solution from converging to a level acceptable for engineering purposes. Next, we compare the results obtained from our nominal solver (no POD-stabilization applied) and available experimental data. Shown in Fig. 5 are the isentropic Mach number values on the surface of the blade. Overall, the predicted values are in good agreement with the experimental data except for a small region upstream of the 7 of 15
8 Figure 2. Flow field near the leading edge. Left: Streamlines. Right: Velocity vectors Log(Global Residual) Iteration Number Figure 3. Convergence history. Note the limit-cycle behavior instability due to shedding. 8 of 15
9 Figure 4. Contour plots for the Eleventh Standard Configuration. Left: Mach number. Right: Static pressure. Note that a single blade passage is used in computations. Multiple passages are shown here for clarity. weak shock. However, as pointed out by Fransson,27 fairly small changes in the inlet flow conditions may result in significantly different flow structures near the shock, which was also observed by other investigators. Figure 5 also reveals small oscillations in surface isentropic Mach number values near the leading edge of the blades. These oscillations are mainly due to the instability caused by a shedding phenomenon near the separation bubble. 2 Surface isentropic Mach number Present solver Experiment Surface location, x/c Figure 5. Surface isentropic Mach number distribution. Note small oscillations caused by the numerical instability near the leading edge. 9 of 15
10 B. Modified Flow Solver Solutions (With POD Stabilization) In this section, we demonstrate how our proposed technique can be used to identify and eliminate the nonlinear limit cycle oscillation found in the example of the previous section. As shown in the previous section, the steady flow computation for the Eleventh Standard test case fails to converge due to a numerical instability caused by a physical phenomenon, that is, the shedding near the separation bubble is clearly seen in the contours of ρe shown in Fig. 6. Here contour plots are generated at every three iterations from the 5001 st to the 5031 st iteration and it is apparent that the limit-cycle behavior has a period spanning approximately 30 iterations. Evidence of shedding can also be seen from the convergence history plot presented in Fig. 7. Figure 6. Contours of ρe for different iteration numbers. Top left to right: Iteration numbers from 5001 to Bottom left to right: Iteration numbers from 5016 to The first step in the POD stabilization algorithm is to assemble the solution snapshots of the POD basis vectors. After the solution goes into limit-cycle at the end of 5000 iterations, we store the computed CFD computations every three iterations for the next 30 iterations. This results in a total of eleven snapshots. We choose the snapshots to span a single cycle of the limit cycle. Shorter periods will not contain all the relevant solution behavior to obtain good POD vectors; longer periods provide no additional useful information. Although not presented here, we also tried fewer or more snapshots and found that ten to 15 snapshots are needed to stabilize the solver considered in this work. Having collected the solution snapshots, the next step is to compute the POD modes using Eq. (17). Next, these POD modes are used to compute estimates of the unstable eigenvalues and eigenvectors using Eq (21). Figure 8 shows the eleven eigenvalues computed in this way. As can be seen, the eigenanalysis reveals two unstable modes, which lie just outside the unit circle and are complex conjugates of one another. Note that we perform the eigenanalyses to determine the unstable modes just once after the stabilization is turned on, and in any event is a very inexpensive operation because of the small size of the problem. Having identified the offending eigenmodes, at each subsequent iteration of the flow solver, we precondition the flow solver using Eq. (12) where the values of β i are computed from Eq. (10) for the two unstable modes. Figure 9 shows the convergence history for the stabilized solver. It can be seen that once the POD stabilization is turned on, the residual drops three orders of magnitude in a few hundred iterations eliminating the numerical instability and fully converges to machine accuracy in around 7500 iterations. Note that for the preconditioning here, we used a total of eleven solution snapshots to compute the five most dominant POD vectors, which in 10 of 15
11 -4.70 Log(Global Residual) Iteration Number Figure 7. Convergence history. Note the limit-cycle behavior instability due to shedding Unstable modes Im(λ) Re(λ) Figure 8. Eigenvalues of the CFD solution. Note that eleven snapshots were used and no order reduction was performed. 11 of 15
12 turn were used to find the two least stable flow solver eigenmodes used to construct the preconditioning. When the total number of POD vectors used was four or fewer, the computations diverged, presumably because the computed flow solve eigen-information was not sufficiently accurate to construct an effective preconditioner. 0.0 POD Stabilization enabled after 5000 iterations -5.0 Log(Global Residual) Iteration Number Figure 9. Convergence history using the POD stabilization technique. Next, we plot the surface isentropic Mach number values when the POD stabilization is applied. As seen in Fig. 10, the stabilization technique eliminates the spatial oscillations near the leading edge of the blades. To further investigate the effectiveness of the proposed method, we turned off the preconditioning at the end of iterations, restarted the solver from the fully-converged solution, and performed computations for an additional iterations. Figure 11 shows the global residual for this case. It is apparent that the numerical instability associated with the shedding phenomenon causes the solution to diverge slowly when the POD technique is disabled. Finally, to re-stabilize the flow solver, we revert back to the preconditioning technique and restart the solution after iterations. Figure 12 shows the convergence histories for the preconditioning constructed using two different sets of solution snapshots. When the snapshot are taken every ten iterations as the solution is diverging (between and iterations), the convergence rate is the slowest. This is not surprising, because the period of the instability is around 15 iterations. When the solution snapshots are taken at more closely spaced intervals, the convergence rate increases because the POD vectors provide a better basis for estimating the eigenvalues and eigenmodes of the system that are used to construct the preconditioner. VII. Conclusions A numerical stabilization method based on a novel POD technique was presented for explicit CFD solvers. Using this scheme, we first take snapshots of the non-converging solution. Using these snapshots, we form POD basis vectors, which we then use to compute approximate eigenvalues and eigenmodes of the unstable CFD modes. Finally, the estimates of the unstable eigenvalues and eigenmodes are used to precondition the explicit scheme, resulting in a stable flow solver. The method was applied to a representative flow solver, and was found to produce fully converged solutions for cases that would otherwise not converge. 12 of 15
13 2 Surface isentropic Mach number Present solver (No POD stabilization) Present solver (With POD stabilization) Experiment Oscillations eliminated with POD stabilization Surface location, x/c Figure 10. Surface isentropic Mach number distribution. Note how POD technique eliminates oscillations caused by the numerical instability near the leading edge. 0.0 POD Stabilization enabled after 5000 iterations Log(Global Residual) POD Stabilization disabled after iters Log(Global Residual) Iteration Number Iteration Number Figure 11. Convergence histories. Left: Solution restarted at the end of iterations and without the POD stabilization technique. Right: Close-up view. 13 of 15
14 0.0 New snapshots (Every 10 iters from ) New snapshots (Every iteration from ) Log(Global Residual) POD stabilization restarted Iteration Number Figure 12. Convergence histories for restarted POD solver References 1 Hall, K. C. and Crawley, E. F., Calculation of Unsteady Flows in Turbomachinery Using the Linearized Euler Equations, AIAA Journal, Vol. 27, No. 6, June 1989, pp Hall, K. C. and Lorence, C. B., Calculation of Three-Dimensional Unsteady Flows in Turbomachinery Using the Linearized Harmonic Euler Equations, Journal of Turbomachinery, Vol. 15, No. 4, Oct. 1993, pp Hall, K. C., Clark, W. S., and Lorence, C. B., A Linearized Euler Analysis of Unsteady Transonic Flows in Turbomachinery, Journal of Turbomachinery, Vol. 116, July 1994, pp Clark, W. S. and Hall, K. C., A Time-Linearized Navier-Stokes Analysis of Stall Flutter, Journal of Turbomachinery, Vol. 122, July 2000, pp Sbardella, L. and Imregun, M., Linearized unsteady viscous turbomachinery flows using hybrid grids, Journal of Turbomachinery, Vol. 123, 2001, pp Hall, K. C. and Ekici, K., Multistage Coupling for Unsteady Flows in Turbomachinery, AIAA Journal, Vol. 43, No. 3, March 2005, pp He, L. and Ning, W., Efficient Approach for Analysis of Unsteady Viscous Flows in Turbomachines, AIAA Journal, Vol. 36, No. 11, 1998, pp Hall, K. C., Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique, presented at the Kerrebrock Symposium, A Symposium in Honor of Professor Jack L. Kerrebrock s 70th Birthday, Massachusetts Inst. of Technology, Cambridge, MA, Jan Hall, K. C., Thomas, J. P., and Clark, W. S., Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique, AIAA Journal, Vol. 40, No. 5, May 2002, pp Ekici, K. and Hall, K. C., Nonlinear Analysis of Unsteady Flows in Multistage Turbomachines Using Harmonic Balance, AIAA Journal, Vol. 45, No. 5, May 2007, pp Ekici, K., Hall, K. C., and Dowell, E. H., Computationally Fast Harmonic Balance Methods for Unsteady Aerodynamic Predictions of Helicopter Rotors, Journal of Computational Physics, Vol. 227, No. 12, 2008, pp Ekici, K. and Hall, K. C., Nonlinear Frequency-Domain Analysis of Unsteady Flows in Turbomachinery with Multiple Excitation Frequencies, AIAA Journal, Vol. 46, No. 8, Aug. 2008, pp Ekici, K. and Hall, K. C., Harmonic Balance Analysis of Limit Cycle Oscillations in Turbomachinery, AIAA Journal, Vol. 49, No. 7, July Huang, H. and Ekici, K., A Harmonic Balance Method for the Analysis of Unsteady Flows in Cascades, AIAA Paper , Ning, W. and He, L., Some Modeling Issues on Trailing-Edge Vortex Shedding, AIAA Journal, Vol. 39, No. 5, 2001, pp Campobasso, M. S. and Giles, M. B., Effects of Flow Instabilities on the Linear Analysis of Turbomachinery Aeroelasticity, Journal of Propulsion and Power, Vol. 19, No. 2, 2003, pp of 15
15 17 Campobasso, M. S. and Giles, M. B., Stabilization of Linear Flow Solver for Turbomachinery Aeroelasticity Using Recursive Projection Method, AIAA Journal, Vol. 42, No. 9, 2004, pp Elliot, J. K., Aerodynamic optimization based on the Euler and Navier-Stokes Equations using Unstructured grids, Ph.D. thesis, Massachusetts Institute of Technology, Cambridge, MA, may Berkooz, G., Holmes, P., and Lumley, J. L., The Propoer Orthogonal Decomposition in the Analysis of Turbulent Flows, Annual Review of Fluid Mechanics, Vol. 25, 1993, pp Willcox, K. and Peraire, J., Balanced model reduction via the proper orthogonal decomposition, AIAA Journal, Vol. 40, 2002, pp Thomas, J. P., Dowell, E. H., and Hall, K. C., Three-dimensional Transonic Aeroelasticity Using Proper Orthogonal Decomposition Based Reduced Order Models, Journal of Aircraft, Vol. 40, No. 3, 2003, pp Spalart, P. R. and Allmaras, S. R., One-Equation Turbulence Model for Aerodynamic Flows, AIAA Paper , Jameson, A., Schmidt, W., and Turkel, E., Numerical Solutions of the Euler Equations by Finite Volume Methods Using Runge-Kutta Time-Stepping Schemes, AIAA Paper , Jameson, A., A Vertex Based Multigrid Algorithm for Three Dimensional Compressible Flow Calculations, ASME Symposium on Numerical Methods for Compressible Flows, Martinelli, L. and Jameson, A., Validation of a Multigrid Method for the Reynolds Averaged Equations, AIAA Paper , Radespiel, R., Rossow, C., and Swanson, R. C., Efficient Cell-Vertex Multigrid Scheme for the Three-Dimensional Navier- Stokes Equations, AIAA Journal, Vol. 28, No. 8, 1990, pp Fransson, T. H., M.Jocker, Boelcs, A., and Ott, P., Viscous and Inviscid Linear/Nonlinear Calculation Versus Quasi 3D Experimental Data for a new Aeroelastic Turbine Standard Configuration, Journal of Turbomachinery, Vol 122, No.1, Sbardella, L. and Imregun, M., Linearized Unsteady Viscous Turbomachinery Flow Using Hybrid Grids, Journal of Turbomachinery, Vol. 123, 2001, pp Cinnella, P., Palma, P. D., Pascazio, G., and Napolitano, M., A Numerical Method for Turbomachinery Aeroelasticity, Journal of Turbomachinery, Vol. 126, 2004, pp of 15
Non-Synchronous Vibrations of Turbomachinery Airfoils
Non-Synchronous Vibrations of Turbomachinery Airfoils 600 500 NSV Frequency,!, hz 400 300 200 F.R. Flutter 100 SFV 0 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 Rotor Speed,!, RPM Kenneth C. Hall,
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 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 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 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 informationGT UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE USING AN EFFICIENT NON-LINEAR HARMONIC BALANCE METHOD
Proceedings of ASME Turbo Expo 213: Turbine Technical Conference and Exposition GT213 June 3-7, 213, San Antonio, Texas, USA GT213-94574 UNSTEADY SIMULATION OF A TWO-STAGE COOLED HIGH PRESSURE TURBINE
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 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 informationLimit Cycle Oscillations of a Typical Airfoil in Transonic Flow
Limit Cycle Oscillations of a Typical Airfoil in Transonic Flow Denis B. Kholodar, United States Air Force Academy, Colorado Springs, CO 88 Earl H. Dowell, Jeffrey P. Thomas, and Kenneth C. Hall Duke University,
More informationUsing Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver
Using Automatic Differentiation to Create a Nonlinear Reduced Order Model Aeroelastic Solver Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708 0300 A novel nonlinear
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 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 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 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 informationComputation of NACA0012 Airfoil Transonic Buffet Phenomenon with Unsteady Navier-Stokes Equations
5th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition 9-2 January 22, Nashville, Tennessee AIAA 22-699 Computation of NACA2 Airfoil Transonic Buffet Phenomenon with
More informationParallel Computation of Forced Vibration for A Compressor Cascade
44th AIAA Aerospace Sciences Meeting and Exhibit 9-2 January 26, Reno, Nevada AIAA 26-628 AIAA Paper 26-628 Parallel Computation of Forced Vibration for A Compressor Cascade Zongjun Hu and Gecheng Zha
More informationNumerical Researches on Aeroelastic Problem of a Rotor due to IGV/Fan Interaction
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 9, Orlando, Florida AIAA 9-865 Numerical Researches on Aeroelastic Problem of a Rotor due to IGV/Fan
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 informationAPPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS
TASK QUARTERLY Vol. 17, Nos 3 4, 2013, pp. 145 154 APPLICATION OF HYBRID CFD/CAA TECHNIQUE FOR MODELING PRESSURE FLUCTUATIONS IN TRANSONIC FLOWS SŁAWOMIR DYKAS, WŁODZIMIERZ WRÓBLEWSKI AND SEBASTIAN RULIK
More informationElliptic Trailing Edge for a High Subsonic Turbine Cascade
Elliptic Trailing Edge for a High Subsonic Turbine Cascade Mahmoud M. El-Gendi 1, Mohammed K. Ibrahim 2, Koichi Mori 3, and Yoshiaki Nakamura 4 1 Graduate School of Engineering, Nagoya University, Nagoya
More informationAPPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM
AIAA 2003-3693 APPLICATION OF SPACE-TIME MAPPING ANALYSIS METHOD TO UNSTEADY NONLINEAR GUST-AIRFOIL INTERACTION PROBLEM Vladimir V. Golubev* and Axel Rohde Embry-Riddle Aeronautical University Daytona
More informationAn Investigation of the Sensitivity of F-16 Fighter Flutter Onset and Limit Cycle Oscillations to Uncertainties
An Investigation of the Sensitivity of F-16 Fighter Flutter Onset and Limit Cycle Oscillations to Uncertainties Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708
More informationA NONLINEAR HARMONIC BALANCE METHOD FOR THE CFD CODE OVERFLOW 2
IFASD-2009-050 A NONLINEAR HARMONIC BALANCE METHOD FOR THE CFD CODE OVERFLOW 2 Chad H. Custer 1, Jeffrey P. Thomas 2, Earl H. Dowell 3, and Kenneth C. Hall 4 Duke University Mechanical Engineering & Materials
More informationTURBINE FORCED RESPONSE PREDICTION USING AN INTEGRATED NONLINEAR ANALYSIS A.I. Sayma, M. Vahdati & M. Imregun Abstract The forced response due to flow
TURBINE FORCED RESPONSE PREDICTION USING AN INTEGRATED NONLINEAR ANALYSIS A.I. Sayma, M. Vahdati & M. Imregun 1 Imperial College of Science Technology and Medicine Mechanical Engineering Department Exhibition
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 informationUsing Automatic Differentiation to Create a Nonlinear Reduced Order Model of a Computational Fluid Dynamic Solver
Using Automatic Differentiation to Create a Nonlinear Reduced Order Model of a Computational Fluid Dynamic Solver Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 2778
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 informationEigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations
NASA/CR-1998-208741 ICASE Report No. 98-51 Eigenmode Analysis of Boundary Conditions for the One-dimensional Preconditioned Euler Equations David L. Darmofal Massachusetts Institute of Technology, Cambridge,
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 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 informationActive Control of Separated Cascade Flow
Chapter 5 Active Control of Separated Cascade Flow In this chapter, the possibility of active control using a synthetic jet applied to an unconventional axial stator-rotor arrangement is investigated.
More informationUnsteady Flow Computation Using a Harmonic Balance Approach Implemented about the OVERFLOW 2 Flow Solver
9th AIAA Computational Fluid Dynamics Conference // June 22-25, 2009, San Antonio, Texas Unsteady Flow Computation Using a Harmonic Balance Approach Implemented about the OVERFLOW 2 Flow Solver Jeffrey
More informationRELIABILITY OF TIME-LINEARIZED FLUTTER PREDICTIONS NEAR THE SURGE LINE
RELIABILITY OF TIME-LINEARIZED FLUTTER PREDICTIONS NEAR THE SURGE LINE Markus May - Dr. Björn Grüber DLR German Aerospace Center, Institute of Aeroelasticity Bunsenstr. 10, 37073 Göttingen, Germany Markus.May@dlr.de
More informationDERIVATED TURBULENCE MODEL TO PREDICT HARMONIC LOADS IN TRANSONIC SEPARATED FLOWS OVER A BUMP
DERIVATED TURBULENCE MODEL TO PREDICT HARMONIC LOADS IN TRANSONIC SEPARATED FLOWS OVER A BUMP M. Philit 1, P. Ferrand 1, S. Labit, J-C Chassaing, S. Aubert 3, T. Fransson 4 1 UMR CNRS 559, Université de
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 informationNumerical Study of Natural Unsteadiness Using Wall-Distance-Free Turbulence Models
Numerical Study of Natural Unsteadiness Using Wall-Distance-Free urbulence Models Yi-Lung Yang* and Gwo-Lung Wang Department of Mechanical Engineering, Chung Hua University No. 707, Sec 2, Wufu Road, Hsin
More informationNON-SYNCHRONOUS VIBRATION OF TURBOMACHINERY AIRFOILS. AFOSR Contract F Final Report
NON-SYNCHRONOUS VIBRATION OF TURBOMACHINERY AIRFOILS AFOSR Contract F49620-03-1-0204 By Robert E. Kielb, Senior Research Scientist and Principal Investigator Kenneth C. Hall, Professor & Chair Meredith
More informationComputation of cascade flutter by uncoupled and coupled methods
International Journal of Computational Fluid Dynamics, Vol. 19, No. 8, November 25, 559 569 Computation of cascade flutter by uncoupled and coupled methods MANI SADEGHI and FENG LIU* Department of Mechanical
More informationAerodynamic force analysis in high Reynolds number flows by Lamb vector integration
Aerodynamic force analysis in high Reynolds number flows by Lamb vector integration Claudio Marongiu, Renato Tognaccini 2 CIRA, Italian Center for Aerospace Research, Capua (CE), Italy E-mail: c.marongiu@cira.it
More informationSENSITIVITY ANALYSIS FOR FORCED RESPONSE IN TURBOMACHINERY USING AN ADJOINT HARMONIC BALANCE METHOD
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds. Crete Island, Greece, 5 10 June
More informationLINEARISED FREQUENCY DOMAIN GUST ANALYSIS OF LARGE CIVIL AIRCRAFT
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
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 informationAn Introduction to the Discontinuous Galerkin Method
An Introduction to the Discontinuous Galerkin Method Krzysztof J. Fidkowski Aerospace Computational Design Lab Massachusetts Institute of Technology March 16, 2005 Computational Prototyping Group Seminar
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 informationDEVELOPMENT OF AN EFFICIENT DESIGN METHOD FOR NON-SYNCHRONOUS VIBRATIONS
DEVELOPMENT OF AN EFFICIENT DESIGN METHOD FOR NON-SYNCHRONOUS VIBRATIONS by Meredith Anne Spiker Department of Mechanical Engineering and Materials Science Duke University Date: Approved: Robert E. Kielb,
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 informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationBlock-Jacobi Implicit Algorithms for the Time Spectral Method
Block-Jacobi Implicit Algorithms for the Time Spectral Method Frédéric Sicot, Guillaume Puigt and Marc Montagnac CERFACS, CFD/AAM, 42 avenue Coriolis F-3057 Toulouse Cedex 0 Many industrial applications
More informationSome effects of large blade deflections on aeroelastic stability
47th AIAA Aerospace Sciences Meeting Including The New Horizons Forum and Aerospace Exposition 5-8 January 29, Orlando, Florida AIAA 29-839 Some effects of large blade deflections on aeroelastic stability
More informationParallel Computations of Unsteady Three-Dimensional Flows in a High Pressure Turbine
Parallel Computations of Unsteady Three-Dimensional Flows in a High Pressure Turbine Dongil Chang and Stavros Tavoularis Department of Mechanical Engineering, University of Ottawa, Ottawa, ON Canada Stavros.Tavoularis@uottawa.ca
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 informationVirtual Aeroelastic Flight Testing for the F-16 Fighter with Stores
Virtual Aeroelastic Flight Testing for the F-16 Fighter with Stores Jeffrey P. Thomas, Earl H. Dowell, and Kenneth C. Hall Duke University, Durham, NC 27708 0300 and Charles M. Denegri Jr. U.S. Air Force
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 informationExplicit algebraic Reynolds stress models for internal flows
5. Double Circular Arc (DCA) cascade blade flow, problem statement The second test case deals with a DCA compressor cascade, which is considered a severe challenge for the CFD codes, due to the presence
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 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 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 informationPRESENTED AT ICAS 2008 ANCHORAGE, ALASKA SEPTEMBER 2008
PRESENTED AT ICAS 2008 ANCHORAGE, ALASKA SEPTEMBER 2008 NONLINEAR AEROELASTICITY OF CURRENT AND FUTURE AEROSPACE VEHICLES Earl Dowell, Kenneth Hall, Robert Kielb, Deman Tang, and Jeffrey Thomas Duke University
More informationA finite-volume algorithm for all speed flows
A finite-volume algorithm for all speed flows F. Moukalled and M. Darwish American University of Beirut, Faculty of Engineering & Architecture, Mechanical Engineering Department, P.O.Box 11-0236, Beirut,
More 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 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 informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS Hierarchy of Mathematical Models 1 / 29 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 2 / 29
More informationGPPS Acoustic Eigenmode Analysis for Swirling Flows in Ducts
Proceedings of Shanghai 07 Global Power and Propulsion Forum 30 th October st November, 07 http://www.gpps.global GPPS-07-003 Acoustic Eigenmode Analysis for Swirling Flows in Ducts Huang Huang AECC Commercial
More informationImprovements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme
Improvements of Unsteady Simulations for Compressible Navier Stokes Based on a RK/Implicit Smoother Scheme Oren Peles and Eli Turkel Department of Applied Mathematics, Tel-Aviv University In memoriam of
More informationDirect Numerical Simulations of Plunging Airfoils
48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition 4-7 January 010, Orlando, Florida AIAA 010-78 Direct Numerical Simulations of Plunging Airfoils Yves Allaneau
More informationNumerical Analysis of Active Cascade Flutter Control with Smart Structure
Proceedings of the International Gas Turbine Congress Tokyo November -7, IGTCTokyo TS-55 Numerical Analysis of Active Cascade Flutter Control with Smart Structure Junichi Kazawa and Toshinori Watanabe
More informationLocal Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver
Local Preconditioning for Low Mach-number Aerodynamic Flow Simulation with a Steady Compressible Navier-Stokes Flow Solver Kurt Sermeus and Eric Laurendeau Advanced Aerodynamics Department Bombardier Aerospace
More informationLAMINAR-TURBULENT TRANSITION FOR ATTACHED AND SEPARATED FLOW
University of Kentucky UKnowledge University of Kentucky Doctoral Dissertations Graduate School 2010 LAMINAR-TURBULENT TRANSITION FOR ATTACHED AND SEPARATED FLOW Qian Zhang University of Kentucky, zqian@engr.uky.edu
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 informationREDUCED ORDER MODELLING OF GUST ANALYSIS USING COMPUTATIONAL FLUID DYNAMICS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More information2. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations
. Getting Ready for Computational Aerodynamics: Fluid Mechanics Foundations We need to review the governing equations of fluid mechanics before examining the methods of computational aerodynamics in detail.
More informationHybrid LES RANS Method Based on an Explicit Algebraic Reynolds Stress Model
Hybrid RANS Method Based on an Explicit Algebraic Reynolds Stress Model Benoit Jaffrézic, Michael Breuer and Antonio Delgado Institute of Fluid Mechanics, LSTM University of Nürnberg bjaffrez/breuer@lstm.uni-erlangen.de
More informationSimulation of Vortex Shedding in a Turbine Stage
THE AMERICAN SOCIETY OF MECHANICAL ENGINEERS 345 E. 47th St_, New York, N.Y. 10017 ^s-gt-242 SThe Society shall not be responsible for statements or opinions advanced in papers or discussion at meetings
More informationTurbulent Boundary Layers & Turbulence Models. Lecture 09
Turbulent Boundary Layers & Turbulence Models Lecture 09 The turbulent boundary layer In turbulent flow, the boundary layer is defined as the thin region on the surface of a body in which viscous effects
More 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 informationINFLUENCE OF GUST MODELLING ON FREE-FLIGHT AEROFOILS
International Forum on Aeroelasticity and Structural Dynamics IFASD 2017 25-28 June 2017 Como, Italy INFLUENCE OF GUST MODELLING ON FREE-FLIGHT AEROFOILS Reik Thormann 1, Giampaolo Pagliuca 1, Sebastian
More informationGPPS NUMERICAL PREDICTION OF UNSTEADY ENDWALL FLOW AND HEAT TRANSFER WITH ONCOMING WAKE
Proceedings of Shanghai 17 Global Power and Propulsion Forum 3 th October 1 st November, 17 http://www.gpps.global GPPS-17-133 NUMERICAL PREDICTION OF UNSTEADY ENDWALL FLOW AND HEAT TRANSFER WITH ONCOMING
More informationPotential of Aeroelastic Tailoring to Improve Flutter Stability of Turbomachinery Compressor Blades
Paper ID: ETC0-0 Proceedings of th European Conference on Turbomachinery Fluid dynamics & Thermodynamics ETC, April -, 0; Stockholm, Sweden Potential of Aeroelastic Tailoring to Improve Flutter Stability
More informationCHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH
82 CHAPTER 4 OPTIMIZATION OF COEFFICIENT OF LIFT, DRAG AND POWER - AN ITERATIVE APPROACH The coefficient of lift, drag and power for wind turbine rotor is optimized using an iterative approach. The coefficient
More informationIMPLEMENTATION OF POD AND DMD METHODS IN APACHE SPARK FRAMEWORK FOR SIMULATION OF UNSTEADY TURBULENT FLOW IN THE MODEL COMBUSTOR
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationEigenmode analysis. for turbomachinery applications
Eigenmode analysis for turbomachinery applications Pierre Moinier, Michael B. Giles Oxford University Computing Laboratory Oxford, United Kingdom September 16, 25 Abstract This paper discusses the numerical
More informationApplication of turbulence. models to high-lift airfoils. By Georgi Kalitzin
Center for Turbulence Research Annual Research Briefs 997 65 Application of turbulence models to high-lift airfoils By Georgi Kalitzin. Motivation and objectives Accurate prediction of the ow over an airfoil
More informationWALL ROUGHNESS EFFECTS ON SHOCK BOUNDARY LAYER INTERACTION FLOWS
ISSN (Online) : 2319-8753 ISSN (Print) : 2347-6710 International Journal of Innovative Research in Science, Engineering and Technology An ISO 3297: 2007 Certified Organization, Volume 2, Special Issue
More informationA STUDY ON AERODYNAMIC DESIGN OF HIGH- EFFICIENCY PROPELLERS AT HIGH-ALTITUDE
A STUDY ON AERODYNAMIC DESIGN OF HIGH- EFFICIENCY PROPELLERS AT HIGH-ALTITUDE Jianhua Xu, Wenping Song & Wenqing Yang National Key Laboratory of Aerodynamic Design and Research, Northwestern Polytechnical
More informationOpen boundary conditions in numerical simulations of unsteady incompressible flow
Open boundary conditions in numerical simulations of unsteady incompressible flow M. P. Kirkpatrick S. W. Armfield Abstract In numerical simulations of unsteady incompressible flow, mass conservation can
More informationINFLUENCE OF AIR COOLING AND AIR-JET VORTEX GENERATOR ON FLOW STRUCTURE IN TURBINE PASSAGE
TASKQUARTERLYvol.19,No2,2015,pp.153 166 INFLUENCE OF AIR COOLING AND AIR-JET VORTEX GENERATOR ON FLOW STRUCTURE IN TURBINE PASSAGE RYSZARD SZWABA, PAWEŁ FLASZYŃSKI, JAN ARTUR SZUMSKI AND PIOTR DOERFFER
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 informationAeroelastic instabilities simulations on turbofan HP compressor blisk at surge-like reversed flow conditions
Aeroelastic instabilities simulations on turbofan HP compressor blisk at surge-like reversed flow conditions Raphaël Murray Kungliga Tekniska Högskolan, SE-100 44 Stockholm, Sweden. This paper presents
More informationContrasting the Harmonic Balance and Linearized Methods for Oscillating-Flap Simulations
Contrasting the Harmonic Balance and Linearized Methods for Oscillating-Flap Simulations Guillaume Dufour, Frédéric Sicot, Guillaume Puigt European Center for Research and Advanced Training in Scientific
More informationFourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics
Fourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics K. Willcox and A. Megretski A new method, Fourier model reduction (FMR), for obtaining stable, accurate, low-order models
More informationManhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza
REPRESENTING PRESENCE OF SUBSURFACE CURRENT TURBINES IN OCEAN MODELS Manhar Dhanak Florida Atlantic University Graduate Student: Zaqie Reza 1 Momentum Equations 2 Effect of inclusion of Coriolis force
More informationGiven the water behaves as shown above, which direction will the cylinder rotate?
water stream fixed but free to rotate Given the water behaves as shown above, which direction will the cylinder rotate? ) Clockwise 2) Counter-clockwise 3) Not enough information F y U 0 U F x V=0 V=0
More informationPhysical Diffusion Cures the Carbuncle Phenomenon
Physical Diffusion Cures the Carbuncle Phenomenon J. M. Powers 1, J. Bruns 1, A. Jemcov 1 1 Department of Aerospace and Mechanical Engineering University of Notre Dame, USA Fifty-Third AIAA Aerospace Sciences
More informationConjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane
Conjugate Heat Transfer Simulation of Internally Cooled Gas Turbine Vane V. Esfahanian 1, A. Shahbazi 1 and G. Ahmadi 2 1 Department of Mechanical Engineering, University of Tehran, Tehran, Iran 2 Department
More informationMultiple Solutions and buffet of Transonic Flow over NACA0012 Airfoil
Fluid Dynamics and Co-located Conferences June 24-27, 23, San Diego, CA 3st AIAA Applied Aerodynamics Conference AIAA 23-326 Multiple Solutions and buffet of Transonic Flow over NACA2 Airfoil Juntao Xiong,
More informationGT Influence of Hot Streak Circumferential Length-Scale in Transonic Turbine Stage
Proceedings of ASME TURBO EXPO 24: International Gas Turbine & Aeroengine Congress & Exhibition June 4-7, 24, Vienna, Austria GT24-5337 Influence of Hot Streak Circumferential Length-Scale in Transonic
More informationImplicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference Method
40th Fluid Dynamics Conference and Exhibit 28 June - 1 July 2010, Chicago, Illinois AIAA 2010-4442 Implicit Large Eddy Simulation of Transitional Flow over a SD7003 Wing Using High-order Spectral Difference
More informationAssessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow
Assessment of Implicit Implementation of the AUSM + Method and the SST Model for Viscous High Speed Flow Simone Colonia, René Steijl and George N. Barakos CFD Laboratory - School of Engineering - University
More information[N175] Development of Combined CAA-CFD Algorithm for the Efficient Simulation of Aerodynamic Noise Generation and Propagation
The 32nd International Congress and Exposition on Noise Control Engineering Jeju International Convention Center, Seogwipo, Korea, August 25-28, 2003 [N175] Development of Combined CAA-CFD Algorithm for
More informationRANS Solutions Using High Order Discontinuous Galerkin Methods
RANS Solutions Using High Order Discontinuous Galerkin Methods Ngoc Cuong Nguyen, Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology, Cambridge, MA 2139, U.S.A. We present a practical
More information