|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
Yinghao Zhou
Join Date: Jun 2024
Location: China
Posts: 3
Rep Power: 2 ![]() |
Hello everyone.
I am using the following equation to calculate the electron potential, and k_eff_ele is a constant; fvScalarMatrix Fai_eleEqn ( fvm::ddt(tau_potential, Fai_ele) - fvm::laplacian(k_eff_ele, Fai_ele) == S_ele ); When k_eff_ele is a smaller value of 30315, The calculated residual can drop below 1e-7, but when k_eff_ele is a larger value of 2126933, the calculation is difficult to converge. I don't know how to modify the code to make it converge. Can anyone help me? The fvSchemes is as follows: ddtSchemes { default Euler; } gradSchemes { default Gauss linear; } divSchemes { default bounded Gauss upwind; div(phi_lamda,lamda) Gauss upwind; div(phi,U) Gauss upwind; div(phi,K) Gauss upwind; div(phi,h) Gauss upwind; div(phi,k) Gauss upwind; div(phi,epsilon) Gauss upwind; div(phi,R) Gauss upwind; div(R) Gauss upwind; div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; div(phi,Yi_h) Gauss limitedLinear 1; div(phi,C) Gauss upwind; } laplacianSchemes { default Gauss linear uncorrected; } interpolationSchemes { default linear; } snGradSchemes { default uncorrected; } The fvSolution is as follows: Fai_ele { solver smoothSolver; smoother symGaussSeidel; tolerance 1e-05; relTol 0.0001; } Fai_eleFinal { solver PBiCGStab; preconditioner diagonal; tolerance 1e-7; relTol 0; } |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Error in New Surface reaction model (Having multiple reactions) | surajkvs | OpenFOAM Programming & Development | 2 | May 23, 2023 21:21 |
Solving for an additional species CO in coalChemistryFoam | N. A. | OpenFOAM Programming & Development | 3 | February 18, 2019 05:58 |
access to member function variable | Rahand | OpenFOAM Programming & Development | 3 | March 29, 2016 15:46 |
NOx calculations slow to converge - faster method? | Antony92from | FLUENT | 0 | January 16, 2016 05:40 |
dieselFoam problem!! trying to introduce a new heat transfer model | vivek070176 | OpenFOAM Programming & Development | 10 | December 23, 2014 23:48 |