• 에너지공학
• /
• 제21권3호
• /
• pp.237-242
• /
• 2012

이 연구에서는 초음속 연소에서 발생하는 열폐색의 다양한 현상을 규명하기 위해 확대 축소 노즐 내부에서 열폐색에 의해 형성되는 2차원의 초음속의 비정상 유동장에 대한 수치해석 결과를 제시한다. 열폐색에 의해 야기되는 이동 충격파를 수치계산하기 위해 TVD 스킴을 이용하며, 노즐의 확대부의 일정영역에 가열을 통하여 열폐색이 발생할 수 있는 조건을 형성하고, 이 때 발생하는 천이현상을 분리부가 있는 경우와 없는 경우에 대해 불시동현상 발생속도, 비추력의 불안정성 등을 통해 비교, 설명한다. 분리부가 있는 경우가 없는 경우에 비해 열폐색에 의해 발생한 경사 충격파가 느린 속도로 상류측으로 이동하여 분리부의 설치가 엔진 불시동의 지연효과가 있음을 제시하였다.

• Coupled systems mechanics
• /
• 제6권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
• /
• 제2권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• 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
• /
• 제20권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.

• 한국해양공학회지
• /
• 제33권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.

• 한국해안·해양공학회논문집
• /
• 제29권4호
• /
• pp.206-216
• /
• 2017

교란운동에너지(TKE)와 레이놀즈 응력의 수직성분($-{\bar{u^{\prime}w^{\prime}}}$)에 대한 한 주기 파장 안에서의 시간변화를 관측자료를 사용하여 분석하였다. 관측자료는 동해에서 온대성저기압이 발달하였던 2017년 1월 14일부터 18일까지 동해안 후정해변에서 측정한 파랑자료를 사용하였다. 이 기간 동안 관측된 모든 파랑자료들 중에서 비슷한 형태를 갖는 수백 개의 규칙파들을 구분하였으며 이 자료를 토대로 Ensemble Average 기법을 사용하여 이 기간 파랑특성을 대표하는 세 개의 평균파를 계산하였다. 그리고 이 평균파를 기준으로 각 파의 요동을 측정하여 한 주기 동안의 교란운동에너지와 레이놀즈 응력을 계산하였다. 이렇게 계산된 자료들을 분석한 결과 교란운동에너지는 파랑의 평균유속과 비슷한 분포를 나타내었으나(즉 유속이 최대값을 나타낼 때 교란운동에너지도 최대값을 나타내었다), $-{\bar{u^{\prime}w^{\prime}}}$는 파랑의 수평유속 방향이 전환되는 '방향전환점'에서 가파르게 증가하는 경향을 나타내었다. 이러한 $-{\bar{u^{\prime}w^{\prime}}}$의 독특한 분포는 Nielsen(1992)에 의해 제안된 난류 convection 현상을 뒷받침하는 발견으로 퇴적물과 같은 물질들의 부유현상이 파랑의 '방향전환점(한 주기 안에서 파랑의 횡단방향 유속 부호가 바뀌는 시점)'에서 촉진될 수 있음을 보여준다. 이렇게 관측된 난류에너지 분포 특성을 CADMAS-SURF 모델을 사용하여 구현해 보았다. 그 결과 교란운동에너지의 경우 모델결과와 관측치 사이에 유사성이 발견되었으나 레이놀즈 응력($-{\bar{u^{\prime}w^{\prime}}}$)의 경우 모델이 '방향전환점'에서의 증가현상을 구현해 내지 못하였다. 이는 CADMAS-SURF와 같은 Reynolds-Averaged Navier-Stokes(RANS) 모델들이 가지는 한계점으로 RANS 모델의 경우 레이놀즈 응력과 같은 난류에너지가 평균유속의 분포에 강한 영향을 받기 때문인 것으로 판명되었다.

• 한국해안·해양공학회논문집
• /
• 제23권5호
• /
• pp.345-357
• /
• 2011

최근 지구환경문제와 에너지원의 다각화를 위한 일환으로 파랑에너지를 이용하는 신재생에너지의 기술개발이 유럽과 일본 등을 중심으로 활발히 추진 및 실용화되고 있다. 특히, 케이슨 내의 공기실에서 파랑에 의한 수면의 상하운동으로 유도되는 공기흐름을 이용하는 진동수주형 파력발전시스템은 가장 효율적인 파랑에너지흡수장치로 알려져 있고, 따라서 상업화에 가장 근접한 파력발전장치 중에 하나이다. 본 연구에서는 진동수주형 파력발전구조물에서 터빈(Wells터빈)에 직접 작용하는 공기흐름속도를 2차원 및 3차원수치실험으로부터 검토하며, 이 때 형상의 변화에 따른 공기의 최대흐름속도를 추정하여 진동수주형 파력발전구조물의 최적형상을 논의한다. 수치해석에서는 기체와 액체의 혼상동적현상을 동일한 지배방정식으로 해석하는 혼상류(2상류)수치모델에 기초한 3차원수치파동수로를 적용하였다. 이로부터 입사주기대에 따라 최적형상의 크기가 상이하게 나타나는 것을 확인할 수 있었고, 최소의 반사율이 발생하는 주기 대에서 공기흐름이 최대로 된다는 것을 알 수 있었다.

• Ocean Systems Engineering
• /
• 제6권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.

• 대한토목학회논문집
• /
• 제28권1B호
• /
• pp.95-109
• /
• 2008

본 연구에서는 LDS 난류응력 모형, Van Rijn의 pick up 함수를 활용하여 일정 경사부에서의 파랑의 이행과 천수, 연이은 쇄파현상, plunging breaker에 후행하는 해저질의 역동적인 부유와 down rush와 후행 파랑에 의한 표사의 재분배를 수치모의 하였다. 이 과정에서 해저질과 소통하는 저면 유체력에 대한 quadratic law를 중심으로 한 기존의 연구 성과들은 정상상태에 기초하여 급속히 가속되고 감속되는 swash 대역의 수리특성을 반영할 수 없다는 결론에 도달하고 이러한 인식에 기초하여 새로운 산출방법이 제시되었다. 새로운 산출방법을 토대로 수치모의하여 비선형 천수과정의 일반적인 특징, 동조 비동조 고차 조화성분으로 전이된 파랑에너지로 인해 상당히 예리하고 왜도된 파형, 파형의 마루로부터 시작되는 물입자 자유낙하, 착수로 인한 커다란 물보라의 형성, 물보라 형성층의 해변으로의 이행, wave finger (Narayanaswamy와 Dalrymple, 2002), swash 대역에서 진행되는 부유사 순환과정, swash 대역에서 처오름으로 인해 부유된 부유사 무리의 off shore 방향으로의 순 이동 등이 비교적 정확히 재현되는 등 상당히 고무적인 결과를 얻을 수 있었다. 이러한 결과는 기존의 Euler 좌표계에서 정의되는 파랑모형과 이동경계 기법의 한계를 뛰어 넘는 것으로 향후 보다 정확한 침식해석이 가능 할 것으로 판단된다.