Neural networks for BEM analysis of steady viscous flows

This paper presents a new neural network‐boundary integral approach for analysis of steady viscous fluid flows. Indirect radial basis function networks (IRBFNs) which perform better than element‐based methods for function interpolation, are introduced into the BEM scheme to represent the variations of velocity and traction along the boundary from the nodal values. In order to assess the effect of IRBFNs, the other features used in the present work remain the same as those used in the standard BEM. For example, Picard‐type scheme is utilized in the iterative procedure to deal with the non‐linear convective terms while the calculation of volume integrals and velocity gradients are based on the linear finite element‐based method. The proposed IRBFN‐BEM is verified on the driven cavity viscous flow problem and can achieve a moderate Reynolds number of 1400 using a relatively coarse uniform mesh. The results obtained such as the velocity profiles along the horizontal and vertical centrelines as well as the properties of the primary vortex are in very good agreement with the benchmark solution. Furthermore, the secondary vortices are also captured by the present method. Thus, it appears that an ability to represent the boundary solution accurately can significantly improve the overall solution accuracy of the BEM. Copyright © 2003 John Wiley & Sons, Ltd.     

A multigrid approach for steady state laminar viscous flows

In this paper, we use the laminar viscous flow in a lid‐driven cavity as an example to describe and verify a numerical scheme for non‐linear partial differential equations. The proposed scheme combines a new analytical method for strongly non‐linear problems, namely the homotopy analysis method, with the multigrid techniques. A family of formulas at different orders is given. At the lowest order, the current approach is the same as the traditional multigrid methods. However, our high‐order scheme needs a fewer number of iterations and less CPU time than the classical ones. Copyright © 2001 John Wiley & Sons, Ltd.     

Meshless Galerkin analysis of Stokes slip flow with boundary integral equations

This paper presents a novel meshless Galerkin scheme for modeling incompressible slip Stokes flows in 2D. The boundary value problem is reformulated as boundary integral equations of the first kind which is then converted into an equivalent variational problem with constraint. We introduce a Lagrangian multiplier to incorporate the constraint and apply the moving least‐squares approximations to generate trial and test functions. In this boundary‐type meshless method, boundary conditions can be implemented exactly and system matrices are symmetric. Unlike the domain‐type method, this Galerkin scheme requires only a nodal structure on the bounding surface of a body for approximation of boundary unknowns. The convergence and abstract error estimates of this new approach are given. Numerical examples are also presented to show the efficiency of the method. Copyright © 2009 John Wiley & Sons, Ltd.     

Some methods of training radial basis neural networks in solving the Navier‐Stokes equations

The purpose of this research is to analyze the application of neural networks and specific features of training radial basis functions for solving 2‐dimensional Navier‐Stokes equations. The authors developed an algorithm for solving hydrodynamic equations with representation of their solution by the method of weighted residuals upon the general neural network approximation throughout the entire computational domain. The article deals with testing of the developed algorithm through solving the 2‐dimensional Navier‐Stokes equations. Artificial neural networks are widely used for solving problems of mathematical physics; however, their use for modeling of hydrodynamic problems is very limited. At the same time, the problem of hydrodynamic modeling can be solved through neural network modeling, and our study demonstrates an example of its solution. The choice of neural networks based on radial basis functions is due to the ease of implementation and organization of the training process, the accuracy of the approximations, and smoothness of solutions. Radial basis neural networks in the solution of differential equations in partial derivatives allow obtaining a sufficiently accurate solution with a relatively small size of the neural network model. The authors propose to consider the neural network as an approximation of the unknown solution of the equation. The Gaussian distribution is used as the activation function.     

Fitted finite element discretization of two‐phase Stokes flow

We propose a novel fitted finite element method for two‐phase Stokes flow problems that uses piecewise linear finite elements to approximate the moving interface. The method can be shown to be unconditionally stable. Moreover, spherical stationary solutions are captured exactly by the numerical approximation. In addition, the meshes describing the discrete interface in general do not deteriorate in time, which means that in numerical simulations, a smoothing or a remeshing of the interface mesh is not necessary. We present several numerical experiments for our numerical method, which demonstrate the accuracy and robustness of the proposed algorithm. Copyright © 2016 John Wiley & Sons, Ltd.     

A new immersed boundary method for compressible Navier–Stokes equations

This paper presents an immersed boundary method for compressible Navier–Stokes equations in irregular domains, based on a local radial basis function approximation. This approach allows one to define a reconstruction of the radial basis functions on each irregular interface cell to treat both the Dirichlet and Neumann boundary conditions accurately on the immersed interfaces. Several numerical examples, including problems with available analytical solutions and the well-documented flow past an airfoil, are presented to test the proposed method. The numerical results demonstrate that the proposed method provides accurate solutions for viscous compressible flows.     

A simple finite element method for Stokes flows with surface tension using unfitted meshes

This paper presents a simple finite element method for Stokes flows with surface tension. The method uses an unfitted mesh that is independent of the interface. Due to the surface force, the pressure has a jump across the interface. Based on the properties of the level set function that implicitly represents the interface, the jump of the pressure is removed, and a new problem without discontinuities is formulated. Then, classical stable finite element methods are applied to solve the new problem. Some techniques are used to show that the method is equivalent to an easy‐to‐implement method that can be regarded as a traditional method with a modified pressure space. However, the matrix of the resulting linear system of equations is the same as that of the traditional method. Optimal error estimates are derived for the proposed method. Finally, some numerical tests are presented to confirm the theoretical results. Copyright © 2015 John Wiley & Sons, Ltd.     

Boundary element analysis of driven cavity flow for low and moderate Reynolds numbers

A boundary element method for steady two‐dimensional low‐to‐moderate‐Reynolds number flows of incompressible fluids, using primitive variables, is presented. The velocity gradients in the Navier–Stokes equations are evaluated using the alternatives of upwind and central finite difference approximations, and derivatives of finite element shape functions. A direct iterative scheme is used to cope with the non‐linear character of the integral equations. In order to achieve convergence, an underrelaxation technique is employed at relatively high Reynolds numbers. Driven cavity flow in a square domain is considered to validate the proposed method by comparison with other published data. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical solution of Navier–Stokes equations using multiquadric radial basis function networks

A numerical method based on radial basis function networks (RBFNs) for solving steady incompressible viscous flow problems (including Boussinesq materials) is presented in this paper. The method uses a ‘universal approximator’ based on neural network methodology to represent the solutions. The method is easy to implement and does not require any kind of ‘finite element‐type’ discretization of the domain and its boundary. Instead, two sets of random points distributed throughout the domain and on the boundary are required. The first set defines the centres of the RBFNs and the second defines the collocation points. The two sets of points can be different; however, experience shows that if the two sets are the same better results are obtained. In this work the two sets are identical and hence commonly referred to as the set of centres. Planar Poiseuille, driven cavity and natural convection flows are simulated to verify the method. The numerical solutions obtained using only relatively low densities of centres are in good agreement with analytical and benchmark solutions available in the literature. With uniformly distributed centres, the method achieves Reynolds number Re = 100 000 for the Poiseuille flow (assuming that laminar flow can be maintained) using the density of , Re = 400 for the driven cavity flow with a density of and Rayleigh number Ra = 1 000 000 for the natural convection flow with a density of . Copyright © 2001 John Wiley & Sons, Ltd.     

Algorithm for analysis of flows in ribbed annuli

Spectral methods for analyses of steady flows in annuli bounded by walls with either axi‐symmetric or longitudinal ribs are developed. The physical boundary conditions are enforced using the immersed boundary conditions concept. In the former case, the Stokes stream function is used to eliminate pressure and to reduce system of field equations to a single fourth‐order partial differential equation. The ribs are assumed to be periodic in the axial direction and this permits representation of the solution in terms of the Fourier expansion. In the latter case, the problem is reduced to the Laplace equations for the flow modifications that can be expressed in terms of the Fourier expansions. The modal functions, which are functions of the radial coordinate, are discretized using Chebyshev polynomials. The problem formulations are closed using either the fixed volume flow rate constraint or the fixed pressure gradient constraint. Various tests have been carried out in order to demonstrate the spectral accuracy of the discretizations, as well as the spectral accuracy of the enforcement of the flow boundary conditions at the ribbed walls using the immersed boundary conditions concept. Special linear solver that takes advantage of the matrix structure has been implemented in order to reduce computational time and memory requirements. It is shown that the algorithm has superior performance when one is interested in the analysis of a large number of geometries, as part of the coefficient matrix that corresponds to the field equation is always the same and one needs to change only the part of the matrix that corresponds to the boundary relations when changing geometry of the flow domain. Copyright © 2011 John Wiley & Sons, Ltd.     

Solving high Reynolds‐number viscous flows by the general BEM and domain decomposition method

In this paper, the domain decomposition method (DDM) and the general boundary element method (GBEM) are applied to solve the laminar viscous flow in a driven square cavity, governed by the exact Navier–Stokes equations. The convergent numerical results at high Reynolds number Re = 7500 are obtained. We find that the DDM can considerably improve the efficiency of the GBEM, and that the combination of the domain decomposition techniques and the parallel computation can further greatly improve the efficiency of the GBEM. This verifies the great potential of the GBEM for strongly non‐linear problems in science and engineering. Copyright © 2004 John Wiley & Sons, Ltd.     

Non‐intrusive reduced‐order modelling of the Navier–Stokes equations based on RBF interpolation

We present a new non‐intrusive model reduction method for the Navier–Stokes equations. The method replaces the traditional approach of projecting the equations onto the reduced space with a radial basis function (RBF) multi‐dimensional interpolation. The main point of this method is to construct a number of multi‐dimensional interpolation functions using the RBF scatter multi‐dimensional interpolation method. The interpolation functions are used to calculate POD coefficients at each time step from POD coefficients at earlier time steps. The advantage of this method is that it does not require modifications to the source code (which would otherwise be very cumbersome), as it is independent of the governing equations of the system. Another advantage of this method is that it avoids the stability problem of POD/Galerkin. The novelty of this work lies in the application of RBF interpolation and POD to construct the reduced‐order model for the Navier–Stokes equations. Another novelty is the verification and validation of numerical examples (a lock exchange problem and a flow past a cylinder problem) using unstructured adaptive finite element ocean model. The results obtained show that CPU times are reduced by several orders of magnitude whilst the accuracy is maintained in comparison with the corresponding high‐fidelity models. Copyright © 2015 John Wiley & Sons, Ltd.     

A boundary element method for steady viscous fluid flow using penalty function formulation

A new boundary element method is presented for steady incompressible flow at moderate and high Reynolds numbers. The whole domain is discretized into a number of eight-noded cells, for each of which the governing boundary integral equation is formulated exclusively in terms of velocities and tractions. The kernels used in this paper are the fundamental solutions of the linearized Navier–Stokes equations with artificial compressibility. Significant attention is given to the numerical evaluation of the integrals over quadratic boundary elements as well as over quadratic quadrilateral volume cells in order to ensure a high accuracy level at high Reynolds numbers. As an illustration, square driven cavity flows are considered for Reynolds numbers up to 1000. Numerical results demonstrate both the high convergence rate, even when using simple (direct) iterations, and the appropriate level of accuracy of the proposed method. Although the method yields a high level of accuracy in the primary vortex region, the secondary vortices are not properly resolved. © 1997 John Wiley & Sons, Ltd.     

A promising boundary element formulation for three‐dimensional viscous flow

In this paper, a new set of boundary‐domain integral equations is derived from the continuity and momentum equations for three‐dimensional viscous flows. The primary variables involved in these integral equations are velocity, traction, and pressure. The final system of equations entering the iteration procedure only involves velocities and tractions as unknowns. In the use of the continuity equation, a complex‐variable technique is used to compute the divergence of velocity for internal points, while the traction‐recovery method is adopted for boundary points. Although the derived equations are valid for steady, unsteady, compressible, and incompressible problems, the numerical implementation is only focused on steady incompressible flows. Two commonly cited numerical examples and one practical pipe flow problem are presented to validate the derived equations. Copyright © 2004 John Wiley & Sons, Ltd.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

A computational methodology for two‐dimensional fluid flows

A weighted residual collocation methodology for simulating two‐dimensional shear‐driven and natural convection flows has been presented. Using a dyadic mesh refinement, the methodology generates a basis and a multiresolution scheme to approximate a fluid flow. To extend the benefits of the dyadic mesh refinement approach to the field of computational fluid dynamics, this article has studied an iterative interpolation scheme for the construction and differentiation of a basis function in a two‐dimensional mesh that is a finite collection of rectangular elements. We have verified that, on a given mesh, the discretization error is controlled by the order of the basis function. The potential of this novel technique has been demonstrated with some representative examples of the Poisson equation. We have also verified the technique with a dynamical core of a two‐dimensional flow in primitive variables. An excellent result has been observed—on resolving a shear layer and on the conservation of the potential and the kinetic energies—with respect to previously reported benchmark simulations. In particular, the shear‐driven simulation at CFL = 2.5 (Courant–Friedrichs–Lewy) and (Reynolds number) exhibits a linear speed up of CPU time with an increase of the time step, Δt. For the natural convection flow, the conversion of the potential energy to the kinetic energy and the conservation of total energy is resolved by the proposed method. The computed streamlines and the velocity fields have been demonstrated. Copyright © 2014 John Wiley & Sons, Ltd.     

A consistent splitting scheme for unsteady incompressible viscous flows I. Dirichlet boundary condition and applications

A well‐recognized approach for handling the incompressibility constraint by operating directly on the discretized Navier–Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the current developments by Guermond and Shen, the possibilities of obtaining accurate pressure and reducing boundary‐layer effect for the pressure are analysed. The present study mainly reports the numerical solutions of an unsteady Navier–Stokes problem based on the so‐called consistent splitting scheme (J. Comput. Phys. 2003; 192 :262–276). At the same time the Dirichlet boundary value conditions are considered. The accuracy of the method is carefully examined against the exact solution for an unsteady flow physics problem in a simply connected domain. The effectiveness is illustrated viz. several computations of 2D double lid‐driven cavity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Local moving least square‐one‐dimensional integrated radial basis function networks technique for incompressible viscous flows

This paper presents a local moving least square‐one‐dimensional integrated radial basis function networks method for solving incompressible viscous flow problems using stream function‐vorticity formulation. In this method, the partition of unity method is employed as a framework to incorporate the moving least square and one‐dimensional integrated radial basis function networks techniques. The major advantages of the proposed method include the following: (i) a banded sparse system matrix which helps reduce the computational cost; (ii) the Kronecker‐ δ property of the constructed shape function which helps impose the essential boundary condition in an exact manner; and (iii) high accuracy and fast convergence rate owing to the use of integration instead of conventional differentiation to construct the local radial basis function approximations. Several examples including two‐dimensional (2D) Poisson problems, lid‐driven cavity flow and flow past a circular cylinder are considered, and the present results are compared with the exact solutions and numerical results from other methods in the literature to demonstrate the attractiveness of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Feedback control of laminar flow separation on NACA23012 airfoil by POD analysis and using perturbed Navier‐Stokes equations

The main purpose of this article is to develop a forced reduced‐order model based on the proper orthogonal decomposition (POD)/Galerkin projection (on isentropic Navier‐Stokes equations) and perturbation method on the compressible Navier‐Stokes equations. The resulting forced reduced‐order model will be used in optimal control of the separated flow over a NACA23012 airfoil at Mach number of 0.2, Reynolds number of 800, and high incidence angle of 24°. The main disadvantage of the POD/Galerkin projection method for control purposes is that controlling parameters do not show up explicitly in the resulting reduced‐order system. The perturbation method and POD/Galerkin projection on the isentropic Navier‐Stokes equations introduce a forced reduced‐order model that can predict the time varying influence of the controlling parameters and the Navier‐Stokes response to external excitations. An optimal control theory based on forced reduced‐order system is used to design a control law for a nonlinear reduced‐order system, which attempts to minimize the vorticity content in the flow field. The test bed is a laminar flow over NACA23012 airfoil actuated by a suction jet at 12% to 18% chord from leading edge and a pair of blowing/suction jets at 15% to 18% and 24% to 30% chord from leading edge, respectively. The results show that wall jet can significantly influence the flow field, remove separation bubbles, and increase the lift coefficient up to 22%, while the perturbation method can predict the flow field in an accurate manner.