• Transactions of the Korean Society of Mechanical Engineers
• /
• v.14 no.3
• /
• pp.716-724
• /
• 1990

Oblique impinging plane jets were investigated experimentally and numerically at Reynolds number 21000. The inclination angle was varied from 90.deg.(normal to the impinging plate) to 60.deg.. The distance H between the nozzle exit and the stagnation point on the impinging plate was fixed at H/D=8. The working fluid was air. The mean velocity components and turbulent quantities were measured by a hot-wire anemometer. And the static pressure distributions on the impinging plate were measured by a Pitot tube. In numerical computation, the governing partial differential equations of elliptic type were solved with conventional k-.epsilon. turbulence model. The measurements show that, after impingement, the jet half width alone the wall increases in both directions, and that similarity for each turbulent quantity such as Reynolds shear stress or turbulent kinetic energy is revealed in the wall jet region. The computed results show some deviation from experimental data in the impingement region, where streamline curvature is significant. However, the computed results agree qualitatively well with measurements.

• Journal of the Korean Society for Advanced Composite Structures
• /
• v.2 no.3
• /
• pp.25-29
• /
• 2011

Some laminate orientations have decreasing values of $D_{16}$, $B_{16}$, $D_{26}$ and $B_{26}$ stiffnesses as the ply number increases. For such plates, the fiber orientations given above behave as specially orthotropic plates and simple formulas developed by the senior author. Most of the bridge and building slabs on girders have large aspect ratios. For such cases further simplification is possible by neglecting the effect of the longitudinal moment terms(Mx) on the relevant partial differential equations of equilibrium. In this paper, the result of the study on the subject problem is presented.

• Structural Engineering and Mechanics
• /
• v.55 no.4
• /
• pp.719-742
• /
• 2015

This paper dealt the free vibration analysis of thick truncated conical composite sandwich shells with transversely flexible cores and simply supported boundary conditions based on a new improved and enhanced higher order sandwich shell theory. Geometries were used in the present work for the consideration of different radii curvatures of the face sheets and the core was unique. The coupled governing partial differential equations were derived by the Hamilton's principle. The in-plane circumferential and axial stresses of the core were considered in the new enhanced model. The first order shear deformation theory was used for the inner and outer composite face sheets and for the core, a polynomial description of the displacement fields was assumed based on the second Frostig's model. The effects of types of boundary conditions, conical angles, length to radius ratio, core to shell thickness ratio and core radius to shell thickness ratio on the free vibration analysis of truncated conical composite sandwich shells were also studied. Numerical results are presented and compared with the latest results found in literature. Also, the results were validated with those derived by ABAQUS FE code.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.1
• /
• pp.61-71
• /
• 1998

This paper describes a dynamic fracture behaviors of structural elements under elastic or elasto-plastic stress waves in two dimensional space. The governing equation of this problem has the type of hyperbolic partial differential equation, which consists of the equation of motions and incremental elasto-plastic constitutive equations. To solve this problem we introduce Zwas' method which is based on the finite difference method. Additionally, in order to deal with the dynamic behavior of elasto-plastic problems, an elasto-plastic loading path in the stress space is proposed to model the plastic yield phenomenon. Based on the result of this computation, the dynamic stress intensity factor at the crack tip of an elastic material is calculated, and the time history of a plastic zone of a elasto-plastic material is to be shown.

• Journal of the Korean Society for Precision Engineering
• /
• v.21 no.5
• /
• pp.73-82
• /
• 2004

This paper presents the study on damage diagnosis of an intelligent cantilevered beams using PZT actuator and PVDF sensor This study provides the theoretical and experimental verification to examine structural damage. Time domain analysis for the non-destructive detection of damage is presented by parameterized partial differential equations and Galerkin approximation techniques. The time histories of the vibration response of structure were used to identify the presence of damage. Furthermore, this systematic approach permits one to use the piezomaterials to both excite and sense the vibration of structures. We also carried out the experimental verification about reliability of theoretical methods fur detecting the damage of a composite beam with PZT actuator and PVDF sensor. Experimental results are presented from tests on cantilevered composite beams which is damaged at different location and different dimensions. The results were compared with the simulation results. Good agreement between the results was found for the time shifts and amplitude difference in transients response of the cantilevered beam.

• Coupled systems mechanics
• /
• v.8 no.2
• /
• pp.129-146
• /
• 2019

The main goal of this work is to develop a numerical simulator to study an isothermal single-phase two-component flow in a naturally fractured oil reservoir, taking into account advection and diffusion effects. We use the Peng-Robinson equation of state with a volume translation to evaluate the properties of the components, and the discretization of the governing partial differential equations is carried out using the Finite Difference Method, along with implicit and first-order upwind schemes. This process leads to a coupled non-linear algebraic system for the unknowns pressure and molar fractions. After a linearization and the use of an operator splitting, the Conjugate Gradient and Bi-conjugated Gradient Stabilized methods are then used to solve two algebraic subsystems, one for the pressure and another for the molar fraction. We studied the effects of fractures in both the flow field and mass transport, as well as in computing time, and the results show that the fractures affect, as expected, the flow creating a thin preferential path for the mass transport.

• Computers and Concrete
• /
• v.27 no.4
• /
• pp.385-393
• /
• 2021

Reinforcement corrosion is the main cause of the durability failure of reinforced concrete (RC) structure. In this paper, a three-dimensional (3D) numerical model of macro-cell corrosion is established to reveal the corrosion mechanisms of steel reinforcement in RC structure. Modified Direct Iteration Method (MDIM) is employed to solve the system of partial differential equations for reinforcement corrosion. Through the sensitivity analysis of electrochemical parameters, it is found that the average corrosion current density is more sensitive to the change of cathodic Tafel slope and anodic equilibrium potential, compared with the other electrochemical parameters. Furthermore, both the anode-to-cathode (A/C) ratio and the anodic length have significant influences on the average corrosion current density, especially when A/C ratio is less than 0.5 and anodic length is less than 35 mm. More importantly, it is demonstrated that the corrosion rate of semi-circumferential corrosion is much larger than that of circumferential corrosion for the same A/C ratio value. The simulation results can give a unique insight into understanding the detailed electrochemical corrosion processes of steel reinforcement in RC structure for application in service life prediction of RC structures in actual civil engineer.

• Steel and Composite Structures
• /
• v.48 no.4
• /
• pp.443-459
• /
• 2023

This paper investigates the stability of a bi-directional functionally graded (BD-FG) cylindrical beam made of imperfect concrete, taking into account size-dependency and the effect of geometry on its stability behavior. Both buckling and dynamic behavior are analyzed using the modified coupled stress theory and the classical beam theory. The BD-FG structure is created by using porosity-dependent FG concrete, with changing porosity voids and material distributions along the pipe radius, as well as uniform and nonuniform radius functions that vary along the beam length. Energy principles are used to generate partial differential equations (PDE) for stability analysis, which are then solved numerically. This study sheds light on the complex behavior of BD-FG structures, and the results can be useful for the design of stable cylindrical microstructures.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.32 no.1
• /
• pp.16-24
• /
• 2022

A computational study of combined thermal and solutal convection (double diffusive convection) in a sealed crystal growth reactor is presented, based on a two-dimensional numerical analysis of the nonlinear and strongly coupled partial differential equations and their associated boundary conditions. The average Nusselt numbers for the source regions are greater than those at the crystal regions for 9.73 × 10³ ≤ Gr_t ≤ 6.22 × 10⁵. The average Nusselt numbers for the source regions varies linearly and increases directly with the thermal Grashof number form 9.73 × 10³ ≤ Gr_t ≤ 6.22 × 10⁵ for aspect ratio, Ar (transport length-to-width) = 1 and 2. Additionally, the average Nusselt numbers for the crystal regions at Ar = 1 are much greater than those at Ar = 2. Also, the occurrence of one unicellular flow structure is caused by both the thermal and solutal convection, which is inherent during the physical vapor transport of Hg₂Br₂. When the aspect ratio of the enclosure increases, the fluid movement is hindered and results in the decrease of thermal buoyancy force.

• Ocean Systems Engineering
• /
• v.12 no.3
• /
• pp.359-384
• /
• 2022

We develop in this work a new well-balanced preserving-positivity path-conservative central-upwind scheme for Saint-Venant-Exner (SVE) model. The SVE system (SVEs) under some considerations, is a nonconservative hyperbolic system of nonlinear partial differential equations. This model is widely used in coastal engineering to simulate the interaction of fluid flow with sediment beds. It is well known that SVEs requires a robust treatment of nonconservative terms. Some efficient numerical schemes have been proposed to overcome the difficulties related to these terms. However, the main drawbacks of these schemes are what follows: (i) Lack of robustness, (ii) Generation of non-physical diffusions, (iii) Presence of instabilities within numerical solutions. This collection of drawbacks weakens the efficiency of most numerical methods proposed in the literature. To overcome these drawbacks a reformulation of the central-upwind scheme for SVEs (CU-SVEs for short) in a path-conservative version is presented in this work. We first develop a finite-volume method of the first order and then extend it to the second order via the averaging essentially non oscillatory (AENO) framework. Our numerical approach is shown to be well-balanced positivity-preserving and shock-capturing. The resulting scheme could be seen as a predictor-corrector method. The accuracy and robustness of the proposed scheme are assessed through a carefully selected suite of tests.