Questions tagged [mesh]

The tag has no usage guidance.

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

Meshing surface of a sphere with a subdomain

I am trying to build a triangle mesh of the surface of a sphere which also includes a subdomain defined by a 'polygon'. Here is a successful example (subdomain defined by the red dots): Note that the ...
0
votes
0answers
23 views

Show deformed mesh in Comsol

I am using COMSOL Multiphysics 5.2a in order to do thermal, mechanical and thermo-mechanical simulations. Due to a boundary load, my components get deformed which I can simulate very well. However, ...
1
vote
1answer
99 views

How to compute gradient of a cell having a boundary face?

In many situations in unstructured mesh solvers, one needs to compute gradient of arbitrary variable $\phi$ such as temperature or velocity at face centers (one of such situations is correction for ...
2
votes
1answer
57 views

Going From Blender Structure defined by triangles to full 3D mesh (Using GMSH?)

I currently have created a model airplane in Blender by drawing a closed volume with triangular planes. I want to do a FEM calculation on this object, meaning I need a fine 3D tetrahedral mesh of this ...
2
votes
1answer
124 views

How to calculate skewness for a mesh?

I am writing a code to calculate mesh quality stats such as: cell volume, face areas and non-orthogonality between faces (basically something like OpenFOAM's ...
2
votes
2answers
136 views

What are the best ways to interpolate a vector field inside (convex) polygons?

I want to interpolate a vector field inside convex polygons in a polygonal mesh. For triangular meshes the scheme uses a piecewise constant interpolation in the triangle, discretized at the center of ...
0
votes
0answers
61 views

FEM Meshing artifact at nodes with fewer neighbors

I wrote a 2D-FEM solver to solve some diffusion process and wanted to verify my code with a test problem. The input was $f(x,y) = x^2+y^2$ and I applied the stiffness matrix on it to get $\Delta f = 4$...
3
votes
1answer
107 views

FEM with elastic inhomogeneous properties leads to mesh-induced anisotropy

I'm solving an elastic homogenization problem and I'm having problems with mesh artifacts. I would like to first give a brief summary of what I do: I have a system with inhomogeneous (but isotropic) ...
2
votes
2answers
102 views

Visualization of Quad/Octree data Structures

I have a simple C++ quadtree implementation. The implementation subdivides node into four children as if there is more than one point in that particular node. At the moment I want to visualize how ...
0
votes
0answers
17 views

Volumetric Mesh Cutting Algorithms For Simulating Sculpture

I am interested in whether there an available algorithms to give me insights on how to cut a volumetric mesh?
0
votes
0answers
47 views

Adding a boundary surface on a STL model in Gmsh

Background info I want to run brain simulations with models from the PHM repository (not relevant to the problem but suits the narrative) which are shipped in STL format. The meshes of the STL have ...
1
vote
0answers
23 views

Finding isocontour lines, populating points on said lines (in MATLAB)

I'm trying to extract data from a potential field. I have an example code of a potential field, using meshgrid in MATLAB. I've pasted the code below. What I'd like to do next is 'convert' this field ...
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 ...
0
votes
0answers
22 views

RANS equations in the ALE formulation

I am simulating an incompressible flow of a newtonian fluid over an oscillating plate using OpenFOAM . As far As I know, in this case, OpenFOAM uses the Arbitrary Lagrangian Eulerian formulation. ...
0
votes
0answers
40 views

Import stl file into Gmsh and create a 3d mesh for CFD

I am trying to develop a case study for teaching CFD. To this end, I am using OpenVSP, Gmsh and SU2. I have created a wing in OpenVSP, exported the CFD mesh as .stl file. Opened the .stl file in Gmsh ...
0
votes
1answer
59 views

Parallel mesh partitioning

When a mesh partitioning takes place and every process works on a part of the mesh is any way to rename the global numbering of nodes(on each process) into a local numbering?Is there any software that ...
1
vote
0answers
68 views

Compute mass matrix in vibrations problem by using finite element method

I have to compute the mass matrix of a Hexahedral mesh. There are 3 methods to compute mass matrix. I'm interested in one method which consists of dividing the mass of the element by the number of ...
2
votes
2answers
324 views

How to determine the Jacobian Ratio for triangle element?

I am trying to implement an algorithm to find the Jacobian ratio for each triangle in mesh as a part of mesh quality check. Let's say that I have vertices of the triangle: $P_1(x_1, y_1, z_1)$, $P_2(...
4
votes
2answers
133 views

Effect of mesh size on solution curves for a 1D problem

I'm interested in studying the effect of mesh size on the behavior of the solution curves of 1D convection-diffusion problem. $$\frac{\partial C}{\partial t} = D\frac{\partial ^2 C}{\partial x^2} - ...
2
votes
1answer
164 views

mesh dependence of numerical adjoint solution

I am solving the steady, two-dimensional adjoint Euler equations, $$A_x^T \partial_x \Psi + A_y^T \partial_y \Psi = 0$$, where $A_x = \partial F_x/\partial U$ and $A_y= \partial F_y/\partial U$ are ...
0
votes
1answer
37 views

Produce vertex displacements from volumetric shrinkage data on unstructured meshes

I was wondering what would be an efficient way to produce compatible displacements for mesh nodes/vertices if the computed data is volume shrinkage of each element/cell in the unstructured mesh? ...
0
votes
0answers
28 views

When semicoarsening is needed in multigrid method?

I am trying to understand multigrid in deep and all the coarsening techiques. So,assuming a 2D grid when semicoarsening is prefered instead of standard coarsening?What is the error's behaviour when a ...
0
votes
0answers
77 views

FEniCS, refinement not 'respecting' domain boundary

Short question: how to ensure that extra points are not included as 'boundary' points after calling the refine function. More details. I am working with a hexahedral mesh in $3$d. Let $X$ be the set ...
0
votes
0answers
23 views

Direct and Inverse efficient mapping of 3D cartesian positions in a 1D array

I saved a sample of $N$ Cartesian locations $\{x_i, y_j, z_k\}$ inside a one-dimensional array $\mathbf{a} = \{a_l\}_{l = 1}^N$. How can I access back (efficiently) the $l$-location of the array $\...
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 ...
0
votes
0answers
47 views

Determining the indices of a VTK mesh

I am trying to assign fixed constraint at specific indices of VTK mesh, however, I only can view STL file in a blender as follows: Upon assigning the favored indices in my scene at the SOFA physics ...
3
votes
2answers
109 views

Find shortest path around a cylinder represented by 3d triangular mesh

Suppose I have a 3d triangular mesh with the topology of a finite cylinder. Let $C$ be a vertex on that mesh. How can I find the shortest path from $C$ to itself that goes around the cylinder? By ...
3
votes
0answers
97 views

Is it possible to retain sets (physical tags) after remeshing the model in Gmsh?

I am currently working on a remeshing algorithm for Abaqus using Gmsh as part of my Bachelors thesis and I seem to be stuck at some point. I am able to remesh .geo and .msh files using the script ...
2
votes
1answer
72 views

Sorting Triangle vertices based on indices for alpha blending (back to front) C++, OpenGL

I have implemented in OpenGL alpha blending for a triangle mesh. I render this mesh using a vector of indices that indicate which vertex from a vertices vector should be drawn. The problem is to ...
0
votes
0answers
81 views

How to read data from simply formatted text file (c++)?

Mesh information like points, faces and cells is to be stored into separate files: e.g. for points file: # points data: x y z N_Points 100 x1 y1 z1 x2 y2 z2 ... cell file: # cells data: ...
0
votes
0answers
52 views

Stokes Equation fails to converge for an ellipse

This might be because of the mesh, but the following code blows up for all values of b not 1. Does anybody have any experience working with the ellipse mesh in Fenics? ...
4
votes
2answers
3k views

Creating 3D Mesh from stl files with gmsh

After long hours of searching for an answer I thought it might be better to ask the community. The problem I have is that I need to convert STL files to mesh files. I know that I therefore need to ...
3
votes
1answer
72 views

Smoothness regularisation of a 2D field on a triangular mesh?

I'm working on an inverse problem where the solution is the values of a 2D scalar field at the vertices of a 2D triangular mesh, such that the field can be defined continuously inside the mesh via ...
0
votes
1answer
174 views

Node re-numbering in 1D mesh GMSH

I'm working with Gmsh to generate 1D meshes of polygon edges. I have noticed that drawing a polygon by "vertices and lines" the program always assigns to the first N nodes of the mesh the coordinates ...
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
142 views

Node renumbering in a 2D mesh

I have a 2D domain which is discretized using Q4 elements. I have the nodal positions and the element connectivity matrix. I would now like to renumber the nodes in such a way that all the interior ...
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 ...
4
votes
2answers
542 views

3D contour mesh computation

I have the value of a function in three dimensions, $f(x, y, z)$, which varies smoothly. I would like to compute a 3D mesh where $f(x, y, z)$ has a particular value. Are there algorithms to do this? ...
4
votes
1answer
175 views

Convergence rate Jacobi/Gauss-Seidel with mesh resolution

In the book A Multigrid Tutorial - Briggs, Henson. McCormick in the beginning of Chapter 3, it is mentioned that ...because the convergence factor behaves as 1-$O(h^{2})$, the coarse grid ...
0
votes
1answer
79 views

Structural mechanics simulation using FLUENT compared to analytical solution

This is a continuation of my previous post, 1D analytical solution vs FEM solution for a bar under compression. For some reason, I cannot comment in it. The analytical solution to the 1-D static ...
0
votes
1answer
108 views

Gmsh exporting wrong mesh DATA [closed]

so hopefully I'll be using gmsh to make meshes out of 2-D cross sections with o thickness. I tried to make a structured mesh with quad elements of a rectangular ...
4
votes
1answer
379 views

Expanding Winding Number algorithm to arcs

I have a problem that I have been attempting to solve for a few days now. I was wondering if I would be able to get some assistance from the community. In order to detect if a point is in a polygon, ...
3
votes
1answer
253 views

Compute mesh of the projection of a 3D surface triangulation

Given a triangulated surface in $\Bbb{R}^3$ we can simply project it on a plane. This will result in a family of triangles which do not form a mesh of the projection for the following reasons: each ...
1
vote
0answers
103 views

Mesh error quantification using error norm for a sharp logistic function [closed]

Consider the following logistic function in 1D for the domain $0$ to $1$. $$ f(x) = \frac{1}{1 + \exp(-50(x-0.5))} $$ Let's say I resolve this on two meshes: first with 11 equispaced points (0.0, 0....
3
votes
1answer
137 views

Discontinuous Galerkin FEM : Control points are mid-points of edges instead of nodes

I am thinking to use discontinuous galerkin FEM (DGFEM) method to estimate discontinuous displacement field $u: \Omega \rightarrow \mathbb{R}^2$ at the crack surface of a material. The domain is ...
1
vote
1answer
497 views

Convert Unstructured Mesh to Structured

I have solved a simple problem on Openfoam, which is a 3D rectangular parallelopiped, filled with rectangular hexahedrons as mesh elements. The structured, mapped mesh was made in ICEM and solved in ...
2
votes
1answer
236 views

Barycentric interpolation equivalent for irregular hexahedra

I have a mesh with irregular hexaedra and I need a fast way to interpolate values at points inside these cells. I know that trilinear interpolation does not work well for large skews. Barycentric ...
1
vote
0answers
147 views

How to choose seed in flood/fill algorithm

I was able to draw to the outer airtight boundary of hole. Now I would like to implement interior flood/fill algorithm. To get an internal seed, only one method comes to my mind which is testing ...
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
2answers
173 views

Sign of integer determinant 4 by 4

I'm in the context of this publication: http://www.gilbertbernstein.com/resources/booleans2009.pdf I applied quantization to my point coordinates: All coordinates are integer lying in [0, 2 power 20]....