search

pdelab / modular-pnp

This repository contains a modularized implementation of the PNP solver in which functionalities are organized into Newton solver techniques, conversion methods between FEniCS and FASP, as well as geometric and parameter functionalities. An additional library specific to strictly PNP functionalities is also within.
https://thepnpsolver.github.io/
3 stars 0 forks source link

PNP benchmark #39

Open maximilianmetti opened 8 years ago

maximilianmetti commented 8 years ago

We need a benchmark to validate the code can solve a physically relevant system.

maximilianmetti commented 8 years ago

we can implement a simple test that verifies Ohm's Law that Voltage = Current * Resistance. This requires computing the current through a cross-sectional surface and comparing it to the boundary conditions for potential. Then, based on this current, we can compute the resistance and see if it matches that predicted by Ohm's law

maximilianmetti commented 8 years ago

Since we do not know the resistance constant of the system in terms of the charge concentrations (at the boundary), we have to compute a Current-Voltage relation and verify that the relationship is nearly linear to indicate that there's a constant (resistance) that describes the relationship. We can do this by solving the system with the BCs for the charges fixed and varying the voltage at the boundary.

maximilianmetti commented 8 years ago

This should move up in priority... It should be the first thing we do once we successfully implement the newton class for the refactor of the code.