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

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.67-73
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler, thin layer Navier-Stokes and full Navier-Stokes ones. are solved using implicit LU-ADI decomposition scheme. The gradient projection method with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.84-85
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for a large scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady computation of the incompressible Navier-Stokes equation, and Smagorinsky LES(Large Eddy Simulation) model is chosen for turbulent flow computation. Both METIS and MPI library are used for domain partitioning and data communication between processors respectively. Tiburon of Hyundai-motor is chosen as the computational model at $Re=7.5{\times}10^{5}$, which is based on the car height. It is confirmed that the drag under road condition is smaller than that of wind tunnel condition.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.171-176
• /
• 1999

In the present paper, three types of the flow solvers were used to investigate the influence on the airfoil shape optimization. The adopted equations, i.e., Euler , thin layer Navier- Stokes and full Navier-Stokes ones, are solved using implicit LU-ADI decomposition scheme. The feasible direction algorithm with the sinusoidal function was used as an optimization algorithm. The present numerical method was applied to the drag minimization problems under the initial shape of NACA0012 airfoils.

• Journal of Mechanical Science and Technology
• /
• v.14 no.8
• /
• pp.867-878
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using three-dimensional Euler equations. Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. Numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, length and the existence of a tunnel on the crossing event. Unsteady aerodynamic loads-a side force and a drag force-acting on the train during the crossing are numerically predicted and analyzed. The side force mainly depends on the nose-shape, and the drag force depends on tunnel existence. Also. a push-pull (i.e.impluse force) force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary to evaluate the stability of high speed multi-car trains. The results also indicate the effectiveness of the present numerical method for simulating the unsteady flow fields induced by bodies in relative motion.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.599-604
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other at the speed of 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using the three-dimensional Euler equations. The Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena properly. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. The numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, the train length and the existence of tunnel when the crossing event occur. Unsteady aerodynamic loads side force and drag force-acting on the train during the crossing are numerically predicted and anlayzed. It is found that the strength of the side force mainly depends on the nose-shape, and that of drag force on tunnel existence. And it is observed that the push-pull like impulsive force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary for the evaluation of the stability of the high speed multi-car train. The results also indicate the effectiveness of the present numerical method for the simulation of unsteady flow field induced by the bodies in the relative motion.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.27 no.8
• /
• pp.1182-1190
• /
• 2003

A parallelized FEM code based on domain decomposition method has been recently developed for large-scale computational fluid dynamics. A 4-step splitting finite element algorithm is adopted for unsteady flow computation of the incompressible Navier-Stokes equation, and Smagorinsky LES model is chosen for turbulent flow computation. Both METIS and MPI Libraries are used for domain partitioning and data communication between processors, respectively. Tiburon model of Hyundai Motor Company is chosen as the computational model at Re=7.5 $\times$ 10$^{5}$ , which is based on the car height. The calculation is carried out under both the wind tunnel condition and the road condition using IBM SP parallel architecture at KISTI Super Computing Center. Compared with the existing experimental data, both the velocity and pressure fields are predicted reasonably well and the drag coefficient is in good agreement. Furthermore, it is confirmed that the drag under the road condition is smaller than that under the wind-tunnel condition.

• Journal of Power System Engineering
• /
• v.17 no.6
• /
• pp.95-101
• /
• 2013

A stereo PIV measurements in a circulating water channel has been performed to investigate the skin friction reduction mechanism of the outer-layer vertical blades first devised by Hutchins. In a recent PIV measurement study, considerable skin friction reduction was achieved as much as 2.73%~7.95% by outer-layer vertical blades array. In the present study, the influence of vertical blades array upon the characteristics of the turbulent coherent structures was analyzed by proper orthogonal decomposition method. It is observed that the vortical structures are cut and deformed by blades array and also the turbulent intensity and the Reynolds stress were weakened by the blades. These phenomena strongly associate the skin-friction drag reduction mechanism in the turbulent boundary layer flow.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.55-58
• /
• 2002

Parallel finite element code has been recently developed for the analysis of the incompressible Wavier-Stokes equations using domain decomposition method. Metis and MPI libraries are used for the domain partitioning of an unstructured mesh and the data communication between sub-domains, respectively. For unsteady computation of the incompressible Navier-Stokes equations, 4-step splitting method is combined with P1P1 finite element formulation. Smagorinsky and dynamic model are implemented for the simulation of turbulent flows. For the validation performance-estimation of the developed parallel code, three-dimensional Laplace equation has been solved. It has been found that the speed-up of 40 has been obtained from the present parallel code fir the bench mark problem. Lastly, the turbulent flows around the MIRA model and Tiburon model have been solved using 32 processors on IBM SMP cluster and unstructured mesh. The computed drag coefficient agrees better with the existing experiment as the mesh resolution of the region increases, where the variation of pressure is severe.

• Korean Journal of Fisheries and Aquatic Sciences
• /
• v.27 no.5
• /
• pp.605-612
• /
• 1994

The methods to determine the hydrodynamic coefficients for the fixed type artificial reefs which were constructed to control ecological system in coastal waters are compared and discussed by model test results. To calculate the wave forces, least square method show good agreement with the experimental results and more stability than maximum force component method or Fourier decomposition method. This modified least square method of weighting the square of measured force turned out to be the most feasible method for maximum force. Using the feasible method, hydrodynamic characteristics for artificial reefs on uniform slopes offshore and breaking zone were studied. They were properly related to Keulegan-Carpenter's number and found larger than previous results. Wave force coefficients for artificial reefs around breaking zone were distributed from 1.5 to 2.5, and the mean value was 2.0. Drag force components were more in evidence than inertia force in maximum force which is important parameter to evaluate stability for high-permeability structures. A formula for the calculation of the maximum force for artificial reefs design is proposed, using structural dimension, water particle velocity and Keulegan-Carpenter's number.