抄録
This paper describes the numerical analysis of room air distribution by the finite element method which can easily deal with any domain, the boundary conditions and so on. The two-equation model of turbulence (κ-ε model) is applied to the governing equations of room air, and the discretization of the basic equations is formulated by the penalty finite element method. which results in the replacement of the equation of continuity by u_<j,j>=-λ^<-1> p where λ>>1 is the penalty parameter. The finite element equations are solved by the partitioning method that the equations are partitioned into the momentum equation of the mean flow and the transport equations for turbulence kinetic energy k and that for turbulence energy dissipation rate ε, and the modified Newton method is employed in the iterative procedure. The accuracy and the stability of the scheme by the influence of the penalty parameter are examined for the two dimensional Poiseuille flow. Though the accuracy of the solutions is improved as the penalty parameter is increased, the large parameter makes the condition of the coefficient matrix ill and the numerical convergence is hard to be obtained for the computer. In this computational experiments the scheme has good accuracy even when λ=10^3. At last the numerical example of the three dimensional room model is carried out and the solutions are confirmed to be fully sufficient. As a result, the finite element method is effective for the prediction of the room air distribution.
