
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E L A S I M U L A T I O N */
/* D ' U N F L U I D E I N S T A T I O N N A I R E B I D I M E N S I O N N E L I D E A L : */
/* */
/* */
/* Definition : */
/* */
/* Il est classiquement defini par le */
/* systeme d'equations differentielles (ou */
/* 'd' denote les derivees ordinaires, et */
/* 'D', les derivees partielles) : */
/* */
/* dx DH */
/* ---- = + ---- */
/* dt Dy */
/* */
/* dy DH */
/* ---- = - ---- */
/* dt Dx */
/* */
/* z = 0 */
/* */
/* avec : */
/* */
/* 1 | 2 2| */
/* H = ------.Log(|z - w |) */
/* 2.pi | | */
/* */
/* ou 'z' et 'w' sont les deux nombres complexes suivants : */
/* */
/* z = x + i.y = [x,y] */
/* w = a + i.b = [a,b] = [1 + epsilon.cos(omega.t),0] */
/* */
/* d'ou : */
/* */
/* 1 | 2 2 2 2 | */
/* H = ------.Log[|(x - y + 2.i.x.y) - (a - b + 2.i.a.b)|] */
/* 2.pi | | */
/* */
/* 1 | 2 2 2 2 | */
/* H = ------.Log[|(x - y ) - (a - b ) + 2.i.(x.y - a.b)|] */
/* 2.pi | | */
/* __________________________________________ */
/* / 2 */
/* 1 / ( 2 2 2 2 ) 2 */
/* H = ------.Log\ / ((x - y ) - (a - b )) + 4.(x.y - a.b) */
/* 2.pi \/ ( ) */
/* */
/* 2 */
/* 1 [{ 2 2 2 2 } 2] */
/* H = ------.Log[{(x - y ) - (a - b )} + 4.(x.y - a.b) ] */
/* 4.pi [{ } ] */
/* */
/* 2 2 */
/* 1 [{ 2 2 } { 2 2 } 2 2 2 2 2] */
/* H = ------.Log[{(x - y )} + {(a - b )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) ] */
/* 4.pi [{ } { } ] */
/* */
/* 2 2 */
/* 1 [{ 2 2 } 2 2 2 2 2 { 2 2 } ] */
/* H = ------.Log[{(x - y )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) + {(a - b )} ] */
/* 4.pi [{ } { } ] */
/* */
/* 3 2 2 2 2 */
/* DH 1 4.x - 4.x.y - 4.x.(a - b ) + 8.x.y - 8.y.a.b */
/* ---- = ------.---------------------------------------------------------------------- */
/* Dx 4.pi 2 2 */
/* { 2 2 } 2 2 2 2 2 { 2 2 } */
/* {(x - y )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) + {(a - b )} */
/* { } { } */
/* */
/* 3 2 2 2 2 */
/* DH 1 4.y - 4.x .y + 4.y.(a - b ) + 8.x .y - 8.x.a.b */
/* ---- = ------.---------------------------------------------------------------------- */
/* Dy 4.pi 2 2 */
/* { 2 2 } 2 2 2 2 2 { 2 2 } */
/* {(x - y )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) + {(a - b )} */
/* { } { } */
/* */
/* ou en simplifiant : */
/* */
/* 3 { 2 2 2 } */
/* DH 1 4.x + 4.x.{y - (a - b )} - 8.y.a.b */
/* ---- = ------.---------------------------------------------------------------------- */
/* Dx 4.pi 2 2 */
/* { 2 2 } 2 2 2 2 2 { 2 2 } */
/* {(x - y )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) + {(a - b )} */
/* { } { } */
/* */
/* 3 { 2 2 2 } */
/* DH 1 4.y + 4.y.{x + (a - b )} - 8.x.a.b */
/* ---- = ------.---------------------------------------------------------------------- */
/* Dy 4.pi 2 2 */
/* { 2 2 } 2 2 2 2 2 { 2 2 } */
/* {(x - y )} - 2.(x - y ).(a - b ) + 4.(x.y - a.b) + {(a - b )} */
/* { } { } */
/* */
/* d'ou : */
/* */
/* dx DH */
/* ---- = F (x,y,z,t) = + ---- */
/* dt x Dy */
/* */
/* dy DH */
/* ---- = F (x,y,z,t) = - ---- */
/* dt y Dx */
/* */
/* z = 0 */
/* */
/* */
/* On va generaliser le nombre complexe 'w' par : */
/* */
/* w = [R + R .cos(R .t),I + I .sin(I .t)] */
/* 1 2 3 1 2 3 */
/* */
/* */
/*************************************************************************************************************************************/