• International Journal of Aeronautical and Space Sciences
• /
• 제2권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.

• 대한기계학회논문집A
• /
• 제27권9호
• /
• pp.1562-1570
• /
• 2003

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. The complexity of these constitutive equations generally requires a stable and accurate numerical method. The radial return mapping is one of the most robust integration scheme currently used. Nonlinear kinematic hardening model of Armstrong-Fredrick type has recovery term and the direction of kinematic hardening increment is not parallel to that of plastic strain increment. In this case, The conventional radial return mapping method cannot be applied directly. In this investigation, we expanded the radial return mapping method to consider the nonlinear kinematic hardening model and implemented this integration scheme into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using Newton method and bisection method. Using dynamic yield condition derived from linearization of flow rule, the integration scheme for elastoplastic and viscoplastic constitutive model was unified. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• International Journal of Aeronautical and Space Sciences
• /
• 제12권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.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2004년도 제22회 춘계학술대회논문집
• /
• pp.426-432
• /
• 2004

Large-Eddy Simulation (LES) is applied for the simulation of compressible flat plate boundary with Reynolds number up to 5 X 10$^{5}$ . Numerical examples include shock/boundary layer interaction and boundary layer transition, aiming future application to the analysis of transonic fan/compressor cascades. The present LES code uses hybrid com-pact/WENO scheme for the spatial discretization and compact diagonalized implicit scheme for the time integration. The present code successfully predicted the bypass transition of subsonic boundary layer. As for supersonic turbulent boundary layer, mean and fluctuation velocity of the attached boundary, as well as the evolution of the friction coefficient and the displacement thickness both upstream and downstream of the separation region are all in good agreement with experiment. The separation point also agreed with the experiment. In the simulation of the shock/laminar boundary layer interaction, the dependence of the transition upon the shock strength is reproduced qualitatively, but the extent of the separation region is overpredicted. These numerical examples show that LES can predict the behavior of boundary layer including transition and shock interaction, which are hardly managed by the conventional Reynolds-averaged Navier-Stokes approach, although there needs to be more effort before achieving quantitative agreement.

• 대한조선학회논문집
• /
• 제38권3호
• /
• pp.1-13
• /
• 2001

사항상태의 선박 주위의 점성유동 계산을 위하여 RANS 방정식에 대한 수치계산법을 개발하였다. 수치계산법은 이산화방법으로 유한체적법에 기초하여, 비점성 대류항에 대하여는 3차 정도의 flux-difference splitting 방법을 사용하고 시간 적분은 Euler 음해법을 사용하였다. 난류모형으로는 Spalart-Allmaras one-equation 모형을 사용하였다. 개발된 수치계산법을 이용하여 선수형상은 같으나 선미형상이 다른 두 VLCC 선형에 대한 조종유체력 및 유동 특성을 계산하고 이를 실험결과와 함께 비교하고 살펴보았다. 계산결과는 구속모형시험과 국부유동계측으로부터 얻은 유체력 및 유동을 잘 예측하고 있을 뿐 아니라 선미형상 차이에서 나타나는 유체력 및 유동특성의 차이도 잘 보여주었다.

• Coupled systems mechanics
• /
• 제8권2호
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Advances in aircraft and spacecraft science
• /
• 제9권5호
• /
• pp.449-462
• /
• 2022

Space-telescopes placed in the Sun-Earth second Lagrange point (L2) observe the sky following a scan strategy that is usually based on a spin-precession motion. Knowing which regions of the sky will be more observed by the instrument is important for the science operations and the instrument calibration. Computing sky observation parameters numerically (discretizing time and the sky) can consume large amounts of time and computational resources, especially when high resolution isrequired.This problem becomesmore critical if quantities are evaluated at detector level instead of considering the instrument entire Field of View (FoV). In previous studies, the authors have derived analytic solutions for quantities that characterize the observation of each point in the sky in terms of observation time according to the scan strategy parameters and the instrument FoV. Analytic solutions allow to obtain results faster than using numerical methods as well as capture detailed characteristics which can be overseen due to discretization limitations. The original approach is based on the analytic expression of the instrument trace over the sky. Such equations are implicit and thusrequiresthe use of numeric solversto compute the quantities.In this work, a new and simpler approach for computing one ofsuch quantities(mean observation time) is presented.The quantity is first computed for pure spin motion and then the effect of the spin axis precession is incorporated under the assumption that the precession motion is slow compared to the spin motion.In this sense, this new approach further simplifies the analytic approach, sparing the use of numeric solvers, which reduces the complexity of the implementation and the computing time.

• 대한기계학회논문집A
• /
• 제24권5호
• /
• pp.1075-1083
• /
• 2000

The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method. The prediction of the inelastic behavior of the structure is an essential part of reliability assessment procedure, because most of the failures are induced by the inelastic deformation, such as creep and plastic deformation. During decades, there has been much progress in understanding of the inelastic behavior of the materials and a lot of inelastic constitutive equations have been developed. These equations consist of the definition of inelastic strain and the evolution of the state variables introduced to quantify the irreversible processes occurred in the material. With respect to the definition of the inelastic strain, the inelastic constitutive models can be categorized into elastoplastic model, unified viscoplastic model and separated viscoplastic model and the different integration methods have been applied to each category. In the present investigation, the generalized integration method applicable for various types of constitutive equations is developed and implemented into ABAQUS by means of UMAT subroutine. The solution of the non-linear system of algebraic equations arising from time discretization with the generalized midpoint rule is determined using line-search technique in combination with Newton method. The strategy to control the time increment for the improvement of the accuracy of the numerical integration is proposed. Several numerical examples are considered to demonstrate the efficiency and applicability of the present method.

• 한국정보과학회논문지:소프트웨어및응용
• /
• 제33권4호
• /
• pp.365-377
• /
• 2006

단백질들은 서로 다른 단백질들과 상호작용 하거나 복합물을 형성함으로써 생물학적으로 중요한 기능을 한다고 알려져 있다. 때문에 대부분의 세포작용에 있어 중요한 역할을 하는 단백질 상호작용의 분석 및 예측에 대한 연구는 여러 연구그룹으로부터 풍부한 데이타가 산출되고 있는 현(現) 게놈시대에서 또 하나의 중요한 이슈가 되고 있다. 본 논문에서는 효모(Saccharomyces cerevisiae)에 대해 공개되어있는 단백질 상호작용 데이타들에서 속성들 간의 연관을 통해 유추 가능한 잠재적 단백질 상호작용들을 예측하기 위한 연관속성 마이닝 방법을 제시한다. 단백질의 속성들 중 연속값을 가지는 속성값들은 최대상호 의존성에 기반을 두어 이산화 하였으며, 정보이론기반 속성선택 알고리즘을 사용하여 단백질들 간의 상호작용 예측을 위해 고려되는 단백질의 속성(attribute) 수 증가에 따른 속성차원문제를 극복하도록 하였다. 속성들 간의 연관성 발견은 데이타마이닝 분야에서 사용되는 연관규칙 발견(association rule discovery) 방법을 사용하였다 논문에서 제안한 방법은 발견된 연관규칙을 통한 단백질 상호작용 예측문제에 있어 최대 약 96.5%의 예측 정확도를 보였으며 속성필터링을 통하여 속성필터링을 하지 않는 기존의 방법에 비해 최대 약 29.4% 연관규칙 발견속도 향상을 보였다.

• 한국항공우주학회지
• /
• 제36권3호
• /
• pp.211-219
• /
• 2008

본 연구에서는 삼차원 비정렬 비점성 유동 해석 코드를 이용하여 전투기 형상의 외부 장착물이 꼬리 날개에 미치는 공력 간섭효과에 대한 연구를 수행하였다. 수치적 기법으로는 격자점 중심(vertex-centered)에 기초한 유한체적법과 시간 적분을 위한 내재적인 point Gauss-Seidel 반복 계산법을 사용하였다. 해석 코드의 검증을 위해 전투기 형상에 대한 정상 유동 계산을 수행하였으며 결과를 실험 데이터와 비교하였다. 외부 장착물의 후류(wake)를 정확히 포착하기 위해 예상되는 후류 영역에 대해 국부적인 격자 조밀화를 수행하였으며 천음속 영역에서의 비정상 유동 해석을 통해 외부 장착물에서 발생하는 후류가 수평꼬리날개에 미치는 간섭효과를 확인하였다. 공력 간섭효과를 감소시키기 위한 대안으로는 뒷전 플랩 꺽임각(trailing edge flap deflection)을 고려하였으며 이에 대한 정량적인 감소효과를 제시하였다.