• Korean Journal of Mathematics
• /
• v.25 no.3
• /
• pp.405-435
• /
• 2017

We deal with a sensitivity analysis of an optimal shape control problem for the stationary Navier-Stokes system. A two-dimensional channel flow of an incompressible, viscous fluid is examined to determine the shape of a bump on a part of the boundary that minimizes the viscous drag. By using the material derivative method and adjoint variables for a shape sensitivity analysis, we derive the shape gradient of the design functional for the model problem.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.1
• /
• pp.20-27
• /
• 2001

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. 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. The capability of the present sensitivity codes to treat complex geometry is successfully demonstrated by analyzing the flows over multi-element airfoils on Chimera overlaid grid systems.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.28 no.1
• /
• pp.9-18
• /
• 2015

In this paper, we verify the optimal topology design for heat conduction problems in steady stated which is obtained numerically using the adjoint design sensitivity analysis(DSA) method. In adjoint variable method(AVM), the already factorized system matrix is utilized to obtain the adjoint solution so that its computation cost is trivial for the sensitivity. For the topology optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of the structure and the allowable volume, respectively. For the experimental validation of the optimal topology design, we compare the results with those that have identical volume but designed intuitively using a thermal imaging camera. To manufacture the optimal design, we apply a simple numerical method to convert it into point cloud data and perform CAD modeling using commercial reverse engineering software. Based on the CAD model, we manufacture the optimal topology design by CNC.

• 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 Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.119-126
• /
• 2004

A continuum-based design sensitivity analysis and topology optimization methods are developed for power flow analysis. Efficient adjoint sensitivity analysis method is employed and further extended to topology optimization problems. Young's moduli of all the finite elements are selected as design variables and parameterized using a bulk material density function. The objective function and constraint are an energy compliance of the system and an allowable volume fraction, respectively. A gradient-based optimization, the modified method of feasible direction, is used to obtain the optimal material layout. Through several numerical examples, we notice that the developed design sensitivity analysis method is very accurate and efficient compared with the finite difference sensitivity. Also, the topology optimization method provides physically meaningful results. The developed is design sensitivity analysis method is very useful to systematically predict the impact on the design variations. Furthermore, the topology optimization method can be utilized in the layout design of structural systems.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.04a
• /
• pp.335-342
• /
• 2002

A continuum-based design sensitivity analysis (DSA) method fur non-shape problems is developed for geometrically nonlinear elastic structures. The non-shape problem is characterized by the design variables that are not associated with the domain of system like sizing, material property, loading, and so on. Total Lagrangian formulation with the Green-Lagrange strain and the second Piola-Kirchhoff stress is employed to describe the geometrically nonlinear structures. The spatial domain is discretized using the 4-node isoparametric plane stress/strain elements. The resulting nonlinear system is solved using the Newton-Raphson iterative method. To take advantage of the derived analytical sensitivity In topology optimization, a fast and efficient design sensitivity analysis method, adjoint variable method, is employed and the material property of each element is selected as non-shape design variable. Combining the design sensitivity analysis method and a gradient-based design optimization algorithm, an automated design optimization method is developed. The comparison of the analytical sensitivity with the finite difference results shows excellent agreement. Also application to the topology design optimization problem suggests a very good insight for the layout design.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.10a
• /
• pp.32-38
• /
• 2000

Wing flap deflection angles of a supersonic transport are optimized to improve transonic cruise performance. For this end, a numerical optimization method is adopted using a three-dimensional unstructured Euler code and a discrete adjoint code. Deflection angles of ten flaps; five for leading edge and five fur railing edge, are employed as design variables. The elliptic equation method is adopted for the interior grid modification during the design process. Interior grid sensitivities are neglected for efficiency. Also tested is the validity of the approximate gradient evaluation method for the present design problem and found that it is applicable for loading edge flap design in cases of no shock waves on the wing surface. The BFGS method is used to minimize the drag with constraints on the lift and upper surface Mach numbers. Two design examples are conducted; one is leading edge flap design, and the other is simultaneous design of leading edge and trailing edge flaps. The latter gave a smaller drag than the former by about two counts.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• v.3B no.1
• /
• pp.23-31
• /
• 2003

A novel 3D shape optimization algorithm is presented for electromagnetic devices carry-ing eddy current. The algorithm integrates the 3D finite element performance analysis and the steepest descent method with design sensitivity and mesh relocation method. For the design sensitivity formula, the adjoint variable vector is defined in complex form based on the 3D finite element method for eddy current problems. A new 3D mesh relocation method is also proposed using the deformation theory of the elastic body under stress to renew the mesh as the shape changes. The design sensitivity f3r the sur-face nodal points is also systematically converted into that for the design variables for the parameterized optimization application. The proposed algorithm is applied to the optimum design of the tank shield model of the transformer and the effectiveness is proved.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.37-43
• /
• 2010

In this paper, a reliability-based design optimization is developed for the topology design of linear structures using a performance measure approach. Spatial domain is discretized using three dimensional Reissner-Mindlin plate elements and design variable is taken as the material property of each element. A continuum based adjoint variable method is employed for the efficient computation of sensitivity with respect to the design and random variables. The performance measure approach of RBDO is employed to evaluate the probabilistic constraints. The topology optimizationproblem is formulated to have probabilistic displacement constraints. The uncertainties such as material property and external loads are considered. Numerical examples show that the developed topology optimization method could effectively yield a reliable design, comparing with the other methods such as deterministic, safety factor, and worst case approaches.