• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2011.04a
• /
• pp.443-444
• /
• 2011

This paper discusses a model order reduction for large order rotor dynamics systems results from the finite element discretization. Typical rotor systems consist of a rotor, built-on parts, and a support system, and require prudent consideration in their dynamic analysis models because they include unsymmetric stiffness, localized nonproportional damping and frequency dependent gyroscopic effects. When the finite element model has a very large number of degrees of freedom because of complex geometry, repeated dynamic analyses to investigate the critical speeds, stability, and unbalanced response are computationally very expensive to finish within a practical design cycle. In this paper, the Krylov-based model order reduction via moment matching significantly speeds up the dynamic analyses necessary to check eigenvalues and critical speeds of a Nelson-Vaugh rotor system. With this approach the dynamic simulation is efficiently repeated via a reduced system by changing a running rotational speed because it can be preserved as a parameter in the process of model reduction. The Campbell diagram by the reduced system shows very good agreement with that of the original system. A 3-D finite element model of the Nelson-Vaugh rotor system is taken as a numerical example to demonstrate the advantages of this model reduction for rotor dynamic simulation.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.2 no.5
• /
• pp.58-65
• /
• 1994

In order to reduce particulate emissions from diesel vehicles, mathematical model is established and analyzed on ceramic wall-flow monolith filter. A wall-flow monolith filter placed in the exhaust stream of a diesel engine can effectively limit the emission of diesel particulates through the monolith. The accumulated particulates can then be periodically combusted inside the monolith by directing hot gas to the monolith while normal engine exhaust is routed around the monolith system. The resulting low flow rates through the monolith require consideration of gas dynamics through the channels as well as particulate combustion to analyze this regeneration process. A mathematical model of the regeneration is formulated as a system of nonlinear partial differential equations describing the conservation of mass, momentum and energy. Numerical solutions are obtained by using a finite difference techniques for the spatial discretization. So we can use filter simulation program for the purpose of filter design and actual filter regeneration

• Structural Engineering and Mechanics
• /
• v.50 no.1
• /
• pp.19-36
• /
• 2014

This paper focuses on a model order reduction (MOR) for large-scale rotordynamic systems by using finite element discretization. Typical rotor-bearing systems consist of a rotor, built-on parts, and a support system. These systems require careful consideration in their dynamic analysis modeling because they include unsymmetrical stiffness, localized nonproportional damping, and frequency-dependent gyroscopic effects. Because of this complex geometry, the finite element model under consideration may have a very large number of degrees of freedom. Thus, the repeated dynamic analyses used to investigate the critical speeds, stability, and unbalanced response are computationally very expensive to complete within a practical design cycle. In this study, we demonstrate that a Krylov subspace-based MOR via moment matching significantly speeds up the rotordynamic analyses needed to check the whirling frequencies and critical speeds of large rotor systems. This approach is very efficient, because it is possible to repeat the dynamic simulation with the help of a reduced system by changing the operating rotational speed, which can be preserved as a parameter in the process of model reduction. Two examples of rotordynamic systems show that the suggested MOR provides a significant reduction in computational cost for a Campbell diagram analysis, while maintaining accuracy comparable to that of the original systems.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2000.11a
• /
• pp.228-231
• /
• 2000

The active vibration control of laminated composite shell has been performed with the optimized sensor/actuator system. PVDF film is used fur the material of sensor/actuator. Finite element method is utilized to model the whole structure including the piezoelectric sensor/actuator system, The distributed selective modal sensor/actuator system is established to prevent the adverse effect of spillover. In the finite element discretization process, the nine-node shell element with five nodal degrees of freedoms is used. Electrode patterns and lamination angles of sensor/actuator are optimized using genetic algorithm. Sensor is designed to minimize the observation spillover, and actuator is designed to minimize the system energy of the control modes under a given initial condition. Modal sensor/actuator profiles are optimized for the first and the second modes suppression of singly curved cantilevered composite shell structure. Discrete LQG method is used as a control law. The real time vibration control with profile optimized sensor/actuator system has been performed. Experimental result shows successful performance of the integrated structure for the active vibration control.

• Advances in aircraft and spacecraft science
• /
• v.3 no.2
• /
• pp.149-170
• /
• 2016

Multi-Disciplinary Optimization (MDO) is widely used to handle the advanced design in several engineering applications. Such applications are commonly simulation-based, in order to capture the physics of the phenomena under study. This framework demands fast optimization algorithms as well as trustworthy numerical analyses, and a synergic integration between the two is required to obtain an efficient design process. In order to meet these needs, an adaptive Computational Fluid Dynamics (CFD) solver and a fast optimization algorithm have been developed and combined by the authors. The CFD solver is based on a high-order discontinuous Galerkin discretization while the optimization algorithm is a high-performance version of the Artificial Bee Colony method. In this work, they are used to address a typical aero-mechanical problem encountered in turbomachinery design. Interesting achievements in the considered test case are illustrated, highlighting the potential applicability of the proposed approach to other engineering problems.

• 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.

• 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.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.147-161
• /
• 2023

Based on the ideas of variational differential quadrature (VDQ) and finite element method (FEM), a numerical approach named as VDQFEM is applied herein to study the large deformations of plate-type structures under static loading with arbitrary shape hole made of functionally graded graphene platelet-reinforced composite (FG-GPLRC) in the context of higher-order shear deformation theory (HSDT). The material properties of composite are approximated based upon the modified Halpin-Tsai model and rule of mixture. Furthermore, various FG distribution patterns are considered along the thickness direction of plate for GPLs. Using novel vector/matrix relations, the governing equations are derived through a variational approach. The matricized formulation can be efficiently employed in the coding process of numerical methods. In VDQFEM, the space domain of structure is first transformed into a number of finite elements. Then, the VDQ discretization technique is implemented within each element. As the last step, the assemblage procedure is performed to derive the set of governing equations which is solved via the pseudo arc-length continuation algorithm. Also, since HSDT is used herein, the mixed formulation approach is proposed to accommodate the continuity of first-order derivatives on the common boundaries of elements. Rectangular and circular plates under various boundary conditions with circular/rectangular/elliptical cutout are selected to generate the numerical results. In the numerical examples, the effects of geometrical properties and reinforcement with GPL on the nonlinear maximum deflection-transverse load amplitude curve are studied.

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

The aim of this work is a numerical comparison (FEM) between lattice pyramidal-core panel and honeycomb core panel for different core thicknesses. By evaluating the mid-span deflection, the shear rigidity and the shear modulus for both core types and different core thicknesses, it is possible to define which core type has got the best mechanical behaviour for each thickness and the evolution of that behaviour as far as the thickness increases. Since a specific base geometry has been used for the lattice pyramidal core, the comparison gives us the opportunity to investigate the unit cell strut angle giving the higher mechanical properties. The presented work considers a detailed FEM modelling of a standard 3-point bending test (ASTM C393/C393M Standard Practice). Detailed FEM modelling addresses to detailed discretization of cores by means of beam elements for lattice core and shell elements for honeycomb core. Facings, instead, have been modelled by using shell elements for both sandwich panels. On lattice core structure, elements of core and facings are directly connected, to better simulate the additive manufacturing process. Otherwise, an MPC-based constraint between facings and core has been used for honeycomb core structure. Both sandwich panels are entirely built of Aluminium alloy. Prior to compare the two models, the FEM sandwich panel model with lattice pyramidal core needs to be validated with 3-point bending test experimental results, in order to ensure a good reliability of the FEM approach and of the comparison. Furthermore, the analytical validation has been performed according to Allen's theory. The FEM analysis is linear static with an increasing midspan load ranging from 50N up to 500N.

• Nuclear Engineering and Technology
• /
• v.55 no.9
• /
• pp.3367-3382
• /
• 2023

Smoothed Particle Hydrodynamics (SPH) is a Lagrangian computational fluid dynamics method that has been widely used in the analysis of physical phenomena characterized by large deformation or multi-phase flow analysis, including free surface. Despite the recent implementation of eddy-viscosity models in SPH methodology, sophisticated turbulent analysis using Lagrangian methodology has been limited due to the lack of computational performance and numerical consistency. In this study, we implement the standard and dynamic Smagorinsky model and dynamic Vreman model as sub-particle scale models based on a weakly compressible SPH solver. The large eddy simulation method is numerically identical to the spatial discretization method of smoothed particle dynamics, enabling the intuitive implementation of the turbulence model. Furthermore, there is no additional filtering process required for physical variables since the sub-grid scale filtering is inherently processed in the kernel interpolation. We simulate lid-driven flow under transition and turbulent conditions as a benchmark. The simulation results show that the dynamic Vreman model produces consistent results with experimental and numerical research regarding Reynolds averaged physical quantities and flow structure. Spectral analysis also confirms that it is possible to analyze turbulent eddies with a smaller length scale using the dynamic Vreman model with the same particle size.