gridap / gridap/Gridap.jl

Higher order triangulation of sphere

Open
#911 2 comments 0 reactions 1 assignee Claimed by @amartinhuertas View on GitHub
Dominant language
Julia
Stars
879
Forks
119
Avg merge
6d 3h
Merged PRs (30d)
4

Description

Hi,

I am trying to solve a Vector-Laplace equation on the sphere using Gridap. I wanted to test the influence of the geometry approximation, reason why I would like to use higher order geometries. I tried to look for possible ways to implement it but still haven't found a feasible way to do it. In particular:

- Generating a surface mesh with gmsh and trying to import the .msh file with the function "GmshDiscreteModel" from GridapGmsh throws the error `ERROR: Only one element type per dimension allowed for the moment. Dimension 3 has 0 different element types`. I attach a .txt file containing the text of the .geo file I am using [sphere.txt](https://github.com/gridap/Gridap.jl/files/11887821/sphere.txt)
- Generating a volume mesh with gmsh, importing it using "GmshDiscreteModel" and then use only the boundary works, but it seems like it is not possible to import higher order .msh files and the error `ERROR: For the moment only for first-order elements` is thrown (Also, in this case how can I access normals?)
- Trying to use the ready made functions in "GridapGeosciences", precisely "https://github.com/gridapapps/GridapGeosciences.jl/blob/master/src/CubedSphereDiscreteModels.jl" and "https://github.com/gridapapps/GridapGeosciences.jl/blob/master/src/CubedSphereTriangulations.jl" seems to successfully complete the Triangulation, but when using "get_normal_vector" it throws `ERROR: MethodError: get_normal_vector(::Gridap.Geometry.BodyFittedTriangulation{2, 3, Gridap.Geometry.UnstructuredDiscreteModel{2, 3, Float64, Gridap.Geometry.Oriented}, Gridap.Geometry.UnstructuredGrid{2, 3, Float64, Gridap.Geometry.Oriented, Nothing}, Gridap.Arrays.IdentityVector{Int64}}) is ambiguous`. I attach a minimal example using the above mentioned files: [test.txt](https://github.com/gridap/Gridap.jl/files/11888063/test.txt).

I was wondering if you know possible issues with the code or a simpler way to implement it.

I am using Julia v1.8.5, Gridap v0.17.17, GridapGmsh v0.6 and FillArrays v0.12.8.

Thanks a lot for the help

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.