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

Passive Scalar Out of Bounds Channel

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By akidess

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2017, 18:22
Default Passive Scalar Out of Bounds Channel
  #1
Senior Member
 
FlyBob91's Avatar
 
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10
FlyBob91 is on a distinguished road
Hello to all,
I'm solving the flow in a very simple geometry: the well known Tjunction channel. It is a microfluids device so the dimensions are of the order of 10^(-6). The flow is completely laminar due to the very low Re number (Re=100). The channel is characterized by two inlets and one outlet. I linked an image to be more clear.
My flow is single phase, water at ambient temperature. Now, after the calculation of the velocity and pressure field with simpleFoam, i run a simulation to calculate the passive scalar (DT=3*10^(-10)). At one inlet i imposed a T=0 while at the other T=1 and then i run scalarTransportFoam. In this manner i can observe how the two flows mix in the outlet channel. The big problem is that my scalar goes out of bound, this means i find value larger 1 and it shouldn't happen. I tried to fix the problem using different schemes for the divergence as suggested in other threads but nothing, the problem is always there.
Is there someone who have experience in this type of problems?

I link also the two cases.
Thanks for help.


Best Regards,
Roberto
Attached Images
File Type: png Schermata del 2017-08-22 23-06-09.png (55.0 KB, 29 views)
Attached Files
File Type: gz Direct.tar.gz (2.2 KB, 2 views)
File Type: gz Scalar.tar.gz (1.8 KB, 6 views)
FlyBob91 is offline   Reply With Quote

Old   August 23, 2017, 03:45
Default
  #2
Senior Member
 
akidess's Avatar
 
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30
akidess will become famous soon enough
Check your mesh.
Tobi likes this.
__________________
*On twitter @akidTwit
*Spend as much time formulating your questions as you expect people to spend on their answer.
akidess is offline   Reply With Quote

Old   August 23, 2017, 05:24
Default
  #3
Senior Member
 
FlyBob91's Avatar
 
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10
FlyBob91 is on a distinguished road
Hi akidess, first thank you for your reply.

Quote:
Originally Posted by akidess View Post
Check your mesh.
What do you mean? Sorry, i don't understand.
Here there is my output after running checkMesh


Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1612+                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1612+
Exec   : checkMesh
Date   : Aug 23 2017
Time   : 10:19:27
Host   : "Nitro"
PID    : 4231
Case   : /home/roby/Scrivania/0.4/Direct
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           306516
    faces:            863850
    internal faces:   810150
    cells:            279000
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     279000
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
    Patch               Faces    Points   Surface topology                  
    walls               45900    46330    ok (non-closed singly connected)  
    inlet1              300      336      ok (non-closed singly connected)  
    inlet2              300      336      ok (non-closed singly connected)  
    movable             6600     6951     ok (non-closed singly connected)  
    outlet              600      651      ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (-0.00015 -0.00165 0) (0.003 0.00165 0.0002)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-6.91724e-18 1.96638e-18 2.57098e-17) OK.
    Max cell openness = 1.29247e-16 OK.
    Max aspect ratio = 1 OK.
    Minimum face area = 1e-10. Maximum face area = 1e-10.  Face area magnitudes OK.
    Min volume = 1e-15. Max volume = 1e-15.  Total volume = 2.79e-10.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 3.4169e-13 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
It seems ok to me. Is there something i forgot?


Best Regards,
Roberto
Attached Files
File Type: gz blockMeshDict.tar.gz (906 Bytes, 1 views)
FlyBob91 is offline   Reply With Quote

Old   August 23, 2017, 05:41
Default
  #4
Senior Member
 
Robert
Join Date: May 2015
Location: Bremen, GER
Posts: 292
Rep Power: 12
RobertHB is on a distinguished road
Hello Roberto,
did you check your fvSchemes settings? Using Gauss limitedLinear 1; for the divergence scheme of the passive scalar usually works fine for me.

Last edited by RobertHB; August 23, 2017 at 06:57.
RobertHB is offline   Reply With Quote

Old   August 23, 2017, 09:51
Default
  #5
Senior Member
 
FlyBob91's Avatar
 
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10
FlyBob91 is on a distinguished road
Quote:
Originally Posted by RobertHB View Post
Hello Roberto,
did you check your fvSchemes settings? Using Gauss limitedLinear 1; for the divergence scheme of the passive scalar usually works fine for me.
Hello, Robert.
I tried the "Gauss limitedLinear01 1" scheme for the divergence but nothing, the solution continues to be out of bounds
FlyBob91 is offline   Reply With Quote

Old   August 25, 2017, 04:31
Default
  #6
Senior Member
 
FlyBob91's Avatar
 
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10
FlyBob91 is on a distinguished road
Little update
I noticed few days ago something strange: of you not copy the phi file in the 0 folder scalarTransportFoam produces a phi file which is different from that produced by simpleFoam. This is very strange. Using the phi file produced by simpleFoam the solution improves but not enough.
Why this difference?


Best Regards,
Roberto
FlyBob91 is offline   Reply With Quote

Old   June 6, 2019, 09:37
Default
  #7
Senior Member
 
Alejandro
Join Date: Jan 2014
Location: Argentina
Posts: 128
Rep Power: 12
ancolli is on a distinguished road
I am facing a similar problem. Did u manage to solve the problem?
ancolli is offline   Reply With Quote

Old   July 17, 2019, 04:42
Default
  #8
Senior Member
 
FlyBob91's Avatar
 
Join Date: Mar 2016
Location: Bergamo
Posts: 157
Rep Power: 10
FlyBob91 is on a distinguished road
Quote:
Originally Posted by ancolli View Post
I am facing a similar problem. Did u manage to solve the problem?
Try with "limitedLinear 1" scheme for the divergence term. It should work.

Best regards.
FlyBob91 is offline   Reply With Quote

Reply


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
problem during mpi in server: expected Scalar, found on line 0 the word 'nan' muth OpenFOAM Running, Solving & CFD 3 August 27, 2018 05:18
Division by zero exception - loop over scalarField Pat84 OpenFOAM Programming & Development 6 February 18, 2017 06:57
Issue symmetryPlane 2.5d extruded airfoil simulation 281419 OpenFOAM Running, Solving & CFD 5 November 28, 2015 14:09
Diverging solution in transonicMRFDyMFoam tsalter OpenFOAM Running, Solving & CFD 30 July 7, 2014 07:20
compressible flow in turbocharger riesotto OpenFOAM 50 May 26, 2014 02:47


All times are GMT -4. The time now is 13:47.