• Proceedings of the Korean Geotechical Society Conference
• /
• 2001.03a
• /
• pp.153-160
• /
• 2001

피에조콘관입시험(PCPT)의 소산시험은 in-situ 상태의 압밀계수(c/sub v/)를 추정하는 방법으로 널리 이용되어왔다. 본 연구에서는 spherical cavity expansion theory 및 axisymmetric uncoupled linear consolidation equation(Gupta & Davidson, 1986)을 이용하여 과압밀점토에서의 초기과잉간극수압의 분포 및 과잉간극수압의 시간에 대한 소산현상을 해석하는 수치해석방법을 제안하여 현장시험결과 및 실내시험결과와 비교 분석하였다. ADIS (alternating direction implicit scheme)를 이용한 FDM 해석을 실시한 결과와 현장시험의 소산곡선은 잘 부합되는 것으로 나타났으며 압밀계수도 실내실험 또는 피에조콘관입 시험에 대한 추정방법으로 산출된 값과 비교적 일치하는 것으로 나타났다.

• Journal of Korea Foundry Society
• /
• v.4 no.1
• /
• pp.21-31
• /
• 1984

IAD( Implicit Alternative Direction ) FDM을 사용하여 금형에서의 L형 순수알미늄 주물에 대하여 그 응고상황의 수치적 해석을 2차원적으로 시도하였다. 계산은 주형비 ( 주형의 부피/주물의 부피 )에 따라 금형 표면에서의 대류현상과 금형/주물계면에서의 Air-gap형성 현상을 고려에 넣은 경우와 그렇지 않은 경우로 나누어 각각 콤퓨터에 의해 계산하여 그 결과를 비교하였다. 또한 금형/주물 계면주위의 온도구배를 구함으로써 Air-gap이 열전달에 미치는 영향을 검토하였다. 이와같은 시도에 의하여 주물의 응고시간, 응고방향 Hot spot의 위치 등 주조에 있어서 매우 중요한 사항들을 알아내는데 Numerical Analysis 이 매우 유용한 방법임을, 특히 금형표면에서의 대류현상과 금형/주물 계면에서의 Air-gap형성현상을 고려에 넣은 경우에 더욱 정확한 결과를 기대할 수 있다는 것이 확인되었다.

• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.3 no.1
• /
• pp.69-77
• /
• 1991

"Ondol" have been used in residential buildings for several thousands years in Korea. The traditional "Ondol" heating system of Korea has changed into the radiant heating system with piping embedded in floors or slabs. This study reports the results of transient experiments performed on a radiant heating system and enclosure. The paper presents some details of the thermal response of slab-heated buildings to varying patterns of heat input. Furthermore, I'll compare the results of experimentation with the ones of the numerical simulation by using the explicit and implicit forms of the finite difference methods (FDM). The study has contributed to testify the feasibility of numerical analysis, and the understanding of the transient behavior of radiant heating panels and enclosure exposed to this type of heating system.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2021.06a
• /
• pp.253-253
• /
• 2021

천수방정식에 대한 불연속 갤러킨 기법 (DG) 모형은 주로 양해법 기반으로 개발되어 적용되어 왔으나, 바닥마찰항의 처리, 과도한 CFL 조건 등의 불리한 점이 지적되어 왔으며, 이로 인하여 실제 적용에서 FDM, FEM 등 다른 고전적인 수치기법과 비교하여 경쟁력을 갖기 어려웠다. 이에 대한 대안으로써, 최근, 불연속 갤러킨 기법에 대한 음해법 기반의 모형이 연구되고 있으며, 다소 복잡한 알고리즘에도 불구하고 적용이 확대되고 있다. 또한, 널리 알려진 바와 같이, 천수방정식의 실제 하도에 대한 적용에 있어 문제점 중 하나는 나비에-스토크스 방정식으로부터 유도할 때 사용된 정수압 가정으로 인하여, 하도의 계단과 같은 불연속 지형에 적용이 용이하지 않다는 것이다. 본 연구에서는 기존에 개발된 불연속 갤러킨 음해법에 불연속 지형의 해석을 위한 표면경사법(surface gradient method)을 결합하여 이러한 문제에 효과적으로 대응할 수 있는 기법을 제시하였다. 개발된 모형의 검증을 위하여, 제방 등 하도 구조물 위의 장주기 조석흐름, 홍수파, 계단 등을 포함하는 댐 붕괴류 모의에 적용하고 실용적인 기능성을 검증하였다. 향후 구조물이 많은 국내 하천에 적용이 가능할 것으로 사료된다.

• Water for future
• /
• v.23 no.3
• /
• pp.319-328
• /
• 1990

The Parameter Identification of 2-demensional estuarine model was carried out using new output ADI-FDM numerical semi-implicit schem transformed in boundary fitted(BF) - coordinate. The hydrodynamic equations which is coupled with the transport equations were used as basic equations in the model. Thompson's equations were used to transform governing equations into rectangular plane equations and his elliptic grid generation scheme was used to generate curvilinear grid system. in BF - coordinates. The parameters to be identified are friction coefficient and disperse coefficient embedded in the governing equations. The numerical output scheme is tidally averaged salinity model in BF - coordinates. The algorithm to optimize norm of error between observations and calculations is the influence coefficinet algorithm associated with least square criterion. The lumped model is conssidered in identification. This paper was concetrated on checking whether the new output scheme might be useful to identify parameters in estuarine salinity model or not. The proposed method was tested through experimental application with hypothetical simple model. The result of the test shows that the proposed method can be used for parameter identification in estuarine model.

• Journal of Korea Water Resources Association
• /
• v.45 no.2
• /
• pp.179-189
• /
• 2012

To simulate the transport of pollutant, a numeric model for the advection-diffusion equation with the decay term is developed. This is finite-difference model using the implicit method (with the weight factor ${\alpha}$) and Gauss-Seidel SOR(successive over-relaxation). This model is compared to the analytical solutions (of simpler dimensional or boundary conditions), and in the condition of Peclet number < 5~20, the result shows stable condition, and Crank-Nicolson method (${\alpha}$=0.5) shows the more accurate results than fully-implicit method (${\alpha}$=1). The mass of advection, diffusion and decay is calculated and the error of mass balance is less than 3%. This model can evaluate the 3-D concentrations of the advection-diffusion and decay problems, but this model uses only the finite-difference method with the fixd grid system, so it can be effectively used in the problems with small Peclet numbers like the pollutant transport in groundwater.

• Journal of Mechanical Science and Technology
• /
• v.16 no.8
• /
• pp.1102-1108
• /
• 2002

In this study, a topological design sensitivity of the ai. bearing surface (ABS) is suggested by using an adjoint variable method. The discrete form of the generalized lubrication equation based on a control volume formulation is used as a compatible condition. A residual function of the slider is considered as an equality constraint function, which represents the slider in equilibrium. The slider thickness parameters at all grid cells are chosen as design variables since they are the topological parameters determining the ABS shape. Then, a complicated adjoint variable equation is formulated to directly handle the highly nonlinear and asymmetric coefficient matrix and vector in the discrete system equation of air-lubricated slider bearings. An alternating direction implicit (ADI) scheme is utilized for the numerical calculation. This is an efficient iterative solver to solve large-scale problem in special band storage. Then, a computer program is developed and applied to a slider model of a sophisticated shape. The simulation results of design sensitivity analysis (DSA) are directly compared with those of FDM at the randomly selected grid cells to show the effectiveness of the proposed approach. The overall distribution of DSA results are reported, clearly showing the region on the ABS where special attention should be given during the manufacturing process.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.3 no.2
• /
• pp.1-4
• /
• 1999

The Helmholtz equation is very important in physics and engineering. However, solution of the Helmholtz equation is in general known as a very difficult phenomenon. For if the ${\omega}$ is negative, the FDM discretized linear system becomes indefinite, whose solution by iterative method requires a very clever preconditioner. In this paper we assume that ${\omega}$ is nonnegative, and determine the optimal ${\rho}$ parameter for the three dimensional ADI iteration for the Helmholtz equation. The ADI(Alternating Direction Implicit) method is also getting new attentions due to the fact that it is very suitable to the vector/parallel computers, for example, as a preconditioner to the Krylov subspace methods. However, classical ADI was developed for two dimensions, and for three dimensions it is known that its convergence behaviour is quite different from that in two dimensions. So far, in three dimensions the so-called Douglas-Rachford form of ADI was developed. It is known to converge for a relatively wide range of ${\rho}$ values but its convergence is very slow. In this paper we determine the necessary conditions of the ${\rho}$ parameter for the convergence and optimal ${\rho}$ for the three dimensional ADI iteration of the Peaceman-Rachford form for the real Helmholtz equation with nonnegative ${\omega}$. Also, we conducted some experiments which is in close agreement with our theory. This straightforward extension of Peaceman-rachford ADI into three dimensions will be useful as an iterative solver itself or as a preconditioner to the the Krylov subspace methods, such as CG(Conjugate Gradient) method or GMRES(m).

• Transactions of the Society of Information Storage Systems
• /
• v.5 no.1
• /
• pp.19-35
• /
• 2009

The present work is an experimental and analytical study on a flexible disk rotating close to a rigid rotating disk in open air. In the analytical study, the air flow in the gap between the flexible disk and the rigid disk is modeled using Navier-Stokes and continuity equations while the flexible disk is modeled using the linear plate theory. The flow equations are discretized using the cell centered finite volume method (FVM) and solved numerically with semi-implicit pressure-linked equations (SIMPLE algorithm). The spatial terms in the disk equation are discretized using the finite difference method (FDM) and the time integration is performed using fourth-order Runge-Kutta method. An experimental test-rig is designed to investigate the dynamics of the flexible disk when rotating close to a co-rotating, a counter-rotating and a fixed rigid disk, which works as a stabilizer. The effects of rotational speed, initial gap height and inlet-hole radius on the flexible disk displacement and its vibration amplitude are investigated experimentally for the different types of stabilizer. Finally, the analytical and experimental results are compared.

• Steel and Composite Structures
• /
• v.49 no.5
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.