• Steel and Composite Structures
• /
• v.31 no.5
• /
• pp.529-539
• /
• 2019

Current paper deals with thermoelastic static and free vibrational behaviors of axisymmetric thick cylinders reinforced with functionally graded (FG) randomly oriented graphene subjected to internal pressure and thermal gradient loads. The heat transfer and mechanical analyses of randomly oriented graphene-reinforced nanocomposite (GRNC) cylinders are facilitated by developing a weak form mesh-free method based on moving least squares (MLS) shape functions. Furthermore, in order to estimate the material properties of GRNC with temperature dependent components, a modified Halpin-Tsai model incorporated with two efficiency parameters is utilized. It is assumed that the distributions of graphene nano-sheets are uniform and FG along the radial direction of nanocomposite cylinders. By comparing with the exact result, the accuracy of the developed method is verified. Also, the convergence of the method is successfully confirmed. Then we investigated the effects of graphene distribution and volume fraction as well as thermo-mechanical boundary conditions on the temperature distribution, static response and natural frequency of the considered FG-GRNC thick cylinders. The results disclosed that graphene distribution has significant effects on the temperature and hoop stress distributions of FG-GRNC cylinders. However, the volume fraction of graphene has stronger effect on the natural frequencies of the considered thick cylinders than its distribution.

• Advances in Computational Design
• /
• v.9 no.1
• /
• pp.39-52
• /
• 2024

If the governing differential equation arising from engineering problems is treated as an analytic, continuous and derivable function, it can be expanded by one point as a series of finite numbers. For the function to be zero for each value of its domain, the coefficients of each term of the same power must be zero. This results in a recursive relationship which, after applying the natural conditions or the boundary conditions, makes it possible to obtain the values of the derivatives of the function with acceptable accuracy. The elastoplastic analysis of an inhomogeneous thick sphere of metallic materials with linear variation of the modulus of elasticity, yield stress and Poisson's ratio as a function of radius subjected to internal pressure is presented. The Beltrami-Michell equation is established by combining equilibrium, compatibility and constitutive equations. Assuming axisymmetric conditions, the spherical coordinate parameters can be used as principal stress axes. Since there is no analytical solution, the natural boundary conditions are applied and the governing equations are solved using a proposed new method. The maximum effective stress of the von Mises yield criterion occurs at the inner surface; therefore, the negative sign of the linear yield stress gradation parameter should be considered to calculate the optimal yield pressure. The numerical examples are performed and the plots of the numerical results are presented. The validation of the numerical results is observed by modeling the elastoplastic heterogeneous thick sphere as a pressurized multilayer composite reservoir in Abaqus software. The subroutine USDFLD was additionally written to model the continuous gradation of the material.

• Journal of the Society of Naval Architects of Korea
• /
• v.50 no.6
• /
• pp.436-449
• /
• 2013

In this paper, a numerical analysis is carried out to study the characteristics of supercavitating flows and the drag of relatively simple two-dimensional and axisymmetric bodies which can be used for supercavity generation device, cavitator, of a high-speed underwater vehicle. In order to investigate the suitability of numerical models, cavity flows around the hemispherical head form and two-dimensional wedge are calculated with combinations of three turbulence models(standard $k-{\epsilon}$, realizable $k-{\epsilon}$, Reynolds stress) and two cavitation models(Schnerr-Sauer, Zwart-Gerber-Belamri). From the results, it is confirmed that the calculated cavity flow is more affected by the turbulence model than the cavitation model. For the calculation of steady state cavity flows, the convergence in case of the realizable $k-{\epsilon}$ model is better than the other turbulence models. The numerical result of the Schnerr-Sauer cavitation model is changed less by turbulence model and more robust than the Zwart-Gerber-Belamri model. Thus the realizable $k-{\epsilon}$ turbulence model and the Schnerr-Sauer cavitation model are applied to calculate supercavitating flows around disks, two dimensional $10^{\circ}$ and $30^{\circ}$ wedges. In case of the disk, the cavitation number dependences of the cavity size and the drag coefficient predicted are similar to either experimental data or Reichardt's semi-empirical equations, but the drag coefficient is overestimated about 3% higher than the Reichardt's equation. In case of the wedges, the cavitation number dependences of the cavity size are similar to experimental data and Newman's linear theory, and the agreement of the cavity length predicted and Newman's linear theory becomes better as decreasing cavitation number. However, the drag coefficients of wedges agree more with experimental data than those of Newman's analytic solution. The cavitation number dependences of the drag coefficients of both the disk and the wedge appear linear and simple formula for estimating the drag of supercavitating disks and wedges are suggested. Consequently, the CFD scheme of this study can be applied for numerical analysis of supercavitating flows of the cavitator and the cavitator design.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.9 no.2
• /
• pp.149-156
• /
• 1999

An experimental simulation on the flow in Czochralski melt using a cold model was carried out to obtain the velocities of fluid flow which affects the oxygen concentration of Czochralski crystal growing system. Low melting point Woods metal with similar Pr number to the silicon melt was adopted as a working fluid. Local flow velocities at numerous positions in the melt were simulataneously measured in three dimension using incorporated magnet probe. The measured velocity field showed a non-axisymmetric pattern dominated by natural convection. The analysis on the correlation between data set of temperatures simultaneously measured at two melt positions showed that the values of correlation coefficients were smaller than those of previous study on the small size of silicon melt and these phenomena are believed to occur because turbulent behavior becomes stronger in large size of the melt.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.2
• /
• pp.12-20
• /
• 1998

Exter ballistics of a typical high-speed projectile is studied through a flow-visualization experiment and an unstructured grid Navier-Srokes computation. Experiment produced a schlieren photograph that adequately shows the characteristic features of this complex flow, namely two kinds of oblique cone shocks and turbulent wake developing into the downstream. A hybrid scheme of finite volume-element method is used to simulate the compressible Reynolds-Averaged Navier-Stok- es solution on unstructured grids. Osher's approximate Riemann solver is used to discretize the cinvection term. Higher-order spatial accuracy is obtained by MUSCL extension and van Albada ty- pe flux limiter is used to stabilize the numerical oscillation near the solution discontinuity. Accurate Gakerkin method is used to discretize the viscous term. Explict fourth-order Runge-Kutta method is used for the time-stepping, which simplifies the application of MUSCL extension. A two-layer k-$\varepsilon$ turbulence model is used to simulate the turbulent wakes accurately. Axisymmetric folw and two-dimensional flow with an angle of attack have been computed. Grid-dependency is also checked by carrying out the computation with doubled meshes. 2-D calculation shows that effect of angle of attack on the flow field is negligible. Axi-symmetric results of the computation agrees well with the flow visualization. Primary oblique shock is represented within 2-3 meshes in numerical results, and the varicose mode of the vortex shedding is clearly captured in the turbulent wake region.

• Journal of the Korean Society of Systems Engineering
• /
• v.13 no.1
• /
• pp.33-39
• /
• 2017

This paper describes the system engineering approach for the heat transfer analysis of plus7 fuel rod for APR1400 using, a commercial software, ANSYS. The fuel rod is composed of fuel pellets, fill gas, end caps, plenum spring and cladding. The heat is transferred from the pellet outward by conduction through the pellet, fill gas and cladding and further by convection from the cladding surface to the coolant in the flow channel. The goal of this paper is to demonstrate the temperature and heat flux change from the fuel centerline to the cladding surface when having maximum fuel centerline temperature at 100% power. This phenomenon is modelled using the ANSYS FEM code and analyzed for steady state temperature distribution across the fuel pellet and clad and the results were compared to the standard values given in APR1400 SSAR. Specifically the applicability of commercial software in the evaluation of nuclear fuel temperature distribution has been accounted. It is note that special codes have been used for fuel rod mechanical analysis which calculates interrelated effects of temperature, pressure, cladding elastic and plastic behavior, fission gas release, and fuel densification and swelling under the time-varying irradiation conditions. To satisfactorily meet this objective we apply system engineering methodologies to formulate the process and allow for verification and validation of the results acquired. The close proximity of the results obtained validated the accuracy of the FEM analysis of the 2D axisymmetric model and 3D model. This result demonstrated the validity of commercial software instead of proprietary in-house code that is more costly to develop and maintain.

• Bulletin of the Korean Space Science Society
• /
• 2003.10a
• /
• pp.63-63
• /
• 2003

There are many actuators and sensors used for attitude control system for KOMPSAT such as Reaction Wheel Assembly, Magnetic Torque Assembly, Dual Thruster Module, Solar array Drive, Three Axis Magnetometer, Conical Earth Sensor, Fine Sun Sensor Assembly, Coarse Sun Sensor Assembly, Gyro Reference Assembly and so on. For KOMPSA T satellite it has been considered using the Magnetic Torquer (MTQ) generating the magnetic dipole moment. In general, the magnetic dipole moment for satellite attitude control system is used for dumping out the excessive reaction wheel momentum so that the reaction wheel speed is not saturated. The objective of this study is to analyze the magnetic field characteristics generated by the Magnetic Torquer using the Maxwell 2D Field Simulator software. Currently, the developing model (DM) of the MTQ is being developed and manufactured at a company under the supervision of KARL MTQ is an electromagnet consisting of a ferromagnetic cylindrical core on which an excitation coil is wound. A current is passed through the coil to produce a dipole momentum in the ferromagnetic core. The configuration of the MTQ will be introduced in the presentation. The 2 dimensional model of the MTQ is drawn as axisymmetric models in RZ plane, and each corresponding material is assigned to the each MTQ object, the core, coil, and background. After the boundary conditions, current sources, and solution parameters are set up, the magnetic field intensities, directions, and other values specified by users can be calculated by using the finite element analysis. The theoretical magnetic field quantities obtained by the Maxwell 2D Simulator can be used for the basis of the development of the MTQ.

• Journal of the Korea Convergence Society
• /
• v.9 no.1
• /
• pp.283-289
• /
• 2018

In automotive disc brake systems, frictional heat is not uniformly dispersed for reasons such as heat flux and thermal deformation. The thermoelastic deformation due to the frictional heat affects the contact pressure distribution and the contact load may be concentrated on the contact portion on the the disc brake surface, resulting in thermoelastic instability. In this study, thermal analysis and thermal deformation analysis considering the contact between disk and pad occurred during braking through 3D axial symmetry model with reference to the experimental equation and Kao's analysis method of contact pressure of disk and pad. ANSYS is used to analyze the thermal and elastic instability problems occurring at the contact surface between the disk and the pad, considering both the thermal and mechanical loads. A 3D axisymmetric model with direct contact between the disk and the pad was constructed to more accurately observe the thermal behavior of the disk by observing the frictional surface temperature, thermal deformation and contact thermal stress of the disk.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.24 no.3
• /
• pp.233-239
• /
• 2004

Smart material is used in various applications such as for glass frame, for medical instruments and for a part of sensors. Smart composite materials ran be applied to a part of aircraft and to the on-line monitoring system for industrial structures, using the shape memory effect. However, it is very difficult to simulate and analyze the shape memory effect in smart composites. In this paper, a two dimensional axisymmetric model was proposed to analyze the smart composite of one fiber and matrix using the finite element method(FEM). The finite element analysis was carried out in two renditions of the room temperature(293K) and a higher temperature (363K). The results we.e compared with the experimental results to confirm the validity of the analysis. In addition, the acoustic emission(AE) technique was used to study the microscopic damage behavior and the effect of pre-strains on TiNi/A16061 shape memory alloy composite.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.3
• /
• pp.247-255
• /
• 2013

The flight control of re-entry vehicles poses a challenge to conventional gain-scheduled flight controllers due to the widely spread aerodynamic coefficients. In addition, a wide range of uncertainties in disturbances must be accommodated by the control system. This paper presents the design of a roll channel controller for a non-axisymmetric reentry vehicle model using the trajectory linearization control (TLC) method. The dynamic equations of a moving mass system and roll control model are established using the Lagrange method. Nonlinear tracking and decoupling control by trajectory linearization can be viewed as the ideal gain-scheduling controller designed at every point along the flight trajectory. It provides robust stability and performance at all stages of the flight without adjusting controller gains. It is this "plug-and-play" feature that is highly preferred for developing, testing and routine operating of the re-entry vehicles. Although the controller is designed only for nominal aerodynamic coefficients, excellent performance is verified by simulation for wind disturbances and variations from -30% to +30% of the aerodynamic coefficients.