No Cover Image

E-Thesis 180 views 87 downloads

Development of a parallel CFD solver with application to arterial flows. / AMARPAL KAPOOR

Swansea University Author: AMARPAL KAPOOR

Abstract

In this research, the finite element method (FEM) was used to solve the nonlinear, incompressible, transient, three dimensional Navier-Stokes equations in their non-conservative form. Linear tetrahedron elements were employed with the elegant, equal order interpolation for both pressure and velocity...

Full description

Published: 2014
Institution: Swansea University
Degree level: Doctoral
Degree name: Ph.D
URI: https://cronfa.swan.ac.uk/Record/cronfa42216
Tags: Add Tag
No Tags, Be the first to tag this record!
first_indexed 2018-08-02T18:54:10Z
last_indexed 2020-09-04T03:03:02Z
id cronfa42216
recordtype RisThesis
fullrecord <?xml version="1.0"?><rfc1807><datestamp>2020-09-03T08:59:18.7915549</datestamp><bib-version>v2</bib-version><id>42216</id><entry>2018-08-02</entry><title>Development of a parallel CFD solver with application to arterial flows.</title><swanseaauthors><author><sid>aa2592aa0bb04bdab92450f24643bf5a</sid><firstname>AMARPAL</firstname><surname>KAPOOR</surname><name>AMARPAL KAPOOR</name><active>true</active><ethesisStudent>false</ethesisStudent></author></swanseaauthors><date>2018-08-02</date><abstract>In this research, the finite element method (FEM) was used to solve the nonlinear, incompressible, transient, three dimensional Navier-Stokes equations in their non-conservative form. Linear tetrahedron elements were employed with the elegant, equal order interpolation for both pressure and velocity. The characteristic based split scheme was formulated in a fully implicit manner to circumvent the time step restrictions of the classical explicit formulations. The monolithic (single step, fully coupled solution procedure for pressures and velocity) form of the CBS scheme was also derived and its suitability was positively demonstrated. Casting the CBS scheme in a monolithic framework, results in the generation of a pressure stabilization term in the mass conservation equation, thereby circumventing the LBB restriction by the elimination of the zero pressure block. An account of all the steps involved in discretizing the Navier-Stokes equations (both in split and monolithic frameworks) was presented in meticulous detail, which included the derivation of the convective and pressure stabilization terms, linearization of the non-linear terms and the consequent derivation of the highly efficient analytical jacobian matrix, along with the temporal and spatial discretizations of the corresponding terms. The monolithic and the split version of the CBS scheme were integrated into a parallel, scalable and extensible Fortran90 software called IFENs. The development of IFENs started during the course of this research and all of its components have been designed and implemented by the author of this thesis. Multi processor parallelism was achieved using the Intel implementation of the most widely used/preferred, Message Passing Interface (MPI) standard. The parallel support needed for the use of a variety of parallel, linear, iterative solvers belonging to the Krylov subspace family (e.g. GMRES and its variants, CG, BiCG, BiCG- stab, etc.), parallel non linear solvers belonging to the Newton-Krylov family (line search newton, trust region newton, nonlinear GMRES, etc.) and parallel preconditioners (incomplete LU, Additive Shwarz Method - ASM, algebraic multigrid, etc.), was provided by the incorporation of PETSc into IFENs. PETSc is a state of the art, non-trivial toolkit, which represents a collection of several parallel libraries useful in high performance scientific computing. Keeping in mind the specific requirements of IFENs, a custom mesh partitioner was implemented. It operated on meshes that were renumbered using bandwidth reducing algorithms like Revere Cuthill Mckee. The possibility of using established domain decomposition libraries like ParMETIS was explored and demonstrated to be counter productive for the demands of this research. After the preliminary testing and validation of the procedures adopted before and during the execution of IFENs, large, high definition domains representative of human arteries (specifically, carotid bifurcations, found in the neck) were considered and the complete incompressible set of Navier-Stokes equations were solved for pressure and velocity fields. During the tenure of this research more than 1000 recorded parallel test cases were executed to test various components of IFENs, as well as various simulations representative of a wide variety of problems. IFENs can easily handle meshes with tens of millions of elements. The largest mesh used for the purpose of this research contained 14.58 million tetrahedrons and 2.489 million nodes, which on average required just 7 minutes per timestep, while executing the classical split framework of the CBS scheme. Results from the simulation of 9 carotid meshes, representative of 4 carotid geometries were presented and found to be in good agreement with the available ultrasound data. The flow fields were analysed and post processed using different techniques for each case. The haemodynamic wall parameters like time averaged wall shear stress and oscillatory shear index were calculated and mapped onto the corresponding boundary nodes. The region in the carotid bifurcation susceptible to the deposition of plaques and consequent stenosis were pointed out and other anomalies were highlighted.</abstract><type>E-Thesis</type><journal/><publisher/><keywords>Computer science.;Biomedical engineering.</keywords><publishedDay>31</publishedDay><publishedMonth>12</publishedMonth><publishedYear>2014</publishedYear><publishedDate>2014-12-31</publishedDate><doi/><url/><notes/><college>COLLEGE NANME</college><CollegeCode>COLLEGE CODE</CollegeCode><institution>Swansea University</institution><degreelevel>Doctoral</degreelevel><degreename>Ph.D</degreename><apcterm/><lastEdited>2020-09-03T08:59:18.7915549</lastEdited><Created>2018-08-02T16:24:28.4641899</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>AMARPAL</firstname><surname>KAPOOR</surname><order>1</order></author></authors><documents><document><filename>0042216-02082018162437.pdf</filename><originalFilename>10797918.pdf</originalFilename><uploaded>2018-08-02T16:24:37.2170000</uploaded><type>Output</type><contentLength>18181210</contentLength><contentType>application/pdf</contentType><version>E-Thesis</version><cronfaStatus>true</cronfaStatus><embargoDate>2018-08-02T00:00:00.0000000</embargoDate><copyrightCorrect>false</copyrightCorrect></document></documents><OutputDurs/></rfc1807>
spelling 2020-09-03T08:59:18.7915549 v2 42216 2018-08-02 Development of a parallel CFD solver with application to arterial flows. aa2592aa0bb04bdab92450f24643bf5a AMARPAL KAPOOR AMARPAL KAPOOR true false 2018-08-02 In this research, the finite element method (FEM) was used to solve the nonlinear, incompressible, transient, three dimensional Navier-Stokes equations in their non-conservative form. Linear tetrahedron elements were employed with the elegant, equal order interpolation for both pressure and velocity. The characteristic based split scheme was formulated in a fully implicit manner to circumvent the time step restrictions of the classical explicit formulations. The monolithic (single step, fully coupled solution procedure for pressures and velocity) form of the CBS scheme was also derived and its suitability was positively demonstrated. Casting the CBS scheme in a monolithic framework, results in the generation of a pressure stabilization term in the mass conservation equation, thereby circumventing the LBB restriction by the elimination of the zero pressure block. An account of all the steps involved in discretizing the Navier-Stokes equations (both in split and monolithic frameworks) was presented in meticulous detail, which included the derivation of the convective and pressure stabilization terms, linearization of the non-linear terms and the consequent derivation of the highly efficient analytical jacobian matrix, along with the temporal and spatial discretizations of the corresponding terms. The monolithic and the split version of the CBS scheme were integrated into a parallel, scalable and extensible Fortran90 software called IFENs. The development of IFENs started during the course of this research and all of its components have been designed and implemented by the author of this thesis. Multi processor parallelism was achieved using the Intel implementation of the most widely used/preferred, Message Passing Interface (MPI) standard. The parallel support needed for the use of a variety of parallel, linear, iterative solvers belonging to the Krylov subspace family (e.g. GMRES and its variants, CG, BiCG, BiCG- stab, etc.), parallel non linear solvers belonging to the Newton-Krylov family (line search newton, trust region newton, nonlinear GMRES, etc.) and parallel preconditioners (incomplete LU, Additive Shwarz Method - ASM, algebraic multigrid, etc.), was provided by the incorporation of PETSc into IFENs. PETSc is a state of the art, non-trivial toolkit, which represents a collection of several parallel libraries useful in high performance scientific computing. Keeping in mind the specific requirements of IFENs, a custom mesh partitioner was implemented. It operated on meshes that were renumbered using bandwidth reducing algorithms like Revere Cuthill Mckee. The possibility of using established domain decomposition libraries like ParMETIS was explored and demonstrated to be counter productive for the demands of this research. After the preliminary testing and validation of the procedures adopted before and during the execution of IFENs, large, high definition domains representative of human arteries (specifically, carotid bifurcations, found in the neck) were considered and the complete incompressible set of Navier-Stokes equations were solved for pressure and velocity fields. During the tenure of this research more than 1000 recorded parallel test cases were executed to test various components of IFENs, as well as various simulations representative of a wide variety of problems. IFENs can easily handle meshes with tens of millions of elements. The largest mesh used for the purpose of this research contained 14.58 million tetrahedrons and 2.489 million nodes, which on average required just 7 minutes per timestep, while executing the classical split framework of the CBS scheme. Results from the simulation of 9 carotid meshes, representative of 4 carotid geometries were presented and found to be in good agreement with the available ultrasound data. The flow fields were analysed and post processed using different techniques for each case. The haemodynamic wall parameters like time averaged wall shear stress and oscillatory shear index were calculated and mapped onto the corresponding boundary nodes. The region in the carotid bifurcation susceptible to the deposition of plaques and consequent stenosis were pointed out and other anomalies were highlighted. E-Thesis Computer science.;Biomedical engineering. 31 12 2014 2014-12-31 COLLEGE NANME COLLEGE CODE Swansea University Doctoral Ph.D 2020-09-03T08:59:18.7915549 2018-08-02T16:24:28.4641899 Faculty of Science and Engineering School of Engineering and Applied Sciences - Uncategorised AMARPAL KAPOOR 1 0042216-02082018162437.pdf 10797918.pdf 2018-08-02T16:24:37.2170000 Output 18181210 application/pdf E-Thesis true 2018-08-02T00:00:00.0000000 false
title Development of a parallel CFD solver with application to arterial flows.
spellingShingle Development of a parallel CFD solver with application to arterial flows.
AMARPAL KAPOOR
title_short Development of a parallel CFD solver with application to arterial flows.
title_full Development of a parallel CFD solver with application to arterial flows.
title_fullStr Development of a parallel CFD solver with application to arterial flows.
title_full_unstemmed Development of a parallel CFD solver with application to arterial flows.
title_sort Development of a parallel CFD solver with application to arterial flows.
author_id_str_mv aa2592aa0bb04bdab92450f24643bf5a
author_id_fullname_str_mv aa2592aa0bb04bdab92450f24643bf5a_***_AMARPAL KAPOOR
author AMARPAL KAPOOR
author2 AMARPAL KAPOOR
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 In this research, the finite element method (FEM) was used to solve the nonlinear, incompressible, transient, three dimensional Navier-Stokes equations in their non-conservative form. Linear tetrahedron elements were employed with the elegant, equal order interpolation for both pressure and velocity. The characteristic based split scheme was formulated in a fully implicit manner to circumvent the time step restrictions of the classical explicit formulations. The monolithic (single step, fully coupled solution procedure for pressures and velocity) form of the CBS scheme was also derived and its suitability was positively demonstrated. Casting the CBS scheme in a monolithic framework, results in the generation of a pressure stabilization term in the mass conservation equation, thereby circumventing the LBB restriction by the elimination of the zero pressure block. An account of all the steps involved in discretizing the Navier-Stokes equations (both in split and monolithic frameworks) was presented in meticulous detail, which included the derivation of the convective and pressure stabilization terms, linearization of the non-linear terms and the consequent derivation of the highly efficient analytical jacobian matrix, along with the temporal and spatial discretizations of the corresponding terms. The monolithic and the split version of the CBS scheme were integrated into a parallel, scalable and extensible Fortran90 software called IFENs. The development of IFENs started during the course of this research and all of its components have been designed and implemented by the author of this thesis. Multi processor parallelism was achieved using the Intel implementation of the most widely used/preferred, Message Passing Interface (MPI) standard. The parallel support needed for the use of a variety of parallel, linear, iterative solvers belonging to the Krylov subspace family (e.g. GMRES and its variants, CG, BiCG, BiCG- stab, etc.), parallel non linear solvers belonging to the Newton-Krylov family (line search newton, trust region newton, nonlinear GMRES, etc.) and parallel preconditioners (incomplete LU, Additive Shwarz Method - ASM, algebraic multigrid, etc.), was provided by the incorporation of PETSc into IFENs. PETSc is a state of the art, non-trivial toolkit, which represents a collection of several parallel libraries useful in high performance scientific computing. Keeping in mind the specific requirements of IFENs, a custom mesh partitioner was implemented. It operated on meshes that were renumbered using bandwidth reducing algorithms like Revere Cuthill Mckee. The possibility of using established domain decomposition libraries like ParMETIS was explored and demonstrated to be counter productive for the demands of this research. After the preliminary testing and validation of the procedures adopted before and during the execution of IFENs, large, high definition domains representative of human arteries (specifically, carotid bifurcations, found in the neck) were considered and the complete incompressible set of Navier-Stokes equations were solved for pressure and velocity fields. During the tenure of this research more than 1000 recorded parallel test cases were executed to test various components of IFENs, as well as various simulations representative of a wide variety of problems. IFENs can easily handle meshes with tens of millions of elements. The largest mesh used for the purpose of this research contained 14.58 million tetrahedrons and 2.489 million nodes, which on average required just 7 minutes per timestep, while executing the classical split framework of the CBS scheme. Results from the simulation of 9 carotid meshes, representative of 4 carotid geometries were presented and found to be in good agreement with the available ultrasound data. The flow fields were analysed and post processed using different techniques for each case. The haemodynamic wall parameters like time averaged wall shear stress and oscillatory shear index were calculated and mapped onto the corresponding boundary nodes. The region in the carotid bifurcation susceptible to the deposition of plaques and consequent stenosis were pointed out and other anomalies were highlighted.
published_date 2014-12-31T03:52:32Z
_version_ 1763752596963590144
score 11.012678

Similar Items