• Journal of applied mathematics & informatics
• /
• v.19 no.1_2
• /
• pp.215-229
• /
• 2005

In this paper, constraint aggregation is combined with the adjoint and multiple shooting strategies for optimal control of differential algebraic equations (DAE) systems. The approach retains the inherent parallelism of the conventional multiple shooting method, while also being much more efficient for large scale problems. Constraint aggregation is employed to reduce the number of nonlinear continuity constraints in each multiple shooting interval, and its derivatives are computed by the adjoint DAE solver DASPKADJOINT together with ADIFOR and TAMC, the automatic differentiation software for forward and reverse mode, respectively. Numerical experiments demonstrate the effectiveness of the approach.

• 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.

• 한국전산유체공학회:학술대회논문집
• /
• 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.

• Nuclear Engineering and Technology
• /
• v.52 no.3
• /
• pp.461-468
• /
• 2020

The improvements of the DeCART/MUSAD code system for uncertainty analysis of HTGR neutronic parameters are presented in this paper. The function for quantifying an uncertainty of critical-spectrumweighted few group cross section was implemented using the generalized adjoint B₁ equation solver. Though the changes between the infinite and critical spectra cause a considerable difference in the contribution by the graphite scattering cross section, it does not significantly affect the total uncertainty. To reduce the number of iterations of the generalized adjoint transport equation solver, the generalized adjoint B₁ solution was used as the initial value for it and the number of iterations decreased to 50%. To reflect the implicit uncertainty, the correction factor was derived with the resonance integral. Moreover, an additional correction factor for the double heterogeneity was derived with the effective cross section of the DH region and it reduces the difference from the complete uncertainty. The code system was examined with the MHTGR-350 Ex.II-2 3D core benchmark. The k_eff uncertainty for Ex.II-2a with only the fresh fuel block was similar to that of the block and the uncertainty for Ex.II-2b with the fresh fuel and the burnt fuel blocks was smaller than that of the fresh fuel block.

• 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.

• International Journal of Automotive Technology
• /
• v.7 no.3
• /
• pp.283-288
• /
• 2006

A mathematical methodology is proposed for designing nozzle hole shapes producing controlled geometric cavitation. The proposed methodology uses an unstructured RANS flow solver, with the ability to compute sensitivity derivatives via an adjoint algorithm. The adjoint formulation for the N-S equations is presented while variation of the turbulence viscosity is not taken into account during the geometry modifications. The sensitivities are calculated in a mode independently of the shape parameterisation. The method is used to develop and evaluate conceptual shapes for nozzle hole cavitation reduction. The localized region at the hole inlet producing cavitation, is parameterised using its radius of curvature, while a cost function is formulated to eliminate the negative pressures present at this location. Sensitivity derivatives are used to assess the dependence of the localized region on the minimum pressure, and to drive the geometry to the targeted shape. The results show that the computer model can provide nozzle hole entry shapes that produce predefined flow characteristics, and thus can be used as an inverse design tool for nozzle hole cavitation control.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.2 no.1
• /
• pp.27-39
• /
• 1998

We describe a fast numerical method for non-separable elliptic equations in self-adjoin form on irregular adaptive domains. One of the most successful results in numerical PDE is developing rapid elliptic solvers for separable EPDEs, for example, Fourier transformation methods for Poisson problem on a square, however, it is known that there is no rapid elliptic solvers capable of solving a general nonseparable problems. It is the purpose of this paper to present an iterative solver for linear EPDEs in self-adjoint form. The scheme discussed in this paper solves a given non-separable equation using a sequence of solutions of Poisson equations, therefore, the most important key for such a method is having a good Poison solver. High performance is achieved by using a fast high-order adaptive Poisson solver which requires only about 500 floating point operations per gridpoint in order to obtain machine precision for both the computed solution and its partial derivatives. A few numerical examples have been presented.

• International Journal of Aeronautical and Space Sciences
• /
• v.8 no.1
• /
• pp.66-78
• /
• 2007

A multi-point aerodynamic shape optimization technique has been developed for helicopter rotor blades in hover based on a continuous adjoint method on unstructured meshes. The Euler flow solver and the continuous adjoint sensitivity analysis were formulated on the rotating frame of reference. The 'objective function and the sensitivity were obtained as a weighted sum of the values at each design point. The blade section contour was modified by using the Hicks-Henne shape functions. The mesh movement due to the blade geometry change was achieved by using a spring analogy. In order to handle the repeated evaluation of the design cycle efficiently, the flow and adjoint solvers were parallelized based on a domain decomposition strategy. A solution-adaptive mesh refinement technique was adopted for the accurate capturing of the wake. Applications were made to the aerodynamic shape optimization of the Caradonna-Tung rotor blades and the UH-60 rotor blades in hover.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.1
• /
• pp.1-10
• /
• 2005

An aerodynamic shape optimization technique has been developed for helicopter rotor blades in hover based on a continuous adjoint method on unstructured meshes. The Euler flow solver and the continuous adjoint sensitivity analysis were formulated on the rotating frame of reference for hovering rotor blades. In order to handle the repeated evaluation of the design cycle efficiently, the flow and adjoint solvers were parallelized using a domain decomposition strategy. A solution-adaptive mesh refinement technique was adopted for the accurate capturing of the tip vortex. Applications were made for the aerodynamic shape optimization of Caradonna-Tung rotor blades and UH60 rotor blades in hover. The results showed that the present method is an effective tool to determine optimum aerodynamic shapes of rotor blades requiring less torque while maintaining the desired thrust level.

• Nuclear Engineering and Technology
• /
• v.51 no.1
• /
• pp.13-30
• /
• 2019

This paper presents the recent improvements in the DeCART code for HTGR analysis. A new 190-group DeCART cross-section library based on ENDF/B-VII.0 was generated using the KAERI library processing system for HTGR. Two methods for the eigen-mode adjoint flux calculation were implemented. An azimuthal angle discretization method based on the Gaussian quadrature was implemented to reduce the error from the azimuthal angle discretization. A two-level parallelization using MPI and OpenMP was adopted for massive parallel computations. A quadratic depletion solver was implemented to reduce the error involved in the Gd depletion. A module to generate equivalent group constants was implemented for the nodal codes. The capabilities of the DeCART code were improved for geometry handling including an approximate treatment of a cylindrical outer boundary, an explicit border model, the R-G-B checker-board model, and a super-cell model for a hexagonal geometry. The newly improved and implemented functionalities were verified against various numerical benchmarks such as OECD/MHTGR-350 benchmark phase III problems, two-dimensional high temperature gas cooled reactor benchmark problems derived from the MHTGR-350 reference design, and numerical benchmark problems based on the compact nuclear power source experiment by comparing the DeCART solutions with the Monte-Carlo reference solutions obtained using the McCARD code.