• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2001년도 추계 학술대회논문집
• /
• pp.84-90
• /
• 2001

Secondary instability of flow past a circular cylinder is examined using direct numerical simulation at Reynolds number 220 and 250. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme, and the pressure and viscosity terms are applied by the 4th order central difference. In spanwise, Navier-Stokes equation is distributed using Spectral Method. The critical Reynolds number for this instability is found to be about Re=190. The secondary instability leads re three-dimensionality with a spanwise wavelength about 4 cylinder diameters at onset (A-mode). Results of three-dimensional effect in wake of a circular cylinder are represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• 한국소성가공학회:학술대회논문집
• /
• 한국소성가공학회 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.

• 대한기계학회:학술대회논문집
• /
• 대한기계학회 2001년도 추계학술대회논문집B
• /
• pp.570-577
• /
• 2001

Three-dimensional time-dependent flow past a circular cylinder is numerically investigated using direct numerical simulation for Reynolds number 280 and 300. The higher-order finite difference scheme is employed for the spatial distributions along with the second order Adams-Bashforth and the first order backward-Euler time integration. In x-y plane, the convection term is applied by the 5th order upwind scheme and the pressure and viscosity terms are applied by the 4th order central difference. And in spanwise, Navier-Stokes equation is distributed using of Spectral Method. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength about 1.0 diameters at onset (B-mode). Present results of three-dimensional effects of in wake of a circular cylinder is represented with spanwise and streamwise vorticity contours as Reynolds numbers.

• 한국추진공학회:학술대회논문집
• /
• 한국추진공학회 2008년 영문 학술대회
• /
• pp.723-728
• /
• 2008

We performed numerical analysis of base drag phenomenon, when a projectile with backward step flies into atmosphere at supersonic speed. We compared with other researchers. From our previous studies that were 2-dimensional simulation, we found out from sophisticated simulations that need dense mesh points to compare base pressure and velocity profile after from base with experimental data. Therefore, we focus on high order spatial disceretization over 3rd order with TVD such as MUSCL TVD 3rd, 5th, and WENO 5th order, and Limiters such as minmod, Triad. Moreover, we enforce to flux averaging schemes such as Roe, RoeM, HLLE, AUSMDV. In present, one dimensional result of Euler tests, there are Sod, Lax, Shu-Osher and interacting blast wave problems. AUSMDV as a flux averaging scheme with MUSCL TVD 5th order as spatial resolution is good agreement with exact solutions than other combinations. We are carrying out the same approaches into 3-dimensional base flow only candidate flux schemes that are Roe, AUSMDV. Additionally, turbulence models are used in 3-dimensional flow, one is Menter s SST DES model and another is Sparlat-Allmaras DES/DDES model in Navier-Stokes equations.

• 천문학회지
• /
• 제56권2호
• /
• pp.287-292
• /
• 2023

This paper investigates the number of scatterings a photon undergoes in random walks before escaping from a medium. The number of scatterings in random walk processes is commonly approximated as τ + τ² in the literature, where τ is the optical thickness measured from the center of the medium. However, it is found that this formula is not accurate. In this study, analytical solutions in sphere and slab geometries are derived for both optically thin and optically thick limits, assuming isotropic scattering. These solutions are verified using Monte Carlo simulations. In the optically thick limit, the number of scatterings is found to be 0.5 τ² and 1.5 τ² in a sphere and slab, respectively. In the optically thin limit, the number of scatterings is ≈ τ in a sphere and ≈ τ (1 - γ - ln τ + τ) in a slab, where γ ≃ 0.57722 is the Euler-Mascheroni constant. Additionally, we present approximate formulas that reasonably reproduce the simulation results well in intermediate optical depths. These results are applicable to scattering processes that exhibit forward and backward symmetry, including both isotropic and Thomson scattering.

• 한국컴퓨터정보학회논문지
• /
• 제20권3호
• /
• pp.89-96
• /
• 2015

IT 기술이 생체역학 분야와 폭넓게 접목되고 있으며 AHRS 센서가 스포츠 모션분석 분야에 소형화 및 가격 경쟁력 측면에서 조명을 받고 있다. 본 논문에서는 피겨스케이트화에 소형의 AHRS 센서를 부착하고, 스핀(spin), 점프, 전/후진, 인/아웃 에지, 토(toe) 등의 기본 동작을 AHRS를 통해 측정한다. AHRS 센서의 측정 오차를 줄이기 위해 Madgwick의 상보필터를 적용하였으며, 짐벌락 현상(Gimbal Lock)을 줄이기 위해 쿼터니언(Quaternion)을 이용하였다. 취득한 9축 궤적 정보에 대해 PCA, ICA, LDA, SVM의 패턴인식 알고리즘을 적용하여 인식정확도 및 실행시간을 구하고, 여러 패턴인식 알고리즘 중에서 어떤 알고리즘이 인식정확도 및 실행시간 측면에서 적용이 가능한지 제시한다. 실험결과, PCA, ICA는 인식정확도가 낮아 사용하기에 부적합하며 LDA, SVM은 인식정확도가 우수하여 피겨스케이팅 기본 동작 인식에 사용이 적합함을 보인다.

• 수산해양기술연구
• /
• 제43권1호
• /
• pp.49-61
• /
• 2007

Otter boards in the trawl are the one of essential equipments for the net mouth to be spread to the horizontal direction. Its performance should be considered in the light of the spreading force to the drag and the stability of towing in the water. Up to the present, studies of the otter boards have focused mainly on the drag and lift force, but not on the stability of otter boards movement in 3 dimensional space. In this study, the otter board is regarded as a rigid body, which has six degrees of freedom motion in three dimensional coordinate system. The forces acting on the otter boards are the underwater weight, the resistance of drag and spread forces and the tension on the warps and otter pendants. The equations of forces were derived and substituted into the governing equations of 6 degrees of freedom motion, then the second order of differential equations to the otter boards were established. For the stable numerical integration of this system, Backward Euler one of implicit methods was used. From the results of the numerical calculation, graphic simulation was carried out. The simulations were conducted for 3 types of otter boards having same area with different aspect ratio(${\lambda}=0.5,\;1.0,\;1.5$). The tested gear was mid-water trawl and the towing speed was 4k't. The length of warp was 350m and all conditions were same to each otter board. The results of this study are like this; First, the otter boards of ${\lambda}=1.0$ showed the longest spread distance, and the ${\lambda}=0.5$ showed the shorted spread distance. Second, the otter boards of ${\lambda}=1.0$ and 1.5 showed the upright at the towing speed of 4k't, but the one of ${\lambda}=0.5$ heeled outside. Third, the yawing angles of three otter boards were similar after 100 seconds with the small oscillation. Fourth, it was revealed that the net height and width are affected by the characteristics of otter boards such as the lift coefficient.

• 한국지반공학회논문집
• /
• 제22권9호
• /
• pp.45-59
• /
• 2006

탄소성 구성방정식은 주로 미분방적식(rate equation)으로 이루어져 있기 때문에 유한요소법 등을 이용한 지반구조물 해석시 미분방정식들에 대한 수치적분을 수행할 수 있는 방법이 필요하다. 구조물의 거동을 해석할시 미분방정식들을 위한 적분방법은 해석결과의 정확성과 유한요소법 모델링의 안전성에 큰 영향을 미치고 있다. 본 논문에서는 최근에 개발되어 사용되고 있는 흙에 관한 구성모델인 "Two-surface soil plasticity model (Manzari and Dafalias 1997)"을 Implicit return-mapping 수치적분방법을 이용하여 실행하는 과정을 제시한다. 본 연구에서 사용된 수치적분방법은 Closest-Point-Projection Method(CPPM) 방법으로 탄성 예측자-소성 교정자(elastic predictor-plastic corrector) 개념을 Implicit Backward Euler방법으로 체계화 시킨 알고리듬이다. 본 연구에서 수행한 "Two-surface soil plasticity model"은 조립토의 비선형거동을 해석하며, Bounding surface 개념 및 비선형 등방경화와 이동경화법칙을 사용하는 모델이다. 본 연구는 CPPM 방법이 정확하고 안정되며 유용한 수치적분을 수행할 수 있는 알고리듬이라는 것을 제시한다. 또한, CPPM 알고리듬은 구성방정식의 해를 반복적으로 해석하는 동안 "Consistent tangent operator $d{\sigma}/d{\varepsilon}$"를 제공하므로, 비선형 유한요소 해석이 2차(quadratic convergence rate)의 수렴 조건을 만족하는데 기여한다는 것을 보여준다.