• 대한기계학회논문집A
• /
• 제24권2호
• /
• pp.472-480
• /
• 2000

A macroscopic yield criterion for porous solids with pressure-sensitive matrices modeled by Coulomb's yield criterion was obtained by generalizing Gurson's yield criterion with consideration of the hydrostatic yield stresses for a spherical thick-walled shell and by fitting the finite element results of a voided cube. The macroscopic yield criterion is valid for negative mean normal stresses as well as for positive mean normal stresses. From the yield criterion, a plastic potential function for the porous solids was derived either for plastic normality flow or for plastic non-normality flow of pressure- sensitive matrices. In addition, the elastic relation, an evolution equation of the plastic flow stress of the matrices and an evolution equation of the void volume fraction were presented to complete a set of constitutive relations. The set of constitutive relations was implemented into a finite element code ABAQUS to analyze the material behavior of rubber-toughened epoxies. The cavitation and the deformation behavior were analyzed around a crack tip under three-point bending and around notch tips under four-point bending. In the numerical analyses, the cavitation of rubber particles was considered via a stress-controlled nucleation model. The numerical results indicate that a reasonable cavitation zone can be obtained with void nucleation controlled by the macroscopic mean normal stress, and a plastic zone is smaller around a notch tip under compression than under tension. These numerical results agree well with corresponding experimental results on the cavitation and plastic zones.

• 한국수자원학회논문집
• /
• 제44권6호
• /
• pp.429-438
• /
• 2011

In order to simulate a free surface flow in a trench channel, a three-dimensional incompressible unsteady Reynolds-averaged Navier-Stokes (RANS) equations are closed with the ${\kappa}-{\epsilon}$ model. The artificial compressibility (AC) method is used. Because the pressure fields can be coupled directly with the velocity fields, the incompressible Navier-Stokes (INS) equations can be solved for the unknown variables such as velocity components and pressure. The governing equations are discretized in a conservation form using a second order accurate finite volume method on non-staggered grids. In order to prevent the oscillatory behavior of computed solutions known as odd-even decoupling, an artificial dissipation using the flux-difference splitting upwind scheme is applied. To enhance the efficiency and robustness of the numerical algorithm, the implicit method of the Beam and Warming method is employed. The treatment of the free surface, so-called interface-tracking method, is proposed using the free surface evolution equation and the kinematic free surface boundary conditions at the free surface instead of the dynamic free surface boundary condition. AC method in this paper can be applied only to the hydrodynamic pressure using the decomposition into hydrostatic pressure and hydrodynamic pressure components. In this study, the boundary-fitted grids are used and advanced each time the free surface moved. The accuracy of our RANS solver is compared with the laboratory experimental and numerical data for a fully turbulent shallow-water trench flow. The algorithm yields practically identical velocity profiles that are in good overall agreement with the laboratory experimental measurement for the turbulent flow.

• 한국대기환경학회지
• /
• 제25권4호
• /
• pp.304-315
• /
• 2009

Numerical simulations were carried out to understand the effect of Sea Surface Temperature (SST) spatial distribution on regional circulation. A three-dimensional non-hydrostatic atmospheric model RAMS, version 6.0, was applied to examine the impact of SST forcing on regional circulation. New Generation Sea Surface Temperature (NGSST) data were implemented to RAMS to compare the results of modeling with default SST data. Several numerical experiments have been undertaken to evaluate the effect of SST for initialization. First was the case with NGSST data (Case NG), second was the case with RAMS monthly data (Case RM) and third was the case with seasonally averaged RAMS monthly data (Case RS). Case NG showed accurate spatial distributions of SST but, the results of RM and RS were $3{\sim}4^{\circ}C$ lower than buoy observation data. By analyzing practical sea surface conditions, large difference in horizontal temperature and wind field for each run were revealed. Case RM and Case RS showed similar horizontal and vertical distributions of temperature and wind field but, Case NG estimated the intensity of sea breeze weakly and land breeze strongly. These differences were due to the difference of the temperature gradient caused by different spatial distributions of SST. Diurnal variations of temperature and wind speed for Case NG indicated great agreement with the observation data and statistics such as root mean squared error, index of agreement, regression were also better than Case RM and Case RS.

• 천문학논총
• /
• 제30권2호
• /
• pp.699-701
• /
• 2015

The high plateaus in west China may provide suitable sites for astronomical observations with the institute's middle-range telescopes and possibly with larger telescopes. Under China-Japan collaborations for site survey in west China, we have been conducting searches for good sites and monitoring their characteristics over several years. As recent results of our site survey show, sites in west Tibet are revealed with a high possibility of good astronomical observations. Weather characteristics at Gar in Ali, Tibet, show its high clear-sky ratios, especially in winter, comparable to Mauna Kea, Hawaii. But it has some wind problem in winter, where stronger wind speeds, over 20m/sec, occur frequently even though the sky is clear. To find calmer sites, we have conducted numerical simulations for the Ali area using the Japan Meteorological Agency NonHydrostatic Model. We have found another site, named ZoZo Hill, near Gar. We will continue to monitor the Gar site to clarify weather characteristics over the whole year and hopefully start to negotiate for site monitoring at ZoZo Hill this year.

• 한국해안해양공학회지
• /
• 제8권2호
• /
• pp.123-136
• /
• 1996

본 연구에서는 시 공간적으로 변화하는 밀도구조를 고려한 3차원 연직함수 전개모형을 제안한다. 정수압을 가정하며 열보존식과 상태방정식의 도입으로 밀도변화가 고려된다 수평방향으로의 변화는 유한 차분 격자상에서 계산되며 유속 및 수온의 연직구조는 선형보간함수를 사용하여 계산된다. 구성된 행렬방정식은 유사변환 개념을 도입하여 시간적분된다. 개발된 모델의 테스트를 위해 간단한 이상해역과 황해역에서 대기와 해양간의 열교환에 따른 수온구조 변동을 실험하였다. 이류효과는 열수송방정식에서만 고려하였으며 연직 와동점성계수와 와동확산계수는 시 공간적으로 일정한 값을 사용하였다. 이상해역에서의 수치실험결과 모델영역의 수심의 차이에 따른 열저장의 차이로 인해 수온의 수평적 구배가 발생하였다. 결과적으로 전향력을 고려하지 않을 경우에는 상층에서는 수심이 증가하는 방향으로 흐름이 발생하고 하층에는 반대방향의 흐름이 유도된 반면 전향력을 고려할 경우에는 수온차에 의한 압력구배력과 전향력이 균형을 이루면서 지형류가 뚜렷하게 나타났다. 황해역에서는 복잡한 흐름이 나타났지만 전체적으로는 지형류의 특성이 우세하게 나타났다.

• 한국환경과학회지
• /
• 제13권8호
• /
• pp.721-731
• /
• 2004

The dispersion of suspended particulates in the coastal complex terrain of mountain-inland basin (city)-sea, considering their recycling was investigated using three-dimensional non-hydrostatic numerical model and lagrangian particle model (or random walk model). Convective boundary layer under synoptic scale westerly wind is developed with a thickness of about I km over the ground in the west of the mountain, while a thickness of thermal internal boundary layer (TIBL) is only confined to less than 200m along the eastern slope of the mountain, below an easterly sea breeze circulation. At the mid of the eastern slop of the mountain, westerly wind confronts easterly sea breeze, which goes to the height of 1700 m above sea level and is finally eastward return flow toward the sea. At this time, particulates floated from the ground surface of the city to the top of TIBL go along the eastern slope of the mountain in the passage of sea breeze, being away the TIBL and reach near the top of the mountain. Then those particulates disperse eastward below the height of sea-breeze circulation and widely spread out over the coastal sea. Total suspended particulate concentration near the ground surface of the city is very low. On the other hand, nighttime radiative cooling produces a shallow nocturnal surface inversion layer (NSIL) of 200 m thickness over the inland surface, but relatively thin thickness less than 100m is found near the mountain surface. As synoptic scale westerly wind should be intensified under the association of mountain wind along the eastern slope of mountain to inland plain and further combine with land-breeze from inland plain toward sea, resulting in strong wind as internal gravity waves with a hydraulic jump motion bounding up to about 1km upper level in the atmosphere in the west of the city and becoming a eastward return flow. Simultaneously, wind near the eastern coastal side of the city was moderate. Since the downward strong wind penetrated into the city, the particulate matters floated near the top of the mountain in the day also moved down along the eastern slope of the mountain, reaching the. downtown and merging in the ground surface inside the NSIL with a maximum ground level concentration of total suspended particulates (TSP) at 0300 LST. Some of them were bounded up from the ground surface to the 1km upper level and the others were forward to the coastal sea surface, showing their dispersions from the coastal NSIL toward the propagation area of internal gravity waves. On the next day at 0600 LST and 0900 LST, the dispersed particulates into the coastal sea could return to the coastal inland area under the influence of sea breeze and the recycled particulates combine with emitted ones from the ground surface, resulting in relatively high TSP concentration. Later, they float again up to the thermal internal boundary layer, following sea breeze circulation.

• 대한토목학회논문집
• /
• 제28권5B호
• /
• pp.575-589
• /
• 2008

본 연구는 지진해일단파(tsunami bore) 혹은 조석단파(tidal bore)와 같은 단파의 동수역학적인 거동특성을 검토할 목적으로, 댐파괴류에서 단파의 형성과 동일한 방법, 즉 수조의 한쪽 끝단에 있는 고수위의 저수조(貯水槽) 게이트를 순간적으로 제거하는 방법으로 단파를 발생시킨다. 이러한 단파의 형성과 전파에 관한 수치시뮬레이션에 이상유(二相流)모델에 기초한 Navier-Stokes식을 적용하며, 이 때 비압축성 및 비혼합성의 액체와 기체흐름을 각각 고려한다. 기체와 액체의 접면을 VOF법으로 추적하고, Navier-Stokes방정식을 수치적으로 풀기 위하여 MCIP법을 적용한다. 1차원인 CIP법을 분할스텝기법을 사용하여 고차원으로 확장한 MCIP법은 수치확산이 매우 작고, 또한 안정된 스킴으로 알려져 있다. 게다가, 난류를 시뮬레이션하기 위하여 그의 유용성이 잘 알려져 있는 LES모델을 사용한다. 단파의 형성과 전파에 관한 수치해석결과를 검증하기 위하여 수리실험을 수행하였으며, 시간경과에 따른 수위변동과 평균유속변동에 대한 수치해석결과 및 실험결과를 비교하여 매우 양호한 상호대응관계를 확인할 수 있었다.

• 한국해안·해양공학회논문집
• /
• 제25권5호
• /
• pp.335-347
• /
• 2013

2008년 2월 일본 홋카이도 서해상의 발달된 저기압에 의해 생성된 폭풍파랑이 동해상 남/남서쪽으로 전파되어 한국과 일본의 동해 해안을 따라 상당한 인명 및 재산 피해를 입혔다. 본 연구는 두 파트로 구성되어 있다. 첫번째 파트에서는 연안역을 따라 상당한 피해를 입은 일본 토야마만에서의 극한 폭풍파랑을 추산하였다. 추산방법으로는 풍파의 성장발달에 중요한 요소인 바람의 강도와 계속 시간의 극한조건을 산정 후, 극한조건을 적용한 동계 온대저기압 상황을 비정역학 기상모델과 스펙트럼 파랑모델을 이용한 수치 실험을 통해 추산하였다. 추산된 토야마만 후시키 토야마에서의 극한 폭풍파랑의 유의파고 및 주기는 각각 6.78 m와 18.28 sec이다. 두 번째 파트에서는 2008년 2월 폭풍파랑으로 인해 북방파제 및 항구에 상당한 피해를 입은 토야마만 후시키항에서의 파랑-구조물 상호작용에 관한 수치실험을 수행하였다. 수치실험은 적합격자세분화 및 wet-dry법이 적용된 비선형천수방정식 모델을 이용하였다. 첫 파트에서 추산된 폭풍파랑 특성은 파랑-구조물 상호작용 수치실험에서 입사파 조건으로 사용되었다. 수치실험 결과, 후시키항의 북방파제가 폭풍파랑에 의해 파손 시, 배후의 만요우부두는 월파 및 월류에 안전하지 못 함이 파악되었다. 또한, 추산 폭풍파랑 상황 하에서 만요우부두의 현 호안시설로는 측면 호안벽으로부터의 월류에 대응하지 못 함이 파악되었다. 두 번째 수치실험결과로부터, wet-dry법이 적용된 적합격자세분화에 의해 세분화된 격자는, 계산부하를 효율적으로 유지하는 동시에, 해안선의 표현 및 해안구조물의 표현에 뛰어남을 확인하였다.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• 제25권1호
• /
• pp.57-81
• /
• 1997

Numerical prediction of nocturnal thermal high in summer of the 1995 near Taegu city located in a basin has been carried out by a non-hydrostatic numerical model over complex terrain through one-way double nesting technique in the Z following coordinate system. Under the prevailing westerly winds, vertical turbulent fluxes of momentum and heat over mountains for daytime hours are quite strong with a large magnitude of more than $120W/\textrm{m}^2$, but a small one of $5W/\textrm{m}^2$ at the surface of the basin. Convective boundary layer (CBL) is developed with a thickness of about 600m over the ground in the lee side of Mt. Hyungje, and extends to the edge of inland at the interface of land sea in the east. Sensible heat flux near the surface of the top of the mountain is $50W/\textrm{m}^2$, but its flux in the basin is almost zero. Convergence of sensible heat flux occurs from the ground surface toward the atmosphere in the lower layer, causing the layer over the mountain to be warmed up, but no convergance of the flux over the basin results from the significant mixing of air within the CBL. As horizontal transport of sensible heat flux from the top of the mountain toward over the basin results in the continuous accumulation of heat with time, enhancing air temperature at the surface of the basin, especially Taegu city to be higher than $39.3^{\circ}C$. Since latent heat fluxes are $270W/\textrm{m}^2$ near the top of the mountain and $300W/\textrm{m}^2$ along the slope of the mountain and the basin, evaporation of water vapor from the surface of the basin is much higher than one from the mountain and then, horizontal transport of latent heat flux is from the basin toward the mountain, showing relative humidity of 65 to 75% over the mountain to be much greater than 50% to 55% in the basin. At night, sensible heat fluxes have negative values of $-120W/\textrm{m}^2$ along the slope near the top of the mountain and $-50W/\textrm{m}^2$ at the surface of the basin, which indicate gain of heat from the lower atmosphere. Nighttime radiative cooling produces a shallow nocturnal surface inversion layer with a thickness of about 100m, which is much lower than common surface inversion layer, and lifts extremely heated air masses for daytime hours, namely, a warm pool of $34^{\circ}C$ to be isolated over the ground surface in the basin. As heat transfer from the warm pool in the lower atmosphere toward the ground of the basin occurs, the air near the surface of the basin does not much cool down, resulting in the persistence of high temperature at night, called nocturnal thermal high or tropical night. High relative humidity of 75% is found at the surface of the basin under the moderate wind, while slightly low relative humidity of 60% is along the eastern slope of the high mountain, due to adiabatic heating by the srong downslope wind. Air temperature near the surface of the basin with high moisture in the evening does not get lower than that during the day and the high temperature produces nocturnal warming situation.