Questions tagged [finite-element-method]

Usage of the Finite Element Method embedded in NDSolve and details on the implementation of the fem in mathematica.

Filter by
Sorted by
Tagged with
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 ...
Snowymint's user avatar
0 votes
1 answer
39 views

Numerical instability due to convection dominated PDE [duplicate]

I am trying to solve this partial diffusion equation shown ...
Snowymint's user avatar
1 vote
1 answer
143 views

Numerical solving diffusion equation in spherical coordinates

Mathematica nicely solves Poisson's equation in spherical coordinates as ...
Rodion Stepanov's user avatar
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" ...
Ulrich Neumann's user avatar
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 ...
AminD's user avatar
  • 21
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 ...
Ariana Fenris's user avatar
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 ...
S. Maths's user avatar
  • 203
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} ...
Daniel Castro's user avatar
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 ...
pp.ch.te's user avatar
  • 191
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 ...
Jole Stock's user avatar
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$): ...
anderstood's user avatar
  • 14.2k
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 $\...
JamesVR's user avatar
  • 335
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 ...
kotozna's user avatar
  • 307
1 vote
1 answer
107 views

Import in Mathematica does not work

I want to import STL files as boundary element mesh. ...
Lewis's user avatar
  • 11
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 ...
Brownian_Motion's user avatar
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 ...
Tayfun's user avatar
  • 31
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 ...
Takashi Miura's user avatar
0 votes
1 answer
71 views

Which approximation functions does NDSolve take when solving second order ODE with the finite element method? Are they linear in parts or some other?

This is the code that I have ...
klaic710's user avatar
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 ...
user94537's user avatar
  • 107
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 ...
user94537's user avatar
  • 107
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: ...
kacper's user avatar
  • 21
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: ...
user94537's user avatar
  • 107
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 ...
Ulrich Neumann's user avatar
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: ...
Peter's user avatar
  • 21
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 ...
Hugh's user avatar
  • 16.1k
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(...
Andrei Nitescu's user avatar
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 ...
icebox207's user avatar
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 ...
Wd Wd's user avatar
  • 11
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 ...
user84456's user avatar
  • 1,524
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 ...
Frida's user avatar
  • 45
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 ...
Fun123's user avatar
  • 47
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: ...
RuiQiang Kang's user avatar
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: <...
Julio Araujo Dos Santos's user avatar
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 ...
Jack LeGrow's user avatar
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 ...
Ashot Matevosyan's user avatar
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 ...
user93372's user avatar
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 ...
A.D's user avatar
  • 41
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 ...
Mam Mam's user avatar
  • 1,753
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 ...
Mam Mam's user avatar
  • 1,753
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 ...
Anthony's user avatar
  • 155
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 ...
Mam Mam's user avatar
  • 1,753
2 votes
0 answers
73 views

Why am I getting a parsing error in NDSolveValue?

I am not sure why I am getting a ...
Charlie's user avatar
  • 506
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}$. ...
Mam Mam's user avatar
  • 1,753
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 ...
Eric Hester's user avatar
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}...
Charlie's user avatar
  • 506
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}\...
Charlie's user avatar
  • 506
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; ...
sra's user avatar
  • 697
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: ...
Alexei Boulbitch's user avatar
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): ...
nikiniki's user avatar
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 ...
Mam Mam's user avatar
  • 1,753

1
2 3 4 5
29