• Journal of computational fluids engineering
• /
• v.21 no.2
• /
• pp.1-11
• /
• 2016

Considerable efforts are required to develop a monotone, robust and stable high-order numerical scheme for solving the hyperbolic system. The discontinuous Galerkin(DG) method is a natural choice, but elimination of the spurious oscillations from the high-order solutions demands a new development of proper limiters for the DG method. There are several available limiters for controlling or removing unphysical oscillations from the high-order approximate solution; however, very few studies were directed to analyze the exact role of the limiters in the hyperbolic systems. In this study, the performance of the several well-known limiters is examined by comparing the high-order($p^1$, $p^2$, and $p^3$) approximate solutions with the exact solutions. It is shown that the accuracy of the limiter is in general problem-dependent, although the Hermite WENO limiter and maximum principle limiter perform better than the TVD and generalized moment limiters for most of the test cases. It is also shown that application of the troubled cell indicators may improve the accuracy of the limiters under some specific conditions.

• Journal of the Society of Naval Architects of Korea
• /
• v.31 no.1
• /
• pp.63-70
• /
• 1994

The three-dimensional incompressible clavier-Stokes equations are solved to simulate the flow field around a Wigley model with free-surface. The IAF(Implicit Approximate Factorization) method is used to show a good success in reducing the computing time. The CPU time is almost an half of that if the IAF method were used. The present method adopts the local linearization and Euler implicit scheme without the pressure-gradient terms for the artificial viscosity. Calculations are carried out at the Reynolds number of $10^6$ and the Froude numbers are 0.25, 0.289 and 0.316. For the approximations of turbulence, the Baldwin-Lomax model is used. The resulting free-surface wave configurations and the velocity vectors are compared with those by the explicit method and experiments.

• Steel and Composite Structures
• /
• v.22 no.1
• /
• pp.161-182
• /
• 2016

A state space differential reproducing kernel (DRK) method is developed for the three-dimensional (3D) analysis of functionally graded material (FGM) axisymmetric circular plates with simply-supported and clamped edges. The strong formulation of this 3D elasticity axisymmetric problem is derived on the basis of the Reissner mixed variational theorem (RMVT), which consists of the Euler-Lagrange equations of this problem and its associated boundary conditions. The primary field variables are naturally independent of the circumferential coordinate, then interpolated in the radial coordinate using the early proposed DRK interpolation functions, and finally the state space equations of this problem are obtained, which represent a system of ordinary differential equations in the thickness coordinate. The state space DRK solutions can then be obtained by means of the transfer matrix method. The accuracy and convergence of this method are examined by comparing their solutions with the accurate ones available in the literature.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.9 no.1
• /
• pp.79-90
• /
• 2005

In this paper we study the numerical solutions of the stochastic differential equations of the form $$du(x,\;t)=f(x,\;t,\;u)dt\;+\;g(x,\;t,\;u)dW(t)\;+\;\sum\limits_{|q|\leq2m}\;A_q(x,\;t)D^qu(x,\;t)dt$$ where $0\;{\leq}\;t\;{\leq}\;T,\;x\;{\in}\;R^{\nu}$, ($R^{nu}$ is the $\nu$-dimensional Euclidean space). Here $u\;{\in}\;R^n$, W(t) is an n-dimensional Brownian motion, $$f\;:\;R^{n+\nu+1}\;{\rightarrow}\;R^n,\;g\;:\;R^{n+\nu+1}\;{\rightarrow}\;R^{n{\times}n},$$, and $$A_q\;:\;R^{\nu}\;{\times}\;[0,\;T]\;{\rightarrow}\;R^{n{\times}n}$$ where ($A_q,\;|\;q\;|{\leq}\;2m$) is a family of square matrices whose elements are sufficiently smooth functions on $R^{\nu}\;{\times}\;[0,\;T]\;and\;D^q\;=\;D^{q_1}_1_{\ldots}_{\ldots}D^{q_{\nu}}_{\nu},\;D_i\;=\;{\frac{\partial}{\partial_{x_i}}}$.

• Journal of Integrative Natural Science
• /
• v.6 no.1
• /
• pp.53-56
• /
• 2013

We are interested in the rate of convergence of solutions of 2D Navier-Stokes equations in a smooth bounded domain as the viscosity tends to zero under Navier friction condition. If the initial velocity is smooth enough($u{\in}W^{2,p}$, p>2), it is known that the rate of convergence is linearly propotional to the viscosity. Here, we consider the rate of convergence for nonsmooth velocity fields when the gradient of the corresponding solution of the Euler equations belongs to certain Orlicz spaces. As a corollary, if the initial vorticity is bounded and small enough, we obtain a sublinear rate of convergence.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.10 no.1
• /
• pp.102-109
• /
• 1986

Based on the Hartenberg-Denavit symbolic equation, which is one of equations for the kinematic analysis of three dimensional (3-D) linkage, a generalized kinematic motion equation is derived utilizing Euler angles and employing the coordinates transformation. The derived equation can feasibly be used for the motion analysis of any type of 3-D linkages as well as 2-D ones. In order to simulate the general motion of 3-D linkgages on digital computer, the generalized equation is programmed through the process of numerical analysis after converting the equation to the type of Newton-Raphson formula and denoting it in matrix form. The feasibility of theoretically derived equation is experimentally proved by comparing the results from the computer with those from experimental setup of three differrent but generally empolyed 3-D linkages.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.165-168
• /
• 2002

This paper describes the dynamics of the weak shock wave propagating inside some kinds of branched pipe bends. Computations are carried out by solving the two-dimensional, compressible, unsteady Euler Equations. The second-order TVD(Total Variation Diminishing) scheme is employed to discretize the governing equations. For computations, two types of branched pipe($90^{\circ}$ branch,$45^{\circ}$ branch) with a diameter of D are used. The incident normal shock wave is assumed at D upstream of the pipe bend entrance, and its Mach number is changed between 1.1 and 2.4. The flow fields are numerically visualized by using the pressure contours and computed schlieren images. The comparison with the experimental data performed for the purpose of validation of computational work. Reflection and diffraction of the propagating shock wave are clarified. The present computations predicted the experimented flow field with a good accuracy.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.456-461
• /
• 2001

This paper depicts the weak shock wave propagating inside some kinds of pipe bends. Computational work is to solve the two-dimensional, compressible, unsteady Euler Equations. The second-order TVD scheme is employed to discretize the governing equations. For the computations, the incident normal shock wave is assumed at the entrance of the pipe bend, and its Mach number is changed between 1.1 and 1.7. The turning angle and radius of the curvature of the pipe bend are changed to investigate the effects on the shock wave structure. The present computational results clearly show the shock wave reflection and diffraction occurring in the pipe bend. In particular, the vortex generation, which occurs at the edge of the bend, and its shedding mechanism are discussed in details.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.18 no.11
• /
• pp.1157-1169
• /
• 2008

This paper presents spectral element formulation which approximates Lamb wave propagation by PZT transducers bonded on a thin plate. A two layer beam model under 2-D plane strain condition is introduced to simulate high-frequency dynamic responses induced by a piezoelectric (PZT) layer rigidly bonded on a base plate. Mindlin-Herrmann and Timoshenko beam theories are employed to represent the first symmetric and anti-symmetric Lamb wave modes on a base plate, respectively. The Euler-Bernoulli beam theory and 1-D linear piezoelectricity are used to model the electro-mechanical behavior of a PZT layer. The equations of motions of a two layer beam model are derived through Hamilton's principle. The necessary boundary conditions associated with the electro-mechanical properties of a PZT layer are formulated in the context of dual functions of a PZT layer as an actuator and a sensor. General spectral shape functions of response field and the associated boundary conditions are obtained through equations of motions converted into frequency domain. Detailed spectrum element formulation for composing the dynamic stiffness matrix of a two layer beam model is presented as well. The validity of the proposed spectral element is demonstrated through numerical examples.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.7 no.3
• /
• pp.110-117
• /
• 1998

A continuum-based configuration design sensitivity analysis method is developed for kinematically driven mechanical systems. The configuration design variable for mechanical systems is defined. The 3-1-3 Euler angle is employed as the orientation design variable. Kinematic admissibility conditions of configuration design change. Direct differentiation method is used to derive the governing equations of the design sensitivity. Numerical examples are presented to demonstrate the validity and effectiveness of the proposed method.