{"title":"ADAPTIVE RBF-FD METHOD FOR POISSON’S EQUATION","authors":"J. Slak, G. Kosec","doi":"10.2495/BE420131","DOIUrl":"https://doi.org/10.2495/BE420131","url":null,"abstract":"Solutions to many physical problems governed by partial differential equations (PDEs) often vary significantly in magnitude throughout the problem domain. Although in some special cases the areas with high error are known in advance, in general the error distribution is unknown beforehand. Adaptive techniques for solving PDEs are a standard way of dealing with this problem, where problematic regions are iteratively refined. A step further is automatic adaptivity, where problematic regions are chosen automatically using an error indicator and then refined, until a certain error threshold is reached. In this paper, we apply a recently published technique for automatic adaptivity for strong form meshless methods and solve the Poisson equation and its generalisations, using the popular RBF-FD method. Both 2D and 3D cases are considered, comparing uniform and adaptive refinement, illustrating the advantages of fully automatic adaptivity.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"16 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"133838617","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"ANALYTICAL 3D BOUNDARY ELEMENT IMPLEMENTATION OF FLAT TRIANGLE AND QUADRILATERAL ELEMENTS FOR POTENTIAL AND LINEAR ELASTICITY PROBLEMS","authors":"N. Dumont, Tatiana Galvão Kurz","doi":"10.2495/BE420011","DOIUrl":"https://doi.org/10.2495/BE420011","url":null,"abstract":"This paper introduces a formulation for 3D potential and linear elasticity problems that end up with the analytical handling of all regular, improper, quasi-singular, singular and hypersingular integrals of an implementation using linear triangle (T3) elements. The extension to flat Q4 and T6 elements is almost straightforward. Results at arbitrarily located internal points are also given analytically. The formulation is based on a generalized transformation to subtriangle coordinates that simplifies the problem’s description and enables the adequate interpretation of all relevant geometric features of a discretized boundary segment, so that it becomes possible to arrive at manageable analytical expressions of all integrals. The paper outlines the main concepts and computational features of the proposed formulation, based on an array with all pre-evaluated integrals required in an implementation. An example of 3D potential problems illustrates all particular cases and the most challenging topological configurations one might deal with in practical applications. The procedure may be easily implemented in a general boundary element code, as the usual numerical quadrature schemes for source points sufficiently far from the integration field remain applicable. There is a work in progress for the implementation of the procedure in the frame of a fast multipole algorithm.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"131222367","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"COUPLED FINITE AND BOUNDARY ELEMENT METHOD FOR SOLVING MAGNETIC HYSTERESIS PROBLEMS","authors":"I. Stupakov, M. Royak, N. Kondratyeva","doi":"10.2495/BE420111","DOIUrl":"https://doi.org/10.2495/BE420111","url":null,"abstract":"This paper is devoted to effective solving of magnetic hysteresis problems in 3D modeling of accelerator magnets. It is proposed to use boundary elements and a reduced scalar magnetic potential in the air and finite elements, in addition to a full scalar magnetic potential in ferromagnetic materials. There is no need to mesh the air domain by using that approach. Also, it allows nonlinear and inhomogeneous magnetization in ferromagnets to be taken into account. Numerical simulation is carried out using three different models of hysteresis: the Preisach model, the model based on dry friction and the Vinokurov model. To verify the accuracy of the developed approach, the simulation results are compared with the magnetic field measurements for the deflecting accelerator magnet.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"41 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"115044013","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"EXACT TIME DEPENDENT CALCULATIONS FOR FLUIDS AND SOLIDS","authors":"E. Kansa","doi":"10.2495/BE420151","DOIUrl":"https://doi.org/10.2495/BE420151","url":null,"abstract":"Meshfree transport calculations have been traditionally performed in the Lagrangian or particle velocity mode connected to an underlying distorting mesh. That is remapped (interpolated) onto a regular mesh. This primitive procedure was necessary with extremely primitive electronic computers that has extremely restrictive memory storage and execution speed. However, the meshless radial basis function approach allows points to move; without diffusive re-interpolation; one can choose to allow the motion at each interior point to convert the partial differential equation into an exact differential, thereby eliminating truncation errors. Given the proper tools to utilize extended arithmetic precision on computers, electronic computers can asymptotically approach the power of ideal academic mathematics and solve numerically problems approaching theoretical methods.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"67 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122020651","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"GLOBAL SINGULAR BOUNDARY METHOD FOR SOLVING 2D NAVIER–STOKES EQUATIONS","authors":"J. Mužík, R. Bulko","doi":"10.2495/BE420181","DOIUrl":"https://doi.org/10.2495/BE420181","url":null,"abstract":"This article presents a numerical algorithm based on the singular boundary method (SBM) for the incompressible steady Navier–Stokes equations formulated using primitive variables. The SBM with the Stokeslet fundamental solution and dual reciprocity (DR) principle has been chosen to solve the nonlinear flow equations. The particular solution of the non-homogeneous Stokes equations is constructed as a linear combination of implicitly local radial basis function. The simple direct iterative scheme was used to handle nonlinearities of Navier–Stokes equations with a variation of the non-homogeneous term of the particular solution. The non-homogeneous term is formed using the nonlinear convective terms of the momentum equations, evaluated using values from previous iterations. It is found that SBM with a localized DR principle gives reasonable results for numerical problems of lid-driven cavity flow up to Re = 3200, and the backward-facing step at Re = 800.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"1 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130119691","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"IMMERSED BOUNDARY METHOD APPLICATION AS A WAY TO BUILD A SIMPLIFIED FLUID-STRUCTURE MODEL","authors":"J. Borges, M. Lourenço, E. Padilla, C. Micallef","doi":"10.2495/BE420221","DOIUrl":"https://doi.org/10.2495/BE420221","url":null,"abstract":"A simplified fluid-structure interaction model, consisting of a cylinder tethered by a spring system interacting dynamically with an incompressible two-dimensional lid-driven cavity flow, is solved using the immersed boundary method. Results show that when the spring forces are weaker than the fluid drag force, the springs stretch freely and the cylinder motion is the direct result of the fluid dynamics action. For higher values of spring forces, the cylinder motion reaches a maximum displacement, and the spring forces induce the cylinder to an oscillatory movement damped by the fluid drag forces. Subsequently the amplitude of the displacement decreases. The cylinder motion is restricted within the mainstream fluid flow, where the maximum displacement reduces as the Reynolds number increases.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"30 2","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"120900588","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"MESH-FREE ANALYSIS OF PLATE BENDING PROBLEMS BY MOVING FINITE ELEMENT APPROXIMATION","authors":"V. Sládek, J. Sládek, M. Repka","doi":"10.2495/BE420191","DOIUrl":"https://doi.org/10.2495/BE420191","url":null,"abstract":"","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"11 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"126561511","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"3D CAUCHY PROBLEM FOR AN ELASTIC LAYER: INTERFACIAL CRACKS DETECTION","authors":"A. Galybin, S. Aizikovich","doi":"10.2495/BE420251","DOIUrl":"https://doi.org/10.2495/BE420251","url":null,"abstract":"This study presents a Cauchy-type problem of 3D elasticity for an elastic layer that can be bonded to an infinite base (half-space) made of dissimilar elastic material. The initial conditions are given on one side of the layer and both stress and displacement vectors are assumed to be known simultaneously. No conditions are specified on the other side. In the case of this side being fully bonded to the base, the stress and displacement vectors are continuous across the interface. This fact introduces certain relationships that have to be imposed on the initial conditions in order to obey continuity. We use these in order to detect a possible appearance of delamination of the interface. By using the double Fourier transform and the general solution of 3D elasticity in terms of harmonic functions, the initial value problem is reduced to a system of Fredholm integral equations of the first kind. Solutions of such systems are usually unstable; therefore, a numerical approach is suggested to overcome this difficulty by using the SVD regularisation. A possibility of delamination detection is discussed.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"35 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132625906","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"NEW BOUNDARY ELEMENT FORMULATION FOR THE SOLUTION OF LAPLACE’S EQUATION","authors":"J. Lobry","doi":"10.2495/BE420061","DOIUrl":"https://doi.org/10.2495/BE420061","url":null,"abstract":"The boundary element method (BEM) is typically used for solving potential problems and has several advantages over the traditional finite element method (FEM). However, the normal derivative of the potential appears explicitly as an unknown, with some inherent ambiguity at the corner nodes. Moreover, the BEM requires time consuming numerical integration (Gaussian quadrature) along the boundary of the domain. In this paper, we propose a new approach of BEM by introducing the weak nodal “cap” flux approach defined in the context of the finite element method. The domain integrals are eliminated at the discrete level by introducing the FEM approximation of the fundamental solutions at every node of the related mesh as basic functions in the Galerkin formulation of the BVP under study. The implementation of this new technique appears to be simpler as no numerical integration on the boundary of the domain is required so that the method leads to a substantially reduced computational burden. Our method is compared to the classical BEM for the numerical solution of the two-dimensional Laplace equation. It is observed that the normal flux presents a better behaviour at corners. A loss of accuracy may occur but it is compensated by a smaller execution time, allowing a finer mesh.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"49 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"134457644","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"NUMERICAL SCHEMES FOR VORTEX SHEET INTENSITY COMPUTATION IN VORTEX METHODS TAKING INTO ACCOUNT THE CURVILINEARITY OF THE AIRFOIL SURFACE LINE","authors":"K. Kuzmina, I. Marchevsky","doi":"10.2495/BE420241","DOIUrl":"https://doi.org/10.2495/BE420241","url":null,"abstract":"In vortex methods, vorticity is the primary computed variable. The problem of the accuracy improvement of vorticity generation simulation at the airfoil surface line in 2D vortex methods is considered. The generated vorticity is simulated by a thin vortex sheet at the airfoil surface line, and it is necessary to determine the intensity of this sheet at each time step. It can be found from the no-slip boundary condition, which leads to a vector boundary integral equation. There are two approaches to satisfy this equation: the first one leads to a singular integral equation of the 1st kind, while the second one leads to a Fredholm-type integral equation of the 2nd kind with bounded kernel for smooth airfoils. Usually, for numerical solution of the boundary integral equation, the airfoil surface line is replaced by a polygon, which consists of straight segments (panels). A discrete analogue of the integral equation can be obtained using the Galerkin method. Different families of basis and projection functions lead to numerical schemes with different complexity and accuracy. For example, a numerical scheme with piecewise-constant basis functions provides the first order of accuracy for vortex sheet intensity, and a numerical scheme with piecewise-linear functions gives the second order of accuracy. However, the velocity field near the airfoil surface line is also of interest. In the case of rectilinear airfoil surface line discretization, the accuracy of velocity field reconstruction has no more than the first order of accuracy for both, piecewise-constant and piecewise-linear numerical schemes. In order to obtain a higher order of accuracy for velocity field reconstruction, it is necessary to take into account the curvilinearity of the airfoil surface line. In this research, we have developed such an approach, which provides the second order of accuracy both, for vortex sheet intensity computation and velocity field reconstruction.","PeriodicalId":429597,"journal":{"name":"Boundary Elements and other Mesh Reduction Methods XLII","volume":"127 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2019-09-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132146910","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}