Farouki_a(n)
- Calculates the value 'a' from the expression 3a² - 2a² = n, where n is the fractional porosity.
Units
- ()
Inputs
- fractional porosity ().
Returns
- return 'a'
Reference
Farouki T. Omar, Cold Regions Science and Technology, 5 (1981) 67-75. The thermal properties of soils in cold regions
Code
float Farouki_a(float fract_por) {
float a;
float nnew = 0.0;
while(fabs(fract_por - nnew) > 0.001) {
a += (fract_por - nnew)*0.25;
nnew = 3*a*a - 2*a*a*a;
}
return a;
}