• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.51-56
• /
• 2000

The Finite Element Solutions Is reported on solid-liquid phase change in porous media with natural convection including freezing. The model is based on volume averaged transport equations, while phase change is assumed to occur over a small temperature range. The FEM (Finite Element Method) algorithm used in this study is 3-step time-splitting method which requires much less execution time and computer storage the velocity-pressure integrated method and the penalty method. And the explicit Lax-Wendroff scheme is applied to nonlinear convective term in the energy equation. For natural convection including melting and solidification the numerical results show reasonable agreement with FDM (Finite Difference Method) results.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.2 s.28
• /
• pp.97-110
• /
• 1998

The numerical damping and dispersion error characteristics associated with difference schemes and a panel shift method used for the calculation of steady free surface flows by a panel method are an analysed in this paper. First, 12 finite difference operators used for the double model flow by Letcher are applied to a two dimensional cylinder with the Kelvin free surface condition and the numerical errors with these schemes are compared with those by the panel shift method. Then, 3-D waves due to a submerged source are calculated by the difference schemes, the panel shift method and also by a higher order boundary element method(HOBEM). Finally, the waves and wave resistance for Wigley's hull are calculated with these three schemes. It is shown that the panel shift method is free of numerical damping and dispersion error and performs better than the difference schemes. However, it can be concluded that the HOBEM also free of the numerical damping and dispersion error is the most stable, accurate and efficient.

• Journal of Industrial Technology
• /
• v.21 no.B
• /
• pp.175-185
• /
• 2001

A numerical model is represented to calculate the reflected waves, the runup of waves and the wave induced velocities on impermeable slopes for the normally incident wave trains of nonlinear monochromatic wave and solitary wave. The finite amplitude shallow water equations with the effects of bottom friction are solved numerically in time domain using an explicit dissipative Lax-Wendroff finite difference method. The numerical model is verified by comparisons with the other numerical results, the measured data and asymptotic results. It is found that the uprushing and downrushing of incident waves may be accurately predicted by the present numerical model. Therefore, the present numerical model can be applicable to swells as well as long waves.

• Journal of applied mathematics & informatics
• /
• v.40 no.1_2
• /
• pp.25-45
• /
• 2022

In this study, the non-standard finite difference method for the numerical solution of singularly perturbed parabolic reaction-diffusion subject to Robin boundary conditions has presented. To discretize temporal and spatial variables, we use the implicit Euler and non-standard finite difference method on a uniform mesh, respectively. We proved that the proposed scheme shows uniform convergence in time with first-order and in space with second-order irrespective of the perturbation parameter. We compute three numerical examples to confirm the theoretical findings.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.178-178
• /
• 2003

• Proceedings of the Korean Geotechical Society Conference
• /
• 2006.03a
• /
• pp.1165-1172
• /
• 2006

In this study, an aspect of settlement, developed from different ground improvement method like suction drain method using vacuum pressure and vertical drain method using overburden pressure, was compared each other. In order to analyze settlement tendency of each method exactly, the finite element analysis program was used. The analyses of vertical settlement and lateral displacement for suction drain method and vertical drain method were conducted independently during the solving stage. The initial condition of drainage zone was fixed with 25m depth and 21m width. After the program analyses, the settlement condition had a different tendency with the ground improvement method. Especially, in the results of vertical drain method, the disparity of settlement between the middle of improved zone and unimproved zone. In the case of suction drain method, however, the difference of settlement was smaller than that of vertical drain method.

• Proceedings of the Korean Geotechical Society Conference
• /
• 1999.03a
• /
• pp.349-356
• /
• 1999

In general, the term soft ground includes clayey soils, which have large compressibility and small shear resistance due to the external load. All process of consolidation in compressible soils can be explained in terms of a transfer of load from an incompressible pore-water to a compressible soil structure. Therefore, one of the most important subjects about the characteristics of the time-dependent consolidation of the clay foundation by the change of load may be the presumption of the final settlement caused by consolidation and the degree of consolidation according to the time. The problems of discontinuous layer interface are very important in the algorithm and programming for the analysis of multi-layered soils using a numerical analysis, finite difference method. Better results can be obtained by the Process for discontinuous layer interface, since it can help consolidation analysis to model the actual ground. The purpose of this paper Provides an efficient computer algorithm based on numerical analysis using finite difference method(F.D.M.) which account for multi-layered soils to determine the degree of consolidation and excess pore pressures relative to time and positions more realistically.

• Journal of the Korean earth science society
• /
• v.36 no.5
• /
• pp.476-483
• /
• 2015

The high-order Laplacian-type filter, which is capable of providing isotropic and sharp cut-off filtering on the spherical domain, is essential in processing geophysical data. In this study, a spherical high-order filter was designed by combining the Fourier method with finite difference-method in the longitude and latitude, respectively. The regular grid system was employed in the filter, which has uniform angular spacing including the poles. The singularity at poles was eliminated by incorporating variable transforms and continuity-matching boundary conditions across poles. The high-order filter was assessed using the Rossby-Haurwitz wave, the observed geopotential, and observed wind field. The performance of the filter was found comparable to the filter based on the Galerkin procedure. The filter, employing the finite difference method, can be designed to give any target order of accuracy, which is an important advantage being unavailable in other methods. The computational complexity is represented with 2n-1 diagonal matrices solver with n being the target order of accuracy. Along with the availability of arbitrary target-order, it is also advantageous that the filter can adopt the reduced grid to increase computational efficiency.

• Geophysics and Geophysical Exploration
• /
• v.5 no.2
• /
• pp.99-107
• /
• 2002

An analytical forward modeling algorithm was developed for the efficient application to the geotechnical engineering in SASW (Spectral Analysis of Surface Waves) method. for the theoretical dispersion curve, the finite difference method using motion stress vector, which was proposed by Aki and Richards, was employed and verified with two earth models. For the stable and fast calculation, it was found that the model size depending on the frequency range is suitable $1.5\~2$ times bigger than the wavelength.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1997.04a
• /
• pp.192-198
• /
• 1997

In the concurrent engineering, the CAD-based design model is necessary for multidisciplinary analysis and for computer-aided manufacturing (CAM). A shape and configuration design velocity field computation of structure has been developed using a computer-aided design (CAD) tool, Pro/ENGINEER. The design Parameterization with CAD tool is to characterize the change in dimensions and movements of geometric control points that govern the shape/orientation of the structural boundary. The boundary velocity is obtained by using a CAD-based finite difference method and the domain velocity field is obtained from finite element analysis (FEA) using the boundary displacement method. In this paper, the continuum configuration DSA for NVH problem, which requires the shape velocity field and the orientation velocity field at the same time, is developed using linear shape functions. For validation of continuum design sensitivity coefficients, design sensitivity coefficients are compared with the coefficients computed using by the finite difference method.