• Wind and Structures
• /
• v.12 no.1
• /
• pp.1-19
• /
• 2009

In this paper, the airflow around an ideal thin plate (hereafter referred to as ITP) with various ratios of central slot is simulated by using the finite-difference-method (FDM)-based Arbitrary-Lagrangian-Eulerian descriptions for the rigid oscillating body. The numerical procedure employs the second-order projection scheme to decouple the governing equations, and the multigrid algorithm with three levels to improve the computational efficiency in evaluating of the pressure equation. The present CFD method is validated through comparing the computed flutter derivatives of the ITP without slot to Theodorsen analytical solutions. Then, the unsteady aerodynamics of the ITP with and without central slot is investigated. It is found that even a smaller ratio of central slot of the ITP has notable effects on pressure distributions of the downstream section, and the pressure distributions on the downstream section will further be significantly affected by the slot ratio and the reduced wind speeds. Continuous increase of $A_2^*$ with the increase of central slot may be the key feature of the slotted ITP. Finally, flutter analyses based on the flutter derivatives of the slotted ITP are performed, and moreover, flutter instabilities of a scaled sectional model of a twin-deck bridge with various ratios of deck slot are investigated. The results confirm that the central slot is effective to improve bridge flutter stabilities, and that the flutter critical wind speeds increase with the increase of slot ratio.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.166-171
• /
• 2000

A two-dimensional/axisymmetric CSCM upwind flux difference splitting Wavier-Stokes method has been developed to study the finite rate chemically react-ing invisicd and viscous hypersonic flows over blunt-body. A upwind method was chosen due to its robustness in capturing the strong bow shock waves. For the nonequilibrium chemically reacting air, NS-I species conservation equations were strongly coupled with flowfield equations through convection and species production terms. The nonequilibrium wall pressure and heat transfer rate distributions along the vehicle were compared with those from equilibrium and perfect gas calculations. The nonequilibrium species distribution shows the reduced concentrations of O and N species when compared with equilibrium species distribution. The solutions resolved strong bow shock waves md heat transfer rate very accurately when compared with central difference schemes.

• The Transactions of the Korean Institute of Electrical Engineers C
• /
• v.51 no.11
• /
• pp.559-566
• /
• 2002

In this paper, we present an accurate and stable method for the solution of the transient electromagnetic response from the conducting wire structures using the time domain integral equation. By using an implicit scheme with the central finite difference approximation for the time domain electric field integral equation, we obtain the transient response from a wire scatterer illuminated by a plane wave and a conducting wire antenna with an impressed voltage source. Also, we consider a wire above a 3-dimensional conducting object. Numerical results are presented, which show the validity of the presented methodology, and compared with a conventional method using backward finite difference approximation.

• Journal of the Korean Mathematical Society
• /
• v.44 no.3
• /
• pp.627-645
• /
• 2007

We investigate the eigenvalues of the semi-circulant preconditioned matrix for the finite difference scheme corresponding to the second-order elliptic operator with the variable coefficients given by $L_vu\;:=-{\Delta}u+a(x,\;y)u_x+b(x,\;y)u_y+d(x,\;y)u$, where a and b are continuously differentiable functions and d is a positive bounded function. The semi-circulant preconditioning operator $L_cu$ is constructed by using the leading term of $L_vu$ plus the constant reaction term such that $L_cu\;:=-{\Delta}u+d_cu$. Using the field of values arguments, we show that the eigenvalues of the preconditioned matrix are clustered at some number. Some numerical evidences are also provided.

• Tribology and Lubricants
• /
• v.12 no.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.

• Tribology and Lubricants
• /
• v.28 no.4
• /
• pp.153-159
• /
• 2012

Tapered roller bearings are widely used in high axial-load and radial-load applications. In this study, a numerical analysis is performed to study a finite line contacts EHL problem between a tapered roller and raceway in tapered roller bearings. Converged solutions are obtained for moderate load and material parameters using a finite difference method with non-uniform grids and the Newton-Raphson method. The contours and sectional plots of pressure distribution and film shape are compared. The pressure distribution and film shapes near both ends of the roller are very different from those in the central part and are transversely asymmetric. The maximum pressure and absolute minimum film thickness always occur at the small end of the roller.

• Tribology and Lubricants
• /
• v.29 no.6
• /
• pp.353-359
• /
• 2013

The rollers and/or races in cylindrical and tapered roller bearings should be profiled to relieve high edge stress concentrations caused by their finite lengths and misalignment. In this study, a numerical analysis was performed to investigate the elastohydrodynamic lubrication (EHL) of a Lundberg profile-type cylindrical roller. A finite difference method with fully nonuniform grids and the Newton-Raphson method were used to present detailed EHL pressure distributions and film shapes, as well as the variations in the minimum and central film thicknesses with the profile modification coefficient. In the Lundberg profile, the maximum pressure and minimum film thickness always occurred near the edges. Proper modification of the Lundberg profile considerably increased the minimum film thickness.

• Structural Engineering and Mechanics
• /
• v.2 no.2
• /
• pp.173-184
• /
• 1994

The objective of this study was to develop a simple and efficient numerical modeling technique for dynamic crack propagation using the finite element method. The study focused on the analysis of a rapidly propagation crack in an elastic body. As already known, discrete crack tip advance with the stationary node procedure results in spurious oscillation in the calculated energy terms. To reduce the spurious oscillation, a simple and efficient moving node procedure is proposed. The procedure does require neither remeshing the discretization nor distorting the original mesh. Two different central difference schemes are also evaluated and compared for dynamic crack propagation problem.

• Journal of Dental Rehabilitation and Applied Science
• /
• v.18 no.4
• /
• pp.301-311
• /
• 2002

The purpose of this study was to compare the distributing pattern of stress on the finite element models with the different vertical bone level of implant fixture. The two kinds of finite element models were designed according to vertical bone level around fixture ($4.0mm{\times}11.5mm$). The cemented crowns for mandibular first and second molars were made. Three- dimensional finite element model was created with the components of the implant and surrounding bone. Vertical loads were applied with force of 200N distributed within 0.5mm radius circle from the center of central fossa and distance 2mm and 4 mm apart from the center of central fossa. Von-Mises stresses were recorded and compared in the supporting bone, fixtures, abutment screws, and crown. The results were as following : (1) In vertical loading at the center circle of central fossa on model 1 and 2, the difference from vertical bone in implant placement did not affect the stress pattern on all components of implant except for crown. (2) With offset distance incerasing and the bone level of implant decreasing, the concentration of stress occured in the buccal side of long crown, around the buccal crestal bone, and on the fixture- abutment interface. As a conclusion, the research showed a tendency to increase the stress on the supporting bone, fixture and screw under the offset loads when the vertical level of bone around fixture was different. Since the same vertical bone bed has more benefits than the different bone bed around fixtures, it is important to prepare a same vertical level of bone bed for the success of implants under occlusal loads.

• Restorative Dentistry and Endodontics
• /
• v.34 no.1
• /
• pp.8-19
• /
• 2009

The purpose of this study was to investigate the effects of four restorative materials under various occlusal loading conditions on the stress distribution at the CEJ of buccal. palatal surface and central groove of occlusal surface of endodontically treated maxillary second premolar, using a 3D finte element analysis. A 3D finite element model of human maxillary second premolar was endodontically treated. After endodontic treatment, access cavity was filled with Amalgam, resin, ceramic or gold of different mechanical properties. A static 500N forces were applied at the buccal (Load-1) and palatal cusp (Load-2) and a static 170N forces were applied at the mesial marginal ridge and palatal cusp simultaneously as centric occlusion (Load-3). Under 3-type Loading condition, the value of tensile stress was analyzed after 4-type restoration at the CEJ of buccal and palatal surface and central groove of occlusal surface Excessive high tensile stresses were observed along the palatal CEJ in Load-1 case and buccal CEJ in Load-2 in all of the restorations. There was no difference in magnitude of stress in relation to the type of restorations. Heavy tensile stress concentrations were observed around the loading point and along the central groove of occlusal surface in all of the restorations. There was slight difference in magnitude of stress between different types of restorations. High tensile stress concentrations around the loading points were observed and there was no difference in magnitude of stress between different types of restorations in Load-3.