• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.10 s.181
• /
• pp.2552-2559
• /
• 2000

Most of ferrous b.c.c weld materials have martensitic transformation during rapid cooling after welding. It is well known that volume expansion due to the phase transformation could influence on the relaxation of welding residual stress. To apply this effect practically, it is necessary to establish a numerical model which is able to estimate the effect of phase transformation on residual stress relaxation quantitatively. For this purpose, the analysis is carried out in two regions, i.e., heating and cooling, because the variation of material properties following a phase transformation in cooling is different in comparison with the case in heating, even at the same temperature. The variation of material properties following phase transformation is considered by the adjustment of specific heat and thermal expansion coefficient, and the distribution of residual stress in analysis is compared with that of experiment by previous study. In this study, simplified numerical procedures considering phase transformation, which based on a commercial finite element package was established through comparing with the experimental data of residual stress distribution by other researcher. To consider the phase transformation effect on residual stress relaxation, the transition of mechanical and thermal property such as thermal expansion coefficient and specific heat capacity was found by try and error method in this analysis. In addition to, since the transformation temperature changes by the kind and control of alloying elements, the steel with many kinds of transformation temperature were selected and the effect of transformation on stress releasement was investigated by the numerical procedures considering phase transformation.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.37 no.2
• /
• pp.187-196
• /
• 2013

A numerical model was developed to predict the density wave oscillation (DWO) in the furnace wall tubes of a fossil-fired once-through boiler. The transient flow fields in the tubes were obtained using a 1D finite volume method in the time domain. A header model was also implemented to simulate the parallel tube connection of the wall tubes. The inlet and outlet mass flow variation in one of the parallel tubes was examined after a heat perturbation to find the DWO. After successful verification with experimental results reported in literature, the developed model was applied to the wall tubes of a 700-MW boiler furnace. In contrast to the simulation of Takitani's experiment, in which the unstable power thresholds tended to rise in the reduced bypass channel flow, no remarkable changes were observed in the power thresholds in the parallel channel modeling of the wall tubes of the boiler furnace.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.7 no.4
• /
• pp.33-38
• /
• 2003

The Physics of the Plume-induced shock and separation Particularly at a high Plume to exit pressure ratio and supersonic speeds up to Mach 3.0 with and without a passive control method, porous extension, were studied using computational techniques. Mass-averaged Navier-Stokes equations with the RNG $\kappa$-$\varepsilon$ turbulence model were solved using a fully implicit finite volume scheme and a 4-stage Runge-Kutta method. The control methodology for plume-afterbody interactions is to use a perforated wall attached at either the nozzle exit or the edge of the missile base. The Effect of porous wall length on plume interference is also investigated The computational results show the main effect of the porous extension on plume-afterbody interactions is to restrain the plume from strongly underexpanding during a change in flight conditions. With control, a change in porous extension length has no significant effect rut plume interference.

• Steel and Composite Structures
• /
• v.35 no.1
• /
• pp.77-92
• /
• 2020

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

• Nuclear Engineering and Technology
• /
• v.19 no.1
• /
• pp.34-41
• /
• 1987

This work concerns with a comparison of the 1-dimensional (or 1-D) load follow analysis method with reference to the detailed 3-dimensional (or 3-D) computations. For this purpose a 1-D two-group finite difference code, HLOFO, and a 3-D coarse-mesh code based on the modified Borresen's method, CMSNAC, are developed. The CMSNAC code is used to obtain the 3-D power peaks and reactivity parameters in response to power swing from 100 to 50 and back to 100% in the 12-3-6-3 load cycle for the BOL of the KORI Unit 1 PWR core. The 3-D result is then compared with the 1-D HLOFO computations, the cross section and buckling inputs of which are obtained by combining the flux-volume weighting scheme with the approximate flux from the auxiliary 3-D computations. It is shown that the 1-D computation has a limited accuracy, yet it is confirmed that it can describe the core axial average behavior which is fairly consistent with the detailed 3-D computation.

• Journal of Energy Engineering
• /
• v.5 no.1
• /
• pp.50-55
• /
• 1996

A numerical simulation on a primary calcinator of porcelain was performed with using Fluent to calculate the heat efficiency by studying velocity vector and temperature profile according to variables such as the location of outlet and porcelain. Control-Volume based Finite Difference Method and Up-wind scheme are used for discretization of differential equation. SIMPLEC Algorithm and standard k-$\varepsilon$ turbulent model are selected to resolve the pressure-velocity coupling and the turbulent. The result of simulation showed that the whole velocity vector field in a calcinator was varied greatly according to the location of outlet. But the whole temperature profile at each zone was still high regardless of the location of outlet because of the radiation. But the temperature of a case with a outlet at sidepart of preheating or cooling zone was little high compared to the case with a outlet on the top of preheating zone. The velocity vector field and temperature profile in a calcinator were almost not affected by the location of porcelain, but the temperature inside a porcelain was much affected according to the place where it was located. The heat efficiency in a calcinator was 44.6% and the gas temperature in the outlet was about 1000 K.

• Tribology and Lubricants
• /
• v.30 no.5
• /
• pp.303-310
• /
• 2014

In order to reduce friction and improve reliability, researchers have applied various surface texturing methods to highly sliding machine elements such as mechanical seals and piston rings. Despite extensive theoretical research on surface texturing, previous numerical results are only applicable to isothermal and iso-viscous conditions. Because the lubricant flow pattern of textured bearing surfaces is much more complicated than that for non-textured bearings, the Navier?Stokes equation is more suitable than the Reynolds equation for the former. This study carries out a thermohydrodynamic (THD) lubrication analysis to investigate the lubrication characteristics of a single micro-dimpled parallel thrust bearing cell. The analysis involves using the continuity, Navier?Stokes, energy, temperature?viscosity relation, and heat conduction equations with the commercial computational fluid dynamics (CFD) code FLUENT. This study discretizes these equations using the finite volume method and solves them using the SIMPLE algorithm. The results include finding the streamlines, pressure and temperature distributions, and variations in the friction force and leakage for various dimple radii and depths. Increasing the dimple radius and decreasing the depth causes a recirculation flow to form because of a strong vortex, and the oil temperature greatly increases compared with the non-textured case. The present numerical scheme and results are applicable to THD analysis of various surface-textured sliding bearings and can lead to further study.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.615-636
• /
• 2014

Aluminium Foam Sandwich (AFS) panels are becoming always more attractive in transportation applications thanks to the excellent combination of mechanical properties, high strength and stiffness, with functional ones, thermo-acoustic isolation and vibration damping. These properties strongly depend on the density of the foam, the morphology of the pores, the type (open or closed cells) and the size of the gas bubbles enclosed in the solid material. In this paper, the vibrational performances of two classes of sandwich panels with an Alulight(R) foam core are studied. Experimental tests, in terms of frequency response function and modal analysis, are performed in order to investigate the effect of different percentage of porosity in the foam, as well as the effect of the random distribution of the gas bubbles. Experimental results are used as a reference for developing numerical models using finite element approach. Firstly, a sensitivity analysis is performed in order to obtain a limit-but-bounded dynamic response, modelling the foam core as a homogeneous one. The experimental-numerical correlation is evaluated in terms of natural frequencies and mode shapes. Afterwards, an update of the previous numerical model is presented, in which the core is not longer modelled as homogeneous. Mass and stiffness are randomly distributed in the core volume, exploring the space of the eigenvectors.

• Wind and Structures
• /
• v.17 no.6
• /
• pp.647-670
• /
• 2013

The paper presents a numerical approach to study of fluid flow characteristics and to predict performance of wind turbines. The numerical model is based on Finite-volume method (FVM) discretization of unsteady Reynolds-averaged Navier-Stokes (URANS) equations. The movement of turbine blades is modeled using moving mesh technique. The turbulence is modeled using commonly used turbulence models: Renormalization Group (RNG) k-${\varepsilon}$ turbulence model and the standard k-${\varepsilon}$ and k-${\omega}$ turbulence models. The model is validated with the experimental data over a large range of tip-speed to wind ratio (TSR) and blade pitch angles. In order to demonstrate the use of numerical method as a tool for designing wind turbines, two dimensional (2-D) and three-dimensional (3-D) simulations are carried out to study the flow through a small scale Darrieus type H-rotor Vertical Axis Wind Turbine (VAWT). The flows predictions are used to determine the performance of the turbine. The turbine consists of 3-symmetrical NACA0022 blades. A number of simulations are performed for a range of approaching angles and wind speeds. This numerical study highlights the concerns with the self-starting capabilities of the present VAWT turbine. However results also indicate that self-starting capabilities of the turbine can be increased when the mounted angle of attack of the blades is increased. The 2-D simulations using the presented model can successfully be used at preliminary stage of turbine design to compare performance of the turbine for different design and operating parameters, whereas 3-D studies are preferred for the final design.

• 한국전산유체공학회:학술대회논문집
• /
• 1996.05a
• /
• pp.58-67
• /
• 1996

This paper describes some computational results of various energy and environmental systems using Patankar's SIMPLE method. The specific topics handled in this study are jet bubbling reactor for flue gas desulfurization, cyclone-type afterburner for incineration, 200m tall stack for 500 MW electric power generation, double skin and heat storage systems of building energy saving for the utilization of solar heating, finally turbulent combustion systems with liquid droplet or pulverized coal particle. A control-volume based finite-difference method with the power-law scheme is employed for discretization. The pressure-velocity coupling is resolved by the use of the revised version of SIMPLE, that is, SIMPLEC. Reynolds stresses are closed using the standard $k-{\varepsilon}$ and RNG $k-{\varepsilon}$ models. Two-phase turbulent combustion of liquid drop or pulverized coal particle is modeled using locally-homogeneous, gas-phase, eddy breakup model. However simple approximate models are incorporated for the modeling of the second phase slip and retardation of ignition without consideration of any detailed particle behavior. Some important results are presented and discussed in a brief note. Especially, in order to make uniform exit flow for the jet bubbling reactor, a well-designed structure of distributor is needed. Further, the aspect ratio in the double skin system appears to be one of important factors to give rise to the visible change of the induced air flow rate. The computational tool employed in this study, in general, appears as a viable method for the design of various engineering system of interest.