• Journal of Ocean Engineering and Technology
• /
• v.25 no.5
• /
• pp.9-14
• /
• 2011

The divergence-free finite elements introduced in this paper are derived from Hermite functions, which interpolate stream functions. Velocity bases are derived from the curl of the Hermite functions. These velocity basis functions constitute a solenoidal function space, and the gradient of the Hermite functions constitute an irrotational function space. The incompressible Navier-Stokes equation is orthogonally decomposed into its solenoidal and irrotational parts, and the decoupled Navier-Stokes equations are then projected onto their corresponding spaces to form appropriate variational formulations. The degrees of the Hermite functions we introduce in this paper are bi-cubis, quartic, and quintic. To verify the accuracy and convergence of the present method, three well-known benchmark problems are chosen. These are lid-driven cavity flow, flow over a backward facing step, and buoyancy-driven flow within a square enclosure. The numerical results show good agreement with the previously published results in all cases.

• 한국전산유체공학회:학술대회논문집
• /
• 2009.11a
• /
• 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 computational fluids engineering
• /
• v.13 no.4
• /
• pp.13-23
• /
• 2008

This paper is an extension of previous study[1] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite (serendipity) cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires four degrees-of-freedom at each element corners. Those degrees-of-freedom are the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational basis functions from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[2].

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.61-65
• /
• 2006

This paper describes a recent development on the divergence free basis function based on a hermite stream function. The well-known cavity problem has been used to compare the accuracy and the convergence of the present method with those of a modified residual method known as one of the stabilized finite element methods. The comparison showed the present method performs better in the accuracy and convergence.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2006.11a
• /
• pp.411-414
• /
• 2006

This paper describes a recent development on the divergence free basis function based on a hermite stream function. The well-known cavity problem has been used to compare the accuracy and the convergence of the present method with those of a modified residual method known as one of the stabilized finite e1ement methods. The comparison showed the present method performs better in the accuracy and convergence.

• 한국전산유체공학회:학술대회논문집
• /
• 2008.03a
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국전산유체공학회:학술대회논문집
• /
• 2008.10a
• /
• pp.33-41
• /
• 2008

This paper is an extension of previous study[9] on a development of a divergence-free element method using a hermite interpolated stream function. Divergence-free velocity bases defined on rectangles derived herein produce pointwise divergence-free flow fields. Hence the explicit imposition of continuity constraint is not necessary and the Galerkin finite element formulation for velocities does not involve the pressure. The divergence-free element of the previous study employed hermite serendipity cubic for interpolation of stream function, and it has been noted a possible discontinuity in variables along element interfaces. This deficiency can be removed by use of a hermite bicubic interpolated stream function, which requires at each element corners four degrees-of-freedom such as the unknown variable, its x- and y-derivatives and its cross derivative. Detailed derivations are presented for both solenoidal and irrotational bases from the hermite bicubic interpolated stream function. Numerical tests are performed on the lid-driven cavity flow, and results are compared with those from hermite serendipity cubics and a stabilized finite element method by Illinca et al[7].

• 한국전산유체공학회:학술대회논문집
• /
• 2007.10a
• /
• pp.220-225
• /
• 2007

This paper is a continuation of the recent development on the hermite-based divergence free basis function and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square cavity with temperature difference across the two sides. The basis functions for the velocities consist of the hermite function and its curl. However, the basis for the temperature are the hermite function and its gradienst. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x- and y-derivatives. Numerical results for the streamlines, the temperatures, the x-velocities and the y-velocities show good agreements with those of De vahl Davis[7].

• Journal of computational fluids engineering
• /
• v.12 no.1
• /
• pp.35-42
• /
• 2007

This paper describes a recent development on the divergence free basis function based on a hermite stream function and verifies its validity by comparing results with those from a modified residual method known as one of stabilized finite element methods. It can be shown that a proper choice of degrees of freedom at a node with a proper arrangement of the hermite interpolation functions can yield solenoidal or divergent free interpolation functions for the velocities. The well-known cavity problem has been chosen for validity of the present algorithm. The comparisons from numerical results between the present and the modified residual showed the present method yields better results in both the velocity and the pressure within modest Reynolds numbers(Re = 1,000).

• Journal of Ocean Engineering and Technology
• /
• v.23 no.5
• /
• pp.1-8
• /
• 2009

This paper is a continuation of the recent development on Hermite-based divergence free element method and deals with a non-isothermal fluid flow thru the buoyancy driven flow in a square enclosure with temperature difference across the two sides. The basis functions for the velocity field consist of the Hermite function and its curl while the basis functions for the temperature field consists of the Hermite function and its gradients. Hence, the number of degrees of freedom at a node becomes 6, which are the stream function, two velocities, the temperature and its x and y derivatives. This paper presents numerical results for Ra = 105, and compares with those from a stabilized finite element method developed by Illinca et al. (2000). The comparison has been done on 32 by 32 uniform elements and the degree of approximation of elements used for the stabilized finite element are linear (Deg. 1) and quadratic (Deg. 2). The numerical results from both methods show well agreements with those of De vahl Davi (1983).