Initializing model with position dependent values

[philomenaweng]

Is there a way to initialize a DFN model with a position-dependent values (i.e. x-dependent porosity or SEI thickness)?

[valentinsulzer]

Porosity, active material fraction, initial concentration: yes. Just pass in a function

def eps(x):
    return ...

parameter_values["Separator porosity"] = eps

The function takes in x in meters relative to the anode current collector.

SEI thickness: not right now but it could be added fairly easily.

[philomenaweng]

thank you!!