• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.189-194
• /
• 2004

As a disk drive becomes widely used in portable environments, one of the important requirements is durability under severe environmental condition, especially, resistance to mechanical shock. An important challenge in the disk recording is to improve disk drive robustness in shock environments. If the system comes In contact with outer shock disturbance, the system gets critical damage in head-gimbal assembly or disk. This paper describes analysis of a HGA(head-gimbal assembly) in micro MO drives to shock loading during both non-operating state and operating state. A finite element model which consists of the disk, suspension, slider and air bearing was used to find structural response of micro MO drives. In the operational case. the air bearing is approximated with four linear elastic springs. The commercially available finite element solver, ANSYS/LS-DYNA, is used to simulate the shock response of the HGA in micro MO drives. In this paper, the mechanical robustness of the suspension is simuiated considering the shock responses of the HGA.

• Journal of ILASS-Korea
• /
• v.23 no.1
• /
• pp.16-21
• /
• 2018

The current study has developed an in-house 3D FEM code in order to model thermoacoustic problems in a gas turbine combustion system and compared calculation results of main instability characteristics with measured ones from a lab-scale partially-premixed combustor. From the comparison of calculation results with the measured data, the current model could successfully capture the harmonic longitudinal instability frequencies and their spatial distributions of the acoustic field as well as the growth rate of self-excited modes.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.29 no.1
• /
• pp.61-65
• /
• 2016

Composites are materials that are widely used in industries such as automobile and aircraft. The composite material is required as a material for using in a high temperature environment as well as acting as a high pressure environment like the nozzle part of the ship. It is important to know the properties of composites. Result values obtained substituting the properties of matrix and fiber numerically have an large error compared with experimental value. In this study we utilize CASADsolver EDISON program for using Finite Element Method. Properties by substituting the fiber and Matrix properties of the composite material properties were compared with those measured in the experiment and calculated by the empirical properties.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1997.10a
• /
• pp.17-23
• /
• 1997

In the two-dimensional Finite Element Method for forming simulation, mesh generation and remeshing process are very significant. In this paper, using the modified splitting mesh generation algorithm, we can overcome the limitation of existing techniques and acquire mesh, which has optimal mesh density. A modified splitting algorithm for automatically generating quadrilateral mesh within a complex domain is described. Unnecessary meshing process for density representation is removed. Especially, during the mesh generation with high gradient density like as shear band representation, the modified mesh density scheme, which will generate quadrilateral mesh with the minimized error, which takes effect on FEM solver, is introduced.

• Journal of the Korean Mathematical Society
• /
• v.50 no.1
• /
• pp.81-93
• /
• 2013

We consider a distributed optimal flux control problem: finding the potential of which gradient approximates the target vector field under an elliptic constraint. Introducing the Lagrange multiplier and a change of variables the Euler-Lagrange equation turns into a coupled equation of an elliptic equation and a reaction diffusion equation. The change of variables reduces iteration steps dramatically when the Gauss-Seidel iteration is considered as a solution method. For the elliptic equation solver we consider the Cell Boundary Element (CBE) method, which is the finite element type flux preserving methods.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2008.11a
• /
• pp.561-562
• /
• 2008

The numerical analysis of sound radiation by vibrating structure is a well known and mature technology used in many industries. Accurate methods based on the boundary or finite element method have been successfully developed over the last two decades and are now available in standard CAE tools. These methods are however known to require significant computational resources which, furthermore, very quickly increase with the frequency of interest. The low speed of most current methods is a main obstacle for a systematic use of acoustic CAE in industrial design processes. In this paper we are going to present a set of innovative techniques that significantly speed-up the calculation of acoustic radiation indicators (acoustic pressure, velocity, intensity and power; contribution vectors). The modeling is based on the well known combination of finite elements and infinite elements but also combines the following ingredients to obtain a very high performance: o a multi-frontal massively parallel sparse direct solver; o a multi-frequency solver based on the Krylov method; o the use of pellicular acoustic modes as a vector basis for representing acoustic excitations; o the numerical evaluation of Green functions related to the specific geometry of the problem under investigation. All these ingredients are embedded in the ACTRAN/AR CAE tool which provides unprecedented performance for acoustic radiation analysis. The method will be demonstrated on several applications taken from various industries.

• 제어로봇시스템학회:학술대회논문집
• /
• 2003.10a
• /
• pp.1349-1353
• /
• 2003

CEMTool is a command style design and analyzing package for scientific and technological algorithm and a matrix based computation language. In this paper, we present a compiler based approach to the implementation of the command mode generalized PDE solver in CEMTool. In contrast to the existing MATLAB PDE Toolbox, our proposed FEM package can deal with the combination of the reserved words such as "laplace" and "convect". Also, we can assign the border lines and the boundary conditions in a very easy way. With the introduction of the lexical analyzer and the parser, our FEM toolbox can handle the general boundary condition and the various PDEs represented by the combination of equations. That is why we need not classify PDE as elliptic, hyperbolic, parabolic equations. Consequently, with our new FEM toolbox, we can overcome some disadvantages of the existing MATLAB PDE Toolbox.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.22 no.2
• /
• pp.101-113
• /
• 2018

The dual singular function method(DSFM) is a numerical algorithm to get optimal solution including corner singularities for Poisson and Helmholtz equations. In this paper, we apply DSFM to solve heat equation which is a time dependent problem. Since the DSFM for heat equation is based on DSFM for Helmholtz equation, it also need to use Sherman-Morrison formula. This formula requires linear solver n + 1 times for elliptic problems on a domain including n reentrant corners. However, the DSFM for heat equation needs to pay only linear solver once per each time iteration to standard numerical method and perform optimal numerical accuracy for corner singularity problems. Because the Sherman-Morrison formula is rather complicated to apply computation, we introduce a simplified formula by reanalyzing the Sherman-Morrison method.

• Journal of computational fluids engineering
• /
• v.10 no.2
• /
• pp.38-47
• /
• 2005

Substructuring methods are often used in finite element structural analyses. In this study a multi-level substructuring(MLSS) algorithm is developed and proposed as a possible candidate for finite element fluid solvers. The present algorithm consists of four stages such as a gathering, a condensing, a solving and a scattering stage. At each level, a predetermined number of elements are gathered and condensed to form an element of higher level. At the highest level, each sub-domain consists of only one super-element. Thus, the inversion process of a stiffness matrix associated with internal degrees of freedom of each sub-domain has been replaced by a sequential static condensation of gathered element matrices. The global algebraic system arising from the assembly of each sub-domain matrices is solved using a well-known iterative solver such as the conjugare gradient(CG) or the conjugate gradient squared(CGS) method. A time comparison with CG has been performed on a 2-D Poisson problem. With one domain the computing time by MLSS is comparable with that by CG up to about 260,000 d.o.f. For 263,169 d.o.f using 8 x 8 sub-domains, the time by MLSS is reduced to a value less than $30\%$ of that by CG. The lid-driven cavity problem has been solved for Re = 3200 using the element interpolation degree(Deg.) up to cubic. in this case, preconditioning techniques usually accompanied by iterative solvers are not needed. Finite element formulation for the incompressible flow has been stabilized by a modified residual procedure proposed by Ilinca et al.[9].

• Journal of the Korean Mathematical Society
• /
• v.59 no.3
• /
• pp.519-548
• /
• 2022

In this paper, a stabilized-penalized collocated finite volume (SPCFV) scheme is developed and studied for the stationary generalized Navier-Stokes equations with mixed Dirichlet-traction boundary conditions modelling an incompressible biological fluid flow. This method is based on the lowest order approximation (piecewise constants) for both velocity and pressure unknowns. The stabilization-penalization is performed by adding discrete pressure terms to the approximate formulation. These simultaneously involve discrete jump pressures through the interior volume-boundaries and discrete pressures of volumes on the domain boundary. Stability, existence and uniqueness of discrete solutions are established. Moreover, a convergence analysis of the nonlinear solver is also provided. Numerical results from model tests are performed to demonstrate the stability, optimal convergence in the usual L² and discrete H¹ norms as well as robustness of the proposed scheme with respect to the choice of the given traction vector.