• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.5 no.1
• /
• pp.55-78
• /
• 2001

Given the anisotropic Poisson equation $-{\nabla}{\cdot}{\mathcal{K}}{\nabla}p=f$, one can convert it into a system of two first order PDEs: the Darcy law for the flux $u=-{\mathcal{K}{\nabla}p$ and conservation of mass ${\nabla}{\cdot}u=f$. A very natural mixed finite volume method for this system is to seek the pressure in the nonconforming P1 space and the Darcy velocity in the lowest order Raviart-Thomas space. The equations for these variables are obtained by integrating the two first order systems over the triangular volumes. In this paper we show that such a method is really a standard finite element method with local recovery of the flux in disguise. As a consequence, we compare two approaches in analyzing finite volume methods (FVM) and shed light on the proper way of analyzing non co-volume type of FVM. Numerical results for Dirichlet and Neumann problems are included.

• Structural Engineering and Mechanics
• /
• v.25 no.2
• /
• pp.181-200
• /
• 2007

A decoupling technique for simulating near-field wave motions in two-phase media is introduced in this paper. First, an equivalent but direct weighted residual method is presented in this paper to solve boundary value problems more explicitly. We applied the Green's theorem for integration by parts on the equivalent integral statement of the field governing equations and then introduced the Neumann conditions directly. Using this method and considering the precision requirement in wave motion simulation, a lumped-mass FEM for two-phase media with clear physical concepts and convenient implementation is derived. Then, considering the innate attenuation character of the wave in two-phase media, an attenuation parameter is introduced into Liao's Multi-Transmitting Formula (MTF) to simulate the attenuating outgoing wave in two-phase media. At last, two numerical experiments are presented and the numerical results are compared with the analytical ones demonstrating that the lumped-mass FEM and the generalized MTF introduced in this paper have good precision.

• Communications of the Korean Mathematical Society
• /
• v.35 no.4
• /
• pp.1143-1158
• /
• 2020

In this paper, we are concerned with the blow-up phenomena for a class of pseudo-parabolic equations with weighted source u_t - △u - △u_t = a(x)f(u) subject to Dirichlet (or Neumann) boundary conditions in any smooth bounded domain Ω ⊂ ℝⁿ (n ≥ 1). Firstly, we obtain the upper and lower bounds for blow-up time of solutions to these problems. Moreover, we also give the estimates of blow-up rate of solutions under some suitable conditions. Finally, three models are presented to illustrate our main results. In some special cases, we can even get some exact values of blow-up time and blow-up rate.

• Journal of applied mathematics & informatics
• /
• v.28 no.3_4
• /
• pp.611-624
• /
• 2010

A Numerical scheme based on collocation method using quartic B-spline functions is designed for the numerical solution of one-dimensional modified equal width wave (MEW) wave equation. Using Von-Neumann approach the scheme is shown to be unconditionally stable. Performance of the method is validated through test problems including single wave, interaction of two waves and use of Maxwellian initial condition. Using error norms $L_2$ and $L_{\infty}$ and conservative properties of mass, momentum and energy, accuracy and efficiency of the suggested method is established through comparison with the existing numerical techniques.

• Nuclear Engineering and Technology
• /
• v.35 no.1
• /
• pp.45-54
• /
• 2003

A hyperbolic two-phase, two-fluid equation system developed in the previous work has been implemented in an existing nuclear safety analysis code, MARS. Although the implicit treatment of interfacial pressure force term introduced in momentum equation of the hyperbolic equation system is required to enhance the numerical stability, it is very difficult to implement in the code because it is not possible to maintain the existing numerical solution structure. As an alternative, two-step approach with stabilizer momentum equations has been selected. The results of a linear stability analysis by Von-Neumann method show the equivalent stability improvement with fully-implicit solution method. To illustrate the applicability, the new solution scheme has been implemented into the best-estimate thermal-hydraulic analysis code, MARS. This paper also includes the comparisons of the simulation results for the perturbation propagation and water faucet problems using both two-step method and the original solution scheme.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.7 no.1
• /
• pp.27-36
• /
• 1970

선박유체역학(船舶流體力學)의 경계조건(境界條件)의 선형화(線型化)를 위한 퍼터베이션(Perturbation)방법(方法)을 따라 배를 가늘고 길게 보느냐, 엷은 수평판(水平板)으로 보느냐 등등에 따라 여러가지 다른 종류(種類)의 수학적(數學的) 정의(定義)에 따르는 이론(理論)을 추궁함으로서, 1) Neumann Ship, 2) Thin Ship, 3) Michell Ship, 4) Slender Ship 등을 이론적(理論的)인 면(面)에서와 수치계산(數値計算) 예(例)로서 비교검토(比較檢討) 하였다. 동시(同時)에 이에 따르는 "직접문제(直接問題)" "간접문제(間接問題)" 등(等)에 관(關)하여도 언급(言及)하고 그 차이점(差異點)을 지적하였다. 이에 의(依)하면 정속전진시(定速前進時)에 L/B가 10인 경우에는 최대광폭의 차(差)가 $5{\sim}8%$정도 되고 또 선수수선반각(船首水線半角)에 상당한 차(差)가 있음으로 같은 이론적(理論的)인 계산치(計算値)에 대(對)하여 저항(抵抗)은 적어도 $10{\sim}15%$정도 차(差)가 생길 것이 아닌가 보게끔 되었다.

• Honam Mathematical Journal
• /
• v.31 no.4
• /
• pp.579-590
• /
• 2009

The precise form of singular functions, singular function representation and the extraction form for the stress intensity factor play an important role in the singular function methods to deal with the domain singularities for the Poisson problems with most common boundary conditions, e.q. Dirichlet or Mixed boundary condition [2, 4]. In this paper we give an elementary step to get the singular functions of the solution for Poisson problem with Neumann boundary condition or Robin boundary condition. We also give singular function representation and the extraction form for the stress intensity with a result showing the number of singular functions depending on the boundary conditions.

• Structural Engineering and Mechanics
• /
• v.35 no.4
• /
• pp.511-533
• /
• 2010

As a truly meshfree method, meshless equilibrium on line method (MELM), for 2D elasticity problems is presented. In MELM, the problem domain is represented by a set of distributed nodes, and equilibrium is satisfied on lines for any node within this domain. In contrary to conventional meshfree methods, test domains are lines in this method, and all integrals can be easily evaluated over straight lines along x and y directions. Proposed weak formulation has the same concept as the equilibrium on line method which was previously used by the authors for enforcement of the Neumann boundary conditions in the strong-form meshless methods. In this paper, the idea of the equilibrium on line method is developed to use as the weak forms of the governing equations at inner nodes of the problem domain. The moving least squares (MLS) approximation is used to interpolate solution variables in this paper. Numerical studies have shown that this method is simple to implement, while leading to accurate results.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.73-80
• /
• 2024

The aim of this paper is to remove the rigid body motion in the interior boundary value problem (BVP) of plane elasticity by solving the interior and exterior BVPs simultaneously. First, we formulate the interior and exterior BVPs simultaneously. The tractions applied on the contour in two problems are the same. After adding and subtracting the two boundary integral equations (BIEs), we will obtain a couple of BIEs. In the coupled BIEs, the properties of relevant integral operators are modified, and those integral operators are generally invertible. Finally, a unique solution for boundary displacement of interior region can be obtained.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.42-57
• /
• 1995

An efficient and accurate scheme has been constructed by taking advantages of the hi-quadratic spline scheme and the higher-order boundary element method selectively depending on computation domains. Boundary surfaces are represented by 8-node boundary elements to describe curved surfaces of a ship and its neighboring free surface more accurately. The variation of the velocity potential complies with the characteristics of the 8-node element on the body surface. But on the free surface, it is assumed to follow that of the hi-quadratic spline scheme. By which, the free surface solution is free from numerical damping and has better numerical dispersion property. As numerical examples, steady and unsteady Neumann-Kelvin problems are considered. Numerical results for a submerged spheroid, Series 60($C_B=0.6$) and a modified support the proposed method. Finally, a new upstream radiation condition is derived using a wave equation operator in order to deal with problems for subcritical reduced frequency. The relevance of this operator has been confirmed in the case of unsteady Kelvin source potential.