• Journal of the Mineralogical Society of Korea
• /
• v.22 no.4
• /
• pp.313-324
• /
• 2009

We explore the effect of particle shape and size on 3-dimensional (3D) network and pore structure of porous earth materials composed of glass beads and silica gel using NMR micro-imaging in order to gain better insights into relationship between structure and the corresponding hydrologic and seismological properties. The 3D micro-imaging data for the model porous networks show that the specific surface area, porosity, and permeability range from 2.5 to $9.6\;mm^2/mm^3$, from 0.21 to 0.38, and from 11.6 to 892.3 D (Darcy), respectively, which are typical values for unconsolidated sands. The relationships among specific surface area, porosity, and permeability of the porous media are relatively well explained with the Kozeny equation. Cube counting fractal dimension analysis shows that fractal dimension increases from ~2.5－2.6 to 3.0 with increasing specific surface area from 2.5 to $9.6\;mm^2/mm^3$, with the data also suggesting the effect of porosity. Specific surface area, porosity, permeability, and cube counting fractal dimension for the natural mongolian sandstone are $0.33\;mm^2/mm^3$, 0.017, 30.9 mD, and 1.59, respectively. The current results highlight that NMR micro-imaging, together with detailed statistical analyses can be useful to characterize 3D pore structures of various porous earth materials and be potentially effective in accounting for transport properties and seismic wave velocity and attenuation of diverse porous media in earth crust and interiors.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.30 no.6
• /
• pp.306-318
• /
• 2018

In order to evaluate the seismic capacity of massive vertical type breakwaters which have intensively been deployed along the coast of South Korea over the last two decades, we carry out the preliminary numerical simulation against the PoHang, GyeongJu, Hachinohe 1, Hachinohe 2, Ofunato, and artificial seismic waves based on the measured time series of ground acceleration. Numerical result shows that significant sliding can be resulted in once non-negligible portion of seismic energy is shifted toward the longer period during its propagation process toward the ground surface in a form of shear wave. It is well known that during these propagation process, shear waves due to the seismic activity would be amplified, and non-negligible portion of seismic energy be shifted toward the longer period. Among these, the shift of seismic energy toward the longer period is induced by the viscosity and internal friction intrinsic in the soil. On the other hand, the amplification of shear waves can be attributed to the fact that the shear modulus is getting smaller toward the ground surface following the descending effective stress toward the ground surface. And the weakened intensity of soil as the number of attacking shear waves are accumulated can also contribute these phenomenon (Das, 1993). In this rationale, we constitute the numerical model using the model by Hardin and Drnevich (1972) for the weakened shear modulus as shear waves go on, and shear wave equation, in the numerical integration of which $Newmark-{\beta}$ method and Modified Newton-Raphson method are evoked to take nonlinear stress-strain relationship into account. It is shown that the numerical model proposed in this study could duplicate the well known features of seismic shear waves such as that a great deal of probability mass is shifted toward the larger amplitude and longer period when shear waves propagate toward the ground surface.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.5 no.1
• /
• pp.21-30
• /
• 1985

The tension leg platform (TLP) is a kind of compliant structures, and is also a type of moored stable platform with a buoyancy exceeding the weight because of having tensioned vertical anchor cables. In this paper, among the various kinds of tension leg structures, Deep Oil Technology (DOT) TLP was analyzed because it has large-displacement portions of the immersed surface such as vertical corner pontoons and small-diameter elongated members such as cross-bracing. It also has results of hydraulic model tests, comparable with theorectical analysis. Because of the vertical axes of symmetry in the three vertical buoyant legs and because there are no larger horizontal buoyant members between these three vertical members, it was decided to develop a numerical algorithm which would predict the dynamic response of the DOT TLP using the previously developed numerical algorithm Floating Vessel Response Simulation (FVRS) for vertically axisymmetric bodies of revolution. In addition, a linearized hydroelastic Morison equation subroutine would be developed to account for the hydrodynamic pressure forces on the small member cross bracing. Interaction between the large buoyant members or small member cross bracings is considered to be negligible and is not included in the analysis. The dynamic response of the DOT TLP in the surge mode is compared with the results of the TLP algorithm for various combinations of diffraction and Morison forces and moments. The results which include the Morison equation are better than the results for diffraction only. This is because the vertically axisymmetric buoyant members are only marginally large enough to consider diffractions effects. The prototype TLP results are expected to be more inertially dominated.

• Korea Journal of Hospital Management
• /
• v.27 no.2
• /
• pp.53-66
• /
• 2022

Aim : The purpose of this study is to investigate the influence of the awareness level of the mandatory employment system for the disabled and the help of the employment increase on the health status of the disabled. Based on this, the purpose of this study is to raise awareness of the compulsory employment system for the disabled and to suggest a way to improve the health of the disabled. Methods : In this study, the first data of the second wave of the 2016-2018 Employment Panel Survey for the Disabled was used to analyze the participation of the disabled in economic activities. a total of 1,648 subjects were analyzed using the chi-square test and generalized estimating equation (GEE). Results : Compared to the group knowing about the compulsory employment system for the disabled, in the group who did not know about the compulsory employment system, the cognitive level was significantly higher for subjective health status (Odds Ratio [OR]: 1.573, 95% Confidence Interval [CI]: 1.252-1.977) and chronic diseases (OR: 1.407, 95% CI: 1.091-1.816), Compared to the group that the compulsory employment system for the disabled is helpful in increasing employment, the group that said it was not helpful at all had depression feel (OR: 2.330, 95% CI: 1.219-4.452) and subjective health status (OR: 2.052, 95% CI: 1.232. -3.416) significantly affected. Conclusion : It was found that the lower the level of awareness and help of the compulsory employment system for the disabled, the negative impact on their health status. Therefore, it is necessary to promote employment growth by raising awareness of the employment system for the disabled. Furthermore, in order to improve the health of the disabled, public relations plans to increase the understanding of the disabled on the system and policies to protect the disabled people who lose their will to work due to difficulties in accessing employment were enacted, and Changes to a disability-friendly working environment should be activated so that participants can participate.

• Journal of the Korean Institute of Landscape Architecture
• /
• v.45 no.6
• /
• pp.115-125
• /
• 2017

This study was to investigate the user's thermal environments under the pergola($L\;7,200{\times}W\;4,200{\times}H\;2,700mn$) covered with Wisteria floribunda(Willd.) DC. according to the variation of leaf area index(LAI). We carried out detailed measurements with two human-biometeorological stations on a popular square Jinju, Korea($N35^{\circ}10^{\prime}59.8^{{\prime}{\prime}}$, $E\;128^{\circ}05^{\prime}32.0^{{\prime}{\prime}}$, elevation: 38m). One of the stations stood under a pergola, while the other in the sun. The measurement spots were instrumented with microclimate monitoring stations to continuously measure air temperature and relative humidity, wind speed, shortwave and longwave radiation from the six cardinal directions at the height of 0.6m so as to calculate the Universal Thermal Climate Index(UTCI) from $9^{th}$ April to $27^{th}$ September 2017. The LAI was measured using the LAI-2200C Plant Canopy Analyzer. The analysis results of 18 day's 1 minute term human-biometeorological data absorbed by a man in sitting position from 10am to 4pm showed the following. During the whole observation period, daily average air temperatures under the pergola were respectively $0.7{\sim}2.3^{\circ}C$ lower compared with those in the sun, daily average wind speed and relative humidity under the pergola were respectively 0.17~0.38m/s and 0.4~3.1% higher compared with those in the sun. There was significant relationship in LAI, Julian day number and were expressed in the equation $y=-0.0004x^2+0.1719x-11.765(R^2=0.9897)$. The average $T_{mrt}$ under the pergola were $11.9{\sim}25.4^{\circ}C$ lower and maximum ${\Delta}T_{mrt}$ under the pergola were $24.1{\sim}30.2^{\circ}C$ when compared with those in the sun. There was significant relationship in LAI, reduction ratio(%) of daily average $T_{mrt}$ compared with those in the sun and was expressed in the equation $y=0.0678{\ln}(x)+0.3036(R^2=0.9454)$. The average UTCI under the pergola were $4.1{\sim}8.3^{\circ}C$ lower and maximum ${\Delta}UTCI$ under the pergola were $7.8{\sim}10.2^{\circ}C$ when compared with those in the sun. There was significant relationship in LAI, reduction ratio(%) of daily average UTCI compared with those in the sun and were expressed in the equation $y=0.0322{\ln}(x)+0.1538(R^2=0.8946)$. The shading by the pergola covered with vines was very effective for reducing daytime UTCI absorbed by a man in sitting position at summer largely through a reduction in mean radiant temperature from sun protection, lowering thermal stress from very strong(UTCI >$38^{\circ}C$) and strong(UTCI >$32^{\circ}C$) down to strong(UTCI >$32^{\circ}C$) and moderate(UTCI >$26^{\circ}C$). Therefore the pergola covered with vines used for shading outdoor spaces is essential to mitigate heat stress and can create better human thermal comfort especially in cities during summer. But the thermal environments under the pergola covered with vines during the heat wave supposed to user "very strong heat stress(UTCI>$38^{\circ}C$)". Therefore users must restrain themselves from outdoor activities during the heat waves.

• Atmosphere
• /
• v.23 no.3
• /
• pp.293-305
• /
• 2013

Nonhydrostatic effects on convectively forced mesoscale flows in two dimensions are numerically investigated using a nondimensional model. An elevated heating that represents convective heating due to deep cumulus convection is specified in a uniform basic flow with constant stability, and numerical experiments are performed with different values of the nonlinearity factor and nonhydrostaticity factor. The simulation result in a linear system is first compared to the analytic solution. The simulated vertical velocity field is very similar to the analytic one, confirming the high accuracy of nondimensional model's solutions. When the nonhydrostaticity factor is small, alternating regions of upward and downward motion above the heating top appear. On the other hand, when the nonhydrostaticity factor is relatively large, alternating updraft and downdraft cells appear downwind of the main updraft region. These features according to the nonhydrostaticity factor appear in both linear and nonlinear flow systems. The location of the maximum vertical velocity in the main updraft region differs depending on the degrees of nonlinearity and nonhydrostaticity. Using the Taylor-Goldstein equation in a linear, steady-state, invscid system, it is analyzed that evanescent waves exist for a given nonhydrostaticity factor. The critical wavelength of an evanescent wave is given by ${\lambda}_c=2{\pi}{\beta}$, where ${\beta}$ is the nonhydrostaticity factor. Waves whose wavelengths are smaller than the critical wavelength become evanescent. The alternating updraft and downdraft cells are formed by the superposition of evanescent waves and horizontally propagating parts of propagating waves. Simulation results show that the horizontal length of the updraft and downdraft cells is the half of the critical wavelength (${\pi}{\beta}$) in a linear flow system and larger than ${\pi}{\beta}$ in a weakly nonlinear flow system.

• Journal of Navigation and Port Research
• /
• v.26 no.1
• /
• pp.106-111
• /
• 2002

The estimation of pile bearing capacity is important since the design details are determined from the result. There are numerous ways of determining the pile design load, but only few of them are chosen in the actual design. According to the recent investigation in Korea, the formulas proposed by Meyerhof based on the SPT N values are most frequently chosen in the design stage. In the study, various static and dynamic formulas have been used in predicting the allowable bearing capacity of a pile. Further, the reliability of these formulas has been verified by comparing the perdicted values with the static and dynamic load test measurements. Also, in most cases, these methods of pile bearing capacity determination do not take the time effect consideration, the actual allowable load as determined from pile load test indicates severe deviation from the design value. The principle results of this study are summarized as follows : As a result of estimate the reliability in criterion of the Davisson method, t was showed that Terzaghi & Peck >Chin>Meyerhof > Modified Meyerhof method was the most reliable method for the prediction of bearing capacity. Comparisons of the various pile-driving formulas showed that Modified Engineering News was the most reliable method. However, a significant error happened between dynamic bearing capacity equation was judged that uncertainty of hammer efficiency, characteristics of variable, time effect etc... was not considered. As a result of considering time effect increased skin friction capacity higher than end bearing capacity. It was found out that it would be possible to increase the skin friction capacity 1.99 times higher than a driving. As a result of considering 7 day's time effect, it was obtained that Engineering news, Modified Engineering News, Hiley, Danish, Gates, CAPWAP(CAse Pile Wave Analysis Program) analysis for relation, repectively, $Q_{u(Restrike)} / Q_{u(EOID)} = 0.98t_{0.1}$ , $0.98t_{0.1}$, $1.17t_{0.1}$, $0.88t_{0.1}$, $0.89t_{0.1}$, $0.97t_{0.1}$.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.21 no.2
• /
• pp.131-136
• /
• 1985

The calculation of the resulting fluid motion is an important problem of ship hydrodynamics. For a partially immersed body the condition of constant pressure at the free surface can be linearized. The resulting linear boundary-value problem for the velocity potential is the Neumann-Kelvin problem. The two-dimensional Neumann-Kelvin problem is studied for the half-immersed circular cylinder by Ursell. Maruo introduced a slender body approach to simplify the Neumann-Kelvin problem in such a way that the integral equation which determines the singularity distribution over the hull surface can be solved by a marching procedure of step by step integration starting at bow. In the present pater for the two-dimensional Neumann-Kelvin problem, it has been suggested that any solution of the problem must have singularities in the corners between the body surface and free surface. There can be infinitely many solutions depending on the singularities in the coroners.

• Proceedings of the Korean Society of Propulsion Engineers Conference
• /
• 2003.05a
• /
• pp.91-93
• /
• 2003

A comprehensive numerical study is carried out to investigate for the understanding of the flow evolution and flame development in a supersonic combustor with normal injection of ncumally injecting hydrogen in airsupersonic flows. The formulation treats the complete conservation equations of mass, momentum, energy, and species concentration for a multi-component chemically reacting system. For the numerical simulation of supersonic combustion, multi-species Navier-Stokes equations and detailed chemistry of H2-Air is considered. It also accommodates a finite-rate chemical kinetics mechanism of hydrogen-air combustion GRI-Mech. 2.11[1], which consists of nine species and twenty-five reaction steps. Turbulence closure is achieved by means of a k-two-equation model (2). The governing equations are spatially discretized using a finite-volume approach, and temporally integrated by means of a second-order accurate implicit scheme (3-5).The supersonic combustor consists of a flat channel of 10 cm height and a fuel-injection slit of 0.1 cm width located at 10 cm downstream of the inlet. A cavity of 5 cm height and 20 cm width is installed at 15 cm downstream of the injection slit. A total of 936160 grids are used for the main-combustor flow passage, and 159161 grids for the cavity. The grids are clustered in the flow direction near the fuel injector and cavity, as well as in the vertical direction near the bottom wall. The no-slip and adiabatic conditions are assumed throughout the entire wall boundary. As a specific example, the inflow Mach number is assumed to be 3, and the temperature and pressure are 600 K and 0.1 MPa, respectively. Gaseous hydrogen at a temperature of 151.5 K is injected normal to the wall from a choked injector.A series of calculations were carried out by varying the fuel injection pressure from 0.5 to 1.5MPa. This amounts to changing the fuel mass flow rate or the overall equivalence ratio for different operating regimes. Figure 1 shows the instantaneous temperature fields in the supersonic combustor at four different conditions. The dark blue region represents the hot burned gases. At the fuel injection pressure of 0.5 MPa, the flame is stably anchored, but the flow field exhibits a high-amplitude oscillation. At the fuel injection pressure of 1.0 MPa, the Mach reflection occurs ahead of the injector. The interaction between the incoming air and the injection flow becomes much more complex, and the fuel/air mixing is strongly enhanced. The Mach reflection oscillates and results in a strong fluctuation in the combustor wall pressure. At the fuel injection pressure of 1.5MPa, the flow inside the combustor becomes nearly choked and the Mach reflection is displaced forward. The leading shock wave moves slowly toward the inlet, and eventually causes the combustor-upstart due to the thermal choking. The cavity appears to play a secondary role in driving the flow unsteadiness, in spite of its influence on the fuel/air mixing and flame evolution. Further investigation is necessary on this issue. The present study features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous works. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the underlying physical mechanisms. Much of the flow unsteadiness is not related to the cavity, but rather to the intrinsic unsteadiness in the flowfield, as also shown experimentally by Ben-Yakar et al. [6], The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The work appears to be the first of its kind in the numerical study of combustion oscillations in a supersonic combustor, although a similar phenomenon was previously reported experimentally. A more comprehensive discussion will be given in the final paper presented at the colloquium.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.31 no.3
• /
• pp.101-114
• /
• 2019

In this study, we carry out the numerical simulation to trace the yearly shoreline change of Mang-Bang beach, which is suffering from erosion problem. We obtain the basic equation (One Line Model for shoreline) for the numerical simulation by assuming that the amount of shoreline retreat or advance is balanced by the net influx of longshore and cross-shore sediment into the unit discretized shoreline segment. In doing so, the energy flux model for the longshore sediment transport rate is also evoked. For the case of cross sediment transport, the modified Bailard's model (1981) by Cho and Kim (2019) is utilized. At each time step of the numerical simulation, we adjust a closure depth according to pertinent wave conditions based on the Hallermeier's analytical model (1978) having its roots on the Shield's parameter. Numerical results show that from 2017.4.26 to 2017.10.15 during which swells are prevailing, a shoreline advances due to the sustained supply of cross-shore sediment. It is also shown that a shoreline temporarily retreats due to the erosion by the yearly highest waves sequentially occurring from mid-October to the end of October, and is followed by gradual recovery of shoreline as high waves subdue and swells prevail. It is worth mentioning that great yearly circulation of shoreline completes when a shoreline retreats due to the erosion by the higher waves occurring from mid-March to the end of March. The great yearly circulation of shoreline mentioned above can also be found in the measured locations of shoreline on 2017.4.5, 2017.9.7, 2017.11.7, 2018.3.14. However, numerically simulated amount of shoreline retreat or advance is more significant than the physically measured one, and it should be noted that these discrepancies become more substantial for the case of RUN II where a closure depth is sustained to be as in the most morphology models like the Genesis (Hanson and Kraus, 1989).