Hi all, I am currently implementing a cell problem in DuMuX and am using an adapted 1p model/Darcy Problem with CCTPFA discretization to do so. After obtaining the solution vector, I need to obtain the gradients of the solution which I will then use to calculate my upscaled conductivities through integration.
In the exchange with Dimitry Pavlov in the mailing list archive (see 2020q2 "How to implement MaterialLawParams that depend on pressure gradient"), Timo suggested using the Jacobian to obtain the gradients (which I could easily make available in my case, as the step is done AFTER solving the cell problem). However I am slightly confused - isn't the Jacobian the partial derivative of the residual with respect to the solution (instead of the partial residual of the solution with respect to the dimension)? At least this is what the handbook seems to imply. Thanks for any help, Helena
_______________________________________________ DuMux mailing list DuMux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
