Questions tagged [finite-element-method]
Usage of the Finite Element Method embedded in NDSolve and details on the implementation of the fem in mathematica.
1,442
questions
0
votes
0
answers
34
views
Solving PDE for Diffusion Equation (Boundary Condition Issue) [duplicate]
I am trying to solve this partial diffusion equation shown
...
0
votes
1
answer
39
views
Numerical instability due to convection dominated PDE [duplicate]
I am trying to solve this partial diffusion equation shown
...
1
vote
1
answer
143
views
Numerical solving diffusion equation in spherical coordinates
Mathematica nicely solves Poisson's equation in spherical coordinates as
...
4
votes
2
answers
252
views
How to mesh a cylinder with helix points?
I would like to mesh a cylinder surface. The mesh should include given cylinderpoints which lie on a helix.
My attempt using "IncludePoints"
...
2
votes
1
answer
133
views
Numerical ground state wavefunction of Schrödinger equation with a Coulomb potential in 2D from NDEigensystem
I want to numerically solve the ground state wave function of the hydrogen atom with the Coulomb potential using the NDEigensystem. Here is the code to get the ground state wave function from the ...
7
votes
4
answers
278
views
Minimal surface bounded between turns of helix
I'm curious to find the shape of a surface bounded between the rungs of a helix, ie the shape of the cloth stretched between the rungs of this child's play tunnel. I'm wondering if we could find it ...
2
votes
0
answers
134
views
Problem with pdetoode for two coupled PDEs
I tried to adapt a code for a single equation to solve the following system using 'pdetoode'
Updated answer
...
3
votes
1
answer
143
views
Solving system of first order PDEs
I am trying to solve the 1st order PDE system
\begin{align}
\xi_u^2+\eta_u^2&= \left(1+\frac{\xi^2+\eta^2}{4} \right)^2\\
\xi_v^2+\eta_v^2&=\left(1+\frac{\xi^2+\eta^2}{4} \right)^2
\end{align}
...
4
votes
1
answer
255
views
Numerical solution of second-order linear hyperbolic PDE
(I'm also searching for analytical solutions to this PDE; check the bountied questions here and here if you have any ideas)
I'm trying to find the numerical solution of the following 2D second-order ...
2
votes
1
answer
196
views
Machine overflow when defining boundary conditions
Recently I have been trying to code Maxwell's equations over a closed surface and have been facing some trouble defining the boundary conditions for the magnetic field.
The equation for the normal of ...
2
votes
1
answer
100
views
Solving Poisson PDE with NDSolve and incomplete BC specifications
When solving the following PDE with a missing BC on the fourth edge ($y=1$):
...
3
votes
4
answers
167
views
NDSolve over derivatives of Heaviside function
I am trying to numerically solve for the strain of a Maxwell material in response to a step stress. The governing equations are $$\dot{\sigma} + \sigma = \dot{\varepsilon}$$ and I want to find $\...
3
votes
1
answer
112
views
Numerically solve for the potential of a point charge in a periodic cubic domain
How can I obtain the 3D numerical solution for the potential (or field) due to a point charge inside a cubic domain with periodic boundary conditions in all directions? I guess I can use ...
1
vote
1
answer
107
views
Import in Mathematica does not work
I want to import STL files as boundary element mesh.
...
1
vote
1
answer
70
views
NDSolveValue::fememrc | Defining initial conditions on a part of a 1D mesh
I am trying to solve a simple PDE in the cylindrical coordinates:
$$ \frac{\partial c}{\partial t} = D \bigg(\frac{1}{r}\frac{\partial c}{\partial r}+\frac{\partial c^2}{\partial r^2}\bigg)$$
The ...
3
votes
1
answer
167
views
Numerical solution of the Richards' equation
I am trying to solve Richards' equation to model fluid flow in soil. The governing partial differential equation, initial condition, and boundary conditions are:
The analytical solution of the ...
6
votes
2
answers
298
views
Simulating buckling
We are trying to implement buckling using a newly implemented FEM solver. However, if we try to reproduce the buckling phenomena using a thin rod, it is just compressed, and we cannot observe the ...
0
votes
1
answer
71
views
2
votes
1
answer
83
views
Post Processing of Solution and Plot of Coupled Partial Differential Equation Over a Semi-Circular Domain
Linked Question: Solution and Plot of Coupled Partial Differential Equation Over a Semi-Circular Domain.II
Many many thanks to @Nasser (https://mathematica.stackexchange.com/users/70/nasser) for his ...
2
votes
1
answer
116
views
Solution and Plot of Coupled Partial Differential Equation Over a Semi-Circular Domain.II
This question is connected with reference to my previous question (Solution of Dimensionless Partial Differential Equation Over a Semi-Circular Domain), where I asked for a help to solve the following ...
1
vote
1
answer
104
views
NDSolveValue result contradicts initial condition
I'm trying to solve the Schroedinger equation in a box with two holes defined by:
...
1
vote
1
answer
85
views
Solve dimensionless PDE in polar coordinate over a semi-circular re
I am trying to solve the following BVP in mathematica
using the following code:
...
2
votes
1
answer
72
views
FEM solution for onedimensional boundary value problem doesn't evaluate [closed]
modified
Based on this, unfortunately closed, question How can we know on how many parts is domain splitted with MaxCellMeasure? [closed] I would like to know why ...
2
votes
1
answer
89
views
Mathematica thinks that an initial condition is a boundary condition
I would like to solve the following reaction-diffusion problem in Mathematica using NDSolve:
...
5
votes
1
answer
120
views
Artifacts when interpolating on an unstructured set of 3D data
When doing an interpolation on an unstructured set of points it is best to use the interpolation method available from the finite element package. However, it is essential that some re-scaling is done ...
3
votes
1
answer
120
views
Second order Poisson ODE
I have a 2nd differential Poisson equation for the electric field in a semiconductor (application for a MOS simulation, hopefully).
I have the formula:
with the condition that after a certain value(...
2
votes
1
answer
83
views
Need help solving cylindrical Laplacian
I'm trying to solve the cylindrical Laplacian for a heated disk in a large cylinder. The cylinder and disk have constant temperatures and I only care about the temperature field between the disk and ...
1
vote
0
answers
47
views
WhenEvent in ParametricNDSolve but not executed
I tried to solve a partial differential equation with parameters, the critical condition of which I used WhenEvent to express, but the solution after bringing in the parameters shows that the ...
2
votes
1
answer
210
views
Solving Schrödinger equation for Dirac comb potential (kicked rotor)
I need to solve the Schrödinger equation for a Dirac delta potential. I could not find the correct way to write the time-dependent potential and how to solve the time-dependent equation for it.
The ...
2
votes
1
answer
64
views
AceGen: Assemble a user specific global array (similar to "Residual")
I have an issue which I know is definitely solvable but I don't know how to implement this in AceGen.
Basically, I want to assemble a matrix $\mathbf{L} \in \mathbb{R}^{n_{dof} \times 6}$ for a global ...
1
vote
1
answer
155
views
How to find eigenvalues?
I have the following system of differetial equations, which are linearized based on steady state solutions. How to discretize them and find eigenvalues or use any other method without discretization ...
1
vote
1
answer
263
views
How to find the analytical or numerical solution of a system of partial differential equations?
I am trying to use mathematica to solve a system of second-order partial differential equations, but I have been unable to solve it. The code of mathematica is as follows:
...
2
votes
2
answers
204
views
Spherical Heat Equation and Convection Boundary Conditions
I'm trying to solve transient heat conduction equation in spherical domain in Mathematica. I made the simplifying hypothesis that temperature varies only with time and radial coordinate.
The code is:
<...
1
vote
1
answer
71
views
Fitting Parameters to the Heat Equation [closed]
As part of my research, I have been trying to use a model of heat conduction through a 2D layer given an input steady-state Gaussian power profile and a heat loss term to the environment. All but one ...
0
votes
0
answers
43
views
NDSolve diverges with Neumann boundary condition [duplicate]
I am solving Poisson-like PDE with the Finite Element Method in Wolfram Mathematica. Only the Neumann boundary condition is imposed on the boundary. Of course, the solution is not unique, most likely ...
4
votes
1
answer
115
views
Finer mesh for selected subregion of a solid 3D cylinder
I am trying to create a finer mesh in a subregion of a solid cylinder in order to improve the resolution there. I don't want to have too many elements for the entire cylinder because it slows down the ...
4
votes
1
answer
111
views
FindRoot error in NDSolve
I have two coupled differential equations as follows
$$
\frac{\partial }{\partial x}U(x,y) =2V(x,y),
$$
$$
\frac{\partial }{\partial y}V(x,y) =V(x,y)U(x,y)+1,
$$
with ...
0
votes
1
answer
131
views
Why do the eigenvalues periodically change with successive increase in the consideration region?
When finding the eigenvalues and eigenfunctions of the system Hc[r, z] using NDEigensystem, the following issue arises: When ...
1
vote
0
answers
72
views
How to find the eigenvalues of a particle in the Coulomb field using NDEigensystem?
As advised in the previous two questions (Why, if I enter an angle into the function, then does the code not work correctly?, let us solve the problem for the Coulomb potential in spherical ...
1
vote
0
answers
94
views
On formulating a Neumann boundary condition
I am attempting to follow this tutorial in the documentation on using FEM to solve PDEs. I am having difficulty understanding how to formulate the Neumann boundary condition for my free-boundary ...
0
votes
0
answers
47
views
Eigenvalues and eigenfunctions of a parameter dependent system
There is the following system $H=-\frac{1}{2}\Delta-a\frac{1}{r}$, where a is a parameter. I'd like to find eigenvalues and eigenfunctions depending on the ...
2
votes
0
answers
73
views
Why am I getting a parsing error in NDSolveValue?
I am not sure why I am getting a ...
0
votes
0
answers
49
views
How can NDSolveValue and ParametricNDSolveValue be used correctly to find eigenvalues and eigenfunctions?
NDEigensystem can be used to find the eigenvalues and eigenfunctions of the following system $H=-\frac{1}{2}\Delta-\frac{1}{r}$.
...
2
votes
0
answers
52
views
Can Mathematica's FEM solve coupled Dirichlet Boundary Conditions?
I am solving a coupled system of PDEs using Mathematica's FEM capabilities.
Specifically, the Navier-Stokes equations with a no-flux stress-free boundary.
To do this, I need to specify a Dirichlet ...
1
vote
1
answer
130
views
Why does NDSolve solution not match physical intuition?
Related to a previous question I posted, I am trying to solve for the electric potential on a box which obeys Ohm's law:
\begin{equation}\nabla\cdot(\overset{\scriptscriptstyle\leftrightarrow}{\sigma}...
1
vote
1
answer
147
views
How to calculate a numerical Fourier transform obtained from NDSolve?
We numerically solve for the electric potential on a box which obeys Ohm's law:
\begin{equation}\nabla\cdot(\overset{\scriptscriptstyle\leftrightarrow}{\sigma} \nabla \Phi) = 0, \label{Eqn:OhmsLaw}\...
1
vote
2
answers
207
views
Normalization of solution of a PDE
I am trying to solve the following PDE:
pde = D[P[x, t], t] + D[J[x, t], x] == 0;
...
3
votes
3
answers
274
views
Area, perimeter and other geometric parameters of a "cucumber"
From the simulation, I have got a region with a shape of a cucumber. Here are the boundary points of this region:
...
1
vote
1
answer
103
views
How to solve system of PDEs (e.g. Navier-Stokes)
I need to solve system of PDEs.
The system of PDEs:
Boundary and initial conditions:
.
I used NDSolve but have some troubles (code below):
...
2
votes
1
answer
93
views
Why does NDEigensystem not show the minimum eigenvalue for a certain parameter range in the cylindrical coordinate system?
In my previous question Why NDEigensystem does not show the minimum eigenvalue?, I asked why the NDEigensystem does not show the minimum eigenvalue for the ...