• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.11a
• /
• pp.938-941
• /
• 2004

Sound generation and radiation from the duct-rotor system are calculated numerically. The wake geometries of a two-bladed rotor are calculated by using a time-marching fiee-wake method without a non-physical model of the far wake. Acoustic free field due to a rotating rotor is obtained by Lowson's equation. Using Kirchhoff source, rotating sources are modeled as stationary ones and can be inserted in the thin body boundary element method. The Kirchhoff source is validated through calculation of acoustic pressure due to a rotating point force. The thin body boundary element method (thin body BEM) is validated through calculation of acoustic radiation of ducted dipole. Using Kirchhoff source and thin body BEM, acoustic radiation of a ducted rotating source is calculated. Acoustic shielding is observed by inserting a duct and shows different phenomena at each major frequency. Acoustic radiation of a real duct-rotor system is also calculated using this method and the ducted acoustic field is significantly different from rotor only.

• The KSFM Journal of Fluid Machinery
• /
• v.3 no.2 s.7
• /
• pp.15-23
• /
• 2000

The present work describes the prediction method for the unsteady flow field and the acoustic pressure field of a ducted axial fan. The prediction method is comprised of time-marching free-wake method, acoustic analogy, and the Kirchhoff-Helmholtz BEM. The predicted sound signal of a rotor is similar to the experiment one. We assume that the rotor rotates with a constant angular velocity and the flow field around the rotor is incompressible and inviscid. Then, a time-marching free-wake method is used to model the fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lowson's method is used to predict the acoustic source. The newly developed Helmholtz-Kirchhoff BEM lot thin body is used to calculate tile sound field of the ducted fan. The ducted fan with 6 blades is analysed and the sound field around the duct is calculated.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2000.11a
• /
• pp.100-105
• /
• 2000

The present work describes the prediction method for the unsteady flow field and the acoustic pressure field of a ducted axial fan. The prediction method is comprised of time-marching free-wake method, acoustic analogy. and the Helmholtz-Kirchhoff BEM. The predicted sound signal of a rotor is similar to the experiment one. We assume that the rotor rotates with a constant angular velocity and the flow field around the rotor is incompressible and inviscid. Then, a time-marching free-wake method is used to model the fan and to calculate the flow field. The force of each element on the blade is calculated by the unsteady Bernoulli equation. Lawson's method is used to predict the acoustic source. The newly developed Helmholtz-Kirchhoff BEM for thin body is used to calculate the sound field of the ducted fan. The ducted fan with 6 blades is analysed and the sound field around the duct is calculated.

• International Journal of Aeronautical and Space Sciences
• /
• v.10 no.2
• /
• pp.63-76
• /
• 2009

A study on the interference effects of overlapping tandem rotors in forward flight is conducted using the time-marching free-wake panel method which adopts field velocity boundary integral formulation. The conventional boundary integral formulation is numerically unstable for the cases when the blade and the wake are in close proximity to each other. In order to avoid this problem, this study applies the field velocity method and modifies the boundary integration formulation. The improved method is used for the parametric study on the advance ratio and the distance between the rotors. These are the parameters that most affect the interference of the tandem rotor in forward flight. Comparison of the aerodynamic performance shows that the horizontal distance between the rotors negligibly influences the overlap-induced power factor for high advance ratio. In addition, it shows that the overlap-induced power factor is inversely proportional to the squared vertical distance between the rotors, and that the overlap-induced power factor increases to a certain extent and decrease back as the advance ratio increases.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.5
• /
• pp.1244-1253
• /
• 1990

Time-Marching methods solving Euler equations are used for calculation of two-dimensional, steady, inviscid flow through a stationary compressor cascade. Calculation method is based on the Denton`s opposed difference scheme. A smoothing in the axial direction is used to increase the stability of solution. The computational grid consists of quadrilateral elements, one of which has four nodes at each corner and the grid points on the upper periodic boundaries are located one pitch away from those on the lower boundaries to satisfy the periodicity condition. Results of calculation show good agreement with other computational and experimental results, proving that the present method of calculation should be applied with confidence for the cascade flow with shock wave.

• Journal of Korea Game Society
• /
• v.20 no.3
• /
• pp.75-84
• /
• 2020

This paper introduces an interactive indirect illumination algorithm which considers indirect visibility. First, a small number of rays are emitted on hemisphere of the current pixel to obtain the first intersection. If this point is directly illuminated by the light source, its illuminated color is collected. Second, in order to approximate the indirect visibility, a 3D ray marching algorithm, which is based on a hierarchy structure, is used to accelerate the ray-voxel intersection. Third, the indirect images are denoised by an edge-avoiding filtering with a local means replacement method.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.19 no.4
• /
• pp.427-435
• /
• 2008

In this paper, a stable marching-on in time(MOT) method with a time domain combined field integral equation(CFIE) is presented to obtain the transient scattering response from arbitrarily shaped three-dimensional conducting bodies. This formulation is based on a linear combination of the time domain electric field integral equation(EFIE) with the magnetic field integral equation(MFIE). The time derivatives in the EFIE and MFIE are approximated using a central finite difference scheme and other terms are averaged over time. This time domain CFIE approach produces results that are accurate and stable when solving for transient scattering responses from conducting objects. Numerical results with the proposed MOT scheme are presented and compared with those obtained from the conventional method and the inverse discrete Fourier transform(IDFT) of the frequency domain CFIE solution.

• The Journal of Korean Institute of Electromagnetic Engineering and Science
• /
• v.22 no.10
• /
• pp.1005-1011
• /
• 2011

In this paper, we present a marching-on-in-degree(MOD) finite difference method(FDM) based on the Helmholtz wave equation for analyzing transient electromagnetic responses in a general dispersive media. The two issues related to the finite difference approximation of the time derivatives and the time consuming convolution operations are handled analytically using the properties of the Laguerre functions. The basic idea here is that we fit the transient nature of the fields, the flux densities, the permittivity with a finite sum of orthogonal Laguerre functions. Through this novel approach, not only the time variable can be decoupled analytically from the temporal variations but also the final computational form of the equations is transformed from finite difference time-domain(FDTD) to a finite difference formulation through a Galerkin testing. Representative numerical examples are presented for transient wave propagation in general Debye, Drude, and Lorentz dispersive medium.

• International Journal of Aeronautical and Space Sciences
• /
• v.6 no.2
• /
• pp.84-96
• /
• 2005

Helical tip vortex is known as stable vortex structure, however the specific frequency component of far wake perturbation induces the vortex pairing in hover and axial flight. It is expected that the tip vortex pairing phenomena may happen in transition flight and very low advance ratio flight so that inflow may be most nonuniform in the low advance ratio flight. The objectives of this paper are that a tip-vortex instability during the transition from hover into very low advance ratio forward flight is numerically predicted to understand a physics by using a time-marching free-wake method. To achieve the objectives, numerical method is firstly validated in typical axial and forward flights cases. Present scheme with trim routine can predict airloads and inflow distribution of forward flight with good accuracy. Then, the transition flight condition is calculated. The rotor used in this wake calculation is a small-scale AH-1G model. By using a tip-vortex trajectory tracking method, the tip-vortex pairing process are clearly observed in transient flight($\mu$=0.03) and disappears at a slightly higher advance ratio($\mu$=0.05). According to the steady flight simulation at $\mu$=0.03, it is confirmed the tip-vortex pairing process is continued in the rear part of rotor disk and not occurs in the front part. Time averaged inflow in this case is predicted as smooth distribution.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.16 no.1
• /
• pp.33-42
• /
• 2003

This paper deals with the FEM analysis of nonlinear sloshing of incompressible, invicid and irrotational flow in two dimensional rectangular tank. We use laplace equation based on potential theory as governing equation. For large amplitude sloshing motion, kinematic and dynamic free surface conditions derived from Bernoulli equation are applied. This problem is solved by FEM using 9-node elements. For the time integration and accurate velocity calculation, we introduce predictor-corrector time marching scheme and least square method. Also, numerical stability in tracking of free surface is obtained by direct calculation of free surface location to time variation. Numerical results of sloshing induced by harmonic excitations, while comparing with those of linear theory and references, prove the accuracy and stability. After verification of our program, we analyze sloshing response characteristics to the fluid height and the excitation amplitude.