• Geosystem Engineering
• /
• v.21 no.6
• /
• pp.309-317
• /
• 2018

We employ a three-dimensional indirect boundary element method (BEM) to simulate temperature change around an underground liquefied natural gas storage cavern. The indirect BEM (IBEM) uses fictitious heat source strength on boundary elements as basic variables which are solved from equations of boundary conditions and then used to compute the temperature change at other points in the considered problem domain. The IBEM requires evaluation of singular integration for temperature change due to heat conduction from a constant heat source on a planar (triangular) region. The singularity can be eliminated by a semi-analytical integration scheme. However, it is found that the semi-analytical integration scheme yields sharp temperature gradient for points close to vertices of triangle. This affects the accuracy of heat flux, if they are evaluated by finite difference method at these points. This difficulty can be overcome by a combination of using a direct numerical integration for these points and the semi-analytical scheme for other points distance away from the vertices. The IBEM and the hybrid integration scheme have been verified with an analytic solution and then used to the application of the underground storage.

• Communications of the Korean Mathematical Society
• /
• v.27 no.4
• /
• pp.753-762
• /
• 2012

We prove some sharp extremal distance results for functions in weighted Bergman spaces on the upper halfplane. We also prove new analogous results in the context of bounded strictly pseudoconvex domains with smooth boundary.

• Bulletin of the Korean Mathematical Society
• /
• v.52 no.4
• /
• pp.1321-1325
• /
• 2015

In this paper we study the semicontinuity of the automorphism groups of domains in multi-dimensional complex space. We give examples to show that known results are sharp (in terms of the required boundary smoothness).

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.14 no.11
• /
• pp.914-922
• /
• 2002

Thermo-acoustic waves can be thermally generated in a compressible flow field by rapid heating and cooling, and chemical reaction near the boundary walls. This mechanism is very important in the space environment in which natural convection does not exist. Also this may be a significant factor for heat transfer when the fluids are close to the thermodynamic critical point. In this study, the generation and transmission characteristics of thermo-acoustic waves in an air-filled confined domain with two-step pulsed heating are studied numerically. The governing equations are discretized using control volume method, and are solved using PISO algorithm and second-order upwind scheme. For the purpose of stable solution, time step was set to the order of $1\times10_-9s,\;and\;grids\;are\;50\times2000$. Results show that temperature and pressure distributions of fluid near the boundary wall subjected to a rapid heating are increased abruptly, and the induced thermo-acoustic wave propagates through the fluid until it decays due to viscous and heat dissipation. Pressure waves have sharp front shape and decay with a long tail in the case of step heating, but these waves have sharp pin shape in the case of pulsed heating.

• Journal of the Korean Mathematical Society
• /
• v.43 no.6
• /
• pp.1143-1158
• /
• 2006

We consider the projectivization of Minkowski space with the analytic continuation of the hyperbolic metric and call this an extended hyperbolic space. We can measure the volume of a domain lying across the boundary of the hyperbolic space using an analytic continuation argument. In this paper we show this method can be further generalized to find the volume of a domain with smooth boundary with suitable regularity in dimension 2 and 3. We also discuss that this volume is invariant under the group of hyperbolic isometries and that this regularity condition is sharp.

• Proceedings of the KIEE Conference
• /
• 2006.07c
• /
• pp.1494-1495
• /
• 2006

A complete finite element analysis method for discharge onset process, which is governed and coupled by charge transport equation and electric field equation, was presented. The charge transport equation of first order was transformed into a second-order one by utilizing the artificial diffusion scheme. The two second-order equations were analyzed by the finite element formulation which is well-developed for second-order ones. The Fowler-Nordheim injection boundary condition was adopted for charge transport equation. After verifying the numerical results by comparing to the analytic solutions using parallel plane electrodes with one carrier system, we extended the result to blade-plane electrodes in 2D xy geometry with three carriers system. Radius of the sharp tip was taken to be 50 ${\mu}m$. When this sharp geometry was solved by utilizing the space discretizing methods, the very sharp tip was found to cause a singularity in electric field and space charge distribution around the tip. To avoid these numerical difficulties in the FEM, finer meshes, a higher order shape function, and artificial diffusion scheme were employed.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.1
• /
• pp.62-72
• /
• 2008

In this study, an analytical model has been developed to predict the build-up of heat field due to a point heat source in the presence of sharp cross-flow depth change. The model has been applied to investigate the effect of the depth change and flow pattern on the heat field. Model results show that, when there is a sharp depth change in cross-flow direction, the heat transport across the boundary of the depth change is enhanced or diminished according to the increasing or decreasing of the horizontal diffusion flux. Including residual components as well as tidal currents give rise to reduce the effect of the horizontal diffusion on the heat transport because of increasing the advection of heat.

• Journal of the Society of Naval Architects of Korea
• /
• v.43 no.5 s.149
• /
• pp.538-549
• /
• 2006

A code is developed to simulate a viscous flow field around a deformable body using the hybrid Cartesian/immersed boundary method. In this method, the immersed boundary(IB) nodes are defined near the body boundary then velocities at the IB nodes are reconstructed based on the interpolation along the normal direction to the body surface. A new method is suggested to define the IB nodes so that a closed fluid domain is guaranteed by a set of IB nodes and the method is applicable to a zero-thickness body such as a sail. To validate the developed code, the vorticity fields are compared with other recent calculations where a cylinder orbits and moves into its own wake. It is shown the code can handle a sharp trailing edge at Reynolds number of $10^5$ under moderate requirements on girds. Finally the developed code is applied to simulate the vortex shedding behind a deforming foil with flapping tail like a fish. It is shown that the acceleration of fluids near the flapping tail contributes to the generation of the thrust for propulsion.

• Journal of Environmental Science International
• /
• v.10 no.6
• /
• pp.417-422
• /
• 2001

The vertical structure of atmosphere was observed In investigate the variation of surface ozone concentration by vertical downward mixing of residual ozone in the atmospheric boundary layer at the Busan coastal area. Airsonde and pilot balloon measurements were made at Gamcheondong and the Kimhae airport for April 26~27, 1996. The vertical potential of potential temperature showed a residual layer between 510m and 1800m from 2100LST April 26 to 0900LST April 27. The downward mixing of ozone in the residual layer of the atmospheric boundary layer was confirmed from vertical profile of mixing ratio near 600m in the morning. The thickness of the sea breeze layer was 900m at 1500LST April 26. Thereafter, it become to be lowered with time A low level jet was measured near 900m at 0300LST on April 27 from a pibal measurement. Early morning sharp increase of surface ozone concentration at the Busan coastal area was caused by vertical downward mixing of ozone concentration rather than by photochemical reaction In the atmospheric boundary layer.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.18 no.3
• /
• pp.701-706
• /
• 1994

The discharge coefficients of critical sonic nozzles were obtained in a high pressure gas flow standard system, which was a gravitational weighing system. The discharge coefficients of critical sonic nozzle farbricated according to ISO specifications are in good agreement with ISO correlation. The discharge coefficients for small inlet radius decrease significantly as the inlet length become short due3 to separation at the sharp-edged inlet. For nozzles having long inlet radius, the effects of inlet length on the discharge coefficients were relatively small, but the effects become significant at the short inlet length. The effect of separation at the sharp-edged inlet is stronger than that of the boundary layer growth. The experimental results support that the shape of critical sonic nozzles suggested by ISO specifications is excellent.