• 한국전산유체공학회:학술대회논문집
• /
• 2004.03a
• /
• 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.

• Journal of the Korean Geotechnical Society
• /
• v.30 no.9
• /
• pp.67-75
• /
• 2014

Methane gas hydrate which is considered energy source for the next generation has an urgent need to develop reliable numerical simulator for coupled THM phenomena in the porous media, to minimize problems arising during the production and optimize production procedures. International collaborations to improve previous numerical codes are in progress, but they still have mismatch in the predicted value and unstable convergence. In this paper, FEM code for fully coupled THM phenomena is developed to analyze methane hydrate dissociation in the porous media. Coupled partial differential equations are derived from four mass balance equations (methane hydrate, soil, water, and hydrate gas), energy balance equation, and force equilibrium equation. Five main variables (displacement, gas saturation, fluid pressure, temperature, and hydrate saturation) are chosen to give higher numerical convergence through trial combinations of variables, and they can analyze the whole region of a phase change in hydrate bearing porous media. The kinetic model is used to predict dissociation of methane hydrate. Developed THM FEM code is applied to the comparative study on a Masuda's laboratory experiment for the hydrate production, and verified for the stability and convergence.

• Geomechanics and Engineering
• /
• v.36 no.4
• /
• pp.355-366
• /
• 2024

The original elastoplastic Hardening Soil model is formulated actually partly under hexagonal pyramidal Mohr-Coulomb failure criterion, and can be only used in specific stress paths. It must be completely generalized under Mohr-Coulomb criterion before its usage in engineering practice. A set of generalized constitutive equations under this criterion, including shear and volumetric yield surfaces and hardening laws, is proposed for Hardening Soil model in principal stress space. On the other hand, a Mohr-Coulumb type yield surface in principal stress space comprises six corners and an apex that make singularity for the normal integration approach of constitutive equations. With respect to the isotropic nature of the material, a technique for processing these singularities by means of Koiter's rule, along with a transforming approach between both stress spaces for both stress tensor and consistent stiffness matrix based on spectral decomposition method, is introduced to provide such an approach for developing generalized Hardening Soil model in finite element analysis code ABAQUS. The implemented model is verified in comparison with the results after the original simulations of oedometer and triaxial tests by means of this model, for volumetric and shear hardenings respectively. Results from the simulation of oedometer test show similar shape of primary loading curve to the original one, while maximum vertical strain is a little overestimated for about 0.5% probably due to the selection of relationships for cap parameters. In simulation of triaxial test, the stress-strain and dilation curves are both in very good agreement with the original curves as well as test data.

• International Journal of Aeronautical and Space Sciences
• /
• v.14 no.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.

• Proceedings of the Korean Powder Metallurgy Institute Conference
• /
• 1997.10a
• /
• pp.7-7
• /
• 1997

During sintering of very porous green bodies, as obtained by compaction of hard powders - such as tungsten carbide or ceramics - or by injection moulding, important shrinkage occurs. Due to heterogeneous green density field, gravity effects, friction on the support, thermal gradients, etc., this shrinkage is often non-uniform, which' may induce significant shape changes. As the ratio of compact dimension to powder size is very high, the mechanics of continuum is relevant to model such phenomena. Thus numerical techniques, such as the finite element method can be used to simulate the sintering process and predict the final shape of the sintered part. Such type of simulation has much been developed in the last decade firstly for hot isostatic pressing and next for die compaction. Finite element modelling has been recently applied to free sintering. The simulation of sintering should be based on constitutive equations describing the thermo-mechanical behaviour of the material under any state of stress and any temperature which may arise within the sintering body. These equations can be drawn either from experimental data or from micromechanical models. The experiments usually consist in free sintering and sinter-forging tests. Indeed applying more complex loading conditions at high temperature under controlled atmosphere is delicate. Micromechanical models describe the constitutive behaviour of aggregates of spheres from the deformation of two-sphere contact either by viscous flow or grain boundary diffusion. Such models are not able to describe complex microstructure and mechanisms as observed in real materials but they can give some basic information on the formulation of constitutive equations. Practically both experimental and theoretical approaches can be coupled to identify the constitutive equations. Such procedure has been performed for modelling the sintering of compacts obtained by die pressing of a mixture of tungsten carbide and cobalt powders. The constitutive behaviour of this material during sintering has been described by a linear viscous constitutive model, whose functions have been fitted from results of free sintering and sinter-forging experiments. This model has next been introduced in ABAQUS finite element code to simulate the sintering of heterogeneous green compacts of various geometries at constant temperature. Examples of simulations are shown and compared with experiments.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Progress in Medical Physics
• /
• v.14 no.1
• /
• pp.1-7
• /
• 2003

In this work we developed a EGS4 control code to calculate the dose distributions for high energy electron beams in water phantom applied longitudinal magnetic field. We reviewed the electron's motion in magnetic field and delivered equations for direction changs of the electron by the external magnetic field. The mathematical results are inserted into the EGS4 code system to account for the presence of external magnetic fields in phantom. The electron pencil beam paths of 6 MeV in water phantom are calculated for magnetic fields of 1-3 T and the dose distributions for a field of 1.0 cm in diameter are calculated for magnetic fields of 0.6-1 T using the code. From the results of path calculations we found that the lateral ranges of the electrons are reduced in the magnetic field of 3 T. For a field of 1 cm diameter and a magnetic field of 1 T, the small dose enhancement near the range of the electrons on the depth dose and the penumbra reduction of 0.15 cm on the beam profile are observed. We discussed and evaluated the results from the theoretical concepts.

• Journal of Ocean Engineering and Technology
• /
• v.33 no.5
• /
• pp.400-410
• /
• 2019

In this study, the numerical code for the 3D nonlinear dynamic analysis of an SLWR (Steel Lazy Wave Riser) was developed using the lumped mass line model in a FORTRAN environment. Because the lumped mass line model is an explicit method, there is no matrix operation. Thus, the numerical algorithm is simple and fast. In the lumped mass line model, the equations of motion for the riser were derived by applying the various forces acting on each node of the line. The applied forces at the node of the riser consisted of the tension, shear force due to the bending moment, gravitational force, buoyancy force, riser/ground contact force, and hydrodynamic force based on the Morison equation. Time integration was carried out using a Runge-Kutta fourth-order method, which is known to be stable and accurate. To validate the accuracy of the developed numerical code, simulations using the commercial software OrcaFlex were carried out simultaneously and compared with the results of the developed numerical code. To understand the nonlinear dynamic characteristics of an SLWR, dynamic simulations of SLWRs excited at the hang-off point and of SLWRs in regular waves were carried out. From the results of these dynamic simulations, the displacements at the maximum bending moments at important points of the design, like the hang-off point, sagging point, hogging points, and touch-down point, were observed and analyzed.

• Journal of the Korea Concrete Institute
• /
• v.16 no.1 s.79
• /
• pp.43-53
• /
• 2004

The current Korean Concrete Design Code(KCI Code) requires the minimum and maximum content of shear s in order to prevent brittle and noneconomic design. However, the required content of the steel reinforcement In KCI Code is quite different to those of the other design codes such as fib-code, Canadian Code, and Japanese Code. Furthermore, since the evaluation equations of the minimum and maximum shear reinforcement for the current KCI Code were based on the experimental results, the equations can not be used for the RC members beyond the experimental application limits. The concrete tensile strength, shear stress, crack inclination, strain perpendicular to the crack, and shear span ratio are strongly related to the lower and upper limits of shear reinforcement. In this research, an evaluation equation for the minimum content of shear reinforcement is theoretical proposed from the Wavier's three principals of the mechanics of materials.

• Journal of computational fluids engineering
• /
• v.20 no.2
• /
• pp.103-108
• /
• 2015

Fluid-body interaction analysis of floating body with six degree-of-freedom motion is presented. In this study, three-dimensional incompressible Navier-Stokes equations are employed as a governing equation. The numerical method is based on a finite-volume approach on a cartesian grid together with a fractional-step method. To represent the body motion, the immersed boundary method for direct forcing is employed. In order to simulate the coupled six degree-of-freedom motion, Euler's equations based on rigid body dynamics are utilized. To represent the complex body shape, level-set based algorithm is utilized. In order to describe the free surface motion, the volume of fluid method utilizing the tangent of hyperbola for interface capturing scheme is employed. This study showed three different continuums(air, water and body) are simultaneously simulated by newly developed code. To demonstrate the applicability of the current approach, two different problems(dam-breaking with stationary obstacle and water entry) are simulated and all results are validated.