No Cover Image

E-Thesis 537 views 190 downloads

A vertex centred finite volume method for solid dynamics. / Miquel Aguirre Font

Swansea University Author: Miquel Aguirre Font

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:29.0413813</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><swanseaauthors><author><sid>e227e8797216b730d8a39ca4f6cf1bb1</sid><ORCID>NULL</ORCID><firstname>Miquel</firstname><surname>Aguirre Font</surname><name>Miquel Aguirre Font</name><active>true</active><ethesisStudent>true</ethesisStudent></author></swanseaauthors><date>2018-08-02</date><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/><journalNumber></journalNumber><paginationStart/><paginationEnd/><publisher/><placeOfPublication/><isbnPrint/><issnPrint/><issnElectronic/><keywords>Computer science.;Mechanics.</keywords><publishedDay>31</publishedDay><publishedMonth>12</publishedMonth><publishedYear>2014</publishedYear><publishedDate>2014-12-31</publishedDate><doi/><url/><notes/><college>COLLEGE NANME</college><department>Engineering</department><CollegeCode>COLLEGE CODE</CollegeCode><institution>Swansea University</institution><degreelevel>Doctoral</degreelevel><degreename>Ph.D</degreename><apcterm/><lastEdited>2018-08-02T16:24:29.0413813</lastEdited><Created>2018-08-02T16:24:29.0413813</Created><path><level id="1">Faculty of Science and Engineering</level><level id="2">School of Engineering and Applied Sciences - Uncategorised</level></path><authors><author><firstname>Miquel</firstname><surname>Aguirre Font</surname><orcid>NULL</orcid><order>1</order></author></authors><documents><document><filename>0042384-02082018162450.pdf</filename><originalFilename>10798092.pdf</originalFilename><uploaded>2018-08-02T16:24:50.1800000</uploaded><type>Output</type><contentLength>13333377</contentLength><contentType>application/pdf</contentType><version>E-Thesis</version><cronfaStatus>true</cronfaStatus><embargoDate>2018-08-02T16:24:50.1800000</embargoDate><copyrightCorrect>false</copyrightCorrect></document></documents><OutputDurs/></rfc1807>
spelling 2018-08-02T16:24:29.0413813 v2 42384 2018-08-02 A vertex centred finite volume method for solid dynamics. e227e8797216b730d8a39ca4f6cf1bb1 NULL Miquel Aguirre Font Miquel Aguirre Font true true 2018-08-02 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. 31 12 2014 2014-12-31 COLLEGE NANME Engineering COLLEGE CODE Swansea University Doctoral Ph.D 2018-08-02T16:24:29.0413813 2018-08-02T16:24:29.0413813 Faculty of Science and Engineering School of Engineering and Applied Sciences - Uncategorised Miquel Aguirre Font NULL 1 0042384-02082018162450.pdf 10798092.pdf 2018-08-02T16:24:50.1800000 Output 13333377 application/pdf E-Thesis true 2018-08-02T16:24:50.1800000 false
title A vertex centred finite volume method for solid dynamics.
spellingShingle A vertex centred finite volume method for solid dynamics.
Miquel Aguirre Font
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 e227e8797216b730d8a39ca4f6cf1bb1
author_id_fullname_str_mv e227e8797216b730d8a39ca4f6cf1bb1_***_Miquel Aguirre Font
author Miquel Aguirre Font
author2 Miquel Aguirre Font
format E-Thesis
publishDate 2014
institution Swansea University
college_str Faculty of Science and Engineering
hierarchytype
hierarchy_top_id facultyofscienceandengineering
hierarchy_top_title Faculty of Science and Engineering
hierarchy_parent_id facultyofscienceandengineering
hierarchy_parent_title Faculty of Science and Engineering
department_str School of Engineering and Applied Sciences - Uncategorised{{{_:::_}}}Faculty of Science and Engineering{{{_:::_}}}School of Engineering and Applied Sciences - Uncategorised
document_store_str 1
active_str 0
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-12-31T03:52:51Z
_version_ 1763752617616343040
score 11.012678

Similar Items