• Honam Mathematical Journal
• /
• v.38 no.1
• /
• pp.47-57
• /
• 2016

In the discrete formulation of the bubble stabilized Legendre Galerkin methods, the system of equations includes the artificial viscosity term as the parameter. We investigate the estimation of this parameter to get the least-squares solution which minimizes the sum of the squares of errors at each node points. Some numerical results are reported.

• Journal of computational fluids engineering
• /
• v.13 no.3
• /
• pp.8-13
• /
• 2008

We present a numerical simulation technique and some preliminary results of the impact and spreading of a droplet containing particles on the solid substrate in 2D. We used the 2nd-order Adams-Bashforth / Crank-Nicholson method to solve the Navier-Stokes equation and employed the level-set method with the continuous surface stress for description of droplet spreading with interfacial tension. The impact velocity has been generated by the instantaneous gravity. The distributed Lagrangian-multipliers method has been combined for the implicit treatment of rigid particles and the discontinuous Galerkin method has been used for the stabilization of the interface advection equation. We investigated the droplet spreading by the inertial force and discussed effects of the presence of particles on the spreading behavior using an example problem. We observed reduced oscillation and spread for the particulate droplet.

• Proceedings of the Korean Institute of Information and Commucation Sciences Conference
• /
• 2000.05a
• /
• pp.143-146
• /
• 2000

This paper presents a study on design for the k-junction of a feeder waveguide with an inductive wall using FDTD method. The feed structure consists of a single waveguide plated on the same layer as radiating waveguide and is characterized by the unit divider, railed a $\pi$-Junction. This $\pi$-Junction with an inductive wall splits part of the power into two branch waveguide through one coupling window, and can excite densely arrayed waveguide at equal phase and amplitude. The power dividing characteristics of a $\pi$ -Junction obtained by FDTD method are compared with one of Galerkin's method of moments. The scattering matrices a $\pi$ -junction calculated by FDTD method are realized.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2001.04a
• /
• pp.335-342
• /
• 2001

In this paper, an improved Element-Free Galerkin (EFG) method is proposed by adding enrichment function to the standard EFG approximation and a discontinuity function is implemented in constructing the shape function across the crack surface. In this method, the singularity and the discontinuity of the crack are efficiently modeled by using initial node distribution to evaluate reliable stress intensity factor, though the standard EFG method requires placing additional nodes near the crack tip. The proposed method enables the initial node distribution to be kept without any additional nodal d.o.f. and expresses the asymptotic stress field near the crack tip successfully. Numerical example verifies the improvement and the effectiveness of the method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.375-382
• /
• 2003

Porous media consist of physically and chemically different materials and have an extremely complicated behavior due to the different material properties of each of its constituents. In addition, the internal structure of porous media has generally a complex geometry that makes the description of its mechanical behavior quite complex. Thus, in order to describe and clarify the deformation behavior of porous media, constitutive models for deformation of porous media coupling several effects such as flow of fluids of thermodynamical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian methods, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of solids and fluids. First of all, governing equations for saturated porous media based on ALE description are derived. Then, weak forms of these equations are obtained in order to implement numerical method using finite element method. Finally, Petrov-Galerkin method Is applied to develop finite element formulation.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.257-270
• /
• 2013

The use of lightweight materials has been steadily increasing in the automotive industry, and presents new challenges to material joining. Among many joining processes, self-piercing riveting (SPR) is particularly promising for joining lightweight materials (such as aluminum alloys) and dissimilar materials (such as steel to Al, and metal to polymer). However, to establish a process window for optimal joint performance, it often requires a long trial-and-error testing of the SPR process. This is because current state of the art in numerical analysis still cannot effectively resolve the problems of severe material distortion and separation in the SPR simulation. This paper presents a coupled meshfree/finite element with a moving boundary algorithm to overcome these numerical difficulties. The simulation results are compared with physical measurements to demonstrate the effectiveness of the present method.

• Journal of the Korean Institute of Telematics and Electronics D
• /
• v.34D no.8
• /
• pp.41-51
• /
• 1997

The device simulator (BANDIS) which can analyze efficiently the electrical characteristics of the semiconductor devices under the three dimensional stationary conditions on the IBM PC was developed. Poisson, electon and hole continuity equations are discretized y te galerkin method using a tetrahedron as af finite element. The frontal solver which has exquisite data structures and advanced input/output functions is dused for the matrix solver which needs the highest cost in the three dimensional device simulation. The discretization method of the continuity equations used in BANDIS are compared with that of the scharfetter-gummel method used in the commercial three-dimensional device. To verify an accuracy and the efficiency of the discretization method, the simulation results of the PN junction diode and the BJT from BANDIS are compared with those of the commercial three-dimensiional device simulator such as DAVINCI. The maximum relative error within 2% and the average number of iterations needed for the convergence is decreased by more than 20%. The total simulation time of the BJT with 25542 nodes is decreased to about 60% compared with that of DAVINCI.

• Structural Engineering and Mechanics
• /
• v.75 no.3
• /
• pp.311-325
• /
• 2020

In applying the spectral stochastic finite element methods to the frequency response analysis, the conventional methods are known to give unstable and inaccurate results near the natural frequencies. To address this issue, a new sensitivity based stabilized formulation for stochastic frequency response analysis is proposed in this paper. The main difference over the conventional spectral methods is that the polynomials of random variables are applied to both numerator and denominator in approximating the harmonic response solution. In order to reflect the resonance behavior of the structure, the denominator polynomials is constructed by utilizing the natural frequency sensitivity and the random mode superposition. The numerator is approximated by applying a polynomial chaos expansion, and its coefficients are obtained through the Galerkin or the spectral projection method. Through various numerical studies, it is seen that the proposed method improves accuracy, especially in the vicinities of structural natural frequencies compared to conventional spectral methods.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.11
• /
• pp.2312-2321
• /
• 2002

The first-order least-squares meshfree method for linear elasticity is presented. The conventional and the compatibility-imposed least-squares formulations are studied on the convergence behavior of the solution and the robustness to integration error. Since the least-squares formulation is a type of mixed formulation and induces positive-definite system matrix, by using shape functions of same order for both primal and dual variables, higher rate of convergence is obtained for dual variables than Galerkin formulation. Numerical examples also show that the presented formulations do not exhibit any volumetric locking for the incompressible materials.