• Korean Journal of Air-Conditioning and Refrigeration Engineering
• /
• v.28 no.7
• /
• pp.293-298
• /
• 2016

A three-dimensional (3D) numerical model of the vertical ground-coupled heat exchanger is useful for analyzing the modern ground source heat pump system. Furthermore, a detailed description of the inner side of the exchanger allows to account for the effects of the thermal capacity. Thus, both methods are included in the proposed numerical model. For the ground portion, a FDM (Finite Difference Method) scheme has been applied using the Cartesian coordinate system. Cylindrical grids are applied for the borehole portion, and the U-tube configuration is adjusted at the grid, keeping the area and distance unchanged. Two sub-models are numerically coupled at each time-step using an iterative method for convergence. The model is validated by a reference 3D model under a continuous heat injection case. The results from a periodic heat injection input show that the proposed thermal capacity model reacts more slowly to the changes, resulting in lower borehole wall temperatures, when compared with a thermal resistance model. This implies that thermal capacity effects may be important factors for system controls.

• Structural Engineering and Mechanics
• /
• v.64 no.5
• /
• pp.621-640
• /
• 2017

This paper reports the numerical investigation conducted to study the influence of Local-Distortional (L-D) interaction mode buckling on post buckling strength erosion in fixed ended lipped channel cold formed steel columns. This investigation comprises of 81 column sections with various geometries and yield stresses that are carefully chosen to cover wide range of strength related parametric ratios like (i) distortional to local critical buckling stress ratio ($0.91{\leq}F_{CRD}/F_{CRL}{\leq}4.05$) (ii) non dimensional local slenderness ratio ($0.88{\leq}{\lambda}_L{\leq}3.54$) (iii) non-dimensional distortional slenderness ratio ($0.68{\leq}{\lambda}_D{\leq}3.23$) and (iv) yield to non-critical buckling stress ratio (0.45 to 10.4). The numerical investigation is carried out by conducting linear and non-linear shell finite element analysis (SFEA) using ABAQUS software. The non-linear SFEA includes both geometry and material non-linearity. The numerical results obtained are deeply analysed to understand the post buckling mechanics, failure modes and ultimate strength that are influenced by L-D interaction with respect to strength related parametric ratios. The ultimate strength data obtained from numerical analysis are compared with (i) the experimental tests data concerning L-D interaction mode buckling reported by other researchers (ii) column strength predicted by Direct Strength Method (DSM) column strength curves for local and distortional buckling specified in AISI S-100 (iii) strength predicted by available DSM based approaches that includes L-D interaction mode failure. The role of flange width to web depth ratio on post buckling strength erosion is reported. Then the paper concludes with merits and limitations of codified DSM and available DSM based approaches on accurate failure strength prediction.

• Proceedings of the KIEE Conference
• /
• 2005.07d
• /
• pp.2897-2899
• /
• 2005

Finite element method (FEM) has been widely used as a useful numerical method that can analyze complex engineering problems in electro-magnetics, mechanics, and others. CEMTool, which is similar to MATLAR, is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present new 3D FEM package in CEMTool environment. In contrast to the existing CEMTool 2D FEM package and MATLAB PDE (Partial Differential Equation) Toolbox, our proposed 3D FEM package can deal with complex 3D models, not a cross-section of 3D models. Consequently, with our new 3D FEM toolbox, we can analyze more diverse engineering Problems which the existing CEMTool 2D FEM package or MATLAB PDE Toolbox can not solve.

• Structural Engineering and Mechanics
• /
• v.88 no.5
• /
• pp.439-449
• /
• 2023

Non-linear vibration characteristics of functionally graded CNT-reinforced composite (FG-CNTRC) cylindrical shell panel on elastic foundation have not been sufficiently examined. In this situation, this study aims at the profound numerical investigation of the non-linear vibration response of FG-CNTRC cylindrical panels on Winkler-Pasternak foundation by introducing an accurate and effective 2-D meshfree-based non-linear numerical method. The large-amplitude free vibration problem is formulated according to the first-order shear deformation theory (FSDT) with the von Karman non-linearity, and it is approximated by Laplace interpolation functions in 2-D natural element method (NEM) and a non-linear partial derivative operator H_NL. The complex and painstaking numerical derivation on the curved surface and the crucial shear locking are overcome by adopting the geometry transformation and the MITC3+ shell elements. The derived nonlinear modal equations are iteratively solved by introducing a three-step iterative solving technique which is combined with Lanczos transformation and Jacobi iteration. The developed non-linear numerical method is estimated through the benchmark test, and the effects of foundation stiffness, CNT volume fraction and functionally graded pattern, panel dimensions and boundary condition on the non-linear vibration of FG-CNTRC cylindrical panels on elastic foundation are parametrically investigated.

• Geomechanics and Engineering
• /
• v.18 no.6
• /
• pp.651-659
• /
• 2019

The IB-SEM numerical method combines the spectral/hp element method and the rigid immersed boundary method. This method avoids the problems of low computational efficiency and errors that are caused by the re-division of the grid when the solids move. Based on the Fourier transformation and the 3D immersed boundary method, the 3D IB-SEM system was established. Then, using the open MPI and the Hamilton HPC service, the computational efficiency was increased substantially. The flows around a cylinder and a sphere were simulated by the system. The surface of the cylinder generates vortices with alternating shedding, and these vortices result in a periodic force acting on the surface of the cylinder. When the shedding vortices enter the flow field behind the cylinder, a recirculation zone is formed. Finally, the three-dimensional pore flow was successfully investigated.

• Proceedings of the KIEE Conference
• /
• 2002.11d
• /
• pp.97-99
• /
• 2002

3D Equivalent magnetic Circuit Network Method (EMCNM) is comparatively the easy way that analyzes 3D models of Electric Machine by using permeance as a distributive magnetic circuit parameter under the existing magnetic equivalent circuit method and Numerical Method. The existing 3D EMCNM could not correctly describe the shape of an analysis target when using rectangular shape element or fan shape element, so it made errors when calculating permeance. Therefore, this paper use the trapezoidal element contained rectangular element. This paper calculate the torque based on the maxwell stress tenser method when the airgap is one layer, three layers, and 5 layers, respectively.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Transactions of Materials Processing
• /
• v.10 no.1
• /
• pp.23-29
• /
• 2001

In this study, the optimized initial hole shape for T-branch forming was proposed to obtain effective welding region. Design variables were determined by approximation analysis using volume constant condition. We performed 3D elastic-plastic FEM(Finite Element Method) analysis to simulate T-branch forming process. The variation of height and thickness of T-branch with various hole shapes was investigated. The optimized initial hole shape equation was obtained by using results for the numerical analysis.

• Geophysics and Geophysical Exploration
• /
• v.9 no.3
• /
• pp.241-249
• /
• 2006

Due to the long tectonic history and the very complex geologic formations in Korea, the anisotropic characteristics of subsurface material may often change very greatly and locally. The algorithms commonly used, however, may not give sufficiently precise computational results of traveltime data particularly for the complex and strong anisotropic model, since they are based on the two-dimensional (2D) earth and/or weak anisotropy assumptions. This study is intended to develope a three-dimensional (3D) modeling algorithm to precisely calculate the first arrival time in the complex anisotropic media. Considering the complex geology of Korea, we assume 3D TTI (tilted transversely isotropy) medium having the arbitrary symmetry axis. The algorithm includes the 2D non-linear interpolation scheme to calculate the traveltimes inside the grid and the 3D traveltime mapping to fill the 3D model with first arrival times. The weak anisotropy assumption, moreover, can be overcome through devising a numerical approach of the steepest descent method in the calculation of minimum traveltime, instead of using approximate solution. The performance of the algorithm developed in this study is demonstrated by the comparison of the analytic and numerical solutions for the homogeneous anisotropic earth as well as through the numerical experiment for the two layer model whose anisotropic properties are greatly different each other. We expect that the developed modeling algorithm can be used in the development of processing and inversion schemes of seismic data acquired in strongly anisotropic environment, such as migration, velocity analysis, cross-well tomography and so on.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1057-1064
• /
• 2005

A conservative finite-volume numerical method using unstructured meshes, which is developed by the authors, is evaluated for its application to several 2-D benchmark problems using a variety of quadrilateral, triangular and hybrid meshes. The present pressure-based numerical method for unstructured mesh clearly demonstrates the same accuracy and robustness as that fur typical structured mesh.