|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Ftab
Join Date: Sep 2011
Posts: 87
Rep Power: 14 ![]() |
Dear all,
As I am trying to solve an ODE for additional variable which looks like this: dX/dt=aX+b I tried with diffusion equation with Zero diffusivity, but Kin_Diff should be between 0 and 1. What should I do? Put Kin_Diff equal to 1e-100? Is there a better way? By the way, the source term was modeled as a subdomain and source. Was it correct? Thanks |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,478
Rep Power: 140 ![]() ![]() ![]() ![]() |
That should work. Yes, you can turn the diffusion off by setting a very low diffusivity.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
Member
Ftab
Join Date: Sep 2011
Posts: 87
Rep Power: 14 ![]() |
Quote:
Thanks a lot Glenn! The simulation gave this error, but ran and converged: Bounds error detected --------------------- Variable: X.Dynamic Diffusivity Locale : wall Bounds error detected --------------------- Variable: Y.Dynamic Diffusivity Locale : wall And this was for D=1e-40 Last edited by ftab; November 15, 2018 at 14:45. |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,478
Rep Power: 140 ![]() ![]() ![]() ![]() |
Looks like 1e-40 is too low and triggers a bounds warning. Try 1e-10. And if you want to be sure do a sensitivity check to ensure that this value is low enough that it is effectively zero (that is, run 1e-9 and compare the results. If no significant difference then it is effectively zero).
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
![]() |
![]() |
![]() |
![]() |
#5 | |
Member
Ftab
Join Date: Sep 2011
Posts: 87
Rep Power: 14 ![]() |
Quote:
Update: This strategy works only for D up to 1e-15, as the truncation error of 1e-16 returns zero and bound error for the Diff coefficient |
||
![]() |
![]() |
![]() |
![]() |
#6 |
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,478
Rep Power: 140 ![]() ![]() ![]() ![]() |
A minor point, but I don't think this is a truncation error. I think simply the bounds of that variable are set to 1e-15 and if you go below that it triggers the bounds error.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum. |
|
![]() |
![]() |
![]() |
Tags |
additional variables, source terms |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Fire and smoke modeling | CFDST | CFX | 84 | April 30, 2019 07:01 |
[Other] OpenFOAM Installation for navalFoam | sachinlb | OpenFOAM Community Contributions | 22 | July 28, 2017 05:26 |
[OpenFOAM.org] Error creating ParaView-4.1.0 OpenFOAM 2.3.0 | tlcoons | OpenFOAM Installation | 13 | April 20, 2016 17:34 |
[swak4Foam] funkySetFields compilation error | tayo | OpenFOAM Community Contributions | 39 | December 3, 2012 05:18 |
DxFoam reader update | hjasak | OpenFOAM Post-Processing | 69 | April 24, 2008 01:24 |