On 01/11/2017 05:14 PM, David Jessop wrote:
Hi Robert,
I get how 'ev_integrate_mat' etc. can be used to evaluate the material values but I want specifically to output these values, say as part of the solution vector. In material_nonlinearity.py there is a post-processing hook that does this, but how would this work in interactive mode?
It is the same. Note that pb.solve() returns a State instance, not a simple vector. So:
state = pb.solve() out = state.create_output_dict()
# Add to out as in post_process()
pb.save_state('foo.vtk', out=out)
Regarding the second part of my question, I want to call the value of "viscosity" on the boundaries within another function (my "flux" function depends on viscosity).
You can call get_viscosity() (talking about nonLinearPoisson.py) directly, no?
r.
Cheers, David
On Wednesday, 11 January 2017 11:35:25 UTC+1, Robert Cimrman wrote:
Hi David,
On 01/11/2017 09:32 AM, David Jessop wrote:
Hi Robert.
How can material properties, such as "viscosity" in my problem be saved to file (say vtk format) for inspection? I'd like to be able to do this at two points in my programme - following the call to set the initial guess for the velocity field and once the solution has converged.
See post_process() in [1]. The term ev_integrate_mat used there allows also obtaining the material evaluated in the quadrature points of all elements, in case you need that.
[1]
http://sfepy.org/doc-devel/examples/linear_elasticity/material_nonlinearity....
Also, can material parameters be evaluated at arbitrary regions, such as the edges of the domain?
You mean manually? What exactly do you need?
Each term evaluates a material where needed - volume integrals in volume quadrature points, surface integrals in the surface (edge or face) quadrature points.