• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.2 s.173
• /
• 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.

• Journal of Korea Water Resources Association
• /
• v.44 no.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.

• Journal of Korean Society for Atmospheric Environment
• /
• v.25 no.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.

• Publications of The Korean Astronomical Society
• /
• v.30 no.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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.8 no.2
• /
• pp.123-136
• /
• 1996

A three-dimensional Galerkin-FEM model which can handle the temporal and spatial variation of density is presented. The hydrostatic approximation is used and density effects are included by means of conservation equation of heat and the equation of state. The finite difference grids are used in the horizontal plane and a set of linear-shape functions is used for the vertical expansion. The similarity transform is introduced to solve resultant matrix equations. The proposed model was first applied to the density-driven circulation in an idealized basin in the presence of the heat exchange between the air and the sea. The advection terms in the momentum equation were ignored, while the convection terms were retained in the heat equation. Coefficients of the vertical eddy viscosity and diffusivity were fixed to be constant. Calculation in a non-rotating idealized basin shows that the difference in heat capacity with depth gives rise to the horizontal gradient of temperature. Consequently, there is a steady new in the upper layer in the direction of increasing depth with compensatory counter flow .in the lower layer. With Coriolis force, geostrophic flow was predominant due to the balance between the pressure gradient and the Coriolis force. As a test in region of irregular topography, the model is applied to the Yellow Sea. Although the resultant flow was very complex, the character of the flow Showed to be geostrophic on the whole.

• Journal of Environmental Science International
• /
• v.13 no.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.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.5B
• /
• pp.575-589
• /
• 2008

In the present work, we investigate the hydrodynamic behavior of a turbulent bore, such as tsunami bore and tidal bore, generated by the removal of a gate with water impounded on one side. The bore generation system is similar to that used in a general dam-break problem. In order to the numerical simulation of the formation and propagation of a bore, we consider the incompressible flows of two immiscible fluids, liquid and gas, governed by the Navier-Stokes equations. The interface tracking between two fluids is achieved by the volume-of-fluid (VOF) technique and the M-type cubic interpolated propagation (MCIP) scheme is used to solve the Navier-Stokes equations. The MCIP method is a low diffusive and stable scheme and is generally extended the original one-dimensional CIP to higher dimensions, using a fractional step technique. Further, large eddy simulation (LES) closure scheme, a cost-effective approach to turbulence simulation, is used to predict the evolution of quantities associated with turbulence. In order to verify the applicability of the developed numerical model to the bore simulation, laboratory experiments are performed in a wave tank. Comparisons are made between the numerical results by the present model and the experimental data and good agreement is achieved.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.25 no.5
• /
• pp.335-347
• /
• 2013

In February 2008, high storm waves due to a developed atmospheric low pressure system propagating from the west off Hokkaido, Japan, to the south and southwest throughout the East Sea (ES) caused extensive damages along the central coast of Japan and along the east coast of Korea. This study consists of two parts. In the first part, we estimate extreme storm wave characteristics in the Toyama Bay where heavy coastal damages occurred, using a non-hydrostatic meteorological model and a spectral wave model by considering the extreme conditions for two factors for wind wave growth, such as wind intensity and duration. The estimated extreme significant wave height and corresponding wave period were 6.78 m and 18.28 sec, respectively, at the Fushiki Toyama. In the second part, we perform numerical experiments on wave-structure interaction in the Fushiki Port, Toyama Bay, where the long North-Breakwater was heavily damaged by the storm waves in February 2008. The experiments are conducted using a non-linear shallow-water equation model with adaptive mesh refinement (AMR) and wet-dry scheme. The estimated extreme storm waves of 6.78 m and 18.28 sec are used for incident wave profile. The results show that the Fushiki Port would be overtopped and flooded by extreme storm waves if the North-Breakwater does not function properly after being damaged. Also the storm waves would overtop seawalls and sidewalls of the Manyou Pier behind the North-Breakwater. The results also depict that refined meshes by AMR method with wet-dry scheme applied capture the coastline and coastal structure well while keeping the computational load efficiently.

• International Union of Geodesy and Geophysics Korean Journal of Geophysical Research
• /
• v.25 no.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.