• 제목/요약/키워드: Navier Stokes Eq.

검색결과 10건 처리시간 0.021초

LES와 Dynamic Smagorinsky 난류모형을 이용한 쇄파역에서의 경계층 Streaming 수치해석 (Numerical Analysis of the Hydraulic Characteristics of a Boundary Layer Streaming over Surf-Zone Using LES and Dynamic Smagorinsky Turbulence Model)

  • 조용준
    • 한국해안·해양공학회논문집
    • /
    • 제32권1호
    • /
    • pp.69-84
    • /
    • 2020
  • 자연 해빈은 해양환경에 따라 침·퇴적을 반복하며 고 파랑에 의해 해빈이 대규모로 침식되더라도 폭풍이 잦아들고 다시 너울이 우세한 해양환경이 회복되는 경우 점진적으로 복원되며, 이러한 해빈 복원은 경계층 streaming을 통해 이루어진다. 이처럼 경계층 streaming은 그 공학적 가치에도 불구하고 해안에서 가용한 표사의 대부분이 공급되는 쇄파 역에서의 경계층 streaming에 대한 우리의 이해는 아직 상당히 부족하다. 이러한 인식에 기초하여 본 연구에서는 쇄파역 경계층 streaming 수리특성을 살펴보기 위해 단조 해안과 사주를 포함한 해안에서의 천수 과정을 수치모의하였다. 수치 모의는 Spatially filtered Navier-Stokes Eq., LES(Large Eddy Simulation), Dynamic Smagorinsky 난류모형으로 구성된 정교한 수치모형에 기초하여 수행되었으며, 이 과정에서 k-ε 난류모형과 LES Turbulence Closure가 모의결과에 미치는 영향도 함께 살펴보았다. 모의결과 해안공학계에 잘 알려진 k-ε 난류모형의 한계로 인해 wall function에 기반한 k-ε 난류모형의 경우 LES와 비교하면 저면 인근 유속이 다소 과다하게 모의 되었다. 또한, 바닥과 가까운 해역에서의 유속이 바닥의 영향으로부터 비교적 자유로운 상층부에서의 유속보다 우월한 Longuet-Higgins(1957)가 이야기하는 전형적인 경계층 streaming이 천수 초입부에서부터 쇄파 역 깊숙이까지 존재하는 것을 확인하였다. 또한, 주기가 상대적으로 긴 경우 경계층 streaming의 세기와 생성범위는 해안 방향으로 확대되며 이러한 경향은 경계층 streaming이 바닥 인근에서 진행되는 마찰로 인한 파랑에너지손실로 결과되며 주기가 긴 경우 천수 과정이 일찍 시작된다는 사실을 상기하면 충분히 수용 가능해 보이며, Longuet-Higgins(1957)의 해석 해에서도 같은 경향을 확인할 수 있다.

불규칙 파랑 비선형 천수 과정 수치해석 - 천수 단계별 파고분포 변화를 중심으로 (Numerical Analysis of Nonlinear Shoaling Process of Random Waves - Centered on the Evolution of Wave Height Distribution at the Varying Stages of Shoaling Process)

  • 김용희;조용준
    • 한국해안·해양공학회논문집
    • /
    • 제32권2호
    • /
    • pp.106-121
    • /
    • 2020
  • 항 외곽시설 신뢰성 설계가 합리적으로 구현하기 위해서는 우리나라 해양환경 특성이 반영된 확률모형이 필요하며 이러한 시각에서 본 연구에서는 천 해역 확률모형 개발을 위한 기초연구의 일부로 불규칙 파랑 천수 과정을 수치 모의하였다. 수치 모의는 자연해안에서 흔히 관측되는 사주가 원빈에 형성된 해안을 대상으로 수행하였으며 파랑모형은 spatially filtered Navier-Stokes Eq., LES[Large Eddy Simulation], one equation dynamic Smagorinsky turbulence closure 등으로 구성하였다. 불규칙 파랑은 우리나라 동해안에서 관측되는 너울 특성을 반영하기 위해 다양한 첨두 증강계수를 지니는 JONSWAP 스펙트럼과 random phase method를 사용하여 모의하였다. 파고분포의 모수는 먼저 수치 모의에서 관측된 자유수면 시계열 자료를 threshold crossing method로 파별 해석[wave by wave analysis]하여 개별 파랑을 특정하고, 이어 이렇게 특정된 파마루와 파곡 빈도 해석결과로부터 산출하였다. 모의결과 현재 천 해역 파고분포를 대표하는 수정 Glukhovskiy 파고분포는 큰 파고와 작은 파고 발생확률은 과다하게, 중간 크기 파고 발생확률은 과소하게 평가하는 것으로 모의 되었으며, 이에 반해 본 논문에서 제시된 파고분포의 경우 일치도가 상당하였다. 또한, 전술한 수정 Glukhovskiy 파고분포와의 간극은 쇄파역에서 제일 현저하게 관측되어 수정 Glukhovskiy 파고분포를 쇄파역 언저리에 거치되는 외곽시설 신뢰성 설계에 적용하는 일은 지양되어야 할 것으로 판단된다.

고속열차의 선두부 형상에 따른 공력특성 변화 (Aerodynamic Characteristics for various front shapes of High Speed Train)

  • 이승철;김사랑;허남건
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 1995년도 추계 학술대회논문집
    • /
    • pp.49-54
    • /
    • 1995
  • A numerical analysis on the effect of the front shape on the aerodynamic characteristics of HST model is made, using FVM based general purpose 3D Navier-Stokes eq. solver, TURBO-3D program. Numerical solutions are compared with each case of different front shape for HST model. The result shows a good quantitative aerodynamic characteristic tendencies for variation of front shape of HST. Thus it may be used as a basis in the design of the shape of real HST.

  • PDF

장파 제어체 Eco-breaker 2의 수리특성 (On the Hydraulic Characteristics of Efficient Long Wave Energy Absorber-Eco-breaker 2)

  • 조용준;김호민
    • 대한토목학회논문집
    • /
    • 제28권5B호
    • /
    • pp.547-558
    • /
    • 2008
  • 최근 폭증하는 해상물동량으로 인해 선박은 날로 대형화되며 현재 우리나라에서는 이러한 추이에 대응하기 위해 초대형 항만건설이 서둘러 진행되고 있으며 그 중 일부는 이미 부분적으로 완성되어 운영 중이거나 계속적인 확장이 예정되어있다. 초대형 항만의 경우 항 정온도 확보를 위해 대형 직립형 방파제가 수반되기 마련이나 비교적 높은 반사율로 인해 장주기의 파랑에 취약한 구조적 문제를 지닌다. 장주기 파랑의 효과적 제어를 위해 본 고에서는 아치형 curtain wall과 유수실로 구성된 된 장파 제어체 -Eco-breaker 2- 가 제시되며 Eco-breaker 2의 수리특성을 수치적으로 규명하였다. 수치모형은 가장 강건한 파랑모형인 Navier-Stokes Eq.과 SPH(Smoothed Particle Hydrodynamics), LES(Large Eddy Simulation) 기법에 기초하여 구성하였다. 수치모형의 검증을 위한 수리모형실험도 병행하였다. 수치모의 결과 Eco-breaker 2는 유수실 내부에 생성되어 아치형 curtain wall 전면 수역과 유수실을 주기적으로 넘나드는 대규모 구조적 와로 장파를 효과적으로 제어할 수 있었다. 또한 대규모 구조적 와의 생성 영역과 세기는 내습하는 파랑의 주기에 종속하며 유수실 체류시간이 입사파 주기의 반을 하회하는 경우 반사계수는 최대 0.18까지 낮아진다. 이러한 결과는 Eco-breaker 2가 환경친화적 항의 개발에 적극적으로 활용될 수 있음을 시사하는 것으로 이를 위해 그 구체적 제원도 함께 제시하였다.

SPH을 활용한 3차원 비선형 파랑모형 개발 (Development of 3-D Nonlinear Wave Driver Using SPH)

  • 조용준;김권수
    • 대한토목학회논문집
    • /
    • 제28권5B호
    • /
    • pp.559-573
    • /
    • 2008
  • Navier-Stokes식, Gaussian 분포형 용출함수를 이용한 내부조파, energy absorbing layer로 삼차원 파랑모형을 새롭게 구성하였다. Navier-Stokes식의 수치적분에는 정교한 수치기법인 SPH(Smoothed Particle Hydrodynamics)가 활용된다. 제안된 파랑모형의 검증은 삼차원 포물형 용기에서의 sloshing현상과 Thacker(1981)의 해석해를 토대로 수행되었다. 초기 수면 형상이 Gaussian hump인 경우와 일방향으로 경사진 경우에 대해 수치모의 하였다. 수치모의 결과 수면이 융기되도록 구속한 외부조건이 해제되면서 시작되는 자유진동의 정성적 거동은 비교적 정확히 모의되었으나 시간이 경과될수록 위상차, 침수선이 퇴각하는 등 초기 수면과는 상당히 다른 결과를 보였다. 최종적인 검증은 쐐기모양 해안에서의 비선형 천수, 굴절거동의 수치모의를 토대로 진행되었다. 수치모의 결과 굴절되는 양이 Hamiltonian ray theory가 제공하는 수치보다 전반적으로 작게 나타났다. 이러한 현상은 이상유체와 선형 이론에 기초한 Hamiltonian ray theory에서 간과된 비선형성, 점성으로 인한 양안과 저면에서의 에너지 감쇄, 쇄파 과정에 유동계에 도입되는 에너지 감쇄, 선행파랑에 의한 down-rush와 조우시 발생하는 도수 등에 기인하는 것으로 판단된다.

Oil boom과 파랑의 비선형상호작용을 고려한 Oil Boom의 누유특성 (Failure Characteristics of Oil Boom Considering the Nonlinear Interaction of Oil Boom with Waves)

  • 조용준;윤대경
    • 한국해안·해양공학회논문집
    • /
    • 제23권3호
    • /
    • pp.193-204
    • /
    • 2011
  • 역동적인 파랑에 노출되는 경우 다양한 failure mode를 쉽게 드러내는 Oil boom의 성능을 개선하기 위해 가장 정교한 파랑모형인 spatially filtered Navier-Stokes 식을 LES (Large Eddy Simulation), 잔차응력에 대한 LDS (Lagrangian Dynamic Smagorinsky 모형), SPH (Smoothed Particle Hydrodynamics) 기법을 활용하여 해석하는 새로운 수치모형이 제언되었다. 이어 부유식 Oil Boom의 누유특성을 규명하기 위해 oil spill, progressive wave, oil boom의 상호작용을 oil boom이 계류삭에 고정되어있는 경우와 oil boom의 excursion이 허용된 경우에 대해 각각 수치모의 하였다. 모의결과 oil boom의 skirt 길이가 수심의 30% 이상이고 excursion이 허용된 경우 oil spill의 차폐 기능은 극대화되는 것으로 밝혀졌다. 이와 더불어 y = 1~2 m 사이에 오일막과 해수의 경계층에서 생성된 와류가 저면으로 확산되면서 시계방향과 반 시계방향의 와류가 엇갈리게 생성되는 coherent eddies가 관측되어 수리실험을 통해 그 존재가 알려진 Kelvin-Helmholz파의 성장과정과 계면으로부터의 일탈과정이 수치모의된 것으로 판단된다.

협수로에서 생성되는 고립파 형태의 항주파와 항주파류 (Solitary Wave-like Ship Induced Waves and Its Associated Currents in a Water Channel of Narrow Width)

  • 조용준;최한림
    • 한국해안·해양공학회논문집
    • /
    • 제27권4호
    • /
    • pp.202-216
    • /
    • 2015
  • 가용한 토지의 부족으로 현재 우리나라 인공운하에서 흔히 관측되는 협수로에서는 고립파 형태의 항주파가 생성되는 것으로 추정된다. 본고에서는 이러한 가설을 확인하기 위해 수치모의를 수행하였다. 수치모의는 삼차원 Navier Stokes 식과 VOF에 기초하여 수행되었으며, 수치모형의 검증은 현재 우리에게 가용한 운하 설계기준 중 가장 빈번히 언급되는 PIANC (1987) 설계안과 본고에서 유도된 해석해에 기초하여 수행되었다. 모의결과, 고립파형태의 항주파가 관측되었으며, 하안 인근에서 계측된 선수파 파고의 경우 수치모의 결과는 PIANC (1987) 설계안을 상회하였으며 상당히 오랜 기간 (이십초 내외) 지속되었다. 선미파의 경우는 하회하였다. 이러한 결과는 광역수로에서 관측되는 항주파의 일반적인 특성과는 상이한 것으로 좁은 수로 폭이 항주파 특성에 영향을 미치는 것으로 판단된다. 또한 선박의 협수로 통과시 선박을 중심으로 양안으로 진행되는 흐름과, 반사로 인해 양안에서 선박으로 진행되는 흐름이 교대로 출현하였으며 이 때 최대 0.90 m/s의 유속이 전 수심대역에서 비교적 균일하게 유지되어, 상당한 쇄굴이 진행될 수 있다고 판단된다.

파형의 왜도가 투과성 방파제 투과율에 미치는 영향 (The Effect of Skewness of Nonlinear Waves on the Transmission Rate through a Porous Wave Breaker)

  • 조용준;강윤구
    • 한국해안·해양공학회논문집
    • /
    • 제29권6호
    • /
    • pp.369-381
    • /
    • 2017
  • 쇄파역에서 흔히 관측되는 왜도된 비선형 파랑은 정현파에 비해 상대적으로 큰 관성력을 지녀 투과율에 미치는 영향이 상당할 것으로 추정되어왔다. 본고에서는 이러한 추론을 확인하기 위해 왜도가 1.02에서 1.032에 분포하는 총 여섯 개의 파랑을 대상으로 한 수치모의를 수행하였다. 수치모의는 OpenFoam 기반 Tool box인 ihFoam을 사용하여 수행하였으며, 왜도된 비선형 파랑은 경사가 1:30인 천퇴를 이용하여 유도하였다. 투수층에서의 유동은 Navier Stokes Eq.에 추가 항력을 도입하여 해석하였으며, 모의결과 왜도된 비선형 파랑은 정현파에 비해 큰 관성력으로 인해 상대적으로 큰 투과율을 지니는 것으로 모의되었다. 또한 Porous층에서의 Damping으로 인한 파속변화 여부를 살펴 보기위해 Damping이 반영된 천수방정식으로부터 분산관계식을 새로이 유도하였으며, 그 결과 투과성 방파제에서의 파형과 마찰력 사이의 위상차가 파속에 영향을 미치는 것으로 보인다.

3차원 공동의 폭변화에 따른 초음속 유동에 대한 수치분석연구 (NUMERICAL ANALYSIS OF THREE DIMENSIONAL SUPERSONIC CAVITY FLOW FOR THE VARIATION OF CAVITY SPANWISE RATIO)

  • 우철훈;김재수;최홍일
    • 한국전산유체공학회:학술대회논문집
    • /
    • 한국전산유체공학회 2006년도 추계 학술대회논문집
    • /
    • pp.181-184
    • /
    • 2006
  • High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation and reattachment, shock and expansion waves. The general cavity flow phenomena include the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity' flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions, The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio(L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyized and compared with the results of Rossiter's Eq.

  • PDF

공동의 폭 변화에 따른 3차원 초음속 공동 유동연구 (NUMERICAL ANALYSIS OF THREE DIMENSIONAL SUPERSONIC CAVITY FLOW FOR THE VARIATION OF CAVITY SPANWISE RATIO)

  • 우철훈;김재수
    • 한국전산유체공학회지
    • /
    • 제11권4호
    • /
    • pp.62-66
    • /
    • 2006
  • High-speed flight vehicle have various cavities. The supersonic cavity flow is complicated due to vortices, flow separation, reattachment, shock waves and expansion waves. The general cavity flow phenomena includes the formation and dissipation of vortices, which induce oscillation and noise. The oscillation and noise greatly affect flow control, chemical reaction, and heat transfer processes. The supersonic cavity flow with high Reynolds number is characterized by the pressure oscillation due to turbulent shear layer, cavity geometry, and resonance phenomenon based on external flow conditions. The resonance phenomena can damage the structures around the cavity and negatively affect aerodynamic performance and stability. In the present study, we performed numerical analysis of cavities by applying the unsteady, compressible three dimensional Reynolds-Averaged Navier-Stokes(RANS) equations with the ${\kappa}-{\omega}$ turbulence model. The cavity model used for numerical calculation had a depth(D) of 15mm cavity aspect ratio (L/D) of 3, width to spanwise ratio(W/D) of 1.0 to 5.0. Based on the PSD(Power Spectral Density) and CSD(Cross Spectral Density) analysis of the pressure variation, the dominant frequency was analyzed and compared with the results of Rossiter's Eq.