|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
exw599
Join Date: Jan 2016
Posts: 6
Rep Power: 10 ![]() |
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. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Join Date: Oct 2013
Posts: 397
Rep Power: 18 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 ![]() |
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.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
exw599
Join Date: Jan 2016
Posts: 6
Rep Power: 10 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#5 | |
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 11 ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#6 | |
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 11 ![]() |
Quote:
are there any simple cases we can refer to? thanks very much! |
||
![]() |
![]() |
![]() |
![]() |
#7 |
New Member
exw599
Join Date: Jan 2016
Posts: 6
Rep Power: 10 ![]() |
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.
|
|
![]() |
![]() |
![]() |
![]() |
#8 | |
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 11 ![]() |
Quote:
![]() |
||
![]() |
![]() |
![]() |
![]() |
#9 | |
New Member
exw599
Join Date: Jan 2016
Posts: 6
Rep Power: 10 ![]() |
Quote:
|
||
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
simpleFoam error - "Floating point exception" | mbcx4jc2 | OpenFOAM Running, Solving & CFD | 12 | August 4, 2015 02:20 |
Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 06:20 |
Unstabil Simulation with chtMultiRegionFoam | mbay101 | OpenFOAM Running, Solving & CFD | 13 | December 28, 2013 13:12 |
calculation stops after few time steps | sivakumar | OpenFOAM Running, Solving & CFD | 7 | March 17, 2013 06:37 |
Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 18:07 |