Questions tagged [finite-element-method]
Tag for the usage of "FiniteElement" Method embedded in NDSolve and implementation of finite element method (fem) in mathematica.
1,221
questions
0
votes
0
answers
39
views
Meshing Error on Imported obj file
I am trying to import a .obj file to be used in a static electric field simulation. I am able to successfully run the simulation with certain other aperture geometries (see here: < https://www....
- 23
0
votes
0
answers
55
views
Simulating the Response of Charged Particles to Apertures with Inhomogeneous Electric Fields Passing through them
I have produced a simulation which allows me to observe a static electric field produced by a charged plate as it passes through some complex apertures in a different, oppositely charged metal plate.
...
- 23
1
vote
2
answers
71
views
Defining Dirichlet Boundary Conditions on an Imported 3D Object with Apertures in it
I am trying to simulate a static electric field (in 2d and 3d) as it passes through a couple of apertures in a 3D "metal sheet" (model imported as a .obj file) given here:
https://www....
- 23
-1
votes
0
answers
50
views
How do I apply bending moments to a THREE-DIMENSIONAL finite element model in Mathematica [closed]
How do I apply bending moments to a THREE-DIMENSIONAL finite element model in Mathematica
- 1
5
votes
1
answer
67
views
NDEigensystem: 1D problem with discontinuous coefficients
I am trying to use NDEigensystem to solve the 1D problem
-cs[x]^2 vx''[x] = w^2 vx[x]
with vx[x] and w the eigenfunction and eigenvalue. The coefficient cs[x] is discontinuous at x = -xp, +xp and vx[x]...
- 85
0
votes
1
answer
118
views
How to solve 2D Heat equation on rectangle with NDSolve?
I am trying to solve the 2D heat equation $$\Delta u =u_t$$ with Dirichlet boundary conditions on the rectangle $(0,a)\times (0,b)$ and with initial condition $u(x,y,0)=g(x,y)$.
This is my code
...
- 125
0
votes
0
answers
52
views
Create a voxelized mesh for a Spherical shell (using a STL mesh file)
I would like to be able to voxelize a spherical shell. I have tried a code that was posted some years ago in voxelize. But it generates the mesh of the entire sphere and not the shell.
I would ...
- 489
6
votes
2
answers
259
views
NDSolve for Laplace equation on disk is not working
I want to solve the laplace equation on the disk using NDSolve.
Particularly I want to solve the problem $$\frac{\partial^2 u }{\partial r^2 } +r^{-1}\frac{\partial u }{\partial r }+r^{-2}\frac{\...
- 173
0
votes
1
answer
88
views
PDE with tangential derivative condition along boundary. NDSolve not working
I'm considering the following PDE
$$
\begin{align}
\left[\cos\theta\frac{\partial}{\partial r}-\frac{\sin\theta}{r}\frac{\partial}{\partial\theta}\right]\left[\sin\theta\frac{\partial}{\partial r}+\...
- 555
8
votes
1
answer
352
views
Droplet of water-alcohol mixture spreading due to evaporation-caused surface tension gradient
I am trying to solve a coupled system of PDEs for 2 functions h[t,r] and c[t,r], with initial conditions of ...
- 81
4
votes
0
answers
100
views
FEM: modelling 2D flow under periodic forcing
Despite many examples of periodic boundary conditions being presented, it remains difficult to obtain a periodic solution. Let us consider a simple example of a flow in a two-dimensional channel under ...
- 783
1
vote
0
answers
48
views
Neumann Boundary Condition for thermal radiation between two bodies
I'm trying to formulate and solve a pde which models thermal radiation between two bodies using a Neumann boundary condition.
I've created a mesh and want to capture the thermal radiation heat ...
2
votes
1
answer
102
views
NDEigensystem returns reversed list of eigenfuctions
I'm considering the two coupled linear differential operators
...
- 783
6
votes
1
answer
151
views
NeumannValue is producing incorrect results
Cross-posted in Wolfram community, an online Mathematica notebook is available there.
The problem under investigation is solving the Laplacian equation to model seepage. See the image below for the ...
- 63
3
votes
1
answer
117
views
Speeding up Boussinesq equations solving?
I am working on Boussinesq equation. The notebook can run perfectly for only 0.6 steps and then the calculation starts running slowly after 0.7. All boundary conditions seemed fine. I am unsure if I ...
- 33
2
votes
1
answer
53
views
Convert boundary condition involving derivative to NeumannValue programmatically
As discussed in e.g.
what causes the error "The dependent variable in the boundary condition needs to be linear" when using NDSolve?
When FiniteElement ...
- 52.5k
2
votes
2
answers
107
views
How to solve 2D eigenvalue problem with robin boundary conditions
I need to solve an eigenvalue problem in 2D as seen in the picture.
I've tried the function NDEigensystem but reading its documentation it seems it has issues with ...
- 21
10
votes
2
answers
815
views
How can I generate a Swiss-Cheese type region and render it transparent?
I am trying to visualize a percolation model that is equivalent to a swiss-cheese where the holes are conducting. At some point, the size of the holes gets large enough for them to touch and a ...
- 101
3
votes
0
answers
87
views
Understanding `NeumannValue`
From this interesting question Reveal the formal PDE of FiniteElement
...
- 37k
5
votes
1
answer
228
views
Neumann boundary conditions are not satisfied
I am solving a system of 3 ODEs and I want to implement neumann BC (d/dx = 0 @x=0). However, the solution does not seem to obey neumann BC when parameter A1/A2/A3 is non zero, as I notice that the ...
- 105
7
votes
1
answer
329
views
Reproducing figure from paper on hydrodynamics
I am trying to solve Eqn. on Page 674 of https://www.nature.com/articles/nphys3667.pdf I shall quote the relevant paragraph :
Next, we proceed to analyse the nonlocal response in a strip with
...
- 345
4
votes
0
answers
75
views
NeumannValue nonlinear pde?
Inspired by the highly interesting question Solving Burger's equation with NDSolve at large time I'll try to understand Finite-Element-solution for nonlinear pde in more detail:
Burger pde
$u_t(x,t)+u(...
- 37k
4
votes
1
answer
143
views
NDEigensystem with FEM: parasitic solutions
I am trying to find the eigensystem of a system of 3 coupled ODEs. Analytically, the system spectrum should have a gap at [-2M,2M] (except for the degenerate state with E=0, strictly at the middle of ...
- 61
9
votes
4
answers
778
views
Solving Burger's equation with NDSolve at large time
I want to solve the following Burger's equation
$$\partial_tu+\partial_x\left(\frac{u^2}{2}\right)=0,~~x\in[0,2\pi],~t>0\\u(x,0)=\frac{1}{3}+\frac{2}{3}\sin(x)$$
with mathematica. Here's my code:
<...
- 193
2
votes
1
answer
63
views
Inactive form for PDE and symmetry
I would like to solve the following non-linear Poisson equation (as a toy problem for a more complicated problem)
...
5
votes
1
answer
175
views
Coupled PDE using NDSolve
The problem is to study the solution of a biharmonic which can be decomposed into a coupled PDE system
$$
\Delta f - \nabla f\cdot\Gamma_{1} = g \\
\Delta g = \nabla f\cdot\Gamma_{2}
$$where $...
- 345
2
votes
1
answer
80
views
Documentation of NDSolve`FEM`FEMNintegrate
Where can I find detailed information concerning the NDSolve`FEM` package?
Thanks!
- 37k
4
votes
2
answers
128
views
Using Neumann boundary conditions for the wave equation
I have the following code to solve the wave equation in 2D:
...
0
votes
0
answers
157
views
Why is regioncentroid not giving correct results?
I am trying to calculate and plot the x-axis of a volume centroid for the region between zx-plane and a function squared $u(t,x,z)$ from time $t=15$ to $25$. The function is a numerical solution to a ...
- 43
3
votes
1
answer
87
views
How to stop simulation after a given model variable achieves certain value?
I would like to stop my simulation after a given model variable (say, plastic strain averaged over entire domain) achieves a given threshold value. Do you know how to do it?
I'm supplying the MWE. How ...
- 517
0
votes
0
answers
53
views
Why does NDSolveValue give such different answers depending on the "Method" used and without indicating an error? [duplicate]
The first two inputs below only differ in the "Method" used in NDSolveValue for the solution of the heat equation in a composite yet they give no error (the error about BCs and ICs is ...
0
votes
1
answer
82
views
If or Piecewise dependent coefficients in NDSolveValue
This:
but the following (placing the partials inside the If) does not. Any suggestions?
5
votes
1
answer
139
views
Solving system of nonlinear PDEs
My aim is to solve a system of nonlinear PDEs arising in nonlinear elasticity. I am new to Mathematica so I started by modifying this example. I tried to change it to neo-Hookean solid. The resulting ...
- 51
8
votes
3
answers
331
views
How to generate a mesh in an area with curves inside
I would like to know how to generate a mesh inside an area, and be able to find the original curve outline inside the area from the mesh.
For example:
There is a rectangle inside which there is a ...
- 83
3
votes
1
answer
116
views
Simulating a gaussian pulse for the wave equation
In the process of investigating diffraction of waves, I am starting with a simple problem consisting in injecting a gaussian pulse at the boundary x=0 of a square domain and then solving the wave ...
- 45
6
votes
1
answer
229
views
Solving 2D Landau-Khalatnikov equation & Poisson equation using finite element approach
I would like to solve Landau-Khalatnikov (LK) equation and Poisson equation using finite element approach to simulate ferroelectric hysteresis loop, electric potential distribution and electric ...
- 63
8
votes
1
answer
213
views
9
votes
2
answers
137
views
How to do boolean operations between a large amount of regions
I want to do boolean operations between a large amount of regions.
For convenience, I consider a simplified problem.
The problem is as follows, digging out a large number of small hemispheres on the ...
- 175
2
votes
1
answer
80
views
NDSolve fails while solving PDE with Dancwerts BC
The following code runs with BC1 (Dirichlet Type). However, it fails for Dancwerts type BC1.
...
- 21
2
votes
1
answer
97
views
Why does SolidMechanicsStrain calculated inside a module give strange results?
I am exploring the FE capabilities of Mathematica and trying a few solid mechanics problems with the aim of bringing them into my teaching. The problem is a simplified version of the coupled heat ...
- 3,564
4
votes
1
answer
218
views
System of 3 second order non linear differential equations
I wish to solve the following system of equations:
$\frac{d^2f}{dr^2}- \frac{1}{r} \frac{df}{dr} = 2 f(r)\phi(r)^2$
$\frac{d^2\phi}{dr^2} + \frac{1}{r}\frac{d\phi}{dr} = \frac{1}{r^2} f(r)^2\phi(r) + \...
- 73
5
votes
1
answer
149
views
How does FindRoot decide if a solution has converged?
I am solving a 1D non-linear differential equation using the finite element method with NDSolve.
From the documentation I understand that the equation is discretized and then solved with FindRoot,
I ...
- 117
6
votes
1
answer
118
views
Solving a Fokker-Planck equation NDSolveValue
I want to solve the following Fokker-Planck equation:
$$\partial_t p = \frac{1}{2}\partial_{xx}\left[\varepsilon^2p\right] + \frac{1}{2}\partial_{yy}\left[\varepsilon^2p\right] - \partial_x\left[(9x-x^...
- 73
6
votes
1
answer
141
views
Fast integration of ElementMeshInterpolation[]-functions
In this simplified example I consider a small FEM-mesh
...
- 37k
5
votes
2
answers
188
views
How to create cylinder pie using `ParametricRegion`?
I'm wondering that this simple definition of a cylinder pie doesn't work
...
- 37k
6
votes
4
answers
198
views
How to order the coordinates generated by ToBoundaryMesh counterclockwise?
The following code works for polygons. A square 4x4 is tried here:
...
1
vote
0
answers
125
views
How to evaluate an Integral which has Interpolating functions
I am trying to evaluate an Integral which has products of Interpolating functions inside, using NIntegrate command. It gives a wrong answer. The code that I used for performing this task is given ...
- 181
0
votes
1
answer
94
views
Coordinates of mesh generated by Mathematica-FEM
MeshCoordinates function yields more nodal coordinates than expected. Can I extract the coordinates of vertices only. Here is an example of a square region with a course mesh to explain the problem. I ...
1
vote
1
answer
77
views
Dirichlet condition from user-defined polygon
I would like to specify a Dirichlet condition from a polyhedron defined by a set of points.
...
- 145
3
votes
1
answer
88
views
The element mesh has insufficient quality
I am quite new at FEM. I was trying to create a triangular mesh of a rectangular region. I wanted to do something very basic: Discretize the domain using squares and then subdivide each square into ...
- 117