• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.5
• /
• pp.1150-1160
• /
• 1993

The behavior of polystyrene in the strip cavity during the packing and cooling stages for an injection molding process is examined numerically. The mathematical model is based on the unified post-filling model and finite element/finite difference methods are used to solve simultaneously the continuity, momentum and energy equations coupled to an equation of state. Simulated results show that the density of the molded parts is lower in the core than at the skin, and that the hotter the melt or the higher the packing pressure, the higher the density in the core. The density variation during the packing stage comes up to 50% compared with the total density variation. Also, the density variation after gate sealing and the effect of cooling rate on the equation of state are negligible.

• International Journal of Automotive Technology
• /
• v.7 no.7
• /
• pp.867-872
• /
• 2006

Fuel cells are devices that convert chemical energy directly into electrical energy. Owing to the high efficiency of the fuel cells, a large number of research work have been done during these years. Among many kinds of the fuel cells, a polymer electrolyte membrane fuel cell is such kind of thing which works under low temperature. Because of the specialty, it stimulated intense global R&D competition. Most of the major world automakers are racing to develop polymer electrolyte membrane fuel cell passenger vehicles. Unfortunately, there are still many problems to be solved in order to make them into the commercial use, such as the thermal and water management in working process of PEMFCs. To solve the difficulites facing the researcher, the analysis of the inner mechanism of PEMFC should be implemented as much as possible and mathematical modeling is an important tool for the research of the fuel cell especially with the combination of experiment. By regarding some of the assumptions and simplifications, using the finite element technique, a two-dimensional electrochemical mode is presented in this paper for the further comparison with experimental data. Based on the principals of the problem, the equations of electronic charge conservation equation, gas-phase continuity equation, and mass balance equation are used in calculating. Finally, modeling results indicate some of the phenomenon in a unit cell, and the relationships between potential and current density.

• Journal of the Society of Naval Architects of Korea
• /
• v.33 no.4
• /
• pp.48-59
• /
• 1996

The numerical step in the unsteady viscous flow analysis can be divided in the space analysis step satisfying continuity equation and the time marching step. In this study the spectral method is applied to solve the pressure Poisson equation in the space analysis step. If the highest order differential term of the pressure Poisson equation is transformed by Fourier series, pressure arid its first derivatives can be expressed by the integral form of Fourier series. So Gibb's phenomena can be eliminated and the spectral method can be applied to non-periodic problems. The numerical analysis of unsteady viscous flow around 2-dimensional circular cylinder and wing is carried out and compared for verification.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2002.10a
• /
• pp.25-32
• /
• 2002

Free vibration analysis of radialiy multi-delaminated beams with through-the-width multi-delamination is performed in the present study. The multiple delaminations are considered to be in a radial manner through the thickness from the top surface of the beam. The natural frequencies of the radially multi-delaminated beams are calculated from a new algorithm that is based on the single compound delaminated beam model. That is, beams with radial multi-delaminations are regarded as the sum of a single compound delamlnated beam that is the single sub-delaminated beam from the top surface of global beam. Each result of frequency equation for the single delaminated beam with unknown boundary conditions obtained through continuity conditions Is updated to the next one, With these sequential operations, the final frequency equation of radially multi-delaminated beams is obtained for both ends boundary conditions of global beam. The numerical results carried out for the beams are compared with those of some references to give the reliance on the proposed algorithm and to investigate the effects of the shape, number, size of multi-delaminations on the natural frequency. Compared with the other previously presented model, the proposed algorithm is more flexible in modeling and formulating as the total array size of frequency equation is always four by four. Therefore, the proposed algorithm will reduce the effort of user in formulating the physical model to the numerical model.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.33 no.2
• /
• pp.69-78
• /
• 2009

Cavitating flow simulation is of practical importance for many hydraulic engineering systems, such as pump, turbine, nozzle, injector, etc. In the present work, a solver for cavitating flow has been developed and applied to simulate the flows past axisymmetric cylinders. Governing equations are the two-phase Navier-Stokes equations, comprised of continuity equation of liquid and vapor phase. The momentum equation is in the mixture phase. The solver employed an implicit, dual time, preconditioned algorithm in curvilinear coordinates. Computations were carried out for three axisymmetric cylinders: hemispherical, ogive, and caliber-0 forebody shape. Then, the present calculations were compared with experiments and other numerical results to validate the present solver. Also, the code has shown its capability to accurately simulate the re-entrant jet phenomena and ventilated cavitation. Hence, it has been found that the present numerical code has successfully accounted for cavitating flows past axisymmetric cylinders.

• Journal of the Korean Mathematical Society
• /
• v.58 no.3
• /
• pp.553-569
• /
• 2021

In this paper, we present and analyze a fully discrete numerical method for solving the time-fractional diffusion wave equation: ∂^β_tu - div(a∇u) = f, 1 < β < 2. We first construct a difference formula to approximate ∂^β_tu by using an interpolation of derivative type. The truncation error of this formula is of O(△t^2+δ-β)-order if function u(t) ∈ C^2,δ[0, T] where 0 ≤ δ ≤ 1 is the Hölder continuity index. This error order can come up to O(△t^3-β) if u(t) ∈ C³ [0, T]. Then, in combinination with the linear finite volume discretization on spatial domain, we give a fully discrete scheme for the fractional wave equation. We prove that the fully discrete scheme is unconditionally stable and the discrete solution admits the optimal error estimates in the H¹-norm and L₂-norm, respectively. Numerical examples are provided to verify the effectiveness of the proposed numerical method.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.49 no.6
• /
• pp.449-456
• /
• 2021

In this paper, we present a strategy to extend solution capability of an existing low Mach number preconditioned compressible solver to incompressible flows with a little modification. To this end, the energy equation that is of the same form of the total energy equation of compressible flows is used. The energy equation is obtained by a linear combination of the thermal energy equation, the continuity equation and the mechanical energy equation. Subsequently, a modified artificial compressibility method in conjunction with a time marching technique is applied to these incompressible governing equations for steady flow solutions. It is found that the Roe average of the common governing equations is equally valid for both the compressible and incompressible flow conditions. The extension of an existing compressible solver to incompressible flows does not affect the original compressible flow analysis. Validity for incompressible flow analysis of the extended solver is examined for various inviscid, laminar and turbulent flows.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.2
• /
• pp.355-365
• /
• 1994

An analytical model of return flow is presented outside the surf zone. The governing equation is derived from the Navier-Stokes equation and the continuity. Each term of the governing equation is evaluated by the ordering analysis. Then the infinitesimal terms, i.e. the turbulent normal stress, the squared vertical velocity of water particle and the streaming velocity, are neglected. The driving forces of return flow are calculated using the linear wave theory for the shallow water approximation. Especially, the space derivative of local wave heights is described considering a shoaling coefficient. The vertical distribution of eddy viscosity is discussed to the customary types which are the constant, the linear function and the exponential function. Each coefficient of the eddy viscosities which sensitively affect the precision of solutions is uniquely decided from the additional boundary condition which the velocity becomes zero at the wave trough level. Also the boundary conditions at the bottom and the continuity relation are used in the integration of the governing equation. The theoretical solutions of present model are compared with the various experimental results. The solutions show a good agreement with the experimental results in the case of constant or exponential function type eddy viscosity.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.14 no.3
• /
• pp.997-1002
• /
• 2013

The article is aiming to investigate the thermal-fluidic characteristics of magnetic fluid in a cavity using GSMAC (generalized-simplified marker and cell method). The transport equations of the magnetic fluid are including the continuity equation, momentum equation and energy equation for natural convection and Maxwell equation and magnetization equation of magnetite nano-sized particles motion. In addition, the heat transfer characteristics such as temperatures and Nusselt numbers and flow characteristics such as streamlines and isotherms of the magnetic fluid were analyzed with the intensity and direction of the magnetic fields. As a result, the thermal-fluidic characteristics of the magnetic fluid in a cavity were could be controlled by the intensity and direction of the magnetic fields.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.4
• /
• pp.775-786
• /
• 2002

A Mixed Volume and Boundary Integral Equation Method is applied for the effective analysis of elastic wave scattering problems and plane elastostatic problems in unbounded solids containing general anisotropic inclusions and voids or isotropic inclusions. It should be noted that this newly developed numerical method does not require the Green's function for anisotropic inclusions to solve this class of problems since only Green's function for the unbounded isotropic matrix is involved in their formulation for the analysis. This new method can also be applied to general two-dimensional elastodynamic and elastostatic problems with arbitrary shapes and number of anisotropic inclusions and voids or isotropic inclusions. In the formulation of this method, the continuity condition at each interface is automatically satisfied, and in contrast to finite element methods, where the full domain needs to be discretized, this method requires discretization of the inclusions only. Finally, this method takes full advantage of the pre- and post-processing capabilities developed in FEM and BIEM. Through the analysis of plane elastostatic problems in unbounded isotropic matrix with orthotropic inclusions and voids or isotropic inclusions, and the analysis of plane wave scattering problems in unbounded isotropic matrix with isotropic inclusions and voids, it will be established that this new method is very accurate and effective for solving plane wave scattering problems and plane elastic problems in unbounded solids containing general anisotropic inclusions and voids/cracks or isotropic inclusions.