• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.30 no.8
• /
• pp.37-45
• /
• 2002

A two-dimensional unsteady, inviscid flow solver has been developed for the simulation of airfoil-vortex interactions on unstructured dynamically adapted meshes. The Euler solver is based on a second-order accurate implicit time integration using a point Gauss-Seidel relaxation scheme and a dual time-step subiteration. A vertex-centered, finite-volume discretization is used in conjunction with the Roe's flux-difference splitting. An unsteady solution-adaptive dynamic mesh scheme is used by adding and deleting mesh points to take account of both spatial and temporal variations of the flow field. The effect of vortex interaction on the aeroelastic displacement of an airfoil attached to the idealized two degree-of-freedom spring system is investigated.

• International Journal of Aeronautical and Space Sciences
• /
• v.2 no.2
• /
• pp.65-72
• /
• 2001

This paper demonstrates an explicit-implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for an existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory. Darcy's law is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full Newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.8
• /
• pp.8-17
• /
• 2005

It is well known that preconditioning methods are efficient for convergence acceleration in the compressible low Mach number flows. In this study, the original Euler equations and three differently nondimensionalized preconditioning methods are implemented in two dimensional inviscid bump flows using the 3rd order MUSCL and DADI schemes as numerical flux discretization and time integration, respectively. The multigrid and local time stepping methods are also used to accelerate the convergence. The test case indicates that a properly modified local preconditioning technique involving concepts of a global preconditioning allows Mach number independent convergence. Besides, an asymptotic analysis for properties of preconditioning methods is added.

• International Journal of Aeronautical and Space Sciences
• /
• v.12 no.4
• /
• pp.318-330
• /
• 2011

This paper investigates the convergence characteristics of the modified artificial compressibility method proposed by Turkel. In particular, a focus is mode on the convergence characteristics due to variation of the preconditioning factor (${\alpha}_u$) and the artificial compressibility (${\beta}$) in conjunction with an upwind method. For the investigations, a code using the modified artificial compressibility is developed. The code solves the axisymmetric incompressible Reynolds averaged Navier-Stokes equations. The cell-centered finite volume method is used in conjunction with Roe's approximate Riemann solver for the inviscid flux, and the central difference discretization is used for the viscous flux. Time marching is accomplished by the approximated factorization-alternate direction implicit method. In addition, Menter's k-${\omega}$ shear stress transport turbulence model is adopted for analysis of turbulent flows. Inviscid, laminar, and turbulent flows are solved to investigate the accuracy of solutions and convergence behavior in the modified artificial compressibility method. The possible reason for loss of robustness of the modified artificial compressibility method with ${\alpha}_u$ >1.0 is given.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.50 no.11
• /
• pp.771-781
• /
• 2022

Sequential convex programming based performance analysis of the designed UAV is performed. The nonlinear optimization problems generated by aerodynamics are approximated to socond order program by discretization and convexification. To improve the performance of the algorithm, the solution of the relaxed problem is used as the initial trajectory. Dive trajectory optimization problem is analyzed through iterative solution procedure of approximated problem. Finally, the maximum final velocity according to the performance of the actuator model was compared.

• Communications of the Korean Mathematical Society
• /
• v.20 no.3
• /
• pp.563-578
• /
• 2005

In [Z. Cai and B. C. Shin, SIAM J. Numer. Anal. 40 (2002), 307-318], we developed the discrete first-order system least squares method for the second-order elliptic boundary value problem by directly approximating $H(div){\cap}H(curl)-type$ space based on the Helmholtz decomposition. Under general assumptions, error estimates were established in the $L^2\;and\;H^1$ norms for the vector and scalar variables, respectively. Such error estimates are optimal with respect to the required regularity of the solution. In this paper, we study solution methods for solving the system of linear equations arising from the discretization of variational formulation which possesses discrete biharmonic term and focus on numerical results including the performances of multigrid preconditioners and the finite element accuracy.

• Journal of applied mathematics & informatics
• /
• v.8 no.1
• /
• pp.9-26
• /
• 2001

In some previous papers the author extended two algorithms proposed by Z. Kovarik for approximate orthogonalization of a finite set of linearly independent vectors from a Hibert space, to the case when the vectors are rows (not necessary linearly independent) of an arbitrary rectangular matrix. In this paper we describe combinations between these two methods and the classical Kaczmarz’s iteration. We prove that, in the case of a consistent least-squares problem, the new algorithms so obtained converge ti any of its solutions (depending on the initial approximation). The numerical experiments described in the last section of the paper on a problem obtained after the discretization of a first kind integral equation ilustrate the fast convergence of the new algorithms. AMS Mathematics Subject Classification : 65F10, 65F20.

• Bulletin of the Korean Mathematical Society
• /
• v.51 no.3
• /
• pp.717-737
• /
• 2014

We study the effect of numerical quadrature in space on semidiscrete and fully discrete piecewise linear finite element methods for parabolic interface problems. Optimal $L^2(L^2)$ and $L^2(H^1)$ error estimates are shown to hold for semidiscrete problem under suitable regularity of the true solution in whole domain. Further, fully discrete scheme based on backward Euler method has also analyzed and optimal $L^2(L^2)$ norm error estimate is established. The error estimates are obtained for fitted finite element discretization based on straight interface triangles.

• Journal of the Korean Mathematical Society
• /
• v.54 no.1
• /
• pp.193-213
• /
• 2017

We investigate an efficient approximation of solution to stochastic Burgers equation driven by an additive space-time noise. We discuss existence and uniqueness of a solution through the Orstein-Uhlenbeck (OU) process. To approximate the OU process, we introduce the Karhunen-$Lo{\grave{e}}ve$ expansion, and sparse grid stochastic collocation method. About spatial discretization of Burgers equation, two separate finite element approximations are presented: the conventional Galerkin method and Galerkin-conservation method. Numerical experiments are provided to demonstrate the efficacy of schemes mentioned above.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.1-28
• /
• 2004

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by an $H^1$-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) of index one. A priori error estimates for semidiscrete scheme are derived for both differ-ential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.