• 한국전산유체공학회:학술대회논문집
• /
• 2004.10a
• /
• pp.145-150
• /
• 2004

In this study the numerical analysis on staging flow with forward ejector is conducted. The forward ejector plays a vital role in staging, which jets out from aftbody. This staging environment needs careful flow analysis for securing staging safety Present study investigates the steady inviscid staging flow phenomena with variation of separation distance. The performance index is forebody base pressure coefficients. The three dominant flow phenomena are observed according to separation distance which could be told as impinging stage, cavity vortex dominancy stage, and pure base flow characteristics stage. Impinging stage shows high thrust for forebody as one might think. However, important point is that cavity vortex dominancy stage can be more favorable for separation than impinging stage as one simply think in certain separation distance.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.95-107
• /
• 1999

The spatial error due to the non-conservative interpolation become first-order when second-order conservative schemes are used, discontinuities are located away from the overlapped regions, and if the length of the overlapped region is not proportional to the grid spacing. Therefore, the solution accuracy is ensured if two domains overlap each other with a fixed grid point and the interpolation is occurred in smooth flow regions. To validate the spatial and temporal accuracy due to the non-conservative interpolation, inviscid and viscous problems are tested.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.48-58
• /
• 1999

The Characteristic Flux Difference Splitting (CFDS) scheme designed to adapt the characteristic boundary conditions at the wall and inflow/outflow boundary planes satisfies Roe's property U, although the CFDS Jacobian matrix is decomposed by a product of elaborate transformation matrices and explicit eigenvalue matrix. When the CFDS algorithm, thus a variant of Roe's scheme, is applied straightforwardly to hypersonic flows over a blunt body, the strong bow shock gradually breaks down near the stagnation point. This numerical instability is widely observed by many researchers employing flux-difference method, known in the literature as the carbuncle phenomenon. Many remedies have been proposed and resulted in partial cures. When the idea of Sanders et al. which identifies the minimum eigenvalues near the discontinuity present is applied to CFDS method, it is shown that the instability problem can be controlled successfully. A few flux splitting methods have also been tested and results are compared against the Nakamori's Mach 8 blunt body flow.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.11a
• /
• pp.42-47
• /
• 1999

In this study, we adopt the point symmetric Gauss-Seidel relaxation algorithm to obtain the steady state solution of the Navier-Stokes equations for the thermal and chemical nonequilibrium flow of air. All of the inviscid, viscous flux Jacobians and thermochemical source Jacobians are included in the implicit part Numerical simulation is performed for the thermal and chemical nonequilibrium flow over blunt body and computational results are presented. The convergence history and CPU time of the present computation are compared with the LU-SGS scheme which employs the approximate Jacobians.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.294-301
• /
• 2011

A two-dimensional hybrid flaw solver has been developed for the accurate and efficient simulation of steady and unsteady flaw fields. The flow solver was cast to accommodate two different topologies of computational meshes. Triangular meshes are adopted in the near-body region such that complex geometric configurations can be easily modeled, while adaptive Cartesian meshes are, utilized in the off-body region to resolve the flaw more accurately with less numerical dissipation by adopting a spatially high-order accurate scheme and solution-adaptive mesh refinement technique. A chimera mesh technique has been employed to link the two flow regimes adopting each mesh topology. Validations were made for the unsteady inviscid vol1ex convection am the unsteady turbulent flaws over an NACA0012 airfoil, and the results were compared with experimental and other computational results.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.458-464
• /
• 2011

A new method identifies coupled fluid-structure system with a reduced set of state variables is presented. Assuming that the structural model is known a priori either from an analysis or a test and using linear transformations between structural and aeroelastic states, it is possible to deduce aerodynamic information from sampled time histories of the aeroelastic system. More specifically given a finite set of structural modes the method extracts generalized aerodynamic force matrix corresponding to these mode shapes. Once the aerodynamic forces are known, an aeroelastic reduced-order model can be constructed in discrete-time, state-space format by coupling the structural model and the aerodynamic system. The resulting reduced-order model is suitable for constant Mach, varying density analysis.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.20-24
• /
• 2011

In this paper, numerical simulation of cavitation flow for modified NACA66 hydrofoil was made by using the multi-phase RANS equation based on pseudo-compressibility. The Homogeneous mixture model comprised of the mixture continuity, mixture momentum and liquid volume fraction equations was utilized. A vertex-centered finite-volume method was used in conjunction 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing The Spalart-Allmaras one equation model was employed for the closure of turbulence. Reasonable agreements were obtained between the calculation results and the experiment for pressure coefficients on the hydrofoil surface.

• 한국전산유체공학회:학술대회논문집
• /
• 2011.05a
• /
• pp.218-224
• /
• 2011

A supersonic inlet at angle of attack has anti-symmetric pressure distribution, and it can make flow instability and structural problem. In this study, numerical analysis of three-dimensional inviscid flow was conducted under various throttle ratio and angle of attack conditions. Throttle ratio was defined as the ratio of the exit area to the smallest cross section area at inlet, and the ratio is controlled from 0 to 2.42. At various angle of attack, the characteristics of steady and unsteady flow around supersonic inlet is observed under different throttling ratios. From these results, pressure recovery curves and pressure history curves were plotted by post processing. Using pressure history data, FFT analysis is also carried out. Through these processes, it shows the tendency of pressure distribution anti-symmetricity and changing dominant frequency as increasing angle of attack.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.190-201
• /
• 1998

The numerical methodology for computing the impact forces of the water entry bodies has been developed. The present method assumed the impact occurs within a very short time interval and the viscous effects do not have time enough to play a significant role in the impact forces, that is, the flow around a water-entry object was assumed as an inviscid potential flow and is solved by the source panel method. The elements fully submerged into the water are routinely treated, but the elements intersected with the effective planar free surface are redefined and reorganized to be amendable to the source panel method. To validate the present code, it has been applied to disk and ogive model and compared with experimental data. Good agreement has been obtained.

• 한국전산유체공학회:학술대회논문집
• /
• 2007.04a
• /
• pp.177-181
• /
• 2007

Present study examines the numerical issues of cell structure simulation for various regimes of detonation phenomena ranging from weakly unstable to highly unstable detonations. Inviscid fluid dynamics equations with $variable-{\gamma}$ formulation and one-step Arrhenius reaction model are solved by a MUSCL-type TVD scheme and 4th order accurate Runge-Kutta time integration scheme. A series of numerical studies are carried out for the different regimes of the detonation phenomena to investigate the computational requirements for the simulation of the detonation wave cell structure by varying the reaction constants and grid resolutions. The computational results are investigated by comparing the solution of steady ZND structure to draw out the minimum grid resolutions and the size of the computational domain for the capturing cell structures of the different regimes of the detonation phenomena.