7 Aug
2014
7 Aug
'14
8:43 p.m.
Hi Robert,
Me again :) but I think this one is easy... in the previous thread about flux I was computing electric field from voltage in the post processing step. I would like to compute energy density also, from electric field. I believe its something like this, in the post_process hook:
# compute the electric field
electric_field = pb.evaluate('ev_grad.i1.Omega( v )', mode='el_avg')
electric_field *= -1.0
out['electric_field'] = Struct(name='Electric field', mode='cell',
data=electric_field, dofs=None)
# compute the energy density
energy_density = pb.evaluate('dw_volume_dot.i1.Omega(p.sigma, ef, ef)',
ef=electric_field) out['energy_density'] = Struct(name='Energy density', mode='cell', data=energy_density, dofs=None)
However, I get an error when its parsing the variables. How do I specify ef variable properly? I attached the relevant files in case its necessary.
Thank you,
Geoff