• The KIPS Transactions:PartA
• /
• v.8A no.4
• /
• pp.503-508
• /
• 2001

Conformal mapping is useful to solve problems in heat conduction, electrostatic potential and fluid flow involving Laplace's equation in two independent variables. Determinations of conformal maps from the unit disk onto a Jordan region eventually requires solving the Theodorsen equation which is in general nonlinear with respect to the boundary correspondence function. H bner's method which has been well known for the efficient method among the many suggestions for the Theodorsen equation, was improved in early study[1, 2]. In this paper Wegmann's method is treated that is more efficient in computation cost rather than H bner's. But we found that a question which is divergent in some difficult problems by numerical experiment of Wegmann's iteration. We analyze theoretically the cause of divergence and propose an improved method by applying a low frequency filter to the Wegmann's method. Numerical experiments by our improved method show convergence for all divergent problems by Wegmann's method.

• Journal of applied mathematics & informatics
• /
• v.25 no.1_2
• /
• pp.67-83
• /
• 2007

The characteristics of steady two-dimensional laminar boundary layer flow of a viscous and incompressible fluid past a moving wedge with suction or injection are theoretically investigated. The transformed boundary layer equations are solved numerically using an implicit finite-difference scheme known as the Keller-box method. The effects of Falkner-Skan power-law parameter (m), suction/injection parameter ($f_0$) and the ratio of free stream velocity to boundary velocity parameter (${\lambda}$) are discussed in detail. The numerical results for velocity distribution and skin friction coefficient are given for several values of these parameters. Comparisons with the existing results obtained by other researchers under certain conditions are made. The critical values of $f_0$, m and ${\lambda}$ are obtained numerically and their significance on the skin friction and velocity profiles is discussed. The numerical evidence would seem to indicate the onset of reverse flow as it has been found by Riley and Weidman in 1989 for the Falkner-Skan equation for flow past an impermeable stretching boundary.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.3
• /
• pp.933-944
• /
• 1996

In this paper, two new techniques, the pattern filling and the refined flow field regeneration, based on the finite element method and Eulerian mesh advancement approach have been developed to analyze incompressible viscous flow with free surfaces. The gorerning equation for flow analysis is Navier-Stokes equation including inertia and gravity effects. The penalty and Newton-Raphson methods are used effectively for finite element formulation. The flow front surface and the volume inflow rate are calculated using the pattern filling technique to select an adequate pattern among five filling patterns at each quadrilateral control volume. By the refined flow field regeneration technique, the new flow field which renders better prediction in flow surface shape is generated and the velocity field at the flow front part is calculated more exactly. Using the new thchniques to be developed, the dam-breaking problem has been analyzed to predict flow phenomenon of fluid and the predicted front positions versus time have been compared with the reported experimental result.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.16 no.2
• /
• pp.42-50
• /
• 2012

The goal of heat transfer studies is the accurate prediction of temperature and heat flux distribution on material boundaries. To this purpose, general-purpose computational fluid dynamics(CFD) code is used : FLUENT. Mass fluxes and pressure ratio are calculated for two types of nozzle. The comparative studies reveal that the computational results are in agreement with the experimental data. Also, heat transfer coefficients from FLUENT for one type of nozzle are very similar and agree well with the experimental data in the diverging part of the nozzle, but the calculated results are large in the converging part. The heat transfer coefficients from Bartz equation are over-predicted. We can consider various reasons for these differences, i.e., laminarization by the highly accelerated flow in the nozzle, turbulent flow model and grid generation.

• Journal of KSNVE
• /
• v.6 no.5
• /
• pp.607-616
• /
• 1996

A combined computational fluid dynamics(CFD)-Kirchhoff method is presented for predicting high-speed impulsive noise generated by a hovering blade. Two types of Kirchhoff integral formula are used; one for the classical linear Kirchhoff formulation and the other for the nonlinear Kirchhoff formulation. An Euler finite difference solver is solved first to obtain the flow field close to the blade, and then this flow field is used as an input to a Kirchhoff formulation to predict the acoustic far-field. These formulas are used at Mach numbers of 0.90 and 0.95 to investigate the effectiveness of the linear and nonlinear Kirchhoff formulas for delocalized flow. During these calculiations, the retarded time equation is also carefully examined, in particular, for the cases of the control surface located outside of the sonic cylinder, where multiple roots are obtained. Predicted results of acoustic far-field pressure with the linear Kirchhoff formulation agree well with experimental data when the control surface is at the certain location(R=1.46), but the correlation is getting worse before or after this specific location of the control surface due to the delocalized nonlinear aerodynamic flow field. Calculations based on the nonlinear Kirchhoff equation using a linear sonic cylinder as a control surface show a reasonable agreement with experimental data in negative amplitudes for both tip Mach numbers of 0.90 and 0.95, except some computational integration problems over a shock. This concliudes that a nonlinear formulation is necessary if the control surface is close to the blade and the flow is delocalized.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.18 no.3
• /
• pp.34-39
• /
• 2014

A study is analyzed on the design factor of center-body diffuser and performed on conceptual design of center-body diffuser with computational fluid dynamic. The flow field of center-body diffuser is calculated using axisymmetric two-dimensional Navier-Stokes equation with $k-{\epsilon}$ turbulencemodel. The center-body diffuser is compared with second throat exhaust diffuser in terms of starting pressure, the degree of vacuum pressure and the design factors. The counter flow jet on cone-tip of the center-body is applied for thermal protection system in the center-body diffuser.

• Water Engineering Research
• /
• v.2 no.3
• /
• pp.179-185
• /
• 2001

The time consuming and expensive nature of experimental research on scouring processes caused by flowing water makes it attractive to develop numerical tools for the predication of the interaction of the fluid flow and the movable bed. In this paper the numerical simulation of scour by a wall jet is presented. The flow is assumed to be two-dimensional, and the alluvium is cohesionless. The solution process, repeated at each time step, involves simulation of a turbulent wall jet flow, solution of the convection-diffusion of sand concentration, and prediction of the bed deformation. For simulation of the jet flow, the governing equations for momentum, mass balance and turbulent parameters are solved by the finite volume method. The SIMPLE scheme with momentum interpolation is used for pressure correction. The convection-diffusion equation is solved for sediment concentration. A boundary condition for concentration at the bed, which takes into account the effect of bed-load, is implemented. The time rate of deposition and scour at the bed is obtained by solving the continuity equation for sediment. The shape and position of the scour hole and deposition of the bed material downstream of the hole appear realistic.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.17 no.4
• /
• pp.43-48
• /
• 2013

A study is analyzed on the adverse pressure gradient and the transient regime of supersonic diffuser with Computational Fluid Dynamic. The flow field of supersonic diffuser is calculated using Axisymmetric two-dimensional Navier-Stokes equation with $k-{\epsilon}$ turbulence model. The transient simulation is compared in terms of mach number and static temperature of vacuum chamber according to pressure variation of rocket engine combustion chamber. Combustion gas flow into the vacuum chamber during operation of the supersonic diffuser. According to this phenomenon, the pressure and the temperature rise in the vacuum chamber were observed. Thus, the protection system will be necessary to prevent the pressure and temperature rise in the transition process during operation of the subsonic diffuser.

• Honam Mathematical Journal
• /
• v.43 no.2
• /
• pp.183-197
• /
• 2021

Recently several authors have extended the Beta function by using its integral representation. However, in many cases no expression of these extended functions in terms of classic special functions is known. In the present paper, we introduce a further extension by defining a family of functions G_r,s : ℝ^*₊ → ℂ, with r, s ∈ ℂ and ℜ(r) > 0. For given r, s, we prove that this function satisfies a second-order linear differential equation with rational coefficients. Solving this ODE, we express G_r,s as a combination of confluent hypergeometric functions. From this we deduce a new integral relation satisfied by Tricomi's function. We then investigate additional specific properties of G_r,1 which take the form of new non trivial integral relations involving exponential and error functions. We discuss the connection between G_r,1 and Stokes' first problem (or Rayleigh problem) in fluid mechanics which consists in determining the flow created by the movement of an infinitely long plate. For $r{\in}{\frac{1}{2}}{\mathbb{N}}^*$, we find additional relations between G_r,1 and Hermite polynomials. In view of these results, we believe the family of extended beta functions G_r,s will find further applications in two directions: (i) for improving our knowledge of confluent hypergeometric functions and Tricomi's function, (ii) and for engineering and physics problems.

• Ocean Systems Engineering
• /
• v.1 no.4
• /
• pp.317-336
• /
• 2011

A time domain finite element based method is employed to analyze wave radiation by multiple cylinders. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a finite element method (FEM). The matrix equation of the FEM is solved through an iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3D) mesh required is generated based on a two-dimensional (2D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary, respectively. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including two-cylinder and four-cylinder cases are provided to show the mutual influence due to cylinders on the first and second waves and forces.