public double Shell_Cd_infiltrate(double Cd_p1,int shell_index,double Cd_p2)
{
//assuming Cd_p2 V > Cd_p1_V
//attention! it could also be negative
//fisk law
double contact_surface =Add_Math.s_circle((1 + shell_index) * shell_radius_unit);
return contact_surface*(Cd_p2-Cd_p1)*0.0015;
//lack of 0 protect
}
public double sum_env_infiltrate_ball(double env_cd,int ball_num,double ball_surface_cd)
{
//undifined
//attention! it could also be negative
//fisk law
double contact_surface = Add_Math.s_circle(ball_radius) * ball_num;//sum surface
double Cd_trans = (env_cd - ball_surface_cd) * 0.0025;
return Cd_trans * contact_surface;
//return 0.01;
}
