• Journal of the Society of Naval Architects of Korea
• /
• v.28 no.2
• /
• pp.159-173
• /
• 1991

This paper describes a potential-based panel method formulated for the analysis of a super-cavitating two-dimensional hydrofoil. The method employs normal dipoles and sources distributed on the foil and cavity surfaces to represent the potential flow around the cavitating hydrofoil. The kinematic boundary condition on the wetted portion of the foil surface is satisfied by requiring that the total potential vanish in the fictitious inner flow region of the foil, and the dynamic boundary condition on the cavity surface is satisfied by requiring thats the potential vary linearly, i.e., the tangential velocity be constant. Green's theorem then results in a potential-based integral equation rather than the usual velocity-based formulation of Hess & Smith type. With the singularities distributed on the exact hydrofoil surface, the pressure distributions are predicted with improved accuracy compared to those of the linearized lilting surface theory, especially near the leading edge. The theory then predicts the cavity shape and cavitation number for an assumed cavity length. To improve the accuracy, the sources and dipoles on the cavity surface are moved to the newly computed cavity surface, where the boundary conditions are satisfied again. This iteration process is repeated until the results are converged. Characteristics of iteration and discretization of the present numerical method are much faster and more stable than the existing nonlinear theories. The theory shows good correlations with the existing theories and experimental results for the super-cavitating flow. In the region of small angles of attack, the present prediction shows and excellent comparison with the Geurst's linear theory. For the long cavity, the method recovers the trends of the Wu's nonlinear theory. In the intermediate regions of the short super-cavitation, the method compares very well with the experimental results of Parkin and also those of Silberman.

• Proceedings of the Acoustical Society of Korea Conference
• /
• spring
• /
• pp.469-474
• /
• 2002

In recent years, modularization of engine parts has increased the application of plastic products in air intake systems. Plastic intake manifolds provide many advantages including reduced weight, contracted cost, and lower intake air temperatures. These manifolds, however, have some weakness when compared with customary aluminium intake manifolds, in that they have low sound transmission loss because of their lower material density. This low transmission loss of plastic intake manifolds causes several problems related to flow noise, especially when the throttle is opened quickly. The physical processes, responsible for this flow noise, include turbulent fluid motion and relative motion of the throttle to the airflow. The former is generated by high-speed airflow in the splits between the throttle valve and the inner-surface of the throttle body and surge-tank, which can be categorized into the quadrupole source. The latter induces the unsteady force on the flow, which can be classified into the dipole source. In this paper, the mechanism of noise generation from the turbulence is only investigated as a preliminary study. Stochastic noise source synthesis method is adopted for the analysis of turbulence-induced, i.e. quadrupole noise by throttle at quick opening state. The method consists of three procedures. The first step corresponds to the preliminary time-averaged Navier-Stokes computation with a $k-\varepsilon$ turbulence model providing mean flow field characteristics. The second step is the synthesis of time-dependent turbulent velocity field associated with quadrupole noise sources. The final step is devoted to the determination of acoustic source terms associated with turbulent velocity. For the first step, we used market available analysis tools such as STAR-CD, the trade names of fluid analysis tools available on the market. The steady state flows at three open angle of throttle valve, i.e. 20, 35 and 60 degree, are numerically analyzed. Then, time-dependent turbulent velocity fields are produced by using the stochastic model and the flow analysis results. Using this turbulent velocity field, the turbulence-originated noise sources, i.e. the self-noise and shear-noise sources are synthesized. Based on these numerical results, it is found that the origin of the turbulent flow and noise might be attributed to the process of formulation and the interaction of two vortex lines formed in the downstream of the throttle valve. These vortex lines are produced by the non-uniform splits between the throttle valve and inner cylinder surface. Based on the analysis, we present the low-noise design of the inner geometry of throttle body.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.1
• /
• pp.151-163
• /
• 2008

A formulation of three-dimensional model of articulated train-b ridge dynamic interaction has been made for the Korean eXpress Train (KTX). Semi-periodic profiles of rail irregularities consisting of elevation, alignment, cross and gauge irregularities have also been proposed using FRA maximum tolerable rail deviations. The effects of rail joints and sleeper step were also included. The resulting system matrices of train and bridge are very spare, and thus, are stored in one-dimensional arrays, yielding a time-efficient solution. A numerical algorithm for computing bridge-train response including an iterative scheme is also formulated. A program simulating train-bridge interaction and solving this problem using the new algorithm is implemented as new modules for the f inite element analysis software named XFINAS. Computed results using the new program are then checked by that of the validated 2-D bridge-train interaction model. This new 3D analysis provides more detailed train responses such as swaying, bouncing, rolling, pitching and yawing accelerations, which are useful inevaluating passenger riding comfort. Train operation safety and derailment could also be directly investigated by relative wheel displacements computed from this program.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.3
• /
• pp.116-126
• /
• 1993

The finite element method(FEM) has been commonly used for structural dynamic analysis. However, the direct global application of FEM to large complex structures such as ships and offshore structures requires considerable computational efforts, and remarkably more in structural dynamic optimization problems. Adoption of the component-mode synthesis method is an efficient means to overcome the above difficulty. Among three classes of the component-mode synthesis method, the free-interface mode method is recognized to have the advantages of better computational efficiency and easier implementation of substructures' experimental results, but the disadvantage of lower accuracy in analytical results. In this paper, an advanced method to improve the accuracy in the application of the free-interface mode method for the vibration analysis of large complex structures is presented. In order to compensate the truncation effect of the higher modes of substructures in the synthesis process, both residual inertia and stiffness effects are taken into account and a frequency shifting technique is introduced in the formulation of the residual compliance of substructures. The introduction of the frequency shrift ins not only excludes cumbersome manipulation of singular matrices for semi-definite substructural systems but gives more accurate results around the specified shifting frequency. Numerical examples of typical structural models including a ship-like two dimensional finite element model show that the analysis results based on the presented method are well competitive in accuracy with those obtained by the direst global FEM analysis for the frequencies which are lower than the highest one employed in the synthesis with remarkably higher computational efficiency and that the presented method is more efficient and accurate than the fixed-interface mode method.

• Journal of the Korea Society for Simulation
• /
• v.28 no.2
• /
• pp.1-14
• /
• 2019

A full scale hydrodynamic simulation that requires an accurate reproduction of shock-induced detonation was conducted for design of an energetic component system. A detailed hydrodynamic analysis SW was developed to validate the reactive flow model for predicting the shock propagation in a train configuration and to quantify the shock sensitivity of the energetic materials. The pyrotechnic device is composed of four main components, namely a donor unit (HNS+HMX), a bulkhead (STS), an acceptor explosive (RDX), and a propellant (BPN) for gas generation. The pressurized gases generated from the burning propellant were purged into a 10 cc release chamber for study of the inherent oscillatory flow induced by the interferences between shock and rarefaction waves. The pressure fluctuations measured from experiment and calculation were investigated to further validate the peculiar peak at specific characteristic frequency (${\omega}_c=8.3kHz$). In this paper, a step-by-step numerical description of detonation of high explosive components, deflagration of propellant component, and deformation of metal component is given in order to facilitate the proper implementation of the outlined formulation into a shock physics code for a full scale hydrodynamic simulation of the energetic component system.

• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.19 no.11
• /
• pp.205-212
• /
• 2018

In order to overcome shortcomings of commercial analysis program for solving certain geotechnical problems, finite element method adopting isoparametric quadrilateral element was selected as a tool for analyzing soil behavior and calculating process was programmed. Two examples were considered in order to verify reliability of the developed program. One of the two examples is the case of acting isotropic confining pressure on finite element and the other is the case of acting shear stress on the sides of the finite element. Isoparametric quadrilateral element was considered as the finite element and displacements in the element can be expressed by node displacements and shape functions in the considered element. Calculating process for determining strain which is defined by derivatives using global coordinates was coded using the Jacobian and the natural coordinates. Four point Gauss rule was adopted to convert double integral which defines stiffness of the element into numerical integration. As a result of executing analysis of the finite element under isotropic confining pressure, calculated stress corresponding to four Gauss points and center of the element were equal to the confining pressure. In addition, according to the analyzed results for the element under shear stress, horizontal stresses and vertical stresses were varied with positions in the element and the magnitudes and distribution pattern of the stresses were thought to be rational.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• Tunnel and Underground Space
• /
• v.32 no.5
• /
• pp.312-326
• /
• 2022

The generalized Hoek-Brown (GHB) criterion, which is recognized as one of the standard failure conditions for rock mass, is specialized for rock engineering applications and covers a wide range of rock mass conditions. Accordingly, many research efforts have been devoted to the incorporation of this criterion into the stability analysis of rock structures. In this study, the slip-line analysis method, which is a kind of elastoplastic analysis method, is combined with the GHB failure criterion to derive analytical equations that can easily calculate the plastic radius and stress distribution in the vicinity of the circular tunnel. In the process of derivation of related formulas, it is assumed that the behavior of rock mass after failure is perfectly plastic and the in-situ stress condition is hydrostatic. In the formulation, it is revealed that the plastic radius can be calculated analytically using the two respective tangential friction angles corresponding to the stress conditions at tunnel wall and elastic-plastic boundary. It is also shown that the plastic radius and stress distribution calculated using the derived analytical equations coincide with the results of Lee & Pietruszczak's numerical method published in 2008. In the latter part of this paper, the influence of the quality of the rock mass on the size of the plastic zone, the stress distribution, and the change of the tangential friction angle was investigated using the derived analytical equations.

• Advances in aircraft and spacecraft science
• /
• v.9 no.5
• /
• pp.415-431
• /
• 2022

Secondary flows have a huge impact on losses generation in modern low pressure gas turbines (LPTs). At design point, the interaction of the blade profile with the end-wall boundary layer is responsible for up to 40% of total losses. Therefore, predicting accurately the end-wall flow field in a LPT is extremely important in the industrial design phase. Since the inlet boundary layer profile is one of the factors which most affects the evolution of secondary flows, the first main objective of the present work is to investigate the impact of two different inlet conditions on the end-wall flow field of the T106A, a well known LPT cascade. The first condition, labeled in the paper as C1, is represented by uniform conditions at the inlet plane and the second, C2, by a flow characterized by a defined inlet boundary layer profile. The code used for the simulations is based on the Discontinuous Galerkin (DG) formulation and solves the Reynolds-averaged Navier-Stokes (RANS) equations coupled with the Spalart Allmaras turbulence model. Secondly, this work aims at estimating the influence of viscosity and turbulence on the T106A end-wall flow field. In order to do so, RANS results are compared with those obtained from an inviscid simulation with a prescribed inlet total pressure profile, which mimics a boundary layer. A comparison between C1 and C2 results highlights an influence of secondary flows on the flow field up to a significant distance from the end-wall. In particular, the C2 end-wall flow field appears to be characterized by greater over turning and under turning angles and higher total pressure losses. Furthermore, the C2 simulated flow field shows good agreement with experimental and numerical data available in literature. The C2 and inviscid Euler computed flow fields, although globally comparable, present evident differences. The cascade passage simulated with inviscid flow is mainly dominated by a single large and homogeneous vortex structure, less stretched in the spanwise direction and closer to the end-wall than vortical structures computed by compressible flow simulation. It is reasonable, then, asserting that for the chosen test case a great part of the secondary flows details is strongly dependent on viscous phenomena and turbulence.

• Structural Engineering and Mechanics
• /
• v.87 no.3
• /
• pp.283-296
• /
• 2023

Free-vibration and buckling analyses of plate problems are investigated with the aid of the strain gradient notation finite element method (SGN-FEM). As SGN-FEM employs physically interpretable polynomials in developing finite elements, parasitic shear sources, which are the cause of shear locking, can be precisely identified and subsequently eliminated. This allows two mutually complementary objectives to be defined in this work, namely, evaluate the efficiency of free-vibration and buckling results provided by corrected models, and study the severity of parasitic shear effects on plate models performance. Parasitic shear are flexural terms erroneously present in shear strain polynomials. It is reviewed here that six parasitic shear terms arise during the formulation of the four-node Mindlin plate element. Two parasitic shear terms have been identified in the in-plane shear strain polynomial while other two have been identified in each of the transverse shear strain polynomials. The element is corrected a-priori, i.e., during development, by simply removing the spurious terms from the shear strain polynomials. The computational implementation of the element in its two versions, namely, containing the parasitic shear terms (PS) and corrected for parasitic shear (SG), allows for assessments of the accuracy of results and of the deleterious effects of parasitic shear in free vibration and buckling analyses. This assessment of the parasitic shear effects is a novelty of this work. Validation of the SG model is done comparing its results with analytical results and results provided by other numerical procedures. Analyses are performed for square plates with different thickness-to-length ratios and boundary conditions. Results for thin plates provided by the PS model do not converge to the correct solutions, which indicates that parasitic shear must be eliminated. That is, analysts should not rely on refinement alone. For thick plates, PS model results can be considered acceptable as deleterious effects are really critical in thin plates. On the other hand, results provided by the SG model converge well for both thin and thick plates. The effectiveness of the SG model is established via high-accuracy results obtained in several examples. It is concluded that corrected SGN-FEM models are efficient alternatives for free-vibration and buckling analysis of Mindlin plate problems, and that precise elimination of parasitic shear is a requirement for sound analyses.