Finite Difference Numerical Methods in EES
Finite Difference Numerical Methods in EES
(OP)
Hey there,
So I'm solving for a hypothetical situation of thermal runaway in Li-ion batteries. The goal is to find the ideal insulating material to ensure that the functional battery does not get above 100 C. Assume that the failing battery is separated from a non-failing battery by a 2 mm thick insulation, and that the battery is 10 mm thick and 200 mm tall. Due to the construction of the battery, the thermal conductivity is anisotropic: 1 W m-1 K-1 and 26 W m-1 K-1 in the x and y directions, respectively. The volumetric heat capacity of the battery is 2.2 kJ L-1 k-1. The failing cell is represented by imposing a 500°C constant surface temperature on the left side of the insulation, and the battery is attached on the top and bottom to cold plates for cooling. The cold plate is represented by a constant surface temperature of 25°C on the top and bottom surfaces of the battery. You are also to assume that the right side of the battery is thermally insulated, and that there is no thermal contact resistance between the insulation and the battery.
Here's the idealized picture.
https://fbcdn-sphotos-c-a.akamaihd.net/hphotos-ak-...
I know that this can be approximated with Ansys and a few other softwares, but I need to solve this using numerical methods.
So far, I have begun doing a nodal analysis to solve it as a 2D finite difference problem. I've found a few basic code formats for this in EES, and would like to be able to build off of them. They're attached to this post. I am thinking that it may be wise to run a separate simulation for the insulation and the battery, so that I can plug in conductivity values for the insulation without having to run the entire system.
If you've seen anything like this before and have some input, it would be great. I'll post updates as I move along in the solution.
Cheers
So I'm solving for a hypothetical situation of thermal runaway in Li-ion batteries. The goal is to find the ideal insulating material to ensure that the functional battery does not get above 100 C. Assume that the failing battery is separated from a non-failing battery by a 2 mm thick insulation, and that the battery is 10 mm thick and 200 mm tall. Due to the construction of the battery, the thermal conductivity is anisotropic: 1 W m-1 K-1 and 26 W m-1 K-1 in the x and y directions, respectively. The volumetric heat capacity of the battery is 2.2 kJ L-1 k-1. The failing cell is represented by imposing a 500°C constant surface temperature on the left side of the insulation, and the battery is attached on the top and bottom to cold plates for cooling. The cold plate is represented by a constant surface temperature of 25°C on the top and bottom surfaces of the battery. You are also to assume that the right side of the battery is thermally insulated, and that there is no thermal contact resistance between the insulation and the battery.
Here's the idealized picture.
https://fbcdn-sphotos-c-a.akamaihd.net/hphotos-ak-...
I know that this can be approximated with Ansys and a few other softwares, but I need to solve this using numerical methods.
So far, I have begun doing a nodal analysis to solve it as a 2D finite difference problem. I've found a few basic code formats for this in EES, and would like to be able to build off of them. They're attached to this post. I am thinking that it may be wise to run a separate simulation for the insulation and the battery, so that I can plug in conductivity values for the insulation without having to run the entire system.
If you've seen anything like this before and have some input, it would be great. I'll post updates as I move along in the solution.
Cheers





RE: Finite Difference Numerical Methods in EES
A 2D problem like this one can be simply solved in Excel. Even a 3D version can be set up, if a quite coarse mesh is acceptable.
prex
http://www.xcalcs.com : Online engineering calculations
http://www.megamag.it : Magnetic brakes and launchers for fun rides
http://www.levitans.com : Air bearing pads
RE: Finite Difference Numerical Methods in EES
There's unlikely a closed form solution, but there are numerical approximations that could be taken, like:
http://www.electronicsprotectionmagazine.com/image...
http://www.electronics-cooling.com/2004/05/simple-...
http://www.electronics-cooling.com/1998/01/calcula...
The articles aren't 100% applicable, so you should do a search on something like "heat spreading
The anisotropy makes it interesting. Possibly, you could scale the horizontal dimension of the battery by 20x, but keep the net heat capacity the same?
TTFN

FAQ731-376: Eng-Tips.com Forum Policies
Need help writing a question or understanding a reply? forum1529: Translation Assistance for Engineers
Of course I can. I can do anything. I can do absolutely anything. I'm an expert!
There is a homework forum hosted by engineering.com: http://www.engineering.com/AskForum/aff/32.aspx
RE: Finite Difference Numerical Methods in EES
RE: Finite Difference Numerical Methods in EES
So, is this homework?
TTFN

FAQ731-376: Eng-Tips.com Forum Policies
Need help writing a question or understanding a reply? forum1529: Translation Assistance for Engineers
Of course I can. I can do anything. I can do absolutely anything. I'm an expert!
There is a homework forum hosted by engineering.com: http://www.engineering.com/AskForum/aff/32.aspx
RE: Finite Difference Numerical Methods in EES
RE: Finite Difference Numerical Methods in EES
TTFN

FAQ731-376: Eng-Tips.com Forum Policies
Need help writing a question or understanding a reply? forum1529: Translation Assistance for Engineers
Of course I can. I can do anything. I can do absolutely anything. I'm an expert!
There is a homework forum hosted by engineering.com: http://www.engineering.com/AskForum/aff/32.aspx
RE: Finite Difference Numerical Methods in EES
To use it: set start=1 to start the calculation. If dt is too large, you go into troubles, then use start=0 to come back to initial conditions.
Then instruct Excel to use iterations, to use manual recalc, and a fair number of iterations (1000?).
Each time you hit the recalc button, when recalculation ends, you see the temperatures at time=t
prex
http://www.xcalcs.com : Online engineering calculations
http://www.megamag.it : Magnetic brakes and launchers for fun rides
http://www.levitans.com : Air bearing pads