A complete methodology for the implementation of XFEM models. Team: - Carlos Hernán Villanueva - Kai Yu
|
|
- Hester Fitzgerald
- 5 years ago
- Views:
Transcription
1 A complete methodology for the implementation of XFEM models Team: - Carlos Hernán Villanueva - Kai Yu
2 Problem Statement QUAD4 DISCONTINUITY
3 Problem Statement Traditional FEM Mesh generates discrete representation of potentially complex geometry Changes in topology or shape Require new mesh Discontinuities must align with element edges
4 Problem Statement GOAL OBJECTIVE Develop a numerical framework Solving problems of complex geometry in 2D and 3D Using XFEM and LSM Decompose a 3D element Arbitrary complex geometry Homogeneous domains capable of integration
5 Background: XFEM Numerical technique Provides local enrichment For solutions to differential equations with discontinuous functions Avoids need to re-mesh Solve shortcomings of FEM
6 Background: Discontinuity A discontinuity can be defined as a rapid change in a field variable (a) Cracks, (b) Shear bands, (c) Interface of two fluids Image from: oduction/figures/fig1.png
7 Background: Level-set function Numerical Scheme Discontinuity of interest is represented as the zero level-set function XFEM XFEM places discontinuity at the boundary layer dividing the grid into negative and positive phases Zero level-set Positive phase Negative phase Image from: res/fig3.png
8 Background: XFEM and LSM LSM XFEM LSM: - Model discontinuity - Update its motion at each calculation XFEM: - Solves the problem - Determines direction of discontinuity
9 Background: Delaunay triangulation Sub-division of geometric objects into integrable domains such as triangles and tetrahedra Circumcircle of any triangle does not contain the vertices of other triangle Image from : ommons/c/c9/delaunay_circumcircles. png
10 Methodology: Setup 4-element 2D mesh. Red areas: phase 2, positive level-set value; blue areas: phase 1 negative level-set value at the node
11 Methodology: Setup 4-element 2D mesh. Red numbers represent the Global Element Id, blue numbers represent the Global Node Id.
12 Methodology: Overview
13 Methodology: Point-to-cell list Point Id Number of cells Cell Ids connected , ,2,3, , , ,
14 Methodology: Overview 4-element 2D mesh. Red numbers represent the Global Element Id, blue numbers represent the Global Node Id.
15 Methodology: Point-to-edge list Global Edge Id Point Ids Cell Ids Local Edge Num
16 Methodology: Edges Edge representation in a QUAD4 element. Green edges represent the local edge index at the element.
17 Methodology: Intersection points Mapping of intersection point to element.
18 Delaunay triangulation Phase 2 Phase 1 Combined
19 Enrichment algorithm Sub-phase (element level) Enrichment (nodal level) Assign degrees of freedom (nodal level)
20 Sub-phase algorithm
21 Enrichment Levels Phase 1 Phase 2 TEMP TEMP_e8 TEMP_e1 TEMP_e9 TEMP_e2 TEMP_e10 TEMP_e3 TEMP_e11 TEMP_e4 TEMP_e12 TEMP_e5 TEMP_e13 TEMP_e6 TEMP_e14
22 Sub-phase algorithm With 27 enrichment levels Triangle number Main Phase Sub-Phase
23 Methodology: Overview 4-element 2D mesh. Red numbers represent the Global Element Id, blue numbers represent the Global Node Id.
24 Enrichment algorithm Nodes/Ele ments C C C C C
25 Enrichment algorithm Nodes/Elem ents 1 U 15 2 C C C U 16 6 C U 17 8 C U 18
26 Methodology: Setup 4-element 2D mesh. Red areas: phase 2, positive level-set value; blue areas: phase 1 negative level-set value at the node
27 Results: Gauss points Gauss points on volume: Phase 1, tetrahedron 1
28 Results: Gauss points on interface and Normal vector Gauss points on interface and normal vector to surface: Phase 1, tetrahedron 1
29 Corroborate results Compare different interface contribution formulations (Stabilized Lagrangian, Nietsche, etc.) and declare a winner By using a combination of different mesh refinements, constraint parameters and property ratios Check how variation of these methods and values affect the results of XFEM Develop a protocol e i e (u + u ) 2 dγ i i dγ i (u 1 u 2 ) 2 dγ (u 2 ) 2 dγ
30 Thermal Tests Mesh size: 20x20,30x30,40x40,50x50 Conductivity ratios: 0.01,0.1,1,10,100,1000 Stabilization factor: 1,10,100,1000 Nominal values are in bold face Vary values of one parameter (mesh size, conductivity ratio, etc) while setting the other parameters to the nominal values Sweep radius of circular inclusion between 2 and 6 using 500 increments
31 Results Mesh refinement sweep Mesh size: variable Conductivity ratios: k2/k1=10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: UMFPACK Pre-conditioner: No scaling Mean Interface and L-2 Errors over Mesh-Refinement Levels Refinement Mean Error Maximum Minimum Std. Dev
32 Results Stabilization factor sweep Mesh size: 30x30 Conductivity ratios: k2/k1=10 Stabilization factor: variable Interface formulation: stabilized Lagrange Solver: UMFPACK Pre-conditioner: Spatial Derivatives Shape Functions Maximum Mean Interface and L-2 Errors of Temperature over Stabilization Factor Levels Refinement Mean Error Maximum Minimum Std. Dev
33 Results Conductivity ratio sweep Mesh size: 30x30 Conductivity ratios: k2/k1=variable Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: UMFPACK Pre-conditioner: Spatial Derivatives Shape Functions Maximum Mean Interface and L-2 Errors of Temperature over Conductivity Ratio Levels Refinement Mean Error E E E E E E E E-04 Maximum E E E E E E E E-04 Minimum E E E E E E E E-04 Std. Dev E E E E E E E E-05
34 Results UMFPACK/GMRES sweep Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: variable Pre-conditioner: Spatial Derivatives Shape Functions Maximum NOTE: UMFPACK and GMRES yielded same condition number results. GMRES pre-conditioner not dumped to Matlab Time UMPFACK, no T scaling UMFPACK, T scaling GMRES, ILU, no T scaling GMRES, ILU, T scaling Condition # Condition # Enrichment Condition # Condition # Enrichment Mean Error E E E E E E+05 Maximum E E E E E E+07 Minimum E E E E E E+00 Std. Dev E E E E E E+06
35 Results UMFPACK/GMRES sweep Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: variable Pre-conditioner: Spatial Derivatives Shape Functions Maximum INTERFACE ERROR
36 Results UMFPACK/GMRES sweep Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: variable Pre-conditioner: Spatial Derivatives Shape Functions Maximum L2 ERROR
37 Results UMFPACK/GMRES sweep Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: variable Pre-conditioner: Spatial Derivatives Shape Functions Maximum CONDITION NUMBER
38 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Shape Functions Maximum Differences in L2 errors Matlab UMPFACK, Shape Functions T scaling femdoc UMFPACK, Shape Functions T scaling DSOL L2 Error Interface Error L2 Error Condition # Enrichment Interface Error L2 Error Condition # Enrichment Mean Error Maximum Minimum Std. Dev E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-05
39 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Shape Functions Maximum
40 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Shape Functions Maximum
41 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Spatial Derivatives Shape Functions Maximum Differences in L2 errors Time Matlab UMFPACK, Spatial Derivatives Shape Functions T scaling femdoc UMFPACK, Spatial Derivatives Shape Functions T scaling DSOL L2 Error Interface Error L2 Error Condition # Enrichment Interface Error L2 Error Condition # Enrichment Mean Error E E E E E E E E E E-05 Maximum Minimum Std. Dev E E E E E E E E E E E E E E E E E E E E E E E E E E E E E E-05
42 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Spatial Derivatives Shape Functions Maximum
43 Results Matlab comparison Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct Pre-conditioner: Spatial Derivatives Shape Functions Maximum
44 Results 2.5D Cylindrical inclusion Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: GMRES Aztec with ILU pre-conditioner Pre-conditioner: No pre-conditioner
45 Results 2.5D Cylindrical inclusion Mesh size: 30x30 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: GMRES Aztec with ILU pre-conditioner Pre-conditioner: No pre-conditioner
46 Results 3D Cylindrical inclusion Mesh size: 50x50 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: GMRES Aztec with ILU pre-conditioner Pre-conditioner: No pre-conditioner
47 Results 3D Cylindrical inclusion Mesh size: 50x50 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: GMRES Aztec with ILU pre-conditioner Pre-conditioner: No pre-conditioner
48 Results 3D Cylindrical inclusion Mesh size: 50x50 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: GMRES Aztec with ILU pre-conditioner Pre-conditioner: No pre-conditioner
49 Results: 2D niv50 thermal
50 Results: 2D niv50 thermal
51 Results: 2D niv50 thermal
52 Results: 2D niv50 thermal Symmetry in interface errors EXTREME VALUES (these extremes appeared in the old enrichment algorithm also): At iterations ~194 and ~305 (194 and 305 are identical) Location of double intersections
53 Results: 2D niv50 thermal Due to the method for computing double intersections and their main phase, the following case occurs: Small phase 1 diamond in center of intersections
54 Results
55 Results
56 Results: 3D niv50 thermal
57 Partial Results Stabilization factor sweep across mesh refinements Mesh size: 20x20, 30x30, 40x40, 50x50 Conductivity ratios: k2/k1 = 10 Stabilization factor: 10 Interface formulation: stabilized Lagrange Solver: Direct UMFPACK Pre-conditioner: Spatial Derivatives Shape Functions Maximum
58 Partial Results Stabilization factor sweep across mesh refinements
59 Questions?
60 Bibliography [1] Zienkiewicz, O.C.; Taylor, R.L.; Zhu, J.Z., Finite Element Method - Its Basis and Fundamentals (6th Edition). [2] Abdelaziz, Y., Hamouine, A., A survey of the extended finite element, Computers & Structures, Volume 86, Issues 11-12, Pages , [3] Fries, T.P., The extended Finite Element Method. < Accessed on September 21 st, [4] Stolarska M., Chopp D.L., Moës N., Belytschko T., Modeling crack growth by level sets in the extended finite element method, International Journal for Numerical Methods in Engineering, Volume 51, Issue 8, Pages , [5] Lee, D.T., Schachter, B.J., Two Algorithms for Constructing a Delaunay Triangulation, International Journal of Computer & Information Sciences, Volume 9, Issue 3, Pages , [6] Hansbo, P., Hansbo, A., A finite element method for the simulation of strong and weak discontinuities in solid mechanics, Computer Methods in Applied Mechanics and Engineering, Volume 193, Issue 33-35, Pages , [7] Joe, B., GEOMPACK3 Computational Geometry in 2D, 3D, ND. < Accessed on November 10 th, [8] Joe, B., Geompack3_prb.f90, a sample problem. < Accessed on November 10 th, [9] Flood Fill. Wikipedia, the free encyclopedia. < Accessed on November 29 th, Images used in this presentation:
A complete methodology for the implementation of XFEM inclusive models
University of Colorado Boulder A complete methodology for the implementation of XFEM inclusive models Author: Carlos Hernan Villanueva Kai Yu Center for Aerospace Structures Department of Mechanical Engineering
More information1 Exercise: Heat equation in 2-D with FE
1 Exercise: Heat equation in 2-D with FE Reading Hughes (2000, sec. 2.3-2.6 Dabrowski et al. (2008, sec. 1-3, 4.1.1, 4.1.3, 4.2.1 This FE exercise and most of the following ones are based on the MILAMIN
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 informationModeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures
University of California, Davis Modeling Discontinuities and their Evolution within Finite Elements: Application to Material Interfaces, 3-D Cracks, and Microstructures N. Sukumar UC Davis Rutgers University,
More informationJune 5, Institute of Structural Analysis Graz University of Technology Lessingstr. 25/II, 8010 Graz, Austria
Higher-order meshing of implicit geometries part I: Integration and interpolation in cut elements arxiv:706.00578v [cs.na] 2 Jun 207 T.P. Fries, S. Omerović, D. Schöllhammer, J. Steidl June 5, 207 Institute
More informationInvestigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models. David M. Trujillo 1. December 2005
Investigation of a Robust Method for Connecting Dissimilar 3D Finite Element Models by David M. Trujillo 1 December 2005 1 Consultant, TRUCOMP, Fountain Valley, California trucomp@earthlink.net Abstract
More informationLS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA.
12 th International LS-DYNA Users Conference FSI/ALE(1) LS-DYNA 980 : Recent Developments, Application Areas and Validation Process of the Incompressible fluid solver (ICFD) in LS-DYNA Part 1 Facundo Del
More informationLecture 2 Unstructured Mesh Generation
Lecture 2 Unstructured Mesh Generation MIT 16.930 Advanced Topics in Numerical Methods for Partial Differential Equations Per-Olof Persson (persson@mit.edu) February 13, 2006 1 Mesh Generation Given a
More information13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY
13.472J/1.128J/2.158J/16.940J COMPUTATIONAL GEOMETRY Lecture 23 Dr. W. Cho Prof. N. M. Patrikalakis Copyright c 2003 Massachusetts Institute of Technology Contents 23 F.E. and B.E. Meshing Algorithms 2
More informationCOMPUTER AIDED ENGINEERING. Part-1
COMPUTER AIDED ENGINEERING Course no. 7962 Finite Element Modelling and Simulation Finite Element Modelling and Simulation Part-1 Modeling & Simulation System A system exists and operates in time and space.
More informationContents. I The Basic Framework for Stationary Problems 1
page v Preface xiii I The Basic Framework for Stationary Problems 1 1 Some model PDEs 3 1.1 Laplace s equation; elliptic BVPs... 3 1.1.1 Physical experiments modeled by Laplace s equation... 5 1.2 Other
More informationImplementation of a XFEM toolbox in Diffpack
Implementation of a XFEM toolbox in Diffpack 1 Md Naim Hossain, 2 Daniel Alves Paladim, 1 Frank Vogel, 2 Prof. tephane P.A. Bordas 1 inutech GmbH, NuremBerg, Germany 2 Cardiff chool of Engineering, Institute
More informationUnstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications
Unstructured Mesh Generation for Implicit Moving Geometries and Level Set Applications Per-Olof Persson (persson@mit.edu) Department of Mathematics Massachusetts Institute of Technology http://www.mit.edu/
More informationGEOMETRY MODELING & GRID GENERATION
GEOMETRY MODELING & GRID GENERATION Dr.D.Prakash Senior Assistant Professor School of Mechanical Engineering SASTRA University, Thanjavur OBJECTIVE The objectives of this discussion are to relate experiences
More informationOutline. Level Set Methods. For Inverse Obstacle Problems 4. Introduction. Introduction. Martin Burger
For Inverse Obstacle Problems Martin Burger Outline Introduction Optimal Geometries Inverse Obstacle Problems & Shape Optimization Sensitivity Analysis based on Gradient Flows Numerical Methods University
More informationHeat Transfer Analysis of a Pipe
LESSON 25 Heat Transfer Analysis of a Pipe 3 Fluid 800 Ambient Temperture Temperture, C 800 500 2 Dia Fluid Ambient 10 20 30 40 Time, s Objectives: Transient Heat Transfer Analysis Model Convection, Conduction
More informationCHAPTER 1. Introduction
ME 475: Computer-Aided Design of Structures 1-1 CHAPTER 1 Introduction 1.1 Analysis versus Design 1.2 Basic Steps in Analysis 1.3 What is the Finite Element Method? 1.4 Geometrical Representation, Discretization
More informationSHAPE SENSITIVITY ANALYSIS INCLUDING QUALITY CONTROL WITH CARTESIAN FINITE ELEMENT MESHES
VI International Conference on Adaptive Modeling and Simulation ADMOS 2013 J. P. Moitinho de Almeida, P. Díez, C. Tiago and N. Parés (Eds) SHAPE SENSITIVITY ANALYSIS INCLUDING QUALITY CONTROL WITH CARTESIAN
More informationTheoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization. M2DO Lab 1,2. 1 Cardiff University
Theoretical Background for OpenLSTO v0.1: Open Source Level Set Topology Optimization M2DO Lab 1,2 1 Cardiff University 2 University of California, San Diego November 2017 A brief description of theory
More informationAdaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics
University of South Carolina Scholar Commons Theses and Dissertations 2015 Adaptivity For Meshfree Point Collocation Methods In Linear Elastic Solid Mechanics Joshua Wayne Derrick University of South Carolina
More informationLevel-set and ALE Based Topology Optimization Using Nonlinear Programming
10 th World Congress on Structural and Multidisciplinary Optimization May 19-24, 2013, Orlando, Florida, USA Level-set and ALE Based Topology Optimization Using Nonlinear Programming Shintaro Yamasaki
More informationCode_Aster. SSNV209 Interface in contact rubbing with X-FEM
Titre : SSNV209 - Interface en contact frottant avec X-FEM Date : 21/07/2015 Page : 1/34 SSNV209 Interface in contact rubbing with X-FEM Summary: This problem corresponds to a quasi-static analysis of
More informationSimilar Pulley Wheel Description J.E. Akin, Rice University
Similar Pulley Wheel Description J.E. Akin, Rice University The SolidWorks simulation tutorial on the analysis of an assembly suggested noting another type of boundary condition that is not illustrated
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 informationJassim Katwan, Osama (2010) Discrete modelling of heat transfer. MSc(R) thesis, University of Glasgow.
Jassim Katwan, Osama (2010) Discrete modelling of heat transfer. MSc(R) thesis, University of Glasgow. http://theses.gla.ac.uk/1602/ Copyright and moral rights for this thesis are retained by the author
More informationCrack Interaction Studies Using XFEM Technique
Journal of Solid Mechanics Vol. 6, No. 4 (214) pp. 41-421 Crack Interaction Studies Using XFEM Technique K. Sharma * Reactor Safety Divison, Bhabha Atomic Research Centre, Trombay, Mumbai, 485, India Received
More informationBarycentric Finite Element Methods
University of California, Davis Barycentric Finite Element Methods N. Sukumar University of California at Davis SIAM Conference on Geometric Design and Computing November 8, 2007 Collaborators and Acknowledgements
More informationAssignment in The Finite Element Method, 2017
Assignment in The Finite Element Method, 2017 Division of Solid Mechanics The task is to write a finite element program and then use the program to analyse aspects of a surface mounted resistor. The problem
More informationPractical Linear Algebra: A Geometry Toolbox
Practical Linear Algebra: A Geometry Toolbox Third edition Chapter 17: Breaking It Up: Triangles Gerald Farin & Dianne Hansford CRC Press, Taylor & Francis Group, An A K Peters Book www.farinhansford.com/books/pla
More informationCommunications in Applied Mathematics and Computational Science
Communications in Applied Mathematics and Computational Science Volume 1 No. 1 2006 A COMPARISON OF THE EXTENDED FINITE ELEMENT METHOD WITH THE IMMERSED INTERFACE METHOD FOR ELLIPTIC EQUATIONS WITH DISCONTINUOUS
More informationEffects of Solvers on Finite Element Analysis in COMSOL MULTIPHYSICS
Effects of Solvers on Finite Element Analysis in COMSOL MULTIPHYSICS Chethan Ravi B.R Dr. Venkateswaran P Corporate Technology - Research and Technology Center Siemens Technology and Services Private Limited
More informationSurgical Cutting on a Multimodal Object Representation
Surgical Cutting on a Multimodal Object Representation Lenka Jeřábková and Torsten Kuhlen Virtual Reality Group, RWTH Aachen University, 52074 Aachen Email: jerabkova@rz.rwth-aachen.de Abstract. In this
More informationAn Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems
An Interface-fitted Mesh Generator and Polytopal Element Methods for Elliptic Interface Problems Long Chen University of California, Irvine chenlong@math.uci.edu Joint work with: Huayi Wei (Xiangtan University),
More informationNodal Integration Technique in Meshless Method
IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 2320-334X, Volume 11, Issue 1 Ver. IV (Feb. 2014), PP 18-26 Nodal Integration Technique in Meshless Method Ahmed MJIDILA
More informationTorsional-lateral buckling large displacement analysis with a simple beam using Abaqus 6.10
Torsional-lateral buckling large displacement analysis with a simple beam using Abaqus 6.10 This document contains an Abaqus tutorial for performing a buckling analysis using the finite element program
More informationModeling Skills Thermal Analysis J.E. Akin, Rice University
Introduction Modeling Skills Thermal Analysis J.E. Akin, Rice University Most finite element analysis tasks involve utilizing commercial software, for which you do not have the source code. Thus, you need
More informationTopology optimization of heat conduction problems
Topology optimization of heat conduction problems Workshop on industrial design optimization for fluid flow 23 September 2010 Misha Marie Gregersen Anton Evgrafov Mads Peter Sørensen Technical University
More informationHow TMG Uses Elements and Nodes
Simulation: TMG Thermal Analysis User's Guide How TMG Uses Elements and Nodes Defining Boundary Conditions on Elements You create a TMG thermal model in exactly the same way that you create any finite
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 informationInfluence of geometric imperfections on tapered roller bearings life and performance
Influence of geometric imperfections on tapered roller bearings life and performance Rodríguez R a, Calvo S a, Nadal I b and Santo Domingo S c a Computational Simulation Centre, Instituto Tecnológico de
More informationThermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA
Thermal Coupling Method Between SPH Particles and Solid Elements in LS-DYNA INTRODUCTION: Jingxiao Xu, Jason Wang LSTC Heat transfer is very important in many industrial and geophysical problems. Many
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 informationChapter 7 Practical Considerations in Modeling. Chapter 7 Practical Considerations in Modeling
CIVL 7/8117 1/43 Chapter 7 Learning Objectives To present concepts that should be considered when modeling for a situation by the finite element method, such as aspect ratio, symmetry, natural subdivisions,
More informationModule: 2 Finite Element Formulation Techniques Lecture 3: Finite Element Method: Displacement Approach
11 Module: 2 Finite Element Formulation Techniques Lecture 3: Finite Element Method: Displacement Approach 2.3.1 Choice of Displacement Function Displacement function is the beginning point for the structural
More informationAn Explanation on Computation of Fracture Mechanics Parameters in ANSYS
University of Tennessee Space Institute Department of Mechanical, Aerospace & Biomedical Engineering Fracture Mechanics Course (ME 524) An Explanation on Computation of Fracture Mechanics Parameters in
More informationMultigrid Pattern. I. Problem. II. Driving Forces. III. Solution
Multigrid Pattern I. Problem Problem domain is decomposed into a set of geometric grids, where each element participates in a local computation followed by data exchanges with adjacent neighbors. The grids
More informationResearch Article Efficient Finite Element Methodology Based on Cartesian Grids: Application to Structural Shape Optimization
Abstract and Applied Analysis Volume 2013, Article ID 953786, 19 pages http://dx.doi.org/10.1155/2013/953786 Research Article Efficient Finite Element Methodology Based on Cartesian Grids: Application
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationcourse outline basic principles of numerical analysis, intro FEM
idealization, equilibrium, solutions, interpretation of results types of numerical engineering problems continuous vs discrete systems direct stiffness approach differential & variational formulation introduction
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationEffectiveness of Element Free Galerkin Method over FEM
Effectiveness of Element Free Galerkin Method over FEM Remya C R 1, Suji P 2 1 M Tech Student, Dept. of Civil Engineering, Sri Vellappaly Natesan College of Engineering, Pallickal P O, Mavelikara, Kerala,
More informationExecutive Summary Sefea Basic Theory
Executive Summary Sefea is one of the newest generations of enriched finite element methods. Developed specifically for low-order 4-node tetrahedron and 3-node triangle in the CAE environment, Sefea achieves
More informationAutomated Modelling of Evolving Discontinuities. Author to whom correspondence should be addressed;
Algorithms 2009, 2, 1008-1030; doi:10.3390/a2031008 Article OPEN ACCESS algorithms ISSN 1999-4893 www.mdpi.com/journal/algorithms Automated Modelling of Evolving Discontinuities Mehdi Nikbakht 1 and Garth
More informationAdvanced Surface Based MoM Techniques for Packaging and Interconnect Analysis
Electrical Interconnect and Packaging Advanced Surface Based MoM Techniques for Packaging and Interconnect Analysis Jason Morsey Barry Rubin, Lijun Jiang, Lon Eisenberg, Alina Deutsch Introduction Fast
More informationThe 3D DSC in Fluid Simulation
The 3D DSC in Fluid Simulation Marek K. Misztal Informatics and Mathematical Modelling, Technical University of Denmark mkm@imm.dtu.dk DSC 2011 Workshop Kgs. Lyngby, 26th August 2011 Governing Equations
More informationStability Analysis of Slope Based on FLAC3D Using Strength Reduction Method. Speaker: Tran, Thi Kim Tu Advisor: Chuen Fa Ni Date: 2017/03/02
Stability Analysis of Slope Based on FLAC3D Using Strength Reduction Method Speaker: Tran, Thi Kim Tu Advisor: Chuen Fa Ni Date: 2017/03/02 1. INTRODUCTION 2. METHODOLOGY 3. RESULTS AND DISCUSSION 4. CONCLUSION
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 informationStream Function-Vorticity CFD Solver MAE 6263
Stream Function-Vorticity CFD Solver MAE 66 Charles O Neill April, 00 Abstract A finite difference CFD solver was developed for transient, two-dimensional Cartesian viscous flows. Flow parameters are solved
More informationGenerative Part Structural Analysis Fundamentals
CATIA V5 Training Foils Generative Part Structural Analysis Fundamentals Version 5 Release 19 September 2008 EDU_CAT_EN_GPF_FI_V5R19 About this course Objectives of the course Upon completion of this course
More informationExample 13 - Shock Tube
Example 13 - Shock Tube Summary This famous experiment is interesting for observing the shock-wave propagation. Moreover, this case uses the representation of perfect gas and compares the different formulations:
More informationOverview and Recent Developments of Dynamic Mesh Capabilities
Overview and Recent Developments of Dynamic Mesh Capabilities Henrik Rusche and Hrvoje Jasak h.rusche@wikki-gmbh.de and h.jasak@wikki.co.uk Wikki Gmbh, Germany Wikki Ltd, United Kingdom 6th OpenFOAM Workshop,
More informationAMS527: Numerical Analysis II
AMS527: Numerical Analysis II A Brief Overview of Finite Element Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao SUNY Stony Brook AMS527: Numerical Analysis II 1 / 25 Overview Basic concepts Mathematical
More informationErrors Estimation of Different Numerical Integration Techniques on Stiffness Matrix and Stress Intensity Factor in XFEM
2013, Scienceline Publication Journal of Civil Engineering and Urbanism Volume 3, Issue 5: 265-278 (2013) (Received: June 20, 2013; Accepted: August 21, 2013; Published: September 30, 2013) ISSN-2252-0430
More informationSTRUCTURAL & MULTIDISCIPLINARY OPTIMIZATION
STRUCTURAL & MULTIDISCIPLINARY OPTIMIZATION Pierre DUYSINX Patricia TOSSINGS Department of Aerospace and Mechanical Engineering Academic year 2018-2019 1 Course objectives To become familiar with the introduction
More informationC. A. D. Fraga Filho 1,2, D. F. Pezzin 1 & J. T. A. Chacaltana 1. Abstract
Advanced Computational Methods and Experiments in Heat Transfer XIII 15 A numerical study of heat diffusion using the Lagrangian particle SPH method and the Eulerian Finite-Volume method: analysis of convergence,
More informationOverview of Unstructured Mesh Generation Methods
Overview of Unstructured Mesh Generation Methods Structured Meshes local mesh points and cells do not depend on their position but are defined by a general rule. Lead to very efficient algorithms and storage.
More informationAutomatic hybrid mesh generation for the boundary face method
Boundary Elements and Other Mesh Reduction Methods XXXVI 139 Automatic hybrid mesh generation for the boundary face method Cheng Huang & Jianming Zhang State Key Laboratory of Advanced Design and Manufacturing
More informationModeling Large Sliding Frictional Contact Along Non-Smooth Discontinuities in X-FEM
Modeling Large Sliding Frictional Contact Along Non-Smooth Discontinuities in X-FEM Seyed Mohammad Jafar TaheriMousavi and Seyedeh Mohadeseh Taheri Mousavi Abstract Modeling large frictional contact is
More informationITU/FAA Faculty of Aeronautics and Astronautics
S. Banu YILMAZ, Mehmet SAHIN, M. Fevzi UNAL, Istanbul Technical University, 34469, Maslak/Istanbul, TURKEY 65th Annual Meeting of the APS Division of Fluid Dynamics November 18-20, 2012, San Diego, CA
More informationhttp://miccom-center.org Topic: Continuum-Particle Simulation Software (COPSS-Hydrodynamics) Presenter: Jiyuan Li, The University of Chicago 2017 Summer School 1 What is Continuum-Particle Simulation?
More informationSTUDY OF FLOW PERFORMANCE OF A GLOBE VALVE AND DESIGN OPTIMISATION
Journal of Engineering Science and Technology Vol. 12, No. 9 (2017) 2403-2409 School of Engineering, Taylor s University STUDY OF FLOW PERFORMANCE OF A GLOBE VALVE AND DESIGN OPTIMISATION SREEKALA S. K.
More informationACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION
Scientific Research of the Institute of Mathematics and Computer Science ACCURACY OF NUMERICAL SOLUTION OF HEAT DIFFUSION EQUATION Ewa Węgrzyn-Skrzypczak, Tomasz Skrzypczak Institute of Mathematics, Czestochowa
More informationTowards a Material Point Method with Powell-Sabin spline basis functions
Towards a Material Point Method with Powell-Sabin spline basis functions MPM with higher-order C 1 -continuous basis functions on triangulations Pascal B.J. de Koster Delft University of Technology Towards
More informationThe Finite Element Method
The Finite Element Method A Practical Course G. R. Liu and S. S. Quek Chapter 1: Computational modeling An overview 1 CONTENTS INTRODUCTION PHYSICAL PROBLEMS IN ENGINEERING COMPUTATIONAL MODELLING USING
More informationCHAPTER-10 DYNAMIC SIMULATION USING LS-DYNA
DYNAMIC SIMULATION USING LS-DYNA CHAPTER-10 10.1 Introduction In the past few decades, the Finite Element Method (FEM) has been developed into a key indispensable technology in the modeling and simulation
More informationStudy of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket
RESEARCH ARTICLE OPEN ACCESS Study of Convergence of Results in Finite Element Analysis of a Plane Stress Bracket Gowtham K L*, Shivashankar R. Srivatsa** *(Department of Mechanical Engineering, B. M.
More informationEvolutionary topology optimization extended finite element method and isolines
Evolutionary topology optimization extended finite element method and isolines Meisam Abdi, Ricky Wildman and Ian Ashcroft Faculty of Engineering, University of Nottingham, University Park, Nottingham,
More informationAn explicit feature control approach in structural topology optimization
th World Congress on Structural and Multidisciplinary Optimisation 07 th -2 th, June 205, Sydney Australia An explicit feature control approach in structural topology optimization Weisheng Zhang, Xu Guo
More informationAnalysis Comparison between CFD and FEA of an Idealized Concept V- Hull Floor Configuration in Two Dimensions
2010 NDIA GROUND VEHICLE SYSTEMS ENGINEERING AND TECHNOLOGY SYMPOSIUM MODELING & SIMULATION, TESTING AND VALIDATION (MSTV) MINI-SYMPOSIUM AUGUST 17-19 DEARBORN, MICHIGAN Analysis Comparison between CFD
More informationGeometric Representations. Stelian Coros
Geometric Representations Stelian Coros Geometric Representations Languages for describing shape Boundary representations Polygonal meshes Subdivision surfaces Implicit surfaces Volumetric models Parametric
More informationA Comparison of the Computational Speed of 3DSIM versus ANSYS Finite Element Analyses for Simulation of Thermal History in Metal Laser Sintering
A Comparison of the Computational Speed of 3DSIM versus ANSYS Finite Element Analyses for Simulation of Thermal History in Metal Laser Sintering Kai Zeng a,b, Chong Teng a,b, Sally Xu b, Tim Sublette b,
More informationModule 3 Mesh Generation
Module 3 Mesh Generation 1 Lecture 3.1 Introduction 2 Mesh Generation Strategy Mesh generation is an important pre-processing step in CFD of turbomachinery, quite analogous to the development of solid
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 informationFinal Report. Discontinuous Galerkin Compressible Euler Equation Solver. May 14, Andrey Andreyev. Adviser: Dr. James Baeder
Final Report Discontinuous Galerkin Compressible Euler Equation Solver May 14, 2013 Andrey Andreyev Adviser: Dr. James Baeder Abstract: In this work a Discontinuous Galerkin Method is developed for compressible
More informationSIMULATION OF AN IMPLANTED PIFA FOR A CARDIAC PACEMAKER WITH EFIELD FDTD AND HYBRID FDTD-FEM
1 SIMULATION OF AN IMPLANTED PIFA FOR A CARDIAC PACEMAKER WITH EFIELD FDTD AND HYBRID FDTD- Introduction Medical Implanted Communication Service (MICS) has received a lot of attention recently. The MICS
More informationX-FEM based modelling of complex mixed mode fatigue crack propagation
X-FEM based modelling of complex mixed mode fatigue crack propagation Hans Minnebo 1, Simon André 2, Marc Duflot 1, Thomas Pardoen 2, Eric Wyart 1 1 Cenaero, Rue des Frères Wright 29, 6041 Gosselies, Belgium
More informationFoundations of Analytical and Numerical Field Computation
Foundations of Analytical and Numerical Field Computation Stephan Russenschuck, CERN-AT-MEL Stephan Russenschuck CERN, TE-MCS, 1211 Geneva, Switzerland 1 Permanent Magnet Circuits 2 Rogowski profiles Pole
More informationA New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes
A New Smoothing Algorithm for Quadrilateral and Hexahedral Meshes Sanjay Kumar Khattri Department of Mathematics, University of Bergen, Norway sanjay@mi.uib.no http://www.mi.uib.no/ sanjay Abstract. Mesh
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationParallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of Earth s Mantle
ICES Student Forum The University of Texas at Austin, USA November 4, 204 Parallel High-Order Geometric Multigrid Methods on Adaptive Meshes for Highly Heterogeneous Nonlinear Stokes Flow Simulations of
More informationDevelopment of a particle difference scheme for dynamic crack propagation problem
Development of a particle difference scheme for dynamic crack propagation problem *Young-Choel Yoon 1) and Kyeong-Hwan Kim ) 1) Dept. of Civil Engineering, Myongji College, Seoul 10-776, Korea ) Dept.
More informationDevelopment of mesh refinement methods at CFD codes for computational fluid mechanics problems DIPLOMA THESIS
NATIONAL TECHNICAL UNIVERSITY OF ATHENS DEPARTMENT OF MECHANICAL ENGINEERING SECTION OF FLUID MECHANICS Development of mesh refinement methods at CFD codes for computational fluid mechanics problems DIPLOMA
More informationOutline of the presentation
Surface Reconstruction Petra Surynková Charles University in Prague Faculty of Mathematics and Physics petra.surynkova@mff.cuni.cz Outline of the presentation My work up to now Surfaces of Building Practice
More informationVisualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium
Visualization of errors of finite element solutions P. Beckers, H.G. Zhong, Ph. Andry Aerospace Department, University of Liege, B-4000 Liege, Belgium Abstract The aim of this paper is to show how to use
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 informationEssay 5 Tutorial for a Three-Dimensional Heat Conduction Problem Using ANSYS
Essay 5 Tutorial for a Three-Dimensional Heat Conduction Problem Using ANSYS 5.1 Introduction The problem selected to illustrate the use of ANSYS software for a three-dimensional steadystate heat conduction
More informationAnsoft HFSS Mesh Menu
Ansoft HFSS After you have seeded your object, you must create the finite element mesh from which the variables and values of your model will be computed. The Mesh menu allows you to: Create or delete
More informationPhysically-Based Modeling and Animation. University of Missouri at Columbia
Overview of Geometric Modeling Overview 3D Shape Primitives: Points Vertices. Curves Lines, polylines, curves. Surfaces Triangle meshes, splines, subdivision surfaces, implicit surfaces, particles. Solids
More informationA methodology for the rigorous verification of plasma simulation codes
A methodology for the rigorous verification of plasma simulation codes Fabio Riva P. Ricci, C. Beadle, F.D. Halpern, S. Jolliet, J. Loizu, J. Morales, A. Mosetto, P. Paruta, C. Wersal École Polytechnique
More informationCHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING
113 CHAPTER 6 EXPERIMENTAL AND FINITE ELEMENT SIMULATION STUDIES OF SUPERPLASTIC BOX FORMING 6.1 INTRODUCTION Superplastic properties are exhibited only under a narrow range of strain rates. Hence, it
More information