No Cover Image

E-Thesis 64 views 24 downloads

A vertex centred finite volume method for solid dynamics. /

Swansea University Author: Aguirre Font, Miquel

Abstract

With the increase of the computational power over the last decades, computational solid dynamics has become a major field of interest in many industrial applications (as for example aerospace, automotive industry, biomedical engineering or manufacturing). Traditionally, these types of simulations ha...

Full description

Published: 2014
Institution: Swansea University
Degree level: Doctoral
Degree name: Ph.D
URI: https://cronfa.swan.ac.uk/Record/cronfa42384
Tags: Add Tag
No Tags, Be the first to tag this record!
first_indexed 2018-08-02T18:54:34Z
last_indexed 2018-08-03T10:10:00Z
id cronfa42384
recordtype RisThesis
fullrecord <?xml version="1.0"?><rfc1807><datestamp>2018-08-02T16:24:50Z</datestamp><bib-version>v2</bib-version><id>42384</id><entry>2018-08-02</entry><title>A vertex centred finite volume method for solid dynamics.</title><alternativeTitle></alternativeTitle><author>Miquel Aguirre Font</author><firstname>Miquel</firstname><surname>Aguirre Font</surname><active>true</active><ORCID/><ethesisStudent>true</ethesisStudent><sid>89fa6381ba826051ade4a373eba5d567</sid><email>378f6a0f803397ec2363ff574d3417ae</email><emailaddr>FyFy4YJxPncboyv9Cm5l2NyvqZQRJmUl2lxhnzSZE7o=</emailaddr><date>2018-08-02</date><deptcode>EEN</deptcode><abstract>With the increase of the computational power over the last decades, computational solid dynamics has become a major field of interest in many industrial applications (as for example aerospace, automotive industry, biomedical engineering or manufacturing). Traditionally, these types of simulations have been carried out using the Finite Element Method (FEM) in conjunction with displacement based formulations, where the displacements are treated as the main problem variables. In the context of solid mechanics, the use of low order (linear) tetrahedral (or triangular) elements is always preferred, primarily due to the complexity of the constitutive models involved (i.e., the number of evaluations needs to be kept as low as possible) as well as the automatic tetrahedral (or triangular) mesh generators available. However, the combination of low order elements with FEM displacement based formulations presents a series of shortcomings, namely: one order of accuracy less for stresses (or strains) than for displacements (or velocities), poor behaviour in bending dominated scenarios, volumetric and shear locking or the appearance of spurious pressure modes. Furthermore, the time integration is usually performed using Newmark integrators, which tend to introduce high frequency noise in the vicinity of sharp gradients. Recently, a new Lagrangian mixed methodology has been presented for the simulation of fast transient dynamics problems. This methodology is in the form of a system of first order conservation laws, where the linear momentum, p, and the deformation gradient tensor, F, are regarded as the two main conservation variables. When thermo-mechanical constitutive models are involved, the formulation is complemented with the first law of thermodynamics (conservation of energy). It has been proven that this formulation circumvents the drawbacks of the low order displacement based FEM methodologies mentioned above. The formulation, presented as a set of conservation laws, allows for standard Computational Fluid Dynamcis (CFD) spatial discretisations. So far, successful implementations have been carried out using cell centred upwind Finite Volume method, two-step Taylor Galerkin, Finite Element Petrov Galerkin (PG), and Hybridizable Discontinuous Galerkin (HDG). The objective of this thesis is to present a new spatial discretisation in order to solve large scale real life problems. To do so, the Jameson-Schmidt-Turkel (JST) scheme, widely know within the CFD community, will be chosen. The JST scheme is a vertex centred finite volume, that combines the use of central differences with an artificial dissipation term. The scheme obtains second order spatial accuracy without the need of linear reconstruction. Furthermore, the artificial dissipation term includes a shock capturing sensor, very suitable in the context of fast dynamics. The scheme can be implemented in an edge-based framework, which combined with the vertex centred storage of the variables results into a computationally efficient scheme. The JST spatial discretisation will be combined with a Total Variation Diminishing (TVD) two-stage Runge-Kutta time integrator. These spatial and temporal discretisations will be adapted to the problem at hand. Specifically, compatibility conditions (involutions) will have to be satisfied by the discrete scheme. Furthermore, numerical corrections will be introduced in order to ensure the conservation of linear and angular momenta. The framework results in a low order computationally efficient solver for solid dynamics, which proves to be very competitive in nearly incompressible scenarios and bending dominated applications. The thesis will present numerical results for one dimension, two dimensions (triangular meshes) and three dimensions (tetrahedral meshes). The problems are chosen in order to prove the order of accuracy, robustness and conservation properties of the algorithm.</abstract><type>E-Thesis</type><journal></journal><volume/><journalNumber></journalNumber><paginationStart></paginationStart><paginationEnd></paginationEnd><publisher></publisher><placeOfPublication></placeOfPublication><isbnPrint></isbnPrint><isbnElectronic/><issnPrint></issnPrint><issnElectronic></issnElectronic><keywords>Computer science.;Mechanics.</keywords><publishedDay>0</publishedDay><publishedMonth>0</publishedMonth><publishedYear>2014</publishedYear><publishedDate>2014-01-01</publishedDate><doi></doi><url></url><notes></notes><college>College of Engineering</college><department>Engineering</department><CollegeCode>CENG</CollegeCode><DepartmentCode>EEN</DepartmentCode><institution>Swansea University</institution><researchGroup></researchGroup><supervisor></supervisor><sponsorsfunders/><grantnumber/><degreelevel>Doctoral</degreelevel><degreename>Ph.D</degreename><lastEdited>2018-08-02T16:24:29Z</lastEdited><Created>2018-08-02T16:24:29Z</Created><path><level id="1">College of Engineering</level><level id="2">Engineering</level></path><authors/><documents><document><filename>0042384-02082018162450.pdf</filename><originalFilename>10798092.pdf</originalFilename><uploaded>2018-08-02T16:24:50Z</uploaded><type>Output</type><contentLength>123456</contentLength><contentType>application/pdf</contentType><version>E-Thesis</version><cronfaStatus>true</cronfaStatus><action>Published to Cronfa</action><actionDate>02/08/2018</actionDate><embargoDate>2018-08-02T16:24:50.18</embargoDate><documentNotes/><copyrightCorrect>false</copyrightCorrect><language/></document></documents></rfc1807>
spelling 2018-08-02T16:24:50Z v2 42384 2018-08-02 A vertex centred finite volume method for solid dynamics. Miquel Aguirre Font Miquel Aguirre Font true true 89fa6381ba826051ade4a373eba5d567 378f6a0f803397ec2363ff574d3417ae FyFy4YJxPncboyv9Cm5l2NyvqZQRJmUl2lxhnzSZE7o= 2018-08-02 EEN With the increase of the computational power over the last decades, computational solid dynamics has become a major field of interest in many industrial applications (as for example aerospace, automotive industry, biomedical engineering or manufacturing). Traditionally, these types of simulations have been carried out using the Finite Element Method (FEM) in conjunction with displacement based formulations, where the displacements are treated as the main problem variables. In the context of solid mechanics, the use of low order (linear) tetrahedral (or triangular) elements is always preferred, primarily due to the complexity of the constitutive models involved (i.e., the number of evaluations needs to be kept as low as possible) as well as the automatic tetrahedral (or triangular) mesh generators available. However, the combination of low order elements with FEM displacement based formulations presents a series of shortcomings, namely: one order of accuracy less for stresses (or strains) than for displacements (or velocities), poor behaviour in bending dominated scenarios, volumetric and shear locking or the appearance of spurious pressure modes. Furthermore, the time integration is usually performed using Newmark integrators, which tend to introduce high frequency noise in the vicinity of sharp gradients. Recently, a new Lagrangian mixed methodology has been presented for the simulation of fast transient dynamics problems. This methodology is in the form of a system of first order conservation laws, where the linear momentum, p, and the deformation gradient tensor, F, are regarded as the two main conservation variables. When thermo-mechanical constitutive models are involved, the formulation is complemented with the first law of thermodynamics (conservation of energy). It has been proven that this formulation circumvents the drawbacks of the low order displacement based FEM methodologies mentioned above. The formulation, presented as a set of conservation laws, allows for standard Computational Fluid Dynamcis (CFD) spatial discretisations. So far, successful implementations have been carried out using cell centred upwind Finite Volume method, two-step Taylor Galerkin, Finite Element Petrov Galerkin (PG), and Hybridizable Discontinuous Galerkin (HDG). The objective of this thesis is to present a new spatial discretisation in order to solve large scale real life problems. To do so, the Jameson-Schmidt-Turkel (JST) scheme, widely know within the CFD community, will be chosen. The JST scheme is a vertex centred finite volume, that combines the use of central differences with an artificial dissipation term. The scheme obtains second order spatial accuracy without the need of linear reconstruction. Furthermore, the artificial dissipation term includes a shock capturing sensor, very suitable in the context of fast dynamics. The scheme can be implemented in an edge-based framework, which combined with the vertex centred storage of the variables results into a computationally efficient scheme. The JST spatial discretisation will be combined with a Total Variation Diminishing (TVD) two-stage Runge-Kutta time integrator. These spatial and temporal discretisations will be adapted to the problem at hand. Specifically, compatibility conditions (involutions) will have to be satisfied by the discrete scheme. Furthermore, numerical corrections will be introduced in order to ensure the conservation of linear and angular momenta. The framework results in a low order computationally efficient solver for solid dynamics, which proves to be very competitive in nearly incompressible scenarios and bending dominated applications. The thesis will present numerical results for one dimension, two dimensions (triangular meshes) and three dimensions (tetrahedral meshes). The problems are chosen in order to prove the order of accuracy, robustness and conservation properties of the algorithm. E-Thesis Computer science.;Mechanics. 0 0 2014 2014-01-01 College of Engineering Engineering CENG EEN Swansea University Doctoral Ph.D 2018-08-02T16:24:29Z 2018-08-02T16:24:29Z College of Engineering Engineering 0042384-02082018162450.pdf 10798092.pdf 2018-08-02T16:24:50Z Output 123456 application/pdf E-Thesis true Published to Cronfa 02/08/2018 2018-08-02T16:24:50.18 false
title A vertex centred finite volume method for solid dynamics.
spellingShingle A vertex centred finite volume method for solid dynamics.
Aguirre Font, Miquel
title_short A vertex centred finite volume method for solid dynamics.
title_full A vertex centred finite volume method for solid dynamics.
title_fullStr A vertex centred finite volume method for solid dynamics.
title_full_unstemmed A vertex centred finite volume method for solid dynamics.
title_sort A vertex centred finite volume method for solid dynamics.
author_id_str_mv 89fa6381ba826051ade4a373eba5d567
author_id_fullname_str_mv 89fa6381ba826051ade4a373eba5d567_***_Aguirre Font, Miquel
author Aguirre Font, Miquel
author2
format E-Thesis
publishDate 2014
institution Swansea University
college_str College of Engineering
hierarchytype
hierarchy_top_id collegeofengineering
hierarchy_top_title College of Engineering
hierarchy_parent_id collegeofengineering
hierarchy_parent_title College of Engineering
department_str Engineering{{{_:::_}}}College of Engineering{{{_:::_}}}Engineering
document_store_str 1
active_str 1
description With the increase of the computational power over the last decades, computational solid dynamics has become a major field of interest in many industrial applications (as for example aerospace, automotive industry, biomedical engineering or manufacturing). Traditionally, these types of simulations have been carried out using the Finite Element Method (FEM) in conjunction with displacement based formulations, where the displacements are treated as the main problem variables. In the context of solid mechanics, the use of low order (linear) tetrahedral (or triangular) elements is always preferred, primarily due to the complexity of the constitutive models involved (i.e., the number of evaluations needs to be kept as low as possible) as well as the automatic tetrahedral (or triangular) mesh generators available. However, the combination of low order elements with FEM displacement based formulations presents a series of shortcomings, namely: one order of accuracy less for stresses (or strains) than for displacements (or velocities), poor behaviour in bending dominated scenarios, volumetric and shear locking or the appearance of spurious pressure modes. Furthermore, the time integration is usually performed using Newmark integrators, which tend to introduce high frequency noise in the vicinity of sharp gradients. Recently, a new Lagrangian mixed methodology has been presented for the simulation of fast transient dynamics problems. This methodology is in the form of a system of first order conservation laws, where the linear momentum, p, and the deformation gradient tensor, F, are regarded as the two main conservation variables. When thermo-mechanical constitutive models are involved, the formulation is complemented with the first law of thermodynamics (conservation of energy). It has been proven that this formulation circumvents the drawbacks of the low order displacement based FEM methodologies mentioned above. The formulation, presented as a set of conservation laws, allows for standard Computational Fluid Dynamcis (CFD) spatial discretisations. So far, successful implementations have been carried out using cell centred upwind Finite Volume method, two-step Taylor Galerkin, Finite Element Petrov Galerkin (PG), and Hybridizable Discontinuous Galerkin (HDG). The objective of this thesis is to present a new spatial discretisation in order to solve large scale real life problems. To do so, the Jameson-Schmidt-Turkel (JST) scheme, widely know within the CFD community, will be chosen. The JST scheme is a vertex centred finite volume, that combines the use of central differences with an artificial dissipation term. The scheme obtains second order spatial accuracy without the need of linear reconstruction. Furthermore, the artificial dissipation term includes a shock capturing sensor, very suitable in the context of fast dynamics. The scheme can be implemented in an edge-based framework, which combined with the vertex centred storage of the variables results into a computationally efficient scheme. The JST spatial discretisation will be combined with a Total Variation Diminishing (TVD) two-stage Runge-Kutta time integrator. These spatial and temporal discretisations will be adapted to the problem at hand. Specifically, compatibility conditions (involutions) will have to be satisfied by the discrete scheme. Furthermore, numerical corrections will be introduced in order to ensure the conservation of linear and angular momenta. The framework results in a low order computationally efficient solver for solid dynamics, which proves to be very competitive in nearly incompressible scenarios and bending dominated applications. The thesis will present numerical results for one dimension, two dimensions (triangular meshes) and three dimensions (tetrahedral meshes). The problems are chosen in order to prove the order of accuracy, robustness and conservation properties of the algorithm.
published_date 2014-01-01T22:15:10Z
_version_ 1639341737407152128
score 10.826822