• International Journal of Aeronautical and Space Sciences
• /
• v.7 no.2
• /
• pp.155-174
• /
• 2006

Parabolized stability equations for compressible flows in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Compressible and incompressible flat plate flow stability under two-dimensional and three¬dimensional disturbances has been investigated to test the present code. Results of the present computation are found to be in good agreement with the multiple scale analysis and DNS data. Stability calculation results by the present PSE code for compressible boundary layer at Mach numbers ranging from 0.02 to 1.5 are also presented and are again seen to be as accurate as the spectral method.

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

An accurate and cost efficient method PSE is used for the stability analysis of 2D or 3D compressible boundary layers. A highly accurate finite difference PSE code has been developed at a general curvilinear coordinate system using an implicit marching procedure to deal with a broad range of transition predictions problems. Evolution of disturbances in compressible flat plate boundary layers are studied for free-stream Mach numbers ranging from 0 to 1.5. The effect of mean-flow nonparallelism is found to be weak on two dimensional waves and strong on three dimensional waves. The maximum amplification rate increases monotonically with Mach number. The present PSE solutions are compared with previous numerical investigations and experimental results and are found to be in good agreement.

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

Nonlinear parabolized stability equations for compressible flow in general curvilinear coordinate system are derived to deal with a broad range of transition prediction problems on complex geometry. A highly accurate finite difference PSE code has been developed using an implicit marching procedure. Blasius flow is tested. The results of the present computation show good agreement with DNS data. Nonlinear interaction can make the T-S fundamental wave more unstable and the onset of its amplitude decay is shifted downstream relative to linear case. For nonlinear calculations, rather small difference in initial amplitude can produce large change during nonlinear region. Compressible secondary instability at Mach number 1.6 is also simulated and showed that 1.1% initial amplitude for primary mode is enough to trigger the secondary growth.

• 한국전산유체공학회:학술대회논문집
• /
• 1998.05a
• /
• pp.110-115
• /
• 1998

An axisymmetric Navier-Stokes procedure has been developed to analyze the pulse jet flow in a ceramic filter unit for the dust dislodging process. Using Baldwin-Lomax turbulence model as a closure relationship, the SIAF(Scalar Implicit Approximate Factorization) algorithm together with the ${\delta}^k-Correction$ iterative time marching scheme is adopted to solve the unsteady compressible Navier-Stokes equations. After some validation tests, the code has been applied to solve the pulse jet flow and examine the effects of geometry and reservoir pressure condition on the pressure level inside the filter unit. To avoid dealing with the uncertainty of such factors as the cohesion of the collected dust and the adhesion of the dust to the medium and also to simplify the analysis, the filter wall is assumed to be impermeable. The results for various test cases are presented.

• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.211-219
• /
• 2007

Here we demonstrate complex transient behavior of viscoelastic liquid described numerically with the Leonov model in straight and contraction channel flow domains. Finite element and implicit Euler time integration methods are employed for spatial discretization and time marching. In order to stabilize the computational procedure, the tensor-logarithmic formulation of the constitutive equation with SUPG and DEVSS algorithms is implemented. For completeness of numerical formulation, the so called traction boundaries are assigned for flow inlet and outlet boundaries. At the inlet, finite traction force in the flow direction with stress free condition is allocated whereas the traction free boundary is assigned at the outlet. The numerical result has illustrated severe forward-backward fluctuations of overall flow rate in inertial straight channel flow ultimately followed by steady state of forward flow. When the flow reversal occurs, the flow patterns exhibit quite complicated time variation of streamlines. In the inertialess flow, it takes much more time to reach the steady state in the contraction flow than in the straight pipe flow. Even in the inertialess case during startup contraction flow, quite distinctly altering flow patterns with the lapse of time have been observed such as appearing and vanishing of lip vortices, coexistence of multiple vortices at the contraction comer and their merging into one.