Questions tagged [adaptive-mesh-refinement]
For questions about AMR, a technique for dynamically updating the precision of a grid to ensure a certain accuracy for a given region.
63
questions
2
votes
1
answer
99
views
AMR framework for efficient simulation of PDEs, potentially with boundary layers
I am aiming at developing a simulation software mostly aimed at reactive fluid flows applications. Here are a simplified list of points to drive the initial choices (code structure, libraries...):
I ...
- 1,764
1
vote
1
answer
88
views
Adaptive mesh refinement with inter-element continuity
I am searching for a library that can perform adaptive mesh refinement (AMR) on large distributed unstructured meshes. For now, the cells are high-order quads/hexa.
I was looking into p4est, which is ...
- 707
0
votes
1
answer
116
views
AMR-Capable meshing software that is not based on quad/octrees
I am looking for AMR/re-meshing software (structured grids would suffice) that is NOT based on quad-octrees, i.e., a fixed refinement rate of 2 but (ideally) something user defined, i.e., ratios of 1....
- 904
5
votes
0
answers
97
views
Grid rearrangement for anisotropic transport simulations
In computational plasma physics, one often faces problems with extremely high transport anisotropy, 1e6 and more, since transport along the magnetic field is much faster than across. To deal with this ...
- 2,370
0
votes
0
answers
36
views
Refluxing step on Finite difference AMR
Hi I am a computer scientist working on MHD code for astrophysics simulation. We use a finite difference scheme where we first solve the spatial derivatives and with them solve the right hand side and ...
- 121
5
votes
0
answers
89
views
How to get an "optimal point" for refinement in FEM adaptive mesh refinement?
Consider the following 1D problem
\begin{align*}
\begin{cases}
\displaystyle
-\frac{d^2u}{dx^2} = f(x), \hspace{0.5cm} x\in (a,b) \\[4mm]
u(a) = u_{a}, \ \ u(b) = u_{b}
\end{cases}
\end{align*}
I ...
- 51
1
vote
1
answer
99
views
numerical integration of integrals in the p-adaptive version of the finite element method
In the p-adaptive version of the finite element method, elements are allowed to have shape functions with arbitrary different polynomial orders. Therefore regarding a 2D problem with quadrilateral ...
- 194
0
votes
1
answer
159
views
Mesh refinement in the Finite Element Method
I need some good references on how to implement programmatically the hp-refinement of meshes in the Finite Element Method in two/three-dimension. I've searched the web a lot and read many articles and ...
- 194
3
votes
1
answer
206
views
Comparison on adaptive mesh refinement on finite elements and finite differences
My current work requires using (Adaptive Mesh Refinement) AMR to resolve multi scale physics. I have a general question whether finite element is better than finite difference in this aspect or not.
I ...
- 119
2
votes
1
answer
138
views
Connectivity of octree as grid
What is the algorithm which transforms the octree/quadtree into a grid? i.e, getting the adjacency or connectivity between the cells in order to be able to perform interpolation, gradient,...etc. I am ...
- 65
3
votes
1
answer
147
views
Proof of R. Verfürth paper on adaptive mesh and bubble functions
I'm studying adaptive meshes, and my professor wrote the following property for a bubble function ( see this scicomp post for the definition I'm using)$b_T$ defined on a triangle $T$.
$$||b_T \phi ||_{...
- 261
4
votes
1
answer
144
views
Computing the residual in a Dual Weighted Residual (DWR) method
I am in the process if computing the Dual-Weighted Residual (DWR) for a linear PDE with a linear functional but I am struggling with the residual part of the calculation.
For example suppose we want ...
- 141
0
votes
1
answer
154
views
Mesh aspect ratio issue with adaptive mesh refinement (AMR)
I am working on implementation of AMR for my finite volume code. Let me use a 2-D mesh to describe my question.
Starting with SINGLE initial cell (let the mesh refine level k = 0) as a root of a quad-...
- 31
1
vote
1
answer
640
views
Finite Difference libray C++
What is the best FD library (or collection of libraries) for C++ codes?
I am looking for some data structure implementation that offers the possibility to do parallel computations on adaptively ...
- 23
1
vote
1
answer
68
views
Recovery of smoothed continuous stresses using the Z^2 error estimator
I have trouble implementing the $Z^2$ error estimator (The superconvergent patch recovery and a posteriori error estimates. Part 1: the recovery technique by Zienkiewicz and Zhu).
For this, I am ...
- 101
2
votes
0
answers
387
views
Implementation of Z^2 error estimator in Abaqus for adaptive mesh refinement
Currently, I am working on a remeshing routine for my simulations (Abaqus 6.14-1) using python scripts. The simulation deals with the Brinell indentation test and as the remeshing software I use Gmsh ...
- 101
3
votes
2
answers
512
views
Data structures of AMR(Adaptive Mesh Refinement) with quadtree
I am graduate student and recently tried to implement quad-tree based AMR. Well I have implemented simple Poisson Solver on node-based quadtree, but I am not sure that my data structure is right. To ...
4
votes
2
answers
1k
views
simple example of an adaptive mesh refinement code
I have been writing an adaptive mesh refinement (amr) code. As a prototype for the code, I have been looking at an adaptive mesh refinement code written by my adviser (written in c). I find looking at ...
4
votes
2
answers
152
views
Good references for dual-weighted residual (DWR) method for goal-oriented adaptive mesh refinement (AMR)
Can anyone help me with good references (books or papers) where I can learn about dual-weighted residual (DWR) method for goal-oriented adaptive mesh refinement (AMR)?
- 63
2
votes
1
answer
66
views
Oscillation term in a posteriori error estimator
Assume that in the a (residual type) posteriori error estimator of some PDE is a term of the form $h_T\|g\|_{L^2(\Omega)}$ involved where $h_T$ is the diameter of an element and $g$ is some known data ...
- 123
4
votes
3
answers
905
views
How can I make sure the flow is divergence-free when I use moving mesh?
I am using projection method and P2/P1 finite element method to solve the incompressible Navier-Stokes equations while the mesh is constantly adapted as the body moves (edge swapping, splitting and ...
3
votes
1
answer
654
views
p-refinement in adaptive methods
Some reference in adaptive techniques say that when the solution $u$ is smooth enough we can use p-refinement instead of h-refinement. And when we have for example singularity, we should use h-...
- 523
2
votes
0
answers
153
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 ...
- 233
1
vote
0
answers
64
views
Non-cubic blocks in Adaptive Mesh Refinement
Most Adaptive Mesh Refinement softwares/logic have cubic blocks. I have explored the BoxLib library and it ONLY supports cubic blocks. To be more precise it supports non-cubic blocks at the coarsest ...
- 467
1
vote
0
answers
273
views
Algorithm for Adaptive Mesh Refinement
I am trying to implement Adaptive Mesh Refinement. I am not a Mathematics/Computational Science person so I will try to write the algorithm in a simpler way. I will be grateful if experts can comment ...
- 467
1
vote
1
answer
269
views
fast adaptive quadrature on equispaced 2-D grid
I need to numerically evaluate 2-D integrals of the form:
$$
\mathcal{I}(\theta) = \int_{0}^{1} \int_0^1 \varphi_\theta(x,y) dx dy
$$
where $\varphi_\theta$ is a family of smooth functions indexed by ...
- 183
1
vote
0
answers
61
views
Square error estimate for adaptive mesh refinement
In a particular implementation(Finite volume advection using upwind) of adaptive mesh refinement the error square estimate for a cell C is given as
$$
\sum_{i = x,y,z} vol * \frac{1}{12} * h^{2} * (\...
- 143
0
votes
1
answer
401
views
How to deal with transition elements in adaptive fem
It is necessary for me to solve a Poisson equation by adaptive finite element method with transition elements technique to get conforming mesh.
For the first local refinement everything is OK, ...
- 523
1
vote
1
answer
490
views
Non-uniform finite difference Adaptive Mesh Refinement
Assuming that the crosses in the figure below are unknowns in a vertex-centered finite difference scheme in an Adaptive Mesh, how can I calculate the double derivate (Laplacian) at the Red x ? The Red ...
- 467
2
votes
2
answers
264
views
Finite element results by different meshes
There are some technique to generate mesh in a domain. My qustion is that:
Is there any difference between the results using different techniques for mesh generation?
If yes which one is better. For ...
- 523
3
votes
1
answer
177
views
Refinement in AMR
Assume I start with an 8x8 coarse mesh (see Fig 1) where the vertices (except boundary vertices) represent the unknown variable.
After iterative approximation - I ...
- 467
2
votes
1
answer
1k
views
Implementing Finite Difference Adaptive Mesh Refinement code
For a start I need to implement a 2-D, vertex centered, finite difference scheme, Adaptive mesh refinement serial code. I have the following doubts before starting:
Is the input to AMR (say 2-D) ...
- 467
2
votes
1
answer
105
views
Strategies for controlling number of new elements in adaptive mesh refinement
I am working on adaptive techniques for solving some elliptic equations.
The technique is based on residual on elements. My problem is that when I use a predefined tolerance for refining elements, the ...
- 523
5
votes
1
answer
217
views
Using finite element error estimators for adaptive mesh refinement
I am in the process of implementing adaptive mesh refinement for a finite element code that solves the Poisson equation. I have had some trouble finding good references on deciding which elements to ...
- 1,869
5
votes
1
answer
633
views
unstructured grid AMR
Are there libraries for conducting parallel AMR on an unstructured grid ? For a finite volume code, polyhedral cells with arbitrarily shaped faces are as easy to handle as hexahedra, and infact ...
- 233
5
votes
2
answers
1k
views
How to implement adaptive mesh refinement using conformal triangles
I am trying to implement adaptive mesh refinement for a finite element code. The code uses (at least for now) linear triangles and so when I do the mesh refinement I want the triangular mesh to remain ...
- 1,869
4
votes
1
answer
411
views
Finite element mesh software
I'm looking for a program to obtain meshes to finite element codes 2D and 3D as complete as possible, preferably in fortran 90 or C/C++.
For example, softwares "Triangle" or "TetGen" generate meshes ...
- 515
0
votes
1
answer
136
views
How to define fluxes for two dimensional convection-diffusion equation?
I want to solve the following differential equation using control volume approach on a Cartesian mesh:
$$\frac{\partial T}{\partial t} + \frac{\partial T}{\partial x} + \frac{\partial T}{\partial y}= \...
- 414
2
votes
0
answers
411
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 ...
- 121
5
votes
1
answer
231
views
Flux at coarse-fine mesh grid interface?
I am trying to solve one dimensional inviscid Burger's equation using adaptive mesh refinement. This is the PDE:
$$\frac{\delta U}{\delta t} + \frac{\delta F}{\delta x} = 0$$
where the flux F of the ...
- 414
4
votes
2
answers
2k
views
How to discretize Burger's equation?
I am trying to solve the very simple one dimensional burgers equation which is:
$$\frac{\delta U}{\delta t} + \frac{\delta F}{\delta x} = 0$$
where the flux F of some variable U is defined as$$ F= \...
- 414
10
votes
0
answers
226
views
Time advance in Adaptive Mesh Refinement method
I am working on solving complex system of 2D PDEs governing the behaviour of plasma in a gas lamp during discharge. Recent tests have shown that because of steep gradients in temperature field and ...
- 131
1
vote
0
answers
133
views
Conservation at grid interface in adaptive mesh refinement
I am using adaptive mesh refinement to solve one dimensional inviscid Burgers equation. However I am facing some difficulty to handle grid interfaces which are not uniform (coarse-fine grid interface)....
- 414
0
votes
1
answer
288
views
Adaptive mesh refinement basic conceptual problem
I am a beginner in adaptive mesh refinement (AMR). After I am done with the first two papers by Dr. Marsha Berger, I was trying to write my own code for a problem which has a parabolic partial ...
- 414
4
votes
1
answer
96
views
How can I compare errors in PDE solvers with non-uniform grids?
Is there a standard approach to testing codes with refined regions? Specifically, I am interested in testing whether the refinement is working correctly.
For the sake of simplicity, let's consider a (...
- 173
4
votes
2
answers
276
views
Adaptive mesh refinement algorithms and the difference between AMR and moving mesh
I'm working on my thesis and a part of it has to do with adaptive mesh refinement. As a computer science major, I'm not too familiar with this field. The best way I can put my knowledge of AMR is: I ...
- 115
3
votes
1
answer
107
views
Fast methods to solve an elliptic PDE if high accuracy is needed only in part of the domain
Does someone know a method to get cheap approximation of harmonic problems (and possibly local approximations)?
Let me explain: I need to compute the solution of an harmonic problem
\begin{equation} ...
- 455
2
votes
0
answers
108
views
Adaptive Mesh vs Uniform Mesh for multiple source/boundary/initial Data
I'm going to ask some beginners' questions.
Adaptive mesh can save many DOFs than a uniform mesh. But it also needs to solve linear systems changing with mesh adaptive process. Is this not problem?
...
- 1,319
4
votes
1
answer
372
views
Meshes in codes that do AMR
In light of this question and some stuff I read online I am wondering if large FE libs (e.g., deal.ii, libmesh etc.) that do AMR keep the entire mesh or possibly a coarse version of the entire mesh on ...
- 1,759
9
votes
3
answers
3k
views
Finite Element Method vs Extended Finite Element Method (FEM vs XFEM)
What are main differences between FEM and XFEM? When should we (not) use XFEM intead of FEM and vice versa? In other words, when I meet a new problem, how I can know to use which one of them?
- 791