Questions tagged [mesh-generation]

The generation of a set of interconnected geometric objects that approximate a desired domain or region in euclidean space.

Filter by
Sorted by
Tagged with
3
votes
1answer
126 views

Mesh with constraints

Is it possible to construct a constrained tetrahedral mesh of a domain using Tetgen or similar software? What I mean by constrained is that there are some nodes or edges that are not free but are ...
3
votes
0answers
534 views

Gmsh for 3D volume with inclusions [closed]

In an attempt to create three-dimensional volumes with inclusions in Gmsh I stumble upon a problem which was non-existent in the two-dimensional case. I'm using the OpenCASCADE geometry kernel ...
3
votes
0answers
222 views

matplotlib contourplot for $\log z$ in the Complex Plane $\mathbb{C}$

I tried using Python's matplotlib on the logarithm and here is what I got, a kind of starburst pattern. Since the angle jumps between $\theta = 0$ and $\theta = 2\pi$, contour assumes there is a ...
2
votes
3answers
138 views

Flux sign and face normal confusion in finite volume method

I implemented a solver for the 2D steady-state heat equation (without heat generation and homogeneous material) $\nabla. (k\nabla T) = 0$, using finite volume method, however, I am having some ...
2
votes
3answers
166 views

For traditional FEM and FVM, why can't we use mesh to represent geometry and use the mesh which represent the geometry to do the computation directly?

Isogeometric analysis [1] has the advantage of integrating geometric and mesh models using NURBS or Spline. At the same time, I would like to ask a question to my friends: for traditional FEM and FVM, ...
2
votes
2answers
85 views

How to find out the difference between a structured and unstructured mesh using the file containing the mesh information?

I have two different mesh files (both are .inp files obtained from Abaqus) that represent the exact same geometries with the same boundary conditions, etc. The only difference is that one of them is ...
2
votes
1answer
88 views

Gmsh meshes flat faces incorrectly for cylindrical faces

I have some C++ code that generates meshes from step files and then analyses these meshes for visibility of the faces from different viewing directions. I currently use CGAL but I would like to switch ...
2
votes
1answer
3k views

Creating a proper quad-mesh in GMSH for an “I”-shaped geometry

I'm trying to make a 2-D quad-only mesh of an "I"-shaped cross-section, with elements of size 5. So far I've managed to make it quad only, but some elements are very distorted when they shouldn't be. ...
2
votes
2answers
616 views

Combined translational and rotational meshing in gmsh

For a flow simulation I am trying to reproduce a specific screw element design for an intermeshing co-rotating twin-screw extruder. I am using gmsh (v2.16) where the element is made from a 2D surface ...
2
votes
2answers
631 views

(Fortran) Integrating/summing over complicated 3D domain

I have some function $F(k_x,k_y,k_z)$ that I wish to numerically integrate over a polygon domain - physically, I am integrating over the first Brillouin Zone (BZ) of the FCC lattice (a truncated ...
2
votes
2answers
1k views

How can I generate shell elements to a mesh

I have a program that generates mesh for given 3d models. The generated mesh must use quadrangle elements. It is required to add "shell elements generation capabilities". The requirement is as ...
2
votes
2answers
1k views

Mesh domain decompositions / mesh partitioning

I have some experience with mpmetis from METIS. It is pretty good software which offers unstructured mesh grid partitioning. But obtained results always minimize edgecuts or total communication volume....
2
votes
1answer
380 views

Filling out a mesh (for Fenics FEM)

I have been using Fenics to model some simple electrostatic problems, but now I need to simulate more complicated geometries. Assuming I have created a mesh of several objects, and know their ...
2
votes
1answer
65 views

Mesh partition quality

I am working on static FE mesh partitioning and in order to achieve a good quality partitioning I want to know how to drecrease interprocessor communication by increasing the connectivity of elements ...
2
votes
1answer
541 views

Two-dimensional mesh in fem: generating P1, P2, P3,… mesh from a P1 mesh

I have a two-dimensional mesh generated by triangle (the mesh generator software is not relevant). This software generates a perfect mesh for approximate the solution by piecewice linear functions (P1 ...
2
votes
1answer
209 views

Mesh generation with mostly the same elements (type/size..)

I'm searching for a way to generate a mesh that consists mostly of the same elements and i mean not only the same type, but the exact same size to ease the FEM-calculation. Is there a particular ...
2
votes
2answers
733 views

CFD mesh generation

Is there any libraries(python is preffered) for generation rectangular 3d mesh (for using with the finite difference method), which will be consistent with boundaries? (I've asked this question at ...
2
votes
2answers
250 views

Algorithm to compute the intersection of meshlines with a boundary

I need a program or an algorithm that computes the intersection of a mesh and a boundary. The mesh is structured orthogonal in nature and the boundary is a circle (for example). This will be used ...
2
votes
0answers
42 views

Isosurface extraction when you know the derivatives

I have a complicated equation that defines a shape in 3D, and I would like to generate a surface mesh. The shape is defined by an isosurface, i.e. the function is positive inside the shape and ...
2
votes
0answers
124 views

scalable parallel mesh/amr on unstructured grid

I am trying to code a scalable parallel AMR for unstructured grid. There seems to be three approaches for this a) Store some global grid info on each processor and partition with parmetis (The ...
2
votes
0answers
342 views

Adaptive Finite Element Method - Laplace

I'm currently attempting to turn my code for solving the laplace equation using finite element approximations into an adaptive one using the dual weighted residual as my error estimator: i.e. my ...
1
vote
1answer
943 views

Algorithm to find boundary faces of mesh

I am given the data about nodes and elements of a tetrahedral mesh. What algorithm do I use to find the boundary faces of the mesh?
1
vote
1answer
353 views

Library for generating Discontinuous Galerkin FEM mesh

There are a number of packages available for generating Continuous Galerkin (CG) FEM meshes (gmsh, tetgen, netgen, etc.), but I have been unable to find one that generates Discontinuous Galerkin (DG) ...
1
vote
2answers
390 views

Which type of meshing is more suited for simulation of electromagnetic metamaterial unit cells?

Electromagnetic metamaterials are resonant, periodic structures (repetition of a unit cell) involving both dielectric and conducting elements, and are used and simulated in frequency ranges from ...
1
vote
3answers
753 views

Create mesh for complicated 3D object for finite element analysis

I see images of steel connections, concrete dams, and other complicated 3D objects in papers which finite element analysis has been performed on them. My questions are: How these objects are created ...
1
vote
1answer
1k views

How can I coarsen a mesh in Gmsh when 'Mesh options' include 'Refine by splitting' but nothing about coarsening?

I am new to Gmsh and I am having trouble creating a circle with a coarse mesh. I use Geometry->Elementary entities->Add->Circle to create the geometry ...
1
vote
1answer
2k views

Determinant of jacobian matrix

I am using an FEM code written in Fortran which I did not design. For a particular problem, the program complains that the determinant of the Jacobian matrix is inferior to zero. I vaguely ...
1
vote
1answer
2k views

Triangular mesh of a 3D surface

I have an irregular grid of points describing this surface (a large subduction fault in South America). The color is depth. Anyway I have 3D coordinates (lon,lat,depth) at irregular intervals. I'm ...
1
vote
2answers
90 views

What are the alternatives for ABAQUS in generating an *.inp file from a CAD model

ABAQUS gives a .inp file (in pre-processing stage) where the information with regard to the preprocessed model is defined, information such as geometry, mesh type, number of elements, boundary ...
1
vote
1answer
202 views

Is there any fundamental difference between meshing for FEM, FVM and FDM?

I am a novice to the field of computational science and have just started studying the FDM and FEM (haven't started on FVM yet). While trying the understand the subject I got this question and trying ...
1
vote
1answer
48 views

stretch elliptic mesh to fit a circle

I've generated a 2D O-mesh around an airfoil. Unfortunately, the mesh O shape is not a perfect circle as you can see in the figure below (in the figure I plotted only the mesh nodes in blue). The red ...
1
vote
1answer
421 views

Global numbering of the vertices in a triangulation consisting of quadratic Lagrange elements

I've asked this question on Computational Science SE too. Consider the following triangulation of a rectangle $(0,a)\times(0,b)$: We have two types of triangles: I enumerate the rectangles(...
1
vote
2answers
517 views

structured to unstructured grid - calculating volume of cell

I have a 3-D structured grid that I am trying to convert over to an unstructured grid. Here is a picture of the 2-D structured grid with the nodes shown. It is axisymmetric about the x axis (y = 0). ...
1
vote
1answer
96 views

Partial differential equation FEM application

I have a PDE which looks like Helmholtz wave equation on one dimensional domain. $$\dfrac{d^2u(x)}{dx^2}+\pi^2u(x)=f(x)$$ where $-\infty <x<\infty $ Also, $f(x)= 1$ for $-0.25<x<0.25$, I ...
1
vote
1answer
50 views

Filling a volume with generalised polyhedra

Given a volume (say, some polyhedron), I need to fill it with smaller polyhedra, such that the space is filled as much as possible. The constraints and relaxations are: (0) For a computation ...
1
vote
1answer
267 views

How to generate a face list from vertices?

I have a little background in writing toy finite volume CFD codes. In 2D Cartesian scenarios, I typically take $x_{\min}$, $x_{\max}$, $y_{\min}$, $y_{\max}$, and the number of points in $x$ and $y$ ...
1
vote
1answer
773 views

A difference running Gmsh in command line and in GUI for 3-D mesh generation [closed]

In Gmsh, I have a test.geo file. The question is: if I run gmsh -3 test.geo in the command line, I get the below error. ...
1
vote
1answer
440 views

Can TETGEN generate triangulation of a 2D point set?

I would like to link my Fortran 90 application to a library which would give me Delaunay triangulation (for a 2D set of points) and tetrahedralization (for a 3D set of points). As suggested in Fastest ...
1
vote
1answer
870 views

Tetrahedral Mesh file format in TETGEN

Ply file formats is used to represent mesh structures in the TETGEN freeware package. All the examples I've seen of 'ply' meshes were used for triangular meshes. ...
1
vote
0answers
45 views

Determining the voxels between two boundary surfaces

Update (Solution) -17 October 2020- I finally did it! You can check the code on my GitLab repo. The function is to do this is the domain_extract which identifies ...
1
vote
0answers
132 views

Duplicate Surface Meshes in different Gmsh files

My objective is to get a conformal surface mesh between two geometries meshed independently in different files. By example, I can start with a simple cube in one .geo file and mesh this geometry. ...
1
vote
0answers
75 views

Metis: how to use and tutorial recommendation

I am new to METIS and trying to use it in my fortran code. I read the manual online. But still, I am not sure about how to implement it my code. I tried the test cases in the graphs directory. For ...
1
vote
0answers
57 views

Efficient initial identification of solid or liquid domains for a block structured Cartesian grid generation system

INTRO Within the last 5 days I was able to generate a block structured Cartesian grid generation system with a combination of Fortran,C++ and Python. I am running intersection tests of the ...
1
vote
0answers
64 views

surface triangulation to boundary representation

Are there any open source tools available which can take a surface triangulation of mesh volumes and convert them to a boundary representation? I am interested in refining multiple region tetrahedral ...
1
vote
0answers
62 views

adaptive / smart grid choice for dynamic programming

I've a read a paper about investing where they used a dynamic programming approach to solve a finite horizon problem, i.e. $$\max_{x_t} E[u(W_T)] $$ where $u$ is a utility function and $W_T$ ...
1
vote
0answers
184 views

How to generate a tet10 mesh from a tet4 mesh

I am writing a FEN solver on tet10 elements. I right now have a CFD grid gen tool (Pointwise) which generates tet4 meshes. From this I have the element connectivity list and the node list. How do I go ...
1
vote
0answers
103 views

Generate discrete set of points in a feasible region

I have two vectors which specify the bounds $x_{min}$ and $x_{max}$ of the sample space. Also, it has to satisfy the linear constraint $Ax \leq b$. How to generate an evenly spaced set of points, ...
1
vote
0answers
84 views

Mesh partitioner with user-defined overlap

I am looking for a mesh partitioner, where I can specify overlap, for example h = 3. I have looked into metis, but I wasn't able to find such a functionality. Is there any other package which ...
1
vote
0answers
31 views

Size of the support

I am working on moving least square approximation. To well-define of this method, the following matrix of size $n \times m$ must be of rank $m$. For $m=3$ $${\bf A }=\begin{pmatrix} 1&x_1&y_1\\...
1
vote
0answers
126 views

Finding boundary intersection points with Cartesian grid?

I would like to find the intersection points e.g. $M_x$ and $M_y$ as in the attached figure. The boundary (solid line) is defined by the Lagrangian points. I am working in C++ (basic - medium ...