• 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].

• 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).

• 한국전산유체공학회:학술대회논문집
• /
• 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].

• Journal of computational fluids engineering
• /
• v.14 no.4
• /
• pp.67-77
• /
• 2009

This paper is a continuation of a recent development on the Hermite-based divergence-free element method and deals with a non-isothermal fluid flow driven by the buoyancy force in a square cavity with temperature difference across the two sides. Two Hermite functions are considered for numerical computations in this paper. One is a cubic function and the other is a quartic function. The degrees-of-freedom of the cubic Hermite function are stream function and its first and second derivatives for the velocity field, and temperature and its first derivatives for the temperature field. The degrees-of-freedom of the quartic Hermite function include two second derivatives and one cross derivative of the stream function in addition to the degrees-of-freedom of the cubic stream function. This paper presents a brief review on the Hermite based divergence-free basis functions and its finite element formulations for the buoyancy driven flow. The present algorithm does not employ any upwinding or a stabilization term. However, numerical values and contour graphs for major flow variables showed good agreements with those by De Vahl Davis[6].

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.4
• /
• pp.55-63
• /
• 1995

This paper describes the procedure to analyze the performance of the end-plated propeller(EPP) by a boundary integral method. The screw blade(SB) and end-plate(EP) are represented by a set of quadrilateral panels, where the source and normal dipole of uniform strength are distributed. The perturbation velocity potential, being the only unknown via the potential-based formulation, is determined by satisfying the flow tangency condition on the blade and the end-plate at the same time. The Kutta condition is satisfied through an iterative process by requiring the null pressure jump across the upper and lower sides of the trailing edges of both the SH and the EP. Sample calculations indicate that the EP increases the loading near the tip of the SB while spreading the trailing vortices along the trailing edge of the EP, thus avoiding the strong tip-vortex formation. Predicted performance of the EPP shows good correlations with the experimental results. The method is therefore considered applicable in designing and analyzing the EPP which may be an alternative for energy-saving propulsive devices.

• Journal of the Society of Naval Architects of Korea
• /
• v.34 no.4
• /
• pp.12-20
• /
• 1997

A potential-based panel method has been developed to investigate the resistance characteristics of a high speed catamaran advancing on the free surface. Normal dipoles and sources are distributed on the body surface while sources are distributed on the free surface. Linearised free surface conditions are used in the present analysis. To avoid the instabilities due to the velocity difference between inner and outer flow of a high speed catamaran, Kutta condition has been applied at the stern. Model test has been carried out not only to validate the numerical results but to confirm the capabilities of a CWC(Circulating Water Channel). It is believed that we can obtain the qualitatively reasonable results in the CWC. Computed results are compared with those of experiments and Insel's experimental values. Since the Kutta condition is applied at the stern, stable solutions are obtained at the high speed range. The present method, using linearised free surface conditions at the high speed range, seems to be a useful tool in the hull form design of a high speed catamaran.

• The Journal of Korean Institute of Communications and Information Sciences
• /
• v.18 no.8
• /
• pp.1059-1068
• /
• 1993

We will put the emphasis on the analysis of arbitrarily shaped microstrip antennas. The most general and rigorous treatment of microstrip antennas is given by the electric field integral equation(EFIE), usally formulated in the spectral domain. In this paper, we use a modification of EFIE, called the mixed potential integral equation(MPIE) , and we solve it in the space domain. This technique uses Green's functions associated with the scalar and vector potential which are calculated by using stratified media theory and are expressed as Sommerfeld integrals. The integral equation is solved by a moment's method using rooftop subsectional basis function. Thus, microstrip patches of any shape can be analysed at any frequency and for any substrate. Numerical results for a rectangular patch and for a L-shaped patch are given and compared with measured values.

• Journal of the Society of Naval Architects of Korea
• /
• v.36 no.2
• /
• pp.9-21
• /
• 1999

The waterjet is very widely used today in propulsion system of high speed vessel but manufacturer is limited because of an efficiency and a difficulty of a manufacture on the inlet configuration of the waterjet. The importancy in designing the inlet configuration of the waterjet lies on the minimization of the cavitation which is largely affecting the efficiency. In this paper, the configuration analysis is carried out to find a optimum shape which is minimizing the cavitation using a two dimensional potential-based panel method with an inlet configuration of a flush type. Also, it is developed a direct design method finding an inlet configuration by a given pressure distributions. The numerically obtained optimum shape using this configuration analysis method show a good agreement compared to the Kashiwadani's results. It is carried out a direct design method over a lip and a ramp of an inlet configuration wish pressure distributions obtained a result of the configuration analysis and the results show a good agreement compared to original configuration.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.38 no.7
• /
• pp.637-646
• /
• 2010

A time-domain panel method based on the potential flow theory and the time-stepping method is developed to predict the steady/unsteady aerodynamic characteristics of FM07, which is the BWB (Blended-wing body) type MAV. In the aerodynamic analyses, we used two types of the initial model(Case I) and the improved model(Case II), which is moved the gravity center toward the rear and has larger aspect ratio. In the steady aerodynamic analyses, it is revealed that improved model has higher lift to drag ratio(L/D) and more stable pitch characteristic than those of the initial model. In the unsteady aerodynamic analyses for sudden acceleration motion similar to the launch phase of MAV, it seemed that there is a rapid increase of the lift coefficient after the launch and unsteady results are good agreed compare with steady results in just a few times. In the analysis for pitch oscillation motion, which is occurred at the cruise condition of the FM07, it shows that unsteady aerodynamic coefficients looped around steady results and the improved model has more sensitive aerodynamic characteristics.