• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 1998년도 춘계 학술대회논문집
• /
• pp.99-109
• /
• 1998

This paper describes the analysis of flow and aeroacoustic field around a car-like body. The governing equations, 3-D unsteady incompressible Navier-Stokes equations, are solved with the iterative time marching scheme. The Chimera grid technique has been applied to efficiently simulate the flow around the side-view mirror, After the flow field analysis has been converged, the aerodynamic noise analysis of the side-view mirror has been performed by solving Ffowcs Williams and Hawkings equation. From the present numerical simulation, the A- and C-pillar vortex are evidently shown and the aerodynamic noise level induced by the side-view mirror is predicted to about 100dB.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 1997년도 제8회 학술강연회논문집
• /
• pp.43-51
• /
• 1997

Steady and unsteady numerical simulations are conducted for the experiments performed to investigate the ram accelerator flow field by using the expansion tube facility in Stanford University. Navier-Stokes equations for chemically reacting flows are analyzed by fully implicit and time accurate numerical methods with Jachimowski's detailed chemistry model for hydrogen-air combustion involving 9 species and 19 reaction steps. Although the steady state assumption shows a good agreement with the experimental schlieren and OH PLIF images for the case of $2H_2$+$O_2$+$17N_2$, it fails in reproducing the combustion region behind the shock intersection point shown in the case of $2H_2$+$O_2$+$12N_2$, mixture. Therefore, an unsteady numerical simulation is conducted for this case and the result shows all the detailed flow stabilization process. The experimental result is revealed to be an instantaneous result during the flow stabilization process. The combustion behind the shock intersection point is the result of a normal detonation formed by the intersection of strong oblique shocks that exist at early stage of the stabilization process. At final stage, the combustion region behind the shock intersection point disappears and the steady state result is retained. The time required for stabilization of the reacting flow in the model ram accelerator is found to be very long in comparison with the experimental test time.

• 한국전산유체공학회지
• /
• 제4권1호
• /
• pp.19-26
• /
• 1999

This paper describes a numerical method for predicting the incompressible unsteady laminar three-dimensional flows with free-surface. The Navier-Stokes equations governing the flows have been discretized by means of finite-difference approximations, and the resulting equations have been solved via the SIMPLE-C algorithm. The free-surface is defined by the motion of a set of marker particles and the interface behaviour was investigated by means of a "Lagrangian" technique. Using the GALA concept of Spalding, the conventional mass continuity equation is modified to form a volumetric or bulk-continuity equation. The use of this bulk-continuity relation allows the hydrodynamic variables to be computed over the entire flow domain including both liquid and gas regions. Thus, the free-surface boundary conditions are imposed implicitly and the problem formulation is greatly simplified. The numerical procedure is validated by comparing the predicted results of a periodic standing waves problems with analytic solutions. The results show that this numerical method produces accurate and physically realistic predictions of three-dimensional free-surface flows.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년도 제30회 춘계학술대회논문집
• /
• pp.275-280
• /
• 2008

Shock tube flow measurement has been often troubled with a finite opening time of diaphragm, but there is no systematic work to investigate its effect on the shock tube flow. In the present study, both the experimental and computational works have been performed on the shock tube flows at low pressure ratios. The computational analysis has been performed using the two-dimensional, unsteady, compressible Navier-Stokes equations, based upon a TVD MUSCL finite difference scheme. It is known that the present computational results reproduce the experimental data with good accuracy and simulate successfully the process of diaphragm opening as a function of time. The concept of an imaginary center is introduced to specify the non-centered expansion wave due to a finite opening time of diaphragm. The results obtained show that the diaphragm opening time is reduced as the initial pressure ratio of shock tube increases, leading to the effect of a finite opening time of diaphragm to be more remarkable at low pressure ratios.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2009년도 춘계학술대회 논문집
• /
• pp.315-319
• /
• 2009

Investigation of the lateral force fluctuations in an axisymmetric overexpanded compressed truncated perfect (CTP) nozzle for the shutdown transient is presented. These nozzles experience side-loads during start-up and shut-down operations, because of the flow separation at nozzle walls. Two types of flow separations such as free shock separation (FSS) and restricted shock separation (RSS) shock structure occur. A two-dimensional unsteady numerical simulation has been carried out over an axisymmetric CTP nozzle to simulate the lateral force fluctuations in nozzle during shutdown process. Reynolds Averaged Navier-Stokes equations are numerically solved using a fully implicit finite volume scheme. Governing equations are solved by coupled implicit scheme. Two equation k-$\omega$ SST turbulence model is selected. Unsteady pressure is measured at four locations along the nozzle wall. Present pressure variation compared well with the experimental data. During shutdown transient, separation pattern varies from FSS to RSS and finally returns to FSS. Several pressure peaks are observed during the RSS separation pattern. These pressure peaks generate lateral force or side loads in rocket nozzle.

• Wind and Structures
• /
• 제16권5호
• /
• pp.411-431
• /
• 2013

A computational study of vortex-induced transverse vibrations of a cylinder with low mass-damping is presented. An Arbitrary Lagrangian-Eulerian (ALE) formulation of the Unsteady Reynolds-Averaged Navier-Stokes equations (URANS), along with the Spalart-Allmaras (SA) one-equation turbulence model, are coupled conservatively with rigid body motion equations of the cylinder mounted on elastic supports in order to study the amplitude and frequency response of a freely vibrating cylinder, its flow-induced motion, Vortex Street, near-wake flow structure, and unsteady loading in a moderate range of Reynolds numbers. The time accurate response of the cylinder from rest to its limit cycle is studied to explore the effects of Reynolds number on the start of large displacements, motion amplitude, and frequency. The computational results are compared with published physical experiments and numerical studies. The maximum amplitudes of displacements computed for various Reynolds numbers are smaller than the experimental values; however, the overall agreement of the results is quite satisfactory, and the upper branch of the limit-cycle displacement amplitude vs. reduced velocity response is captured, a feature that was missed by other studies. Vortex shedding modes, lock-in phenomena, frequency response, and phase angles are also in agreement with experiments.

• Wind and Structures
• /
• 제39권1호
• /
• pp.31-45
• /
• 2024

This paper examines the flow characteristics around an inclined prism with a U-shaped cross-section ("U-profile") and investigates the connection between the flow and flow-induced vibrations. The study employs a combined approach that involves wind tunnel experiments and computational fluid dynamics (CFD) using an unsteady Reynolds-averaged Navier-Stokes (RANS) turbulence model. Distinct vortex formation patterns are observed in the flow field surrounding the stationary inclined profile. When the cavity of the profile faces away from the incoming flow, large vortices develop behind the profile. Conversely, when the cavity is oriented towards the oncoming flow, these vortices form within the cavity. Notably, due to the slow movement of these large vortices through the cavity, the frequency at which vortices are shed in the negative inclination case is lower compared to the positive inclination, where they form in the wake. Wind tunnel experiments reveal an intermittent transition between the two vortex formation patterns at zero inclination. Large vortices sporadically emerge both in the cavity and behind the profile. The simulation results demonstrate that when these large vortices occur at a frequency close to the structure's natural frequency, they induce prominent pitch vibrations. This phenomenon is also sought after and presented in coupled vibration experiments. Additionally, the simulations indicate that when the natural frequency of the structure is considerably lower than the vortex shedding frequency, this type of vibration can be observed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.119-125
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory and solitary wave theory are employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with each other. The marker-density method is better then the hybrid method. Also they are compared with other existing experimental results. The Agreement between the experimental data and the computation results is good.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2002년도 학술대회지
• /
• 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.

• Journal of Ship and Ocean Technology
• /
• 제8권2호
• /
• pp.61-69
• /
• 2004

Unsteady nonlinear wave motions on the free surface over a plane beach of constant slope are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier-Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. The second-order Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with the plane beach of constant slope in surf zone. The results are compared with other existing experimental results. Agreement between the experimental data and the computation results is good.