• Journal of computational fluids engineering
• /
• v.18 no.2
• /
• pp.66-71
• /
• 2013

In this paper, the supersonic flows around space launch vehicles have been numerically simulated by using a 3-D RANS flow solver. The focus of the study was made for investigating plume-induced flow separation(PIFS). For this purpose, a vertex-centered finite-volume method was utilized in conjunction with 2nd-order Roe's FDS to discretize the inviscid fluxes. The viscous fluxes were computed based on central differencing. The Spalart-Allmaras model was employed for the closure of turbulence. The Gauss-Seidel iteration was used for time integration. To validate the flow solver, calculation was made for the 0.04 scale model of the Saturn-5 launch vehicle at the supersonic flow condition without exhaust plume, and the predicted results were compared with the experimental data. Good agreements were obtained between the present results and the experiment for the surface pressure coefficient and the Mach number distribution inside the boundary layer. Additional calculations were made for the real scale of the Saturn-5 configuration with exhaust plume. The flow characteristics were analyzed, and the PIFS distances were validated by comparing with the flight data. The KSLV-1 is also simulated at the several altitude conditions. In case of the KSLV-1, PIFS was not observed at all conditions, and it is expected that PIFS is affected by the nozzle position.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.11
• /
• pp.3538-3548
• /
• 1996

In this study, as a part of basic study for developing the simulation program for the assemssment of reliability of electronic EMC packaging which covers from EMC mixing step to thermal analysis, comparison between a measured and predicted values of material properties of EMC and finitde element analysis of thermal stress are performed. For the experimental testing specimens of fifty, sixty hive and eighty percent filler($13\mu m$, spherical silica) weight fraction are fabricated using tranfer molding. Coefficient of thermal expansion, elastic modulus and thermla conductivity are measured using these specimens and then these measured values are compared with the predicted values by various equations ( such as dilute suspension method. self consistent method, generalized self consistent method, Hashin-Shtrikman's bounds. Shapery's bounds, Nielsen's method and others). Measured values are distributed within the upper and lower bounds of equations. Measured elastic modulus and coefficient of thermal expansion approaches closely the perdicted values with self consisten mehtod and upper bound of Shaperys equation respectively. However small differences of thermal conductivity between the different filler volume fraction are obserbed. FEM analysis indicates that firstly stress is concentrated at the corner section of EMC and secondly EMC with eighty percent filler weight fraction shows less thermal stress when package is cooling down and relatively high thermal stress when package is heating up.

• 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.