Hello there. Actually, i’m ending to solve my problem, but i’m founding some issues to plot the solutions: from __future__ import print_function from petsc4py import PETSc from fenics import * import matplotlib.pyplot as plt from IPython.display import clear_output import numpy as np import time fr…

Could you please state what you expect to get, and show what plots you are currently getting? As far as I can tell from looking at your code F1 and F2 are residuals. These can be associated with the finite element space of the test function u1 and v1, and plotted over space for each time step. Is …

Hello, sorry for this. what I’m getting is: [image] [image] [image] There is a lot of lines and I don’t know why, it should be the single accumulated energy at each time step, That is, the total elastic, fracture and total energy at each time step. Also, it seems like that is not adding both e…

[image] janduigamer13: elastic_energy = assemble(F1) fracture_energy = assemble(F2) elastic_energy_values.append(elastic_energy) fracture_energy_values.append(fracture_energy) Start by considering these variables. What do you see when you print them? I would expect these to be a …

Ok, I’ve been thinking for a while and I don’t know if it’s correct, but with your idea, I’ve got: [image] [image] [image] I summed all values and, as a good news, I’ve got only positive terms: Cumulative Elastic Energy: 0.011589729018203784 Cumulative Fracture Energy: 2.1286995202790905e-15 …

[image] janduigamer13: Where can I found more about get_local() functions? See Efficient way to convert data from numpy to fenics Function and documentation .

Brief doubt with plotting fenics solution

Legacy DOLFIN

violetus December 19, 2023, 10:05pm 6

See Efficient way to convert data from numpy to fenics Function and documentation.

Topic		Replies	Views
How to extract and plot the eigen-function dolfinx	16	960	October 16, 2023
Error in NewFrac FEniCSx tutorial : ValueError: too many values to unpack (expected 2) General	5	568	September 19, 2024
Why the solution is not appearing	6	161	May 4, 2024
Plotting energy on mesh interval General	9	318	June 8, 2023
How can I plot the result of the velocity using plot function in vedo.dolfin instead of using paraview	3	543	June 2, 2022

Brief doubt with plotting fenics solution

Related topics