[exw599]

Hello everyone,

I was using porousSimpleFoam to solve the fluid velocity field of the flow passing a U-shape porous bed. I also want to take the porous subdomain as a solid chunk and solve the stress field using solidDisplacementFoam.

I know you could write additional equations in your own solver to couple different physics, but I just can not find a way to ONLY solve the additional equations in PART of the whole domain (since porous domain is only a fraction of the entire domain).

Any suggestions and helps would be appreciated.

[chriss85]

You should take a look at how multi-region solvers are written. In short, you load different meshes and their fields and solve the equations separately on them. You might need boundary conditions to couple the values between the regions.

[akidess]

A simple alternative is to multiply your equation system with an indicator field which is 1 only where you want the equations to be solved.

[exw599]

Thank you both for your replies, I have checked the multi region solvers but it seems to decompose the whole domain and solve different equations on each subdomain. Maybe a more practical approach for me is to output the results from a sub-region and make it serve as the initial condition for another solver. I am searching for the method to output VectorField only from a sub-region.

[cfdopenfoam]

Originally Posted by exw599

Thank you both for your replies, I have checked the multi region solvers but it seems to decompose the whole domain and solve different equations on each subdomain. Maybe a more practical approach for me is to output the results from a sub-region and make it serve as the initial condition for another solver. I am searching for the method to output VectorField only from a sub-region.

hi, exw599, have you got a good way to do this? i now have exactly the same situation.

[cfdopenfoam]

Originally Posted by akidess

A simple alternative is to multiply your equation system with an indicator field which is 1 only where you want the equations to be solved.

hi, Anton. thank you for your hints here. it's a good approach, but i am wondering how the boundary conditions could be imposed because they only need to be specified at the boundaries of the subdomain.

are there any simple cases we can refer to?

thanks very much!

[exw599]

Hi in the end I simply solve two different equations on using two different solvers with the results from first solver serving as initial condition for the second solver. You can use mapfield function to map field value between different meshes as long as they have intersections.

[cfdopenfoam]

Originally Posted by exw599

Hi in the end I simply solve two different equations on using two different solvers with the results from first solver serving as initial condition for the second solver. You can use mapfield function to map field value between different meshes as long as they have intersections.

thank you for your enlightment. this is very helpful and can be a very nice solution for steady case. unfortunately i work on an unsteady case and i do very long time simulation. it's nearly impossible to write required fields at each time step and then read as solved parameters for the second solver.

[exw599]

Originally Posted by cfdopenfoam

thank you for your enlightment. this is very helpful and can be a very nice solution for steady case. unfortunately i work on an unsteady case and i do very long time simulation. it's nearly impossible to write required fields at each time step and then read as solved parameters for the second solver.

Hi that depends on the ratio of the relaxation times between the two coupling dynamics, I suppose if one is significantly smaller than the other you could safely ignore the dynamics with smaller relaxation time and solve it using steady state solver instead and solve full dynamics for the other. Hope this helps.