• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2009년 추계학술대회논문집
• /
• pp.49-55
• /
• 2009

This paper evaluates performances of a recently developed divergence-free finite element method based on Hermite interpolated stream functions. Velocity bases are derived from Hermite interpolated stream functions to form divergence-free basis functions. These velocity basis functions constitute a solenoidal function space, and the simple gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into a solenoidal and an irrotational parts, and the decoupled Navier-Stokes equations are projected onto their corresponding spaces to form proper variational formulations. To access accuracy and convergence of the present algorithm, three test problems are selected. They are lid-driven cavity flow, flow over a backward-facing step and buoyancy-driven flow within a square enclosure. Hermite interpolation functions from cubic to quintic are chosen to run the test problems. Numerical results are shown. In all cases it has shown that the present method has performed well in accuracies and convergences. Moreover, the present method does not require an upwinding or a stabilized term.

• Journal of Mechanical Science and Technology
• /
• 제15권6호
• /
• pp.813-820
• /
• 2001

Solutions of inviscid rotational flows near the corners of an arbitrary angle and within a triangle of arbitrary shapes are presented. The corner-flow solutions has a rotational component as a particular solution. The addition of irrotatoinal components yields a general solution, which is indeterminate unless the far-field condition is imposed. When the corner angle is less than 90$^{\circ}$the flow asymptotically becomes rotational. For the corner angle larger than 90$^{\circ}$it tends to become irrotational. The general solution for the corner flow is then applied to rotational flows within a triangle (Method I). The error level depends on the geometry, and a parameter space is presented by which we can estimate the error level of solutions. On the other hand, Method II employing three separate coordinate systems is developed. The error level given by Method II is moderate but less dependent on the geometry.