Algorithms for construction of Element Partition Trees for Direct Solver executed over h refined grids with B-splines and C 0 separators
|
|
- Claud Atkins
- 6 years ago
- Views:
Transcription
1 Available online at ScienceDirect Procedia Computer Science 108C (2017) International Conference on Computational Science, ICCS 2017, June 2017, Zurich, Switzerland Algorithms for construction of Element Partition Trees for Direct Solver executed over h refined grids with B-splines and C 0 separators Bartosz Janota 1, Maciej Paszyński 1 AGH University of Science and Technology, Kraków, Poland Abstract We propose a way of performing isogeometric finite element method (IGA-FEM) computations over h refined grids with B-spline basis functions. Namely, we propose to use the B-spline basis functions defined over patches of elements with C 0 separators between the refinement levels. Our definition of the B-splines and C 0 separators allows introduction of arbitrary order B-splines over 2D grids refined towards singularities. We also present an algorithm for construction of element partition trees (EPT) over h refined grids with modified B-splines. The EPT allows generating an ordering which gives a linear computational cost of the multi-frontal solver over 2D grids refined towards a point or an edge singularity. We present the algorithm for transforming the EPT into an ordering. We also verify the linear computational cost of the proposed method on grids with point and edge singularity. We compare our method to h-adaptive finite element method (h-fem) computations with Lagrange polynomials The Authors. Published by Elsevier B.V. Peer-review under responsibility of the scientific committee of the International Conference on Computational Science Keywords: multi-frontal direct solver, element partition trees, b-spline, ordering 1 Introduction The main idea of the isogeometric analysis [10] (isogeometric finite element method IGA-FEM) is to connect computational methods with computer aided design tools, such as AutoCAD. Before the IGA-FEM has been invented, the geometry of the computational domain had often been defined in Computer Aided Design (CAD) tools using B-splines or their rationalized version called Non-Uniform Rational B-splines (NURBS). The computations, however, had been performed in finite element method (FEM) tools, using hierarchical Lagrange polynomials, and so the mesh and the basis functions were converted from the CAD world to the FEM world. Since the invention of the IGA-FEM, both steps are integrated. The IGA-FEM has multiple applications in shear deformable shell theory [5], phase field models [18, 19], phase-separation simulations with application to cancer growth simulations [21, 22], wind turbine aerodynamics [24], incompressible hyper-elasticity [20], turbulent flow The Authors. Published by Elsevier B.V. Peer-review under responsibility of the scientific committee of the International Conference on Computational Science /j.procs
2 858 Bartosz Janota et al. / Procedia Computer Science 108C (2017) simulations [8], transport of drugs in cardiovascular applications [23] or the blood flow simulations itself [3, 2, 7]. We introduce rationalized B-splines and C 0 separators between same refinement levels elements, which allows extending the IGA-FEM with B-splines into grids refined toward a point or an edge singularities. Alternative approach is the application of T-splines [28], which requires the introduction of the T-junction extensions to obtain the analysis suitable grids, but in this paper we focus on the B-spline basis functions. We also propose an algorithm for construction of orderings to be used by the multi-frontal direct solvers, based on the concepts introduced in [15, 16] for standard finite element method. When comparing classical FEM and IGA-FEM on regular grids, the classical FEM is only C 0 on the borders of elements, while C (p 1) inside particular elements, whereas IGA-FEM ensures C p 1 regular solutions if one uses B-splines of the order p. Of course, there are some limitations, such as the IGA-FEM solution increases the computational cost, as described in [11]. The application of B-spline basis functions results in O(N 1.5 p 3 ) computational complexity on homogeneous 2D meshes and O(N 2 p 3 ) on homogeneous 3D meshes, where N is the number of unknowns over the regular mesh and p is the B-spline order. The computational complexity is p 3 times higher in comparison to standard FEM with Lagrange polynomials. The computational complexity can be reduced down to O(Np 2 )in2d or O(N 4 /3p 2 )in3d, by using shared-memory [13] or distributed memory solvers [14]. In the second case, the communication complexity is also of the same order as the computational complexity. However, when we switch into the computational costs, the constants in front of the computational complexity come into the picture. As it is shown in [17], the computational cost of IGA-FEM for the fixed number of elements becomes lower than the cost of FEM. However, when we utilize patches of higher order elements with C 0 separators between them, the resulting method (called riga) is cheaper than both IGA-FEM and FEM. In [17] the higher order elements are grouped into patches over the regular grids. In this paper we rather focus on adaptive, h-refined grids. In this paper, we propose to identify patches of elements with same refinement levels and introduce C 0 separators between them. 2 B-splines over refined meshes Our idea of the introduction of B-splines over the 2D mesh with singularities is to introduce C 0 separators between the refinement levels surrounding the singularities and to rationalize the B-splines. The rationlization is necessary to fulfill the partition of unity property. In this case, we only have C 0 continuity between the refinement levels. However, we still have higher order continuity inside the patches of elements, located inside the refinement levels. In the following considerations, for simplicity, we will focus on two-dimensional meshes with point and edge singularity, presented in Figure 1. We truncate the B-spline basis functions with the C 0 separators defined over the mesh, and rationalize them, to obtain the partition of unity feature Ri k(x) = B n i,k(x)ω i j=1 B j,k(x)ω j. The drawback of introductions the C 0 sepparators is that some of the supports of basis functions have rectangular shape, and some other are non-rectangular ones. In Figures 2 and 3 we can see some examples of basis functions spanned over the 3 rd depth level 2D mesh. In the case of non-rectangular supports, the green elements belong to the support of the basis functions located on the yellow point. The same situation exists when we want to span basis functions from the blue points, except the most center ones. Each non-rectangular basis function is built from two pieces, as it is illustrated in Figure 4. Red basis function is a 3/4 of the regular basis function and this is the first piece. The second one is
3 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Figure 1: 2D mesh, refined towards the point singularity (a) Regular 4x4 mesh, 1 st level depth, no separators (b) Regular 4x4 mesh, 4 th level depth, double lines yet mean C 0 separators Figure 2: 3 rd depth level 2D mesh with some sample basis functions combined from two regular basis functions (dark and light blue) and the black one. The black basis function is a sum of two regular basis functions and one linear function. Thanks to them, we have continuity on the edges between red basis function and all the others and a zero value in the most central corner of that basis function. The most important thing is that the red value between two blue basis functions is zero so that basis function perfectly fits into 3. 3 Element Partition Tress Element Partition Tree (EPT) can be built based on the mesh, as shown in Section (2), and the basis functions spanned over that mesh. EPT is an essential base for the ordering (see Section 3.1). Ordering resulting from EPT may turn the not-nice matrix into the one that is multi-diagonal or almost multi-diagonal, and the resulting cost of the factorization will be linear or quasi-linear. On Algorithm 1 you can see EPT generation algorithm. The algorithm cuts layer by layer, getting deeper into the mesh. Thanks to that, we are sure, that the new sub-tree (created from cut part) is not too deep. We create those sub-trees recursively. As you can see on Figure 5, EPT is deep and unbalanced. EPT can be also generated for other refined grids, like for the edge singularity case. EPT is going to be used as a base for ordering.
4 860 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Figure 3: Non-rectangular support Figure 4: Non-rectangular support basis function Figure 5: Complete EPT 4
5 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Algorithm 1 EPT generation procedure GenerateEptNodes(rootNode) for i =1todepth do There is an example I refer to, see 5 ProcessEptNode(left) red on 5 ProcessEptNode(right) blue on 5 ProcessEptNode(top) green on 5 ProcessEptNode(bottom) purple on 5 end for ProcessEptNode(innerMost16Nodes) orange on 5 end procedure procedure ProcessEptNode(node) AddEptNode(node) if node can be split then ProcessEptNode(nodeF irsthalf) ProcessEptNode(nodeSecondHalf) end if end procedure procedure AddEptNode(node) EPT EPT + node end procedure 3.1 Ordering Figure 6: Original FEM matrix for depth =8 Figure 6 presents the original FEM matrix A generated for the regular mesh with point singularity (depth = 8) with the natural ordering of elements. This matrix is not multi-diagonal at all (it is sparse but the non-zero entries are spread around the whole matrix and it has some blocks of non-zero entries far from the diagonal, which results in a large number of additional non-zero entries generated during the factorization), so the cost of the elimination is non-linear.
6 862 Bartosz Janota et al. / Procedia Computer Science 108C (2017) We want to reduce the bandwidth as much as possible and to achieve that we need to permute the rows somehow to make numerical factorization more efficient. We can do that thanks to that very basic matrix transformation: B = P 1 AP (1) where P is the permutation matrix based on elimination ordering and B is the outcome matrix that gives as the linear cost of the multi-frontal computations. The permutation matrix P and the ordering can be obtained by post-order traversal of EPT. This is illustrated in Figure 7 for the exemplary EPT for 1/4 of 5. The ordering is generated by Algorithm 2. If we meet any B-spline that has not been added yet, we should add it to the ordering list. Main of that algorithm is to eliminate each of the basis functions as deep in the EPT as it is possible, so this is why we eliminate 1, 2, 3, 4, 5, 6 and 7 in the leaves, and not earlier. 9 and 11 could be eliminated in the root node, but they exist deeper as well, so we add them to the ordering list once we meet them there. There are some basis functions which do not exist deeper in the EPT than in the root, so, in the last step we add 8, 10 and 12 to the ordering list. Figure 7: EPT ordering, EPT for 1/4 of 5, last step of 2 Figure 8: 8 th, 12 th and 10 th basis function added, we skip all of the remaining basis functions, because they are already added, final ordering list: [1, 2, 9, 3, 4, 5, 7, 6, 11, 8, 12, 10]
7 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Algorithm 2 Ordering generation procedure TraverseEptPostOrder(currentNode) TraverseEptPostOrder(leftSon) traverse the left subtree TraverseEptPostOrder(rightSon) traverse the right subtree if orderinglist does not contain any B-splines spanned over currentnode then orderinglist orderinglist + newly encountered B-splines end if end procedure Our matrix A after the application of the generated ordering is presented in Figure 9. We can see that the matrix now is multiple-diagonal, and the cost of its factorization will be linear. Figure 9: Multi-diagonal matrix, after permutation of 6 with a use of ordering (7) 4 Results evaluation Our theoretical estimation of computational cost can be proved with the results shown below. To do that, we measured multi-frontal direct solvers computational cost and we express that cost in number of floating-point operations (FLOPS). We are using ordering (obtained from EPT) to control MUMPS solver. We run our experiment for the two-dimensional mesh refined towards point and edge singularity. We compare the mesh with B-spline basis functions and C 0 separators with the mesh with Lagrange polynomials. In the first case, we generate the ordering by using the EPT and submit it to MUMPS solver. In the second case we use the automatic ordering proposed by MUMPS, and the grid obtained from hp2d code. We consider the square grids and 2 nd order B-splines.
8 864 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Figure 10: FLOPS estimate for point singularity with respect to the degrees of freedom (log scale), for B-splines with C 0 separators, ordering from EPT, passed to MUMPS solver, versus Lagrange polynomials and MUMPS with automatic ordering Figure 11: FLOPS estimate for edge singularity with respect to the degrees of freedom (log scale), for B-splines with C 0 separators, ordering from EPT, passed to MUMPS solver, versus Lagrange polynomials and MUMPS with automatic ordering The experiments are presented in Figures 10 and 11. For the case of point singularity, both costs are linear with respect to the number of degrees of freedom, and the costs are almost identical. For the case of edge singularity, the B-splines with C 0 separators with EPT ordering give linear cost, while Lagrange polynomials with MUMPS with automatic ordering gives quasi-linear ordering.
9 Bartosz Janota et al. / Procedia Computer Science 108C (2017) Conclusions In this paper we have shown that introducing C 0 separators between refinement levels and modifying B-splines by rationalizing and truncating them with the separators, results in a linear computational cost reduction when compared to classical FEM with Lagrange polynomials. Namely, for point singularity the costs are comparable, and for the edge singularity in 2D the B-splines give linear cost while Lagrange polynomials give quasi-linear cost. We proved that rationalized B-splines span over patches of elements separated with C 0 separators are the keys. They can be used as the basis functions for grid refinements both for point and edge singularity. For the generation of the ordering with our B-splines and C 0 separators we performed a postorder traversal of the element partition tree. We showed that our approach allows to reduce the computational cost when comparing to classical FEM with Lagrange polynomials. The future work will involve implementation in three-dimensions [25] and parallel computations in two and three-dimensions [26, 27]. 6 Acknowledgements This work is supported by National Science Centre, Poland grant no. 2016/21/B/ST6/ References [1] Y. Bazilevs, L. Beirao da Veiga, J.A. Cottrell, T.J.R. Hughes, and G. Sangalli, Isogeometric analysis: Approximation, stability and error estimates for h-refined meshes, Mathematical Methods and Models in Applied Sciences, 16 (2006) [2] Y. Bazilevs, V.M. Calo, J.A. Cottrell, T.J.R. Hughes, A. Reali, G. Scovazzi, Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows, Computer Methods in Applied Mechanics and Engineering 197 (2007) [3] Y. Bazilevs, V.M. Calo, Y. Zhang, T.J.R. Hughes: Isogeometric fluid-structure interaction analysis with applications to arterial blood flow, Computational Mechanics 38 (2006). [4] Y. Bazilevs, V. M. Calo, J. A. Cottrell, J. A. Evans, S. Lipton, M. A. Scott, T. W. Sederberg, Isogeometric analysis using T-splines, Computer Methods in Applied Mechanics and Engineering, 199 (2010) [5] D.J. Benson, Y. Bazilevs, M.-C. Hsu and T.J.R. Hughes, A large-deformation, rotation-free isogeometric shell, Computer Methods in Applied Mechanics and Engineering, 200 (2011) [6] M. Bubak, J. Kitowski, K. Wiatr, EScience on Distributed Computing Infrastructure: Achievements of PLGrid Plus Domain-specific Services and Tools, vol (2014) Springer [7] V.M. Calo, N. Brasher, Y. Bazilevs, T.J.R. Hughes, Multiphysics Model for Blood Flow and Drug Transport with Application to Patient-Specific Coronary Artery Flow, Computational Mechanics, 43(1) (2008) [8] K. Chang, T.J.R. Hughes, V.M. Calo, Isogeometric variational multiscale large-eddy simulation of fully-developed turbulent flow over a wavy wall, Computers and Fluids, 68 (2012) [9] B. Janota, P. Lipski, and M. Paszyński, Computational complexity of isogeometric analysis with T-splines and B-splines over 2D grids refined towards singularities, Procedia Computer Science, proceedings of International Conference on Computational Science ICCS, (2016). [10] J. Cottrell, and T. Hughes, and Y. Bazilevs, Isogeometric Analysis: Toward Integration of CAD and FEA, (2009).
10 866 Bartosz Janota et al. / Procedia Computer Science 108C (2017) [11] N. Collier, D. Pardo, L. Dalcin, M. Paszyński, and V. Calo, The cost of continuity: A study of the performance of isogeometric finite elements using direct solvers, Computer Methods in Applied Mechanics and Engineering, (2012), 213, [12] L. Piegl, and W. Tiller, The NURBS Book (Second Edition), Springer-Verlag New York, Inc., (1997). [13] M. Woźniak, K. Kuźnik, M. Paszyński, V. Calo, and D. Pardo, Computational cost estimates for parallel shared memory isogeometric multi-frontal solvers, Computers & Mathematics with Applications, 67(10) (2014) [14] M. Woźniak, M. Paszyński, D. Pardo, L. Dalcin, V. Calo, Computational cost of isogeometric multi-frontal solvers on parallel distributed memory machines, Computers Methods in Applied Mechanics and Engineering, 284 (2015) [15] H. AbouEisha, M. Moshkov, V. Calo, M. Paszyński, D. Goik, and K. Jopek, Dynamic Programming Algorithm for Generation of Optimal Elimination Trees for Multi-frontal Direct Solver Over H-refined Grids, Procedia Computer Science, 29 (2014) 29, [16] M. Paszyński, Fast Solvers for Mesh-Based Computations, CRC Press, Taylor & Francis, (2016). [17] D. Garcia, D. Pardo, L. Dalcin, M. Paszyński, N. Collier, V. Calo, The value of continuity: Refined isogeometric analysis and fast direct solvers, Computer Methods in Applied Mechanics and Engineering (2017) DOI: /j.cma [18] L. Dedè,T.J.R. Hughes, S. Lipton, V.M. Calo, Structural topology optimization with isogeometric analysis in a phase field approach, USNCTAM2010, 16th US National Congree of Theoretical and Applied Mechanics. [19] L. Dedè, M. J. Borden, T.J.R. Hughes, Isogeometric analysis for topology optimization with a phase field model, ICES REPORT 11-29, The Institute for Computational Engineering and Sciences, The University of Texas at Austin (2011). [20] R. Duddu, L. Lavier, T.J.R. Hughes, V.M. Calo, A finite strain Eulerian formulation for compressible and nearly incompressible hyper-elasticity using high-order NURBS elements, International Journal of Numerical Methods in Engineering, 89(6) (2012) [21] H. Gómez, V.M. Calo, Y. Bazilevs, T.J.R. Hughes, Isogeometric analysis of the Cahn-Hilliard phase-field model, Computer Methods in Applied Mechanics and Engineering 197 (2008) [22] H. Gómez, T.J.R. Hughes, X. Nogueira, V.M. Calo, Isogeometric analysis of the isothermal Navier- Stokes-Korteweg equations. Computer Methods in Applied Mechanics and Engineering 199 (2010) [23] S. Hossain, S.F.A. Hossainy, Y. Bazilevs, V.M. Calo, T.J.R. Hughes, Mathematical modeling of coupled drug and drug-encapsulated nanoparticle transport in patient-specific coronary artery walls, Computational Mechanics, doi: /s , (2011). [24] M.-C. Hsu, I. Akkerman, Y. Bazilevs, High-performance computing of wind turbine aerodynamics using isogeometric analysis, Computers and Fluids, 49(1) (2011) [25] A. Paszyńska, E. Grabska, M. Paszyński, A Graph Grammar Model of the hp Adaptive Three Dimensional Finite Element Method. Part II., Fundamenta Informaticae, 114(2) (2012) [26] M. Paszyński, A. Paszyńska, Graph transformations for modeling parallel hp-adaptive Finite Element Method, Lecture Notes in Computer Science, 4967 (2008) [27] M. Paszyński, On the Parallelization of Self-Adaptive hp-finite Element Methods Part I. Composite Programmable Graph Grammar Model, Fundamenta Informaticae, 93(4) (2009) [28] L. Beirao da Veiga, A. Buffa, G. Sangalli, R. Vazquez, Analysis-suitable T-splines of arbitrary degree: definition and properties, Mathematical Models and Methods in Applied Sciences, 23(11) 2013
Fast solvers for mesh-based computations
Fast solvers for mesh-based computations Maciej Paszyński Department of Computer Science AGH University, Krakow, Poland home.agh.edu.pl/paszynsk Collaborators: Anna Paszyńska (Jagiellonian University,Poland)
More informationIsogeometric Analysis of Fluid-Structure Interaction
Isogeometric Analysis of Fluid-Structure Interaction Y. Bazilevs, V.M. Calo, T.J.R. Hughes Institute for Computational Engineering and Sciences, The University of Texas at Austin, USA e-mail: {bazily,victor,hughes}@ices.utexas.edu
More informationOptimization of Execution Time, Energy Consumption, and Accuracy during Finite Element Method Simulations
Optimization of Execution Time, Energy Consumption, and Accuracy during Finite Element Method Simulations Maciej Woźniak, Marcin Łoś, Leszek Siwik, Darusz Król, Maciej Paszyński Department of Computer
More informationHYBRID DIRECT AND ITERATIVE SOLVER FOR H ADAPTIVE MESHES WITH POINT SINGULARITIES
HYBRID DIRECT AND ITERATIVE SOLVER FOR H ADAPTIVE MESHES WITH POINT SINGULARITIES Maciej Paszyński Department of Computer Science, AGH University of Science and Technology, Kraków, Poland email: paszynsk@agh.edu.pl
More informationAPPLICATION OF PROJECTION-BASED INTERPOLATION ALGORITHM FOR NON-STATIONARY PROBLEM
Computer Science 7 (3) 26 http://dx.doi.org/.7494/csci.26.7.3.297 Maciej Woźniak Maciej Paszyński APPLICATION OF PROJECTION-BASED INTERPOLATION ALGORITHM FOR NON-STATIONARY PROBLEM Abstract In this paper,
More informationHypergraph Grammars in hp-adaptive Finite Element Method
Available online at www.sciencedirect.com Procedia Computer Science 18 (2013 ) 1545 1554 2013 International Conference on Computational Science Hypergraph Grammars in hp-adaptive Finite Element Method
More informationThe NURBS and GeoPDEs packages Octave software for research on IGA
The NURBS and GeoPDEs packages Octave software for research on IGA Rafael Vázquez IMATI Enrico Magenes, Pavia Consiglio Nazionale delle Ricerche R. Vázquez (IMATI-CNR Italy) The NURBS and GeoPDEs packages
More informationAN ISOGEOMETRIC REISSNER-MINDLIN SHELL WITH LAGRANGE BASIS
11th World Congress on Computational Mechanics (WCCM XI) 5th European Conference on Computational Mechanics (ECCM V) 6th European Conference on Computational Fluid Dynamics (ECFD VI) E. Oñate, J. Oliver
More informationIsogeometric Collocation Method
Chair for Computational Analysis of Technical Systems Faculty of Mechanical Engineering, RWTH Aachen University Isogeometric Collocation Method Seminararbeit By Marko Blatzheim Supervisors: Dr. Stefanie
More informationA direct solver with reutilization of previously-computed LU factorizations for h-adaptive finite element grids with point singularities
Maciej Paszynski AGH University of Science and Technology Faculty of Computer Science, Electronics and Telecommunication Department of Computer Science al. A. Mickiewicza 30 30-059 Krakow, Poland tel.
More informationICES REPORT September M.A. Scott, T.J.R. Hughes, T.W. Sederberg, M.T. Sederberg
ICES REPORT 14-33 September 2014 An integrated approach to engineering design and analysis using the Autodesk T-spline plugin for Rhino3d by M.A. Scott, T.J.R. Hughes, T.W. Sederberg, M.T. Sederberg The
More informationRecent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA
Recent Developments in Isogeometric Analysis with Solid Elements in LS-DYNA Liping Li David Benson Attila Nagy Livermore Software Technology Corporation, Livermore, CA, USA Mattia Montanari Nik Petrinic
More informationA NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION OF FLOW DOMAINS
6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 11 15 June 2018, Glasgow, UK A NURBS-BASED APPROACH FOR SHAPE AND TOPOLOGY OPTIMIZATION
More informationInteractive modelling and simulation using blending techniques. Aleksander Pedersen
Interactive modelling and simulation using blending techniques Aleksander Pedersen NUC R&D Group in mathematical and geometrical modeling, numerical simulations, programming and visualization, Narvik University
More informationParallel multi-frontal solver for isogeometric finite element methods on GPU
Parallel multi-frontal solver for isogeometric finite element methods on GPU Maciej Paszyński Department of Computer Science, AGH University of Science and Technology, Kraków, Poland email: paszynsk@agh.edu.pl
More informationFrequency Domain Analysis for Isogeometric Element in LS-DYNA
Frequency Domain Analysis for Isogeometric Element in LS-DYNA Liping Li 1, Yun Huang 1, Zhe Cui 1, Stefan Hartmann 2, David J. Benson 3 1 Livermore Software Technology Corporation, Livermore, CA, USA 2
More informationHigh-performance computing of wind turbine aerodynamics using isogeometric analysis
Iowa State University From the SelectedWorks of Ming-Chen Hsu October, 2011 High-performance computing of wind turbine aerodynamics using isogeometric analysis Ming-Chen Hsu, University of California,
More informationMODELLING STIFFENED LIGHTWEIGHT STRUCTURES WITH ISOGEOMETRIC ANALYSIS VIA MORTAR METHODS
ECCOMAS Congress 2016 VII European Congress on Computational Methods in Applied Sciences and Engineering M. Papadrakakis, V. Papadopoulos, G. Stefanou, V. Plevris (eds.) Crete Island, Greece, 5 10 June
More informationLOCAL AND GLOBAL APPROACHES TO FRACTURE MECHANICS USING ISOGEOMETRIC ANALYSIS METHOD
LOCAL AND GLOBAL APPROACHES TO FRACTURE MECHANICS USING ISOGEOMETRIC ANALYSIS METHOD 1* ABDOLGHAFOOR KHADEMALRASOUL, 2 REZA NADERI 1 Ph.D. Student., Department of Civil Engineering, Shahroud University
More informationDynamic Programming Algorithm for Generation of Optimal Elimination Trees for Multi-Frontal Direct Solver over h-refined Grids
Procedia Computer Science Volume 29, 2014, Pages 947 959 ICCS 2014. 14th International Conference on Computational Science Dynamic Programming Algorithm for Generation of Optimal Elimination Trees for
More informationPetIGA. A Framework for High Performance Isogeometric Analysis. Santa Fe, Argentina. Knoxville, United States. Thuwal, Saudi Arabia
PetIGA A Framework for High Performance Isogeometric Analysis Lisandro Dalcin 1,3, Nathaniel Collier 2, Adriano Côrtes 3, Victor M. Calo 3 1 Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationModelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin
Modelling of Isogeometric Analysis for Plane Stress Problem Using Matlab GUI Norliyana Farzana Binti Zulkefli, Ahmad Razin Bin Zainal Abidin Faculty of Civil Engineering, Universiti Teknologi Malaysia,
More informationA POINTWISE APPROACH FOR ENFORCEMENT OF ESSENTIAL BOUNDARY CONDITIONS IN THE ISOGEOMETRIC ANALYSIS * Y. BAZARGAN-LARI **
IJST, Transactions of Mechanical Engineering, Vol. 38, No. M1 +, pp 167-179 Printed in The Islamic Republic of Iran, 2014 Shiraz University A POINTWISE APPROACH FOR ENFORCEMENT OF ESSENTIAL BOUNDARY CONDITIONS
More informationSingularities in Isogeometric Analysis
Paper 45 Singularities in Isogeometric Analysis Civil-Comp Press, 01 Proceedings of the Eighth International Conference on Engineering Computational Technology, B.H.V. Topping, (Editor), Civil-Comp Press,
More informationFEAP - - A Finite Element Analysis Program
FEAP - - A Finite Element Analysis Program Version 8.5 Isogeometric User Manual Robert L. Taylor & Sanjay Govindjee Department of Civil and Environmental Engineering University of California at Berkeley
More informationICES REPORT Volumetric T-spline Construction Using Boolean Operations
ICES REPORT 13-19 July 2013 Volumetric T-spline Construction Using Boolean Operations by L. Liu, Y. Zhang, T.J.R. Hughes, M.A. Scott, T.W. Sederberg The Institute for Computational Engineering and Sciences
More informationMAFELAP Mini-Symposium: Advanced FEM methodologies and isogeometric analysis. Organisers: Matthias Möller and Bernd Simeon
MAFELAP 2016 Conference on the Mathematics of Finite Elements and Applications 14 17 June 2016 Mini-Symposium: Advanced FEM methodologies and isogeometric analysis Organisers: Matthias Möller and Bernd
More informationThe effect of parameterization on isogeometric analysis of free-form curved beams
The effect of parameterization on isogeometric analysis of free-form curved beams Hosseini, SF, Moetakef-Imani, B, Hadidi-Moud, S & Hassani, B Author post-print (accepted) deposited by Coventry University
More informationAvailable online at ScienceDirect. Procedia Engineering 136 (2016 ) Dynamic analysis of fuel tank
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 136 (2016 ) 45 49 The 20 th International Conference: Machine Modeling and Simulations, MMS 2015 Dynamic analysis of fuel tank
More informationA Study of Isogeometric Analysis for Scalar Convection-Diffusion Equations
A Study of Isogeometric Analysis for Scalar Convection-Diffusion Equations Volker John a,b,, Liesel Schumacher c, a Weierstrass Institute for Applied Analysis and Stochastics, Leibniz Institute in Forschungsverbund
More informationA generalized finite element formulation for arbitrary basis functions: from isogeometric analysis to XFEM
A generalized finite element formulation for arbitrary basis functions: from isogeometric analysis to XFEM D.J. Benson 1, Y. Bazilevs 1, E. De Luycker 1, M.-C. Hsu 1, M. Scott 2, T.J.R. Hughes 2, and T.
More informationDevelopment of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics
Development of an Integrated Computational Simulation Method for Fluid Driven Structure Movement and Acoustics I. Pantle Fachgebiet Strömungsmaschinen Karlsruher Institut für Technologie KIT Motivation
More informationCurrent Status of Isogeometric Analysis in LS-DYNA
Current Status of Isogeometric Analysis in LS-DYNA Stefan Hartmann Developer Forum, September 24 th, 2013, Filderstadt, Germany Development in cooperation with: D.J. Benson: Professor of Applied Mechanics,
More informationADVANTAGES OF IDENTICAL SPLINE-BASED INTERFACES FOR FLUID-STRUCTURE INTERACTION ABSTRACT
ADVANTAGES OF IDENTICAL SPLINE-BASED INTERFACES FOR FLUID-STRUCTURE INTERACTION Norbert Hosters, Michel Make, Thomas Spenke, Stefanie Elgeti and Marek Behr Chair for Computational Analysis of Technical
More informationConstruction of polynomial spline spaces over quadtree and octree T-meshes
Construction of polynomial spline spaces over quadtree and octree T-meshes M. Brovka J.I. López J.M. Escobar J.M. Cascón, R. Montenegro October, 2014 Abstract We present a new strategy for constructing
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions,
More informationAvailable online at ScienceDirect. Procedia Engineering 111 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 111 (2015 ) 115 120 XXV R-S-P seminar, Theoretical Foundation of Civil Engineering (24RSP) (TFoCE 2015) Numerical simulation
More informationParallel hp-finite Element Simulations of 3D Resistivity Logging Instruments
Parallel hp-finite Element Simulations of 3D Resistivity Logging Instruments M. Paszyński 1,3, D. Pardo 1,2, L. Demkowicz 1, C. Torres-Verdin 2 1 Institute for Computational Engineering and Sciences 2
More informationAvailable online at ScienceDirect. Procedia Engineering 155 (2016 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 155 (2016 ) 332 341 International Symposium on "Novel structural skins - Improving sustainability and efficiency through new
More informationComputergrafik. Matthias Zwicker. Herbst 2010
Computergrafik Matthias Zwicker Universität Bern Herbst 2010 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling Piecewise Bézier curves Each segment
More informationEXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY
EXACT BUCKLING SOLUTION OF COMPOSITE WEB/FLANGE ASSEMBLY J. Sauvé 1*, M. Dubé 1, F. Dervault 2, G. Corriveau 2 1 Ecole de technologie superieure, Montreal, Canada 2 Airframe stress, Advanced Structures,
More informationRevised Sheet Metal Simulation, J.E. Akin, Rice University
Revised Sheet Metal Simulation, J.E. Akin, Rice University A SolidWorks simulation tutorial is just intended to illustrate where to find various icons that you would need in a real engineering analysis.
More informationApproximation of 3D-Parametric Functions by Bicubic B-spline Functions
International Journal of Mathematical Modelling & Computations Vol. 02, No. 03, 2012, 211-220 Approximation of 3D-Parametric Functions by Bicubic B-spline Functions M. Amirfakhrian a, a Department of Mathematics,
More informationAn algorithm for tensor product approximation of three-dimensional material data for implicit dynamics simulations
An algorithm for tensor product approximation of three-dimensional material data for implicit dynamics simulations Krzysztof Podsiadło, Marcin Łoś, Leszek Siwik, and Maciej Woźniak AGH University of Science
More informationNodal Basis Functions for Serendipity Finite Elements
Nodal Basis Functions for Serendipity Finite Elements Andrew Gillette Department of Mathematics University of Arizona joint work with Michael Floater (University of Oslo) Andrew Gillette - U. Arizona Nodal
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationComputergrafik. Matthias Zwicker Universität Bern Herbst 2016
Computergrafik Matthias Zwicker Universität Bern Herbst 2016 Today Curves NURBS Surfaces Parametric surfaces Bilinear patch Bicubic Bézier patch Advanced surface modeling 2 Piecewise Bézier curves Each
More informationAdaptive Isogeometric Analysis by Local h-refinement with T-splines
Adaptive Isogeometric Analysis by Local h-refinement with T-splines Michael Dörfel 1, Bert Jüttler 2, Bernd Simeon 1 1 TU Munich, Germany 2 JKU Linz, Austria SIMAI, Minisymposium M13 Outline Preliminaries:
More informationICES REPORT February The Institute for Computational Engineering and Sciences The University of Texas at Austin Austin, Texas 78712
ICES REPORT 3-03 February 03 Isogeometric Collocation: Cost Comparison with Galerkin Methods and Extension to Adaptive Hierarchical NURBS Discretizations by Dominik Schillinger, John A. Evans, Alessandro
More informationComputer Graphics Curves and Surfaces. Matthias Teschner
Computer Graphics Curves and Surfaces Matthias Teschner Outline Introduction Polynomial curves Bézier curves Matrix notation Curve subdivision Differential curve properties Piecewise polynomial curves
More informationICES REPORT January D. Schillinger, L. Dede, M.A. Scott, J.A. Evans, M.J. Borden, E. Rank, T.J.R. Hughes
ICES REPORT 12-05 January 2012 An Isogeometric Design-through-analysis Methodology based on Adaptive Hierarchical Refinement of NURBS, Immersed Boundary Methods, and T-spline CAD Surfaces by D. Schillinger,
More informationA substructure based parallel dynamic solution of large systems on homogeneous PC clusters
CHALLENGE JOURNAL OF STRUCTURAL MECHANICS 1 (4) (2015) 156 160 A substructure based parallel dynamic solution of large systems on homogeneous PC clusters Semih Özmen, Tunç Bahçecioğlu, Özgür Kurç * Department
More informationCHAPTER 4. Numerical Models. descriptions of the boundary conditions, element types, validation, and the force
CHAPTER 4 Numerical Models This chapter presents the development of numerical models for sandwich beams/plates subjected to four-point bending and the hydromat test system. Detailed descriptions of the
More information3D NURBS-ENHANCED FINITE ELEMENT METHOD
7th Workshop on Numerical Methods in Applied Science and Engineering (NMASE 8) Vall de Núria, 9 a 11 de enero de 28 c LaCàN, www.lacan-upc.es 3D NURBS-ENHANCED FINITE ELEMENT METHOD R. Sevilla, S. Fernández-Méndez
More informationA NURBS Enhanced extended Finite Element Approach for Unfitted CAD Analysis
A NURBS Enhanced extended Finite Element Approach for Unfitted CAD Analysis Grégory Legrain To cite this version: Grégory Legrain. A NURBS Enhanced extended Finite Element Approach for Unfitted CAD Analysis.
More informationComparison between Bezier and Hermite cubic interpolants in elastic spline formulations
Acta Mech DOI 10.1007/s00707-013-1020-1 Sayed-Mahdi Hashemi-Dehkordi Pier Paolo Valentini Comparison between Bezier and Hermite cubic interpolants in elastic spline formulations Received: 2 August 2013
More informationAbaqus User Element implementation of NURBS based Isogeometric Analysis
1 / 22 Abaqus User Element implementation of NURBS based Isogeometric Analysis 1, A. Duval 1, F. Maurin 1, H. Al Akhras 1 1 Université de Lyon, CNRS INSA-Lyon, Laboratoire de Mécanique des Contacts et
More informationIsogeometric fluid structure interaction analysis with emphasis on non-matching discretizations, and with application to wind turbines
Isogeometric fluid structure interaction analysis with emphasis on non-matching discretizations, and with application to wind turbines Yuri Bazilevs a,, Ming-Chen Hsu a, and Michael A. Scott b a Department
More informationRevision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering. Introduction
Revision of the SolidWorks Variable Pressure Simulation Tutorial J.E. Akin, Rice University, Mechanical Engineering Introduction A SolidWorks simulation tutorial is just intended to illustrate where to
More informationModule 1: Introduction to Finite Element Analysis. Lecture 4: Steps in Finite Element Analysis
25 Module 1: Introduction to Finite Element Analysis Lecture 4: Steps in Finite Element Analysis 1.4.1 Loading Conditions There are multiple loading conditions which may be applied to a system. The load
More informationAvailable online at ScienceDirect. Procedia Engineering 90 (2014 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 90 (2014 ) 288 293 10th International Conference on Mechanical Engineering, ICME 2013 Parallelization of enriched free mesh
More informationSimulation of rotation and scaling algorithm for numerically modelled structures
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Simulation of rotation and scaling algorithm for numerically modelled structures To cite this article: S K Ruhit et al 2018 IOP
More informationICES REPORT January D. Schillinger, L. Dede, M.A. Scott, J.A. Evans, M.J. Borden, E. Rank, T.J.R. Hughes
ICES REPORT 12-05 January 2012 An Isogeometric Design-through-analysis Methodology based on Adaptive Hierarchical Refinement of NURBS, Immersed Boundary Methods, and T-spline CAD Surfaces by D. Schillinger,
More informationAvailable online at ScienceDirect. Procedia Engineering 99 (2015 )
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 99 (2015 ) 575 580 APISAT2014, 2014 Asia-Pacific International Symposium on Aerospace Technology, APISAT2014 A 3D Anisotropic
More informationMODELING OF A MICRO-GRIPPER COMPLIANT JOINT USING COMSOL MULTIPHYSICS SIMULATION
MODELING OF A MICRO-GRIPPER COMPLIANT JOINT USING COMSOL MULTIPHYSICS SIMULATION Mihăiţă Nicolae ARDELEANU, Veronica DESPA, Ioan Alexandru IVAN Valahia University from Targoviste E-mail: mihai.ardeleanu@valahia.ro,
More informationRELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS
International Conference on Economic Engineering and Manufacturing Systems Braşov, 26 27 November 2009 RELIABILITY OF THE FEM CALCULATIONS OF THE FRACTURE MECHANICS PARAMETERS Galina TODOROVA, Valentin
More informationAvailable online at ScienceDirect. Energy Procedia 69 (2015 )
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 69 (2015 ) 1885 1894 International Conference on Concentrating Solar Power and Chemical Energy Systems, SolarPACES 2014 Heliostat
More informationParameterization adaption for 3D shape optimization in aerodynamics
Parameterization adaption for 3D shape optimization in aerodynamics Badr Abou El Majd Computer Science and Decision Aiding Laboratory Hassan II University B.P. 5366 Maârif, Casablanca, Morocco b.abouelmajd@fsac.ac.ma
More informationIsogeometric Finite Element Analysis Using T- Splines
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2009-08-12 Isogeometric Finite Element Analysis Using T- Splines Jingang Li Brigham Young University - Provo Follow this and additional
More informationTo mesh or not to mesh: that is the question
To mesh or not to mesh: that is the question Stefan Kollmannsberger Computation in Engineering Computation in Engineering Ordinarius: Prof. Dr.rer.nat. Ernst Rank Numerical Mechanics: Dr.-Ing. Stefan Kollmannsberger,
More informationAN APPROACH FOR LOAD BALANCING FOR SIMULATION IN HETEROGENEOUS DISTRIBUTED SYSTEMS USING SIMULATION DATA MINING
AN APPROACH FOR LOAD BALANCING FOR SIMULATION IN HETEROGENEOUS DISTRIBUTED SYSTEMS USING SIMULATION DATA MINING Irina Bernst, Patrick Bouillon, Jörg Frochte *, Christof Kaufmann Dept. of Electrical Engineering
More informationMethod of Finite Elements I
Institute of Structural Engineering Page 1 Held by Prof. Dr. E. Chatzi, Dr. P. Steffen Assistants: Adrian Egger (HIL E 13.3), Harry Mylonas (HIL H33.1), Konstantinos Tatsis (HIL H33.1) Lectures homepage:
More informationVOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL
VOLCANIC DEFORMATION MODELLING: NUMERICAL BENCHMARKING WITH COMSOL The following is a description of the model setups and input/output parameters for benchmarking analytical volcanic deformation models
More informationNURBS: Non-Uniform Rational B-Splines AUI Course Denbigh Starkey
NURBS: Non-Uniform Rational B-Splines AUI Course Denbigh Starkey 1. Background 2 2. Definitions 3 3. Using NURBS to define a circle 4 4. Homogeneous coordinates & control points at infinity 9 5. Constructing
More informationA numerical study of Isogeometric Analysis for cavity flow and flow around a cylindrical obstacle
A numerical study of Isogeometric Analysis for cavity flow and flow around a cylindrical obstacle Babak S. Hosseini 1, Matthias Möller 2, and Stefan Turek 1 1 TU Dortmund, Institute of Applied Mathematics
More informationLocal Modification of Subdivision Surfaces Based on Curved Mesh
Local Modification of Subdivision Surfaces Based on Curved Mesh Yoshimasa Tokuyama Tokyo Polytechnic University tokuyama@image.t-kougei.ac.jp Kouichi Konno Iwate University konno@cis.iwate-u.ac.jp Junji
More informationIsogeometric Analysis (IGA) Part I (MS62) and Part II (MS78)
Isogeometric Analysis (IGA) Part I (MS62) and Part II (MS78) July 11, 4:00-6:00 PM Challenges in Isogeometric Analysis (IGA) Tor Dokken, Oslo, Norway Locally Refined B-splines, Tom Lyche, University of
More informationCOORDINATE MEASUREMENTS OF COMPLEX-SHAPE SURFACES
XIX IMEKO World Congress Fundamental and Applied Metrology September 6 11, 2009, Lisbon, Portugal COORDINATE MEASUREMENTS OF COMPLEX-SHAPE SURFACES Andrzej Werner 1, Malgorzata Poniatowska 2 1 Faculty
More informationMotivation Patch tests Numerical examples Conclusions
Weakening the tight coupling between geometry and simulation in isogeometric analysis: from sub- and super- geometric analysis to Geometry Independent Field approximation (GIFT) Elena Atroshchenko, Gang
More informationConstrained modification of the cubic trigonometric Bézier curve with two shape parameters
Annales Mathematicae et Informaticae 43 (014) pp. 145 156 http://ami.ektf.hu Constrained modification of the cubic trigonometric Bézier curve with two shape parameters Ede Troll University of Debrecen
More informationNon-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, Politecnico di Milano, February 3, 2017, Lesson 1
Non-Linear Finite Element Methods in Solid Mechanics Attilio Frangi, attilio.frangi@polimi.it Politecnico di Milano, February 3, 2017, Lesson 1 1 Politecnico di Milano, February 3, 2017, Lesson 1 2 Outline
More informationIsogeometric Simulation of Turbine Blades for Aircraft Engines
Isogeometric Simulation of Turbine Blades for Aircraft Engines David Großmann 1, Bert Jüttler 2, Helena Schlusnus 3, Johannes Barner 1, Anh-Vu Vuong 4 1 MTU Aero Engines GmbH, Munich, Germany 2 Johannes
More informationSolving Partial Differential Equations on Overlapping Grids
**FULL TITLE** ASP Conference Series, Vol. **VOLUME**, **YEAR OF PUBLICATION** **NAMES OF EDITORS** Solving Partial Differential Equations on Overlapping Grids William D. Henshaw Centre for Applied Scientific
More informationGenerating Airplane Wings for Numerical Simulation and Manufacturing
Generating Airplane Wings for Numerical Simulation and Manufacturing Karl-Heinz Brakhage Philipp Lamby Institute of Geometry and Applied Mathematics RWTH Aachen, University of Technology D-52056 Aachen,
More informationAPPLICATIONS OF CFD AND PCA METHODS FOR GEOMETRY RECONSTRUCTION OF 3D OBJECTS
Mathematical Modelling and Analysis 2005. Pages 123 128 Proceedings of the 10 th International Conference MMA2005&CMAM2, Trakai c 2005 Technika ISBN 9986-05-924-0 APPLICATIONS OF CFD AND PCA METHODS FOR
More informationSolving a Two Dimensional Unsteady-State. Flow Problem by Meshless Method
Applied Mathematical Sciences, Vol. 7, 203, no. 49, 242-2428 HIKARI Ltd, www.m-hikari.com Solving a Two Dimensional Unsteady-State Flow Problem by Meshless Method A. Koomsubsiri * and D. Sukawat Department
More informationMAHARAVO H. RANDRIANARIVONY
MAHARAVO H. RANDRIANARIVONY EDUCATION: Technical University of Chemnitz, Germany Ph. D. in Computer Science. 2006 Supervisor: Prof. Guido Brunnett. Technical University of Chemnitz, Germany M. Sc. in Mathematics
More informationIntroduction to Finite Element Analysis using ANSYS
Introduction to Finite Element Analysis using ANSYS Sasi Kumar Tippabhotla PhD Candidate Xtreme Photovoltaics (XPV) Lab EPD, SUTD Disclaimer: The material and simulations (using Ansys student version)
More informationImprovement of Reduced Order Modeling based on POD
Author manuscript, published in "The Fifth International Conference on Computational Fluid Dynamics (28)" Improvement of Reduced Order Modeling based on POD M. Bergmann 1, C.-H. Bruneau 2, and A. Iollo
More informationFree-Form Shape Optimization using CAD Models
Free-Form Shape Optimization using CAD Models D. Baumgärtner 1, M. Breitenberger 1, K.-U. Bletzinger 1 1 Lehrstuhl für Statik, Technische Universität München (TUM), Arcisstraße 21, D-80333 München 1 Motivation
More informationIntroduction to the Finite Element Method (3)
Introduction to the Finite Element Method (3) Petr Kabele Czech Technical University in Prague Faculty of Civil Engineering Czech Republic petr.kabele@fsv.cvut.cz people.fsv.cvut.cz/~pkabele 1 Outline
More informationMODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
A. N. Johnson et al., Int. J. Comp. Meth. and Exp. Meas., Vol. 3, No. 3 (2015) 269 278 MODELING MIXED BOUNDARY PROBLEMS WITH THE COMPLEX VARIABLE BOUNDARY ELEMENT METHOD (CVBEM) USING MATLAB AND MATHEMATICA
More informationAlgebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations
Algebraic Flux Correction Schemes for High-Order B-Spline Based Finite Element Approximations Numerical Analysis group Matthias Möller m.moller@tudelft.nl 7th DIAM onderwijs- en onderzoeksdag, 5th November
More informationINTERSECTION OF CURVES FACTS, COMPUTATIONS, APPLICATIONS IN BLOWUP*
South Bohemia Mathematical Letters Volume 24, (2016), No. 1, 10-16. INTERSECTION OF CURVES FACTS, COMPUTATIONS, APPLICATIONS IN BLOWUP* PAVEL CHALMOVIANSKÝ abstrakt. We deal with application of intersection
More informationSecond International Workshop on Scientific Computing and Applications. Kananaskis, Canada, May 28 - June 1, 2000
Second International Workshop on Scientific Computing and Applications. Kananaskis, Canada, May 28 - June 1, 2000 Program May 28 (Sunday) 19:00-21:00 Registration and reception Session Chairman: Y. Wong
More informationParallel Mesh Multiplication for Code_Saturne
Parallel Mesh Multiplication for Code_Saturne Pavla Kabelikova, Ales Ronovsky, Vit Vondrak a Dept. of Applied Mathematics, VSB-Technical University of Ostrava, Tr. 17. listopadu 15, 708 00 Ostrava, Czech
More informationMESH INDEPENDENCE AND CONVERGENCE TEST IN THE OIL FILM CHARACTERISTICS OF A TILTING PAD THRUST BEARING
Mesh independence and convergence test in the oil film characteristics MESH INDEPENDENCE AND CONVERGENCE TEST IN THE OIL FILM CHARACTERISTICS OF A TILTING PAD THRUST BEARING D.V.Srikanth a, K.K.Chaturvedi
More information3D Resistivity Logging Instruments
Parallel hp-finite Element Simulations of 3D Resistivity Logging Instruments M.Paszyński (1), D.Pardo (2), L.Demkowicz (3), C.Torres-Verdin (2) (1) Department of Modeling and Information Technology, AGH
More informationFinite Element Method. Chapter 7. Practical considerations in FEM modeling
Finite Element Method Chapter 7 Practical considerations in FEM modeling Finite Element Modeling General Consideration The following are some of the difficult tasks (or decisions) that face the engineer
More information