The Stack Overflow podcast is back! Listen to an interview with our new CEO.

Questions tagged [finite-element-method]

Tag for the usage of "FiniteElement" Method embedded in NDSolve and implementation of finite element method (fem) in mathematica.

38 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
573 views

From where to learn finite element method?

Can anyone reccomend me a book or site for learning finite element method with mathematica, besides wolfram language official site, https://reference.wolfram.com/language/FEMDocumentation/tutorial/...
297 views

Create vtk file of ParaView from Mathematica [FEM post-processing]

I'm doing FEM in Mathematica and want to visualize the results in Paraview (e.g., create a video of the displacement field). Now I have the element mesh, the displacement and stress fields. How can I ...
496 views

benchmark of Mathematica's FEM?

Does anyone have any numbers on how mathematica compares to other commercial (eg ANSYS) or free FEM softwares (eg. FreeFEM++, FEniCS, elmer)? If this is too vague say solving the diffusion equation in ...
407 views

FEM Mesh for pipes and tubes

I work with pipe systems and would like to use finite elements within NDSolve for regions inside or exterior to a pipe. How to do I make a mesh of the region ...
316 views

Stress Operators for Finite Element Analysis

For stress analysis in the finite element method you need a stress operator. For two dimensional plane stress it can be found here. For two dimensional plain strain it can be found here. I need the ...
292 views

How to use periodic boundary conditions on the derivative of u[x,y,z]?

I have a structural mechanical problem, where i try to calculate the strain field in the modelled material. I simplified the problem in the following example: There is a material with two phases, ...
420 views

Symbolic Weak Form

Usually I write the weak form by hand for my FEM code, but it's a little annoying and mechanic sometimes. So, I wonder, is there any way to generate the symbolic weak form in Mathematica? For ...
79 views

Direction of NeumannValue on boundary

When solving partial differential equations NeumannValue is used to specify the flux across the boundary. For normal fluxes the detailed notes state They appear on the boundary ∂Ω of the region Ω ...
315 views

How to increase boundary mesh quality?

Update: I tracked down the current issue to errors generated by numerically integrating over the boundary mesh. What's the proper way to integrate a function over the boundary of an implicit region? ...
145 views

NDSolve post-processing: Calculate the flow over a FEM-boundary

I'm trying to calculate the heat loss from buried pipes with NDSolve`FEM. For this I set Dirichlet conditions at the model top and bottom and temperature dependant Neumann conditions at the contact ...
52 views

Interface points of NDEigensystem

When solving an eigenvalue problem with "NDEigensystem", e.g. a 1D Eigenvalue problem with the interval composed of different materials, which should be solved by the pure numerical method such as FEM,...
143 views

Postprocessing of FEM: extremely slow speed of NIntegrate

I've gone through the documentation for NDSolve FEM. There are quite a lot tips and techniques for accelerating NDSolve, and only quite few tips for post processing. Somehow I found out that the ...
183 views

Reducing a FEM 3D Interpolation Function to a 2D Interpolation Function

I'm trying to solve the heat equation in 3D using FEM and back out a slice in the xz plane. My problem is when I define a function that is a slice in the xz plane it uses the full 3D interpolation ...
124 views

Stress analysis in a disk with circular holes

Writing the following code: ...
128 views

Coupled parabolic differential equations with time delay

Is it possible for NDSolve to solve delay partial differential equations with simple Neumann boundary conditions? An example I have is as below: ...
953 views

FDTD Electromagnetic Simulations with Mathematica

Mathematica has some Finite Element capabilities (as explained here). Do you think mathematica is a realistic/sensible option for electromagnetic simulations? (aka making a movie of electric fields ...
268 views

How to Implement a “Function Object” suitable for NDSolve and Finite Element Method?

I have a bunch of methods that shoud return an expression that behave like a function, i.e. an expression that, when supplied with suitable arguments, evaluate to something. A simple ...
258 views

How can I improve the solution of my PDE?

I want to improve the last code to find a better solution for umi. For now I'm using only MaxStepFraction. Also, when I put no ...
503 views

221 views

Puzzling NDSolve[] behavior for PDE (smooth solution, inconsistent with boundary conditions)

Consider the following: NDSolve[{D[z[x, y], x, x] + D[z[x, y], y, y] == 0, z[x, 0] == Sin[x], z[0, y] == Cos[y]}, z[x, y], x, y] {{z[x, y] -> ...
85 views

Problems with the eigenvalues calculated using NDEigenvalue

I'm trying to solve a Sturm-Liouville problem like $\qquad -\psi''(z)+(\frac{1}{z}+2\,z)\psi'(z)=\lambda\,\psi(z)$ using NDEigensystem in order to learn how to ...
48 views

AceFEM: Manipulation of SMTData[“MatrixGlobal”] after building but before solving

I want to manipulate the global stiffness matrix before I solve the system of equations. Unfortunately, the command SMTNewtonIteration[] builds up the global matrix and the global vector and solve the ...
110 views

PeriodicBoundaryConditions an a square

I hope the answer is simple, but I can't see the solution: I want to solve a PDE for u[x,y] under Periodic boundary conditions on say a square size 10 x 10. I.e. <...
81 views

How to draw a triangle shape function plot with shape function area being (s, t,1 - s - t)

I want to draw a dynamic graph showing the varying shape of a triangular element with the change in the point P[x, y] as shown below. I want to find the shape ...