CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

solve scalarTransport in subdomain

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2016, 11:00
Default solve scalarTransport in subdomain
  #1
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
dear all,

i want to solve scalarTransport in subdomain during runTime, after the flow fields of the whole domain are obtained by a flow model. is it possible to do it in this way:

Code:
#include "createSubMesh.H"
#include "createSubFields.H"  // the scalar fields are registered on the subMesh

// time marching, after U solved over the whole domain
solve
(
// scalar equations...
);
or

is there any other alternative to implement multiphysics modeling in subdomain? Two threads related to my problem but not fix it: http://www.cfd-online.com/Forums/ope...-equation.html
http://www.cfd-online.com/Forums/ope...subdomain.html

(i work with OF-2.3.0)

thanks for your interest. your suggestion will make my day.

/karekle
cfdopenfoam is offline   Reply With Quote

Old   May 16, 2016, 23:48
Default
  #2
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
greetings to all,

i think i should describe my problem more specifically, with the attached image. i hope this can help make myself clear.

after some searching, it seems that modeling multiphysics on subdomain is still a remained problem. there could be several reasons.

its difficult for one to only work with portion of the mesh to solve matrix equations.
to model physics on subdomain, one might re-create sub-mesh or split mesh into regions, and have to define new boundary conditions on the subdomain.
what makes things more difficult is that one may need to solve multiphysics on subdomain during runTime, which means one can not develop a new solver to handle the scalarTransport calculations after the flow simulation ends.

i am now still struggling on this way and, i think this problem deserves efforts. i will keep updating this thread based on my knowledge and put my solution here if i find one. also, i hope senior OF-users could kindly offer your comments, hints or suggestions here, which will receive highly appreciation!

thank you for your time.

/Karekle
Attached Images
File Type: png subdomain.png (70.0 KB, 56 views)
cfdopenfoam is offline   Reply With Quote

Old   May 17, 2016, 11:07
Default
  #3
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
sorry that I missed sth important here.

the problem, both the flow and transport are transient, which means that we can not use mapfields to map flow fields from the whole domain to the subdomain and start another transport simulation. however, this can obviously be done in a steady case. we may need to do mapfields at every time step during runTime and solve the transport on the subdomain.
cfdopenfoam is offline   Reply With Quote

Old   May 17, 2016, 22:22
Default
  #4
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,205
Rep Power: 31
arjun will become famous soon enougharjun will become famous soon enough
Solving in a subdomain or solving in some regions while not solving in others could be a big design change.
I am not familiar with openform so I can not really tell you how difficult or easy it is to do in openfoam.

If I had to do it in my solver FVUS wilkdkatze, it was as easy as setting up the case. Because this is base of design of FVUS.

In my case I would just create two regions and connect them by GGI (inteface). I would solve flow on two regions and energy on only 1.

If the region where T is solved changing with iterations or time, then FVUS has concept of partial region. Once should be able to define partion region and use this to create solver for Engery.

(Note FVUS at the moment only active on immersed boundary as partion region but other option for partial region creation like from volume of fraction are around the corner).

Good luck.



Quote:
Originally Posted by cfdopenfoam View Post
sorry that I missed sth important here.

the problem, both the flow and transport are transient, which means that we can not use mapfields to map flow fields from the whole domain to the subdomain and start another transport simulation. however, this can obviously be done in a steady case. we may need to do mapfields at every time step during runTime and solve the transport on the subdomain.
arjun is offline   Reply With Quote

Old   May 17, 2016, 22:52
Default
  #5
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
Quote:
Originally Posted by arjun View Post
Solving in a subdomain or solving in some regions while not solving in others could be a big design change.

If I had to do it in my solver FVUS wilkdkatze, it was as easy as setting up the case. Because this is base of design of FVUS.

In my case I would just create two regions and connect them by GGI (inteface). I would solve flow on two regions and energy on only 1.

If the region where T is solved changing with iterations or time, then FVUS has concept of partial region. Once should be able to define partion region and use this to create solver for Engery.

(Note FVUS at the moment only active on immersed boundary as partion region but other option for partial region creation like from volume of fraction are around the corner).

Good luck.
thank you for your hints, Arjun.

I just googled your solver, FVUS wilkdkatze, but got little information on that. if you are happy to tell more, i am very much looking forward to it.
cfdopenfoam is offline   Reply With Quote

Old   May 17, 2016, 23:13
Default
  #6
Senior Member
 
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,205
Rep Power: 31
arjun will become famous soon enougharjun will become famous soon enough
Quote:
Originally Posted by cfdopenfoam View Post
thank you for your hints, Arjun.

I just googled your solver, FVUS wilkdkatze, but got little information on that. if you are happy to tell more, i am very much looking forward to it.

On youtube,

https://www.youtube.com/channel/UCs6...M6JtRa8bLuVr6A

I am in process of releasing it. You could beta test it if you want.
Feedback help me, make it better.
arjun is offline   Reply With Quote

Old   May 17, 2016, 23:20
Default
  #7
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
Quote:
Originally Posted by arjun View Post
On youtube,

https://www.youtube.com/channel/UCs6...M6JtRa8bLuVr6A

I am in process of releasing it. You could beta test it if you want.
Feedback help me, make it better.
thank you. i will have a try.
cfdopenfoam is offline   Reply With Quote

Old   May 19, 2016, 08:16
Default
  #8
Member
 
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10
cfdopenfoam is on a distinguished road
hi all.

A little progress on the proposed problem.

I modified the scalarTransportFoam and developed a very simple solver subregionTransport, following the ideas mentioned above. I also setup a pitzDaily case to test the solver. unfortunately it did not work well. But I think we can work together to fix it. I do not know how to attach the solver and case here. if interested in my problem, plz leave a message. I will send you the files.
cfdopenfoam is offline   Reply With Quote

Old   June 27, 2017, 10:47
Default
  #9
Member
 
Ashish Kumar
Join Date: Jun 2015
Posts: 33
Rep Power: 9
ashish.svm is on a distinguished road
Quote:
Originally Posted by cfdopenfoam View Post
hi all.

A little progress on the proposed problem.

I modified the scalarTransportFoam and developed a very simple solver subregionTransport, following the ideas mentioned above. I also setup a pitzDaily case to test the solver. unfortunately it did not work well. But I think we can work together to fix it. I do not know how to attach the solver and case here. if interested in my problem, plz leave a message. I will send you the files.
Hi,

I am facing a similar problem. I want to solve an equation of a sub-domain of the original mesh. If you have resolved the problem and if you can share the techniques used, then please post here.

Thanks a lot in advance.
ashish.svm is offline   Reply With Quote

Old   March 28, 2020, 11:33
Default any luck?
  #10
Member
 
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 16
r2d2 is on a distinguished road
Hi, I am also interested in a correct setup/solver for solving some scalar transport on just a part of the domain using the global velocity field, as you mentioned in the two messages at the begin of the thread. Did you manage to get it work?
r2d2 is offline   Reply With Quote

Reply

Tags
modeling, multiphysics, openfoam, subdomain

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Solve Flow or VOF simultaneously ? Ramsey FLUENT 1 February 16, 2011 13:16
Regarding the adaptation of scalarTransport Foam as well as general openFOAM logic. HFonten1 OpenFOAM 0 August 10, 2010 13:05
Linearized NS euqations: how to solve them?(problem with Matrix operations..) matteoL OpenFOAM Running, Solving & CFD 0 November 18, 2009 06:58
Subdomain or Interfaces (CHT) sandeep_tu CFX 8 July 14, 2009 11:06
set temperature on subdomain, in transient mode nicolas CFX 0 August 16, 2006 06:42


All times are GMT -4. The time now is 03:50.