/* funzione per interpolare l'energia */
Double_t en ( Double_t *x, Double_t *par ) {
// par[1] = .186;
// par[0] = 1.68;
/* variabile ausiliaria: par[0] = $\hbar c$ */
Double_t l = TMath::Power( B * par[0] / x[0], 1./3. );
/*
* termine tra parentesi (ho sfruttato le funzioni note per non
* dover riscrivere un'altra funzione per calcolare $l d\beta/dl$)
*/
Double_t tmp = l * betad( l, 1. ) - beta( l, 1. );
/* moltiplico per il coefficiente */
tmp = tmp * l * l / ( par[1] * B );
/* primo termine fuori dalle parentesi */
tmp += l * betad( l, 0. );
/* incremento di uno */
tmp += 1.;
/* 'par[2]' è la densità d'energia $e_0$ */
return par[2] - tmp / ( 3. * B * x[0] * x[0] );
}
tmp<volSymmTensorField> Foam::relativeVelocityModel::tauDm() const
{
volScalarField betac(alphac_*rhoc_);
volScalarField betad(alphad_*rhod_);
// Calculate the relative velocity of the continuous phase w.r.t the mean
volVectorField Ucm(betad*Udm_/betac);
return tmp<volSymmTensorField>
(
new volSymmTensorField
(
"tauDm",
betad*sqr(Udm_) + betac*sqr(Ucm)
)
);
}
