• 한국전산유체공학회:학술대회논문집
• /
• 2006.05a
• /
• pp.311-314
• /
• 2006

The supersonic flow around tandem cavities was investigated by three- dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes(RANS) equation with the $\kappa-\omega$ thrbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split using van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge-Kutta method. The aspect ratio of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two-dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the fire cavity flow cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Energy Engineering
• /
• v.21 no.3
• /
• pp.237-242
• /
• 2012

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with a isolator. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, unstart phenomena, fluctuations of specific thrust caused by thermal choking and effects of isolator. The adverse pressure gradient caused by heat addition brings about separation of the wall boundary layers and formation of the oblique shock wave that proceed to upstream. The proceeding speed of the oblique shock wave to upstream direction for the convergent-divergent nozzle with isolator is lower than that for the nozzle without isolator.

• Coupled systems mechanics
• /
• v.6 no.4
• /
• pp.439-464
• /
• 2017

The hydro-elastic system consisting of a pre-stretched highly elastic plate, compressible Newtonian viscous fluid, and the rigid wall is considered and it is assumed that on the plate a lineal-located time-harmonic force acts. It is required to investigate the dynamic behavior of this system and determine how the problem parameters and especially the pre-straining of the plate acts on this behavior. The elasticity relations of the plate are described through the harmonic potential and linearized (with respect to perturbations caused by external time-harmonic force) form of these relations is used in the present investigation. The plane-strain state in the plate is considered and the motion of that is described within the scope of the three-dimensional linearized equations of elastic waves in elastic bodies with initial stresses. The motion of the fluid is described by the linearized Navier-Stokes equations and it is considered the plane-parallel flow of this fluid. The Fourier transform with respect to the space coordinate is applied for a solution to the corresponding boundary-value problem. Numerical results on the frequency response of the interface normal stress and normal velocity and the influence of the initial stretching of the plate on this response are presented and discussed. In particular, it is established that the initial stretching of the plate can decrease significantly the absolute values of the aforementioned quantities.

• Current Optics and Photonics
• /
• v.2 no.3
• /
• pp.261-268
• /
• 2018

Stimulated polariton scattering (SPS) from the $A_1$ transverse optical (TO) modes of $BaTiO_3$ bulk crystal generating a terahertz (THz) wave with the noncollinear phase-matching (NPM) condition is theoretically investigated. To our best knowledge, this is the first report on THz wave generation from $BaTiO_3$ bulk crystal via SPS. Phase-matching (PM) characteristics in the NPM configuration are analyzed. Effective parametric gain lengths for the Stokes and THz waves in the NPM configuration are calculated. The effective parametric gain coefficient and absorption coefficient of the THz wave in $BaTiO_3$ are theoretically simulated. The THz phonon flux densities generated via SPS in $BaTiO_3$ are theoretically calculated by solving the coupled wave equations under the NPM condition. The PM characteristics and THz-wave parametric gain characteristics in $BaTiO_3$ are compared to those in $MgO:LiNbO_3$. The results of the analysis indicate that $BaTiO_3$ is an attractive optical crystal for efficient THz wave generation via SPS.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2004.03a
• /
• pp.671-686
• /
• 2004

The shock-induced combustion ramjet (shcramjet) is a hypersonic airbreathing propulsion concept which over-comes the drawbacks of the long, massive combustors present in the scramjet by using a standing oblique detonation wave (a coupled shock-combustion front) as a means of nearly instantaneous heat addition. A novel shcramjet combustor design that makes use of wedge-shaped flameholders to avoid detonation wave-wall interactions is proposed and analyzed with computational fluid dynamics (CFD) simulations in this study. The laminar, two-dimensional Navier-Stokes equations coupled with a non-equilibrium hydrogen-air combustion model based on chemical kinetics are used to represent the physical system. The equations are solved with the WARP (window-allocatable resolver for propulsion) CFD code (see: Parent, B. and Sislian, J. P., “The Use of Domain Decomposition in Accelerating the Convergence of Quasihyperbolic Systems”, J. of Comp. Physics, Vol. 179, No. 1,2002, pages 140-169). The solver was validated with experimental results found in the literature. A series of steady-state numerical simulations was conducted using WARP and it was deter-mined by means of thrust potential calculations that this combustor design is a viable one for shcramjet propulsion: assuming a shcramjet flight Mach number of twelve at an altitude of 36,000 m, the geometrical dimensions used for the combustor give rise to an operational range for combustor inlet Mach numbers between six and eight. Different shcramjet flight Mach numbers would require different combustor dimensions and hence a variable geometry system in or-der to be viable.

• Journal of Mechanical Science and Technology
• /
• v.20 no.8
• /
• pp.1256-1265
• /
• 2006

The supersonic flows around tandem cavities were investigated by two-dimensional and three-dimensional numerical simulations using the Reynolds-Averaged Navier-Stokes (RANS) equation with the k- ω turbulence model. The flow around a cavity is characterized as unsteady flow because of the formation and dissipation of vortices due to the interaction between the freestream shear layer and cavity internal flow, the generation of shock and expansion waves, and the acoustic effect transmitted from wake flow to upstream. The upwind TVD scheme based on the flux vector split with van Leer's limiter was used as the numerical method. Numerical calculations were performed by the parallel processing with time discretizations carried out by the 4th-order Runge- Kutta method. The aspect ratios of cavities are 3 for the first cavity and 1 for the second cavity. The ratio of cavity interval to depth is 1. The ratio of cavity width to depth is 1 in the case of three dimensional flow. The Mach number and the Reynolds number were 1.5 and $4.5{\times}10^5$, respectively. The characteristics of the dominant frequency between two- dimensional and three-dimensional flows were compared, and the characteristics of the second cavity flow due to the first cavity flow was analyzed. Both two dimensional and three dimensional flow oscillations were in the 'shear layer mode', which is based on the feedback mechanism of Rossiter's formula. However, three dimensional flow was much less turbulent than two dimensional flow, depending on whether it could inflow and outflow laterally. The dominant frequencies of the two dimensional flow and three dimensional flows coincided with Rossiter's 2nd mode frequency. The another dominant frequency of the three dimensional flow corresponded to Rossiter's 1st mode frequency.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.1
• /
• pp.76-84
• /
• 2019

For decades, solitary waves have commonly been used to simulate tsunami conditions in numerical studies. However, the main component of a tsunami waveform acts at completely different spatial and temporal distributions than a solitary waveform. Thus, this study applied a 2-D numerical wave tank that included a non-reflected tsunami generation system based on Navier-Stokes equations (LES-WASS-2D) to directly simulate the run-up of a tsunami-like solitary wave on a slope. First, the waveform and velocity due to the virtual depth factor were applied to the numerical wave tank to generate a tsunami, which made it possible to generate the wide waveform of a tsunami, which was not reproduced with the existing solitary wave approximation theory. Then, to validate the applied numerical model, the validity and effectiveness of the numerical wave tank were verified by comparing the results with the results of a laboratory experiment on a tsunami run-up on a smooth impermeable 1:19.85 slope. Using the numerical results, the run-up characteristics due to a tsunami-like solitary wave on an impermeable slope were also discussed in relation to the volume ratio. The maximum run-up heights increased with the ratio of the tsunami waveform. Therefore, the tsunami run-up is highly likely to be underestimated compared to a real tsunami if the solitary wave of the approximation theory is applied in a tsunami simulation in a coastal region.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.4
• /
• pp.206-216
• /
• 2017

The temporal distribution of the turbulence kinetic energy (TKE) and the vertical component of Reynolds stresses ($-{\bar{u^{\prime}w^{\prime}}}$) was measured during one wave period under high wave energy conditions. The wave data were obtained at Hujeong Beach in the east coast of Korea at January 14~18 of 2017 when an extratropical cyclone was developed in the East Sea. Among the whole thousands of waves measured during the period, hundreds of regular waves that had with similar pattern were selected for the analysis in order to give three representing mean wave patterns using the ensemble average technique. The turbulence properties were then estimated based on the selected wave data. It is interesting to find out that $-{\bar{u^{\prime}w^{\prime}}}$ has one clear peak near the time of flow reversal while TKE has two peaks at the corresponding times of maximum cross-shore velocity magnitudes. The distinguished pattern of Reynolds stress indicates that vertical fluxes of such properties as suspended sediments may be enhanced at the time when the horizontal flow direction is reversed to disturb the flows, supporting the turbulence convection process proposed by Nielsen (1992). The characteristic patterns of turbulence properties are examined using the CADMAS-SURF Reynolds-Averaged Navier-Stokes (RANS) model. Although the model can reasonably simulate the distribution of TKE pattern, it fails to produce the $-{\bar{u^{\prime}w^{\prime}}}$ peak at the time of flow reversal, which indicates that the application of RANS model is limited in the prediction of some turbulence properties such as Reynolds stresses.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.23 no.5
• /
• pp.345-357
• /
• 2011

Recently, as part of diversifying energy sources and earth environmental issues, technology development of new renewable energy using wave energy is actively promoted and commercialized around Europe and Japan etc. In particular, OWC(Oscillating Water Column) wave power generation system using air flow induced by vertical movement of the water surface by waves in an air-chamber within caisson is known as the most efficient wave energy absorption device and therefore, is one of the wave power generation apparatus the closest to commercialization. This study examines air flow velocity, which operates turbine(Wells turbine) directly in oscillating water column type wave power generation structure from two-and three-dimensional numerical experiments and discusses optimal shape of oscillating water column type wave power generation structure by estimating the maximum flow rate of air according to change in shape. The three-dimensional numerical wave flume was applied in interpretation for this study which is the model for the immiscible two-phase flow based on the Navier-Stokes Equation. From this, it turned out that size of optimal shape appears differently according to the incident wave period and air flow is maximized at the period where minimum reflection ratio occurs.

• Ocean Systems Engineering
• /
• v.6 no.1
• /
• pp.23-60
• /
• 2016

The major objective of this study was to develop further understanding of 3D nearshore hydrodynamics under a variety of wave and tidal forcing conditions. The main tool used was a comprehensive 3D numerical model - combining the flow module of Delft3D with the WAVE solver of XBeach - of nearshore hydro- and morphodynamics that can simulate flow, sediment transport, and morphological evolution. Surf-swash zone hydrodynamics were modeled using the 3D Navier-Stokes equations, combined with various turbulence models (${\kappa}-{\varepsilon}$, ${\kappa}-L$, ATM and H-LES). Sediment transport and resulting foreshore profile changes were approximated using different sediment transport relations that consider both bed- and suspended-load transport of non-cohesive sediments. The numerical set-up was tested against field data, with good agreement found. Different numerical experiments under a range of bed characteristics and incident wave and tidal conditions were run to test the model's capability to reproduce 3D flow, wave propagation, sediment transport and morphodynamics in the nearshore at the field scale. The results were interpreted according to existing understanding of surf and swash zone processes. Our numerical experiments confirm that the angle between the crest line of the approaching wave and the shoreline defines the direction and strength of the longshore current, while the longshore current velocity varies across the nearshore zone. The model simulates the undertow, hydraulic cell and rip-current patterns generated by radiation stresses and longshore variability in wave heights. Numerical results show that a non-uniform seabed is crucial for generation of rip currents in the nearshore (when bed slope is uniform, rips are not generated). Increasing the wave height increases the peaks of eddy viscosity and TKE (turbulent kinetic energy), while increasing the tidal amplitude reduces these peaks. Wave and tide interaction has most striking effects on the foreshore profile with the formation of the intertidal bar. High values of eddy viscosity, TKE and wave set-up are spread offshore for coarser grain sizes. Beach profile steepness modifies the nearshore circulation pattern, significantly enhancing the vertical component of the flow. The local recirculation within the longshore current in the inshore region causes a transient offshore shift and strengthening of the longshore current. Overall, the analysis shows that, with reasonable hypotheses, it is possible to simulate the nearshore hydrodynamics subjected to oceanic forcing, consistent with existing understanding of this area. Part II of this work presents 3D nearshore morphodynamics induced by the tides and waves.