Engineering Papers⌕ Search

Engineering topics

Baysal, Oktay

Publications and source records attributed to Baysal, Oktay.

At least 37 records · Page 2

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A discrete sensitivity analysis algorithm had previously been developed and applied to two-dimensional aerodynamic optimization problems, where the computational domains were discretized by using single grids. The sparse, unsymmetric systems of linear equations resulting from this algorithm were solved by a direct matrix inversion matrix. However, for large two-dimensional problems and, practically, all three-dimensional problems, direct inversion methods become inapplicable, primarily due to the prohibitive computer storage needed. In an attempt to alleviate such hindrances, the sensitivity analysis with domain decomposition (SADD) scheme was developed. This scheme divides the computational domain into smaller and nonoverlapping subdomains (multiblock grids) that are solved separately. Then, the final solution is constructed from the subdomain solutions. As the number of grid points in the interface boundaries of the subdomains becomes large, the computer memory required to store the effective coefficient matrix of these interface points starts to increase. Presented in this Technical Note is a preconditioned iterative procedure to overcome this particular problem.

Eleshaky, Mohamed E.↗

Three-dimensional aerodynamic shape optimization of supersonic delta wings

A recently developed three-dimensional aerodynamic shape optimization procedure AeSOP(sub 3D) is described. This procedure incorporates some of the most promising concepts from the area of computational aerodynamic analysis and design, specifically, discrete sensitivity analysis, a fully implicit 3D Computational Fluid Dynamics (CFD) methodology, and 3D Bezier-Bernstein surface parameterizations. The new procedure is demonstrated in the preliminary design of supersonic delta wings. Starting from a symmetric clipped delta wing geometry, a Mach 1.62 asymmetric delta wing and two Mach 1. 5 cranked delta wings were designed subject to various aerodynamic and geometric constraints.

Burgreen, Greg W.↗

Improving the efficiency of aerodynamic shape optimization

The computational efficiency of an aerodynamic shape optimization procedure that is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid-point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit methodology to calculate the highly converged flow solutions that are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. Practically identical optimization results are obtained that are independent of the method used to represent the surface. A substantial factor of 8 decrease in computational time for the optimization process is achieved by implementing both of the design procedure improvements.

Burgreen, Greg W.↗

Design of 3-D Nacelle near Flat-Plate Wing Using Multiblock Sensitivity Analysis (ADOS)

One of the major design tasks involved in reducing aircraft drag is the integration of the engine nacelles and airframe. With this impetus, nacelle shapes with and without the presence of a flat-plate wing nearby were optimized. This also served as a demonstration of the 3-D version of the recently developed aerodynamic design optimization methodology using sensitivity analysis, ADOS. The required flow analyses were obtained by solving the three-dimensional, compressible, thin-layer Navier-Stokes equations using an implicit, upwind-biased, finite volume scheme. The sensitivity analyses were performed using the preconditioned version of the SADD scheme (sensitivity analysis on domain decomposition). In addition to demonstrating the present method's capability for automatic optimization, the results offered some insight into two important issues related to optimizing the shapes of multicomponent configurations in close proximity. First, inclusion of the mutual interference between the components resulted in a different shape as opposed to shaping an isolated component. Secondly, exclusion of the viscous effects compromised not only the flow physics but also the optimized shapes even for isolated components.

Eleshaky, Mohamed E.↗

Dynamic-overlapped-grid simulation of aerodynamically determined relative motion

Currently, there is a need to develop a means of analyzing and studying unsteady flowfields which involve multiple component configurations with at least one of the components in relative motion with respect to the others. Two of the important phenomena that such analyses can help to understand are the unsteady aerodynamic interference and the boundary-induced component of the flowfield. With this motivation, a computational method is developed which couples the governing equations of the unsteady flowfield and the rigid-body dynamics in six degrees-of-freedom. These equations are solved on composite meshes of overlapped subdomain grids which can move with respect to each other. Initially, several measures that reduce the numerical error are studied and compared with the exact solution of a moving normal shock in a tube. It is concluded that a second-order accurate method, for spatial and temporal discretizations as well as for the moving subdomain interpolations, is needed as a minimum measure. Furthermore, the CFL numbers should be restricted to unity. Then, the method is used to simulate the flowfield history and predict the aerodynamically determined trajectory of a store dropped from its initial position under a wing.

Yen, Guan-Wei↗

Airfoil shape optimization using sensitivity analysis on viscous flow equations

An aerodynamic shape optimization method has previously been developed by the authors using the Euler equations and has been applied to supersonic-hypersonic nozzle designs. This method has also included a flowfield extrapolation (or flow prediction) method based on the Taylor series expansion of an existing CFD solution. The present paper reports on the extension of this method to the thin-layer Navier-Stokes equations in order to account for the viscous effects. Also, to test the method under highly nonlinear conditions, it has been applied to the transonic flows. Initially, the success of the flow prediction method is tested. Then, the overall method is demonstrated by optimizing the shapes of two supercritical transonic airfoils at zero angle of attack. The first one is shape optimized to achieve a minimum drag while obtaining a lift above a specified value. Whereas, the second one is shape optimized for a maximum lift while attaining a drag below a specified value. The results of these two cases indicate that the present method can produce successfully optimized aerodynamic shapes.

Eleshaky, Mohamed E.↗

Aerodynamic shape optimization using preconditioned conjugate gradient methods

In an effort to further improve upon the latest advancements made in aerodynamic shape optimization procedures, a systematic study is performed to examine several current solution methodologies as applied to various aspects of the optimization procedure. It is demonstrated that preconditioned conjugate gradient-like methodologies dramatically decrease the computational efforts required for such procedures. The design problem investigated is the shape optimization of the upper and lower surfaces of an initially symmetric (NACA-012) airfoil in inviscid transonic flow and at zero degree angle-of-attack. The complete surface shape is represented using a Bezier-Bernstein polynomial. The present optimization method then automatically obtains supercritical airfoil shapes over a variety of freestream Mach numbers. Furthermore, the best optimization strategy examined resulted in a factor of 8 decrease in computational time as well as a factor of 4 decrease in memory over the most efficient strategies in current use.

Burgreen, Greg W.↗

Preconditioned domain decomposition scheme for three-dimensional aerodynamic sensitivity analysis

A preconditioned domain decomposition scheme is introduced for the solution of the 3D aerodynamic sensitivity equation. This scheme uses the iterative GMRES procedure to solve the effective sensitivity equation of the boundary-interface cells in the sensitivity analysis domain-decomposition scheme. Excluding the dense matrices and the effect of cross terms between boundary-interfaces is found to produce an efficient preconditioning matrix.

Eleshaky, Mohammed E.↗

Viscous simulation method for unsteady flows past multicomponent configurations

The present algorithm for the numerical simulation of flows about complex configurations (whose multiple and nonsimilar components have arbitraty geometries) employs a hybridization of the domain decomposition techniques for grid generation as well as to reduce computer-memory requirements. A fully vectorized, finite-volume, upwind-biased, approximately factored multigrid method is used to solve 3D Reynolds-averaged unsteady and compressible Navier-Stokes equations simulating supersonic flows past an ogive-nose-cylinder near or within a cavity. The time-averaged surface pressures obtained compare favorably with the wind tunnel data.

Fouladi, Kamran↗

Tangential, semisubmerged, and internal store carriage and separation

Accounts are presented of tangential, semisubmerged, and internal store carriage drag, as well as of the results to date of stores-separation investigations employing both computational and experimental methods. It is demonstrated that CFD can contribute to the data base needed for internal stores carriage design and trade studies. Attention is given to the case of an internally carried store separating from its bay at supersonic speeds, where various kinds of interference are encountered. A code is developed for the simulation of these unsteady flows through time-accurate computations; computation results are noted to closely reproduce experiment data.

Wilcox, Floyd J., Jr.↗

Transonic solutions of a wing/pylon/finned store using hybrid domain decomposition

Transonic Euler calculations about a complex multicomponent configuration are presented. The 3D Euler equations are solved utilizing an upwind-biased, alternating direction implicit, approximately factored, multigrid algorithm. Computational results are compared to experimental data of the finned store in a carriage position.

Newman, James C., III↗

Computing high-speed flows past an oscillating cylinder near a vertical wall

A computational method to simulate unsteady flows involving moving rigid boundaries and interference has been developed. The method is used to solve inviscid equations governing the fluid flow and the dynamic equations governing the motion of rigid bodies. A second-order accurate, upwind-biased, and alterating-direction-implicit method is used to solve the governing equations of the flow. A kinematic domain decomposition (KDD) procedure is extended to treat 3D problems with a high degree of accuracy and generality. The method under consideration is applied to both transonic and supertransonic flows. Both cases involve flow past a cylinder which is forced to pitch sinusoidally near a vertical wall. Benefits of the proposed approach include accurate calculation of the flow around 3D moving multiple bodies with interference; reduction of a numerical error; in particular, the dispersion error which strongly affects wave propagation; and minimization of the phase error which is accumulated according to the time advance procedure.

Yen, Guan-Wei↗

Navier-Stokes computations of cavity aeroacoustics with suppression devices

Effectiveness of two devices to suppress the cavity acoustics was computationally investigated. Two dimensional, computational simulations were performed for the transonic, turbulent flows past a cavity, which was first equipped with a rear face ramp and then with a spoiler. The Reynolds-averaged, unsteady, compressible, full Navier-Stokes equations were solved time accurately by a second-order accurate, implicit, upwind, finite-volume method. The effect of turbulence was included through a Baldwin-Lomax model with modifications for the multiple-wall effects and for the highly vortical flow with a shear layer. The results included instantaneous and time-averaged flow properties, and time-series analyses of the pressure inside the cavity, which compared favorably with the available experimental data. These results were also contrasted with the computed aeroacoustics of the same cavity (length-to-depth ratio of 4.5), but without a device, to demonstrate the suppression effectiveness.

Baysal, Oktay↗

Improving the efficiency of aerodynamic shape optimization procedures

The computational efficiency of an aerodynamic shape optimization procedure which is based on discrete sensitivity analysis is increased through the implementation of two improvements. The first improvement involves replacing a grid point-based approach for surface representation with a Bezier-Bernstein polynomial parameterization of the surface. Explicit analytical expressions for the grid sensitivity terms are developed for both approaches. The second improvement proposes the use of Newton's method in lieu of an alternating direction implicit (ADI) methodology to calculate the highly converged flow solutions which are required to compute the sensitivity coefficients. The modified design procedure is demonstrated by optimizing the shape of an internal-external nozzle configuration. A substantial factor of 8 decrease in computational time for the optimization process was achieved by implementing both of the design improvements.

Burgreen, Greg W.↗

Aerodynamic shape optimization via sensitivity analysis on decomposed computational domains

Direct and iterative method considered to be most applicable to large systems of linear equations arising in discrete sensitivity analysis are assessed. Based on a single-domain grid, computations are performed using a banded matrix solver and an iterative solver, the generalized minimum residual (GMRES) method. The banded matrix solver is found to be generally the most economical method for those applications where the number of right-hand sides is large (i.e., a large number of design variables or a large number of adjoint vectors). For systems of equations that are too large to be solved by direct methods, an approach is proposed whereby the computational domain is divided into small subdomains, and each subdomain is solved separately.

Eleshaky, Mohamed E.↗

Aerodynamic sensitivity analysis methods for the compressible Euler equations

This study presents a mathematical formulation developed for aerodynamic sensitivity coefficients based on a discretized form of the compressible 2D Euler equations. A brief motivating introduction to the aerodynamic sensitivity analysis and the reasons behind an integrated flow/sensitivity analysis for design algorithms are presented. The finite difference approach and the quasi-analytical approach are used to determine the aerodynamic sensitivity coefficients. A new flow prediction concept, which is an outcome of the direct method in the quasi-analytical approach, is developed and illustrated with an example. Surface pressure coefficient distributions of a nozzle-afterbody configuration obtained from the predicted flowfield solution are compared successfully with their corresponding values obtained from a flowfield analysis code and the experimental data.

Baysal, Oktay↗