Hdglab: an Open-Source Implementation of the Hybridisable Discontinuous Galerkin Method in MATLAB
Total Page:16
File Type:pdf, Size:1020Kb

Load more
Recommended publications
-
Arxiv:1911.09220V2 [Cs.MS] 13 Jul 2020
MFEM: A MODULAR FINITE ELEMENT METHODS LIBRARY ROBERT ANDERSON, JULIAN ANDREJ, ANDREW BARKER, JAMIE BRAMWELL, JEAN- SYLVAIN CAMIER, JAKUB CERVENY, VESELIN DOBREV, YOHANN DUDOUIT, AARON FISHER, TZANIO KOLEV, WILL PAZNER, MARK STOWELL, VLADIMIR TOMOV Lawrence Livermore National Laboratory, Livermore, USA IDO AKKERMAN Delft University of Technology, Netherlands JOHANN DAHM IBM Research { Almaden, Almaden, USA DAVID MEDINA Occalytics, LLC, Houston, USA STEFANO ZAMPINI King Abdullah University of Science and Technology, Thuwal, Saudi Arabia Abstract. MFEM is an open-source, lightweight, flexible and scalable C++ library for modular finite element methods that features arbitrary high-order finite element meshes and spaces, support for a wide variety of dis- cretization approaches and emphasis on usability, portability, and high-performance computing efficiency. MFEM's goal is to provide application scientists with access to cutting-edge algorithms for high-order finite element mesh- ing, discretizations and linear solvers, while enabling researchers to quickly and easily develop and test new algorithms in very general, fully unstructured, high-order, parallel and GPU-accelerated settings. In this paper we describe the underlying algorithms and finite element abstractions provided by MFEM, discuss the software implementation, and illustrate various applications of the library. arXiv:1911.09220v2 [cs.MS] 13 Jul 2020 1. Introduction The Finite Element Method (FEM) is a powerful discretization technique that uses general unstructured grids to approximate the solutions of many partial differential equations (PDEs). It has been exhaustively studied, both theoretically and in practice, in the past several decades [1, 2, 3, 4, 5, 6, 7, 8]. MFEM is an open-source, lightweight, modular and scalable software library for finite elements, featuring arbitrary high-order finite element meshes and spaces, support for a wide variety of discretization approaches and emphasis on usability, portability, and high-performance computing (HPC) efficiency [9]. -
Codesaturne Practical User's Guide
EDF R&D Fluid Dynamics, Power Generation and Environment Department Single Phase Thermal-Hydraulics Group 6, quai Watier F-78401 Chatou Cedex Tel: 33 1 30 87 75 40 Fax: 33 1 30 87 79 16 JUNE 2017 Code Saturne documentation Code Saturne version 5.0.0 practical user's guide contact: [email protected] http://code-saturne.org/ c EDF 2017 Code Saturne EDF R&D Code Saturne version 5.0.0 practical user's documentation guide Page 1/142 ABSTRACT Code Saturne is a system designed to solve the Navier-Stokes equations in the cases of 2D, 2D ax- isymmetric or 3D flows. Its main module is designed for the simulation of flows which may be steady or unsteady, laminar or turbulent, incompressible or potentially dilatable, isothermal or not. Scalars and turbulent fluctuations of scalars can be taken into account. The code includes specific modules, referred to as \specific physics", for the treatment of Lagrangian particle tracking, semi-transparent radiative transfer, gas combustion, pulverised coal combustion, electricity effects (Joule effect and elec- tric arcs) and compressible flows. Code Saturne relies on a finite volume discretisation and allows the use of various mesh types which may be hybrid (containing several kinds of elements) and may have structural non-conformities (hanging nodes). The present document is a practical user's guide for Code Saturne version 5.0.0. It is the result of the joint effort of all the members in the development team. It presents all the necessary elements to run a calculation with Code Saturne version 5.0.0. -
MFEM: a Modular Finite Element Methods Library
MFEM: A Modular Finite Element Methods Library Robert Anderson1, Andrew Barker1, Jamie Bramwell1, Jakub Cerveny2, Johann Dahm3, Veselin Dobrev1,YohannDudouit1, Aaron Fisher1,TzanioKolev1,MarkStowell1,and Vladimir Tomov1 1Lawrence Livermore National Laboratory 2University of West Bohemia 3IBM Research July 2, 2018 Abstract MFEM is a free, lightweight, flexible and scalable C++ library for modular finite element methods that features arbitrary high-order finite element meshes and spaces, support for a wide variety of discretization approaches and emphasis on usability, portability, and high-performance computing efficiency. Its mission is to provide application scientists with access to cutting-edge algorithms for high-order finite element meshing, discretizations and linear solvers. MFEM also enables researchers to quickly and easily develop and test new algorithms in very general, fully unstructured, high-order, parallel settings. In this paper we describe the underlying algorithms and finite element abstractions provided by MFEM, discuss the software implementation, and illustrate various applications of the library. Contents 1 Introduction 3 2 Overview of the Finite Element Method 4 3Meshes 9 3.1 Conforming Meshes . 10 3.2 Non-Conforming Meshes . 11 3.3 NURBS Meshes . 12 3.4 Parallel Meshes . 12 3.5 Supported Input and Output Formats . 13 1 4 Finite Element Spaces 13 4.1 FiniteElements....................................... 14 4.2 DiscretedeRhamComplex ................................ 16 4.3 High-OrderSpaces ..................................... 17 4.4 Visualization . 18 5 Finite Element Operators 18 5.1 DiscretizationMethods................................... 18 5.2 FiniteElementLinearSystems . 19 5.3 Operator Decomposition . 23 5.4 High-Order Partial Assembly . 25 6 High-Performance Computing 27 6.1 Parallel Meshes, Spaces, and Operators . 27 6.2 Scalable Linear Solvers . -
FELICITY: a MATLAB/C++ TOOLBOX for DEVELOPING FINITE ELEMENT METHODS and SIMULATION MODELING\Ast
SIAM J. SCI.COMPUT. \bigcircc 2018 Society for Industrial and Applied Mathematics Vol. 40, No. 2, pp. C234{C257 FELICITY: A MATLAB/C++ TOOLBOX FOR DEVELOPING FINITE ELEMENT METHODS AND SIMULATION MODELING\ast SHAWN W. WALKERy Abstract. This paper describes a MATLAB/C++ finite element toolbox, called FELICITY, for simulating various types of systems of partial differential equations (e.g., coupled elliptic/parabolic problems) using the finite element method. It uses MATLAB in an object-oriented way for high-level manipulation of data structures in finite element codes, while utilizing a domain-specific language (DSL) and code generation to automate low-level tasks such as matrix assembly (via the MATLAB mex interface). We describe the fundamental functionality of the toolbox's MATLAB interface, such as using higher order Lagrange (simplicial) meshes, defining finite element spaces, allocating degrees- of-freedom, assembling discrete bilinear and linear forms, and interpolation over meshes. Moreover, we describe in-depth how automatic code generation is implemented in FELICITY. Two example problems and their implementation are provided to demonstrate the ability of FELICITY to solve coupled problems with interacting subdomains of different co-dimension. Future improvements are also discussed. Key words. finite elements, coupled systems, geometric flows, code generation, MATLAB, open source software AMS subject \bfc \bfl \bfa \bfs \bfifi\bfc\bfa\bft\bfo\bfn. 68N30, 65N30, 65M60, 68N19, 68N20 DOI. 10.1137/17M1128745 1. Introduction. The development of numerical methods to solve partial dif- ferential equations (PDEs) continues to advance to address new domain areas and problems of increasing complexity. With this, the number of software packages has increased to address the needs for researching new methods and simulating large prob- lems. -
Mofem: an Open Source, Parallel Finite Element Library
MoFEM: An open source, parallel finite element library Kaczmarczyk, Ł., Ullah, Z., Lewandowski, K., Meng, X., Zhou, X-Y., Athanasiadis, I., Nguyen, H., Chalons- Mouriesse, C-A., Richardson, E. J., Miur, E., Shvarts, A. G., Wakeni, M., & Pearce, C. J. (2020). MoFEM: An open source, parallel finite element library. The Journal of Open Source Software, 5(45). https://doi.org/doi.org/10.21105/joss.01441 Published in: The Journal of Open Source Software Document Version: Publisher's PDF, also known as Version of record Queen's University Belfast - Research Portal: Link to publication record in Queen's University Belfast Research Portal Publisher rights © 2020 The Authors. This is an open access article published under a Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution and reproduction in any medium, provided the author and source are cited. General rights Copyright for the publications made accessible via the Queen's University Belfast Research Portal is retained by the author(s) and / or other copyright owners and it is a condition of accessing these publications that users recognise and abide by the legal requirements associated with these rights. Take down policy The Research Portal is Queen's institutional repository that provides access to Queen's research output. Every effort has been made to ensure that content in the Research Portal does not infringe any person's rights, or applicable UK laws. If you discover content in the Research Portal that you believe breaches copyright or violates any law, please contact [email protected]. Download date:06. -
Developments in the Modeling & Simulation Program at EDF
Developments in the Modeling & Simulation Program at EDF. Potential Collaboration Topics CASL Industry Council Meeting. Charleston, SC. April 4-5 2017 Didier Banner Presentation outline EDF’s M&S tools and software policy Current trends in numerical simulation -------------------- On-going CASL – CEA –EDF collaboration Potential collaboration on the NESTOR data EDF Key figures • French NPP fleet • 58 operating reactors, from 900 MW to 1450 MW • 157 to 205 fuel assemblies per reactor • Fuel cycles - 12 or 18 months • Fuel assemblies renewal from 1/4th to 1/3rd • Some estimated costs* • One day of outage: ~1 M€ • Total fuel cost: ~5 €/MWh • Major retrofit in France: ~50 b€ Including post-Fukushima program: ~10 b€ EDF R&D KEY FIGURES Use of Modelling &Simulation - examples Resistance to impact Tightness of the (projectiles) Seismic Analysis containment vessel Environmental impacts Behaviour of turbines Dismantling Waste Storage Tightness of the primary loop Control of nuclear Behaviour of the reactions pressure vessel EDF Modeling and Simulation policy Models Specific studies: i.e FSI interaction,irradiation, turbulence,.. Codes i.e CFD (Saturne), Neutronics (Cocagne),Mechanics (Aster) Platforms Interoperability, Users’s experience --------------------- Development Strategy - examples EDF Open-Source CFD (Saturne), Mechanics(Aster), Free Surface Flow EDF developments-not open source Neutronics, Electromagnetics, … Codevelopment/Partnership Two-phase flow (Neptune), Fast transient dynamics,.. Commercial Software: Ansys, Abaqus, EDF Modeling -
Development of a Coupling Approach for Multi-Physics Analyses of Fusion Reactors
Development of a coupling approach for multi-physics analyses of fusion reactors Zur Erlangung des akademischen Grades eines Doktors der Ingenieurwissenschaften (Dr.-Ing.) bei der Fakultat¨ fur¨ Maschinenbau des Karlsruher Instituts fur¨ Technologie (KIT) genehmigte DISSERTATION von Yuefeng Qiu Datum der mundlichen¨ Prufung:¨ 12. 05. 2016 Referent: Prof. Dr. Stieglitz Korreferent: Prof. Dr. Moslang¨ This document is licensed under the Creative Commons Attribution – Share Alike 3.0 DE License (CC BY-SA 3.0 DE): http://creativecommons.org/licenses/by-sa/3.0/de/ Abstract Fusion reactors are complex systems which are built of many complex components and sub-systems with irregular geometries. Their design involves many interdependent multi- physics problems which require coupled neutronic, thermal hydraulic (TH) and structural mechanical (SM) analyses. In this work, an integrated system has been developed to achieve coupled multi-physics analyses of complex fusion reactor systems. An advanced Monte Carlo (MC) modeling approach has been first developed for converting complex models to MC models with hybrid constructive solid and unstructured mesh geometries. A Tessellation-Tetrahedralization approach has been proposed for generating accurate and efficient unstructured meshes for describing MC models. For coupled multi-physics analyses, a high-fidelity coupling approach has been developed for the physical conservative data mapping from MC meshes to TH and SM meshes. Interfaces have been implemented for the MC codes MCNP5/6, TRIPOLI-4 and Geant4, the CFD codes CFX and Fluent, and the FE analysis platform ANSYS Workbench. Furthermore, these approaches have been implemented and integrated into the SALOME simulation platform. Therefore, a coupling system has been developed, which covers the entire analysis cycle of CAD design, neutronic, TH and SM analyses. -
Overview of Research in the Scidac Center for Simulation of Fusion Relevant RF Actuators P
Overview of Research in the SciDAC Center for Simulation of Fusion Relevant RF Actuators P. T. Bonoli and D. L. Green on behalf of the RF-SciDAC Team 60th Annual Meeting of the APS Division of Plasma Physics November 5–9, 2018 Portland, Oregon Poster BP11.00071 RF-SciDAC : Center for Simulation of Fusion Relevant RF Actuators http://rfscidac4.org Abstract An overview is given of research which has the goal of developing a predictive simulation capability of the self-consistent interaction of RF power with the scrape-off layer (SOL), including the effects of plasma sheaths, ponderomotive forces near an antenna, and turbulence and transport. This will make it possible to answer critical questions related to how RF power modifies properties of the SOL, and how, in turn, the SOL affects the propagation and absorption of RF waves. Targeted problems include the impact of high power RF systems on plasma facing materials, including high-Z impurity sputtering and transport induced by large RF-induced sheath potentials, localized thermal loads, and antenna damage. In order to address the high geometric fidelity required to describe the 3D magnetic field and 3D solid geometry of the RF launching structures and the surrounding vacuum vessel the open-source scalable Modular Finite Element Framework (MFEM) is being utilized. The Parallel Unstructured Mesh Infrastructure (PUMI) is also being used in order to move to global solution domains, which are presently precluded by existing meshing approaches. *Work supported by US DoE contract numbers DE-SC0018090, DE-SC0018319, DE- SC0018275, FWP 3ERAT952, FWP 2017-LLNL-SCW1619, and Work Proposal 3203. -
Modeling and Distributed Computing of Snow Transport And
Modeling and distributed computing of snow transport and delivery on meso-scale in a complex orography Modelitzaci´oi computaci´odistribu¨ıda de fen`omens de transport i dip`osit de neu a meso-escala en una orografia complexa Thesis dissertation submitted in fulfillment of the requirements for the degree of Doctor of Philosophy Programa de Doctorat en Societat de la Informaci´oi el Coneixement Alan Ward Koeck Advisor: Dr. Josep Jorba Esteve Distributed, Parallel and Collaborative Systems Research Group (DPCS) Universitat Oberta de Catalunya (UOC) Rambla del Poblenou 156, 08018 Barcelona – Barcelona, 2015 – c Alan Ward Koeck, 2015 Unless otherwise stated, all artwork including digital images, sketches and line drawings are original creations by the author. Permission is granted to copy, distribute and/or modify this document under the terms of the Creative Commons BY-SA License, version 4.0 or ulterior at the choice of the reader/user. At the time of writing, the license code was available at: https://creativecommons.org/licenses/by-sa/4.0/legalcode Es permet la lliure c`opia, distribuci´oi/o modificaci´od’aquest document segons els termes de la Lic`encia Creative Commons BY-SA, versi´o4.0 o posterior, a l’escollida del lector o usuari. En el moment de la redacci´o d’aquest text, es podia accedir al text de la llic`encia a l’adre¸ca: https://creativecommons.org/licenses/by-sa/4.0/legalcode 2 In memoriam Alan Ward, MA Oxon, PhD Dublin 1937-2014 3 4 Acknowledgements A long-term commitment such as this thesis could not prosper on my own merits alone. -
Improved Support for Parallel Adaptive Simulation in CEED WBS 2.2.6.06, Milestone CEED-MS29
ECP Milestone Report Improved Support for Parallel Adaptive Simulation in CEED WBS 2.2.6.06, Milestone CEED-MS29 Mark Shephard Valeria Barra Jed Brown Jean-Sylvain Camier Veselin Dobrev Yohan Dudouit Paul Fischer Tzanio Kolev David Medina Misun Min Cameron Smith Morteza H. Siboni Jeremy Thompson Tim Warburton July 12, 2019 DOCUMENT AVAILABILITY Reports produced after January 1, 1996, are generally available free via US Department of Energy (DOE) SciTech Connect. Website http://www.osti.gov/scitech/ Reports produced before January 1, 1996, may be purchased by members of the public from the following source: National Technical Information Service 5285 Port Royal Road Springfield, VA 22161 Telephone 703-605-6000 (1-800-553-6847) TDD 703-487-4639 Fax 703-605-6900 E-mail [email protected] Website http://www.ntis.gov/help/ordermethods.aspx Reports are available to DOE employees, DOE contractors, Energy Technology Data Exchange representatives, and International Nuclear Information System representatives from the following source: Office of Scientific and Technical Information PO Box 62 Oak Ridge, TN 37831 Telephone 865-576-8401 Fax 865-576-5728 E-mail [email protected] Website http://www.osti.gov/contact.html This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. -
Femengine: Finite Element Method C++ Code Based on Functional And
2019 Ivannikov Ispras Open Conference (ISPRAS) FEMEngine: finite element method C++ code based on functional and template metaprogramming Alexey Gurin∗, Alexey Baykin∗, Timofey Polyansky‡, Anton Krivtsov§¶ ∗Laboratory of Digital and Intelligent Hydrocarbon Production Systems Lavrentyev’s Institute of Hydrodynamics of the SB RAS, Novosibirsk, Russia Email: [email protected] ‡Department of Physics Novosibirsk State University, Novosibirsk, Russia §Higher School of Theoretical Mechanics Peter the Great St. Petersburg Polytechnic University (SPbPU), St. Petersburg, Russia ¶Laboratory ”Discrete models in mechanics” Institute for Problems in Mechanical Engineering, St. Petersburg, Russia acheived by the metaprogramming technique. The FEniCS Abstract—The paper discusses the problems of a finite element method programming. Modern C++ functional programming and Form Compiler, written in Python as a part of the FEniCS template metaprogramming approach to finite element analysis Project, generates and compile an efficient C++ code for a is presented. This approach simplifies the implementation of matrix assembly on the fly. This form of metaprogramming is an effective assembly of the stiffness matrix for a problem universal, but very difficult to implement and maintain. defined by a weak form. This method is tested by a solution During the last decade, the C++ language has been of the Poisson equation on an unstructured 3D tetrahedral mesh using FEM C++ library FEMEngine developed by authors. The greatly improved. Three new standards are introduced: C++11, function which calculates the matrix is generated by higher order C++14, C++17. The ability to program in a functional functions during the compilation stage. The performance of the paradigm has been added and the template metaprogramming computation is analyzed by studying of a disassembled code and technique has been simplified. -
(OSCAE.Initiative) at Universiti Teknologi Malaysia
Open Source Computer Aided Engineering Inititive (OSCAE.Initiative) at Universiti Teknologi Malaysia Abu Hasan Abdullah Faculty of Mechanical Engineering Universiti Teknologi Malaysia March 1, 2017 Abstract Open Source Computer Aided Engineering Initiative (OSCAE.Initiative) is a public statement of principles relating to open source software for Computer Aided Engineering. It aims to pro- mote the use of open source software in engineering discipline with a goal that within the next five years, open source software will become very common tools for conducting CAE analyses. Towards this end, a small OSCAE.Initiative lab to be used for familiarizing and developmental testing of related software packages has been setup at the Marine Technology Centre, Universiti Teknology Malaysia. This paper describes the infrastructure setup, computing hardware components, multiplat- form operating environments and various categories of software packages offered at this lab which has already been identified as a blueprint for another two labs sharing a kindred spirit. 1 Overview The OSCAE.Initiative lab at Universiti Teknologi Malaysia, is located on the second floor of the Marine Technology Centre (MTC). It was established to promote the use of open source software in engineering discipline with a goal that open source software will become the tools of choice for conducting CAE analyses and simulations. Currently, the lab provides • general purpose computing facilities for educational and academic use, and • specialized (CAE) lab environments, software, and support for students, researchers and staff of MTC. Access is by terminals throughout the lab. While users are not required to have their own computers, it is recognized that many do, and facilities are provided to transfer data to and from the OSCAE.Initiative systems so that personal computers and workstations at OSCAE.Initiative lab can complement each other.