I am trying to code it and most of it is done for now.
Here is my take, during the initialization of model, you first solve the poisson equation for this variable (Ue as you called it in other thread). After it is solved the field comes to be smooth. dUedn is no longer as sharp as it is in the start for cplus and cminus. So basically in my model init cplus, cminus from user. using cplus and cminus, solve poisson equation for Ue to initialize it. Iterations 1. Solve Ue poisson equation 2. Solve Cplus and Cminus transport. This is my rough outline. |
|
All times are GMT -4. The time now is 16:37. |