• Journal of Mechanical Science and Technology
• /
• v.20 no.12
• /
• pp.2230-2241
• /
• 2006

This paper presents a sole application of boundary element method to the conjugate heat transfer problem of thermally developing laminar flow in a thick walled pipe when the fluid velocities are fully developed. Due to the coupled mechanism of heat conduction in the solid region and heat convection in the fluid region, two separate solutions in the solid and fluid regions are sought to match the solid-fluid interface continuity condition. In this method, the dual reciprocity boundary element method (DRBEM) with the axial direction marching scheme is used to solve the heat convection problem and the conventional boundary element method (BEM) of axisymmetric model is applied to solve the heat conduction problem. An iterative and numerically stable BEM solution algorithm is presented, which uses the coupled interface conditions explicitly instead of uncoupled conditions. Both the local convective heat transfer coefficient at solid-fluid interface and the local mean fluid temperature are initially guessed and updated as the unknown interface thermal conditions in the iterative solution procedure. Two examples imposing uniform temperature and heat flux boundary conditions are tested in thermally developing region and compared with analytic solutions where available. The benchmark test results are shown to be in good agreement with the analytic solutions for both examples with different boundary conditions.

• Coupled systems mechanics
• /
• v.9 no.1
• /
• pp.5-28
• /
• 2020

Despite a widely-held belief that the finite element method is the method for the solution of solid mechanics problems, which has for 30 years dissuaded solid mechanics scientists from paying any attention to the finite volume method, it is argued that finite volume methods can be a viable alternative. It is shown that it is simple to understand and implement, strongly conservative, memory efficient, and directly applicable to nonlinear problems. A number of examples are presented and, when available, comparison with finite element methods is made, showing that finite volume methods can be not only equal to, but outperform finite element methods for many applications.

• Proceedings of the KSME Conference
• /
• 2008.11a
• /
• pp.442-447
• /
• 2008

This paper presents a new approach to simulate fluid-solid interaction problems involving non-matching interfaces. The coupling between fluid and solid domains with dissimilar finite element meshes consisting of 4-node quadrilateral elements is achieved by using the interface element method (IEM). Conditions of compatibility between fluid and solid meshes are satisfied exactly by introducing the interface elements defined on interfacing regions. Importantly, a consistent transfer of loads through matching interface element meshes guarantees the present method to be an efficient approach of the solution strategy to fluid-solid interaction problems. An arbitrary Lagrangian-Eulerian (ALE) description is adopted for the fluid domain, while for the solid domain an updated Lagrangian formulation is considered to accommodate finite deformations of an elastic structure. The stabilized equal order velocity-pressure elements for incompressible flows are used in the motion of fluids. Fully coupled equations are solved simultaneously in a single computational domain. Numerical results are presented for fluid-solid interaction problems involving nonmatching interfaces to demonstrate the effectiveness of the methodology.

• Journal of Korean Association for Spatial Structures
• /
• v.8 no.4
• /
• pp.81-90
• /
• 2008

This paper presents a finite element formulation for the three-dimensional hierarchical solid element using Integrals of Legendre polynomials. The proposed hexahedral solid element is composed of four different modes including vertex, edge, face, and internal mode, respectively. The eigenvalue and patch test have been carried out to confirm the zero-energy mode and constant strain condition. In addition to these, a posteriori error estimation has been studied for the p-adaptive finite element analysis that is based on a smoothing technique to compute a post-processed solution from the finite element solution. The uniform p-refinement and non-uniform p-refinement are compared in terms of convergence rate as the number of degree of freedom is increased. The simple cantilever beam is tested to show the performance of the proposed solid element.

• Structural Engineering and Mechanics
• /
• v.76 no.3
• /
• pp.379-393
• /
• 2020

In this paper, we propose an automatic procedure to improve the accuracy of finite element solutions using enriched 2D solid finite elements (4-node quadrilateral and 3-node triangular elements). The enriched elements can improve solution accuracy without mesh refinement by adding cover functions to the displacement interpolation of the standard elements. The enrichment scheme is more effective when used adaptively for areas with insufficient accuracy rather than the entire model. For given meshes, an error for each node is estimated, and then proper degrees of cover functions are applied to the selected nodes. A new error estimation method and cover function selection scheme are devised for the proposed adaptive enrichment scheme. Herein, we demonstrate the proposed enrichment scheme through several 2D problems.

• Structural Engineering and Mechanics
• /
• v.28 no.6
• /
• pp.695-713
• /
• 2008

The aim of this contribution is to present a new link/beam finite element suitable for electrothermo-structural analysis of uni-axially graded materials. Continuous polynomial variation of geometry and material properties will be considered. The element matrix and relations for solution of Joule's heat (and its distribution to the element nodes) have been established in the sense of a sequence method of a coupled problem solution. The expression for the solution of nodal forces caused by a continuously distributed temperature field has also been derived. The theoretical part of this contribution is completed by numerical validation, which proves the high accuracy and effectiveness of the proposed element. The results of the performed experiments are compared with those obtained using the more expensive multiphysical link element and solid element of the FEM program Ansys. The proposed finite element could be used not only in the multiphysical analysis of the current paths and actuators but also in analysis of other 1D construction parts made of composite or uni-axially graded materials.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.12 no.4
• /
• pp.642-651
• /
• 1988

The present paper develops finite element procedures to calculate displacements, strains and stresses in initially stressed elastic solids subjected to static or time-dependent loading conditions. As a point of departure, we employ Hamilton's principle to obtain nonlinear equations of motion characterizing the displacement in a solid. The equations of motion reduce to linear equations of motion if incremental stresses are assumed to be infinitesimal. In the case of linear problem, finite element solutions are obtained by Newmark's direct integration method and by modal analysis. An analytic solution is referred to compare with the linear finite element solution. In the case of nonlinear problem, finite element solutions are obtained by Newton-Raphson iteration method and compared with the linear solution. Finally, the effect of the order of Gauss-Legendre numerical integration on the nonlinear finite element solution, has been investigated.

• Nuclear Engineering and Technology
• /
• v.54 no.8
• /
• pp.2801-2808
• /
• 2022

The synergistic effect of ZrC nanoparticle pining and Re solution in W matrix on the thermal stability of tungsten was studied by investigating the evolution of the microstructure, hardness and tensile properties after annealing in a temperature range of 1000-1700 ℃. The results of metallography, electron backscatter diffraction pattern and Vickers micro-hardness indicate that the rolled W-1wt%Re-0.5 wt% ZrC alloy has a higher recrystallization temperature (1600 ℃-1700 ℃) than that of the rolled pure W (1200 ℃), W-0.5 wt%ZrC (1300 ℃), W-0.5 wt%HfC (1400-1500 ℃) and W-K-3wt%Re alloy fabricated by the same technology. The molecular dynamics simulation results indicated that solution Re atoms in W matrix can slow down the self-diffusion of W atoms and form dragging effect to delay the growth of W grain, moreover, the diffusion coefficient decrease with increasing Re content. In addition, the ZrC nanoparticles can pin the grain boundaries and dislocations effectively, preventing the recrystallization. Therefore, synergistic effect of solid solution Re element and dispersed ZrC nanoparticles significantly increase recrystallization temperature.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.22 no.3
• /
• pp.284-291
• /
• 2002

An attempt was made to evaluate the degree of aging degradation in thermally aged 2.25Cr-1Mo steel by electrical resistivity measurement. Artificial aging was performed to simulate the microstructural degradation in 2.25Cr-1Mo steel arising from long time exposure at $540^{\circ}C$. Microstructural parameter (amount of solid solution element), mechanical property (ductile-brittle transition temperature) and electrical resistivity were measured to investigate the mutual relationship among these parameters. Depletion of solid solution element(Mo and Cr) in matrix was detected after aging. The ductile-brittle transition temperature(DBTT) increased rapidly in the initial stage of aging and then saturated afterward. On the other hand, the electrical resistivity decreased rapidly in the beginning and then saturated in the later stage of aging.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.16 no.1
• /
• pp.29-36
• /
• 2008

In this paper, the orbital forming analysis of an automotive hub bearing was studied to predict forming performances using the explicit finite element method. To find an efficient solution technique for the orbital forming, axisymmetric finite element models and 3D solid element models were solved and numerically compared. The time scaling and mass scaling techniques were introduced to reduce the excessive computational time caused by small element size in case of the explicit finite element method. It was found from the numerical simulations on the orbital forming that the axisymmetric element models showed the similar results to the 3D solid element models in forming loads whereas the deformations at the inner race of bearing were quite different. Finally the strains at the inner race of bearing and the forming forces to the peen were measured for the same product of the numerical model by test, and were compared with the 3D solid element results. It was founded that the test results were in good agreements with the numerical ones.