OF 15 FLUENT Mesh with interfaces

hjasak · March 12, 2009, 12:49

As I said: if you think it is legal that in your geometry you have uncovered faces on the GGI interface, please use bridging: it will do what is necessary.

If you think the interface should be completely covered, set bridging to off. If you get the message, there is something wrong with your setup and you need to investigate by hand.

Clear?

Hrv

bastil · March 12, 2009, 13:05

Ok I see, thanks. In this case it is legal so I switched and I get:

Attempt to cast type wall to type lduInterface
From function refCast<to>(From&)

Seems like there is something wrong with my bc?

Regards

bastil · March 12, 2009, 13:26

Great, works. Thanks Hrv for your help once more. One ggi was still of type "wall".

Regards.

hjasak · March 12, 2009, 13:34

Happens to me as well - reading the error message (however obscure) will usually help.

Enjoy,

Hrv

bastil · March 12, 2009, 16:16

Hrv,

finally I want to guide you to my add-on:

http://www.cfd-online.com/OpenFOAM_Discus/messages/126/11626.html?1236852411#POST3342 1

It helps do identify the problems.

Regards

bastil · March 13, 2009, 04:52

Hrv, all:

I still have troubles in runnig it on a cluster. On my workstation it runs on one cpu as well as on two, but on the cluster I get floating point errors.
I guess this has to do with parallelisation, I am unsure about how to do it. My first guess for decomposeParDict was:

preservePatches (patch1 patch1_shadow);
preservePatches (patch2 patch2_shadow);
preservePatches (patch3 patch3_shadow);

But obviously this is not the right way decause decomposePar complaines. So what I do now is:

preservePatches (patch1 patch1_shadow
patch2 patch2_shadow patch3 patch3_shadow);

Feedback now looks as expected but I get floating point exceptions. And ideas?

Regards

hjasak · March 13, 2009, 05:07

Well, unless you finished off the parallelisation while I wasn't looking, this should not run (ie. it is dependent on the decomposition). There is still 1-2 bits to implement - you are basically trying to run in parallel without parallel action on GGI.

Wait for it a bit more, I working on this.

Hrv

bastil · March 13, 2009, 05:41

Ok so just choosing preserve patches does not make sure it will run? And if it runes I had luck, right?

Thanks once more and regards.

BastiL

jens_klostermann · March 13, 2009, 05:44

Hi BastiL,

Your problem sounds to me like you have cyclics (cyclicGgi) in more than 1 direction, maybe even in all directions?

One "quick and dirty" solution for now is using preserveFaceZones (instead of preservePatches). preserveFaceZones should be used in that way that there exist a thin cell layer at the cyclics. (via faceSet and setsToZones).

Jens

bastil · March 13, 2009, 05:51

Jens,

I don't have cyclics just "traditional" non-conformal interfaces. However, your waorkaround may help me. Can you please give me some details about how to create these FaceZones? Thanks.

BastiL

bastil · March 13, 2009, 08:56

Ok so what I tried is to create faceZones using setSets. I put patch1 and patch1_shadow into one face zone. So I get 5 face zones for my 10 patches that form 5 ggis. However, I get the same behaviour as before. How do you handle this, Jens?
Seems like there are some other things going on regarding parallelisation.

Regards BastiL

kev4573 · March 31, 2009, 17:36

Bastil, did you get your decomposition to preserve both sides of the entire GGI interface?

Kevin

March 12, 2009, 12:49	As I said: if you think it is	#21
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,907 Rep Power: 33	As I said: if you think it is legal that in your geometry you have uncovered faces on the GGI interface, please use bridging: it will do what is necessary. If you think the interface should be completely covered, set bridging to off. If you get the message, there is something wrong with your setup and you need to investigate by hand. Clear? Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

March 12, 2009, 13:05	Ok I see, thanks. In this case	#22
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Ok I see, thanks. In this case it is legal so I switched and I get: Attempt to cast type wall to type lduInterface From function refCast<to>(From&) Seems like there is something wrong with my bc? Regards

March 12, 2009, 13:26	Great, works. Thanks Hrv for y	#23
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Great, works. Thanks Hrv for your help once more. One ggi was still of type "wall". Regards.

March 12, 2009, 13:34	Happens to me as well - readin	#24
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,907 Rep Power: 33	Happens to me as well - reading the error message (however obscure) will usually help. Enjoy, Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

March 12, 2009, 16:16	Hrv, finally I want to guid	#25
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Hrv, finally I want to guide you to my add-on: http://www.cfd-online.com/OpenFOAM_Discus/messages/126/11626.html?1236852411#POST3342 1 It helps do identify the problems. Regards

March 13, 2009, 04:52	Hrv, all: I still have trou	#26
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Hrv, all: I still have troubles in runnig it on a cluster. On my workstation it runs on one cpu as well as on two, but on the cluster I get floating point errors. I guess this has to do with parallelisation, I am unsure about how to do it. My first guess for decomposeParDict was: preservePatches (patch1 patch1_shadow); preservePatches (patch2 patch2_shadow); preservePatches (patch3 patch3_shadow); But obviously this is not the right way decause decomposePar complaines. So what I do now is: preservePatches (patch1 patch1_shadow patch2 patch2_shadow patch3 patch3_shadow); Feedback now looks as expected but I get floating point exceptions. And ideas? Regards

March 13, 2009, 05:07	Well, unless you finished off	#27
hjasak Senior Member Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,907 Rep Power: 33	Well, unless you finished off the parallelisation while I wasn't looking, this should not run (ie. it is dependent on the decomposition). There is still 1-2 bits to implement - you are basically trying to run in parallel without parallel action on GGI. Wait for it a bit more, I working on this. Hrv __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

March 13, 2009, 05:41	Ok so just choosing preserve p	#28
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Ok so just choosing preserve patches does not make sure it will run? And if it runes I had luck, right? Thanks once more and regards. BastiL

March 13, 2009, 05:44	Hi BastiL, Your problem sou	#29
jens_klostermann Senior Member Jens Klostermann Join Date: Mar 2009 Posts: 117 Rep Power: 17	Hi BastiL, Your problem sounds to me like you have cyclics (cyclicGgi) in more than 1 direction, maybe even in all directions? One "quick and dirty" solution for now is using preserveFaceZones (instead of preservePatches). preserveFaceZones should be used in that way that there exist a thin cell layer at the cyclics. (via faceSet and setsToZones). Jens

March 13, 2009, 05:51	Jens, I don't have cyclics	#30
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Jens, I don't have cyclics just "traditional" non-conformal interfaces. However, your waorkaround may help me. Can you please give me some details about how to create these FaceZones? Thanks. BastiL

March 13, 2009, 08:56	Ok so what I tried is to creat	#31
bastil Senior Member BastiL Join Date: Mar 2009 Posts: 530 Rep Power: 20	Ok so what I tried is to create faceZones using setSets. I put patch1 and patch1_shadow into one face zone. So I get 5 face zones for my 10 patches that form 5 ggis. However, I get the same behaviour as before. How do you handle this, Jens? Seems like there are some other things going on regarding parallelisation. Regards BastiL

March 31, 2009, 17:36		#32
kev4573 Senior Member Kevin Smith Join Date: Mar 2009 Posts: 104 Rep Power: 17	Bastil, did you get your decomposition to preserve both sides of the entire GGI interface? Kevin

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Partioning of Mesh having nonconformal interfaces	Saji	FLUENT	0	July 15, 2006 05:58
Grid interfaces in Fluent	Omar Qazi	FLUENT	1	January 28, 2006 19:55
Fluent - discontinuous grid interfaces	Simon Archen	FLUENT	2	November 16, 2001 10:52
Sliding mesh interfaces	Karl Kevala	FLUENT	1	February 20, 2001 05:31
Periodic grid-interfaces in FLUENT 5	tom	Main CFD Forum	8	May 13, 1999 21:27