• 제목/요약/키워드: Implicit euler

검색결과 74건 처리시간 0.023초

A Study on the Nozzle-Rotor Interactions of Partial Admission Supersonic Turbines

  • Seong, Young-sik;Han, Seong-hoon;Kim, Kui-soon;Park, Chang-kyu
    • 한국추진공학회:학술대회논문집
    • /
    • 한국추진공학회 2004년도 제22회 춘계학술대회논문집
    • /
    • pp.391-397
    • /
    • 2004
  • The performance characteristics of partial admission supersonic turbines are analyzed by using the commercial CFD program FLUENT6.0. The governing equations were discretized with Euler implicit method in time and 2nd-order upwind scheme of FVM in space. The k-$\varepsilon$ turbulence model was utilized to describe the turbulent flow field. In order to investigate the nozzle--rotor interactions and the effect of partial admission, the flows in supersonic turbine rotor cascades with a nozzle are computed. Extensive computations of partial admission supersonic turbines provide the shock structures and flow patterns in the nozzle and rotor. It is clearly shown that the nozzle flow is highly affected by the shocks or expansion waves propagated from the rotor leading edge. And the rotor flow is also affected by the shocks or wakes originated from the nozzle.

  • PDF

Time-dependent simplified spherical harmonics formulations for a nuclear reactor system

  • Carreno, A.;Vidal-Ferrandiz, A.;Ginestar, D.;Verdu, G.
    • Nuclear Engineering and Technology
    • /
    • 제53권12호
    • /
    • pp.3861-3878
    • /
    • 2021
  • The steady-state simplified spherical harmonics equations (SPN equations) are a higher order approximation to the neutron transport equations than the neutron diffusion equation that also have reasonable computational demands. This work extends these results for the analysis of transients by comparing of two formulations of time-dependent SPN equations considering different treatments for the time derivatives of the field moments. The first is the full system of equations and the second is a diffusive approximation of these equations that neglects the time derivatives of the odd moments. The spatial discretization of these methodologies is made by using a high order finite element method. For the time discretization, a semi-implicit Euler method is used. Numerical results show that the diffusive formulation for the time-dependent simplified spherical harmonics equations does not present a relevant loss of accuracy while being more computationally efficient than the full system.

COMPUTATIONAL METHOD FOR SINGULARLY PERTURBED PARABOLIC REACTION-DIFFUSION EQUATIONS WITH ROBIN BOUNDARY CONDITIONS

  • GELU, FASIKA WONDIMU;DURESSA, GEMECHIS FILE
    • Journal of applied mathematics & informatics
    • /
    • 제40권1_2호
    • /
    • pp.25-45
    • /
    • 2022
  • In this study, the non-standard finite difference method for the numerical solution of singularly perturbed parabolic reaction-diffusion subject to Robin boundary conditions has presented. To discretize temporal and spatial variables, we use the implicit Euler and non-standard finite difference method on a uniform mesh, respectively. We proved that the proposed scheme shows uniform convergence in time with first-order and in space with second-order irrespective of the perturbation parameter. We compute three numerical examples to confirm the theoretical findings.

부분입사형 초음속 터빈의 노즐과 익렬의 상호작용에 관한 연구 (A Study on the Nozzle-Rotor Interactions of Partial Admission Supersonic Turbines)

  • 성영식;한성훈;김귀순;박창규
    • 한국추진공학회지
    • /
    • 제8권3호
    • /
    • pp.53-60
    • /
    • 2004
  • 노즐과 익렬의 상호작용과 부분입사 효과를 연구하기 위해 노즐과 터빈의 익렬을 같이 계산하였고, 충격파의 구조와 유동의 현상을 밝혀내기 위해 많은 계산을 하였다. 지배방정식은 시간에 관해서는 오일러 음해법으로, 공간에 관해서는 유한체적법을 사용하여 2차의 상류차분법으로 이산화하였으며 $\kappa$-$\varepsilon$ 난류 모델을 사용하였다. 본 논문에서는 익렬의 앞전에서 발생한 충격파와 확장파가 노즐 유동에 영향을 주는 것과 익렬내부의 유동이 노즐에서 발생한 충격파의 영향을 받는 것을 확인하였다.

Numerical result of complex quick time behavior of viscoelastic fluids in flow domains with traction boundaries

  • Kwon, Young-Don
    • Korea-Australia Rheology Journal
    • /
    • 제19권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.

1단 축류 터빈의 비정상 내부유동특성에 관한 2차원 해석(I) (Two-Dimensional Analysis of Unsteady Flow through One Stage of Axial Turbine (I))

  • 박준염;엄인식;백제현
    • 대한기계학회논문집B
    • /
    • 제23권11호
    • /
    • pp.1371-1378
    • /
    • 1999
  • Flow through turbomachinery has a very complex structure and Is Intrinsically unsteady. In addition, trend to highly loaded turbomachinery makes the flow extremely complex due to the interaction between rotor and stator. In this study, flows through UTRC LSRR turbine are numerically analyzed using 2 dimensional Navier-Stokes equations. The convective terms of the governing equations are discretized using the Van-Leer's FVS(Flux vector splitting) with an upwind TVD scheme. The conventional central differencing is used to discretize the diffusion terms on the finite volume. The accurate unsteady motion is achieved by using a 2nd order accurate, 3-point Euler implicit scheme. The quasi-conservative zonal scheme is used for calculating the flow variables on the zonal interface between the rotor and stator. The axial gap between stator and rotor has been configured in two variations, 15% and 65% of average chord length. The analysis program is validated using experimental results and the effect of axial gap is examined. The numerical analysis results are presented by time averaged pressure coefficient and pressure magnitude coefficient and compared with experimental results.

비정렬 동적격자를 이용한 블레이드-와류 간섭에 따른 공탄성 변위예측 (Prediction of Aeroelastic Displacement Under Close BVI Using Unstructured Dynamic Meshes)

  • 조규원;오우섭;권오준;이인
    • 한국항공우주학회지
    • /
    • 제30권8호
    • /
    • pp.37-45
    • /
    • 2002
  • 본 연구에서는 이차원에서 비정상 비점성 유동해석을 위한 비정렬 동적 편자 기법을 개발하였다. 유동해석 기법은 시간에 대해 2차의 정확도를 갖는 내재적인 시간적분법을 사용하였으며, 격자중심의 유한 체적법과 Roe의 풍상차분법을 이용하여 공간에 대한 차분화를 하였다. 시간과 공간에 대한 정확도를 증가시키기 위해서는 해에 따라 원하는 위치에 격자점들을 임의로 추가할 수 있는 비정상 동적 적응격자 기법을 사용하였다. 이를 이용하여 이차원의 2자유도를 갖는 스프링 에어포일 시스템의 와류와의 간섭현상에 따른 공탄성적 변위를 예측하였다.

다결정재 소성변형의 탄소성 해석을 위한 접선강성 개발 (A New Tangent Stiffness for Anisotropic Elasto-Viscoplastic Analysis of Polycrystalline Deformations)

  • 윤종헌;허훈;이용신
    • 한국소성가공학회:학술대회논문집
    • /
    • 한국소성가공학회 2006년도 춘계학술대회 논문집
    • /
    • pp.349-352
    • /
    • 2006
  • The plastic deformation of polycrystalline materials is induced by changes of the microstructure when the loading is beyond the critical state of stress. Constitutive models for the crystal plasticity have the common objective which relates microscopic single crystals in the crystallographic texture to the macroscopic continuum point. In this paper, a new consistent tangent stiffness for the anisotropic elasto-viscoplastic analysis of polycrystalline deformation is developed, which can be used in the finite element analysis for the slip-dominated large deformation of polycrystalline materials. In order to calculate the consistent tangent stiffness, the state function is defined based on the consistency condition between the elastic and plastic stress. The rate of shearing increment($\Delta{\gamma}^{\alpha}$) is calculated with satisfying the consistency condition. The consistency condition becomes zero when the trial resolved shear stress($\tau^{{\alpha}^*}$) becomes resolved shear stress($\tau^{\alpha}$) at every step. Iterative method is utilized to calculate the rate of shearing increment based on the implicit backward Euler method. The consistent tangent stiffness can be formulated by differentiating the rate of shearing increment with total strain increment after the instant rate of shearing increment converges. The proposed tangent stiffness is applied to the ABAQUS/Standard by implementing in the ABAQUS/UMAT.

  • PDF

점성 및 충격파효과를 고려한 천음속 터빈 케스케이드의 유체유발 진동해석 (Flow-induced Vibration of Transonic Turbine Cascades Considering Viscosity and Shock Wave Effects)

  • 오세원;박웅;김동현
    • 한국소음진동공학회논문집
    • /
    • 제16권9호
    • /
    • pp.937-948
    • /
    • 2006
  • In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

점성 및 충격파 효과를 고려한 천음속 터빈 케스케이드의 유체유발 진동해석 (Flow-Induced Vibration of Transonic Turbine Cascades Considering Viscosity and Shock Wave Effects)

  • 오세원;김동현;박웅
    • 한국소음진동공학회:학술대회논문집
    • /
    • 한국소음진동공학회 2006년도 춘계학술대회논문집
    • /
    • pp.793-802
    • /
    • 2006
  • In this study, a fluid/structure coupled analysis system for simulating complex flow-induced vibration (FIV) phenomenon of cascades has been developed. The flow is modeled using Euler and Wavier-Stokes equations with different turbulent models. The fluid domains are modeled using the unstructured grid system with dynamic deformations due to the motion of structural boundary. The Spalart-Allmaras (S-A) and the SST ${\kappa}-{\omega}$ turbulent models are used to predict the transonic turbulent flows. A fully implicit time marching scheme based on the Newmark direct integration method is used in order to solve the coupled governing equations for viscous flow-induced vibration phenomena. For the purpose of validation for the developed FIV analysis system, comparison results for computational analyses of steady and unsteady aerodynamics and flutter analyses are presented in the transonic flow region. In addition, flow-induced vibration analyses for the isolated cascade and multi-blades cascade models have been conducted to show the physical fluid-structure interaction effects in the time domain.

  • PDF