|
[Sponsors] |
[swak4Foam] reactingMultiPhaseEulerFoam problems with groovyBC |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
July 1, 2016, 11:50 |
reactingMultiPhaseEulerFoam problems with groovyBC
|
#1 |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 11 |
Dear All,
Would it be possible please to help me with the below issue. Maybe some of you already stumbled in this problem and solved it. I am dealing with simulations of fluidized beds of particles. In particular I have to perform simulations where the air flow from the bed distribution grid is variable with time and from point to point of the grid. With twoPhaseEulerFoam version 2.4.0 I was able to use groovyBC to specify the time and position dependence of the velocity at the distribution grid. Also I was able to compute the minimum fluidization velocity for comparison with experimental data by increasing the velocity at the distribution grid over time. To do this I used as a boundary condition in the U.air file: Code:
inlet { aliases { alphaPart alpha.particles; alphaAir alpha.air; } type groovyBC; variables ( "t=time();" "u0=0.3000;" "u1=0.6000;" "t0=5.0;" "t1=7.0;" "uf=max(u0,min((u1-u0)/(t1-t0)*(t-t0)+u0,u1));" "inVel=uf/min(1.0-alphaPart,0.38);" ); valueExpression "vector(0,inVel,0)"; value $internalField; } Code:
"inVel=uf/min(1.0-alphaPart,0.38);" I tried to do the same with reactingMultiphaseEulerFoam from release 3.0.1 of OF since I would like to expand my simulation scope including multiple particles phases with different diameter and properties and mass exchange between the phases. Unfortunately I run into a problem that I was not able to overcome. The problem is readily apparent even in the simplest of cases. At first I thought that I made some mistakes in the case setup for the new solver but in the end I was able to reproduce it using the fluidized bed tutorial from reactingTwophaseEulerFoam with small modification to adapt it to reactingMultiphaseEulerFoam. The case is quite simple. Setting a surface velocity of 0.1 m/s using InterstitialInletVelocity Code:
inlet { type interstitialInletVelocity; inletVelocity uniform (0 0.1 0); alpha alpha.air; value $internalField; } If I try the same using the groovyBC boundary conditions, that worked for twoPhaseEulerFoam version 2.4.0, I basically obtain plugs of particles that are shot out of the computational domain and the simulation crash with a segfault after 0.15 s. See attachments test-02* This strange behavior happens also if I specify a surface velocity of 0 m/s as shown in the attachments test-03* Code:
inlet { aliases { alphaPart alpha.particles; alphaAir alpha.air; } type groovyBC; variables ( "t=time();" "u0=0.0000;" "u1=0.0000;" "t0=0.0;" "t1=2.0;" "uf=max(u0,min((u1-u0)/(t1-t0)*(t-t0)+u0,u1));" "inVel=uf/min(1.0-alphaPart,0.38);" ); valueExpression "vector(0,inVel,0)"; value $internalField; } What am I missing? The groovyBC worked fine with 2.4. Have I to modify my boundary condition somehow to meet some requirement of reactingMultiphaseEulerFoam 3.0.1? It is possible at all to use groovyBC with vertion 3.0.1 to do the same thing as in twoPhaseEulerFoam 2.4.0 Do I have to write my own BC to implement time a space variation of velocity at the inlet grid? Any help from the community is appreciated Best Regards Luca P.S. I have attached the test file I have derived from the original tutorial in tutorials/multiphase/reactingTwoPhaseEulerFoam/RAS/fluidisedBed/ To meet the file size requirements I have removed the polymesh directory, but the mesh is the same as in the tutorial. The pictures are in the following posts. |
|
July 1, 2016, 11:52 |
|
#2 |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 11 |
Test 01 images
|
|
July 1, 2016, 11:52 |
|
#3 |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 11 |
Test 02 Images
|
|
July 1, 2016, 11:53 |
|
#4 |
New Member
Join Date: Jul 2015
Posts: 7
Rep Power: 11 |
Test 03 Images
|
|
December 13, 2016, 07:46 |
|
#5 |
New Member
Sebastien L
Join Date: Mar 2016
Posts: 1
Rep Power: 0 |
Hello, just seen your message. Even if i cannot really help you, i'm a bit surprised because you say that your are modeling particles, and you are using a laminar phase with no viscosity for your model, not a granular phase.
In fact the reactingmultiphaseeulerfoam is not equivalent to the reactingtwophaseeulerfoam because there is no granular therory included into it. It treats the phases a bit differently, with sharpning on interphase between fluids but keeping euler-euler modeling of phases inside the fluid with individual transport equation for phases but common pressure. |
|
Tags |
boundary condition, groovybc, time, twophaseeulerfoam, variable |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[swak4Foam] groovyBC and funkySetFields married and got a kid named swak4Foam | gschaider | OpenFOAM Community Contributions | 169 | August 10, 2023 10:01 |
[swak4Foam] having problems configuring parabolic inlet using groovyBC | asharma | OpenFOAM Community Contributions | 18 | February 12, 2018 09:40 |
[swak4Foam] GroovyBC for 2D wave flume! | Hisham | OpenFOAM Community Contributions | 13 | January 20, 2012 06:04 |
[swak4Foam] Wall heat transfer using groovyBC (XiFoam solver) | usergk | OpenFOAM Community Contributions | 7 | February 4, 2011 14:36 |
[swak4Foam] groovyBC: problems compiling: "flex: not found" and "undefined reference to ..." | sega | OpenFOAM Community Contributions | 12 | February 17, 2010 10:30 |