• Journal of Ocean Engineering and Technology
• /
• 제26권6호
• /
• pp.7-13
• /
• 2012

This paper presents a numerical study of the wave loads acting on offshore structures using a Cartesian-grid-based flow simulation method. Finite volume discretization with a volume-of-fluid (VOF) method is adopted to solve two-phase Navier-Stokes equations. Among the many variations of the VOF method, the CICSAM scheme is applied. The body boundary conditions are satisfied using a porosity function, and wave generation is carried out by using transient (wave or damping) zone approaches. In order to validate the present numerical method, three different basic offshore structures, including a sphere, Pinkster barge, and Wigley model, are numerically investigated. First, diffraction and radiation problems are solved using the present numerical method. The wave exciting and drift forces from the diffraction problems are compared with potential-based solutions. The added mass and wave damping forces from the radiation problems are also compared with the potential results. Next, the wave-induced motion responses of the structures are calculated and compared with the existing experimental data. The comparison results are fairly good, showing the validity of the present numerical method.

• Journal of Energy Engineering
• /
• 제2권3호
• /
• pp.268-275
• /
• 1993

A series of parametric investigations are performed in order to resolve the flow characteristic of a Quebec city stoker incinerator. The parameters considered in this study are five internal configurations of the Quebec city stoker itself and its modified ones, primary air velocity, the injection velocity and angle of the secondary air, and the reduction of the stoker exit area. A control-volume based finite-difference method by Patankar together with the power-law scheme is employed for discretization. The resolution of the pressure-velocity coupling is made by the use of SIMPLEC algorithm. The standard, two equation, k-$\varepsilon$ model is incorporated for the closure of turbulence. The size of recirculation region, turbulent viscosity, the mass fraction of the secondary air and pressure drop are calculated in order to analyze the characteristics of flow field. The results are physically acceptable and discussed in detail. The flow field of the Quebec city stoker shows the strong recirculation zone together with the high turbulence intensity over the upper part of the incinerator.

• Journal of Energy Engineering
• /
• 제5권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.

• Nuclear Engineering and Technology
• /
• 제51권6호
• /
• pp.1487-1503
• /
• 2019

The methods and performance of a pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT employing the drift-flux model are presented. This code aims at providing an accurate yet fast core thermal-hydraulics solution capability to high-fidelity multiphysics core analysis systems targeting massively parallel computing platforms. The four equation drift-flux model is adopted for two-phase calculations, and numerical solutions are obtained by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models involving turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analyses. ESCOT is parallelized by a domain decomposition scheme that involves both radial and axial decomposition to enable highly parallelized execution. The ESCOT solutions are validated through the applications to various experiments which include CNEN $4{\times}4$, Weiss et al. two assemblies, PNNL $2{\times}6$, RPI $2{\times}2$ air-water, and PSBT covering single/two-phase and unheated/heated conditions. The parameters of interest for validation include various flow characteristics such as turbulent mixing, spacer grid pressure drop, cross-flow, reverse flow, buoyancy effect, void drift, and bubble generation. For all the validation tests, ESCOT shows good agreements with measured data in the extent comparable to those of other subchannel-scale codes: COBRA-TF, MATRA and/or CUPID. The execution performance is examined with a mini-sized whole core consisting of 89 fuel assemblies and for an OPR1000 core. It turns out that it is about 1.5 times faster than a subchannel code based on the two-fluid three field model and the axial domain decomposition scheme works as well as the radial one yielding a steady-state solution for the OPR1000 core within 30 s with 104 processors.

• Steel and Composite Structures
• /
• 제38권5호
• /
• pp.477-496
• /
• 2021

The main objective of this paper is to study vibration of sandwich open cylindrical panel with damaged core and FG face sheets based on three-dimensional theory of elasticity. The structures are made of a damaged isotropic core and two external face sheets. These skins are strengthened at the nanoscale level by randomly oriented Carbon nanotubes (CNTs) and are reinforced at the microscale stage by oriented straight fibers. These reinforcing phases are included in a polymer matrix and a three-phase approach based on the Eshelby-Mori-Tanaka scheme and on the Halpin-Tsai approach, which is developed to compute the overall mechanical properties of the composite material. Three complicated equations of motion for the panel under consideration are semi-analytically solved by using 2-D differential quadrature method. Several parametric analyses are carried out to investigate the mechanical behavior of these multi-layered structures depending on the damage features, through-the-thickness distribution and boundary conditions. It is seen that for the large amount of power-law index "P", increasing this parameter does not have significant effect on the non-dimensional natural frequency parameters of the FG sandwich curved panel. Results indicate that by increasing the value of isotropic damage parameter "D" up to the unity (fully damaged core) the frequency would tend to become zero. One can dictate the fiber variation profile through the radial direction of the sandwich panel via the amount of "P", "b" and "c" parameters. It should be noticed that with increase of volume fraction of fibers, the frequency parameter of the panels does not increase necessarily, so by considering suitable amounts of power-law index "P" and the parameters "b" and "c", one can get dynamic characteristics similar or better than the isotropic limit case for laminated FG curved panels.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제27권4호
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• Polymer(Korea)
• /
• 제29권6호
• /
• pp.549-556
• /
• 2005

Effects of the dimensionless variables on the heat transport phenomena in the extrusion process of a single screw extruder have been studied numerically. Based on the understanding of the solids conveying related to the geometrical structure and characteristics of the screw, the heat balance equation for the solids conveying zone was established and normalized. The finite volume method and power-law scheme were applied to derive a discretized equation and the equation was solved using the alternating direction iterative method with relaxation. Effects of the dimensionless parameters, Biot and Peclet numbers, that define the heat transfer characteristics of the solids conveying zone have been investigated with respect to the temperature of the feeding zone and the length of the solids conveying zone. As the Biot number is increased, the heat loss by cooling dominates to decrease the temperature of the barrel but it has little effects on the temperature of the solids bed and the length of the solids conveying zone. On the other hand, if the Peclet number is increased, the convection term dominates to decrease the temperature of the solids bed and it results in an increase in the length of the solids conveying zone.

• International Journal of Aeronautical and Space Sciences
• /
• 제14권1호
• /
• pp.30-45
• /
• 2013

The present work focuses on the unsteady aerodynamics and aeroelastic properties of a small-medium sized wind-turbine blade operating under ideal conditions. A tapered/twisted blade representative of commercial blades used in an experiment setup at the National Renewable Energy Laboratory is considered. The aerodynamic loads are computed using Computational Fluid Dynamics (CFD) techniques. For this purpose, FLUENT$^{(R)}$, a commercial finite-volume code that solves the Navier-Stokes and the Reynolds-Averaged Navier-Stokes (RANS) equations, is used. Turbulence effects in the 2D simulations are modeled using the Wilcox k-w model for validation of the CFD approach. For the 3D aerodynamic simulations, in a first approximation, and considering that the intent is to present a methodology and workflow philosophy more than highly accurate turbulent simulations, the unsteady laminar Navier-Stokes equations were used to determine the unsteady loads acting on the blades. Five different blade pitch angles were considered and their aerodynamic performance compared. The structural dynamics of the flexible wind-turbine blade undergoing significant elastic displacements has been described by a nonlinear flap-lag-torsion slender-beam differential model. The aerodynamic quasi-steady forcing terms needed for the aeroelastic governing equations have been predicted through a strip-theory based on a simple 2D model, and the pertinent aerodynamic coefficients and the distribution over the blade span of the induced velocity derived using CFD. The resulting unsteady hub loads are achieved by a first space integration of the aeroelastic equations by applying the Galerkin's approach and by a time integration using a harmonic balance scheme. Comparison among two- and three- dimensional computations for the unsteady aerodynamic load, the flap, lag and torsional deflections, forces and moments are presented in the paper. Results, discussions and pertinent conclusions are outlined.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• 제23권8호
• /
• pp.1040-1047
• /
• 1999

The characteristics of flame propagation in inert particle-laden $H_2$/Air premixed gas are numerically investigated on this study. The 2nd order TVD scheme is applied to numerical analysis of governing equations and multi-step chemical reaction model and detailed transport properties are sued to solve chemical reaction terms. Radiation heat transfer is computed by applying the finite volume method to a radiative transfer equation. The burning velocities against the mole fractions of hydrogen agree well with results performed by different workers. The inert particles play significant roles in the flame propagation on account of momentum and heat transfer between gas and particles. Gas temperature, pressure and flame propagation speed are decreased as the loading ratio of particle is increased. Also the products behind flame zone contain lots of water vapor whose absorption coefficient is much larger than that of unburned gas. Thus, the radiation effect of gas and particles must be considered simultaneously for the flame propagation in a mixture of $H_2$/Air and inert particles. As a result, it is founded that because the water vapor emits much radiation and this emitted radiation is released at boundaries as radiant heat loss as well as reabsorbed by gas and particles, flame propagation speed and flame structure are altered with radiation effect.

• 한국전산유체공학회:학술대회논문집
• /
• 한국전산유체공학회 2004년도 춘계 학술대회논문집
• /
• pp.146-153
• /
• 2004

The Numerical study has been carried out to investigate the effects of chemical reaction and thermal radiation on the rocket plume flow-field at various altitudes. The theoretical formulation is based on the Navier-Stokes equations for compressible flows along with the infinitely fast chemistry and thermal radiation. The governing equations were solved by a finite volume fully-implicit TVD(Total Variation Diminishing) code which uses Roe's approximate Riemann solver and MUSCL(Monotone Upstream-centered Schemes for Conservation Laws) scheme. LU-SGS (Lower Upper Symmetric Gauss Seidel) method is used for the implicit solution strategy. An equilibrium chemistry module for hydrocarbon mixture with detailed thermo-chemical properties and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. In this study, kerosene-fueled rocket was assumed operating at O/F ratio of 2.34 with a nozzle expansion ratio of 6.14. Flight conditions considered were Mach number zero at ground level, Mach number 1.16 at altitude 5.06km and Mach number 2.9 at altitude 17.34km. Numerical results gave the understandings on the detailed plume structures at different altitude conditions. The diffusive effect of the thermal radiation on temperature field and the effect of chemical recombination during the expansion process could be also understood. By comparing the results from frozen flow and infinitely fast chemistry assumptions, the excess temperature of the exhaust gas resulting from the chemical recombination seems to be significant and cannot be neglected in the view point of performance, thermal protection and flow physics.