• Journal of the Society of Naval Architects of Korea
• /
• v.47 no.1
• /
• pp.11-19
• /
• 2010

A code is developed to simulate incompressible free surface flows using the Roe's flux-difference splitting scheme. An interface of two fluids is considered as a moving contact discontinuity. The continuities of pressure and normal velocity across the interface are enforced by the conservation law in the integral sense. The fluxes are computed using the Roe's flux-difference splitting scheme for two incompressible fluids. The interface can be identified based on the computed density distribution. However, no additional treatment is required along the interface during the whole computations. Complicated time evolution of the interface including topological change can be captured without any difficulties. The developed code is applied to simulate the Rayleigh-Taylor instability of two incompressible fluids in the density ratio of 7.2:1 and the broken dam problem of water-air. The present results are compared with other available results and good agreements are achieved for the both cases.

• Proceedings of The Korean Society of Agricultural and Forest Meteorology Conference
• /
• 2019.08a
• /
• pp.118-119
• /
• 2019

• Journal of The Korean Astronomical Society
• /
• v.32 no.2
• /
• pp.137-147
• /
• 1999

Recent studies show the importance of understanding three-dimensional magnetic reconnect ion on the solar surface. For this purpose, I consider non-coplanar magnetic reconnection, a simple case of three-dimensional reconnect ion driven by a collision of two straight flux tubes which are not on the same plane initially. The relative angle e between the two tubes characterizes such reconnection, and can be regarded as a measure of magnetic shear. The observable characteristics of non-coplanar reconnection are compared between the two cases of small and large angles. An important feature of the non-coplanar reconnect ion is that magnetic twist can be produced via the re-ordering of field lines. This is a consequence of the conversion of mutual helicity into self helicities by reconnection. It is shown that the principle of energy conservation when combined with the production of magnetic twist puts a low limit on the relative angle between two flux tubes for reconnect ion to occur. I provide several observations supporting the magnetic twist generation by reconnection, and discuss its physical implications for the origin of magnetic twist on the solar surface and the problem of coronal heating.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Journal of Advanced Marine Engineering and Technology
• /
• v.22 no.5
• /
• pp.693-701
• /
• 1998

This study presents result of comparison between two other numerical method method of char-acteristics(MOC) and Lax-Wendroff method(LWM) applied at wave action analysis of intake and exhaust pipe in terms of calculated pressure velocity and emitted noise in the time and the fre-quency domain by means of fast Fourier transform analysis. Particularly FCT(Flux Corrected Transport)scheme is appended to LWM to protest unaceptable overshoots occurring near discon-tinuity. The final conclusion of this study is that MOC should be replaced by a second order finite difference approach because of larger contributions due to high frequency components than the results from the method of characteristics. Clear benefits we can get by change are faster calcula-tion higher accuracy conservation of mass and consistent calculation method.

• Journal of Power System Engineering
• /
• v.3 no.3
• /
• pp.14-21
• /
• 1999

The wave nature of heat conduction has been developed in situations involving extreme thermal gradients, very short times, or temperatures near absolute zero. Under the excitation of a periodic surface heating in a finite medium, the hyperbolic and parabolic heat conduction equations and the damped wave equations in heat flux are presented for comparative analysis by using the Green's function with the integral transform technique. The Kummer transformation is also utilized to accelerate the rate of convergence of these solutions. On the other hand, the temperature distributions are obtained through integration of the energy conservation law with respect to time. For hyperbolic heat conduction, the heat flux distribution does not exist throughout all the region in a finite medium within the range of very short times(${\xi}<{\eta}_l$). It is shown that due to the thermal relaxation time, the hyperbolic heat conduction equation has thermal wave characteristics as the damped wave equation has wave nature.

• 한국전산유체공학회:학술대회논문집
• /
• 1999.05a
• /
• pp.186-191
• /
• 1999

An Euler solver is developed to predict accurate aerodynamic data such as lift coefficient, drag coefficient, and moment coefficient. The conservation law form of the compressible Euler equations are used in the generalized curvilinear coordinates system. The Euler solver uses a finite volume method and the second order Roe's flux difference splitting scheme with min-mod flux limiter to calculate the fluxes accurately. An implicit scheme which includes the boundary conditions is implemented to accelerate the convergence rate. The multi-block grid is integrated into the flow solver for complex geometry. The flowfields are analyzed around NACA 0012 airfoil in the cases of $M_{\infty}=0.75,\;\alpha=2.0\;and\;M_{\infty}=0.80,\;\alpha=1.25$. The numerical results are compared with other numerical results from the literature. The final goal of this research is to prepare a robust and an efficient Navier-Stokes solver eventually.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.2
• /
• pp.7-12
• /
• 2008

A far field solution of the slowly varying force on an offshore structure by gravity ocean waves was shown as a function of the reflection and transmission of the body disturbed waves. The solution was obtained from the conservation of the momentum flux, which simply describes various wave forces, while making it unnecessary to compute complicated integration over a control surface. The solution was based on the assumption that the frequency difference of the bichromatic incident waves is small and its second order term is negligible. The final solution is expressed in term of the reflection and transmission waves, i.e. their amplitudes and phase angles. Consequently, it shows that not only the amplitudes but also the phase differences make critical contributions to the slowly varying force. In a limiting case, the slowly varying force solution gives the one of the mean drift force, which is only dependent on the reflection wave amplitude. An approximation is also suggested in a case where only the mean drift force information is available.

• Ocean and Polar Research
• /
• v.46 no.1
• /
• pp.93-103
• /
• 2024

Estuaries are connecting pathways where terrestrial carbon is transported to the ocean and environments where various biogeochemical cycles occur. It is essential to estimate the carbon flux across the land-sea continuum to accurately determine the global carbon budget. Additionally, understanding the carbon characteristics of estuarine environments provides valuable information for watershed management and coastal ecosystem conservation. This paper introduces research results in Korea regarding the organic carbon cycle in estuarine environments. In particular, it focuses on research results concerning organic carbon characteristics using stable and radioisotopes, and, based on this, suggests directions for future study.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.589-602
• /
• 2005

We present the numerical algorithm for the model for high-strain rate deformation in hyperelastic-viscoplastic materials based on a fully conservative Eulerian formulation by Plohr and Sharp. We use a hyperelastic equation of state and the modified Steinberg and Lund's rate dependent plasticity model for plasticity. A two-dimensional approximate Riemann solver is constructed in an unsplit manner to resolve the complex wave structure and combined with the second order TVD flux. Numerical results are also presented.