柱坐标系和球坐标系下NS方程的直接推导.doc

Derivation of 3D Euler and Navier-Stokes Equations in Cylindrical Coordinates Dingxi Wang School of Engineering, Durham Contents 1. Derivation of 3D Euler Equation in Cylindrical coordinates 2. Derivation of Euler Equation in Cylindrical coordinates moving at in tangential direction 3. Derivation of 3D Navier-Stokes Equation in Cylindrical Coordinates 1. Derivation of 3D Euler Equation in Cylindrical coordinates Euler Equation in Cartesian coordinates (1.1) Where ? Conservative flow variables ? Inviscid/convective flux in x direction ? Inviscid/convective flux in y direction ? inviscid/convective flux in z direction And their specific definitions are as follows ,,, ? Total enthalpy Some relationship We want to perform the following coordinates transformation Because According to Cramer’s ruler, we have (1.2.1) (1.2.2) Where Similar to the above (1.2.3) (1.2.4) In addition, the following relations hold between cylindrical coordinate and Cartesian coordinate , ,,,, (1.3) (1.4.1) (1.4.2) Derivation Multiplying the both side of equation (1.1) by and applying equalities (1.4.1) and (1.4.2) gives, (1.5) Differentiating the following w.r.t. time gives , , (1.6.1) (1.6.2) Expanding the term and applying the relationships (1.6) yields, (1.7.1) Expanding the term and applying the relationships (1.6) yields, (1.7.2) Substituting relationships (1.7) into equation (1.5) and rearranging gives, (1.8) As we can see from expressions (1.7), the momentum equations in radial and tangential directions contain velocities in Cartesian coordinate; we need to replace them with corresponding variables in cylindrical coordinate. Writing down the momentum equations in radial and tangential directions as follows, (1.9.1) (1.9.2) Multiplying (1.9.1) by and (1.9.2) by, then summing up and applying expressions (1.6) and rearran