• Journal of the Korea Academia-Industrial cooperation Society
• /
• v.20 no.8
• /
• pp.8-14
• /
• 2019

Bird-strike is one of the most important design factors for safety in the aviation industry. Bird-strikes have been the cause of significant damage to aircraft and rotorcraft structures and the loss of life. This study used DYTRAN software to simulate the transient response of an Euler-Lagrangian composite helicopter blade that has been impacted by a bird. The Arbitrary Lagrangian Eulerian (ALE) method and a suitable equation of state were applied to model the bird. ALE was applied to the bird-strike analysis due to the large difference between the properties of the blade and bird. The debris of the bird was assumed to be a fluid and applied as Euler elements after the collision. Through the analysis of bird impacts, the leading-edge of the rotor blade (50.8 mm) was used to identify a positive margin of 1.18 based on the TSAI-FILL criteria. The results are assessed to be sufficiently reliable and may be evaluated to replace tests with various analysis conditions. The structural stability of the rotor blade could be assessed by applying various load conditions and different modeling methods in the future.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.1B
• /
• pp.95-109
• /
• 2008

We numerically analyzed the nonlinear shoaling, a plunging breaker and its accompanying energetic suspension of sediment at a bed, and a redistribution of suspended sediments by a down rush of preceding waves and the following plunger using SPH with a Gaussian kernel function, Lagrangian Dynamic Smagorinsky model (LDS), Van Rijn's pick up function. In that process, we came to the conclusion that the conventional model for the tractive force at a bottom like a quadratic law can not accurately describe the rapidly accelerating flow over a swash zone, and propose new methodology to accurately estimate the bottom tractive force. Using newly proposed wave model in this study, we can successfully duplicate severely deformed water surface profile, free falling water particles, a queuing splash after the landing of water particles on the free surface and a wave finger due to the structured vortex on a rear side of wave crest (Narayanaswamy and Dalrymple, 2002), a circulation of suspended sediments over a swash zone, net transfer of sediments clouds suspended over a swash zone toward the offshore, which so far have been regarded very difficult features to mimic in the computational fluid mechanics.

• Atmosphere
• /
• v.23 no.1
• /
• pp.1-11
• /
• 2013

The angular momentum transport of an idealized axisymmetric vortex in the developing stage was investigated using the Weather Research and Forecast (WRF) model. The balanced axisymmetric vortex was constructed based on an empirical function for tangential wind, and the temperature, geopotential, and surface pressure were obtained from the balanced equation. The numerical simulation was carried out for 6 days on the f-plane with the Sea Surface Temperature (SST) set as constant. The weak vortex at initial time was intensified with time, and reached the strength of tropical cyclone in a couple of days. The Absolute Angular Momentum (AAM) was transported along with the secondary circulation of the vortex. Total AAM integrated over a cylinder of radius of 2000 km decreased with simulation time, but total kinetic energy increased rapidly. From the budget analysis, it was found that the surface friction is mainly responsible for the decrease of total AAM. Also, contribution of the surface friction to the AAM loss was about 90% while that of horizontal advection was as small as 8%. The trajectory of neutral numerical tracers following the secondary circulation was presented for the Lagrangian viewpoint of the transports of absolute angular momentum. From the analysis using the trajectory of tracers it was found that the air parcel was under the influence of the surface friction continuously until it leaves the boundary layer near the core. Then the air parcel with reduced amount of angular momentum compared to its original amount was transported from boundary layer to upper level of the vortex and contributed to form the anti-cyclone. These results suggest that the tropical cyclone loses angular momentum as it develops, which is due to the dissipation of angular momentum by the surface friction.

• Proceedings of KOSOMES biannual meeting
• /
• 2009.06a
• /
• pp.89-94
• /
• 2009

This study aims to examine the the characteristics if the sea water circulation in related to the change of the water exchange rate according to the the development phases of Busan Port. To clarify the characteristics, the water exchange caused by construction if new-port and river discharge conditions was examined by the numerical experiments using the Lagrangian particle tracking model based on 2-D shallow water equation. This study deals with the charge of flow field and water exchange after redevelopment using numerical simulation technique, based on the field observation and 4 rivers discharge flow. As a result if numerical simulation, Variation if current velocity in Busan North Port is almost nothing before and after construction. Water exchange of channel in area of redevelopment is a little. Futhermore, the tidal exchange had a tendency to be small both before and after redevelopment by 60%.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.38 no.5
• /
• pp.74-84
• /
• 1996

Dynamic loading of structures often causes excursions of stresses well into the inelastic range and the influence of geometric changes on the dynamic response is also significant in many cases. Therefore, both material and geometric nonlinearity effects should be considered in case that a dynamic load acts on the structure. For developing a program to analyze the dynamic response of an axisymmetric shell in this study, the material nonlinearity effect on the dynamic response was formulated by the elasto-viscoplastic model highly corresponding to the real behavior of the material. Also, the geometrically nonlinear behavior is taken into account using a total Lagrangian coordinate system, and the equilibrium equation of motion was numerically solved by a central difference scheme. A complete finite element program has been developed and the results obtained by it are compared with those in the references 1 and 2. The results are in good agreement with each other. As a case study of its application, the developed program was applied to a dynamic response analysis of a nuclear reinforced concrete containment structure. The results obtained from the' numerical examples are summarized as follows : 1. The dynamic magnification factor of the displacement and the stress were unrelated with the concrete strength. 2. As shown by the results that the displacement dynamic magnification factor were form 1.7 to 2.3 and the stress dynamic magnification factor from 1.8 to 2.5, the dynamic magnification factor of stress were larger than that of displacement. 3. The dynamic magnification factor of stress on the exterior surface was larger than that on the interior surface of the structure.

• Journal of the Korean Crystal Growth and Crystal Technology
• /
• v.13 no.3
• /
• pp.145-151
• /
• 2003

Using computational fluid dynamics (CFD) method, we investigated three-dimensional fluid flow field and particle movement with respect to the injected gas flow rate variation in typical cyclone separator system. The results of numerical investigation were deduced by coupling the analysis of fluid flow field with Wavier-stokes equation and the tracking of the particle trajectory with Langrangian approach. It was shown that the increasing of injected gas flow rate resulted in the increasing of pressure loss in the separator. This change of inner pressure had an effect on an aspect of the fluid flow in the separator. Particle movement was determined by fluid flow in the separator and was fully depended on a diameter of particles under the fixed flow rate. Increasing of injected gas flow rate was led to an increasing of the trace of particle, so the particles moved to the lower part of the separator. For this reason, the minimum diameters of the particles were decreased and increased the separation rate under the fixed particle diameter. In conclusion, the changes of injected gas flow rate have an important factor to the fluctuation of the fluid flow field and particle trajectory in the separator.

• Journal of the Korea Convergence Society
• /
• v.9 no.6
• /
• pp.53-60
• /
• 2018

In this paper, we described the unknown parameter identification using Dual Extended Kalman Filter for precise control of 2-link manipulator. 2-link manipulator has highly non-linear characteristic with changed parameter thought tasks. The parameter kinds of mass and inertia of system is important to handle with the manipulator robustly. To solve the control problem by estimating the state and unknown parameters of the system through the proposed method. In order to verify the performance of proposed method, we simulate the implementation using Matlab and compare with results of RLS algorithm. At the results, proposed method has a better performance than those of RLS and verify the estimation performance in the parameter estimation.

• Steel and Composite Structures
• /
• v.35 no.1
• /
• pp.77-92
• /
• 2020

The present paper outlined a procedure for geometrically nonlinear dynamic analysis of functionally graded graphene platelets-reinforced (GPLR-FG) nanocomposite cylinder subjected to mechanical shock loading. The governing equation of motion for large deformation problems is derived using meshless local Petrov-Galerkin (MLPG) method based on total lagrangian approach. In the MLPG method, the radial point interpolation technique is employed to construct the shape functions. A micromechanical model based on the Halpin-Tsai model and rule of mixture is used for formulation the nonlinear functionally graded distribution of GPLs in polymer matrix of composites. Energy dissipation in analyses of the structure responding to dynamic loads is considered using the Rayleigh damping. The Newmark-Newton/Raphson method which is an incremental-iterative approach is implemented to solve the nonlinear dynamic equations. The results of the proposed method for homogenous material are compared with the finite element ones. A very good agreement is achieved between the MLPG and FEM with very fine meshing. In addition, the results have demonstrated that the MLPG method is more effective method compared with the FEM for very large deformation problems due to avoiding mesh distortion issues. Finally, the effect of GPLs distribution on strength, stiffness and dynamic characteristics of the cylinder are discussed in details. The obtained results show that the distribution of GPLs changed the mechanical properties, so a classification of different types and volume fraction exponent is established. Indeed by comparing the obtained results, the best compromise of nanocomposite cylinder is determined in terms of mechanical and dynamic properties for different load patterns. All these applications have shown that the present MLPG method is very effective for geometrically nonlinear analyses of GPLR-FG nanocomposite cylinder because of vanishing mesh distortion issue in large deformation problems. In addition, since in proposed method the distributed nodes are used for discretization the problem domain (rather than the meshing), modeling the functionally graded media yields to more accurate results.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.12
• /
• pp.1738-1746
• /
• 2002

In the FSI (Fluid-Structure Interaction) problems, two different governing equations are to be solved together. One is fur the fluid and the other for the structure. Furthermore, a kinematic constraint should be imposed along the boundary between the fluid and the structure. We use the combined formulation, which incorporates both the fluid and structure equations of motion into a single coupled variational equation so that it is not necessary to calculate the fluid force on the surface of structure explicitly when solving the equations of motion of the structure. A two-dimensional channel flow divided by a Bernoulli-Euler beam is considered and the dynamic response of the beam under the influence of channel flow is studied. The Navier-Stokes equations are solved using a P2P1 Galerkin finite element method with ALE (Arbitrary Lagrangian-Eulerian) algorithm. The internal structural damping effect is not considered in this study and numerical results are compared with a previous work fer steady case. In addition to the Reynolds number, two non-dimensional parameters, which govern this fluid-structure system, are proposed. It is found that the larger the dynamic viscosity and density of the fluid are, the larger the damping of the beam is. Also, the added mass is found to be linearly proportional to the density of the fluid.

• Journal of the Society of Naval Architects of Korea
• /
• v.51 no.5
• /
• pp.409-418
• /
• 2014

Simulation of flow past a complex marine structure requires a fine resolution in the vicinity of the structure, whereas a coarse resolution is enough far away from it. Therefore, a lot of grid cells may be wasted, when a simple Cartesian grid system is used for an Immersed Boundary Method (IBM). To alleviate this problems while maintaining the Cartesian frame work, we adopted an Adaptive Mesh Refinement (AMR) scheme where the grid system dynamically and locally refines as needed. In this study, We implemented a moving IBM and an AMR technique in our basic 3D incompressible Navier-Stokes solver. A Volume Of Fluid (VOF) method was used to effectively treat the free surface, and a recently developed Lagrangian Dynamic Subgrid-scale Model (LDSM) was incorporated in the code for accurate turbulence modeling. To capture vortex induced vibration accurately, the equation for the structure movement and the governing equations for fluid flow were solved at the same time implicitly. Also, We have developed an interface by using AutoLISP, which can properly distribute marker particles for IBM, compute the geometrical information of the object, and transfer it to the solver for the main simulation. To verify our numerical methodology, our results were compared with other authors' numerical and experimental results for the benchmark problems, revealing excellent agreement. Using the verified code, we investigated the following cases. (1) simulating flow around a floating sphere. (2) simulating flow past a marine structure.