Question about the demo thermal_mechanical_full

[hw-wu]

Hello everyone,

I've been exploring this code with a keen interest in solving coupling problems. The code successfully passed the benchmark tests and have been able to run demos.hyperelasticity.example without any issues. However, I've encountered a stumbling block while trying to work with the thermal_mechanical_full demo. I'm facing a FileNotFoundError as detailed below:

FileNotFoundError: [Errno 2] No such file or directory: '/home/myname/Desktop/jax-fem/demos/thermal_mechanical_full/input/numpy/points.npy'

I'm wondering if there's something I might miss. Any guidance or suggestions would be greatly appreciated!

[tianjuxue]

@xwpken Please also push the input data files to Github if they are not too large. Thanks!

[xwpken]

@hw-wu Hello, we have uploaded the necessary data files. If you have installed the FEniCS, you may also run the 'fenics.py' in the same folder to genereate these files.

[hw-wu]

@xwpken I really appreciate your reply. It seems that FEniCS is used to generate input file... However, when I installed mshr, it said

File "", line 198, in _run_module_as_main
File "", line 88, in _run_code
File "/home/hongwei/Desktop/jax-fem/demos/thermal_mechanical_full/fenics.py", line 16, in
domain = Rectangle(Point(0., 0.), Point(L, L)) - Circle(Point(0., 0.), R, 100)
^^^^^^^^^
NameError: name 'Rectangle' is not defined

And after searching, I found Rectangle() and Circle() functions have been deprecated (https://fenicsproject.discourse.group/t/rectangle-and-circle-functions-deprecated/3845).
So we should use an old version of mshr to generate it?

[tianjuxue]

@xwpken I really appreciate your reply. It seems that FEniCS is used to generate input file... However, when I installed mshr, it said

File "", line 198, in _run_module_as_main File "", line 88, in _run_code File "/home/hongwei/Desktop/jax-fem/demos/thermal_mechanical_full/fenics.py", line 16, in domain = Rectangle(Point(0., 0.), Point(L, L)) - Circle(Point(0., 0.), R, 100) ^^^^^^^^^ NameError: name 'Rectangle' is not defined

And after searching, I found Rectangle() and Circle() functions have been deprecated (https://fenicsproject.discourse.group/t/rectangle-and-circle-functions-deprecated/3845). So we should use an old version of mshr to generate it?

You don't need to run FEniCS yourself. I think Weipeng already uploaded the necessary input data to the latest repo. Please pull from Github.

[hw-wu]

Thanks for your instant reply. Yeah, you are right. The example.py can run now.. but I am trying to understand how to generate the mesh like @xwpken did in the demo...

[xwpken]

Thanks to Prof Xue for his explanation. For using 'fenics.py' to generate data files, you can first check the version of FEniCS. For my Linux machine, the version of FEniCS (including mshr) is 2019.1.0. In the latest FEniCSx, many APIs have changed, so running errors may occur.

[hw-wu]

Thanks for your reply. It can run now. But I also have a question. It seems that the problem is solved twice by FEniCs and JAX-FEM? I am wondering whether it is possibile to use other mesh generation tool like gmsh?