I am trying to implement a chemical solver in an implicit approximate factorization code. I need a chemical source term Jacobian. A derivation has been done for this in my group in the past, however the derivation assumed that the fluxes are written as such:
where rho_1 is the density of species 1, etc.
The code I am using now writes the fluxes in this way:
Would anyone know how to "convert" -so to speak- the Jacobian for the first type of flux definition to a Jacobian for the second type of flux definition?
Thanks very much in advance for your replies
Re: Chemical Jacobians
Assuming Q1 is the old Q vector and Q2 is the new Q vector and W is the sorce term vector, you then already have d_W/d_Q1 and you want to compute d_W/d_Q2
d_W/d_Q2 = d_W/d_Q1 * d_Q1/d_Q2
you have to compute the transformation matrix d_Q1/d_Q2 and post multiply the existing jacobian with it. THat should be all.
|All times are GMT -4. The time now is 21:59.|