# 3D-Airfoil

**URL:** <https://pyfr.discourse.group/t/3d-airfoil/272>\
**Category:** General\
**Created:** [20 February 2020 15:14 UTC](https://pyfr.discourse.group/t/3d-airfoil/272 "2020-02-20T15:14:19Z")\
**Posts on this page:** 5\
**Page:** 1

<div class="post-metadata">

**Author:** ![Stian\_Hjorteland](https://avatars.discourse-cdn.com/v4/letter/s/b782af/32.png) [@Stian\_Hjorteland](https://pyfr.discourse.group/u/Stian_Hjorteland)\
**Post date:** [20 February 2020 15:14 UTC](https://pyfr.discourse.group/t/3d-airfoil/272/1 "2020-02-20T15:14:19Z")

</div>

Hi guys.

As a part of my master thesis i am simulating flow on an airfoil (90 deg) in openfoam. I also want to compare my results with pyfr (3rd order). Somehow, my setup at the present moment wont start to simulate. I am using Ubuntu on a Windows computer (through Oracle Virtual Box).  
My attempt on a setup for the .ini is attached. It is worth mentioning that the mesh is an O-grid from pointwise with 1000000 cells. The .pyfrm file is therefore made by pointwise when exporting and should be fine which is why i think the problem is in the .ini file.

I would be grateful for any tips as i am a new user of PyFR and find it a bit different from openfoam.

Best regards Stian

[Airfoil-DU96-w-180.ini](https://pyfr.discourse.group/uploads/short-url/b908JndOwiVbYvQqxCPaLVX6bSB.ini) (2.15 KB)

---

<div class="post-metadata">

**Author:** ![Giorgio\_Giangaspero](https://avatars.discourse-cdn.com/v4/letter/g/ce73a5/32.png) [@Giorgio\_Giangaspero](https://pyfr.discourse.group/u/Giorgio_Giangaspero)\
**Post date:** [21 February 2020 09:32 UTC](https://pyfr.discourse.group/t/3d-airfoil/272/2 "2020-02-21T09:32:59Z")

</div>

Hi Stian

do you get any error messages?  
best

---

<div class="post-metadata">

**Author:** ![Stian\_Hjorteland](https://avatars.discourse-cdn.com/v4/letter/s/b782af/32.png) [@Stian\_Hjorteland](https://pyfr.discourse.group/u/Stian_Hjorteland)\
**Post date:** [21 February 2020 09:58 UTC](https://pyfr.discourse.group/t/3d-airfoil/272/3 "2020-02-21T09:58:24Z")

</div>

Hi, yes:

msk600@virtual:~/AirfoilPyFR$ export OMP\_NUM\_THREAD=1  
msk600@virtual:~/AirfoilPyFR$ pyfr run -b openmp -p Airfoil-DU96-w-180.pyfrm Airfoil-DU96-w-180.ini  
Traceback (most recent call last):  
File “/usr/bin/pyfr”, line 11, in   
load\_entry\_point(‘pyfr==1.5.0’, ‘console\_scripts’, ‘pyfr’)()  
File “/usr/share/pyfr/pyfr/ **main**.py”, line 110, in main  
args.process(args)  
File “/usr/share/pyfr/pyfr/ **main**.py”, line 235, in process\_run  
args, NativeReader(args.mesh), None, Inifile.load(args.cfg)  
File “/usr/share/pyfr/pyfr/ **main**.py”, line 216, in \_process\_common  
solver = get\_solver(backend, rallocs, mesh, soln, cfg)  
File “/usr/share/pyfr/pyfr/solvers/ **init**.py”, line 14, in get\_solver  
return get\_integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)  
File “/usr/share/pyfr/pyfr/integrators/ **init**.py”, line 43, in get\_integrator  
return integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)  
File “/usr/share/pyfr/pyfr/integrators/std/controllers.py”, line 12, in **init**  
super(). **init** (\*args, \*\*kwargs)  
File “/usr/share/pyfr/pyfr/integrators/std/steppers.py”, line 165, in **init**  
super(). **init** (\*args, \*\*kwargs)  
File “/usr/share/pyfr/pyfr/integrators/std/steppers.py”, line 8, in **init**  
super(). **init** (\*args, \*\*kwargs)  
File “/usr/share/pyfr/pyfr/integrators/std/base.py”, line 12, in **init**  
super(). **init** (\*args, \*\*kwargs)  
File “/usr/share/pyfr/pyfr/integrators/base.py”, line 59, in **init**  
self.system = systemcls(backend, rallocs, mesh, initsoln, nreg, cfg)  
File “/usr/share/pyfr/pyfr/solvers/base/system.py”, line 34, in **init**  
eles, elemap = self.\_load\_eles(rallocs, mesh, initsoln, nreg, nonce)  
File “/usr/share/pyfr/pyfr/solvers/base/system.py”, line 107, in \_load\_eles  
eles.set\_backend(self.backend, nreg, nonce)  
File “/usr/share/pyfr/pyfr/util.py”, line 48, in **call**  
return proxylist(x(\*args, \*\*kwargs) for x in self)  
File “/usr/share/pyfr/pyfr/util.py”, line 48, in   
return proxylist(x(\*args, \*\*kwargs) for x in self)  
File “/usr/share/pyfr/pyfr/solvers/navstokes/elements.py”, line 12, in set\_backend  
super().set\_backend(backend, nscalupts, nonce)  
File “/usr/share/pyfr/pyfr/solvers/baseadvecdiff/elements.py”, line 26, in set\_backend  
super().set\_backend(backend, nscal\_upts, nonce)  
File “/usr/share/pyfr/pyfr/solvers/baseadvec/elements.py”, line 26, in set\_backend  
super().set\_backend(backend, nscal\_upts, nonce)  
File “/usr/share/pyfr/pyfr/solvers/base/elements.py”, line 179, in set\_backend  
for i in range(nscal\_upts)]  
File “/usr/share/pyfr/pyfr/solvers/base/elements.py”, line 179, in   
for i in range(nscal\_upts)]  
File “/usr/share/pyfr/pyfr/backends/base/backend.py”, line 20, in newfn  
m = fn(self, \*args, \*\*kwargs)  
File “/usr/share/pyfr/pyfr/backends/base/backend.py”, line 133, in matrix  
return self.matrix\_cls(self, ioshape, initval, extent, aliases, tags)  
File “/usr/share/pyfr/pyfr/backends/base/types.py”, line 110, in **init**  
aliases, tags)  
File “/usr/share/pyfr/pyfr/backends/base/types.py”, line 65, in **init**  
backend.malloc(self, extent)  
File “/usr/share/pyfr/pyfr/backends/base/backend.py”, line 71, in malloc  
obj.onalloc(data, 0)  
File “/usr/share/pyfr/pyfr/backends/openmp/types.py”, line 22, in onalloc  
self.\_set(self.\_initval)  
File “/usr/share/pyfr/pyfr/backends/openmp/types.py”, line 31, in \_set  
self.data[:, :self.ncol] = self.\_pack(ary)  
File “/usr/share/pyfr/pyfr/backends/base/types.py”, line 88, in \_pack  
mode=‘constant’)  
File “/usr/lib/python3/dist-packages/numpy/lib/arraypad.py”, line 1363, in pad  
newmat = narray.copy()  
MemoryError  
msk600@virtual:~/AirfoilPyFR$

Any ideas?

Best regards Stian

fredag 21. februar 2020 10.33.36 UTC+1 skrev Giorgio Giangaspero følgende:

---

<div class="post-metadata">

**Author:** ![turbvs](https://avatars.discourse-cdn.com/v4/letter/t/f14d63/32.png) [@turbvs](https://pyfr.discourse.group/u/turbvs)\
**Post date:** [21 February 2020 12:13 UTC](https://pyfr.discourse.group/t/3d-airfoil/272/4 "2020-02-21T12:13:51Z")

</div>

Hi Stian,

A 1 million elements mesh with polynomial order 3 is not a small problem, of course, depending on your machine.

You might be exceeding the memory limit. Have you had a look at that?

BR,  
Vishal

---

<div class="post-metadata">

**Author:** ![Stian\_Hjorteland](https://avatars.discourse-cdn.com/v4/letter/s/b782af/32.png) [@Stian\_Hjorteland](https://pyfr.discourse.group/u/Stian_Hjorteland)\
**Post date:** [2 March 2020 11:08 UTC](https://pyfr.discourse.group/t/3d-airfoil/272/5 "2020-03-02T11:08:39Z")

</div>

fredag 21. februar 2020 13.14.05 UTC+1 skrev Vishal Saini følgende:
