Finding computeQpOffDiagJacobian for sink and recombination #27660
-
I am trying to add sink and recombination to the Multi-phase poly-crystalline Allen-Cahn Kernel. The top equation is the equation already used by MOOSE. The second equation is from: https://www.sciencedirect.com/science/article/pii/S0927025622004839#! and it adds vacancy production, vacancy-interstial recombination, and sink terms to the grand potential equation in terms of chemical potential. From SuspectabilityTimeDirivative, matdiffusion, maskedbodyforce(deals with h_m*s_v^o term), and coupledswitchingtimedirivative I was able to recreate the equation from image 1:
I am creating a new kernal that adds the missing term:
Here is the residual in moose i created:
I copied below my .c file, I am struggling to find documentation on finding my off diagonal Jacobian equation. I have been following the steps provided by the INL: https://mooseframework.inl.gov/getting_started/examples_and_tutorials/tutorial01_app_development/step04_weak_form.html#procedure p.s. Am I missing an override anywhere?
|
Beta Was this translation helpful? Give feedback.
Replies: 1 comment 5 replies
-
Hello
what is n here? I see _qp (which seems correct) in your kernel For the Jacobian:
For the off-diagonal
All of these should be override. The nice thing about override is that if you made a mistake in the signature (for example the return type), it wont compile thanks to override
NB: do not forget to check your Jacobian with the Jacobian testing tools in MOOSE. See JacobianTest on the website, or troubleshooting failed solves |
Beta Was this translation helpful? Give feedback.
Hello
It s missing this header : JvarMapKernelInterface.h
Does it exist in your version of moose?
Can you find it in either the framework or the modules headers?
Guillaume