• Coupled systems mechanics
• /
• v.8 no.1
• /
• pp.71-93
• /
• 2019

Finite element simulations of solid mechanics problems often involve the use of Non-Confirming Meshes (NCM) to increase accuracy in capturing nonlinear behavior, including damage and plasticity, in part of a solid domain without an undue increase in computational costs. In the presence of material nonlinearity and plasticity, higher-order variables are often needed to capture nonlinear behavior and material history on non-conforming interfaces. The most popular formulations for coupling non-conforming meshes are dual methods that involve the interpolation of a traction field on the interface. These methods are subject to the Ladyzhenskaya-Babuska-Brezzi (LBB) stability condition, and are therefore limited in their implementation with the higher-order elements needed to capture nonlinear material behavior. Alternatively, the enriched discontinuous Galerkin approach (EDGA) (Haikal and Hjelmstad 2010) is a primal method that provides higher order kinematic fields on the interface, and in which interface tractions are computed from local finite element estimates, therefore facilitating its implementation with nonlinear material models. The inclusion of higher-order interface variables, however, presents the issue of preserving material history at integration points when a increase in integration order is needed. In this study, the enriched discontinuous Galerkin approach (EDGA) is extended to the case of small-deformation plasticity. An interface-driven Gauss-Kronrod integration rule is proposed to enable adaptive enrichment on the interface while preserving history-dependent material data at existing integration points. The method is implemented using classical J2 plasticity theory as well as the pressure-dependent Drucker-Prager material model. We show that an efficient treatment of interface variables can improve algorithmic performance and provide a consistent approach for coupling non-conforming meshes in inelasticity.

• Nuclear Engineering and Technology
• /
• v.49 no.7
• /
• pp.1369-1378
• /
• 2017

The present paper reports the development of a computational code based on the Adaptive Group of Ink Drop Spread (AGIDS) for reconstruction of the neutron noise sources in reactor cores. AGIDS algorithm was developed as a fuzzy inference system based on the active learning method. The main idea of the active learning method is to break a multiple input-single output system into a single input-single output system. This leads to the ability to simulate a large system with high accuracy. In the present study, vibrating absorber-type neutron noise source in an International Atomic Energy Agency-two dimensional reactor core is considered in neutron noise calculation. The neutron noise distribution in the detectors was calculated using the Galerkin finite element method. Linear approximation of the shape function in each triangle element was used in the Galerkin finite element method. Both the real and imaginary parts of the calculated neutron distribution of the detectors were considered input data in the developed computational code based on AGIDS. The output of the computational code is the strength, frequency, and position (X and Y coordinates) of the neutron noise sources. The calculated fraction of variance unexplained error for output parameters including strength, frequency, and X and Y coordinates of the considered neutron noise sources were $0.002682{\sharp}/cm^3s$, 0.002682 Hz, and 0.004254 cm and 0.006140 cm, respectively.

• Structural Engineering and Mechanics
• /
• v.14 no.1
• /
• pp.71-84
• /
• 2002

A finite element formulation for the time-domain analysis of linear transient elastodynamic problems is presented based on the weak form obtained by applying the Galerkin's method to the integro-differential equations which contain the initial conditions implicitly and does not include the inertia terms. The weak form is extended temporally under the assumptions of the constant and linear time variations of field variables, since the time-stepping algorithms such as the Newmark method and the Wilson ${\theta}$-method are not necessary, obtaining two kinds of implicit finite element equations which are tested for numerical accuracy and convergency. Three classical examples having finite and infinite domains are solved and numerical results are compared with the other analytical and numerical solutions to show the versatility and accuracy of the presented formulation.

• Journal of Mechanical Science and Technology
• /
• v.20 no.1
• /
• pp.110-124
• /
• 2006

The multi scale wavelet-Galerkin method implemented in an adaptive manner has an advantage of obtaining accurate solutions with a substantially reduced number of interpolation points. The method is becoming popular, but its numerical efficiency still needs improvement. The objectives of this investigation are to present a new numerical scheme to improve the performance of the multi scale adaptive wavelet-Galerkin method and to give detailed implementation procedure. Specifically, the subdomain technique suitable for multiscale methods is developed and implemented. When the standard wavelet-Galerkin method is implemented without domain subdivision, the interaction between very long scale wavelets and very short scale wavelets leads to a poorly-sparse system matrix, which considerably worsens numerical efficiency for large-sized problems. The performance of the developed strategy is checked in terms of numerical costs such as the CPU time and memory size. Since the detailed implementation procedure including preprocessing and stiffness matrix construction is given, researchers having experiences in standard finite element implementation may be able to extend the multi scale method further or utilize some features of the multiscale method in their own applications.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1999.03b
• /
• pp.253-256
• /
• 1999

The Streamline Upwind Petrov-Galerkin(SUPG) finite element method is used to solve the two-dimensional laminar and turbulent flow. The flow is simulated by averaged Navier-Stokes equations with a penalty function approach and the lograithmic(k-$\varepsilon$) turbulent model is employed to take into account its turbulent behavior. The near-wall viscous sub-layer model is employed to approach the dominant viscous effects in the near wall zones. To find a good-enough initial guess of the Newton-Raphson iteration solving Nonlinear Matrix the Incremental method is considered for momentum and the Incomplete logarithmic turbu-lent equations for Turbulence. The validation of our method is investigated in comparision with published experimental data.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.10a
• /
• pp.375-382
• /
• 2003

Porous media consist of physically and chemically different materials and have an extremely complicated behavior due to the different material properties of each of its constituents. In addition, the internal structure of porous media has generally a complex geometry that makes the description of its mechanical behavior quite complex. Thus, in order to describe and clarify the deformation behavior of porous media, constitutive models for deformation of porous media coupling several effects such as flow of fluids of thermodynamical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian methods, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of solids and fluids. First of all, governing equations for saturated porous media based on ALE description are derived. Then, weak forms of these equations are obtained in order to implement numerical method using finite element method. Finally, Petrov-Galerkin method Is applied to develop finite element formulation.

• Journal of Korea Water Resources Association
• /
• v.32 no.3
• /
• pp.237-246
• /
• 1999

A finite element model is studied to simulate unsteady free surface flow based on dynamic wave equation and collocation method. The collocation method is used in conjunction with Hermite polynomials, and resulting matrix equations are solved by skyline method. The model is verified by applying to hydraulic jump, nonlinear disturbance propagation and dam-break flow in a horizontal frictionless channel. The computed results are compared with those by Bubnov-Galerkin and Petrov-Galerkin methods. It is also applied to the North Han River to simulate the floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge hydrographs. The suggested model has proven to be one of the promising scheme for simulating the gradually and rapidly varied unsteady flow in open channels.

• Journal of Wetlands Research
• /
• v.13 no.3
• /
• pp.395-409
• /
• 2011

The objective of this study is to develop an accurate and robust two-dimensional finite element method for turbulence simulation in open channels. The model is based on Streamline Upwind/Petrov-Galerkin finite element method and Boussinesq's eddy viscosity theory. The method developed in the study is depth-averaged mixing length model which assumes anisotropic and local equilibrium state of turbulence. The model calibration and validation were performed by comparing with analytical solutions and observed data. Several numerical simulations were carried out, which examined the performance of the turbulence model for the purpose of sensitivity analysis. The uniform channels that appear horizontal flow and vertical flow were carried out. The model was also applied to the Han river was in for the applicability test. The results were compared with the observed data. The suggested model displayed reasonable flow distribution compare to the observed data in natural river flow. As a result of this study, the two-dimensional finite element model provides a reliable results for flow distribution based on the turbulence simulation in open channels.

• Proceedings of the KSME Conference
• /
• 2001.11b
• /
• pp.475-480
• /
• 2001

The flowfield of transverse jet in a supersonic air stream subjected to shock wave turbulent boundary layer interactions is simulated numerically by Generalized Taylor Galerkin(GTG) finite element methods. Effects of turbulence are taken into account with a two-equation $(k-\varepsilon)$ model with a compressibility correction. Injection pressures and slot widths are varied in the present study. Pressure, separation extents, and penetration heights are compared with experimental data. Favorable comparisons with experimental measurements are demonstrated.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.26 no.3
• /
• pp.375-383
• /
• 2002

The flowfield of transverse jet in a supersonic air stream subjected to shock wave turbulent boundary layer interactions is simulated numerically by Generalized Taylor Galerkin(GTG) finite element methods. Effects of turbulence are taken into account with a two-equation (k-$\varepsilon$) model with a compressibility correction. Injection pressures and slot widths are varied in the present study. Pressure, separation extents, and penetration heights are compared with experimental data. Favorable comparisons with experimental measurements are demonstrated.