• Journal of applied mathematics & informatics
• /
• 제14권1_2호
• /
• pp.97-113
• /
• 2004

In this paper we shall study moving boundary problems, and we introduce an approach for solving a wide range of them by using calculus of variations and optimization. First, we transform the problem equivalently into an optimal control problem by defining an objective function and artificial control functions. By using measure theory, the new problem is modified into one consisting of the minimization of a linear functional over a set of Radon measures; then we obtain an optimal measure which is then approximated by a finite combination of atomic measures and the problem converted to an infinite-dimensional linear programming. We approximate the infinite linear programming to a finite-dimensional linear programming. Then by using the solution of the latter problem we obtain an approximate solution for moving boundary function on specific time. Furthermore, we show the path of moving boundary from initial state to final state.

• 한국지하수토양환경학회:학술대회논문집
• /
• 한국지하수토양환경학회 2001년도 총회 및 춘계학술발표회
• /
• pp.77-80
• /
• 2001

The element-free Galerkin (EFG) method is one of meshless methods, which is an efficient method of modeling problems of fluid or solid mechanics with complex boundary shapes and large changes in boundary conditions. This paper discusses the theory of the EFG method and its applications to modeling of groundwater flow. In the EFG method, shape functions are constructed based on the moving least square (MLS) approximation, which requires only set of nodes. The EFG method can eliminate time-consuming mesh generation procedure with irregular shaped boundaries because it does not require any elements. The coupled EFG-FEM technique was introduced to treat Dirichlet boundary conditions. A computer code EFGG was developed and tested for the problems of steady-state and transient groundwater flow in homogeneous or heterogeneous aquifers. The accuracy of solutions by the EFG method was similar to that by the FEM. The EFG method has the advantages in convenient node generation and flexible boundary condition implementation.

• Structural Engineering and Mechanics
• /
• 제21권3호
• /
• pp.315-332
• /
• 2005

The Element-free Galerkin Method has become a very popular tool for the simulation of mechanical problems with moving boundaries. The internally applied Moving Least Squares interpolation uses in general Gaussian or cubic weighting functions and has compact support. Due to the approximative character of this interpolation the obtained shape functions do not fulfill the interpolation conditions, which causes additional numerical effort for the application of the boundary conditions. In this paper a new weighting function is presented, which was designed for meshless shape functions to fulfill these essential conditions with very high accuracy without any additional effort. Furthermore this interpolation gives much more stable results for varying size of the influence radius and for strongly distorted nodal arrangements than existing weighting function types.

• Ocean Systems Engineering
• /
• 제6권3호
• /
• pp.245-264
• /
• 2016

The iterative boundary element method (IBEM) developed originally before for cavitating two-dimensional (2-D) and three-dimensional (3-D) hydrofoils moving under free surface is modified and applied to the case of 2-D (two-dimensional) airfoils and 3-D (three-dimensional) wings over water. The calculation of the steady-state flow characteristics of an inviscid, incompressible fluid past 2-D airfoils and 3-D wings above free water surface is of practical importance for air-assisted marine vehicles such as some racing boats including catamarans with hydrofoils and WIG (Wing-In-Ground) effect crafts. In the present paper, the effects of free surface both on 2-D airfoils and 3-D wings moving steadily over free water surface are investigated in detail. The iterative numerical method (IBEM) based on the Green's theorem allows separating the airfoil or wing problems and the free surface problem. Both the 2-D airfoil surface (or 3-D wing surface) and the free surface are modeled with constant strength dipole and constant strength source panels. While the kinematic boundary condition is applied on the airfoil surface or on the wing surface, the linearized kinematic-dynamic combined condition is applied on the free surface. The source strengths on the free surface are expressed in terms of perturbation potential by applying the linearized free surface conditions. No radiation condition is enforced for downstream boundary in 2-D airfoil and 3-D wing cases and transverse boundaries in only 3-D wing case. The method is first applied to 2-D NACA0004 airfoil with angle of attack of four degrees to validate the method. The effects of height of 2-D airfoil from free surface and Froude number on lift and drag coefficients are investigated. The method is also applied to NACA0015 airfoil for another validation with experiments in case of ground effect. The lift coefficient with different clearance values are compared with those of experiments. The numerical method is then applied to NACA0012 airfoil with the angle of attack of five degrees and the effects of Froude number and clearance on the lift and drag coefficients are discussed. The method is lastly applied to a rectangular 3-D wing and the effects of Froude number on wing performance have been investigated. The numerical results for wing moving under free surface have also been compared with those of the same wing moving above free surface. It has been found that the free surface can affect the wing performance significantly.

• 한국해양공학회지
• /
• 제21권6호
• /
• pp.53-58
• /
• 2007

A particle method, recognized as one of gridless methods, has been developed to investigate non-linear free-surface motions interacting with structures. This method is more feasible and effective than conventional grid methods for solving flow fieldswith complicated boundary shapes. The method consists of particle interaction models representing pressure gradient, diffusion, incompressibility, and the free-surface boundary conditions without grids. In the present study, broken dam problems with various viscosity values are simulated to validate the developed method.

• Interaction and multiscale mechanics
• /
• 제6권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.

• Structural Engineering and Mechanics
• /
• 제25권4호
• /
• pp.467-480
• /
• 2007

Two methods were developed for analyzing problems with two adjacent sub-domains modeled by different kinds of elements in finite element method. Each sub-domain can be defined independently without the consideration of equivalent division with common nodes used for the interface. These two methods employ an individual interface to accomplish the compatibility. The MLSA method uses the moving least square approximation which is the basic formulation for Element Free Galerkin Method to formulate the interface. The displacement field assumed by this method does not pass through nodes on the common boundary. Therefore, nodes can be chosen freely for this method. The results show that the MLSA method has better approximation than traditional methods.

• 한국해양공학회:학술대회논문집
• /
• 한국해양공학회 2006년 창립20주년기념 정기학술대회 및 국제워크샵
• /
• pp.382-385
• /
• 2006

A high-order spectral/boundary-element method is newly adapted as an efficient numerical tool. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved fly using the high-order spectral method and body potential is solved fly using the high-order boundary element method. Through the combination of these two methods, the wave-making problems fly a submerged sphere moving with the large amplitude oscillation are solved in time-domain. With the example calculations, nonlinear effects on free-surface profiles and hydrodynamic forces are shown and discussed.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 춘계학술대회논문집
• /
• pp.165-169
• /
• 2009

The effect of the Dirichlet boundary condition for the redistance equation of level set method on the solutionof sloshing problem is investigated by adopting four Dirichlet boundary conditions. For the solution of the incompressible Navier-Stokes equations, P1P1 four-step fractional finite element method is employed and a least-square finite element method is used for the solutions of the two hyperbolic type equations of level set method; advection and redistance equation. ALE (Arbitrary Lagrangian Eulerian) method is used to deal with a moving computational domain. It has been shown that the free surface motion in a sloshing tank is strongly dependent on the type of the Dirichlet boundary condition and the results of broken dam and sloshing problems using various Dirichlet boundary conditions are discussed and compared with the existing experimental results.

• KIEE International Transaction on Electrical Machinery and Energy Conversion Systems
• /
• 제4B권4호
• /
• pp.180-183
• /
• 2004

This paper presents a description of the point collocation method and its application to the electromagnetic field computation. The interpolation scheme is based on the fast moving least square reproducing kernel approximation. In the method, the integration cell is not required and the essential boundary conditions can be enforced directly. Numerical simulations on 1-D and 2-D problems are carried out to validate the method. It is found that computational efficiency is higher than the general mesh-free methods.