# 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
Filter by
Sorted by
Tagged with
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....
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. ...
1 vote
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....
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
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]...
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 ...
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 ...
259 views

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 ...
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 ...
1 vote
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 ...
102 views

### NDEigensystem returns reversed list of eigenfuctions

I'm considering the two coupled linear differential operators ...
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 ...
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 ...
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 ...
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 ...
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 ...
87 views

### Understanding NeumannValue

From this interesting question Reveal the formal PDE of FiniteElement ...
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 ...
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 ...
75 views

80 views

### Documentation of NDSolveFEMFEMNintegrate

Where can I find detailed information concerning the NDSolveFEM package? Thanks!
128 views

### Using Neumann boundary conditions for the wave equation

I have the following code to solve the wave equation in 2D: ...
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 ...
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 ...
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 ...
82 views

### If or Piecewise dependent coefficients in NDSolveValue

This: but the following (placing the partials inside the If) does not. Any suggestions?
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 ...
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 ...
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 ...
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 ...
213 views

### Quarter a given triangular meshregion

For a given , examplary simple, triangular meshregion ...
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 ...
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. ...
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 ...
218 views

141 views

### Fast integration of ElementMeshInterpolation[]-functions

In this simplified example I consider a small FEM-mesh ...
188 views

### How to create cylinder pie using ParametricRegion?

I'm wondering that this simple definition of a cylinder pie doesn't work ...
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
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 ...
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