• ETRI Journal
• /
• v.28 no.1
• /
• pp.107-110
• /
• 2006

The dimensions of the human body vary by age, sex, and race. The internal structure and outer dimensions of a body exposed to an electromagnetic field is important for accurate dosimetry. The average physical size of Korean adult males between the ages 18 to 24 was investigated, and a male volunteer was selected whose physical condition is within the physical standards, ${\pm}5%$. Magnetic resonance images and partially computerized tomography images of the volunteer were acquired. The intervals between the transverse images were 1 mm for the head and 3 mm for the rest of the body. About 30 different tissues were manually classified by an anatomist on the raw images, and the segmented images were implemented in the form of a text file appropriate for numerical formulation.

• 한국전산유체공학회:학술대회논문집
• /
• 2003.10a
• /
• pp.191-192
• /
• 2003

A numerical method for a wave diffraction problem in three-dimensional channels is developed. The physical models are various shapes of channel connected to the open sea. When a ship or an offshore structure is moored in various configurations of channel connected to an open sea, the prediction of the hydrodynamic force exerting on the moored ship could be important for the prediction of its motion. It is assumed that the fluid is inviscid and incompressible and its motion is irrotational. From the continuity equation, the Laplace equation can be obtained as the governing equation. The surface tension at free surface is neglected, and wave amplitude is assumed to be small compared to the wave length. Then the free surface condition can be linearized. The numerical method used here is the localized finite element method based on a variational formulation

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.10a
• /
• pp.75-78
• /
• 2003

In order to describe the Bauschinger and transient behavior of orthotropic fiber-reinforced composites, a combined isotropic-kinematic hardening law based on the non-linear kinematic hardening rule was considered here, in particular, based on the Chaboche type law. In this modified constitutive law, the anisotropic evolution of the back-stress was properly accounted for. Also, to represent the orthotropy of composite materials, Hill's 1948 quadratic yield function and the orthotropic elasticity constitutive equations were utilized. Furthermore, the numerical formulation to update the stresses was also developed based on the incremental deformation theory for the boundary value problems. Numerical examples confirmed that the new law based on the anisotropic evolution of the back-stress complies well with the constitutive behavior of highly anisotropic materials such as fiber-reinforced composites.

• Journal of Advanced Marine Engineering and Technology
• /
• v.27 no.2
• /
• pp.289-298
• /
• 2003

Meshfree approximations exhibit significant Potential to solve partial differential equations. Meshfree methods have been successfully applied to various problems which the traditional finite element methods have difficulties to handle including the quasi-static and dynamic fracture, large deformation problems, contact problems, and strain localization problems. Reproducing Kernel Particle Method (RKPM) is used in this research fur to its built-in feature of multi-resolution. the sound mathematical foundation and good numerical performance. A formulation of RKPM is reviewed and numerical examples are given to verify the accuracy of the proposed meshfree method for largely deformed elasto-plastic material.

• Proceedings of the KSME Conference
• /
• 2001.11a
• /
• pp.341-348
• /
• 2001

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, it will be established that this new method is very accurate and effective for solving plane elastic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.23 no.10
• /
• pp.1296-1309
• /
• 1999

A three dimensional compressible Navier-Stokes code is developed to analyze flowfields and performance of a vaned diffuser in a centrifugal compressor. It employs scalar implicit approximate factorization, finite volume formulation, second order upwind differencing and a two-equation $q-{\omega}$ turbulence model based on the integration to the wall. Pressure recovery and loss coefficients of a vaned diffuser are evaluated using a developed computer code. The simulated three dimensional flows show how through flow structure affects pressure recovery performance and loss coefficients of a vane for design and off-design inlet flow angles. Development of complex three dimensional flow over the inlet region and leading edge are very influential to the overall flow and performance.

• Korea-Australia Rheology Journal
• /
• v.18 no.3
• /
• pp.161-167
• /
• 2006

This study presents results on the numerical simulation of Newtonian and non-Newtonian flow in a channel obstructed by an asymmetric array of obstacles for clarifying the descriptive ability of current non-Newtonian constitutive equations. Jones and Walters (1989) have performed the corresponding experiment that clearly demonstrates the characteristic difference among the flow patterns of the various liquids. In order to appropriately account for flow properties, the Navier-Stokes, the Carreau viscous and the Leonov equations are employed for Newtonian, shear thinning and extension hardening liquids, respectively. Making use of the tensor-logarithmic formulation of the Leonov model in the computational scheme, we have obtained stable solutions up to relatively high Deborah numbers. The peculiar characteristics of the non-Newtonian liquids such as shear thinning and extension hardening seem to be properly illustrated by the flow modeling. In our opinion, the results show the possibility of current constitutive modeling to appropriately describe non-Newtonian flow phenomena at least qualitatively, even though the model parameters specified for the current computation do not precisely represent material characteristics.

• The Korean Journal of Rheology
• /
• v.2 no.1
• /
• pp.33-45
• /
• 1990

This study analyzes the planar 4:1 contraction flow of viscoelastic fluids with retardation time using finite volume method. To consider separately the elasticity effect of the viscoelastic fluid without shear thinn-ing effect, Oldroyd B liquid model is adopted for the numerical simulation. Instead of the stream function-vorticity formulation, SIMPLER algorithm with staggered grid system which incorporates primitive variable has been introduced in discretizing the momentum equations. An upwind corrected scheme has been used in discetizing the constitutive equations for the non-Newtonian part of the stress. The size of the corner vortex is shown to be slightly influenced by the Weissenberg number. However as the Weissenberg number is increased the chang-ing of the vortex shape agrees qualitatively well with some experimental studies.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 1993.04a
• /
• pp.52-56
• /
• 1993

An efficient method is presented for determining transient responses of locally nonlinear structures using substructure eigenproperties and Lagrange multiplier technique. The method is based upon the mode synthesis formulation procedure, but does not construct the equations of motion of the combined whole structure compared with the conventional methods. For modal bases of each linear substructure, either fixed or free interface modes can be employed. The transient analysis is based upon the recurrence discrete-time state equations and offers the simplicity of the Euler integration method without requiring small time increment and iterative solution procedure. Numerical examples reveal that the method is very accurate and efficient in calculating transient responses compared with the direct numerical integration method.

• Structural Engineering and Mechanics
• /
• v.5 no.3
• /
• pp.221-237
• /
• 1997

This paper presents a new uniaxial material model for rate-sensitive analysis addressing both the transient and steady-state responses. The new model adopts visco-plastic theory for the rate-sensitive response, and employs a three-parameter representation of the overstress as a function of the strain-rate. The third parameter is introduced in the new model to control its transient response characteristics, and to provide flexibility in fitting test data on the variation of overstress with strain-rate. Since the governing visco-plastic differential equation cannot be integrated analytically due to its inherent nonlinearity, a new single-step numerical integration procedure is proposed, which leads to high levels of accuracy almost independent of the size of the integration time-step. The new model is implemented within the nonlinear analysis program ADAPTIC, which is used to provide several verification examples and comparison with other experimental and numerical results. The companion paper extends the three-parameter model to trilinear static stress-strain relationships for steel and concrete, and presents application examples of the proposed models.