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: 17	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.

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

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: 31	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: 17	Thanks! I will chase that up. Tom