
[Sponsors] 
May 16, 2016, 11:00 
solve scalarTransport in subdomain

#1 
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
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... ); …… is there any other alternative to implement multiphysics modeling in subdomain? Two threads related to my problem but not fix it: http://www.cfdonline.com/Forums/ope...equation.html http://www.cfdonline.com/Forums/ope...subdomain.html (i work with OF2.3.0) thanks for your interest. your suggestion will make my day. /karekle 

May 16, 2016, 23:48 

#2 
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
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 recreate submesh 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 OFusers could kindly offer your comments, hints or suggestions here, which will receive highly appreciation! thank you for your time. /Karekle 

May 17, 2016, 11:07 

#3 
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
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. 

May 17, 2016, 22:22 

#4  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,205
Rep Power: 31 
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:


May 17, 2016, 22:52 

#5  
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
Quote:
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. 

May 17, 2016, 23:13 

#6  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 1,205
Rep Power: 31 
Quote:
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. 

May 17, 2016, 23:20 

#7  
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
Quote:


May 19, 2016, 08:16 

#8 
Member
Karelke Yu
Join Date: Dec 2014
Posts: 96
Rep Power: 10 
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. 

June 27, 2017, 10:47 

#9  
Member
Ashish Kumar
Join Date: Jun 2015
Posts: 33
Rep Power: 9 
Quote:
I am facing a similar problem. I want to solve an equation of a subdomain 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. 

March 28, 2020, 11:33 
any luck?

#10 
Member
Radu Mustata
Join Date: Mar 2009
Location: Zaragoza, Spain
Posts: 99
Rep Power: 16 
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?


Tags 
modeling, multiphysics, openfoam, subdomain 
Thread Tools  Search this Thread 
Display Modes  


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 