• Korea-Australia Rheology Journal
• /
• v.19 no.4
• /
• pp.233-242
• /
• 2007

We performed a numerical investigation to find out the optimal choice of the spatial discretization in the distributed-Lagrangian-multiplier/fictitious-domain (DLM/FD) method for the solid/fluid interaction problem. The elastic solid bar attached on the bottom in a pressure-driven channel flow of a Newtonian fluid was selected as a model problem. Our formulation is based on the scheme of Yu (2005) for the interaction between flexible bodies and fluid. A fixed regular rectangular discretization was applied for the description of solid and fluid domain by using the fictitious domain concept. The hydrodynamic interaction between solid and fluid was treated implicitly by the distributed Lagrangian multiplier method. Considering a simplified problem of the Stokes flow and the linearized elasticity, two numerical factors were investigated to clarify their effects and to find the optimum condition: the distribution of Lagrangian multipliers and the solid/fluid interfacial condition. The robustness of this method was verified through the mesh convergence and a pseudo-time step test. We found that the fluid stress in a fictitious solid domain can be neglected and that the Lagrangian multipliers are better to be applied on the entire solid domain. These results will be used to extend our study to systems of elastic particle in the Stokes flow, and of particles in the viscoelastic fluid.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.189-194
• /
• 2001

In this paper, non-reacting and reacting flowfields were computed using a preconditioned Navier-Stokes solver. The preconditioning technique of Merkle et al. and TVD scheme or Chakravarthy and Osher was employed and the results obtained using developed code have a good agreement with the previous results and experimental data. The preconditioned Wavier-Stokes equation set with low Reynolds number $\kappa-\epsilon$ equation and species continuity equations, are discretized with strongly implicit manner and time integrated with LU-SSOR scheme. For the purpose of treating unsteady problem the duel-time stepping scheme was employed. For the validation of the code in incompressible flow regime, steady driven square cavity flow was considered and calculation result shows reasonably good agreement with the result of incompressible code. Shock wave/boundary layer interaction problem was considered to show the shock capturing performance of preconditioned-TVD scheme. To validate unsteady flow, acoustic oscillation problem was calculated, and supersonic premix flame of $H_2$-air reaction problem which is calculated with turbulence model, 9-species/18-reaction step reaction model, shows reasonable agreement with the previous results. As a result, the preconditioning method has an advantage to calculate incompressible and compressible flow through one code and preconditioned solver easily developed from standard compressible code with minor efforts. But additional computational time and computer memory is required due to preconditioning matrix.

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

• Journal of the Korean Society for Industrial and Applied Mathematics
• /
• v.14 no.2
• /
• pp.125-140
• /
• 2010

This paper develops a least-squares approach to the solution of the optimal control problem for the Navier-Stokes equations. We recast the optimality system as a first-order system by introducing velocity-flux variables and associated curl and trace equations. We show that a least-squares principle based on $L^2$ norms applied to this system yields optimal discretization error estimates in the $H^1$ norm in each variable.

• 한국전산유체공학회:학술대회논문집
• /
• 2001.05a
• /
• pp.149-154
• /
• 2001

Parallel performance of a Myrinet based PC-cluster was tested and compared with a conventional Fast-Ethernet system. A preconditioned Navier-Stokes code was parallelized with domain decomposition technique, and used for the parallel performance test. Speed-up ratio was examined as a major performance parameter depending on the number of processor and the network topology. As was expected, Myrinet system shows a superior parallel performance to the Fast-Ethernet system even with a single network adpater for a dual processor SMP machine. A test for the dependency on problem size also shows that network communication speed is a crucial factor for parallelized computational fluid dynamics analysis and the Myrinet system is a plausible candidate for high performance parallel computing system.

• Proceedings of the KSME Conference
• /
• 2002.08a
• /
• pp.229-232
• /
• 2002

For the purpose of Thermal Protection Material design problem, a numerical analysis of axisymmetric high temperature supersonic impinging jet flows of exhaust gas from combustor on curved surfaces has been accomplished. A modifed CSCM Upwind Navier-Stokes method which is able to cure the carbuncle Phenomena has been developed to study strong shock wave structure and thermodynamic wall properties such as pressure and heat transfer rate on various curved surfaces. The results show that the maximum heat transfer rate which is the most important parameter affecting thermo-chemical surface ablation on the plate did not occur at the center of jet impingement, but rather on a circle slightly away from the center of impingement and the shear stress distribution alone the wall is similar to the wall heat transfer late distribution.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.14 no.1
• /
• pp.11-19
• /
• 2010

An approach to reduce cancellation errors of the temperature preconditioned Navier-Stokes equations is proposed. This approach is also applied to the conventional preconditioning methods. Adiabatic laminar viscous flows around a circular cylinder are calculated at different Mach numbers. It is shown that a redefinition of total enthalpy for reducing magnitude of the enthalpy remarkably alleviates cancellation problems of the temperature preconditioning.

• Journal of the Korean Mathematical Society
• /
• v.57 no.4
• /
• pp.915-933
• /
• 2020

A two-level finite element method based on the Newton iterative method is proposed for solving the Navier-Stokes/Darcy model. The algorithm solves a nonlinear system on a coarse mesh H and two linearized problems of different loads on a fine mesh h = O(H^4-𝜖). Compared with the common two-grid finite element methods for the considered problem, the presented two-level method allows for larger scaling between the coarse and fine meshes. Moreover, we prove the stability and convergence of the considered two-level method. Finally, we provide numerical experiment to exhibit the effectiveness of the presented method.

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

Time marching algorithms applied to compressible Navier-Stokes equation have a convergence problem at low Mach number. It is mainly due to the eigenvalue stiffness and pressure singularity as Mach number approaches to zero. Among the several methods to overcome the shortcomings of time marching scheme, time derivative preconditioning method have been used successfully. In this numerical analysis, we adopted a preconditioner of K.H. Chen and developed a two-dimensional, axisymmetric Navier-Stokes program. The steady state driven cavity flow and backward facing step flow problems were computed to confirm the accuracy and the robustness of preconditioned algorithm for low Mach number flows. And the transonic and supersonic flows insice the JPL axisymmetric nozzle internal flow is exampled to investigate the effects of preconditioning at high Mach number flow regime. Test results showed excellent agreement with the experimental data.

• Journal of the Korean Mathematical Society
• /
• v.46 no.5
• /
• pp.1007-1025
• /
• 2009

First-order least-squares method of a distributed optimal control problem for the incompressible Stokes equations is considered. An optimality system for the optimal solution are reformulated to the equivalent first-order system by introducing the vorticity and then the least-squares functional corresponding to the system is defined in terms of the sum of the squared $H^{-1}$ and $L^2$ norms of the residual equations of the system. Finite element approximations are studied and optimal error estimates are obtained. Resulting linear system of the optimality system is symmetric and positive definite. The V-cycle multigrid method is applied to the system to test computational efficiency.