Meshless Local Petrov-Galerkin (MLPG) Method for Incompressible Viscous Fluid Flows

2006 ◽  
Author(s):  
Mohammad Haji Mohammadi

In this paper, the truly Meshless Local Petrov-Galerkin (MLPG) method is extended for computation of unsteady incompressible flows, governed by the Navier–Stokes equations (NSE), in vorticity-stream function formulation. The present method is a truly meshless method based only on a number of randomly located nodes. The formulation is based on two equations including stream function Poisson equation and vorticity advection-dispersion-reaction eq. (ADRE). The meshless method is based on a local weighted residual method with the Heaviside step function and quartic spline as the test functions respectively over a local subdomain. Moving Least Square approximation (MLS) is employed in shape function construction for approximation of a gauss point. Due to dissatisfaction of kronecker delta property in MLS approximation, the penalty method is employed to enforce the essential boundary conditions. In order to overcome instability and numerical errors encountering in convection dominant flows, a new upwinding scheme is used to stabilize the convection operator in the streamline direction (as is done in SUPG). In this upwinding technic, instead of moving subdomains the weight function is shifted in the direction of flow. The efficiency, accuracy and robustness are demonstrated by some test problems, including the standard driven cavity together with the driven cavity flow in an L shaped cavity and flow in a Z shaped channel. The comparison of computational results shows that the developed method is capable of accurate resolution of flow fields in complex geometries.

2019 ◽  
Vol 286 ◽  
pp. 07003
Author(s):  
M. Rammane ◽  
S. Mesmoudi ◽  
A. Tri ◽  
B. Braikat ◽  
N. Damil

In this work, we propose to investigate numerically the incompressible flows by the Asymptotic Numerical Method (ANM) with the Moving Least Square (MLS). The mathematical formulation is based on theNavier-Stokes equations written in a strongly formulation to avoid all difficulties of the numerical integration. The used algorithm is developed to investigate the effective of the ANM with the MLS in the strongly formulation.


2016 ◽  
Vol 2016 ◽  
pp. 1-8 ◽  
Author(s):  
Hongfen Gao ◽  
Gaofeng Wei

Combining the finite covering technical and complex variable moving least square, the complex variable meshless manifold method can handle the discontinuous problem effectively. In this paper, the complex variable meshless method is applied to solve the problem of elastic dynamics, the complex variable meshless manifold method for dynamics is established, and the corresponding formula is derived. The numerical example shows that the numerical solutions are in good agreement with the analytical solution. The CVMMM for elastic dynamics and the discrete forms are correct and feasible. Compared with the traditional meshless manifold method, the CVMMM has higher accuracy in the same distribution of nodes.


2021 ◽  
Vol 8 (3) ◽  
pp. 418-424
Author(s):  
Syed Fazuruddin ◽  
Seelam Sreekanth ◽  
G. Sankara Sekhar Raju

Incompressible 2-D Navier-stokes equations for various values of Reynolds number with and without partial slip conditions are studied numerically. The Lid-Driven cavity (LDC) with uniform driven lid problem is employed with vorticity - Stream function (VSF) approach. The uniform mesh grid is used in finite difference approximation for solving the governing Navier-stokes equations and developed MATLAB code. The numerical method is validated with benchmark results. The present work is focused on the analysis of lid driven cavity flow of incompressible fluid with partial slip conditions (imposed on side walls of the cavity). The fluid flow patterns are studied with wide range of Reynolds number and slip parameters.


2012 ◽  
Vol 446-449 ◽  
pp. 3633-3638
Author(s):  
Yu Ling Jiao ◽  
Guang Wei Meng ◽  
Xu Xi Qin

moving least square meshless method is a numerical approximation based on points that do not generate the grid of cells, as long as the node information. Basis function and weight function meshless method for the calculation of accuracy have a significant impact. In order to compare the order of the base functions and powers of the radius of influence domain function meshless method for computational accuracy and efficiency , this paper selected first, second and third basis function and spline-type weight function in a different influence domain radius, respectively construct the field function. Mindlin plate element is derived based on the format of the plate bending problem meshless discrete equations. Programming examples are calculated with elastic plate bending problems non-grid solutions, and analysis and comparison of their accuracy and efficiency, results show that the meshless method using elastic plate bending problem is feasible and effective.


2012 ◽  
Vol 2012 ◽  
pp. 1-17 ◽  
Author(s):  
K. Poochinapan

The aim of this paper is to study the properties of approximations to nonlinear terms of the 2D incompressible Navier-Stokes equations in the stream function formulation (time-dependent biharmonic equation). The nonlinear convective terms are numerically solved by using the method with internal iterations, compared to the ones which are solved by using explicit and implicit schemes (operator splitting scheme Christov and Marinova; (2001)). Using schemes and algorithms, the steady 2D incompressible flow in a lid-driven cavity is solved up to Reynolds number Re =5000 with second-order spatial accuracy. The schemes are thoroughly validated on grids with different resolutions. The result of numerical experiments shows that the finite difference scheme with internal iterations on nonlinearity is more efficient for the high Reynolds number.


2013 ◽  
Vol 394 ◽  
pp. 128-133
Author(s):  
Yuan Ding Wang ◽  
Jun Jie Tan ◽  
Xiao Wei Cai ◽  
Deng Feng Ren

Large Eddy Simulation (LES) based on the least square meshless method was proposed in the present paper to simulate the classical turbulent flow around a stationary 2D circular cylinder. The subgrid scale model of Smagorinsky-Lily was employed to close the Navier-Stokes equations filtered by Favre filter. The Reynolds number is 3900 which means that the flow is subcritical and the wake is fully turbulent but the cylinder boundary is still laminar. Results obtained in this paper were evaluated by comparison with published experimental results and other numerical results. The results obtained in the present work show better agreement with the experimental values than other two-dimensional LES results .


2013 ◽  
Vol 2013 ◽  
pp. 1-14 ◽  
Author(s):  
Mohamed A. El-Beltagy ◽  
Mohamed I. Wafa

A two-dimensional stochastic solver for the incompressible Navier-Stokes equations is developed. The vorticity-stream function formulation is considered. The polynomial chaos expansion was integrated with an unstructured node-centered finite-volume solver. A second-order upwind scheme is used in the convection term for numerical stability and higher-order discretization. The resulting sparse linear system is solved efficiently by a direct parallel solver. The mean and variance simulations of the cavity flow are done for random variation of the viscosity and the lid velocity. The solver was tested and compared with the Monte-Carlo simulations and with previous research works. The developed solver is proved to be efficient in simulating the stochastic two-dimensional incompressible flows.


2017 ◽  
Vol 7 (4) ◽  
pp. 714-727 ◽  
Author(s):  
Dongdong He ◽  
Kejia Pan

AbstractIncompressible flows with zero Reynolds number can be modeled by the Stokes equations. When numerically solving the Stokes flow in stream-vorticity formulation with high-order accuracy, it will be important to solve both the stream function and velocity components with the high-order accuracy simultaneously. In this work, we will develop a fifth-order spectral/combined compact difference (CCD) method for the Stokes equation in stream-vorticity formulation on the polar geometries, including a unit disk and an annular domain. We first use the truncated Fourier series to derive a coupled system of singular ordinary differential equations for the Fourier coefficients, then use a shifted grid to handle the coordinate singularity without pole condition. More importantly, a three-point CCD scheme is developed to solve the obtained system of differential equations. Numerical results are presented to show that the proposed spectral/CCD method can obtain all physical quantities in the Stokes flow, including the stream function and vorticity function as well as all velocity components, with fifth-order accuracy, which is much more accurate and efficient than low-order methods in the literature.


2014 ◽  
Vol 2014 ◽  
pp. 1-11 ◽  
Author(s):  
Fen Li ◽  
Xiaolin Li

The interpolating boundary element-free method (IBEFM) is developed in this paper for boundary-only analysis of unilateral problems which appear in variational inequalities. The IBEFM is a direct boundary only meshless method that combines an improved interpolating moving least-square scheme for constructing interpolation functions with boundary integral equations (BIEs) for representing governing equations. A projection operator is used to formulate the BIEs and then the formulae of the IBEFM are obtained for unilateral problems. The convergence of the developed meshless method is derived mathematically. The capability of the method is also illustrated and assessed through some numerical experiments.


Sign in / Sign up

Export Citation Format

Share Document