• Water Engineering Research
• /
• v.6 no.1
• /
• pp.31-38
• /
• 2005

Groundwater recharge is defined in an addition of water to groundwater reservoir. Recently, many people have been moving to the Edwards aquifer and urban and agricultural industry have been expending. Hydrologists and water planning managers concern about insufficient groundwater amounts and irrigation water price variability. In this paper, I focus on estimates of local recharge volumes and quantify preferential flow through GIS technique. Chloride Mass Balance (CMB) and hydrochemical components have been widely applied to recharge rate and evaluate flow paths. The CMB method is based on relationship between wet-dry chloride deposition data and Rainfall data. These data are manipulated using ArcGIS. Especially, hydrochemical concentration distribution is good index for groundwater residence times or flow paths such as $[Mg^{2+}]/[Ca^{2+}],[Cl]$ and log$([Ca^{2+}]+[Mg^{2+}])/[Na^+]$. Well information such as hydrological-hydrochemical data are imported into ArcGIS and manipulated by interpolation techniques. For each potentiometric surface and water quality, point data are converted to spatial data through each Kriging and Inverse Distance Weighted (IDW) techniques.

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

• Journal of computational fluids engineering
• /
• v.5 no.2
• /
• pp.38-46
• /
• 2000

It has been indicated that the rotor/stator interaction has distinct causes of unsteadiness, such as the viscous vortex shedding, wake/stator interaction and potential rotor/stator interaction. In this paper, the mechanism of unsteady potential interaction in one stage axial compressor is numerically investigated for blade row ratio 1:1 and 2:3 at design point and for blade row ratio 2:3 at off-design point in two-dimensional view point. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting(FVS) and Cubic spline interpolation is applied on zonal interface. In this study the flow unsteadiness due to potential interaction are found to be larger in blade row ratio 2:3 than in 1:1. The total pressure rise in blade row ratio 2:3 is closer to the real value in design point than that in 1:1. The change of unsteady pressure amplitude according to the variation of stator exit pressure is very small.

• 한국전산유체공학회:학술대회논문집
• /
• 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.19 no.2
• /
• pp.49-57
• /
• 2014

In the present study, 2-D gas-kinetic flow solver on unstructured meshes was developed for flows from continuum to transitional regimes. The gas-kinetic BGK scheme is based on numerical solutions of the BGK simplification of the Boltzmann transport equation. In the initial reconstruction, the unstructured version of the linear interpolation is applied to compute left and right states along a cell interface. In the gas evolution step, the numerical fluxes are computed from the evaluation of the time-dependent gas distribution function around a cell interface. Two-dimensional compressible flow calculations were performed to verify the accuracy and robustness of the current gas-kinetic approach. Gas-kinetic BGK scheme was successfully applied to two-dimensional steady and unsteady flow simulations with strong contact discontinuities. Exemplary hypersonic viscous simulations have been conducted to analyze the performances of the gas-kinetic scheme. The computed results show fair agreement with other standard particle-based approaches for both continuum part and transitional part.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.25 no.11
• /
• pp.1518-1526
• /
• 2001

In this paper, the mechanism of unsteady potential interaction and wake interaction in one stage axial turbine is numerically investigated at design point in two-dimensional viewpoint. The numerical technique used is the upwind scheme of Van-Leer's Flux Vector Splitting (FVS) and Cubic spline interpolation is applied on zonal interface between stator and rotor. The inviscid analysis is used to embody the influence of potential interaction only and viscous analysis is used to embody the influences of both potential interaction and wake interaction at the same time. The potential-flow disturbance from the stator into a rotor passage and the periodic blockage effect of rotor produce the unsteady pressure on the blade surface in inviscid analysis. After the wake is cut by rotor, two counterrotating votical patterns flanking the wake centerline in the passage are generated. So, these phenomena magnify the unsteady pressure in viscous analysis than that in inviscid analysis. The resulting unsteady forces on the rotor, generated by the combined interaction of the two effects by potential and wake interaction, are discussed.

• Proceedings of the KSME Conference
• /
• 2000.04b
• /
• pp.599-604
• /
• 2000

In order to study unsteady aerodynamic loads on high speed trains passing by each other at the speed of 350km/h, three-dimensional flow fields around trains during the crossing event are numerically simulated using the three-dimensional Euler equations. The Roe's FDS with MUSCL interpolation is employed to simulate wave phenomena properly. An efficient moving grid system based on domain decomposition techniques is developed to analyze the unsteady flow field induced by the restricted motion of a train on a rail. The numerical simulations of the trains passing by on the double-track are carried out to study the effect of the train nose-shape, the train length and the existence of tunnel when the crossing event occur. Unsteady aerodynamic loads side force and drag force-acting on the train during the crossing are numerically predicted and anlayzed. It is found that the strength of the side force mainly depends on the nose-shape, and that of drag force on tunnel existence. And it is observed that the push-pull like impulsive force successively acts on each car and acts in different directions between the neighborhood cars. The maximum change of the impulsive force reaches about 3 tons. These aerodynamic force data are absolutely necessary for the evaluation of the stability of the high speed multi-car train. The results also indicate the effectiveness of the present numerical method for the simulation of unsteady flow field induced by the bodies in the relative motion.

• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.916-921
• /
• 2003

In most industrial applications, the geometrical complexity is combined with the moving boundaries. These problems considerably increase the computational difficulties since they require, respectively, regeneration and deformation of the grid. As a result, engineering flow simulation is restricted. In order to solve this kind of problems the immersed boundary method was developed. In this study, the immersed boundary method is applied to the numerical simulation of stationary, rotating and oscillating cylinders in the 2-dimensional square cavity. No-slip velocity boundary conditions are given by imposing feedback forcing term to the momentum equation. Besides, this technique is used with a second-order accurate interpolation scheme in order to improve the accuracy of flow near the immersed boundaries. The governing equations for the mass and momentum using the immersed boundary method are discretized on the non-staggered grid by using the finite volume method(FVM). This study presents the possibility of the immersed boundary method to apply to the complex flow experienced in the industrial applications.

• International Journal of Fluid Machinery and Systems
• /
• v.10 no.1
• /
• pp.63-75
• /
• 2017

Many axial fans have circular arc blades with constant thickness. It is still a challenging task to calculate their performance, i.e. to predict how large their pressure rise and pressure losses are. For this task a need for cascade data exists. Therefore, the designer needs a method which works quickly for design purposes. In the present contribution a design method for such cascades consisting of circular arc blades with constant thickness is described. It is based on a singularity method which is combined with a CFD-data-based flow loss model. The flow loss model uses CFD-data to predict the total pressure losses. An interpolation method for the CFD-data are applied and described in detail. Data of measurements are used to validate the CFD-data and parameter variations are conducted. The parameter variations include the variation of the camber angle, pitch chord ratio and the Reynolds number. Additionally, flow patterns of two dimensional cascades consisting of circular arc blades with constant thickness are shown.

• Journal of Biomedical Engineering Research
• /
• v.23 no.2
• /
• pp.101-107
• /
• 2002

This paper presents a segmentation algorithm to extract endocardial contour and epicardial contour of left ventricle in MR Cardiac images. The algorithm is based on a generalized gradient vector flow(GGVF) snake and a prediction of initial contour(PIC). Especially. the proposed algorithm uses physical characteristics of endocardial and epicardial contours, cross profile correlation matching(CPCM), and a mixed interpolation model. In the experiment, the proposed method is applied to short axis MR cardiac image set, which are obtained by Siemens, Medinus, and GE MRI Systems. The experimental results show that the proposed algorithm can extract acceptable epicardial and endocardial walls. We calculate quantitative parameters from the segmented results, which are displayed graphically. The segmented left vents role is visualized volumetrically by surface rendering. The proposed algorithm is implemented on Windows environment using Visual C ++.