CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > SU2

MULTIZONE with species transport bombing

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 10, 2022, 09:04
Default MULTIZONE with species transport bombing
  #1
New Member
 
Robert Tramel
Join Date: Nov 2022
Posts: 3
Rep Power: 2
DoctorRobert is on a distinguished road
I am trying to run an analysis that involves multizone fluid-fluid flow with species. To create a simple test case I took the .su2 mesh file for the incompressible flow/species transport tutorial and duplicated and shifted one mesh to create a simple multizone case. The problem works beautifully without species transport but when I add that in the code bombs in the setup phase and I get the following error message. Can anyone help me out here?

Error in "virtual void CMultizoneOutput::SetMultizoneHistoryOutputFields( const COutput* const*, const CConfig* const*)":
-------------------------------------------------------------------------
Per Surface output group unknown: SPECIES_COEFF_SURF
------------------------------ Error Exit -------------------------------
DoctorRobert is offline   Reply With Quote

Old   November 10, 2022, 12:28
Default
  #2
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 269
Rep Power: 15
bigfootedrockmidget is on a distinguished road
Hi,


We have not implemented/tested this keyword for multizone yet.

You could first try by removing the SPECIES_COEFF_SURF keyword from the config file. Do you need to optimize a species related objective?
bigfootedrockmidget is offline   Reply With Quote

Old   November 10, 2022, 17:37
Default
  #3
New Member
 
Robert Tramel
Join Date: Nov 2022
Posts: 3
Rep Power: 2
DoctorRobert is on a distinguished road
I never included the term in the cfg files because I didn't know it was needed. Should I be defining it? I am not trying to optimize anything species related. I'm trying multizone because I want to use a body force in certain zones as a crude but simple way to mimic a porous media.
DoctorRobert is offline   Reply With Quote

Old   November 11, 2022, 02:46
Default
  #4
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 269
Rep Power: 15
bigfootedrockmidget is on a distinguished road
I did a quick test. If I run the CHT case:
https://su2code.github.io/tutorials/Static_CHT/
and add species transport, I do not get this message. I just add this to the flow_cylinder.cfg fluid setup:
Code:
% -------------------- SCALAR TRANSPORT ---------------------------------------%
%
KIND_SCALAR_MODEL= SPECIES_TRANSPORT
DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
DIFFUSIVITY_CONSTANT= 0.001
%
CONV_NUM_METHOD_SPECIES= SCALAR_UPWIND
MUSCL_SPECIES= NO
SLOPE_LIMITER_SPECIES = NONE
%
TIME_DISCRE_SPECIES= EULER_IMPLICIT
%
SPECIES_INIT= 1.0, 0.0
SPECIES_CLIPPING= YES
SPECIES_CLIPPING_MIN= 0.0, 0.0
SPECIES_CLIPPING_MAX= 1.0, 1.0

Maybe you can check what the difference in setup is.


Can you share your setup somewhere? If you can create a bug report on github with a case that reproduces this error, we can look into it.

https://github.com/su2code/SU2/issues
bigfootedrockmidget is offline   Reply With Quote

Old   November 11, 2022, 20:11
Default
  #5
New Member
 
Robert Tramel
Join Date: Nov 2022
Posts: 3
Rep Power: 2
DoctorRobert is on a distinguished road
The problem comes when you specify MARKER_ANALYZE to compute mass flux . However when you remove that there is no communication of species information across the interface. as can be seen in the following images of velocity, species_0, and species_1
Attached Images
File Type: jpg velocity.jpg (28.7 KB, 8 views)
File Type: jpg species_0.jpg (27.2 KB, 6 views)
File Type: jpg species_1.jpg (25.3 KB, 6 views)
DoctorRobert is offline   Reply With Quote

Old   November 16, 2022, 08:33
Default
  #6
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 269
Rep Power: 15
bigfootedrockmidget is on a distinguished road
Hi,



OK, there were actually 2 problems here. One is that the species mass fractions were not transported across the zone interface. The other is that for MARKER_ANALYZE, we did not add species information, but later in the code we assumed that it was present.



Both are now fixed in the feature branch, but it will be a while before it ends up in the master version. In the mean time, you could clone the SU2 git repository and then do a git checkout feature_species_multizone and compile the fixed branch yourself.
Thanks for spotting the problem!
bigfootedrockmidget is offline   Reply With Quote

Old   November 19, 2022, 18:19
Default
  #7
Senior Member
 
bigfoot
Join Date: Dec 2011
Location: Netherlands
Posts: 269
Rep Power: 15
bigfootedrockmidget is on a distinguished road
By the way, in the develop branch, species transport has been extended with species dependent properties. the keyword PASSIVE_SCALAR has therefore also been replaced by SPECIES_TRANSPORT
If you clone the code from github, you can find some examples in the Testcases/species_transport directory
bigfootedrockmidget is offline   Reply With Quote

Reply

Tags
incompressible, multizone, species transport

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
Modelling Biomass Combustion via Species Transport Racheal FLUENT 39 January 8, 2022 07:42
UDF for Sorption in Species Transport through Porous Medium Kosch Fluent UDF and Scheme Programming 6 October 21, 2020 17:46
Direct numerical simulation of species transport equation with phase change Pmaroul Main CFD Forum 2 October 12, 2018 16:02
A question about species box of species transport model tjushang FLUENT 2 June 25, 2014 07:56
user defined scalar vs species transport gvj_mech FLUENT 0 March 27, 2014 14:33


All times are GMT -4. The time now is 21:31.