This paper presents a review of the methodology, problems and progress in computational aeroacoustics(CAA). The nature, characteristics, and objectives of aeroacoustics problems are quite different from the commonly encountered CFD problems. In this paper, computational methods that are designed especially for CAA applications are introduced. The potential offered by CAA, the numerical issues which need to be given careful attention, and some of the recent progress in solving aeroacoustic problem are discussed

A coupled general circulation model (GCM) of the global atmosphere, oceans and lands is used for finding the future climate at times of doubled carbon dioxide concentration (DCDC) of the atmosphere. Two runs of the model were made in order to find the future change. Global changes at times of DCDC may be characterized by a global warmig of $1.4^{\circ}C$, a $3\%$ global precipitation increase, and an increase in the surface available water (SAW) over the global land among others. The estimated increase in SAW over the land implies that river discharge is likely to increase with increased chance of severe weather as a result of the future global warming.

Pressure-flow control characteristics of a commercially available cerebrospinal flow(CSF) control shunt valve was tested using fluid-solid interaction analysis. Pre-stress of the valve diaphragm(membrane) was computed for proper valve opening. The results were ir good agreements with the valve specification listed in the commercially available CSF control valve. The results of the study can be effectively used to design variety of CSF control shunt valves.

In this study, the three-dimensional blood flow within the sac of KTAH(Korean artificial heart) is simulated using fluid-structure interaction model. The numerical method employed in this study is the finite element commercial package ADINA. The thrombus formation is one of the most critical problems in KTAH. High fluid shear stress or stagnated flow are believed to be the main causes of these disastrous phenomenon. We solved the fluid-structure interaction between the 3D blood flow in the sac and the surrounding sac material. The sac material is assumed as linear elastic material and the blood as incompressible viscous fluid. Numerical solutions show that high shear stress region and stagnated flow are found near the upper part of the sac and near the comer of the outlet during diastole stage.

A Numerical study of laminar vortex-shedding past a circular cylinder has been performed widely by many researchers. Many factors, such as numerical technique and domain size, number and shape of grid, affected predicting vortex shedding and Strouhal number. In the present study, the effect of convection scheme, time discretization methods and grid dependence were investigated. The present paper presents the finite volume solution of unsteady flow past circular cylinder at Re=200, 400. The Strouhal number was predicted using UDS, CDS, Hybrid, Power-law, LUDS, QUICK scheme for convection term, implicit and crank-nicolson methods for time discretization. The grid dependence was investigated using H-type mesh and O-type mesh. It also studied that the effect of mesh size of the nearest adjacent grid of circular cylinder. The effect of convection scheme is greater than the effect of time discretization on predicting Strouhal. It has been found that the predicted Strouhal number changed with mesh size and shape.

Two-dimensional Navier-Stokes codes are used to simulate the shock tunnel in Konkuk university. In order to design experiments in impulse facilities properly and to interpret data from such facilities, it is necessary to understand how the flow approaches steady state. This is done by determining the transient flow field and flow establishment time around a given model. This will be accomplished by developing appropriate CFD codes which solve the Navier-Stokes equations, and simulating the starting process and resulting unsteady viscous flow phenomena. The starting process in a shock tunnel consists of multiple shock interactions and contact discontinuities, which are difficult to solve with the classical shock capturing schemes. A recently developed high resolution scheme is adapted for resolving the unsteady phenomena of those multiple shock interactions and contact surfaces during the starting process. The bifurcation phenomenon due to the interactions of the reflected shock from the end of the shock tube with the boundary layer generated by the incident shock becomes of particular interest. By comparing with the experiment results, the accuracy of the numerical analysis is validated and it is demonstrated that the properties which can hardly be obtained through the experiment can be estimated.

We study a conceptual numerical model on shock-vortex interaction setting an impulsive shock in a compressible vertex. Navier-Stokes equations are solved for the investigation of interactive structure and acoustic wave propagation. The rotationally symmetric vortex enforces two compression-expansion pairs resultantly forming a quadrupolar shape. These compressive and expansive waves cylindrically propagate to the far field and turn to acoustic waves. Using a fine uniform Cartesian grid system and a TVD-high resolution method, the flow data irl: precisely obtained to extend our interest to the sound source.

The Interactions of shock wave with turbulent boundary layers in high-speed flows cause complex flowfields which result in increased adverse pressure gradients, skin friction and temperature. Accurate and reliable prediction of such phenomena is needed in designing high-speed propulsion systems. Such analyses of the complex flowfields require sophisticated numerical scheme that can resolve interactions between shock wave and boundary layers accurately. Therefore the purpose of the present. article is to introduce an accurate and efficient mixed explicit-implicit generalized Galerkin finite element method. To demonstrate the validity of the theory and numerical procedure, several benchmark cases are investigated.

An aerodynamic design method has been developed by using a three-dimensional unstructured Euler code and an adjoint code with a discrete approach. The resulting adjoint code is applied to a wing design problem of super-sonic transport with a wing-body-nacelle configuration. Hicks-Henne shape functions are adopted far the surface geometry perturbation, and the elliptic equation method is employed fer the interior grid modification during the design process. Interior grid sensitivities are neglected except those for design parameters associated with nacelle translation. The Sequential Quadratic Programming method is used to minimize the drag with constraints on the lift and airfoil thickness. Successful design results confirm validity and efficiency of the present design method.

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.

Small-sized fans for cooling electronic components are usually produced by injection molding in unified form including blade and hub. However, in case of middle and large-sized fans, they deal ill various products with varied numbers of blade or stagger angles after designing a single blade by molding or aluminum die casting. At this occasion, it is necessary to study design method for high-performance axial fans can be operated with various numbers of blades or stagger angle using unique blade for specified conditions. Therefore, the goal in this study is developing a optimal design method which improves performance of axial fans within the large range of operating by single blade.

An iterative time marching procedure for solving incompressible internal flow has been applied to the flow around a rotor-stator. This procedure solves three-dimensional incompressible Reynolds-averaged Navier-Stokes equation on a moving, time-deforming, non-orthogonal body-fitted grid using second-order accurate schemes for the time derivatives and third/second-order accurate schemes for the spatial derivatives. To handle rotationg geometry, the multiblock technique is applied and the overall flow domain is subdivided into two blocks. In each block, a grid is generated and flowfield is solved independently of the other blocks. The boundary data for each block is provided by the neighboring blocks using bilinear interpolation technique.

Ground effect of three-dimensional wing is studied. LU-factored Implicit upwind TVD scheme and Baldwin-Lomax turbulence model are used for this calculation. To investigate ground effect, NACA 4415 wing at M=0.5 calculated. Two different angles of attack and three cases of flight height are calculated. As increasing angle of attack, the ground effect becomes strong. In case of NACA 4415 wing in ground effect, strength of wing tip vortex becomes stronger than that of free flight.

The three-dimensional unsteady compressible Full Navier-Stokes equation solver with sliding multi-block method has been applied to analyze three dimensional characteristics of the flow field and compression wave around the high speed train which Is entering into a tunnel. The numerical scheme of AF + ADI was used to efficiently solve Navier-Stokes equations in the curvilinear coordinate system. The vortex formation around the nose region was found and the generation of compression wave due to the blockage effects was observed ahead of the train in the form of plane wave. The three dimensional characteristics of the flow field compared to the analytic results were discussed in detail. The variation of pressure of tunnel wall surface and velocity profile of the train are identified as the train enters into a tunnel. The changes in aerodynamic forces and streamlines of each specific sections are also discussed and presented.

A method of two and three dimensional orthogonal grid generation with control of spacing by using the covariant Laplace equation is Presented. An important feature of the methodology is its ability to control effectively the grid spacing especially near the boundaries still maintaining good orthogonality in whole field. The method is based on the concept of decomposition of the global transformation into consecutive transformation of an approximate conformal mapping and au auxiliary orthogonal mapping to have linear and uncoupled equations. Control of cell spacing is based on the concept of reference arc length, and orthogonal correction is performed in the auxiliary domain. It is concluded that the methodology can successfully generate well controlled orthogonal grids around bodies of 2 and 3 dimensional configurations.

연소실 안으로 분출되는 스월 유동의 vortex breakdown mechanism 에 대한 연구를 하였다. 3 차원 유한 체적기법과 Runge-Kutta 시간 적분법이 적용되었으며, 난류모델은 dynamic large eddy simulation (DLES)이 적용되었다. 계산 시간의 효율성과 기억용량을 효과적으로 사용하기 위하여 message passing interface (MPI) 병렬계산 기법이 적용되었다. 스월 난류 유동에 있어서 vortex breakdown 거동을 가시적으로 표착 하였는데. 이는 스월 유동에 의한 난류 응력 증대, 난류 생성/소산율 증대 및 혼합율 증대에 대한 실험적 근거를 뒷받침하는 매우 중요한 결과이다. 또한 평균 속도와 레이놀스 응력에 대한 계산 결과도 실험 결과와 비교하였다.

Two-dimensional steady flow of two immiscible, compressible fluids are considered when the temperature of each layer is constant. Both upper and lower fluids are bounded by two horizontal rigid boundaries with symmetric obstruction of compact support at the tourer boundary. By using asymptotic method, we derive the forced K-dV equation governing interfacial wave. Various solutions and numerical results are presented.

An S-duct flow is subjected to an entrance flow of Mach 0.6. The duct turns $30^{\circ}$ and reverses its turn by $30^{\circ}$ followed by a straight section. Such an internal flow induces a secondary flow due to curvature effect. Goal of this paper is to show the sensitivity of outflow boundary conditions on the quality of numerical solutions as well as to show curvature effect on the flow field. The often-used Baldwin-Lomax turbulence model is shown to be less functional on the concave region when the secondary flow has its strong Influence.

Transient aerodynamic response of an airfoil to a rapidly deploying spoiler is numerically investigated using a two-dimensional turbulent compressible Navier-Stokes flow model. The spoiler moving relative to a stationary airfoil is treated by an overset grid bounded by a 'dynamic domain-dividing line' the concept of which is developed first..in this paper. The fluid-dynamic mechanism of the adverse lift due to the rapidly deploying spoiler is analyzed. Also the effect of spoiler deploying rate on the initial behavior of the aerodynamic response is expounded, which is of interest in view of active control technology and controller design for the spoiler. The results of present computation about the stationary as well as moving spoilers are relatively in good agreement with the existing experimental data.

Two and Three dimensional dynamic adaptation code is developed for transient computations. This code involves mesh refinement and coarsening to either add points in high gradient regions of flow or remove points where they are not needed, for high spatial accuracy. Temporary cell algorithm is used to maintain the original grid quality. To show the assessment of the accuracy and efficiency, two dimensional study and unsteady flows are computed. Also, three dimensional steady computations are made to assess the refinement using temporary cell algorithm. The result shows the high spatial accuracy primarily in discontinuity regions in steady and unsteady computation.

For the analysis of compressible flow with real gas effect, characteristic form of Roe's Riemann solver was derived again using real gas equation of state and variable specific heal and it was extended to multi component reactive system. From this study, it is known that some correction should be made for the use of existing numerical algorithm. 1) Sonic speed and characteristic variable should be corrected with real gas effect. 2) Roe's average was applicable only with the assumption of constant properties. 3) Artificial damping term and characteristic variables should be corrected but their influence may not be significant.

In this paper, numerical study has been done for the improvement of the superdetonative ram accelerator performance and for the design optimization of the system. The objective function to optimize the premixture composition is the ram tube length required to accelerate projectile from initial velocity $V_o$ to target velocity $V_e$. The premixture is composed of $H_2,\;O_2,\;N_2$ and the mole numbers of these species are selected at design variables. RSM(Response Surface Methodology) which is widely used for the complex optimization problems is selected as the optimization technique. In particular, to improve the non-linearity of the response and to consider the accuracy and efficiency of the solution, design space stretching technique has been applied. Separate sub-optimization routine is introduced to determine the stretching position and clustering parameters which construct the optimum regression model. Two step optimization technique has been applied to obtain the optimal system. With the application of stretching technique, we can perform system optimization with a small number of experimental points, and construct precise regression model for highly non-linear domain. The error to compared with analysis result is only $0.01\%$ and it is demonstrated that present method can be applied more practical design optimization problems with many design variables.

A two-dimensional/axisymmetric CSCM upwind flux difference splitting Wavier-Stokes method has been developed to study the finite rate chemically react-ing invisicd and viscous hypersonic flows over blunt-body. A upwind method was chosen due to its robustness in capturing the strong bow shock waves. For the nonequilibrium chemically reacting air, NS-I species conservation equations were strongly coupled with flowfield equations through convection and species production terms. The nonequilibrium wall pressure and heat transfer rate distributions along the vehicle were compared with those from equilibrium and perfect gas calculations. The nonequilibrium species distribution shows the reduced concentrations of O and N species when compared with equilibrium species distribution. The solutions resolved strong bow shock waves md heat transfer rate very accurately when compared with central difference schemes.

In this numerical work, three dimensional supersonic laminar flow and heat transfer of a blunt body(sphere-cone) at Mach 5 is simulated. The effects of angle of attack and the spin rate on the now and heat transfer are analysed. To solve the three dimensional compressible Wavier-Stokes equation, a finite volume method with the modified LDFSS scheme is employed for spatial discretization, and a point SGS implicit method is used for time integration. It is found that the heat transfer rate increases at the windward side and decreases at the leeward side with the angle of attack. The heat transfer rate at all surfaces slightly increases with the spin rate.

An efficient Newton-GMRES algorithm is presented for computing two-dimensional steady compressible viscous flows on unstructured hybrid meshes. The scheme is designed on cell-centered finite volume method which accepts general polygonal meshes. Steady-state solution is obtained with pseudo-transient continuation strategy. The preconditioned, restarted general minimum residual(GMRES) method is employed in matrix-free form to solve the linear system arising at each Newton iteration. The incomplete LU fartorization is employed for the preconditioning of linear system. The Spalart-Allmars one equation turbulence model is fully coupled with the flow equations to simulate turbulence effect. The accuracy, efficiency and robustness of the presently developed method are demonstrated on various test problems including laminar and turbulent flows over flat plate and airfoils.

A well-known pressure correction method, a SIMPLE algorithm is extended to treat compressible flows. Collocated grids are used and density is linked to pressure via an equation of state. The influence of pressure on density in the case of compressible flows is implicitly incorporated into the extended SIMPLE algorithm. The first-order Upwind and high-order Quick scheme are compared with respect to an accuracy and convergence time at all speeds. The extended method is verified on a number of test cases and the results we compared with other numerical results available in the literature. The calculated results show that the Quick scheme improves accuracy at all speed and also reduces the calculation time at supersonic flows, compared with the Upwind scheme.