• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.161-171
• /
• 1999

Aerodynamic sensitivity analysis codes are developed via the hand-differentiation using a direct differentiation method and an adjoint method respectively from discrete two-dimensional compressible Navier-Stokes equations. Unlike previous other researches, Baldwin-Lomax algebraic turbulence model is also differentiated by hand to obtain design sensitivities with respect to design variables of interest in turbulent flows. Discrete direct sensitivity equations and adjoint equations are efficiently solved by the same time integration scheme adopted in the flow solver routine. The required memory for the adjoint sensitivity code is greatly reduced at the cost of the computational time by allowing the large banded flux jacobian matrix unassembled. Direct sensitivity code results are found to be exactly coincident with sensitivity derivatives obtained by the finite difference. Adjoint code results of a turbulent flow case show slight deviations from the exact results due to the limitation of the algebraic turbulence model in implementing the adjoint formulation. However, current adjoint sensitivity code yields much more accurate sensitivity derivatives than the adjoint code with the turbulence eddy viscosity being kept constant, which is a usual assumption for the prior researches.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.04a
• /
• pp.57-65
• /
• 2009

An efficient and high-fidelity design approach for wing-body shape optimization is presented. Depending on the size of design space and the number of design of variable, aerodynamic shape optimization process is carried out via different optimization strategies at each design stage. In the first stage, global optimization techniques are applied to planform design with a few geometric design variables. In the second stage, local optimization techniques are used for wing surface design with a lot of design variables to maintain a sufficient design space with a high DOF (Degree of Freedom) geometric change. For global optimization, Kriging method in conjunction with Genetic Algorithm (GA) is used. Asearching algorithm of EI (Expected Improvement) points is introduced to enhance the quality of global optimization for the wing-planform design. For local optimization, a discrete adjoint method is adopted. By the successive combination of global and local optimization techniques, drag minimization is performed for a multi-body aircraft configuration while maintaining the baseline lift and the wing weight at the same time. Through the design process, performances of the test models are remarkably improved in comparison with the single stage design approach. The performance of the proposed design framework including wing planform design variables can be efficiently evaluated by the drag decomposition method, which can examine the improvement of various drag components, such as induced drag, wave drag, viscous drag and profile drag.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.66-72
• /
• 2000

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1102-1108
• /
• 2002

In this study, a topological design sensitivity of the ai. bearing surface (ABS) is suggested by using an adjoint variable method. The discrete form of the generalized lubrication equation based on a control volume formulation is used as a compatible condition. A residual function of the slider is considered as an equality constraint function, which represents the slider in equilibrium. The slider thickness parameters at all grid cells are chosen as design variables since they are the topological parameters determining the ABS shape. Then, a complicated adjoint variable equation is formulated to directly handle the highly nonlinear and asymmetric coefficient matrix and vector in the discrete system equation of air-lubricated slider bearings. An alternating direction implicit (ADI) scheme is utilized for the numerical calculation. This is an efficient iterative solver to solve large-scale problem in special band storage. Then, a computer program is developed and applied to a slider model of a sophisticated shape. The simulation results of design sensitivity analysis (DSA) are directly compared with those of FDM at the randomly selected grid cells to show the effectiveness of the proposed approach. The overall distribution of DSA results are reported, clearly showing the region on the ABS where special attention should be given during the manufacturing process.

• Proceedings of the KIEE Conference
• /
• 2008.07a
• /
• pp.655-656
• /
• 2008

This research presents a topology optimization to minimize eddy current loss maintaining mechanical robustness of structure part in electrical machines A design sensitivity equation for the topology optimization is derived by employing the discrete system equations combined with the adjoint variable method. As a numerical example, frame design of a C-core actuator is performed by the proposed method.

• Proceedings of the KIEE Conference
• /
• 2006.07b
• /
• pp.651-652
• /
• 2006

This research presents a topology optimization for manipulating the main heat flow in coupled magneto-thermal systems. The heat generated by eddy currents is considered in the design domain assuming an adiabatic boundary. For a practical optimization, the convection condition is considered in the topological process of the thermal field. Topology design sensitivity is derived by employing the discrete system equations combined with the adjoint variable method. As numerical examples, a simple iron and a C-core design heated-up by eddy currents demonstrate the strength of the proposed approach to solve the coupled problem.

• Proceedings of the KIEE Conference
• /
• 2007.07a
• /
• pp.138-139
• /
• 2007

This research provides machine designers with some intuition to consider both, magnetic and heat transfer effects. A topological multi-objective function includes magnetic energy and heat inflow rate to the system, which equals to the total heat dissipation by conduction and convection. For the thermal field regarding the heat inflow, introduced as a reaction force, topology design sensitivity is derived by employing discrete equations. The adjoint variable method is used to avoid numerous sensitivity evaluations. As a numerical example, a C-core design excited by winding current demonstrates the strength of the multi-physical approach.