fargOCA issueshttps://gitlab.oca.eu/DISC/fargOCA/-/issues2022-06-29T16:31:11Zhttps://gitlab.oca.eu/DISC/fargOCA/-/issues/544Viscosity alpha specification2022-06-29T16:31:11ZAlain O' MiniussiViscosity alpha specificationAs discussed in #524
In this boundary condition @morby we would like to change the vector alpha (viscosity.value in VT_ALPHA) with time with something like:
$`\alpha(Rmed[i]) = (\alpha_M-\alpha_m)*\exp(-t/T)+\alpha_m`$
this is not p...As discussed in #524
In this boundary condition @morby we would like to change the vector alpha (viscosity.value in VT_ALPHA) with time with something like:
$`\alpha(Rmed[i]) = (\alpha_M-\alpha_m)*\exp(-t/T)+\alpha_m`$
this is not possible at the moment since alpha is defined as const
```
reals const& alpha = disk.physic().viscosity.alpha();
```
I think that the way to do is to introduce another viscosity type,
and provide $`\alpha_M`$ and $`\alpha_m`$ as parameters.
The reason is that we will also have to change the vector alpha as a function of the evolution with time of the disk temperature and density.
Therefore, we will need also to change the call to the viscosity field.
But may be there is a better way to do.
@alainm : what do you suggest ?Viscosity ReorganisationAlain O' MiniussiAlain O' Miniussi