×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Jobs

Finite Difference Numerical Methods in EES

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

RE: Finite Difference Numerical Methods in EES

What's 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

The wording you used suggests a homework problem, but your other postings do not.

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

(OP)
@prex EES is Engineering Equation Solver. It's a program used frequently for thermodynamic and fluids calculations. The coarsest mesh possible is 2mm because it is necessary to have a nod on the outside of the insulation and one on the contact surface with the battery. If you have done anything similar to this in Excel, could I see an example? I'm much more comfortable with excel myself.

RE: Finite Difference Numerical Methods in EES

Upon rereading your OP and this last reply, I'm a bit confused. You seem to proposing to reinvent the wheel, as it were, since you're talking about doing a finite-element solution, which is exactly what ANSYS is.

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

(OP)
@IRstuff I realized after your first comment that this should definitely be in the homework forum, so I went ahead and threw it in there as well. Thanks for the heads up on that. I just kept it in here as well because of prez's reply for convenience. Do I need to take it down from this forum?

RE: Finite Difference Numerical Methods in EES

No, that's fine. Engineering.com's homework forum is still sparsely populated, and there may be some utility to having this thread staying here.

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

You have it. It's a relatively rough setup (in 2D), that can easily be made better (e.g. dx different from dy, use of symmetry,...). It also possibly contains errors.
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

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members!


Resources