Source Term for Joule Heating

TomStarke · August 9, 2010, 04:03

Hi,

I am currently setting up a simulation for Joule heating. So far I have produced a field of the electrical potential V using the laplacian solver. I guess the next step would be to map this field to a simulation that solves the heat equation.

I guess, a solution would be to add a source term to the laplacian solver of the form |grad V|^2 / rho, (rho is the resistivity and a scalar field variable), and then map the potential field V to that and solve.

Here my question: how do I modify the solver equation to include the source term |grad V|^2 / rho?

I guess it's not too tricky but I don't know how to implement the square.

Thanks for your help.

akidess · August 9, 2010, 04:49

I guess the sqr() function is what you are looking for? The programmers guide has a good overview of functions like that in section 1.4.1 / P-23ff (referring to OF1.5).

TomStarke · August 9, 2010, 20:13

Thanks! I will chase that up.

Tom

August 9, 2010, 04:03	Source Term for Joule Heating	#1
TomStarke New Member Join Date: Aug 2010 Posts: 5 Rep Power: 16	Hi, I am currently setting up a simulation for Joule heating. So far I have produced a field of the electrical potential V using the laplacian solver. I guess the next step would be to map this field to a simulation that solves the heat equation. I guess, a solution would be to add a source term to the laplacian solver of the form \|grad V\|^2 / rho, (rho is the resistivity and a scalar field variable), and then map the potential field V to that and solve. Here my question: how do I modify the solver equation to include the source term \|grad V\|^2 / rho? I guess it's not too tricky but I don't know how to implement the square. Thanks for your help.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
momentum source term	zwdi	FLUENT	14	June 27, 2017 16:40
DxFoam reader update	hjasak	OpenFOAM Post-Processing	69	April 24, 2008 02:24
DecomposePar links against liblamso0 with OpenMPI	jens_klostermann	OpenFOAM Bugs	11	June 28, 2007 18:51
UDF Scalar Code: HT 1	Greg Perkins	FLUENT	8	October 20, 2000 13:40
UDFs for Scalar Eqn - Fluid/Solid HT	Greg Perkins	FLUENT	0	October 11, 2000 04:43

August 9, 2010, 04:49		#2
akidess Senior Member Anton Kidess Join Date: May 2009 Location: Germany Posts: 1,377 Rep Power: 30	I guess the sqr() function is what you are looking for? The programmers guide has a good overview of functions like that in section 1.4.1 / P-23ff (referring to OF1.5).

August 9, 2010, 20:13		#3
TomStarke New Member Join Date: Aug 2010 Posts: 5 Rep Power: 16	Thanks! I will chase that up. Tom