A coupled problem with different equations in different geometries

ywem · September 6, 2018, 13:23

Hi,
I am solving a problem consisting of two coupled equaitons. These two equations are solved in two adjacent geometries, respectivley. Coupling takes place at the interface between them and the interface moves according to the solutions of the equations. Some kind Foamers advised using dynamic mesh and it can deal with one equation in one geometry. But how can I deal with the case when two different equations concern on different geometries.
Dose anyone has some ideas about the method to solve this problem using OpenFOAM?
Many thanks.

mAlletto · September 7, 2018, 04:29

The are also dynamic mesh and overset mesh technologies available for the volume of fluid method. Maybe this helps

simrego · September 7, 2018, 05:48

Just some ideas, maybe helps, maybe just ignore them:

1,
What if you use the interFoam for a base, but use the alpha to blend the equations? I'm not sure if it is possible, but I think it is. Something like:
eqn1 = "dummy stuff";
eqn2 = "dummy stuff";
Solve ( alpha*eqn1 + (1-alpha)*eqn2 == 0 );
The problem with this is to figure out the equation which will change the alpha field, and maybe at the interface there will be some weird stuff...

2,
Or you can check the chtMultiRegionFoam solver. This solver is for two types of regions (solids, and fluids). But you can rewrite it as you wish. You can define the equations for each type of regions, where each region has it's own geometry (mesh). The communication between them is done with a coupled boundary condition, but maybe you will have to write your own BC for your problem, and you can also add the dynamicMesh functionality to it.

ywem · September 7, 2018, 11:30

Dear Alletto and simrego,
Thank you very much for both of the suggestions. Hopefully they will work.
Best wishes,

September 6, 2018, 13:23	A coupled problem with different equations in different geometries	#1
ywem Member Wenming Yang Join Date: Jun 2018 Posts: 42 Rep Power: 8	Hi, I am solving a problem consisting of two coupled equaitons. These two equations are solved in two adjacent geometries, respectivley. Coupling takes place at the interface between them and the interface moves according to the solutions of the equations. Some kind Foamers advised using dynamic mesh and it can deal with one equation in one geometry. But how can I deal with the case when two different equations concern on different geometries. Dose anyone has some ideas about the method to solve this problem using OpenFOAM? Many thanks.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[ANSYS Meshing] Help with element size	sandri_92	ANSYS Meshing & Geometry	14	November 14, 2018 08:54
[Commercial meshers] Problem in import mesh with two coupled region from gambit	rapierrz	OpenFOAM Meshing & Mesh Conversion	2	October 30, 2017 13:20
three coupled equations	ptub	Main CFD Forum	5	February 21, 2016 13:33
coupled equations: prediction correction method	zaynab	Main CFD Forum	0	November 12, 2015 06:49
problem with gradient as source term in coupled equations	Cyp	OpenFOAM Programming & Development	0	August 16, 2010 14:12

September 7, 2018, 04:29		#2
mAlletto Senior Member Michael Alletto Join Date: Jun 2018 Location: Bremen Posts: 616 Rep Power: 16	The are also dynamic mesh and overset mesh technologies available for the volume of fluid method. Maybe this helps

September 7, 2018, 05:48		#3
simrego Senior Member anonymous Join Date: Jan 2016 Posts: 416 Rep Power: 14	Just some ideas, maybe helps, maybe just ignore them: 1, What if you use the interFoam for a base, but use the alpha to blend the equations? I'm not sure if it is possible, but I think it is. Something like: eqn1 = "dummy stuff"; eqn2 = "dummy stuff"; Solve ( alphaeqn1 + (1-alpha)eqn2 == 0 ); The problem with this is to figure out the equation which will change the alpha field, and maybe at the interface there will be some weird stuff... 2, Or you can check the chtMultiRegionFoam solver. This solver is for two types of regions (solids, and fluids). But you can rewrite it as you wish. You can define the equations for each type of regions, where each region has it's own geometry (mesh). The communication between them is done with a coupled boundary condition, but maybe you will have to write your own BC for your problem, and you can also add the dynamicMesh functionality to it.

September 7, 2018, 11:30		#4
ywem Member Wenming Yang Join Date: Jun 2018 Posts: 42 Rep Power: 8	Dear Alletto and simrego, Thank you very much for both of the suggestions. Hopefully they will work. Best wishes,