• 한국해양공학회지
• /
• 제10권3호
• /
• pp.62-73
• /
• 1996

Numerical simulation is made of the three-dimensional wave breaking motion about a part of a floating offshore structure containing a circular cylinder mounted vertically onto a lower hull in regular periodic gravity wave generated by a numerical wave maker. TUMMAC-VIII finite-difference method is newly developed for such a problem. By use of density-function technique the three-dimensional wave breaking motion is approximately implenented in the framework of rectangular grid system. A porosity technique is devised for the implementation of the no-slip bydy boundary conditions. The generation of breaking waves by the interaction of incident waves with the structure is well simulated and interesting features of breaking waves are revealed with containing degree of quantitative and qualitative accuracy.

• 대한수학회지
• /
• 제57권5호
• /
• pp.1239-1266
• /
• 2020

In this paper, some novel discrete formulations for stabilizing the mixed finite element method Q1-Q0 (bilinear velocity and constant pressure approximations) are introduced and discussed for the generalized Stokes problem. These are based on stabilizing discontinuous pressure approximations via local jump operators. The developing idea consists in a reduction of terms in the local jump formulation, introduced earlier, in such a way that stability and convergence properties are preserved. The computer implementation aspects and numerical evaluation of these stabilized discrete formulations are also considered. For illustrating the numerical performance of the proposed approaches and comparing the three versions of the local jump methods alongside with the global jump setting, some obtained results for two test generalized Stokes problems are presented. Numerical tests confirm the stability and accuracy characteristics of the resulting approximations.

• Selected Papers of The Society of Naval Architects of Korea
• /
• 제2권1호
• /
• pp.63-78
• /
• 1994

This paper deals with the treatment of the open boundary in two-dimensional free-surface wave problems. Two numerical schemes are investigated for the implementation of the open boundary condition. One is to add the artificial damping term to the dynamic free-surface boundary condition, in which the determination of suitable damping coefficient and the damping zone is the most important. The other is a modified Orlanski's method, which is known to be very useful for the uni-directional waves. Using these two schemes, numerical tests have been conducted for a few typical free-surface wave problems. To obtain the numerical solution of the free-surface boundary value problem, the fundamental source-distribution method is used and the fully nonlinear free-surface boundary conditions are applied. The computed results are presented in comparison with those of others for the proof of practicality of these two schemes.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• 제22권2호
• /
• pp.125-136
• /
• 2018

In this article, a homotopy perturbation transform method (HPTM) and the Laplace transform combined with Taylor expansion method are presented for solving Volterra integral equations with a convolution kernel. The (HPTM) is innovative in Laplace transform algorithm and makes the calculation much simpler while in the Laplace transform and Taylor expansion method we first convert the integral equation to an algebraic equation using Laplace transform then we find its numerical inversion by power series. The numerical solution obtained by the proposed methods indicate that the approaches are easy computationally and its implementation very attractive. The methods are described and numerical examples are given to illustrate its accuracy and stability.

• Journal of Mechanical Science and Technology
• /
• 제16권4호
• /
• pp.522-531
• /
• 2002

In order to study the thermal-hydraulic behavior of the cable-in-conduit-conductor (CICC), a numerical model has been developed. In the model, the high heat transfer approximation between superconducting strands and supercritical helium is adopted. The strong coupling of heat transfer at the front of normal zone generates a contact discontinuity in temperature and density. In order to obtain the converged numerical solutions, a moving mesh method is used to capture the contact discontinuity in the short front region of the normal zone. The coupled equation is solved using the finite element method with the artificial viscosity term. Details of the numerical implementation are discussed and the validation of the code is performed for comparison of the results with thse of GANDALF and QSAIT.

• 대한조선학회논문집
• /
• 제41권1호
• /
• pp.23-30
• /
• 2004

In the Lagrangian vortex particle method based on the vorticity-velocity formulation for solving the incompressible Navier-Stokes equations, a numerical scheme for calculating pressure fields is presented. Implementation of the numerical method is directly connected with the well-established surface panel methods, just by dealing with the dynamic coupling among vorticity field. Assuming the vorticity and the velocity fields are to be calculated in time domain analysis, the pressure calculation for a complete set of solution at present time step is performed in a similar way to the one used in the Eulerian description. For a validation of the present method, we illustrate the early development of the viscous flow about an impulsive started circular cylinder for Reynolds number 550. The comparative study with the Eulerian finite Volume method provides an extensive understanding and application of the mesh-free Lagrangian vortex methods for numerical simulation of viscous flows around arbitrary bodies of general shape.

• Computers and Concrete
• /
• 제4권1호
• /
• pp.67-82
• /
• 2007

The discrete crack-concept is applied to study the 3D propagation of tensile-dominated failure in plain concrete. To this end the Partition of Unity Finite Element Method (PUFEM) is utilized and the strong discontinuity approach is followed. A consistent linearized implementation of the PUFEM is combined with a predictor-corrector algorithm to track the crack path, which leads to a robust numerical description of concrete cracking. The proposed concept is applied to study concrete failure during the PCT3D test and the predicted numerical results are compared to experimental data. The proposed numerical concept provides a clear interface for constitutive models and allows an investigation of their impact on concrete cracking under 3D conditions, which is of significant scientific interests to interpret results from 3D experiments.

• Coupled systems mechanics
• /
• 제7권6호
• /
• pp.669-690
• /
• 2018

This paper deals with numerical modeling of dynamic failure phenomena in rate-sensitive brittle and/or ductile materials. To this end, a two-dimensional continuum viscodamage-embedded discontinuity model, which is based on our previous work (see Do et al. 2017), is developed. More specifically, the pre-peak nonlinear and rate-sensitive hardening response of the material behavior, representing the fracture-process zone creation, is described by a rate-dependent continuum damage model. Meanwhile, an embedded displacement discontinuity model is used to formulate the post-peak response, involving the macro-crack creation accompanied by exponential softening. The numerical implementation in the context of the finite element method exploiting the second-order mid-point scheme is discussed in detail. In order to show the performance of the model several numerical examples are included.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.72-78
• /
• 2014

A hybrid particle-mesh method as the combination between the Vortex-In-Cell (VIC) method and penalization method has been achieved in recent years. The VIC method, which is based on the vorticity-velocity formulation, offers particle-mesh algorithms to numerically simulate flows past a solid body. The penalization method is used to enforce boundary conditions at a body surface with a decoupling between body boundaries and computational grids. The main advantage of the hybrid particle-mesh method is an efficient implementation for solid boundaries of arbitrary complexity on Cartesian grids. However, a numerical simulation of flows in large domains is still not too easy. In this study, a multi-domain approach is thus proposed to further reduce computation cost and easily implement it. We validate the implementation by numerical simulations of an incompressible viscous flow around an impulsively started circular cylinder.

• 한국전산유체공학회지
• /
• 제19권2호
• /
• pp.58-65
• /
• 2014

The adaptive wavelet method is studied for the enhancement of computational efficiency of three-dimensional flows. For implementation of the method for three-dimensional Euler equation, wavelet decomposition process is introduced based on the previous two-dimensional adaptive wavelet method. The order of numerical accuracy of an original solver is preserved by applying modified thresholding value. In order to assess the efficiency of the proposed algorithm, the method is applied to the computation of flow field around ONERA-M6 wing in transonic regime with 4th and 6th order interpolating polynomial respectively. Through the application, it is confirmed that the three-dimensional adaptive wavelet method can reduce the computational time while conserving the numerical accuracy of an original solver.