• Journal of the Korean Mathematical Society
• /
• v.57 no.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.

• Bulletin of the Korean Mathematical Society
• /
• v.53 no.3
• /
• pp.853-874
• /
• 2016

Based on a particular overlapping domain decomposition technique, a parallel finite element discretization algorithm for the generalized Stokes equations is proposed and investigated. In this algorithm, each processor computes a local approximate solution in its own subdomain by solving a global problem on a mesh that is fine around its own subdomain and coarse elsewhere, and hence avoids communication with other processors in the process of computations. This algorithm has low communication complexity. It only requires the application of an existing sequential solver on the global meshes associated with each subdomain, and hence can reuse existing sequential software. Numerical results are given to demonstrate the effectiveness of the parallel algorithm.

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.21 no.3
• /
• pp.109-124
• /
• 2017

The numerical solution of the Stokes equation with discontinuous viscosity and singular force term is challenging, due to the discontinuity of pressure, non-smoothness of velocity, and coupled discontinuities along interface.In this paper, we give an efficient algorithm to solve this problem by employing spectral Legendre and Chebyshev approximations.First, we present the algorithm for a problem defined in rectangular domain with straight line interface. Then it is generalized to a domain with smooth curve boundary and interface by employing spectral element method. Numerical experiments demonstrate the accuracy and efficiency of our algorithm and its spectral convergence.

• Journal of applied mathematics & informatics
• /
• v.37 no.5_6
• /
• pp.507-519
• /
• 2019

The ground for this paper is to examine the generalized Stokes' first and second issues for an incompressible couple pressure liquid under isothermal conditions. Exact solutions for each problem are acquired by using the Laplace transform (LT) with respect to the time variable t and the sine Fourier transform (FT) with respect to the y-variable. Further, a comparison is given of the obtained results and the results of Devakar and Lyengar [1] and by using the four inverse Laplace transform algorithms (Stehfest's, Tzou's, Talbot, Fourier series) in the space time domain utilizing a numerical methodology. Moreover, velocity profiles are plotted and considered for various occasions and distinctive estimations of couple stress parameters. At the end, the outcomes are exhibited by graphs and in tabular forms.

• Honam Mathematical Journal
• /
• v.43 no.2
• /
• pp.183-197
• /
• 2021

Recently several authors have extended the Beta function by using its integral representation. However, in many cases no expression of these extended functions in terms of classic special functions is known. In the present paper, we introduce a further extension by defining a family of functions G_r,s : ℝ^*₊ → ℂ, with r, s ∈ ℂ and ℜ(r) > 0. For given r, s, we prove that this function satisfies a second-order linear differential equation with rational coefficients. Solving this ODE, we express G_r,s as a combination of confluent hypergeometric functions. From this we deduce a new integral relation satisfied by Tricomi's function. We then investigate additional specific properties of G_r,1 which take the form of new non trivial integral relations involving exponential and error functions. We discuss the connection between G_r,1 and Stokes' first problem (or Rayleigh problem) in fluid mechanics which consists in determining the flow created by the movement of an infinitely long plate. For $r{\in}{\frac{1}{2}}{\mathbb{N}}^*$, we find additional relations between G_r,1 and Hermite polynomials. In view of these results, we believe the family of extended beta functions G_r,s will find further applications in two directions: (i) for improving our knowledge of confluent hypergeometric functions and Tricomi's function, (ii) and for engineering and physics problems.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.1 s.232
• /
• pp.80-86
• /
• 2005

This study presents numerical solutions of the two-dimensional Navier-Stokes equations for supersonic unsteady flow in a convergent-divergent nozzle with heat addition. The TVD scheme in generalized coordinates is employed in order to calculate the moving shock waves caused by thermal choking. We discuss on transient characteristics, start and unstart phenomena, fluctuations of specific thrust caused by thermal choking and viscous effects. We prove that the control of separation of boundary layer is the most important key problem to prevent the thermal choking.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.