• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.10
• /
• pp.89-101
• /
• 1995

선형 탄성 등방성 물체 내에 있는 일반적인 복합모드 크랙 문제들을 해석하기 위한 이중 경계적분방정식의 일반식과 계산해법이 제시되었다. 크랙면이 포함된 물체 해석에 있어서 유일한 해를 얻기 위하여, 한 면상의 점에는 변위 경계적분방정식이 적용되었고 마주하고 있는 상대면 상의 점에는 인력 경계적분방정식이 적용되었다. 인력 및 변위 경계적분방정식의 강특이해 및 초특이해 적분항들은 수치해법을 적용하기 전에 정상화되었다. 정상화과정 중 보정되는 강특이적분항이 상대 크랙면 상의 특이해 요소를 따라 직접 적분되는 것을 격리시키기 위하여, 특이해 적분 경로를 완만한 곡면으로 우회시킨 가상의 비특이해 보조경계로 대치하여 적분값을 계산하였다. 제시된 해법의 정확성과 효율성을 예시하기 위하여, 2차원 및 3차원 크랙 문제의 변형 후 모습과 응력강도계수 계산 결과를 보였다.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.20 no.6
• /
• pp.553-563
• /
• 2008

If an arbitrary topography is approximated to a number of vertical steps, both variational approximation and eigenfunction expansion method can be used to compute linear wave transformation over the bottom. In this study a scatterer method associated with variational approximation is proposed to calculate reflection and transmission coefficients. Present method may be shown to be more simple and direct than the successive-application-matrix method by O'Hare and Davies. And Several numerical examples are given which are in good agreement with existing results.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.29 no.1A
• /
• pp.31-43
• /
• 2009

In this study, the numerical method is presented, which can consider the various train types and can solve the equations of motion for a vehicle-bridge interaction analysis by non-iteration procedure through formulating the coupled equations of motion. The coupled equations of motion for the vehicle-bridge interaction are solved by the Newmark ${\beta}$ of a direct integration method, and by composing the effective stiffness matrix and the effective force vector according to a analysis step, those can be solved with the same manner of the solving procedure of equilibrium equations in static analysis. Also, the effective stiffness matrix is reconstructed by the Skyline method for increasing the analysis effectiveness. The Cholesky's matrix decomposition scheme is applied to the analysis procedure for minimizing the numerical errors that can be generated in directly calculating the inverse matrix. The equations of motion for the conventional trains are derived, and the numerical models of the conventional trains are idealized by a set of linear springs and dashpots with 16 degrees of freedom. The bridge models are simplified by the 3 dimensional space frame element which is based on the Euler-Bernoulli theory. The rail irregularities of vertical and lateral directions are generated by the PSD functions of the Federal Railroad Administration (FRA). The results of the vehicle-bridge interaction analysis are verified by the experimental results for the railway plate girder bridges of a span length with 12 m, 18 m, and the experimental and analytical data are applied to the low pass filtering scheme, and the basis frequency of the filtering is a 2 times of the 1st fundamental frequency of a bridge bending.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.2 no.1
• /
• pp.1-9
• /
• 1978

냉각용 Fin의 설계문제를 일반적인 최적설계문제로 바꾸어서 일반화된 Steepest Descent 방법에 의한 수치적 방법을 도입하여 해결하였다. 보다 실제적인 문제를 다룰 수 있도록 여러가지 제한조건을 고려한 Fin의 최적곡선 모양의 해를 얻었으며 이 방법의 유용성을 보였다. 사다리꼴의 Fin 설계예에서 위 방법을 이용한 해와 직접 계산에 의한 열전달량의 등고선 그림으로부터 구한 해와 일치함을 보였다.

• Tunnel and Underground Space
• /
• v.26 no.4
• /
• pp.266-273
• /
• 2016

We have numerically modeled thermal evolution of Ulleung Island after an emplacement of magma chamber. The disk-shape magma chamber is assumed to locate at 2.9 km beneath the island and has a diameter and a thickness of 10 km and 300 (or 600) m, respectively. The geothermal gradients evaluated from the numerical modeling coincide well with the range of the geotherms (${\sim}95^{\circ}C/km$) observed from the well logging. Although there are limitations in the application of the numerical results directly to the interpretation of the observed geotherms, we believe that an existence of a hot magma chamber in molten or in solidified state is the most plausible explanation for the observed geotherms.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.14 no.1
• /
• pp.19-33
• /
• 2002

A numerical method to efficiently secure necessary design informations of the hydraulic characteristics of rubble mound breakwater was attempted here. The method combines the exterior wave field with the interior wave field which is formulated incorporating porous media flow inside the breakwaters. An approximate method based on the long wave assumption was used for the exterior wave field while a boundary element method was used for the interior wave field. A hydraulic experiment was also performed to verify the validity of the numerical analysis. The numerical results were compared with experimental data and results from existing formulae. They generally agreed in both reflection and transmission coefficients. The calculated pore pressures also showed a similar pattern with experimental data, even if they gave some significant differences in their values fur some cases. The main cause of such differences can be attributed to the strongly nonlinear wave field occurring on the frontal slope of the breakwater. The direct input of dynamic pressures(measured from hydraulic experiment) into the numerical method was suggested as a promising method to enhance the predictability of pore pressures.

• Journal of the Geological Society of Korea
• /
• v.54 no.6
• /
• pp.683-694
• /
• 2018

Subduction has been the focal point of numerical studies for decades because it plays an important role in the Earth's mass and energy circulations and generates earthquakes and arc volcanoes which are closely related to the human lives. Among the studies on subduction, numerical modeling has been broadly applied to the quantitative studies on the subducting slab in the mantle which cannot be directly observed. In this study, we benchmark the numerical subduction modeling using a finite element package, COMSOL $Multiphysics^{(R)}$ and the results are consistent with the previously reported benchmark results.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.29 no.6
• /
• pp.350-362
• /
• 2017

The Eulerian nearshore current model is more advantageous than the Lagrangian model in the way that numerical results from the Eulerian model can be directly compared with the measurements by the stationary equipment. It is because the wave mass flux is not included in the computed mass flux of Euleran nearshore current model. In addition, the Eulerian model can simulate the longshore currents with depth varying parabolic profile. However, the numerical models proposed by different researcher have different forms of the wave stress terms. For example, wave stresses in Newberger and Allen's (2007) model is constant over the depth, while those of Chun (2012) are vertically distributed. In the present study, these wave stress terms were compared against Hamilton et al.'s (2001) laboratory experiments to see the effects of different wave stress terms performed on the computation of nearshore currents.

• Economic and Environmental Geology
• /
• v.27 no.1
• /
• pp.41-47
• /
• 1994

A fast Hankel transform (FHT) algorithm (Anderson, 1982) is applied to numerical evaluation of many Green's tensor integrals encountered in three-dimensional electromagnetic modeling using integral equations. Efficient computation of Hankel transforms is obtained by a combination of related and lagged convolutions which are available in the FHT. We express Green's tensor integrals for a layered half-space, and rewrite those to a form of related functions so that the FHT can be applied in an efficient manner. By use of the FHT, a complete or full matrix of the related Hankel transform can be rapidly and accurately calculated for about the same computation time as would be required for a single direct convolution. Computing time for a five-layer half-space shows that the FHT is about 117 and 4 times faster than conventional direct and multiple lagged convolution methods, respectively.

• Journal of Energy Engineering
• /
• v.11 no.2
• /
• pp.130-135
• /
• 2002

Bypass air conditioning systems are divided into three types; outdoor air bypass, mixed air bypass and return air bypass system. Among bypass air conditioning systems, a return air bypass system is more effective than other two systems because it doesn't induce unconditioned outdoor air into conditioned room. The numerical study on the bypass air conditioning system shows this system can maintain indoor RH(Relative Humidity) less than a conventional CAV (Constant Air Volume) air conditioning system by adjusting face and bypass dampers at part load. A simulator was built to compare results of a numerical experiment and those of a simulator experiment. The results of the simulator experiment was nearly same as those of the numerical experiment; when a design sensible load (the ratio of sensible load to total sensible load) was 70 percent (at this time, RSHF=0.7), indoor relative humidity (in case of both numerical experiments and simulator experiments) was maintained below 60% specified by ASHRAE STANDARD 62-1999. The bypass air conditioning system is expected to be applied to many buildings where the Percentage of latent loads or air change tate is high.