• Journal of applied mathematics & informatics
• /
• v.14 no.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.

• Interaction and multiscale mechanics
• /
• v.2 no.4
• /
• pp.333-352
• /
• 2009

We introduce a radial basis collocation method to solve axially moving beam problems which involve $2^{nd}$ order differentiation in time and $4^{th}$ order differentiation in space. The discrete equation is constructed based on the strong form of the governing equation. The employment of multiquadrics radial basis function allows approximation of higher order derivatives in the strong form. Unlike the other approximation functions used in the meshfree methods, such as the moving least-squares approximation, $4^{th}$ order derivative of multiquadrics radial basis function is straightforward. We also show that the standard weighted boundary collocation approach for imposition of boundary conditions in static problems yields significant errors in the transient problems. This inaccuracy in dynamic problems can be corrected by a statically condensed semi-discrete equation resulting from an exact imposition of boundary conditions. The effectiveness of this approach is examined in the numerical examples.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.86-87
• /
• 2003

Unstructured grid system is suitable for flows of complex geometries. For problems with moving boundary walls, the grid system must be changed and deformed with time if we use a body fitted grid system. In this paper, a new moving-grid finite-volume method on unstructured grid system is proposed and developed for unsteady compressible flows with shock waves. To assure geometric conservation laws on moving grid system, a control volume on the space-time unified domain is adopted for estimating numerical flux. The method is described and applied for two-dimensional flows.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.28 no.6
• /
• pp.619-627
• /
• 2004

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method. The results agree well with previous numerical and experimental results. This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications. The usefulness of this method will be confirmed when we solve the complex geometries and moving bodies.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.25 no.5
• /
• pp.439-446
• /
• 2012

This paper presents an implicit moving least squares(MLS) difference method for improving the solution accuracy of 1-D free boundary problems, which implicitly updates the topology change of moving interface. The conventional MLS difference method explicitly updates the moving interface; it requires no iterative solution procedure but results in the loss of accuracy. However, the newly developed implicit scheme makes the total system nonlinear involving iterative solution procedure, but numerical verification show that it dramatically elevates the solution accuracy with moderate computation increase. Through numerical experiments for melting problems having moving singularity, it is verified that the proposed method can achieve the second order accuracy.

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.93-100
• /
• 2009

The gridless (or meshfree) methods, such as MPS, SPH, FPM an so forth, are feasible and robust for the problems with moving boundary and/or complicated boundary shapes, because these methods do not need to generate a grid system. In this study, a gridless solver, which is based on the combination of moving least square interpolations on a cloud of points with point collocation for evaluating the derivatives of governing equations, is presented for two-dimensional unsteady incompressible Navier-Stokes problem in the low Reynolds number. A MAC-type algorithm was adopted and the Poission equation for the pressure was solved successively in the moving least square sense. Some typical problems were solved by the presented solver for the validation and the results obtained were compared with analytic solutions and the numerical results by conventional CFD methods, such as a FVM.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.4
• /
• pp.315-322
• /
• 2009

This paper presents a novel numerical method based on the extended moving least squares finite difference method(MLS FDM) for solving 1-D Stefan problem. The MLS FDM is employed for easy numerical modelling of the moving boundary and Taylor polynomial is extended using wedge function for accurate capturing of interfacial singularity. Difference equations for the governing equations are constructed by implicit method which makes the numerical method stable. Numerical experiments prove that the extended MLS FDM show high accuracy and efficiency in solving semi-infinite melting, cylindrical solidification problems with moving interfacial boundary.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.1
• /
• pp.39-48
• /
• 2013

This paper develops a 2-D moving least squares(MLS) difference method for Stefan problem by extending the 1-D version of the conventional method. Unlike to 1-D interfacial modeling, the complex topology change in 2-D domain due to arbitrarily moving boundary is successfully modelled. The MLS derivative approximation that drives the kinetics of moving boundary is derived while the strong merit of MLS Difference Method that utilizes only nodal computation is effectively conserved. The governing equations are differentiated by an implicit scheme for achieving numerical stability and the moving boundary is updated by an explicit scheme for maximizing numerical efficiency. Numerical experiments prove that the MLS Difference Method shows very good accuracy and efficiency in solving complex 2-D Stefan problems.

• Proceedings of the Korean Society of Soil and Groundwater Environment Conference
• /
• 2001.04a
• /
• 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.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.4 no.3
• /
• pp.146-155
• /
• 1992

It has been conventional to treat the land boundary as a fixed one in numerical modeling of tidal flows, particularly in the finite element scheme. However conventional models using the fixed land boundary result in unrealistic tidal flows in inter-tidal zones which exist over wide coastal area in Korea. In this study, a 2-dimensional hydrodynamic model, using finite element method for moving boundary problems was developed. The performance of the model was tested in a rectangular channel with an open boundary at one end and a moving boundary at the other end. The model was applied to calculate the tidal currents in Maro Hae, located in the southwestern part of Korea where wide tidal flats develop. The behavior of tidal currents in the Udolmok and near the tidal flats in the study area was satisfactory when compared with the observed data. Variation of tidal currents due to the construction of Kochunam sea-dyke which barrages large area of tidal flat was presented. The results of this study confirm the efficiency of moving boundary treatment in coastal numerical models.