Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method

Detalhes bibliográficos
Autor(a) principal: Santos, Milene Francisco dos
Data de Publicação: 2022
Tipo de documento: Dissertação
Idioma: eng
Título da fonte: Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)
Texto Completo: http://hdl.handle.net/10316/103048
Resumo: Dissertação de Mestrado em Matemática apresentada à Faculdade de Ciências e Tecnologia
id RCAP_bac1fcfd27e0b0f1bb26571d3d4caf5b
oai_identifier_str oai:estudogeral.uc.pt:10316/103048
network_acronym_str RCAP
network_name_str Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)
repository_id_str 7160
spelling Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin MethodOptimização de ordem em domínios curvos: método de Galerkin descontínuoquação de Helmholtzmétodo de Galerkin descontínuoomínios com fronteira curvarescontrução polinomialHelmholtz's equationdiscontinuous Galerkin methodcurved boundary domainspolynomial reconstructionDissertação de Mestrado em Matemática apresentada à Faculdade de Ciências e TecnologiaO problema que motivou a escolha do tema desta tese prende-se com descrição do comportamento das ondas electromagnéticas na córnea humana, a fim de compreender as razões que conduzem à opacidade da córnea. As equações de Maxwell, que descrevem a propagação de campos electromagnéticos, constituem o modelo matemático natural para o nosso estudo. Neste trabalho, optámos por considerar um modelo simplificado dado pela equação de Helmholtz, assumindo a variação harmónica no tempo dos campos electromagnéticos. Para resolver a equação de Helmholtz recorremos ao método dos elementos finitos descontínuos Galerkin (DG). Uma vez que estamos interessados em resolver a equação num domínio com fronteira curva, que pretende mimetizar a córnea humana, as malhas poligonais nãomalhas poligonais não se ajustam exactamente ao domínio físico, o que conduz a uma redução da precisão do método numérico.Nesta tese propomos duas abordagens para lidar com a redução da precisão do método da DG num domínio com fronteira curva, na resolução do problema de Helmholtz em duas dimensões com condições de fronteira de Dirichlet homogéneas. O primeiro método designa-se por DG-ROD (Reconstruction for Off-site Data) e baseia-se na reconstrução polinomial da condição de fronteira imposta no domínio computacional, tendo em conta a condição de fronteira imposta no domínio físico. Os testes numéricos realizados evidenciam a diminuição do erro e o aumento da ordem de convergência do método, em relação ao método DG clássico. O segundo método proposto designa-se DG-NM (Nelder-Mead) com controlo do passo e baseia-se na alteração da condição de fronteira imposta no domínio computacional através da resolução de um problema de minimização sem restrições. Esse problema de minimização é resolvido com uma variante do método de Nelder-Mead. Os testes numéricos realizados evidenciam a diminuição do erro em relação ao método DG clássico e uma diminuição do número de iterações em relação ao método de NM clássico. Ambos os métodos sugeridos nesta tese têm a vantagem de não utilizar malhas com elementos curvos nem transformações não lineares do elemento da malha para o elemento de referência, comparativamente a outras alternativas para lidar com domínios com fronteira curva.The problem that motivated the choice of the subject of this thesis is related to the description of the behaviour of electromagnetic waves in the human cornea, in order to understand the reasons that lead to the opacity of the cornea. Maxwell's equations, which describe the propagation of electromagnetic fields, are the natural mathematical model for our study. In this work, we have chosen to consider a simplified model given by the Helmholtz equation, assuming the harmonic variation in time of the electromagnetic fields. To solve the Helmholtz equation we use the discontinuous finite element method Galerkin (DG). Since we are interested in solving the equation in a domain with a curved boundary, which intends to mimic the human cornea, polygonal meshes do not polygonal meshes do not fit exactly into the physical domain, which leads to a reduction in the accuracy of the numerical method.In this thesis we propose two approaches to deal with the reduction in the accuracy of the DG method in a domain with a curved boundary, in solving the Helmholtz problem in two dimensions with homogeneous Dirichlet boundary conditions. The first method is called DG-ROD (Reconstruction for Off-site Data), is based on a polynomial reconstruction of the boundary condition imposed on the computational domain that takes into account the boundary condition imposed on the physical domain. The numerical tests show a reduction of the error and an increase in the order of convergence of the method, in relation to the classical DG method. The second method proposed, called DG-NM (Nelder-Mead) with step size control, is based on changing the boundary condition imposed on the computational domain by solving an unconstrained minimisation problem. This minimisation problem is solved with a variant of the Nelder-Mead method. The numerical tests evidence a decrease in the error in relation to the classical DG method and a decrease in the number of iterations in relation to the classical NM method. Both methods suggested in this thesis have the advantage of not requiring the generation of curved meshes to adjust the boundary nor complex nonlinear transformations to map the curved elements to the reference one, in relation to other alternatives to deal with curved boundary domains.Outro - FEDER -- Fundo Europeu de Desenvolvimento Regional, through COMPETE 2020 -- Programa Operacional Fatores de Competitividade, and the National Funds through FCT -- Fundação para a Ciência e a Tecnologia, project no. POCI-01-0145-FEDER-028118. PTDC/MAT-APL/28118/20172022-07-22info:eu-repo/semantics/publishedVersioninfo:eu-repo/semantics/masterThesishttp://hdl.handle.net/10316/103048http://hdl.handle.net/10316/103048TID:203080114engSantos, Milene Francisco dosinfo:eu-repo/semantics/openAccessreponame:Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)instname:Agência para a Sociedade do Conhecimento (UMIC) - FCT - Sociedade da Informaçãoinstacron:RCAAP2022-10-19T20:38:16Zoai:estudogeral.uc.pt:10316/103048Portal AgregadorONGhttps://www.rcaap.pt/oai/openaireopendoar:71602024-03-19T21:19:55.180335Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos) - Agência para a Sociedade do Conhecimento (UMIC) - FCT - Sociedade da Informaçãofalse
dc.title.none.fl_str_mv Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
Optimização de ordem em domínios curvos: método de Galerkin descontínuo
title Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
spellingShingle Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
Santos, Milene Francisco dos
quação de Helmholtz
método de Galerkin descontínuo
omínios com fronteira curva
rescontrução polinomial
Helmholtz's equation
discontinuous Galerkin method
curved boundary domains
polynomial reconstruction
title_short Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
title_full Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
title_fullStr Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
title_full_unstemmed Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
title_sort Order Optimisation in Curved Boundary Domains: Discontinuous Galerkin Method
author Santos, Milene Francisco dos
author_facet Santos, Milene Francisco dos
author_role author
dc.contributor.author.fl_str_mv Santos, Milene Francisco dos
dc.subject.por.fl_str_mv quação de Helmholtz
método de Galerkin descontínuo
omínios com fronteira curva
rescontrução polinomial
Helmholtz's equation
discontinuous Galerkin method
curved boundary domains
polynomial reconstruction
topic quação de Helmholtz
método de Galerkin descontínuo
omínios com fronteira curva
rescontrução polinomial
Helmholtz's equation
discontinuous Galerkin method
curved boundary domains
polynomial reconstruction
description Dissertação de Mestrado em Matemática apresentada à Faculdade de Ciências e Tecnologia
publishDate 2022
dc.date.none.fl_str_mv 2022-07-22
dc.type.status.fl_str_mv info:eu-repo/semantics/publishedVersion
dc.type.driver.fl_str_mv info:eu-repo/semantics/masterThesis
format masterThesis
status_str publishedVersion
dc.identifier.uri.fl_str_mv http://hdl.handle.net/10316/103048
http://hdl.handle.net/10316/103048
TID:203080114
url http://hdl.handle.net/10316/103048
identifier_str_mv TID:203080114
dc.language.iso.fl_str_mv eng
language eng
dc.rights.driver.fl_str_mv info:eu-repo/semantics/openAccess
eu_rights_str_mv openAccess
dc.source.none.fl_str_mv reponame:Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)
instname:Agência para a Sociedade do Conhecimento (UMIC) - FCT - Sociedade da Informação
instacron:RCAAP
instname_str Agência para a Sociedade do Conhecimento (UMIC) - FCT - Sociedade da Informação
instacron_str RCAAP
institution RCAAP
reponame_str Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)
collection Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos)
repository.name.fl_str_mv Repositório Científico de Acesso Aberto de Portugal (Repositórios Cientìficos) - Agência para a Sociedade do Conhecimento (UMIC) - FCT - Sociedade da Informação
repository.mail.fl_str_mv
_version_ 1799134092856393728