I want to simulate the coupled equations of photodesorption of Oxygen molecules from palladium surface due to laser irradiation in femtosecond range . These are the equations .
Desorption.pdf
In fipy I substitute the adsorbate temperature T_ads in other equations to solve directly and I got four equations , 1 - electrons temperature, 2 - phonons temperature, 3 - Adsorbate temperature , 4 - Desorption yield.
Here is the first code
then I thought that due to very small change in desorption ,this happens. Problem arises when I want to do parametric variation with respect of neu_pfac. After calculation I see that yield is linearly dependent on neu_pfac but from the equations I can see that it can not be. I could not correct the code because I did not understand where are the errors.
I want to simulate the coupled equations of photodesorption of Oxygen molecules from palladium surface due to laser irradiation in femtosecond range . These are the equations . Desorption.pdf In fipy I substitute the adsorbate temperature T_ads in other equations to solve directly and I got four equations , 1 - electrons temperature, 2 - phonons temperature, 3 - Adsorbate temperature , 4 - Desorption yield. Here is the first code
in this case I see that there is no change in desorption yield. So I put these h= 0.5 - Theta , in the equations and I can see the change.
then I thought that due to very small change in desorption ,this happens. Problem arises when I want to do parametric variation with respect of neu_pfac. After calculation I see that yield is linearly dependent on neu_pfac but from the equations I can see that it can not be. I could not correct the code because I did not understand where are the errors.