BCs #134
BCs #134
Conversation
RemDelaporteMathurin
left a comment
RemDelaporteMathurin
left a comment
There was a problem hiding this comment.
First pass on the Concentration part.
| name: python3 | ||
| --- | ||
|
|
||
| # Concentration # |
There was a problem hiding this comment.
I think instead of Concentration + Particle Flux + Surface reactions, we can have everything under "Hydrogen Transport"?
|
|
||
| ## Defining fixed concentration ## | ||
|
|
||
| BCs need to be assigned to surfaces using FESTIM's `SurfaceSubdomain` class. To define the concentration of a specific species, we can use `FixedConcentrationBC`: |
There was a problem hiding this comment.
| BCs need to be assigned to surfaces using FESTIM's `SurfaceSubdomain` class. To define the concentration of a specific species, we can use `FixedConcentrationBC`: | |
| BCs need to be assigned to surfaces using FESTIM's `SurfaceSubdomain` class. To set the concentration of a specific species on a boundary, we can use `FixedConcentrationBC`: |
|
|
||
| $$ | ||
|
|
||
| BC = 10 + x^2 + T(t+2) |
There was a problem hiding this comment.
| BC = 10 + x^2 + T(t+2) | |
| c = 10 + x^2 + T(t+2) |
| BC = 10 + x^2 + T(t+2) | ||
|
|
||
| $$ | ||
|
|
There was a problem hiding this comment.
| To do so, we define a python function. |
| ```{code-cell} ipython3 | ||
| my_custom_value = lambda x, t, T: 10 + x[0]**2 + T *(t + 2) | ||
|
|
||
| my_bc = FixedConcentrationBC(subdomain=boundary, value=my_custom_value, species=H) | ||
| ``` |
There was a problem hiding this comment.
Here we should add a note block showing that this is equivalent to
def my_custom_value(x, t, T):
return 10 + x[0]**2 + T *(t + 2)| my_bc = FixedConcentrationBC(subdomain=boundary, value=my_custom_value, species=H) | ||
| ``` | ||
|
|
||
| ## Choosing a solubility law ## |
There was a problem hiding this comment.
Same comment as above. We need to show more about the integration. For example, show that this is equivalent to a FixedConcentrationBC but it's just providing convenience
|
|
||
| $$ S_{ext} = \varphi \cdot f(x) $$ | ||
|
|
||
| $$\varphi = 1\cdot 10^{13} \quad \mathrm{m}^{-2}\mathrm{s}^{-1}$$ |
There was a problem hiding this comment.
| $$\varphi = 1\cdot 10^{13} \quad \mathrm{m}^{-2}\mathrm{s}^{-1}$$ | |
| $$\varphi = 10^{13} \quad \mathrm{m}^{-2}\mathrm{s}^{-1}$$ |
| my_bc = HenrysBC(subdomain=3, H_0=1.5, E_H=0.2, species=H, pressure=pressure_value) | ||
| ``` | ||
|
|
||
| ## Plasma implantation approximation ## |
There was a problem hiding this comment.
Maybe this one could go in a section "Hydrogen transport: advanced"
| my_bc = HenrysBC(subdomain=3, H_0=1.5, E_H=0.2, species=H, pressure=pressure_value) | ||
| ``` | ||
|
|
||
| ## Plasma implantation approximation ## |
There was a problem hiding this comment.
For this application, because it's a "real-world" application, I would use real numbers.
Typically, the implantation range is a few nanometers for example.
| my_bc = HenrysBC(subdomain=3, H_0=1.5, E_H=0.2, species=H, pressure=pressure_value) | ||
| ``` | ||
|
|
||
| ## Plasma implantation approximation ## |
There was a problem hiding this comment.
There's actually an issue here: this is not a boundary condition approximation, you set an actual volumetric source.
Boundary conditions chapter: