• Kyungpook Mathematical Journal
• /
• 제60권3호
• /
• pp.629-645
• /
• 2020

A uniformly convergent numerical method is developed for solving singularly perturbed 1-D parabolic convection-diffusion problems. The developed method applies a non-standard finite difference method for the spatial derivative discretization and uses the implicit Runge-Kutta method for the semi-discrete scheme. The convergence of the method is analyzed, and it is shown to be first order convergent. To validate the applicability of the proposed method two model examples are considered and solved for different perturbation parameters and mesh sizes. The numerical and experimental results agree well with the theoretical findings.

• 설비공학논문집
• /
• 제28권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
• /
• 제64권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.

• Structural Engineering and Mechanics
• /
• 제88권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.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2005년도 제36회 하계학술대회 논문집 D
• /
• 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.

• Geomechanics and Engineering
• /
• 제18권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.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 2002년도 추계학술대회 논문집 전기기기 및 에너지변환시스템부문
• /
• 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
• /
• 제4권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.

• 지구물리와물리탐사
• /
• 제9권3호
• /
• pp.241-249
• /
• 2006

한국의 지질 환경은 암석 분포가 매우 다양하고 복잡한 구조 활동의 영향을 받아 지하매질의 이방성 특성이 국부적으로 심하게 변화한다. 기존의 이방성 주시 모델링의 경우 지질 모델을 2차원으로 단순화시킴으로써 이러한 복잡한 지질 환경을 제대로 고려할 수 없었다. 또한 약 이방성 가정을 사용하여 실제로 나타날 수 있는 지하 매질의 심각한 이방성 영향을 주시 모델링에서 고려할 수 없었다. 이에 이 연구에서는 보다 실제적이고 복잡한 3차원 횡등방성 매질(transversely isotropic media)에서 q-P파의 초동 주시 양상을 모사할 수 있는 주시 모델링 알고리듬을 개발하였다. 이 알고리듬에서는 2차원 비선형 주시 내삽(2D nonlinear traveltime interpolation) 기법과 주시의 3차원 격자 채움법(mapping)을 이용한 직접 전파법(direct calculation)을 통해 급격한 물성의 변화에도 주시 계산이 가능하도록 하였다. 또한, 최소 주시 계산과정에서 수치 미분을 통한 최대 경사법(steepest descent method)을 사용하여 약 이방성 가정을 극복하였다. 개발된 알고리듬은 해석해와 비교하여 그 타당성을 검증하였고 3차원 2층구조에 대한 주시 계산을 수행하여 물성이 급격히 변화하는 모델에 대해서도 안정적으로 주시 계산이 이루어짐을 확인하였다. 이 연구에서 개발한 3차원 주시 모델링 알고리듬은 향후 구조보정이나 토모그래피 알고리듬 개발에 사용될 수 있을 것으로 기대한다.

• 소성∙가공
• /
• 제10권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.