• Tribology and Lubricants
• /
• 제12권3호
• /
• pp.115-122
• /
• 1996

In this paper, a finite difference method and the Newton-Raphson method are used to evaluate the Hamrock and Dowson's EHL film thickness formulas in elliptical contact problems. The minimum and central film thicknesses are compared with the Hamrock and Dowson's numerical results for various dimensionless parameters and with their film thickness formulas. The results of present analysis are more accurate and physically reasonable. The minimum film thickness formula is similar with the Hamrock and Dowson's results, however, the central film thickness formula shows large differences. Therefore, the Hamrock and Dowson's central film thickness formula should be replaced by following equation. $H_{c} = 4.88U^{0.68}G^{0.44}W^{0.096}(1-0.58e^{-0.60k})$ More accurate film thickness formula for general elliptical contact problems can be expected using present numerical methods and further research should be required.

• Korea-Australia Rheology Journal
• /
• 제21권4호
• /
• pp.289-297
• /
• 2009

We numerically analyzed flow characteristics of the polymer melt in the screw equipment using a proper modeling and investigated design parameters which have influence on the mixing performance as the capability of the screw equipment. We considered the non-Newtonian and non-isothermal flow in a single rotor equipment to investigate the mixing performance with respect to screw dimensions as shape parameter of the single rotor equipment and screw speed as process parameter. We used Bird-Carreau-Yasuda model as a viscous model of the polymer melt and the particle tracking method to investigate the mixing performance in the screw equipment and considered four mixing performance indexes: residence time distribution, deformation rate, total strain and particle standard deviation as a new mixing performance index. We compared these indexes to determine design parameters and object function. On basis of the analysis results, we carried out the optimal design by using the response surface method and design of experiments. In conclusion, the differences of results between the optimal value and numerical analysis are about 5.0%.

• International Journal of Fluid Machinery and Systems
• /
• 제2권3호
• /
• pp.254-259
• /
• 2009

In this paper it was attempted to treat the hydrodynamic journal bearing as a time-based nonlinear reaction source in each step of rotor rotation in order to observe the bearing effect more realistically and accurately in stead of the conventional method of simple linearized stiffness and damping. Lubrication analysis based on finite element method is employed to calculate the hydrodynamic reaction of bearing and Newmark's method was used to calculate the rotor dynamics in the time domain. Simulation for an industrial electrical motor showed remarkable results with differences compared to those by the conventional method in the dynamic behavior of the rotor.

• Journal of applied mathematics & informatics
• /
• 제6권1호
• /
• pp.175-184
• /
• 1999

Along with the computation and analysis for nonlinear system being more and more involved in the fields such as automation control electronic technique and electrical power system the nonlin-ear theory has become quite a attractive field for academic research. In this paper we derives the solutions for state equation of nonlinear system by using the inverse operator expression of the so-lutions is obtained. An actual computation example is given giving a comparison between IOM and Runge-kutta method. It has been proved by our investigation that IOM has some distinct advantages over usual approximation methods in that it is computationally con-venient rapidly convergent provides accurate solutions not requiring perturbation linearization or the massive computation inherent in discrietization methods such as finite differences. So the IOM pro-vides an effective method for the solution of nonlinear system is of potential application valuable in nonlinear computation.

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

A finite-difference method based on conservative supra characteristic method type upwind flux difference splitting has been developed to study the nonequilibrium chemically reacting inviscid flow. For nonequilibrium air, NS-1 species equations were strongly coupled with flowfield equations through convection and species production terms. Inviscid nonequilibrium chemically reacting air mixture flows over Blunt-body were solved to demonstrate the capability of the current method. At low altitude flight conditions the nonequilibrium air models predicted almost the same temperature, density and pressure behind the shock as equilibrium flow: however, at high altitudes they showed substantial differences due to nonequilibrium chemistry effect. The new nonequilibrium chemically reacting upwind flux difference splitting mettled can be extended to viscous flow and multi-dimensional flow conditions.

• Nuclear Engineering and Technology
• /
• 제55권2호
• /
• pp.734-748
• /
• 2023

A comprehensive investigation of the Ronen method is performed in homogeneous and heterogeneous slab problems from the Sood benchmark, considering isotropic and linearly-anisotropic problems. Three finite differences implementations are exercised and compared. The results are compared to reference solutions using one and two energy groups. The validation is performed for the criticality eigenvalue and the fundamental neutron flux distribution. The results demonstrate the significantly improved accuracy achievable by the Ronen method using a broad set of problems. For standard convergence tolerances, the maximal deviation in criticality eigenvalue is less than ten pcm, and the maximal deviation in the spatial distribution of the flux is less than 2%, always located near sharp interfaces or vacuum boundaries.

• Computers and Concrete
• /
• 제18권3호
• /
• pp.319-336
• /
• 2016

In this study, both two- and three-dimensional (2D and 3D) finite-volume-based models were developed to analyze the heat transfer mechanisms through the porous structures of cellular concretes under steady-state heat transfer conditions and to investigate the differences between the 2D and 3D modeling results. The 2D and 3D reconstructed pore networks were generated from the microstructural information measured by 3D images captured by X-ray computerized tomography (X-CT). The computed effective thermal conductivities based on the 2D and 3D calculations performed on the reconstructed porous structures were found to be nearly identical to those evaluated from the 2D cross-sectional images and the 3D X-CT images, respectively. In addition, the 3D computed effective thermal conductivity was found to agree better with the measured values, in comparison with the 2D reconstruction and real cross-sectional images. Finally, the thermal conductivities computed for different reconstructed porous 3D structures of cellular concretes were compared with those obtained from 2D computations performed on 2D reconstructed structures. This comparison revealed the differences between 2D and 3D image-based modeling. A correlation was thus derived between the results of the 3D and 2D models.

• 마이크로전자및패키징학회지
• /
• 제10권2호
• /
• pp.19-26
• /
• 2003

광통신용 광학부품의 신뢰성 특성은 솔더 조인트의 열 사이클에 따른 소성(Plastic)과 크립(Creep) 변형에 가장 큰 영향을 받는다. 열 사이클에 따른 소성과 크립 변형 증가로 인해 정렬 틀어짐이 발생하며 이는 광손실 변화의 주요인이 된다. 또한, 소성과 크립 변형량이 증가 또는 계속 누적이 될 경우 솔더의 피로수명 한계로 인해 제품 불량 발생의 원인이 된다. 이러한 열적 사이클에 따른 광부품의 신뢰성을 확보하기 위해 본 논문에서는 유한요소해석법(FEM)을 적용하였다. 소성과 크립 변형의 변화량을 유한요소해석으로 계산하고 이를 크립 피로 파괴(Creep-Fatigue) 수명 예측 모델에 적용하여 그 수명을 예측하였다. 솔더와 모재와의 계면 또는 솔더 내부에서 생성되는 온도에 따른 소성과 크립 변형을 파악하기 위해 텔코디아(Telcordia)의 광부품 신뢰성 온도 사이클(-40 to 75)을 적용하였다. 승온과 냉각 속도의 변화에 따른 영향을 검토하기 위해 1/min, 10/min 및 50/min으로 변화를 주고 유지 시간을 1시간으로 고정할 경우의 결과를 비교 분석하였다.

• 대한치과보철학회지
• /
• 제45권2호
• /
• pp.274-282
• /
• 2007

Statement of problem: Main consideration was given to the stresses at the site of implant entry into the cortical bone at the alveolar crest. As a suspectible factor affecting the occurrence of stress concentrations, the contact angle between the implant and the alveolar crest bone was addressed. Purpose: The purpose of this study is to evaluate angles between the alveolar crest bone and the implant effect on the implant crestal area induced stresses using a finite element method. Material and methods: Cylindrically shaped, standard size ITI implants entering into alveolar crest with four different contact angles of 0, 15, 30, and 45 deg. with the long axis of the implant were axisymmetrically modelled. Alterations of stresses around the implants were computed and compared at the cervical cortical bone. Results and conclusion: The results demonstrated that regardless of the difference of the implant/alveolar crest bone contact angles, stress concentration occurred at the cervical bone and the angle differences led to insignificant variations in stress level.

• Steel and Composite Structures
• /
• 제15권5호
• /
• pp.481-505
• /
• 2013

This paper focuses on thermal post-buckling analysis of functionally graded beams with temperature dependent physical properties by using the total Lagrangian Timoshenko beam element approximation. Material properties of the beam change in the thickness direction according to a power-law function. The beam is clamped at both ends. In the case of beams with immovable ends, temperature rise causes compressible forces and therefore buckling and post-buckling phenomena occurs. It is known that post-buckling problems are geometrically nonlinear problems. Also, the material properties (Young's modulus, coefficient of thermal expansion, yield stress) are temperature dependent: That is the coefficients of the governing equations are not constant in this study. This situation suggests the physical nonlinearity of the problem. Hence, the considered problem is both geometrically and physically nonlinear. The considered highly non-linear problem is solved considering full geometric non-linearity by using incremental displacement-based finite element method in conjunction with Newton-Raphson iteration method. In this study, the differences between temperature dependent and independent physical properties are investigated for functionally graded beams in detail in post-buckling case. With the effects of material gradient property and thermal load, the relationships between deflections, critical buckling temperature and maximum stresses of the beams are illustrated in detail in post-buckling case.