Adaptive Numerical Modeling of Engineering Problems Using Hierarchical Fup Basis Functions and Control Volume IsoGeometric Analysis

2021 ◽  
Author(s):  
◽  
Grgo Kamber ◽  

The main objective of this thesis is to utilize the powerful approximation properties of Fup basis functions for numerical solutions of engineering problems with highly localized steep gradients while controlling spurious numerical oscillations and describing different spatial scales. The concept of isogeometric analysis (IGA) is presented as a unified framework for multiscale representation of the geometry and solution. This fundamentally high-order approach enables the description of all fields as continuous and smooth functions by using a linear combination of spline basis functions. Classical IGA usually employs Galerkin or collocation approach using B-splines or NURBS as basis functions. However, in this thesis, a third concept in the form of control volume isogeometric analysis (CV-IGA) is used with Fup basis functions which represent infinitely smooth splines. Novel hierarchical Fup (HF) basis functions is constructed, enabling a local hp-refinement such that they can replace certain basis functions at one resolution level with new basis functions at the next resolution level that have a smaller length of the compact support (h-refinement), but also higher order (p-refinement). This hp-refinement property enables spectral convergence which is significant improvement in comparison to the hierarchical truncated B-splines which enable h-refinement and polynomial convergence. Thus, in domain zones with larger gradients, the algorithm uses smaller local spatial scales, while in other region, larger spatial scales are used, controlling the numerical error by the prescribed accuracy. The efficiency and accuracy of the adaptive algorithm is verified with some classic 1D and 2D benchmark test cases with application to the engineering problems with highly localized steep gradients and advection-dominated problems.

Author(s):  
Mohammad Tamsir ◽  
Neeraj Dhiman ◽  
F.S. Gill ◽  
Robin

This paper presents an approximate solution of 3D convection diffusion equation (CDE) using DQM based on modified cubic trigonometric B-spline (CTB) basis functions. The DQM based on CTB basis functions are used to integrate the derivatives of space variables which transformed the CDE into the system of first order ODEs. The resultant system of ODEs is solved using SSPRK (5,4) method. The solutions are approximated numerically and also presented graphically. The accuracy and efficiency of the method is validated by comparing the solutions with existing numerical solutions. The stability analysis of the method is also carried out.


Water ◽  
2018 ◽  
Vol 10 (12) ◽  
pp. 1787 ◽  
Author(s):  
Luka Malenica ◽  
Hrvoje Gotovac ◽  
Grgo Kamber ◽  
Srdjan Simunovic ◽  
Srikanth Allu ◽  
...  

A novel numerical model for groundwater flow in karst aquifers is presented. A discrete-continuum (hybrid) approach, in which a three-dimensional matrix flow is coupled with a one-dimensional conduit flow, was used. The laminar flow in the karst matrix is described by a variably saturated flow equation to account for important hydrodynamic effects in both the saturated and unsaturated zones. Turbulent conduit flow for both free surface and pressurized flow conditions was captured via the noninertia wave equation, whereas the coupling of two flow domains was established through an exchange term proportional to head differences. The novel numerical approach based on Fup basis functions and control-volume formulation enabled us to obtain smooth and locally conservative numerical solutions. Due to its similarity to the isogeometric analysis concept (IGA), we labeled it as control-volume isogeometric analysis (CV-IGA). Since realistic verification of the karst flow models is an extremely difficult task, the particular contribution of this work is the construction of a specially designed 3D physical model ( dimensions: 5.66 × 2.95 × 2.00 m) in order to verify the developed numerical model under controlled laboratory conditions. Heterogeneous porous material was used to simulate the karst matrix, and perforated pipes were used as karst conduits. The model was able to capture many flow characteristics, such as the interaction between the matrix and conduit, rainfall infiltration through the unsaturated zone, direct recharge through sinkholes, and both free surface and pressurized flow in conduits. Two different flow experiments are presented, and comparison with numerical results confirmed the validity of the developed karst flow model under complex laboratory conditions.


2020 ◽  
Vol 92 (10) ◽  
pp. 1437-1461 ◽  
Author(s):  
Grgo Kamber ◽  
Hrvoje Gotovac ◽  
Vedrana Kozulić ◽  
Luka Malenica ◽  
Blaž Gotovac

2013 ◽  
Vol 2013 ◽  
pp. 1-9
Author(s):  
Xinkang Li ◽  
Jifa Zhang ◽  
Yao Zheng

Isogeometric analysis (IGA) based on nonuniform rational B-splines (NURBS) is developed for static analysis of beams and plates using the third-order shear deformation theory (TSDT). TSDT requires C1-continuity of generalized displacements; quadratic, cubic, and quartic NURBS basis functions are utilized to satisfy this requirement. Due to the noninterpolatory nature of NURBS basis functions, a penalty method is presented to enforce the essential boundary conditions. A series of numerical examples of thick and thin beams and plates with different boundary conditions are presented. Results are compared with analytical solutions and other numerical methods. First-order shear deformation theory (FSDT) is also employed and compared with TSDT in the stress analysis. The effects of aspect ratios and boundary conditions are discussed as well.


Author(s):  
Prithvi Bhat ◽  
Zeger Bontinck ◽  
Jacopo Corno ◽  
Sebastian Schöps ◽  
Herbert DeGersem

Purpose This paper aims to propose the use of isogeometric analysis (IGA) for the simulation of electrical machines to represent their geometries exactly and obtain numerical solutions of high accuracy and regularity. Design/methodology/approach IGA makes use of non-uniform rational b-splines to parametrise the domain and approximate the solution spaces. Dealing with the different stator and rotor topologies, the computational domain is split into two non-overlapping parts on which Maxwell’s equations are solved independently and are interconnected by a classical Schwarz domain decomposition scheme. The results are compared with the conventional polynomial finite element method (FEM). Findings The new methodology is reliable and efficient. The obtained solutions of the fields are in good agreement with the ones obtained by the FEM approach. IGA offers a better accuracy than FEM. Originality/value The application of IGA combined with domain decomposition to the model of an electric machine is a new and original contribution.


2021 ◽  
Vol 381 ◽  
pp. 113779
Author(s):  
Wenbin Hou ◽  
Kai Jiang ◽  
Xuefeng Zhu ◽  
Yuanxing Shen ◽  
Ping Hu

2014 ◽  
Vol 986-987 ◽  
pp. 1418-1421
Author(s):  
Jun Shan Li

In this paper, we propose a meshless method for solving the mathematical model concerning the leakage problem when the pressure is tested in the gas pipeline. The method of radial basis function (RBF) can be used for solving partial differential equation by writing the solution in the form of linear combination of radius basis functions, that is, when integrating the definite conditions, one can find the combination coefficients and then the numerical solution. The leak problem is a kind of inverse problem that is focused by many engineers or mathematical researchers. The strength of the leak can find easily by the additional conditions and the numerical solutions.


2017 ◽  
Vol 21 (3) ◽  
pp. 835-866 ◽  
Author(s):  
Meng Wu ◽  
Bernard Mourrain ◽  
André Galligo ◽  
Boniface Nkonga

AbstractMotivated by the magneto hydrodynamic (MHD) simulation for Tokamaks with Isogeometric analysis, we present splines defined over a rectangular mesh with a complex topological structure, i.e., with extraordinary vertices. These splines are piecewise polynomial functions of bi-degree (d,d) and parameter continuity. And we compute their dimension and exhibit basis functions called Hermite bases for bicubic spline spaces. We investigate their potential applications for solving partial differential equations (PDEs) over a physical domain in the framework of Isogeometric analysis. For instance, we analyze the property of approximation of these spline spaces for the L2-norm; we show that the optimal approximation order and numerical convergence rates are reached by setting a proper parameterization, although the fact that the basis functions are singular at extraordinary vertices.


2016 ◽  
Vol 311 ◽  
pp. 503-536 ◽  
Author(s):  
Xuefeng Zhu ◽  
Ping Hu ◽  
Zheng-Dong Ma

Sign in / Sign up

Export Citation Format

Share Document