Method of regularized sources for axisymmetric Stokes flow problems

2016 ◽  
Vol 26 (3/4) ◽  
pp. 1226-1239 ◽  
Author(s):  
Kai Wang ◽  
Shiting Wen ◽  
Rizwan Zahoor ◽  
Ming Li ◽  
Božidar Šarler

Purpose – The purpose of this paper is to find solution of Stokes flow problems with Dirichlet and Neumann boundary conditions in axisymmetry using an efficient non-singular method of fundamental solutions that does not require an artificial boundary, i.e. source points of the fundamental solution coincide with the collocation points on the boundary. The fundamental solution of the Stokes pressure and velocity represents analytical solution of the flow due to a singular Dirac delta source in infinite space. Design/methodology/approach – Instead of the singular source, a non-singular source with a regularization parameter is employed. Regularized axisymmetric sources were derived from the regularized three-dimensional sources by integrating over the symmetry coordinate. The analytical expressions for related Stokes flow pressure and velocity around such regularized axisymmetric sources have been derived. The solution to the problem is sought as a linear combination of the fields due to the regularized sources that coincide with the boundary. The intensities of the sources are chosen in such a way that the solution complies with the boundary conditions. Findings – An axisymmetric driven cavity numerical example and the flow in a hollow tube and flow between two concentric tubes are chosen to assess the performance of the method. The results of the newly developed method of regularized sources in axisymmetry are compared with the results obtained by the fine-grid second-order classical finite difference method and analytical solution. The results converge with a finer discretization, however, as expected, they depend on the value of the regularization parameter. The method gives accurate results if the value of this parameter scales with the typical nodal distance on the boundary. Originality/value – Analytical expressions for the axisymmetric blobs are derived. The method of regularized sources is for the first time applied to axisymmetric Stokes flow problems.

2019 ◽  
Vol 29 (5) ◽  
pp. 1786-1814 ◽  
Author(s):  
A.A. Avramenko ◽  
N.P. Dmitrenko ◽  
I.V. Shevchuk ◽  
A.I. Tyrinov ◽  
V.I. Shevchuk

Purpose The paper aims to consider heat transfer in incompressible flow in a rotating flat microchannel with allowance for boundary slip conditions of the first and second order. The novelty of the paper encompasses analytical and numerical solutions of the problem, with the latter based on the lattice Boltzmann method (LBM). The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). An increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes a similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which takes into account the second-order boundary conditions. Design/methodology/approach The paper is based on the comparisons of an analytical solution and a numerical solution, which employs the lattice Boltzmann method. Both mathematical approaches used the first-order and second-order slip boundary conditions. The results obtained using both methods agree well with each other. Findings The analytical solution of the problem includes relations for the velocity and temperature profiles and for the Nusselt number depending on the rotation rate of the microchannel and slip velocity. It was demonstrated that the velocity profiles at high rotation rates transform from parabolic to M-shaped with a minimum at the channel axis. The temperature profiles tend to become uniform (i.e. almost constant). The increase in the channel rotation rate contributes to the increase in the Nusselt number. An increase in the Prandtl number causes the similar effect. The trend caused by the effect of the second-order slip boundary conditions depends on the closure hypothesis. It is shown that heat transfer in a flat microchannel can be successfully modeled using the LBM methodology, which considers the second-order boundary conditions. Originality/value The novelty of the paper encompasses analytical and numerical solutions of the problem, whereas the latter are based on the LBM.


2015 ◽  
Vol 7 (1) ◽  
pp. 13-30 ◽  
Author(s):  
Wenzhen Qu ◽  
Wen Chen

AbstractIn this paper, an improved singular boundary method (SBM), viewed as one kind of modified method of fundamental solution (MFS), is firstly applied for the numerical analysis of two-dimensional (2D) Stokes flow problems. The key issue of the SBM is the determination of the origin intensity factor used to remove the singularity of the fundamental solution and its derivatives. The new contribution of this study is that the origin intensity factors for the velocity, traction and pressure are derived, and based on that, the SBM formulations for 2D Stokes flow problems are presented. Several examples are provided to verify the correctness and robustness of the presented method. The numerical results clearly demonstrate the potentials of the present SBM for solving 2D Stokes flow problems.


2017 ◽  
Vol 17 (8) ◽  
pp. 128-133
Author(s):  
A.A. Bolshakov

In the paper the approximate solution for a problem of a rectangular plate on a two-parameter elastic base is suggested. The double series of beam functions satisfying elastic support boundary conditions are constructed. The analytical expressions for series function coefficients are obtained.


Author(s):  
Anis Younes ◽  
Ahmed Makradi ◽  
Ali Zidane ◽  
Qian Shao ◽  
Lyazid Bouhala

Purpose – The purpose of this paper is to develop an efficient non-iterative model combining advanced numerical methods for solving buoyancy-driven flow problems. Design/methodology/approach – The solution strategy is based on two independent numerical procedures. The Navier-Stokes equation is solved using the non-conforming Crouzeix-Raviart (CR) finite element method with an upstream approach for the non-linear convective term. The advection-diffusion heat equation is solved using a combination of Discontinuous Galerkin (DG) and Multi-Point Flux Approximation (MPFA) methods. To reduce the computational time due to the coupling, the authors use a non-iterative time stepping scheme where the time step length is controlled by the temporal truncation error. Findings – Advanced numerical methods have been successfully combined to solve buoyancy-driven flow problems on unstructured triangular meshes. The accuracy of the results has been verified using three test problems: first, a synthetic problem for which the authors developed a semi-analytical solution; second, natural convection of air in a square cavity with different Rayleigh numbers (103-108); and third, a transient natural convection problem of low Prandtl fluid with horizontal temperature gradient in a rectangular cavity. Originality/value – The proposed model is the first to combine advanced numerical methods (CR, DG, MPFA) for buoyancy-driven flow problems. It is also the first to use a non-iterative time stepping scheme based on local truncation error control for such coupled problems. The developed semi analytical solution based on Fourier series is also novel.


2013 ◽  
Vol 2013 ◽  
pp. 1-11 ◽  
Author(s):  
S. R. Mahmoud

The problem of generalized magneto-thermoelastic diffusion in an infinite rotating nonhomogeneity medium subjected to certain boundary conditions is studied. The chemical potential is also assumed to be a known function of time at the boundary of the cavity. The analytical expressions for the displacements, stresses, temperature, concentration, and chemical potential are obtained. Comparison was made between the results obtained in the presence and absence of diffusion. The results indicate that the effect of nonhomogeneity, rotation, magnetic field, relaxation time, and diffusion is very pronounced.


2015 ◽  
Vol 56 (10) ◽  
pp. 103101 ◽  
Author(s):  
Quy-Dong To ◽  
Van-Huyen Vu ◽  
Guy Lauriat ◽  
Céline Léonard

1975 ◽  
Vol 67 (4) ◽  
pp. 787-815 ◽  
Author(s):  
Allen T. Chwang ◽  
T. Yao-Tsu Wu

The present study further explores the fundamental singular solutions for Stokes flow that can be useful for constructing solutions over a wide range of free-stream profiles and body shapes. The primary singularity is the Stokeslet, which is associated with a singular point force embedded in a Stokes flow. From its derivatives other fundamental singularities can be obtained, including rotlets, stresslets, potential doublets and higher-order poles derived from them. For treating interior Stokes-flow problems new fundamental solutions are introduced; they include the Stokeson and its derivatives, called the roton and stresson.These fundamental singularities are employed here to construct exact solutions to a number of exterior and interior Stokes-flow problems for several specific body shapes translating and rotating in a viscous fluid which may itself be providing a primary flow. The different primary flows considered here include the uniform stream, shear flows, parabolic profiles and extensional flows (hyper-bolic profiles), while the body shapes cover prolate spheroids, spheres and circular cylinders. The salient features of these exact solutions (all obtained in closed form) regarding the types of singularities required for the construction of a solution in each specific case, their distribution densities and the range of validity of the solution, which may depend on the characteristic Reynolds numbers and governing geometrical parameters, are discussed.


Sign in / Sign up

Export Citation Format

Share Document