Skip to content

Commit

Permalink
Merge pull request #4 from cbm755/deps
Browse files Browse the repository at this point in the history 
Update README with installation instructions
  • Loading branch information
@hanrach
hanrach committed Aug 16, 2021
2 parents 00f817e + a05fb4f commit 8b1fa99
Showing 1 changed file with 8 additions and 1 deletion.
9 changes: 8 additions & 1 deletion README.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -4,4 +4,11 @@
This repository contains the solver for pseudo-2D model of Li-ion battery based on Finite Difference Method with automatic differentiation via [JAX](https://github.com/google/jax).

The model is taken from the P2D model outlined in the [paper](http://web.mit.edu/braatzgroup/Torchio_JElectSoc_2016.pdf) by Torchio et al.
Run `examples/main_reorder.py` to get solution.

## Getting started

* `pip install numpy scipy` (or get these from your OS)
* `pip install jax[cpu]`
* `pip install scikit-umfpack`

Then run `examples/main_reorder.py`.

2 comments on commit 8b1fa99

@jsfeng-fudan
Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Hi,@hanrach
I study to use the code to compute the heat during charge and discharge, for exchange, in naive method, the heat can be collected by : Qohm, Qrxn and Qrev,
in p2d_main_fn.py,

图片

    print("-------------------------------- positive electrode ------------------------------")
    print("heat generated  by Ohm mechanism in positive electrode : ")
    print(peq.Qohm(phis_ne,phis_pe,phie_ne,phie_pe,uvec_pe, uvec_sep, uvec_ne,Tvec_pe))
    print("heat generated  by chemiacal reaction in positive electrode :")
    print(peq.Qrxn(jvec_pe,eta_pe))
    print("heat generated  by entropy change in positive electrode :")
    pcmax =51554
    print(peq.Qrev(jvec_pe,Tvec_pe,cmat_pe,cmat_pe,51554))

I can obtain the Ohm heat and reaction heat as below
图片

but I cannot arrive at the reversible heat. The Qrev() function prototype is Qrev(j,Tc,cM,cMp,cmax). I dont know the meaning of cM and cMp, I can usig the umat_pe to replace them.
can you help me to solve this problem.

by the way, the Ohm heat is so large ,is it reasonable?

@cbm755
Copy link
Collaborator

@cbm755 cbm755 commented on 8b1fa99 Sep 27, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'm not sure about the Ohm heat. I recall the problem has a wide variety of scales... I filed #7 to help track this down

Please sign in to comment.