Galerkin finite element methods for wave problems
Indian Academy of Sciences (India)
these methods. Keywords. hp-Finite element method; continuous Galerkin methods; wave solutions; Gibbs' phenomenon. 1. Introduction. Galerkin methods belong to the class of solution methods for PDEs where the solution residue is minimized giving rise to the well-known weak formulation of problems. In this approach,.
Finite element or Galerkin type semidiscrete schemes
Durgun, K.
1983-01-01
A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.
Discontinuous Galerkin finite element methods for hyperbolic differential equations
van der Vegt, Jacobus J.W.; van der Ven, H.; Boelens, O.J.; Boelens, O.J.; Toro, E.F.
2002-01-01
In this paper a suryey is given of the important steps in the development of discontinuous Galerkin finite element methods for hyperbolic partial differential equations. Special attention is paid to the application of the discontinuous Galerkin method to the solution of the Euler equations of gas
Galerkin finite element methods for wave problems
Indian Academy of Sciences (India)
We compare here the accuracy, stability and wave propagation properties of a few Galerkin methods. The basic Galerkin methods with piecewise linear basis functions (called G1FEM here) and quadratic basis functions (called G2FEM) have been compared with the streamwise-upwind Petrov Galerkin (SUPG) method for ...
An implicit discontinuous Galerkin finite element model for water waves
van der Vegt, Jacobus J.W.; Ambati, V.R.; Bokhove, Onno
2005-01-01
We discuss a new higher order accurate discontinuous Galerkin finite element method for non-linear free surface gravity waves. The algorithm is based on an arbitrary Lagrangian Eulerian description of the flow field using deforming elements and a moving mesh, which makes it possible to represent
An implicit discontinuous Galerkin finite element model for water waves
van der Vegt, Jacobus J.W.; Tomar, S.K.; Yao, Z.H.; Yuan, M.W.; Zhong, W.X.
2004-01-01
An overview is given of a discontinuous Galerkin finite element method for linear free surface water waves. The method uses an implicit time integration method which is unconditionally stable and does not suffer from the frequently encountered mesh dependent saw-tooth type instability at the free
Space-time discontinuous Galerkin finite element methods
van der Vegt, Jacobus J.W.; Deconinck, H.; Ricchiuto, M.
2006-01-01
In these notes an introduction is given to space-time discontinuous Galerkin (DG) finite element methods for hyperbolic and parabolic conservation laws on time dependent domains. the space-time DG discretization is explained in detail, including the definition of the numerical fluxes and
Galerkin finite element methods for wave problems
Indian Academy of Sciences (India)
piecewise linear and quadratic basis functions for solving the one-dimensional wave equation. 2. One-dimensional wave equation. First, we briefly discuss the Galerkin method that employs piecewise quadratic polynomials for the basis or interpolating functions. We will call this as G2FEM for ease of reference. Here, one ...
Massey, Thomas Christopher
2002-01-01
A Flexible Galerkin Finite Element Method (FGM) is a hybrid class of finite element methods that combine the usual continuous Galerkin method with the now popular discontinuous Galerkin method (DGM). A detailed description of the formulation of the FGM on a hyperbolic partial differential equation, as well as the data structures used in the FGM algorithm is presented. Some hp-convergence results and computational cost are included. Additionally, an a posteriori error estimate f...
Finite element and discontinuous Galerkin methods for transient wave equations
Cohen, Gary
2017-01-01
This monograph presents numerical methods for solving transient wave equations (i.e. in time domain). More precisely, it provides an overview of continuous and discontinuous finite element methods for these equations, including their implementation in physical models, an extensive description of 2D and 3D elements with different shapes, such as prisms or pyramids, an analysis of the accuracy of the methods and the study of the Maxwell’s system and the important problem of its spurious free approximations. After recalling the classical models, i.e. acoustics, linear elastodynamics and electromagnetism and their variational formulations, the authors present a wide variety of finite elements of different shapes useful for the numerical resolution of wave equations. Then, they focus on the construction of efficient continuous and discontinuous Galerkin methods and study their accuracy by plane wave techniques and a priori error estimates. A chapter is devoted to the Maxwell’s system and the important problem ...
Discontinuous Galerkin finite element method for shallow two-phase flows
Rhebergen, Sander; Bokhove, Onno; van der Vegt, Jacobus J.W.
We present a discontinuous Galerkin finite element method for two depth-averaged two-phase flow models. One of these models contains nonconservative products for which we developed a discontinuous Galerkin finite element formulation in Rhebergen et al. (2008) J. Comput. Phys. 227, 1887-1922. The
Finite element, discontinuous Galerkin, and finite difference evolution schemes in spacetime
International Nuclear Information System (INIS)
Zumbusch, G
2009-01-01
Numerical schemes for Einstein's vacuum equation are developed. Einstein's equation in harmonic gauge is second-order symmetric hyperbolic. It is discretized in four-dimensional spacetime by finite differences, finite elements and interior penalty discontinuous Galerkin methods, the latter being related to Regge calculus. The schemes are split into space and time and new time-stepping schemes for wave equations are derived. The methods are evaluated for linear and nonlinear test problems of the Apples-with-Apples collection.
International Nuclear Information System (INIS)
Guerreiro, J.N.C.; Loula, A.F.D.
1988-12-01
The mixed Petrov-Galerkin finite element formulation is applied to transiente and steady state creep problems. Numerical analysis has shown additional stability of this method compared to classical Galerkin formulations. The accuracy of the new formulation is confirmed in some representative examples of two dimensional and axisymmetric problems. (author) [pt
van der Vegt, Jacobus J.W.; van der Ven, H.
1998-01-01
A new discretization method for the three-dimensional Euler equations of gas dynamics is presented, which is based on the discontinuous Galerkin finite element method. Special attention is paid to an efficient implementation of the discontinuous Galerkin method that minimizes the number of flux
Dual and primal mixed Petrov-Galerkin finite element methods in heat transfer problems
International Nuclear Information System (INIS)
Loula, A.F.D.; Toledo, E.M.
1988-12-01
New mixed finite element formulations for the steady state heat transfer problem are presented with no limitation in the choice of conforming finite element spaces. Adding least square residual forms of the governing equations of the classical Galerkin formulation the original saddle point problem is transformed into a minimization problem. Stability analysis, error estimates and numerical results are presented, confirming the error estimates and the good performance of this new formulation. (author) [pt
Discontinuous Galerkin finite element methods for (non)conservative partial differential equations
Rhebergen, Sander
2010-01-01
The first research topic in this thesis is the development of discontinuous Galerkin (DG) finite element methods for partial differential equations containing nonconservative products, which are present in many two-phase flow models. For this, we combine the theory of Dal Maso, LeFloch and Murat, in
Applications of Taylor-Galerkin finite element method to compressible internal flow problems
Sohn, Jeong L.; Kim, Yongmo; Chung, T. J.
1989-01-01
A two-step Taylor-Galerkin finite element method with Lapidus' artificial viscosity scheme is applied to several test cases for internal compressible inviscid flow problems. Investigations for the effect of supersonic/subsonic inlet and outlet boundary conditions on computational results are particularly emphasized.
Space-time discontinuous Galerkin finite element method for inviscid gas dynamics
van der Ven, H.; van der Vegt, Jacobus J.W.; Bouwman, E.G.; Bathe, K.J.
2003-01-01
In this paper an overview is given of the space-time discontinuous Galerkin finite element method for the solution of the Euler equations of gas dynamics. This technique is well suited for problems which require moving meshes to deal with changes in the domain boundary. The method is demonstrated
Two fluid space-time discontinuous Galerkin finite element method. Part II: Applications
Sollie, W.E.H.; van der Vegt, Jacobus J.W.
2009-01-01
The numerical method for two fluid flow computations presented in Sollie, Bokhove \\& van der Vegt, Two Fluid Space-Time Discontinuous Galerkin Finite Element Method. Part I: Numerical Algorithm is applied to a number of one and two dimensional single and two fluid test problems, including a magma -
The Galerkin Finite Element Method for A Multi-term Time-Fractional Diffusion equation
Jin, Bangti; Lazarov, Raytcho; Liu, Yikan; Zhou, Zhi
2014-01-01
We consider the initial/boundary value problem for a diffusion equation involving multiple time-fractional derivatives on a bounded convex polyhedral domain. We analyze a space semidiscrete scheme based on the standard Galerkin finite element method using continuous piecewise linear functions. Nearly optimal error estimates for both cases of initial data and inhomogeneous term are derived, which cover both smooth and nonsmooth data. Further we develop a fully discrete scheme based on a finite...
Directory of Open Access Journals (Sweden)
Haitao Che
2011-01-01
Full Text Available We investigate a H1-Galerkin mixed finite element method for nonlinear viscoelasticity equations based on H1-Galerkin method and expanded mixed element method. The existence and uniqueness of solutions to the numerical scheme are proved. A priori error estimation is derived for the unknown function, the gradient function, and the flux.
International Nuclear Information System (INIS)
Hughes, T.J.R.; Hulbert, G.M.; Franca, L.P.
1988-10-01
Galerkin/least-squares finite element methods are presented for advective-diffusive equations. Galerkin/least-squares represents a conceptual simplification of SUPG, and is in fact applicable to a wide variety of other problem types. A convergence analysis and error estimates are presented. (author) [pt
Liu, Meilin
2011-07-01
A discontinuous Galerkin finite element method (DG-FEM) with a highly-accurate time integration scheme is presented. The scheme achieves its high accuracy using numerically constructed predictor-corrector integration coefficients. Numerical results show that this new time integration scheme uses considerably larger time steps than the fourth-order Runge-Kutta method when combined with a DG-FEM using higher-order spatial discretization/basis functions for high accuracy. © 2011 IEEE.
The use of Galerkin finite-element methods to solve mass-transport equations
Grove, David B.
1977-01-01
The partial differential equation that describes the transport and reaction of chemical solutes in porous media was solved using the Galerkin finite-element technique. These finite elements were superimposed over finite-difference cells used to solve the flow equation. Both convection and flow due to hydraulic dispersion were considered. Linear and Hermite cubic approximations (basis functions) provided satisfactory results: however, the linear functions were computationally more efficient for two-dimensional problems. Successive over relaxation (SOR) and iteration techniques using Tchebyschef polynomials were used to solve the sparce matrices generated using the linear and Hermite cubic functions, respectively. Comparisons of the finite-element methods to the finite-difference methods, and to analytical results, indicated that a high degree of accuracy may be obtained using the method outlined. The technique was applied to a field problem involving an aquifer contaminated with chloride, tritium, and strontium-90. (Woodard-USGS)
Wheeler, Mary
2013-11-16
We study the numerical approximation on irregular domains with general grids of the system of poroelasticity, which describes fluid flow in deformable porous media. The flow equation is discretized by a multipoint flux mixed finite element method and the displacements are approximated by a continuous Galerkin finite element method. First-order convergence in space and time is established in appropriate norms for the pressure, velocity, and displacement. Numerical results are presented that illustrate the behavior of the method. © Springer Science+Business Media Dordrecht 2013.
The Galerkin finite element method for a multi-term time-fractional diffusion equation
Jin, Bangti
2015-01-01
© 2014 The Authors. We consider the initial/boundary value problem for a diffusion equation involving multiple time-fractional derivatives on a bounded convex polyhedral domain. We analyze a space semidiscrete scheme based on the standard Galerkin finite element method using continuous piecewise linear functions. Nearly optimal error estimates for both cases of initial data and inhomogeneous term are derived, which cover both smooth and nonsmooth data. Further we develop a fully discrete scheme based on a finite difference discretization of the time-fractional derivatives, and discuss its stability and error estimate. Extensive numerical experiments for one- and two-dimensional problems confirm the theoretical convergence rates.
A High Order Filter with Galerkin Finite Element Method for the Spherical Local domain Model
Lee, C. H.; Cheong, H. B.; Kang, H. G.
2017-12-01
A High Order Filter with Galerkin Finite Element Method for the Spherical Local domain ModelChung-Hui Lee1 and Hyeong-Bin Cheong and Hyun-Gyu KangDepartment of Environmental Atmospheric Sciences, Pukyong National University, Busan, Korea (1 chlee@pukyong.ac.kr) A high-order filter with Galerkin finite element method is constructed by applying a two dimensional finite element method with quadrilateral basis functions to the spherical limited area domain. The quadrilateral basis function is defined as four shape-functions on separate four grid-boxes which share the same gridpoint. A first-order derivative is represented with an algebraic equation consisting of nine point stencil. Helmholtz equation on a sphere is the basic component of the high order filter and the filtering is performed by solving this equation with two dimensional finite element method. As the theory describes, for spherical Laplacian operator and first-order derivative, the convergence rates of the error were revealed to be second-order and fourth-order, respectively. In addition, since the convergence rate of errors for the filter in this study was the same as the filter with Fourier finite element method, the accuracy of the filter is comparable to the filter based on the Fourier finite element method. The high-order filter was applied to the WRF (Weather Research and Forecasting) as hyper-viscosity and its performance was compared with those of the built-in viscosity scheme of the WRF model. As a result of the tropical cyclone simulation, the forecast error for the high-order filter and the built-in viscosity were similar for the minimum pressure and track prediction. However, for the precipitation and rainfall distribution, the prediction with high-order filter appeared closer to observations than those with built-in viscosity.
Numerical modelling of MPA-CVD reactors with the discontinuous Galerkin finite element method
Houston, Paul; Sime, Nathan
2017-07-01
In this article we develop a fully self consistent mathematical model describing the formation of a hydrogen plasma in a microwave power assisted chemical vapour deposition (MPA-CVD) reactor employed for the manufacture of synthetic diamond. The underlying multi-physics model includes constituent equations for the background gas mass average velocity, gas temperature, electromagnetic field energy and plasma density. The proposed mathematical model is numerically approximated based on exploiting the discontinuous Galerkin finite element method. We demonstrate the practical performance of this computational approach on a variety of CVD reactor geometries for a range of operating conditions.
Liu, Meilin
2012-08-01
A discontinuous Galerkin finite element method (DG-FEM) with a highly accurate time integration scheme for solving Maxwell equations is presented. The new time integration scheme is in the form of traditional predictor-corrector algorithms, PE CE m, but it uses coefficients that are obtained using a numerical scheme with fully controllable accuracy. Numerical results demonstrate that the proposed DG-FEM uses larger time steps than DG-FEM with classical PE CE) m schemes when high accuracy, which could be obtained using high-order spatial discretization, is required. © 1963-2012 IEEE.
Huang, Chih-Hsu; Lin, Chou-Ching K; Ju, Ming-Shaung
2015-02-01
Compared with the Monte Carlo method, the population density method is efficient for modeling collective dynamics of neuronal populations in human brain. In this method, a population density function describes the probabilistic distribution of states of all neurons in the population and it is governed by a hyperbolic partial differential equation. In the past, the problem was mainly solved by using the finite difference method. In a previous study, a continuous Galerkin finite element method was found better than the finite difference method for solving the hyperbolic partial differential equation; however, the population density function often has discontinuity and both methods suffer from a numerical stability problem. The goal of this study is to improve the numerical stability of the solution using discontinuous Galerkin finite element method. To test the performance of the new approach, interaction of a population of cortical pyramidal neurons and a population of thalamic neurons was simulated. The numerical results showed good agreement between results of discontinuous Galerkin finite element and Monte Carlo methods. The convergence and accuracy of the solutions are excellent. The numerical stability problem could be resolved using the discontinuous Galerkin finite element method which has total-variation-diminishing property. The efficient approach will be employed to simulate the electroencephalogram or dynamics of thalamocortical network which involves three populations, namely, thalamic reticular neurons, thalamocortical neurons and cortical pyramidal neurons. Copyright © 2014 Elsevier Ltd. All rights reserved.
Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman’s Equations
Iliev, Oleg P.
2010-01-01
We present a two-scale finite element method for solving Brinkman\\'s equations with piece-wise constant coefficients. This system of equations model fluid flows in highly porous, heterogeneous media with complex topology of the heterogeneities. We make use of the recently proposed discontinuous Galerkin FEM for Stokes equations by Wang and Ye in [12] and the concept of subgrid approximation developed for Darcy\\'s equations by Arbogast in [4]. In order to reduce the error along the coarse-grid interfaces we have added a alternating Schwarz iteration using patches around the coarse-grid boundaries. We have implemented the subgrid method using Deal.II FEM library, [7], and we present the computational results for a number of model problems. © 2010 Springer-Verlag Berlin Heidelberg.
Johnson, Ryan; Kercher, Andrew; Schwer, Douglas; Corrigan, Andrew; Kailasanath, Kazhikathra
2017-11-01
This presentation focuses on the development of a Discontinuous Galerkin (DG) method for application to chemically reacting flows. The in-house code, called Propel, was developed by the Laboratory of Computational Physics and Fluid Dynamics at the Naval Research Laboratory. It was designed specifically for developing advanced multi-dimensional algorithms to run efficiently on new and innovative architectures such as GPUs. For these results, Propel solves for convection and diffusion simultaneously with detailed transport and thermodynamics. Chemistry is currently solved in a time-split approach using Strang-splitting with finite element DG time integration of chemical source terms. Results presented here show canonical unsteady reacting flow cases, such as co-flow and splitter plate, and we report performance for higher order DG on CPU and GPUs.
Uranus, H.P.; Hoekstra, Hugo; van Groesen, Embrecht W.C.
A Galerkin finite element scheme furnished with 1st-order Bayliss-Gunzburger-Turkel-like boundary conditions is formulated to compute both the guided and leaky modes of anisotropic channel waveguides of non-magnetic materials with diagonal permittivity tensor. The scheme is formulated using
Uranus, H.P.; Hoekstra, Hugo; van Groesen, Embrecht W.C.
2003-01-01
A Galerkin finite element scheme furnished with 1st-order Bayliss-Gunzberger-Turkel-like boundary conditions is formulated to compute both the guided and leaky modes of anisotropic channel waveguides of non-magnetic material with diagonal permitivity tensor. The scheme is formulated using
A nodal discontinuous Galerkin finite element method for nonlinear elastic wave propagation.
Bou Matar, Olivier; Guerder, Pierre-Yves; Li, YiFeng; Vandewoestyne, Bart; Van Den Abeele, Koen
2012-05-01
A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation.
Directory of Open Access Journals (Sweden)
Jinfeng Wang
2014-01-01
Full Text Available We discuss and analyze an H1-Galerkin mixed finite element (H1-GMFE method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H1-GMFE method. Based on the discussion on the theoretical error analysis in L2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H1-norm. Moreover, we derive and analyze the stability of H1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure.
Wang, Jinfeng; Zhao, Meng; Zhang, Min; Liu, Yang; Li, Hong
2014-01-01
We discuss and analyze an H 1-Galerkin mixed finite element (H 1-GMFE) method to look for the numerical solution of time fractional telegraph equation. We introduce an auxiliary variable to reduce the original equation into lower-order coupled equations and then formulate an H 1-GMFE scheme with two important variables. We discretize the Caputo time fractional derivatives using the finite difference methods and approximate the spatial direction by applying the H 1-GMFE method. Based on the discussion on the theoretical error analysis in L 2-norm for the scalar unknown and its gradient in one dimensional case, we obtain the optimal order of convergence in space-time direction. Further, we also derive the optimal error results for the scalar unknown in H 1-norm. Moreover, we derive and analyze the stability of H 1-GMFE scheme and give the results of a priori error estimates in two- or three-dimensional cases. In order to verify our theoretical analysis, we give some results of numerical calculation by using the Matlab procedure. PMID:25184148
Directory of Open Access Journals (Sweden)
L. Jones Tarcius Doss
2012-01-01
Full Text Available A quadrature-based mixed Petrov-Galerkin finite element method is applied to a fourth-order linear ordinary differential equation. After employing a splitting technique, a cubic spline trial space and a piecewise linear test space are considered in the method. The integrals are then replaced by the Gauss quadrature rule in the formulation itself. Optimal order a priori error estimates are obtained without any restriction on the mesh.
Kou, Jisheng
2013-06-20
We analyze a combined method consisting of the mixed finite element method for pressure equation and the discontinuous Galerkin method for saturation equation for the coupled system of incompressible two-phase flow in porous media. The existence and uniqueness of numerical solutions are established under proper conditions by using a constructive approach. Optimal error estimates in L2(H1) for saturation and in L∞(H(div)) for velocity are derived. Copyright © 2013 John Wiley & Sons, Ltd.
Abedi, Reza; Mudaliar, Saba
2017-12-01
We present an asynchronous spacetime discontinuous Galerkin (aSDG) method for time domain electromagnetics in which space and time are directly discretized. By using differential forms we express Maxwell's equations and consequently their discontinuous Galerkin discretization for arbitrary domains in spacetime. The elements are discretized with electric and magnetic basis functions that are discontinuous across all inter-element boundaries and can have arbitrary high and per element spacetime orders. When restricted to unstructured grids that satisfy a specific causality constraint, the method has a local and asynchronous solution procedure with linear solution complexity in terms of the number of elements. We numerically investigate the convergence properties of the method for 1D to 3D uniform grids for energy dissipation, an error relative to the exact solution, and von Neumann dissipation and dispersion errors. Two dimensional simulations demonstrate the effectiveness of the method in resolving sharp wave fronts.
Directory of Open Access Journals (Sweden)
Xinming Zhang
2009-01-01
Full Text Available A wavelet Galerkin finite-element method is proposed by combining the wavelet analysis with traditional finite-element method to analyze wave propagation phenomena in fluid-saturated porous medium. The scaling functions of Daubechies wavelets are considered as the interpolation basis functions to replace the polynomial functions, and then the wavelet element is constructed. In order to overcome the integral difficulty for lacking of the explicit expression for the Daubechies wavelets, a kind of characteristic function is introduced. The recursive expression of calculating the function values of Daubechies wavelets on the fraction nodes is deduced, and the rapid wavelet transform between the wavelet coefficient space and the wave field displacement space is constructed. The results of numerical simulation demonstrate that the method is effective.
Memon, Sajid
2012-01-01
In this article, a posteriori error estimates are derived for mixed finite element Galerkin approximations to second order linear parabolic initial and boundary value problems. Using mixed elliptic reconstructions, a posteriori error estimates in L∞(L2)- and L2(L2)-norms for the solution as well as its flux are proved for the semidiscrete scheme. Finally, based on a backward Euler method, a completely discrete scheme is analyzed and a posteriori error bounds are derived, which improves upon earlier results on a posteriori estimates of mixed finite element approximations to parabolic problems. Results of numerical experiments verifying the efficiency of the estimators have also been provided. © 2012 Society for Industrial and Applied Mathematics.
Adib, Arash; Poorveis, Davood; Mehraban, Farid
2018-03-01
In this research, two equations are considered as examples of hyperbolic and elliptic equations. In addition, two finite element methods are applied for solving of these equations. The purpose of this research is the selection of suitable method for solving each of two equations. Burgers' equation is a hyperbolic equation. This equation is a pure advection (without diffusion) equation. This equation is one-dimensional and unsteady. A sudden shock wave is introduced to the model. This wave moves without deformation. In addition, Laplace's equation is an elliptical equation. This equation is steady and two-dimensional. The solution of Laplace's equation in an earth dam is considered. By solution of Laplace's equation, head pressure and the value of seepage in the directions X and Y are calculated in different points of earth dam. At the end, water table is shown in the earth dam. For Burgers' equation, least-square method can show movement of wave with oscillation but Galerkin method can not show it correctly (the best method for solving of the Burgers' equation is discrete space by least-square finite element method and discrete time by forward difference.). For Laplace's equation, Galerkin and least square methods can show water table correctly in earth dam.
Uranus, H.P.; Hoekstra, Hugo; van Groesen, Embrecht W.C.
A simple high-order Galerkin finite element scheme is formulated to compute both the guided and leaky modes of anisotropic planar waveguides with a diagonal permitivity tensor. Transparent boundary conditions derived from the Sommerfeld radiation conditions are used to model the fields at the
Czech Academy of Sciences Publication Activity Database
Dolejší, V.; Kůs, Pavel
2008-01-01
Roč. 73, č. 12 (2008), s. 1739-1766 ISSN 0029-5981 Keywords : backward difference formula * discontinuous Galerkin method * adaptive choice of the time step Subject RIV: BA - General Mathematics Impact factor: 2.229, year: 2008 http://onlinelibrary.wiley.com/doi/10.1002/nme.2143/abstract
Directory of Open Access Journals (Sweden)
Alberto Zingaro
2018-02-01
Full Text Available Microfluidics is a multidisciplinary area founding applications in several fields such as the aerospace industry. Microelectromechanical systems (MEMS are mainly adopted for flow control, micropower generation and for life support and environmental control for space applications. Microflows are modeled relying on both a continuum and molecular approach. In this paper, the compressible Navier–Stokes (CNS equations have been adopted to solve a two-dimensional unsteady flow for a viscous micro shock-channel problem. In microflows context, as for the most gas dynamics applications, the CNS equations are usually discretized in space using finite volume method (FVM. In the present paper, the PDEs are discretized with the nodal discontinuous Galerkin finite element method (DG–FEM in order to understand how the method performs at microscale level for compressible flows. Validation is performed through a benchmark test problem for microscale applications. The error norms, order of accuracy and computational cost are investigated in a grid refinement study, showing a good agreement and increasing accuracy with reference data as the mesh is refined. The effects of different explicit Runge–Kutta schemes and of different time step sizes have also been studied. We found that the choice of the temporal scheme does not really affect the accuracy of the numerical results.
van der Vegt, Jacobus J.W.; Izsak, F.; Bokhove, Onno
2007-01-01
Abstract. A detailed a priori error estimate is provided for a continuous-discontinuous Galerkin ��?nite element method for the generalized 2D vorticity dynamics equations. These equations describe several types of geophysical ﬂows, including the Euler equations. The algorithm consists of a
International Nuclear Information System (INIS)
Pereira, L.C.M.
1991-05-01
A methodology for thermo-hydraulic analysis of pressurized water reactor cores using a discontinuous Petrov-Galerkin formulation and a frontal method for solution of linear systems with hard disk storage is presented. this procedure aims at achieving significant economy of computer memory. the complete development of the subchannel equations is presented, beginning in the conservation equations. The cast in variational and finite element formulations is also presented. (author). 30 refs, 14 tabs
Directory of Open Access Journals (Sweden)
Maria Carla Piastra
2018-02-01
Full Text Available In Electro- (EEG and Magnetoencephalography (MEG, one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017. It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages, be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM% of 1.5% and mean magnitude errors (MAG% of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented
Piastra, Maria Carla; Nüßing, Andreas; Vorwerk, Johannes; Bornfleth, Harald; Oostenveld, Robert; Engwer, Christian; Wolters, Carsten H
2018-01-01
In Electro- (EEG) and Magnetoencephalography (MEG), one important requirement of source reconstruction is the forward model. The continuous Galerkin finite element method (CG-FEM) has become one of the dominant approaches for solving the forward problem over the last decades. Recently, a discontinuous Galerkin FEM (DG-FEM) EEG forward approach has been proposed as an alternative to CG-FEM (Engwer et al., 2017). It was shown that DG-FEM preserves the property of conservation of charge and that it can, in certain situations such as the so-called skull leakages , be superior to the standard CG-FEM approach. In this paper, we developed, implemented, and evaluated two DG-FEM approaches for the MEG forward problem, namely a conservative and a non-conservative one. The subtraction approach was used as source model. The validation and evaluation work was done in statistical investigations in multi-layer homogeneous sphere models, where an analytic solution exists, and in a six-compartment realistically shaped head volume conductor model. In agreement with the theory, the conservative DG-FEM approach was found to be superior to the non-conservative DG-FEM implementation. This approach also showed convergence with increasing resolution of the hexahedral meshes. While in the EEG case, in presence of skull leakages, DG-FEM outperformed CG-FEM, in MEG, DG-FEM achieved similar numerical errors as the CG-FEM approach, i.e., skull leakages do not play a role for the MEG modality. In particular, for the finest mesh resolution of 1 mm sources with a distance of 1.59 mm from the brain-CSF surface, DG-FEM yielded mean topographical errors (relative difference measure, RDM%) of 1.5% and mean magnitude errors (MAG%) of 0.1% for the magnetic field. However, if the goal is a combined source analysis of EEG and MEG data, then it is highly desirable to employ the same forward model for both EEG and MEG data. Based on these results, we conclude that the newly presented conservative DG
A Moving Discontinuous Galerkin Finite Element Method for Flows with Interfaces
2017-12-07
cuboid elements. In general, the method works in conjunction with standard local grid operations, including edge collapse, to ensure that degenerate...Shock: Line Plots............................................................................................ 31 iv TABLES 1 Intersecting Oblique...standard local grid operations including edge collapse, cf. [62], while in order to maintain a sufficient level of refinement, remaining cells can be refined
Directory of Open Access Journals (Sweden)
Seyed Abolfazl Hosseini
2016-02-01
Full Text Available In the present paper, development of the three-dimensional (3D computational code based on Galerkin finite element method (GFEM for solving the multigroup forward/adjoint diffusion equation in both rectangular and hexagonal geometries is reported. Linear approximation of shape functions in the GFEM with unstructured tetrahedron elements is used in the calculation. Both criticality and fixed source calculations may be performed using the developed GFEM-3D computational code. An acceptable level of accuracy at a low computational cost is the main advantage of applying the unstructured tetrahedron elements. The unstructured tetrahedron elements generated with Gambit software are used in the GFEM-3D computational code through a developed interface. The forward/adjoint multiplication factor, forward/adjoint flux distribution, and power distribution in the reactor core are calculated using the power iteration method. Criticality calculations are benchmarked against the valid solution of the neutron diffusion equation for International Atomic Energy Agency (IAEA-3D and Water-Water Energetic Reactor (VVER-1000 reactor cores. In addition, validation of the calculations against the P1 approximation of the transport theory is investigated in relation to the liquid metal fast breeder reactor benchmark problem. The neutron fixed source calculations are benchmarked through a comparison with the results obtained from similar computational codes. Finally, an analysis of the sensitivity of calculations to the number of elements is performed.
Miller, R.E.
1977-01-01
A steady-state simulation model was applied to the shallow hydrothermal system in the East Mesa area of Imperial Valley, Calif. The steady-state equations of flow and heat transport were solved by use of a Galerkin, finite-element method. A solution was obtained by iterating between the temperature and pressure equations, using updated densities and viscosities. Temperature and pressure were obtained for each node, and corresponding head values were calculated. The simulated temperature and pressure patterns correlated well with the observed patterns. Additional data, mainly from test drilling, would be required for construction of a similar model of the deep hydrothermal system.
Liu, G. R.; Nguyen-Xuan, H.; Nguyen-Thoi, T.; Xu, X.
2009-06-01
A carefully designed procedure is presented to modify the piecewise constant strain field of linear triangular FEM models, and to reconstruct a strain field with an adjustable parameter α. A novel Galerkin-like weakform derived from the Hellinger-Reissner variational principle is proposed for establishing the discretized system equations. The new weak form is very simple, possesses the same good properties of the standard Galerkin weakform, and works particularly well for strain construction methods. A superconvergent alpha finite element method (S αFEM) is then formulated by using the constructed strain field and the Galerkin-like weakform for solid mechanics problems. The implementation of the S αFEM is straightforward and no additional parameters are used. We prove theoretically and show numerically that the S αFEM always achieves more accurate and higher convergence rate than the standard FEM of triangular elements (T3) and even more accurate than the four-node quadrilateral elements (Q4) when the same sets of nodes are used. The S αFEM can always produce both lower and upper bounds to the exact solution in the energy norm for all elasticity problems by properly choosing an α. In addition, a preferable- α approach has also been devised to produce very accurate solutions for both displacement and energy norms and a superconvergent rate in the energy error norm. Furthermore, a model-based selective scheme is proposed to formulate a combined S αFEM/NS-FEM model that handily overcomes the volumetric locking problems. Intensive numerical studies including singularity problems have been conducted to confirm the theory and properties of the S αFEM.
Sirenko, Kostyantyn
2013-01-01
A scheme that discretizes exact absorbing boundary conditions (EACs) to incorporate them into a time-domain discontinuous Galerkin finite element method (TD-DG-FEM) is described. The proposed TD-DG-FEM with EACs is used for accurately characterizing transient electromagnetic wave interactions on two-dimensional waveguides. Numerical results demonstrate the proposed method\\'s superiority over the TD-DG-FEM that employs approximate boundary conditions and perfectly matched layers. Additionally, it is shown that the proposed method can produce the solution with ten-eleven digit accuracy when high-order spatial basis functions are used to discretize the Maxwell equations as well as the EACs. © 1963-2012 IEEE.
Quezada de Luna, M.; Farthing, M.; Guermond, J. L.; Kees, C. E.; Popov, B.
2017-12-01
The Shallow Water Equations (SWEs) are popular for modeling non-dispersive incompressible water waves where the horizontal wavelength is much larger than the vertical scales. They can be derived from the incompressible Navier-Stokes equations assuming a constant vertical velocity. The SWEs are important in Geophysical Fluid Dynamics for modeling surface gravity waves in shallow regimes; e.g., in the deep ocean. Some common geophysical applications are the evolution of tsunamis, river flooding and dam breaks, storm surge simulations, atmospheric flows and others. This work is concerned with the approximation of the time-dependent Shallow Water Equations with friction using explicit time stepping and continuous finite elements. The objective is to construct a method that is at least second-order accurate in space and third or higher-order accurate in time, positivity preserving, well-balanced with respect to rest states, well-balanced with respect to steady sliding solutions on inclined planes and robust with respect to dry states. Methods fulfilling the desired goals are common within the finite volume literature. However, to the best of our knowledge, schemes with the above properties are not well developed in the context of continuous finite elements. We start this work based on a finite element method that is second-order accurate in space, positivity preserving and well-balanced with respect to rest states. We extend it by: modifying the artificial viscosity (via the entropy viscosity method) to deal with issues of loss of accuracy around local extrema, considering a singular Manning friction term handled via an explicit discretization under the usual CFL condition, considering a water height regularization that depends on the mesh size and is consistent with the polynomial approximation, reducing dispersive errors introduced by lumping the mass matrix and others. After presenting the details of the method we show numerical tests that demonstrate the well
Kumar, Dinesh; Rai, K N
2016-12-01
Hyperthermia is a process that uses heat from the spatial heat source to kill cancerous cells without damaging the surrounding healthy tissues. Efficacy of hyperthermia technique is related to achieve temperature at the infected cells during the treatment process. A mathematical model on heat transfer in multilayer tissues in finite domain is proposed to predict the control temperature profile at hyperthermia position. The treatment technique uses dual-phase-lag model of heat transfer in multilayer tissues with modified Gaussian distribution heat source subjected to the most generalized boundary condition and interface at the adjacent layers. The complete dual-phase-lag model of bioheat transfer is solved using finite element Legendre wavelet Galerkin approach. The present solution has been verified with exact solution in a specific case and provides a good accuracy. The effect of the variability of different parameters such as lagging times, external heat source, metabolic heat source and the most generalized boundary condition on temperature profile in multilayer tissues is analyzed and also discussed the effective approach of hyperthermia treatment. Furthermore, we studied the modified thermal damage model with regeneration of healthy tissues as well. For viewpoint of thermal damage, the least thermal damage has been observed in boundary condition of second kind. The article concludes with a discussion of better opportunities for future clinical application of hyperthermia treatment. Copyright Â© 2016 Elsevier Ltd. All rights reserved.
quadratic spline finite element method
Directory of Open Access Journals (Sweden)
A. R. Bahadir
2002-01-01
Full Text Available The problem of heat transfer in a Positive Temperature Coefficient (PTC thermistor, which may form one element of an electric circuit, is solved numerically by a finite element method. The approach used is based on Galerkin finite element using quadratic splines as shape functions. The resulting system of ordinary differential equations is solved by the finite difference method. Comparison is made with numerical and analytical solutions and the accuracy of the computed solutions indicates that the method is well suited for the solution of the PTC thermistor problem.
Stable Galerkin versus equal-order Galerkin least-squares elements for the stokes flow problem
International Nuclear Information System (INIS)
Franca, L.P.; Frey, S.L.; Sampaio, R.
1989-11-01
Numerical experiments are performed for the stokes flow problem employing a stable Galerkin method and a Galerkin/Least-squares method with equal-order elements. Error estimates for the methods tested herein are reviewed. The numerical results presented attest the good stability properties of all methods examined herein. (A.C.A.S.) [pt
Continuum damage growth analysis using element free Galerkin ...
Indian Academy of Sciences (India)
This paper presents an elasto-plastic element free Galerkin formulation based on Newton–Raphson algorithm for damage growth analysis. Isotropic ductile damage evolution law is used. A study has been carried out in this paper using the proposed element free Galerkin method to understand the effect of initial damage ...
Shakib, Farzin; Hughes, Thomas J. R.
1991-01-01
A Fourier stability and accuracy analysis of the space-time Galerkin/least-squares method as applied to a time-dependent advective-diffusive model problem is presented. Two time discretizations are studied: a constant-in-time approximation and a linear-in-time approximation. Corresponding space-time predictor multi-corrector algorithms are also derived and studied. The behavior of the space-time algorithms is compared to algorithms based on semidiscrete formulations.
Directory of Open Access Journals (Sweden)
Wrana B.
2011-03-01
Full Text Available Artykuł podejmuje zagadnienie analizy rozchodzenia sie fal naprezeniowych w gruncie w ujeciu metody elementów skonczonych bazujac na sformułowaniu rozwiazania ciagłego w przestrzeni i nieciagłego w dziedzinie czasu Galerkina (space and time-discontinous Galerkin TDG finite element method. W tym sformułowaniu zarówno przemieszczenia jak i predkosci sa wielkosciami nieznanymi wzajemnie niezaleznymi aproksymowanymi ciagłymi funkcjami kształtu w przestrzeni i nieciagłymi funkcjami kształtu w czasie. Do opisu zachowania sie gruntu w pełni nasyconego woda zastosowano sformułowanie u-p w ujeciu metody elementów skonczonych. Grunt traktowany jest, jako osrodek dwufazowy składajacy sie ze szkieletu i wody w porach. Zastosowane sformułowanie uwzglednia tłumienie osrodka przez uwzglednienie dyssypacji energii proporcjonalnej do predkosci wody wzgledem szkieletu. W artykule przedstawiono porównanie proponowanej metody rozwiazania numerycznego w dziedzinie czasu do metod obecnie stosowanych, takich jak: metoda róznicy centralnej, metoda Houbolta, Wilsona θ, HHT-α oraz najczesciej stosowanej metody Newmarka. Z porównania wynika, ze proponowana metoda jest metoda stabilna o małym błedzie numerycznego rozwiazania.
Finite elements methods in mechanics
Eslami, M Reza
2014-01-01
This book covers all basic areas of mechanical engineering, such as fluid mechanics, heat conduction, beams, and elasticity with detailed derivations for the mass, stiffness, and force matrices. It is especially designed to give physical feeling to the reader for finite element approximation by the introduction of finite elements to the elevation of elastic membrane. A detailed treatment of computer methods with numerical examples are provided. In the fluid mechanics chapter, the conventional and vorticity transport formulations for viscous incompressible fluid flow with discussion on the method of solution are presented. The variational and Galerkin formulations of the heat conduction, beams, and elasticity problems are also discussed in detail. Three computer codes are provided to solve the elastic membrane problem. One of them solves the Poisson’s equation. The second computer program handles the two dimensional elasticity problems, and the third one presents the three dimensional transient heat conducti...
Electrical machine analysis using finite elements
Bianchi, Nicola
2005-01-01
OUTLINE OF ELECTROMAGNETIC FIELDSVector AnalysisElectromagnetic FieldsFundamental Equations SummaryReferencesBASIC PRINCIPLES OF FINITE ELEMENT METHODSIntroductionField Problems with Boundary ConditionsClassical Method for the Field Problem SolutionThe Classical Residual Method (Galerkin's Method)The Classical Variational Method (Rayleigh-Ritz's Method)The Finite Element MethodReferencesAPPLICATIONS OF THE FINITE ELEMENT METHOD TO TWO-DIMENSIONAL FIELDSIntroductionLinear Interpolation of the Function fApplication of the Variational MethodSimple Descriptions of Electromagnetic FieldsAppendix: I
The Blended Finite Element Method for Multi-fluid Plasma Modeling
2016-07-01
Briefing Charts 3. DATES COVERED (From - To) 07 June 2016 - 01 July 2016 4. TITLE AND SUBTITLE The Blended Finite Element Method for Multi-fluid Plasma...BLENDED FINITE ELEMENT METHOD FOR MULTI-FLUID PLASMA MODELING Éder M. Sousa1, Uri Shumlak2 1ERC INC., IN-SPACE PROPULSION BRANCH (RQRS) AIR FORCE RESEARCH...MULTI-FLUID PLASMA MODEL 2 BLENDED FINITE ELEMENT METHOD Blended Finite Element Method Nodal Continuous Galerkin Modal Discontinuous Galerkin Model
Analysis of Tube Drawing Process – A Finite Element Approach ...
African Journals Online (AJOL)
In this paper the effect of die semi angle on drawing load in cold tube drawing has been investigated numerically using the finite element method. The equation governing the stress distribution was derived and solved using Galerkin finite element method. An isoparametric formulation for the governing equation was utilized ...
International Nuclear Information System (INIS)
Lee, Byeong Hae
1992-02-01
This book gives descriptions of basic finite element method, which includes basic finite element method and data, black box, writing of data, definition of VECTOR, definition of matrix, matrix and multiplication of matrix, addition of matrix, and unit matrix, conception of hardness matrix like spring power and displacement, governed equation of an elastic body, finite element method, Fortran method and programming such as composition of computer, order of programming and data card and Fortran card, finite element program and application of nonelastic problem.
Continuum damage growth analysis using element free Galerkin ...
Indian Academy of Sciences (India)
using the proposed element free Galerkin method to understand the effect of initial damage and its growth on structural ... Effect of material discontinuity on damage growth analysis is also presented. Keywords. Damage ... 1996), Lagrange multiplier technique (Dolbow & Belytschko 1998), penalty method (Liu. 2002), full ...
Finite element analysis of nonlinear creeping flows
International Nuclear Information System (INIS)
Loula, A.F.D.; Guerreiro, J.N.C.
1988-12-01
Steady-state creep problems with monotone constitutive laws are studied. Finite element approximations are constructed based on mixed Petrov-Galerkin formulations for constrained problems. Stability, convergence and a priori error estimates are proved for equal-order discontinuous stress and continuous velocity interpolations. Numerical results are presented confirming the rates of convergence predicted in the analysis and the good performance of this formulation. (author) [pt
Superconvergence in Galerkin finite element methods
Wahlbin, Lars B
1995-01-01
This book is essentially a set of lecture notes from a graduate seminar given at Cornell in Spring 1994. It treats basic mathematical theory for superconvergence in the context of second order elliptic problems. It is aimed at graduate students and researchers. The necessary technical tools are developed in the text although sometimes long proofs are merely referenced. The book gives a rather complete overview of the field of superconvergence (in time-independent problems). It is the first text with such a scope. It includes a very complete and up-to-date list of references.
Peridynamic Multiscale Finite Element Methods
Energy Technology Data Exchange (ETDEWEB)
Costa, Timothy [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Bond, Stephen D. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Littlewood, David John [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Moore, Stan Gerald [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
2015-12-01
art of local models with the flexibility and accuracy of the nonlocal peridynamic model. In the mixed locality method this coupling occurs across scales, so that the nonlocal model can be used to communicate material heterogeneity at scales inappropriate to local partial differential equation models. Additionally, the computational burden of the weak form of the peridynamic model is reduced dramatically by only requiring that the model be solved on local patches of the simulation domain which may be computed in parallel, taking advantage of the heterogeneous nature of next generation computing platforms. Addition- ally, we present a novel Galerkin framework, the 'Ambulant Galerkin Method', which represents a first step towards a unified mathematical analysis of local and nonlocal multiscale finite element methods, and whose future extension will allow the analysis of multiscale finite element methods that mix models across scales under certain assumptions of the consistency of those models.
Ablative Thermal Response Analysis Using the Finite Element Method
Dec John A.; Braun, Robert D.
2009-01-01
A review of the classic techniques used to solve ablative thermal response problems is presented. The advantages and disadvantages of both the finite element and finite difference methods are described. As a first step in developing a three dimensional finite element based ablative thermal response capability, a one dimensional computer tool has been developed. The finite element method is used to discretize the governing differential equations and Galerkin's method of weighted residuals is used to derive the element equations. A code to code comparison between the current 1-D tool and the 1-D Fully Implicit Ablation and Thermal Response Program (FIAT) has been performed.
Generalized multiscale finite element method. Symmetric interior penalty coupling
Efendiev, Yalchin R.
2013-12-01
Motivated by applications to numerical simulations of flows in highly heterogeneous porous media, we develop multiscale finite element methods for second order elliptic equations. We discuss a multiscale model reduction technique in the framework of the discontinuous Galerkin finite element method. We propose two different finite element spaces on the coarse mesh. The first space is based on a local eigenvalue problem that uses an interior weighted L2-norm and a boundary weighted L2-norm for computing the "mass" matrix. The second choice is based on generation of a snapshot space and subsequent selection of a subspace of a reduced dimension. The approximation with these multiscale spaces is based on the discontinuous Galerkin finite element method framework. We investigate the stability and derive error estimates for the methods and further experimentally study their performance on a representative number of numerical examples. © 2013 Elsevier Inc.
Finite elements and approximation
Zienkiewicz, O C
2006-01-01
A powerful tool for the approximate solution of differential equations, the finite element is extensively used in industry and research. This book offers students of engineering and physics a comprehensive view of the principles involved, with numerous illustrative examples and exercises.Starting with continuum boundary value problems and the need for numerical discretization, the text examines finite difference methods, weighted residual methods in the context of continuous trial functions, and piecewise defined trial functions and the finite element method. Additional topics include higher o
Solution of Differential Equation by Means of Finite Element Method
Hayashi, Daigoro; 林, 大五郎
1989-01-01
The finite element method may be defined as the modern "Method of Weighted Residuals" (MWR). This paper describes how to solve the differential equations which are essential in order to explain quantiatively a number of valuable geological and geodynamic problems.The methods to solve linear differential equation, non-linear equation, non-linear non-steady equation, Laplace equation and incompreeible New tonian flow problem are explained by means of the Galerkin finite element method.
Stabilized Finite Elements in FUN3D
Anderson, W. Kyle; Newman, James C.; Karman, Steve L.
2017-01-01
A Streamlined Upwind Petrov-Galerkin (SUPG) stabilized finite-element discretization has been implemented as a library into the FUN3D unstructured-grid flow solver. Motivation for the selection of this methodology is given, details of the implementation are provided, and the discretization for the interior scheme is verified for linear and quadratic elements by using the method of manufactured solutions. A methodology is also described for capturing shocks, and simulation results are compared to the finite-volume formulation that is currently the primary method employed for routine engineering applications. The finite-element methodology is demonstrated to be more accurate than the finite-volume technology, particularly on tetrahedral meshes where the solutions obtained using the finite-volume scheme can suffer from adverse effects caused by bias in the grid. Although no effort has been made to date to optimize computational efficiency, the finite-element scheme is competitive with the finite-volume scheme in terms of computer time to reach convergence.
Numerical simulation of fluid-structure interactions with stabilized finite element method
Sváček, Petr
2016-03-01
This paper is interested to the interactions of the incompressible flow with a flexibly supported airfoil. The bending and the torsion modes are considered. The problem is mathematically described. The numerical method is based on the finite element method. A combination of the streamline-upwind/Petrov-Galerkin and pressure stabilizing/Petrov-Galerkin method is used for the stabilization of the finite element method. The numerical results for a three-dimensional problem of flow over an airfoil are shown.
Finite Element Method for Linear Multiterm Fractional Differential Equations
Directory of Open Access Journals (Sweden)
Abdallah A. Badr
2012-01-01
Full Text Available We consider the linear multiterm fractional differential equation (fDE. Existence and uniqueness of the solution of such equation are discussed. We apply the finite element method (FEM to obtain the numerical solution of this equation using Galerkin approach. A comparison, through examples, between our techniques and other previous numerical methods is established.
Finite Element Method for Linear Multiterm Fractional Differential Equations
Badr, Abdallah A.
2012-01-01
We consider the linear multiterm fractional differential equation (fDE). Existence and uniqueness of the solution of such equation are discussed. We apply the finite element method (FEM) to obtain the numerical solution of this equation using Galerkin approach. A comparison, through examples, between our techniques and other previous numerical methods is established.
(ajst) finite element analysis of a fluid-structure
African Journals Online (AJOL)
liquid flow. The fluid-structure interaction is found to be governed by Poisson's ratio. In this steady finite element method based on Galerkin formulation is applied. Numerical results show a good similarity with those of the literature obtained by the characteristics method. Key words : Fluid-structure interaction, flexible pipe, ...
Finite element approximations of the stokes flow problem based upon various variational principles
International Nuclear Information System (INIS)
Franca, L.P.; Hughers, T.J.R.; Stenberg, R.
1989-05-01
Finite element methods are constructed by adding to the usual Galerkin method terms that are mesh-dependent least-squares forms of the Euler-Lagrange equations. The methods are consistent and possess additional stability compared to the Galerkin method. Finite element interpolations, which are unstable in the Galerkin approach, are now convergent. The methodology is applied to the velocity-pressure formulation, a.k.a., Herrmann's formulation, to the stress-velocity formulation, a.k.a., Hellinger-Reissner's formulation and to a new formulation based on augmented stress, pressure and velocity [pt
Weiser, Martin
2016-01-01
All relevant implementation aspects of finite element methods are discussed in this book. The focus is on algorithms and data structures as well as on their concrete implementation. Theory is covered as far as it gives insight into the construction of algorithms. Throughout the exercises a complete FE-solver for scalar 2D problems will be implemented in Matlab/Octave.
Element free Galerkin formulation of composite beam with longitudinal slip
Energy Technology Data Exchange (ETDEWEB)
Ahmad, Dzulkarnain; Mokhtaram, Mokhtazul Haizad [Department of Civil Engineering, Universiti Selangor, Bestari Jaya, Selangor (Malaysia); Badli, Mohd Iqbal; Yassin, Airil Y. Mohd [Faculty of Civil Engineering, Universiti Teknologi Malaysia, Skudai, Johor (Malaysia)
2015-05-15
Behaviour between two materials in composite beam is assumed partially interact when longitudinal slip at its interfacial surfaces is considered. Commonly analysed by the mesh-based formulation, this study used meshless formulation known as Element Free Galerkin (EFG) method in the beam partial interaction analysis, numerically. As meshless formulation implies that the problem domain is discretised only by nodes, the EFG method is based on Moving Least Square (MLS) approach for shape functions formulation with its weak form is developed using variational method. The essential boundary conditions are enforced by Langrange multipliers. The proposed EFG formulation gives comparable results, after been verified by analytical solution, thus signify its application in partial interaction problems. Based on numerical test results, the Cubic Spline and Quartic Spline weight functions yield better accuracy for the EFG formulation, compares to other proposed weight functions.
Finite element solution algorithm for incompressible fluid dynamics
Baker, A. J.
1974-01-01
A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the transient motion of a viscous incompressible fluid, i.e., hydrodynamics. Dependent variable transformation renders the differential equation description uniformly elliptic. The finite element algorithm is established using the Galerkin criterion on a local basis within the Method of Weighted Residuals. It is unconstrained with respect to system linearity, computational mesh uniformity or solution domain closure regularity. The finite element matrices are established using a linear 'natural coordinate function' description. Computational solutions using the COMOC computer program illustrate the various features of the algorithm including recirculating flows.
Directory of Open Access Journals (Sweden)
Samira Hosseini
Full Text Available Abstract One of the main drawbacks of Element Free Galerkin (EFG method is its dependence on moving least square shape functions which don’t satisfy the Kronecker Delta property, so in this method it’s not possible to apply Dirichlet boundary conditions directly. The aim of the present paper is to discuss different aspects of three widely used methods of applying Dirichlet boundary conditions in EFG method, called Lagrange multipliers, penalty method, and coupling with finite element method. Numerical simulations are presented to compare the results of these methods form the perspective of accuracy, convergence and computational expense. These methods have been implemented in an object oriented programing environment, called INSANE, and the results are presented and compared with the analytical solutions.
Directory of Open Access Journals (Sweden)
M.H.R. Ghoreishy
2008-02-01
Full Text Available This research work is devoted to the footprint analysis of a steel-belted radial tyre (185/65R14 under vertical static load using finite element method. Two models have been developed in which in the first model the tread patterns were replaced by simple ribs while the second model was consisted of details of the tread blocks. Linear elastic and hyper elastic (Arruda-Boyce material models were selected to describe the mechanical behavior of the reinforcing and rubbery parts, respectively. The above two finite element models of the tyre were analyzed under inflation pressure and vertical static loads. The second model (with detailed tread patterns was analyzed with and without friction effect between tread and contact surfaces. In every stage of the analysis, the results were compared with the experimental data to confirm the accuracy and applicability of the model. Results showed that neglecting the tread pattern design not only reduces the computational cost and effort but also the differences between computed deformations do not show significant changes. However, more complicated variables such as shape and area of the footprint zone and contact pressure are affected considerably by the finite element model selected for the tread blocks. In addition, inclusion of friction even in static state changes these variables significantly.
A finite element primer for beginners the basics
Zohdi, Tarek I
2014-01-01
The purpose of this primer is to provide the basics of the Finite Element Method, primarily illustrated through a classical model problem, linearized elasticity. The topics covered are:(1) Weighted residual methods and Galerkin approximations,(2) A model problem for one-dimensional?linear elastostatics,(3) Weak formulations in one dimension,(4) Minimum principles in one dimension,(5) Error estimation in one dimension,(5) Construction of Finite Element basis functions in one dimension,(6) Gaussian Quadrature,(7) Iterative solvers and element by element data structures,(8) A model problem for th
International Nuclear Information System (INIS)
Pereira, Luis Carlos Martins
1998-06-01
New Petrov-Galerkin formulations on the finite element methods for convection-diffusion problems with boundary layers are presented. Such formulations are based on a consistent new theory on discontinuous finite element methods. Existence and uniqueness of solutions for these problems in the new finite element spaces are demonstrated. Some numerical experiments shows how the new formulation operate and also their efficacy. (author)
International Nuclear Information System (INIS)
Franca, L.P.; Toledo, E.M.; Loula, A.F.D.; Garcia, E.L.M.
1988-12-01
A new finite element method is employed to approximate axisymmetric shell problems. This formulation enhances stability and accuracy, from thin to moderately thick shells, compared to the correspondent Galerkin finite element approximations. Numerical results illustrate the good performance of the present method on some typical pressure vessels aplications. (author) [pt
International Nuclear Information System (INIS)
Al-Akhrass, Dina
2014-01-01
Simulations in solid mechanics exhibit several difficulties, as dealing with incompressibility, with nonlinearities due to finite strains, contact laws, or constitutive laws. The basic motivation of our work is to propose efficient finite element methods capable of dealing with incompressibility in finite strain context, and using elements of low order. During the three last decades, many approaches have been proposed in the literature to overcome the incompressibility problem. Among them, mixed formulations offer an interesting theoretical framework. In this work, a three-field mixed formulation (displacement, pressure, volumetric strain) is investigated. In some cases, this formulation can be condensed in a two-field (displacement - pressure) mixed formulation. However, it is well-known that the discrete problem given by the Galerkin finite element technique, does not inherit the 'inf-sup' stability condition from the continuous problem. Hence, the interpolation orders in displacement and pressure have to be chosen in a way to satisfy the Brezzi-Babuska stability conditions when using Galerkin approaches. Interpolation orders must be chosen so as to satisfy this condition. Two possibilities are considered: to use stable finite element satisfying this requirement, or to use finite element that does not satisfy this condition, and to add terms stabilizing the FE Galerkin formulation. The latter approach allows the use of equal order interpolation. In this work, stable finite element P2/P1 and P2/P1/P1 are used as reference, and compared to P1/P1 and P1/P1/P1 formulations stabilized with a bubble function or with a VMS method (Variational Multi-Scale) based on a sub-grid-space orthogonal to the FE space. A finite strain model based on logarithmic strain is selected. This approach is extended to three and two field mixed formulations with stable or stabilized elements. These approaches are validated on academic cases and used on industrial cases. (author)
Generalized multiscale finite element methods. nonlinear elliptic equations
Efendiev, Yalchin R.
2013-01-01
In this paper we use the Generalized Multiscale Finite Element Method (GMsFEM) framework, introduced in [26], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation. With this convention, we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin (CG) or discontinuous Galerkin (DG) global formulations. Here, we use Symmetric Interior Penalty Discontinuous Galerkin approach. Both methods yield a predictable error decline that depends on the respective coarse space dimension, and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples. © 2014 Global-Science Press.
The dimension split element-free Galerkin method for three-dimensional potential problems
Meng, Z. J.; Cheng, H.; Ma, L. D.; Cheng, Y. M.
2018-02-01
This paper presents the dimension split element-free Galerkin (DSEFG) method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares (IMLS) approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
Optical Finite Element Processor
Casasent, David; Taylor, Bradley K.
1986-01-01
A new high-accuracy optical linear algebra processor (OLAP) with many advantageous features is described. It achieves floating point accuracy, handles bipolar data by sign-magnitude representation, performs LU decomposition using only one channel, easily partitions and considers data flow. A new application (finite element (FE) structural analysis) for OLAPs is introduced and the results of a case study presented. Error sources in encoded OLAPs are addressed for the first time. Their modeling and simulation are discussed and quantitative data are presented. Dominant error sources and the effects of composite error sources are analyzed.
International Nuclear Information System (INIS)
Bailey, T S; Adams, M L; Yang, B; Zika, M R
2005-01-01
We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses piecewise linear weight and basis functions in the finite element approximation, and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We show that this new PWL method gives solutions comparable to those from Palmer's finite-volume method. However, since the PWL method produces a symmetric positive definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids
International Nuclear Information System (INIS)
Bailey, T.S.; Adams, M.L.; Yang, B.; Zika, M.R.
2005-01-01
We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses piecewise linear weight and basis functions in the finite element approximation, and it can be applied on arbitrary polygonal (2-dimensional) or polyhedral (3-dimensional) grids. We show that this new PWL method gives solutions comparable to those from Palmer's finite-volume method. However, since the PWL method produces a symmetric positive definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids. (authors)
High speed inviscid compressible flow by the finite element method
Zienkiewicz, O. C.; Loehner, R.; Morgan, K.
1984-01-01
The finite element method and an explicit time stepping algorithm which is based on Taylor-Galerkin schemes with an appropriate artificial viscosity is combined with an automatic mesh refinement process which is designed to produce accurate steady state solutions to problems of inviscid compressible flow in two dimensions. The results of two test problems are included which demonstrate the excellent performance characteristics of the proposed procedures.
A code for obtaining temperature distribution by finite element method
International Nuclear Information System (INIS)
Bloch, M.
1984-01-01
The ELEFIB Fortran language computer code using finite element method for calculating temperature distribution of linear and two dimensional problems, in permanent region or in the transient phase of heat transfer, is presented. The formulation of equations uses the Galerkin method. Some examples are shown and the results are compared with other papers. The comparative evaluation shows that the elaborated code gives good values. (M.C.K.) [pt
Convergence of a residual based artificial viscosity finite element method
Nazarov, Murtazo
2013-02-01
We present a residual based artificial viscosity finite element method to solve conservation laws. The Galerkin approximation is stabilized by only residual based artificial viscosity, without any least-squares, SUPG, or streamline diffusion terms. We prove convergence of the method, applied to a scalar conservation law in two space dimensions, toward an unique entropy solution for implicit time stepping schemes. © 2012 Elsevier B.V. All rights reserved.
Sparse adaptive finite elements for radiative transfer
International Nuclear Information System (INIS)
Widmer, G.; Hiptmair, R.; Schwab, Ch.
2008-01-01
The linear radiative transfer equation, a partial differential equation for the radiation intensity u(x,s), with independent variables x element of D is contained in R n in the physical domain D of dimension n=2,3, and angular variable s element of S 2 :={y element of R 3 :|y|=1}, is solved in the n+2-dimensional computational domain DxS 2 . We propose an adaptive multilevel Galerkin finite element method (FEM) for its numerical solution. Our approach is based on (a) a stabilized variational formulation of the transport operator, (b) on so-called sparse tensor products of two hierarchic families of finite element spaces in H 1 (D) and in L 2 (S 2 ), respectively, and (c) on wavelet thresholding techniques to adapt the discretization to the underlying problem. An a priori error analysis shows, under strong regularity assumptions on the solution, that the sparse tensor product method is clearly superior to a discrete ordinates method, as it converges with essentially optimal asymptotic rates while its complexity grows essentially only as that for a linear transport problem in R n . Numerical experiments for n=2 on a set of example problems agree with the convergence and complexity analysis of the method and show that introducing adaptivity can improve performance in terms of accuracy vs. number of degrees even further
Statistical finite element analysis.
Khalaji, Iman; Rahemifar, Kaamran; Samani, Abbas
2008-01-01
A novel technique is introduced for tissue deformation and stress analysis. Compared to the conventional Finite Element method, this technique is orders of magnitude faster and yet still very accurate. The proposed technique uses preprocessed data obtained from FE analyses of a number of similar objects in a Statistical Shape Model framework as described below. This technique takes advantage of the fact that the body organs have limited variability, especially in terms of their geometry. As such, it is well suited for calculating tissue displacements of body organs. The proposed technique can be applied in many biomedical applications such as image guided surgery, or virtual reality environment development where tissue behavior is simulated for training purposes.
Adaptive finite element methods for differential equations
Bangerth, Wolfgang
2003-01-01
These Lecture Notes discuss concepts of `self-adaptivity' in the numerical solution of differential equations, with emphasis on Galerkin finite element methods. The key issues are a posteriori error estimation and it automatic mesh adaptation. Besides the traditional approach of energy-norm error control, a new duality-based technique, the Dual Weighted Residual method for goal-oriented error estimation, is discussed in detail. This method aims at economical computation of arbitrary quantities of physical interest by properly adapting the computational mesh. This is typically required in the design cycles of technical applications. For example, the drag coefficient of a body immersed in a viscous flow is computed, then it is minimized by varying certain control parameters, and finally the stability of the resulting flow is investigated by solving an eigenvalue problem. `Goal-oriented' adaptivity is designed to achieve these tasks with minimal cost. At the end of each chapter some exercises are posed in order ...
Finite element analysis of rotating beams physics based interpolation
Ganguli, Ranjan
2017-01-01
This book addresses the solution of rotating beam free-vibration problems using the finite element method. It provides an introduction to the governing equation of a rotating beam, before outlining the solution procedures using Rayleigh-Ritz, Galerkin and finite element methods. The possibility of improving the convergence of finite element methods through a judicious selection of interpolation functions, which are closer to the problem physics, is also addressed. The book offers a valuable guide for students and researchers working on rotating beam problems – important engineering structures used in helicopter rotors, wind turbines, gas turbines, steam turbines and propellers – and their applications. It can also be used as a textbook for specialized graduate and professional courses on advanced applications of finite element analysis.
International Nuclear Information System (INIS)
Thornton, E.A.; Ramakrishnan, R.
1986-06-01
Prediction of compressible flow phenomena using the finite element method is of recent origin and considerable interest. Two shock capturing finite element formulations for high speed compressible flows are described. A Taylor-Galerkin formulation uses a Taylor series expansion in time coupled with a Galerkin weighted residual statement. The Taylor-Galerkin algorithms use explicit artificial dissipation, and the performance of three dissipation models are compared. A Petrov-Galerkin algorithm has as its basis the concepts of streamline upwinding. Vectorization strategies are developed to implement the finite element formulations on the NASA Langley VPS-32. The vectorization scheme results in finite element programs that use vectors of length of the order of the number of nodes or elements. The use of the vectorization procedure speeds up processing rates by over two orders of magnitude. The Taylor-Galerkin and Petrov-Galerkin algorithms are evaluated for 2D inviscid flows on criteria such as solution accuracy, shock resolution, computational speed and storage requirements. The convergence rates for both algorithms are enhanced by local time-stepping schemes. Extension of the vectorization procedure for predicting 2D viscous and 3D inviscid flows are demonstrated. Conclusions are drawn regarding the applicability of the finite element procedures for realistic problems that require hundreds of thousands of nodes
A combined continuous-discontinuous finite element method for convection-diffusion problems
Devloo, PRB; Forti, T; Gomes, SM
2007-01-01
Discontinuous Galerkin (DGM) method combines the advantages of stability of finite volume method and the accuracy of continuous finite element method (FEIVI). Applications of the DGM are particularly valuable where the solution presents high-gradients or discontinuities, such as boundary layers and shock problems. A disadvantage of the DGM is the higher computational cost when comparing to classic finite element method, due to the increased number of degrees of freedom. With this motivation, ...
Solution of Finite Element Equations
DEFF Research Database (Denmark)
Krenk, Steen
An important step in solving any problem by the finite element method is the solution of the global equations. Numerical solution of linear equations is a subject covered in most courses in numerical analysis. However, the equations encountered in most finite element applications have some special...... features that justify the development of specialized solution algorithms....
Massively Parallel Finite Element Programming
Heister, Timo
2010-01-01
Today\\'s large finite element simulations require parallel algorithms to scale on clusters with thousands or tens of thousands of processor cores. We present data structures and algorithms to take advantage of the power of high performance computers in generic finite element codes. Existing generic finite element libraries often restrict the parallelization to parallel linear algebra routines. This is a limiting factor when solving on more than a few hundreds of cores. We describe routines for distributed storage of all major components coupled with efficient, scalable algorithms. We give an overview of our effort to enable the modern and generic finite element library deal.II to take advantage of the power of large clusters. In particular, we describe the construction of a distributed mesh and develop algorithms to fully parallelize the finite element calculation. Numerical results demonstrate good scalability. © 2010 Springer-Verlag.
Block-iterative finite element computations for incompressible flow problems
International Nuclear Information System (INIS)
Tezduyar, T.E.; Liou, J.; Glowinski, R.; Nguyen, T.; Poole, S.
1988-01-01
A block-iterative finite element procedure is presented for two-dimensional fluid dynamics computations on multiply-connected domains based on the vorticity-stream function formulation of the incompressible Navier-Stokes equations. The difficulty associated with the convection term in the vorticity transport equation is addressed by using a streamline-upwind/Petrov-Galerkin scheme. Element-by-element preconditioned iteration techniques with high degree of vectorization and high computational speed are employed to solve the linear equation system for each block. The authors conclude that performance evaluations show the potential of these techniques to be used for large-scale computations
Galerkin and weighted Galerkin methods for a forward-backward heat equation
Lu, H.
1997-01-01
Galerkin and weighted Galerkin methods are proposed for the numerical solution of parabolic partial differential equations where the diffusion coefficient takes different signs. The approach is based on a simultaneous discretization of space and time variables by using continuous finite element
The finite element response Matrix method
International Nuclear Information System (INIS)
Nakata, H.; Martin, W.R.
1983-01-01
A new method for global reactor core calculations is described. This method is based on a unique formulation of the response matrix method, implemented with a higher order finite element method. The unique aspects of this approach are twofold. First, there are two levels to the overall calculational scheme: the local or assembly level and the global or core level. Second, the response matrix scheme, which is formulated at both levels, consists of two separate response matrices rather than one response matrix as is generally the case. These separate response matrices are seen to be quite beneficial for the criticality eigenvalue calculation, because they are independent of k /SUB eff/. The response matrices are generated from a Galerkin finite element solution to the weak form of the diffusion equation, subject to an arbitrary incoming current and an arbitrary distributed source. Calculational results are reported for two test problems, the two-dimensional International Atomic Energy Agency benchmark problem and a two-dimensional pressurized water reactor test problem (Biblis reactor), and they compare well with standard coarse mesh methods with respect to accuracy and efficiency. Moreover, the accuracy (and capability) is comparable to fine mesh for a fraction of the computational cost. Extension of the method to treat heterogeneous assemblies and spatial depletion effects is discussed
Generalized multiscale finite element method for elasticity equations
Chung, Eric T.
2014-10-05
In this paper, we discuss the application of generalized multiscale finite element method (GMsFEM) to elasticity equation in heterogeneous media. We consider steady state elasticity equations though some of our applications are motivated by elastic wave propagation in subsurface where the subsurface properties can be highly heterogeneous and have high contrast. We present the construction of main ingredients for GMsFEM such as the snapshot space and offline spaces. The latter is constructed using local spectral decomposition in the snapshot space. The spectral decomposition is based on the analysis which is provided in the paper. We consider both continuous Galerkin and discontinuous Galerkin coupling of basis functions. Both approaches have their cons and pros. Continuous Galerkin methods allow avoiding penalty parameters though they involve partition of unity functions which can alter the properties of multiscale basis functions. On the other hand, discontinuous Galerkin techniques allow gluing multiscale basis functions without any modifications. Because basis functions are constructed independently from each other, this approach provides an advantage. We discuss the use of oversampling techniques that use snapshots in larger regions to construct the offline space. We provide numerical results to show that one can accurately approximate the solution using reduced number of degrees of freedom.
Finite Element Based Formulation of Lattice Boltzmann Equation
International Nuclear Information System (INIS)
Jo, Jong Chull; Roh, Kyung Wan; Kwon, Young W.; Kwon, Young W.
2008-01-01
The Lattice Boltzmann Method (LBM) has been developed for application to thermal-fluid problems. Recently, the technique was also applied to fluid-structure interaction problems. Most of those studies considered a regular shape of lattice or mesh like square and cubic grids. In order to apply the LBM to more practical cases, it is necessary to be able to solve complex or irregular shapes of problem domains. There have been different kinds of approaches to address the problems. The most common technique was using the finite volume formulation of the lattice Boltzmann equation. Another approach was a point-wise interpolation technique for irregular grids. Other techniques were based on the finite element method. Generally, the finite element method is very powerful for solving two or three-dimensional complex or irregular shapes of domains using the isoparametric element formulation which is based on a mathematical mapping from a regular shape of element in an imaginary domain to a more general and irregular shape of element in the physical domain. In addition, there are variety of choices of finite elements such as triangular or quadrilateral shapes in 2-D, or tetrahedral, triangular prism, or general six-sided solids in 3-D. As a result, the present study presents a new finite element formulation for the lattice Boltzmann equation using the general weighted residual technique. Among the weighted residual formulations, the collocation method, Galerkin method or method of moments are used to develop the finite element based LBM
Programming the finite element method
Smith, I M; Margetts, L
2013-01-01
Many students, engineers, scientists and researchers have benefited from the practical, programming-oriented style of the previous editions of Programming the Finite Element Method, learning how to develop computer programs to solve specific engineering problems using the finite element method. This new fifth edition offers timely revisions that include programs and subroutine libraries fully updated to Fortran 2003, which are freely available online, and provides updated material on advances in parallel computing, thermal stress analysis, plasticity return algorithms, convection boundary c
On symmetric pyramidal finite elements
Czech Academy of Sciences Publication Activity Database
Liu, L.; Davies, K. B.; Yuan, K.; Křížek, Michal
2004-01-01
Roč. 11, 1-2 (2004), s. 213-227 ISSN 1492-8760 R&D Projects: GA AV ČR IAA1019201 Institutional research plan: CEZ:AV0Z1019905 Keywords : mesh generation * finite element method * composite elements Subject RIV: BA - General Mathematics Impact factor: 0.108, year: 2004
Finite element methods for engineers
Fenner, Roger T
2013-01-01
This book is intended as a textbook providing a deliberately simple introduction to finite element methods in a way that should be readily understandable to engineers, both students and practising professionals. Only the very simplest elements are considered, mainly two dimensional three-noded “constant strain triangles”, with simple linear variation of the relevant variables. Chapters of the book deal with structural problems (beams), classification of a broad range of engineering into harmonic and biharmonic types, finite element analysis of harmonic problems, and finite element analysis of biharmonic problems (plane stress and plane strain). Full Fortran programs are listed and explained in detail, and a range of practical problems solved in the text. Despite being somewhat unfashionable for general programming purposes, the Fortran language remains very widely used in engineering. The programs listed, which were originally developed for use on mainframe computers, have been thoroughly updated for use ...
Finite element analysis of a fluid-structure interaction in flexible pipe ...
African Journals Online (AJOL)
The obtained mathematical system is constituted of four non-linear hyperbolic partial differential equations describing the wave propagation in both pipe wall and liquid flow. The fluid-structure interaction is found to be governed by Poisson's ratio. In this steady finite element method based on Galerkin formulation is applied.
Time-integration methods for finite element discretisations of the second-order Maxwell equation
Sarmany, D.; Bochev, Mikhail A.; van der Vegt, Jacobus J.W.
This article deals with time integration for the second-order Maxwell equations with possibly non-zero conductivity in the context of the discontinuous Galerkin finite element method (DG-FEM) and the H(curl)-conforming FEM. For the spatial discretisation, hierarchic H(curl)-conforming basis
Time-integration methods for finite element discretisations of the second-order Maxwell equation
Sarmany, D.; Bochev, Mikhail A.; van der Vegt, Jacobus J.W.
This article deals with time integration for the second-order Maxwell equations with possibly non-zero conductivity in the context of the discontinuous Galerkin finite element method DG-FEM) and the $H(\\mathrm{curl})$-conforming FEM. For the spatial discretisation, hierarchic
International Nuclear Information System (INIS)
Bailey, Teresa S.; Adams, Marvin L.; Yang, Brian; Zika, Michael R.
2008-01-01
We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses recently introduced piecewise linear weight and basis functions in the finite element approximation and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We first demonstrate some analytical properties of the PWL method and perform a simple mode analysis to compare the PWL method with Palmer's vertex-centered finite-volume method and with a bilinear continuous finite element method. We then show that this new PWL method gives solutions comparable to those from Palmer's. However, since the PWL method produces a symmetric positive-definite coefficient matrix, it should be substantially more computationally efficient than Palmer's method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids
Membrane finite element method for simulating fluid flow in porous medium
Directory of Open Access Journals (Sweden)
Mei-li Zhan
2009-06-01
Full Text Available A new membrane finite element method for modeling fluid flow in a porous medium is presented in order to quickly and accurately simulate the geo-membrane fabric used in civil engineering. It is based on discontinuous finite element theory, and can be easily coupled with the normal Galerkin finite element method. Based on the saturated seepage equation, the element coefficient matrix of the membrane element method is derived, and a geometric transform relation for the membrane element between a global coordinate system and a local coordinate system is obtained. A method for the determination of the fluid flux conductivity of the membrane element is presented. This method provides a basis for determining discontinuous parameters in discontinuous finite element theory. An anti-seepage problem regarding the foundation of a building is analyzed by coupling the membrane finite element method with the normal Galerkin finite element method. The analysis results demonstrate the utility and superiority of the membrane finite element method in fluid flow analysis of a porous medium.
The Iris biometric feature segmentation using finite element method
Directory of Open Access Journals (Sweden)
David Ibitayo LANLEGE
2015-05-01
Full Text Available This manuscript presents a method for segmentation of iris images based on a deformable contour (active contour paradigm. The deformable contour is a novel approach in image segmentation. A type of active contour is the Snake. Snake is a parametric curve defined within the domain of the image. Snake properties are specified through a function called energy functional. This means they consist of packets of energy which expressed as partial Differential Equations. The partial Differential Equation is the controlling engine of the active contour since this project, the Finite Element Method (Standard Galerkin Method implementation for deformable model is presented.
A finite element method for extended KdV equations
Directory of Open Access Journals (Sweden)
Karczewska Anna
2016-09-01
Full Text Available The finite element method (FEM is applied to obtain numerical solutions to a recently derived nonlinear equation for the shallow water wave problem. A weak formulation and the Petrov–Galerkin method are used. It is shown that the FEM gives a reasonable description of the wave dynamics of soliton waves governed by extended KdV equations. Some new results for several cases of bottom shapes are presented. The numerical scheme presented here is suitable for taking into account stochastic effects, which will be discussed in a subsequent paper.
Finite element method for simulation of the semiconductor devices
International Nuclear Information System (INIS)
Zikatanov, L.T.; Kaschiev, M.S.
1991-01-01
An iterative method for solving the system of nonlinear equations of the drift-diffusion representation for the simulation of the semiconductor devices is worked out. The Petrov-Galerkin method is taken for the discretization of these equations using the bilinear finite elements. It is shown that the numerical scheme is a monotonous one and there are no oscillations of the solutions in the region of p-n transition. The numerical calculations of the simulation of one semiconductor device are presented. 13 refs.; 3 figs
FINITE ELEMENT ANALYSIS OF STRUCTURES
Directory of Open Access Journals (Sweden)
PECINGINA OLIMPIA-MIOARA
2015-05-01
Full Text Available The application of finite element method is analytical when solutions can not be applied for deeper study analyzes static, dynamic or other types of requirements in different points of the structures .In practice it is necessary to know the behavior of the structure or certain parts components of the machine under the influence of certain factors static and dynamic . The application of finite element in the optimization of components leads to economic growth , to increase reliability and durability organs studied, thus the machine itself.
An enriched finite element method to fractional advection-diffusion equation
Luan, Shengzhi; Lian, Yanping; Ying, Yuping; Tang, Shaoqiang; Wagner, Gregory J.; Liu, Wing Kam
2017-08-01
In this paper, an enriched finite element method with fractional basis [ 1,x^{α }] for spatial fractional partial differential equations is proposed to obtain more stable and accurate numerical solutions. For pure fractional diffusion equation without advection, the enriched Galerkin finite element method formulation is demonstrated to simulate the exact solution successfully without any numerical oscillation, which is advantageous compared to the traditional Galerkin finite element method with integer basis [ 1,x] . For fractional advection-diffusion equation, the oscillatory behavior becomes complex due to the introduction of the advection term which can be characterized by a fractional element Peclet number. For the purpose of addressing the more complex numerical oscillation, an enriched Petrov-Galerkin finite element method is developed by using a dimensionless fractional stabilization parameter, which is formulated through a minimization of the residual of the nodal solution. The effectiveness and accuracy of the enriched finite element method are demonstrated by a series of numerical examples of fractional diffusion equation and fractional advection-diffusion equation, including both one-dimensional and two-dimensional, steady-state and time-dependent cases.
On Using Particle Finite Element for Hydrodynamics Problems Solving
Directory of Open Access Journals (Sweden)
E. V. Davidova
2015-01-01
Full Text Available The aim of the present research is to develop software for the Particle Finite Element Method (PFEM and its verification on the model problem of viscous incompressible flow simulation in a square cavity. The Lagrangian description of the medium motion is used: the nodes of the finite element mesh move together with the fluid that allows to consider them as particles of the medium. Mesh cells deform when in time-stepping procedure, so it is necessary to reconstruct the mesh to provide stability of the finite element numerical procedure.Meshing algorithm allows us to obtain the mesh, which satisfies the Delaunay criteria: it is called \\the possible triangles method". This algorithm is based on the well-known Fortune method of Voronoi diagram constructing for a certain set of points in the plane. The graphical representation of the possible triangles method is shown. It is suitable to use generalization of Delaunay triangulation in order to construct meshes with polygonal cells in case of multiple nodes close to be lying on the same circle.The viscous incompressible fluid flow is described by the Navier | Stokes equations and the mass conservation equation with certain initial and boundary conditions. A fractional steps method, which allows us to avoid non-physical oscillations of the pressure, provides the timestepping procedure. Using the finite element discretization and the Bubnov | Galerkin method allows us to carry out spatial discretization.For form functions calculation of finite element mesh with polygonal cells, \
Directory of Open Access Journals (Sweden)
Guang Wei Meng
2015-01-01
Full Text Available A new method using the enriched element-free Galerkin method (EEFGM to model functionally graded piezoelectric materials (FGPMs with cracks was presented. To improve the solution accuracy, extended terms were introduced into the approximation function of the conventional element-free Galerkin method (EFGM to describe the displacement and electric fields near the crack. Compared with the conventional EFGM, the new approach requires smaller domain to describe the crack-tip singular field. Additionally, the domain of the nodes was not affected by the crack. Therefore, the visibility method and the diffraction method were no longer needed. The mechanical response of FGPM was discussed, when its material parameters changed exponentially in a certain direction. The modified J-integrals for FGPM were deduced, whose results were compared with the results of the conventional EFGM and the analytical solution. Numerical example results illustrated that this method is feasible and precise.
Directory of Open Access Journals (Sweden)
Ping Zhang
2014-01-01
Full Text Available The variational multiscale element free Galerkin method is extended to simulate the Stokes flow problems in a circular cavity as an irregular geometry. The method is combined with Hughes’s variational multiscale formulation and element free Galerkin method; thus it inherits the advantages of variational multiscale and meshless methods. Meanwhile, a simple technique is adopted to impose the essential boundary conditions which makes it easy to solve problems with complex area. Finally, two examples are solved and good results are obtained as compared with solutions of analytical and numerical methods, which demonstrates that the proposed method is an attractive approach for solving incompressible fluid flow problems in terms of accuracy and stability, even for complex irregular boundaries.
On cell entropy inequality for discontinuous Galerkin methods
Jiang, Guangshan; Shu, Chi-Wang
1993-01-01
We prove a cell entropy inequality for a class of high order discontinuous Galerkin finite element methods approximating conservation laws, which implies convergence for the one dimensional scalar convex case.
Parallel, adaptive finite element methods for conservation laws
Biswas, Rupak; Devine, Karen D.; Flaherty, Joseph E.
1994-01-01
We construct parallel finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. A posteriori estimates of spatial errors are obtained by a p-refinement technique using superconvergence at Radau points. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We compare results using different limiting schemes and demonstrate parallel efficiency through computations on an NCUBE/2 hypercube. We also present results using adaptive h- and p-refinement to reduce the computational cost of the method.
An adaptive finite element method for steady and transient problems
International Nuclear Information System (INIS)
Benner, R.E. Jr.; Davis, H.T.; Scriven, L.E.
1987-01-01
Distributing integral error uniformly over variable subdomains, or finite elements, is an attractive criterion by which to subdivide a domain for the Galerkin/finite element method when localized steep gradients and high curvatures are to be resolved. Examples are fluid interfaces, shock fronts and other internal layers, as well as fluid mechanical and other boundary layers, e.g. thin-film states at solid walls. The uniform distribution criterion is developed into an adaptive technique for one-dimensional problems. Nodal positions can be updated simultaneously with nodal values during Newton iteration, but it is usually better to adopt nearly optimal nodal positions during Newton iteration upon nodal values. Three illustrative problems are solved: steady convection with diffusion, gradient theory of fluid wetting on a solid surface and Buckley-Leverett theory of two phase Darcy flow in porous media
Jiang, Lijian
2010-08-01
In this paper, we discuss a numerical multiscale approach for solving wave equations with heterogeneous coefficients. Our interest comes from geophysics applications and we assume that there is no scale separation with respect to spatial variables. To obtain the solution of these multiscale problems on a coarse grid, we compute global fields such that the solution smoothly depends on these fields. We present a Galerkin multiscale finite element method using the global information and provide a convergence analysis when applied to solve the wave equations. We investigate the relation between the smoothness of the global fields and convergence rates of the global Galerkin multiscale finite element method for the wave equations. Numerical examples demonstrate that the use of global information renders better accuracy for wave equations with heterogeneous coefficients than the local multiscale finite element method. © 2010 IMACS.
Steady-state solution of the PTC thermistor problem using a quadratic spline finite element method
Directory of Open Access Journals (Sweden)
Bahadir A. R.
2002-01-01
Full Text Available The problem of heat transfer in a Positive Temperature Coefficient (PTC thermistor, which may form one element of an electric circuit, is solved numerically by a finite element method. The approach used is based on Galerkin finite element using quadratic splines as shape functions. The resulting system of ordinary differential equations is solved by the finite difference method. Comparison is made with numerical and analytical solutions and the accuracy of the computed solutions indicates that the method is well suited for the solution of the PTC thermistor problem.
A spectral-element discontinuous Galerkin lattice Boltzmann method for nearly incompressible flows
Min, Misun; Lee, Taehun
2011-01-01
We present a spectral-element discontinuous Galerkin lattice Boltzmann method for solving nearly incompressible flows. Decoupling the collision step from the streaming step offers numerical stability at high Reynolds numbers. In the streaming step, we employ high-order spectral-element discontinuous Galerkin discretizations using a tensor product basis of one-dimensional Lagrange interpolation polynomials based on Gauss-Lobatto-Legendre grids. Our scheme is cost-effective with a fully diagonal mass matrix, advancing time integration with the fourth-order Runge-Kutta method. We present a consistent treatment for imposing boundary conditions with a numerical flux in the discontinuous Galerkin approach. We show convergence studies for Couette flows and demonstrate two benchmark cases with lid-driven cavity flows for Re = 400-5000 and flows around an impulsively started cylinder for Re = 550-9500. Computational results are compared with those of other theoretical and computational work that used a multigrid method, a vortex method, and a spectral element model.
Solid finite elements through three decades
Venkatesh, DN; Shrinivasa, U
1994-01-01
conventionally, solid finite elements have been looked upon as just generalizations of two-dimensional finite elements. In this article we trace their development starting from the days of their inception. Keeping in tune with our perceptions on developing finite elements, without taking recourse to any extra variational techniques, we discuss a few of the techniques which have been applied to solid finite elements. Finally we critically examine our own work on formulating solid finite elemen...
Bause, Markus; Radu, Florin A; Köcher, Uwe
2017-01-01
Variational time discretization schemes are getting of increasing importance for the accurate numerical approximation of transient phenomena. The applicability and value of mixed finite element methods in space for simulating transport processes have been demonstrated in a wide class of works. We consider a family of continuous Galerkin-Petrov time discretization schemes that is combined with a mixed finite element approximation of the spatial variables. The existence and uniqueness of the semidiscrete approximation and of the fully discrete solution are established. For this, the Banach-Nečas-Babuška theorem is applied in a non-standard way. Error estimates with explicit rates of convergence are proved for the scalar and vector-valued variable. An optimal order estimate in space and time is proved by duality techniques for the scalar variable. The convergence rates are analyzed and illustrated by numerical experiments, also on stochastically perturbed meshes.
An element-free Galerkin (EFG) method for generalized Fisher equations (GFE)
International Nuclear Information System (INIS)
Shi Ting-Yu; Ge Hong-Xia; Cheng Rong-Jun
2013-01-01
A generalized Fisher equation (GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance. The exact mathematical result of the GFE has been widely used in population dynamics and genetics, where it originated. Many researchers have studied the numerical solutions of the GFE, up to now. In this paper, we introduce an element-free Galerkin (EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics. Compared with other numerical methods, the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. In comparison with the traditional method, numerical solutions show that the new method has higher accuracy and better convergence. Several numerical examples are presented to demonstrate the effectiveness of the method
Finite element method for incompressible two-fluid model using a fractional step method
International Nuclear Information System (INIS)
Uchiyama, Tomomi
1997-01-01
This paper presents a finite element method for an incompressible two-fluid model. The solution algorithm is based on the fractional step method, which is frequently used in the finite element calculation for single-phase flows. The calculating domain is divided into quadrilateral elements with four nodes. The Galerkin method is applied to derive the finite element equations. Air-water two-phase flows around a square cylinder are calculated by the finite element method. The calculation demonstrates the close relation between the volumetric fraction of the gas-phase and the vortices shed from the cylinder, which is favorably compared with the existing data. It is also confirmed that the present method allows the calculation with less CPU time than the SMAC finite element method proposed in my previous paper. (author)
Automation of finite element methods
Korelc, Jože
2016-01-01
New finite elements are needed as well in research as in industry environments for the development of virtual prediction techniques. The design and implementation of novel finite elements for specific purposes is a tedious and time consuming task, especially for nonlinear formulations. The automation of this process can help to speed up this process considerably since the generation of the final computer code can be accelerated by order of several magnitudes. This book provides the reader with the required knowledge needed to employ modern automatic tools like AceGen within solid mechanics in a successful way. It covers the range from the theoretical background, algorithmic treatments to many different applications. The book is written for advanced students in the engineering field and for researchers in educational and industrial environments.
Nonlinear, finite deformation, finite element analysis
Nguyen, Nhung; Waas, Anthony M.
2016-06-01
The roles of the consistent Jacobian matrix and the material tangent moduli, which are used in nonlinear incremental finite deformation mechanics problems solved using the finite element method, are emphasized in this paper, and demonstrated using the commercial software ABAQUS standard. In doing so, the necessity for correctly employing user material subroutines to solve nonlinear problems involving large deformation and/or large rotation is clarified. Starting with the rate form of the principle of virtual work, the derivations of the material tangent moduli, the consistent Jacobian matrix, the stress/strain measures, and the objective stress rates are discussed and clarified. The difference between the consistent Jacobian matrix (which, in the ABAQUS UMAT user material subroutine is referred to as DDSDDE) and the material tangent moduli ( C e ) needed for the stress update is pointed out and emphasized in this paper. While the former is derived based on the Jaumann rate of the Kirchhoff stress, the latter is derived using the Jaumann rate of the Cauchy stress. Understanding the difference between these two objective stress rates is crucial for correctly implementing a constitutive model, especially a rate form constitutive relation, and for ensuring fast convergence. Specifically, the implementation requires the stresses to be updated correctly. For this, the strains must be computed directly from the deformation gradient and corresponding strain measure (for a total form model). Alternatively, the material tangent moduli derived from the corresponding Jaumann rate of the Cauchy stress of the constitutive relation (for a rate form model) should be used. Given that this requirement is satisfied, the consistent Jacobian matrix only influences the rate of convergence. Its derivation should be based on the Jaumann rate of the Kirchhoff stress to ensure fast convergence; however, the use of a different objective stress rate may also be possible. The error associated
Planas Badenas, Ramon
2005-01-01
The solution of the stationary Stokes problem through the finite element method using linear elements leads to instabilities in the obtained solution. There are two sources of difficulties to solve numerically this problem. The first one is related to the incompressibility condition and the second one relies in the analysis of the behavior of a bulid with high viscosity. This work is devoted to the analysis, and implementation in a Matlab code, of the Galerkin Least-Squares stabilization of t...
ANSYS duplicate finite-element checker routine
Ortega, R.
1995-01-01
An ANSYS finite-element code routine to check for duplicated elements within the volume of a three-dimensional (3D) finite-element mesh was developed. The routine developed is used for checking floating elements within a mesh, identically duplicated elements, and intersecting elements with a common face. A space shuttle main engine alternate turbopump development high pressure oxidizer turbopump finite-element model check using the developed subroutine is discussed. Finally, recommendations are provided for duplicate element checking of 3D finite-element models.
Energy Technology Data Exchange (ETDEWEB)
Pereira, Luis Carlos Martins
1998-06-15
New Petrov-Galerkin formulations on the finite element methods for convection-diffusion problems with boundary layers are presented. Such formulations are based on a consistent new theory on discontinuous finite element methods. Existence and uniqueness of solutions for these problems in the new finite element spaces are demonstrated. Some numerical experiments shows how the new formulation operate and also their efficacy. (author)
Entropy conservative finite element schemes
Tadmor, E.
1986-01-01
The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.
A finite element approach for predicting nozzle admittances
Sigman, R. K.; Zinn, B. T.
1983-01-01
A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation by using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. The boundary conditions are described for both straight and curved sonic lines. A two dimensional FEM with linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.
Theoretical determination of nozzle admittances using a finite element approach
Sigman, R. K.; Zinn, B. T.
1980-01-01
A finite element method is used to predict the admittances of axisymmetric nozzles. It is assumed that the flow in the nozzle is isentropic and irrotational, and the disturbances are small so that linear analyses apply. An approximate, two dimensional compressible model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear wave equation. The differential form of the acoustic equation is transformed to an integral equation using Galerkin's method, and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. A two-dimensional FEM using linear elements is used to solve the acoustic equation. A one dimensional FEM is also used to solve the reduced equation of Crocco, and the solution verifies the sufficiency of the boundary residual formulation. Comparison between computed admittances and experimental data is shown to be quite good.
Vertical discretization with finite elements for a global hydrostatic model on the cubed sphere
Yi, Tae-Hyeong; Park, Ja-Rin
2017-06-01
A formulation of Galerkin finite element with basis-spline functions on a hybrid sigma-pressure coordinate is presented to discretize the vertical terms of global Eulerian hydrostatic equations employed in a numerical weather prediction system, which is horizontally discretized with high-order spectral elements on a cubed sphere grid. This replaces the vertical discretization of conventional central finite difference that is first-order accurate in non-uniform grids and causes numerical instability in advection-dominant flows. Therefore, a model remains in the framework of Galerkin finite elements for both the horizontal and vertical spatial terms. The basis-spline functions, obtained from the de-Boor algorithm, are employed to derive both the vertical derivative and integral operators, since Eulerian advection terms are involved. These operators are used to discretize the vertical terms of the prognostic and diagnostic equations. To verify the vertical discretization schemes and compare their performance, various two- and three-dimensional idealized cases and a hindcast case with full physics are performed in terms of accuracy and stability. It was shown that the vertical finite element with the cubic basis-spline function is more accurate and stable than that of the vertical finite difference, as indicated by faster residual convergence, fewer statistical errors, and reduction in computational mode. This leads to the general conclusion that the overall performance of a global hydrostatic model might be significantly improved with the vertical finite element.
Mulder, W.A.; Zhebel, E.; Minisini, S.
2013-01-01
We analyse the time-stepping stability for the 3-D acoustic wave equation, discretized on tetrahedral meshes. Two types of methods are considered: mass-lumped continuous finite elements and the symmetric interior-penalty discontinuous Galerkin method. Combining the spatial discretization with the
Shock capturing in discontinuous Galerkin spectral elements via the entropy viscosity method
Hackl, Jason; Shringarpure, Mrugesh; Fischer, Paul; Balachandar, Sivaramakrishnan
2017-11-01
We present a 3D discontinuous Galerkin spectral element solver for compressible flows with shock waves using artificial viscosity to regularize the solution for representation by nested tensor products of high-order Lagrange polynomials. The viscosity is constructed from a smoothed evaluation of the residual of an entropy inequality, localizing the artificial viscosity around shock waves and other flow features that would otherwise not be representable in spectral elements without thermodynamic violations due to Gibbs oscillations. Applied to the Guermond-Popov (2014) stress tensor, this smoothed, continuous artificial viscosity is easily integrated with the non-symmetric numerical fluxes of Baumann and Oden (1999). The method is implemented on top of nek5000, leveraging an outstanding high-performance spectral element code to solve shocked flows over curved surfaces. The interaction of a Mach 3 shock with a sphere is shown to demonstrate this capability. Supported by the Department of Energy Predictive Science Academic Alliance Program Contract DE-NA0002378.
A finite element method for the thermochemical decomposition of polymeric materials. I - Theory
Sullivan, R. M.; Salamon, N. J.
1992-01-01
The governing differential equations are developed to model the thermomechanical behavior of chemically decomposing, polymeric materials. These equations account for thermal and gaseous diffusion through a poroelastic, transversely isotropic solid. The Bubnov-Galerkin finite element method is applied to the governing equations to cast the coupled set into a single matrix equation. A method for solving these equations simultaneously at each time step is discussed.
Development of three-dimensional transport code by the double finite element method
International Nuclear Information System (INIS)
Fujimura, Toichiro
1985-01-01
Development of a three-dimensional neutron transport code by the double finite element method is described. Both of the Galerkin and variational methods are adopted to solve the problem, and then the characteristics of them are compared. Computational results of the collocation method, developed as a technique for the vaviational one, are illustrated in comparison with those of an Ssub(n) code. (author)
Application of a circulation model in bays, using the finite element method
International Nuclear Information System (INIS)
Soares, R.
1984-01-01
The circulation of water was studied in different areas in 'Baia de Sepetiba', in the State of Rio de Janeiro, Brazil. The method applied on the mathematical studies was Galerkin's method and ths originated a system of equations which described all the water motions. The Finite Element method used, had great sensitivity to modifications of input data. Comparison between computed and measured data was made in order to verify the conclusions. (M.A.C.) [pt
Physically-based Surface Texture Synthesis Using a Coupled Finite Element System.
Bajaj, Chandrajit; Zhang, Yongjie; Xu, Guoliang
2008-01-01
This paper describes a stable and robust finite element solver for physically-based texture synthesis over arbitrary manifold surfaces. Our approach solves the reaction-diffusion equation coupled with an anisotropic diffusion equation over surfaces, using a Galerkin based finite element method (FEM). This method avoids distortions and discontinuities often caused by traditional texture mapping techniques, especially for arbitrary manifold surfaces. Several varieties of textures are obtained by selecting different values of control parameters in the governing differential equations, and furthermore enhanced quality textures are generated by fairing out noise in input surface meshes.
Transient finite element modeling of functional electrical stimulation.
Filipovic, Nenad D; Peulic, Aleksandar S; Zdravkovic, Nebojsa D; Grbovic-Markovic, Vesna M; Jurisic-Skevin, Aleksandra J
2011-03-01
Transcutaneous functional electrical stimulation is commonly used for strengthening muscle. However, transient effects during stimulation are not yet well explored. The effect of an amplitude change of the stimulation can be described by static model, but there is no differency for different pulse duration. The aim of this study is to present the finite element (FE) model of a transient electrical stimulation on the forearm. Discrete FE equations were derived by using a standard Galerkin procedure. Different tissue conductive and dielectric properties are fitted using least square method and trial and error analysis from experimental measurement. This study showed that FE modeling of electrical stimulation can give the spatial-temporal distribution of applied current in the forearm. Three different cases were modeled with the same geometry but with different input of the current pulse, in order to fit the tissue properties by using transient FE analysis. All three cases were compared with experimental measurements of intramuscular voltage on one volunteer.
Finite element analysis of the neutron transport equation in spherical geometry
International Nuclear Information System (INIS)
Kim, Yong Ill; Kim, Jong Kyung; Suk, Soo Dong
1992-01-01
The Galerkin formulation of the finite element method is applied to the integral law of the first-order form of the one-group neutron transport equation in one-dimensional spherical geometry. Piecewise linear or quadratic Lagrange polynomials are utilized in the integral law for the angular flux to establish a set of linear algebraic equations. Numerical analyses are performed for the scalar flux distribution in a heterogeneous sphere as well as for the criticality problem in a uniform sphere. For the criticality problems in the uniform sphere, the results of the finite element method, with the use of continuous finite elements in space and angle, are compared with the exact solutions. In the heterogeneous problem, the scalar flux distribution obtained by using discontinuous angular and spatical finite elements is in good agreement with that from the ANISN code calculation. (Author)
Domain decomposition methods for mortar finite elements
Energy Technology Data Exchange (ETDEWEB)
Widlund, O.
1996-12-31
In the last few years, domain decomposition methods, previously developed and tested for standard finite element methods and elliptic problems, have been extended and modified to work for mortar and other nonconforming finite element methods. A survey will be given of work carried out jointly with Yves Achdou, Mario Casarin, Maksymilian Dryja and Yvon Maday. Results on the p- and h-p-version finite elements will also be discussed.
A first course in finite elements
Fish, Jacob
2007-01-01
Developed from the authors, combined total of 50 years undergraduate and graduate teaching experience, this book presents the finite element method formulated as a general-purpose numerical procedure for solving engineering problems governed by partial differential equations. Focusing on the formulation and application of the finite element method through the integration of finite element theory, code development, and software application, the book is both introductory and self-contained, as well as being a hands-on experience for any student. This authoritative text on Finite Elements:Adopts
A blended continuous–discontinuous finite element method for solving the multi-fluid plasma model
Energy Technology Data Exchange (ETDEWEB)
Sousa, E.M., E-mail: sousae@uw.edu; Shumlak, U., E-mail: shumlak@uw.edu
2016-12-01
The multi-fluid plasma model represents electrons, multiple ion species, and multiple neutral species as separate fluids that interact through short-range collisions and long-range electromagnetic fields. The model spans a large range of temporal and spatial scales, which renders the model stiff and presents numerical challenges. To address the large range of timescales, a blended continuous and discontinuous Galerkin method is proposed, where the massive ion and neutral species are modeled using an explicit discontinuous Galerkin method while the electrons and electromagnetic fields are modeled using an implicit continuous Galerkin method. This approach is able to capture large-gradient ion and neutral physics like shock formation, while resolving high-frequency electron dynamics in a computationally efficient manner. The details of the Blended Finite Element Method (BFEM) are presented. The numerical method is benchmarked for accuracy and tested using two-fluid one-dimensional soliton problem and electromagnetic shock problem. The results are compared to conventional finite volume and finite element methods, and demonstrate that the BFEM is particularly effective in resolving physics in stiff problems involving realistic physical parameters, including realistic electron mass and speed of light. The benefit is illustrated by computing a three-fluid plasma application that demonstrates species separation in multi-component plasmas.
Interaction of high-speed compressible viscous flow and structure by adaptive finite element method
International Nuclear Information System (INIS)
Limtrakarn, Wiroj; Dechaumphai, Pramote
2004-01-01
Interaction behaviors of high-speed compressible viscous flow and thermal-structural response of structure are presented. The compressible viscous laminar flow behavior based on the Navier-Stokes equations is predicted by using an adaptive cell-centered finite-element method. The energy equation and the quasi-static structural equations for aerodynamically heated structures are solved by applying the Galerkin finite-element method. The finite-element formulation and computational procedure are described. The performance of the combined method is evaluated by solving Mach 4 flow past a flat plate and comparing with the solution from the finite different method. To demonstrate their interaction, the high-speed flow, structural heat transfer, and deformation phenomena are studied by applying the present method to Mach 10 flow past a flat plate
FEBio: finite elements for biomechanics.
Maas, Steve A; Ellis, Benjamin J; Ateshian, Gerard A; Weiss, Jeffrey A
2012-01-01
In the field of computational biomechanics, investigators have primarily used commercial software that is neither geared toward biological applications nor sufficiently flexible to follow the latest developments in the field. This lack of a tailored software environment has hampered research progress, as well as dissemination of models and results. To address these issues, we developed the FEBio software suite (http://mrl.sci.utah.edu/software/febio), a nonlinear implicit finite element (FE) framework, designed specifically for analysis in computational solid biomechanics. This paper provides an overview of the theoretical basis of FEBio and its main features. FEBio offers modeling scenarios, constitutive models, and boundary conditions, which are relevant to numerous applications in biomechanics. The open-source FEBio software is written in C++, with particular attention to scalar and parallel performance on modern computer architectures. Software verification is a large part of the development and maintenance of FEBio, and to demonstrate the general approach, the description and results of several problems from the FEBio Verification Suite are presented and compared to analytical solutions or results from other established and verified FE codes. An additional simulation is described that illustrates the application of FEBio to a research problem in biomechanics. Together with the pre- and postprocessing software PREVIEW and POSTVIEW, FEBio provides a tailored solution for research and development in computational biomechanics.
Finite element coiled cochlea model
Isailovic, Velibor; Nikolic, Milica; Milosevic, Zarko; Saveljic, Igor; Nikolic, Dalibor; Radovic, Milos; Filipović, Nenad
2015-12-01
Cochlea is important part of the hearing system, and thanks to special structure converts external sound waves into neural impulses which go to the brain. Shape of the cochlea is like snail, so geometry of the cochlea model is complex. The simplified cochlea coiled model was developed using finite element method inside SIFEM FP7 project. Software application is created on the way that user can prescribe set of the parameters for spiral cochlea, as well as material properties and boundary conditions to the model. Several mathematical models were tested. The acoustic wave equation for describing fluid in the cochlea chambers - scala vestibuli and scala timpani, and Newtonian dynamics for describing vibrations of the basilar membrane are used. The mechanical behavior of the coiled cochlea was analyzed and the third chamber, scala media, was not modeled because it does not have a significant impact on the mechanical vibrations of the basilar membrane. The obtained results are in good agreement with experimental measurements. Future work is needed for more realistic geometry model. Coiled model of the cochlea was created and results are compared with initial simplified coiled model of the cochlea.
Theoretical determination of nozzle admittances using a finite element method
Sigman, R. K.; Zinn, B. T.
1979-01-01
A finite element method (FEM) is used to predict the admittances of axisymmetric nozzles. The flow in the nozzle is assumed to be isentropic and the disturbances are assumed to be small so that linear analyses apply. An approximate two dimensional compressible flow model is used to describe the steady flow in the nozzle. The propagation of acoustic disturbances is governed by the complete linear acoustic wave equation. This partial differential wave equation is transformed to an integral equation using Galerkin's method and Green's theorem is applied so that the acoustic boundary conditions can be introduced through the boundary residuals. A two dimensional finite element method using linear triangular elements is used to solve the integral acoustic equation. A one dimensional FEM is used to solve the reduced nozzle acoustic equation developed by Crocco and the solution is used to verify the sufficiency of the boundary residual formation. It is shown that agreement between predicted values of the admittance and experimental data is quite good.
Discontinuous Galerkin Method for Hyperbolic Conservation Laws
Mousikou, Ioanna
2016-11-11
Hyperbolic conservation laws form a special class of partial differential equations. They describe phenomena that involve conserved quantities and their solutions show discontinuities which reflect the formation of shock waves. We consider one-dimensional systems of hyperbolic conservation laws and produce approximations using finite difference, finite volume and finite element methods. Due to stability issues of classical finite element methods for hyperbolic conservation laws, we study the discontinuous Galerkin method, which was recently introduced. The method involves completely discontinuous basis functions across each element and it can be considered as a combination of finite volume and finite element methods. We illustrate the implementation of discontinuous Galerkin method using Legendre polynomials, in case of scalar equations and in case of quasi-linear systems, and we review important theoretical results about stability and convergence of the method. The applications of finite volume and discontinuous Galerkin methods to linear and non-linear scalar equations, as well as to the system of elastodynamics, are exhibited.
Simplified Discontinuous Galerkin Methods for Systems of Conservation Laws with Convex Extension
Barth, Timothy J.
1999-01-01
Simplified forms of the space-time discontinuous Galerkin (DG) and discontinuous Galerkin least-squares (DGLS) finite element method are developed and analyzed. The new formulations exploit simplifying properties of entropy endowed conservation law systems while retaining the favorable energy properties associated with symmetric variable formulations.
Finite-Element Software for Conceptual Design
DEFF Research Database (Denmark)
Lindemann, J.; Sandberg, G.; Damkilde, Lars
2010-01-01
and research. Forcepad is an effort to provide a conceptual design and teaching tool in a finite-element software package. Forcepad is a two-dimensional finite-element application based on the same conceptual model as image editing applications such as Adobe Photoshop or Microsoft Paint. Instead of using...
Non-linear finite element modeling
DEFF Research Database (Denmark)
Mikkelsen, Lars Pilgaard
The note is written for courses in "Non-linear finite element method". The note has been used by the author teaching non-linear finite element modeling at Civil Engineering at Aalborg University, Computational Mechanics at Aalborg University Esbjerg, Structural Engineering at the University...
Nonlinear finite element modeling of corrugated board
A. C. Gilchrist; J. C. Suhling; T. J. Urbanik
1999-01-01
In this research, an investigation on the mechanical behavior of corrugated board has been performed using finite element analysis. Numerical finite element models for corrugated board geometries have been created and executed. Both geometric (large deformation) and material nonlinearities were included in the models. The analyses were performed using the commercial...
Least-squares finite element methods
Bochev, Pavel
2009-01-01
Since their emergence, finite element methods have taken a place as one of the most versatile and powerful methodologies for the approximate numerical solution of Partial Differential Equations. This book presents the theory and practice of least-square finite element methods, their strengths and weaknesses, successes, and open problems
Finite element and finite difference methods in electromagnetic scattering
Morgan, MA
2013-01-01
This second volume in the Progress in Electromagnetic Research series examines recent advances in computational electromagnetics, with emphasis on scattering, as brought about by new formulations and algorithms which use finite element or finite difference techniques. Containing contributions by some of the world's leading experts, the papers thoroughly review and analyze this rapidly evolving area of computational electromagnetics. Covering topics ranging from the new finite-element based formulation for representing time-harmonic vector fields in 3-D inhomogeneous media using two coupled sca
Books and monographs on finite element technology
Noor, A. K.
1985-01-01
The present paper proviees a listing of all of the English books and some of the foreign books on finite element technology, taking into account also a list of the conference proceedings devoted solely to finite elements. The references are divided into categories. Attention is given to fundamentals, mathematical foundations, structural and solid mechanics applications, fluid mechanics applications, other applied science and engineering applications, computer implementation and software systems, computational and modeling aspects, special topics, boundary element methods, proceedings of symmposia and conferences on finite element technology, bibliographies, handbooks, and historical accounts.
Convergence of Crank-Nicolson-Galerkin discrete scheme for ...
African Journals Online (AJOL)
We studied the maximum-norm error estimate for the Galerkin finite element discretization in time of a stochastic wave equation by the Crank-Nicolson time stepping finite difference method. The error estimate was obtained by using the notions of rational function and resolvent estimates.
Finite element analysis of piezoelectric materials
International Nuclear Information System (INIS)
Lowrie, F.; Stewart, M.; Cain, M.; Gee, M.
1999-01-01
This guide is intended to help people wanting to do finite element analysis of piezoelectric materials by answering some of the questions that are peculiar to piezoelectric materials. The document is not intended as a complete beginners guide for finite element analysis in general as this is better dealt with by the individual software producers. The guide is based around the commercial package ANSYS as this is a popular package amongst piezoelectric material users, however much of the information will still be useful to users of other finite element codes. (author)
Simulation of incompressible flows with heat and mass transfer using parallel finite element method
Directory of Open Access Journals (Sweden)
Jalal Abedi
2003-02-01
Full Text Available The stabilized finite element formulations based on the SUPG (Stream-line-Upwind/Petrov-Galerkin and PSPG (Pressure-Stabilization/Petrov-Galerkin methods are developed and applied to solve buoyancy-driven incompressible flows with heat and mass transfer. The SUPG stabilization term allows us to solve flow problems at high speeds (advection dominant flows and the PSPG term eliminates instabilities associated with the use of equal order interpolation functions for both pressure and velocity. The finite element formulations are implemented in parallel using MPI. In parallel computations, the finite element mesh is partitioned into contiguous subdomains using METIS, which are then assigned to individual processors. To ensure a balanced load, the number of elements assigned to each processor is approximately equal. To solve nonlinear systems in large-scale applications, we developed a matrix-free GMRES iterative solver. Here we totally eliminate a need to form any matrices, even at the element levels. To measure the accuracy of the method, we solve 2D and 3D example of natural convection flows at moderate to high Rayleigh numbers.
On higher order pyramidal finite elements
Czech Academy of Sciences Publication Activity Database
Liu, L.; Davies, K.B.; Křížek, Michal; Guan, L.
2011-01-01
Roč. 3, č. 2 (2011), s. 131-140 ISSN 2070-0733 R&D Projects: GA AV ČR(CZ) IAA100190803 Institutional research plan: CEZ:AV0Z10190503 Keywords : pyramidal polynomial basis functions * finite element method * composite elements * three-dimensional mortar elements Subject RIV: BA - General Mathematics Impact factor: 0.750, year: 2011
Finite element methods a practical guide
Whiteley, Jonathan
2017-01-01
This book presents practical applications of the finite element method to general differential equations. The underlying strategy of deriving the finite element solution is introduced using linear ordinary differential equations, thus allowing the basic concepts of the finite element solution to be introduced without being obscured by the additional mathematical detail required when applying this technique to partial differential equations. The author generalizes the presented approach to partial differential equations which include nonlinearities. The book also includes variations of the finite element method such as different classes of meshes and basic functions. Practical application of the theory is emphasised, with development of all concepts leading ultimately to a description of their computational implementation illustrated using Matlab functions. The target audience primarily comprises applied researchers and practitioners in engineering, but the book may also be beneficial for graduate students.
Advanced finite element method in structural engineering
Long, Yu-Qiu; Long, Zhi-Fei
2009-01-01
This book systematically introduces the research work on the Finite Element Method completed over the past 25 years. Original theoretical achievements and their applications in the fields of structural engineering and computational mechanics are discussed.
Multigrid methods for mortar finite elements
Wohlmuth, Barbara
2000-01-01
Multigrid methods for mortar finite elements / R. Krause ; B. Wohlmuth. - In: Multigrid methods VI / Erik Dick ... (ed.). - Berlin u.a. : Springer, 2000. - S. 136-142 (Lecture notes in computational science and engineering ; 14)
FINITE ELEMENT MODEL FOR PREDICTING RESIDUAL ...
African Journals Online (AJOL)
direction (σx) had a maximum value of 375MPa (tensile) and minimum value of ... These results shows that the residual stresses obtained by prediction from the finite element method are in fair agreement with the experimental results.
Discrete mechanics Based on Finite Element Methods
Chen, Jing-bo; Guo, Han-Ying; Wu, Ke
2002-01-01
Discrete Mechanics based on finite element methods is presented in this paper. We also explore the relationship between this discrete mechanics and Veselov discrete mechanics. High order discretizations are constructed in terms of high order interpolations.
ANSYS mechanical APDL for finite element analysis
Thompson, Mary Kathryn
2017-01-01
ANSYS Mechanical APDL for Finite Element Analysis provides a hands-on introduction to engineering analysis using one of the most powerful commercial general purposes finite element programs on the market. Students will find a practical and integrated approach that combines finite element theory with best practices for developing, verifying, validating and interpreting the results of finite element models, while engineering professionals will appreciate the deep insight presented on the program's structure and behavior. Additional topics covered include an introduction to commands, input files, batch processing, and other advanced features in ANSYS. The book is written in a lecture/lab style, and each topic is supported by examples, exercises and suggestions for additional readings in the program documentation. Exercises gradually increase in difficulty and complexity, helping readers quickly gain confidence to independently use the program. This provides a solid foundation on which to build, preparing readers...
Finite elements in CAD and ADINA
International Nuclear Information System (INIS)
Bathe, K.J.
1986-01-01
The use of finite element methods in computer-aided-design - CAD - is discussed. Some current capabilities are presented and important future developments are outlined. The discussion focusses on the use of the ADINA program in CAD applications. (orig.)
Finite element approximation of the Isaacs equation
Salgado, Abner J.; Zhang, Wujun
2015-01-01
We propose and analyze a two-scale finite element method for the Isaacs equation. The fine scale is given by the mesh size $h$ whereas the coarse scale $\\varepsilon$ is dictated by an integro-differential approximation of the partial differential equation. We show that the method satisfies the discrete maximum principle provided that the mesh is weakly acute. This, in conjunction with weak operator consistency of the finite element method, allows us to establish convergence of the numerical s...
Finite Element Model of Gear Induction Hardening
Hodek, J; Zemko, M; Shykula, P
2015-01-01
International audience; This paper presents a finite element model of a gear induction hardening process. The gear was surface-heated by an induction coil and quickly cooled by spraying water. The finite element model was developed as a three-dimensional model. The electromagnetic field, temperature field, stress distribution and microstructure distribution were examined. Temperature and microstructural characteristics were measured and used. The gear material data was obtained in part by mea...
The finite element method in electromagnetics
Jin, Jianming
2014-01-01
A new edition of the leading textbook on the finite element method, incorporating major advancements and further applications in the field of electromagnetics The finite element method (FEM) is a powerful simulation technique used to solve boundary-value problems in a variety of engineering circumstances. It has been widely used for analysis of electromagnetic fields in antennas, radar scattering, RF and microwave engineering, high-speed/high-frequency circuits, wireless communication, electromagnetic compatibility, photonics, remote sensing, biomedical engineering, and space exploration. The
Surgery simulation using fast finite elements
DEFF Research Database (Denmark)
Bro-Nielsen, Morten
1996-01-01
This paper describes our recent work on real-time surgery simulation using fast finite element models of linear elasticity. In addition, we discuss various improvements in terms of speed and realism......This paper describes our recent work on real-time surgery simulation using fast finite element models of linear elasticity. In addition, we discuss various improvements in terms of speed and realism...
Error Estimates for a Semidiscrete Finite Element Method for Fractional Order Parabolic Equations
Jin, Bangti
2013-01-01
We consider the initial boundary value problem for a homogeneous time-fractional diffusion equation with an initial condition ν(x) and a homogeneous Dirichlet boundary condition in a bounded convex polygonal domain Ω. We study two semidiscrete approximation schemes, i.e., the Galerkin finite element method (FEM) and lumped mass Galerkin FEM, using piecewise linear functions. We establish almost optimal with respect to the data regularity error estimates, including the cases of smooth and nonsmooth initial data, i.e., ν ∈ H2(Ω) ∩ H0 1(Ω) and ν ∈ L2(Ω). For the lumped mass method, the optimal L2-norm error estimate is valid only under an additional assumption on the mesh, which in two dimensions is known to be satisfied for symmetric meshes. Finally, we present some numerical results that give insight into the reliability of the theoretical study. © 2013 Society for Industrial and Applied Mathematics.
Quadrature representation of finite element variational forms
DEFF Research Database (Denmark)
Ølgaard, Kristian Breum; Wells, Garth N.
2012-01-01
This chapter addresses the conventional run-time quadrature approach for the numerical integration of local element tensors associated with finite element variational forms, and in particular automated optimizations that can be performed to reduce the number of floating point operations. An alter......This chapter addresses the conventional run-time quadrature approach for the numerical integration of local element tensors associated with finite element variational forms, and in particular automated optimizations that can be performed to reduce the number of floating point operations...
Mirzaee, H.; Ryan, J.K.; Kirby, R.M.
2011-01-01
The discontinuous Galerkin (DG) methods provide a high-order extension of the finite volume method in much the same way as high-order or spectral/hp elements extend standard finite elements. However, lack of inter-element continuity is often contrary to the smoothness assumptions upon which many
A simple finite element method for boundary value problems with a Riemann–Liouville derivative
Jin, Bangti
2016-02-01
© 2015 Elsevier B.V. All rights reserved. We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order α∈(3/2,2) on the unit interval (0,1). The standard Galerkin finite element approximation converges slowly due to the presence of singularity term xα-^{1} in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and ^{L2}(D) error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.
Finite element analysis of tibial fractures
DEFF Research Database (Denmark)
Wong, Christian Nai En; Mikkelsen, Mikkel Peter W; Hansen, Leif Berner
2010-01-01
INTRODUCTION: Fractures of the tibial shaft are relatively common injuries. There are indications that tibial shaft fractures share characteristics in terms of site, type and local fracture mechanisms. In this study, we aimed to set up a mathematical, computer-based model using finite element...... of bony healing. The biomechanical results are the basis for fracture healing, biomechanical fall analysis and stability analysis of osteosynthesis. MATERIAL AND METHODS: A finite element model of the bony part of the lower leg was generated on the basis of computed tomography data from the Visible Human...... Project. The data consisted of 21,219 3D elements with a cortical shell and a trabecular core. Three types of load of torsion, a direct lateral load and axial compression were applied. RESULTS: The finite element linear static analysis resulted in relevant fracture localizations and indicated relevant...
Quadrilateral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Benzley, Steven E
2012-10-16
Techniques for coarsening a quadrilateral mesh are described. These techniques include identifying a coarsening region within the quadrilateral mesh to be coarsened. Quadrilateral elements along a path through the coarsening region are removed. Node pairs along opposite sides of the path are identified. The node pairs along the path are then merged to collapse the path.
A finite element solution of transonic flow
Tatum, K. E.
1978-01-01
The use of finite elements is explored in a field in which its use has previously not been deemed very feasible, that of transonic flow. The specific problem chosen is that of steady small-disturbance transonic flow. The nonlinear equations are formulated with an artificial viscosity term added to yield the proper domain of dependence and directional bias in supersonic regions and across imbedded shock waves. Justification is given for the problem and means of solution chosen, and the potential advantages of the finite element procedure over standard finite difference procedures are discussed. Several possible improvements on the method as presently derived are stated. Computational mesh requirements and certain mesh variations are described. Some results equivalent to finite difference calculations are given as a sample solution.
A strongly conservative finite element method for the coupling of Stokes and Darcy flow
Kanschat, G.
2010-08-01
We consider a model of coupled free and porous media flow governed by Stokes and Darcy equations with the Beavers-Joseph-Saffman interface condition. This model is discretized using divergence-conforming finite elements for the velocities in the whole domain. Discontinuous Galerkin techniques and mixed methods are used in the Stokes and Darcy subdomains, respectively. This discretization is strongly conservative in Hdiv(Ω) and we show convergence. Numerical results validate our findings and indicate optimal convergence orders. © 2010 Elsevier Inc.
Finite element procedures for time-dependent convection-diffusion-reaction systems
Tezduyar, T. E.; Park, Y. J.; Deans, H. A.
1988-01-01
New finite element procedures based on the streamline-upwind/Petrov-Galerkin formulations are developed for time-dependent convection-diffusion-reaction equations. These procedures minimize spurious oscillations for convection-dominated and reaction-dominated problems. The results obtained for representative numerical examples are accurate with minimal oscillations. As a special application problem, the single-well chemical tracer test (a procedure for measuring oil remaining in a depleted field) is simulated numerically. The results show the importance of temperature effects on the interpreted value of residual oil saturation from such tests.
International Nuclear Information System (INIS)
Franca, L.P.; Stenberg, R.
1989-06-01
Stability conditions are described to analyze a variational formulation emanating from a variational principle for linear isotropic elasticity. The variational principle is based on four dependent variables (namely, the strain tensor, augmented stress, pressure and displacement) and is shown to be valid for any compressibility including the incompressible limit. An improved convergence error analysis is established for a Galerkin-least-squares method based upon these four variables. The analysis presented establishes convergence for a wide choice of combinations of finite element interpolations. (author) [pt
Calculation of two-dimensional thermal transients by the method of finite elements
International Nuclear Information System (INIS)
Fontoura Rodrigues, J.L.A. da.
1980-08-01
The unsteady linear heat conduction analysis throught anisotropic and/or heterogeneous matter, in either two-dimensional fields with any kind of geometry or three-dimensional fields with axial symmetry is presented. The boundary conditions and the internal heat generation are supposed time - independent. The solution is obtained by modal analysis employing the finite element method under Galerkin formulation. Optionally, it can be used with a reduced resolution method called Stoker Economizing Method wich allows a decrease on the program processing costs. (Author) [pt
Geometrically unfitted finite element methods and applications
Burman, Erik; Larson, Mats; Olshanskii, Maxim
2017-01-01
This book provides a snapshot of the state of the art of the rapidly evolving field of integration of geometric data in finite element computations. The contributions to this volume, based on research presented at the UCL workshop on the topic in January 2016, include three review papers on core topics such as fictitious domain methods for elasticity, trace finite element methods for partial differential equations defined on surfaces, and Nitsche’s method for contact problems. Five chapters present original research articles on related theoretical topics, including Lagrange multiplier methods, interface problems, bulk-surface coupling, and approximation of partial differential equations on moving domains. Finally, two chapters discuss advanced applications such as crack propagation or flow in fractured poroelastic media. This is the first volume that provides a comprehensive overview of the field of unfitted finite element methods, including recent techniques such as cutFEM, traceFEM, ghost penalty, and aug...
Verification of Orthogrid Finite Element Modeling Techniques
Steeve, B. E.
1996-01-01
The stress analysis of orthogrid structures, specifically with I-beam sections, is regularly performed using finite elements. Various modeling techniques are often used to simplify the modeling process but still adequately capture the actual hardware behavior. The accuracy of such 'Oshort cutso' is sometimes in question. This report compares three modeling techniques to actual test results from a loaded orthogrid panel. The finite element models include a beam, shell, and mixed beam and shell element model. Results show that the shell element model performs the best, but that the simpler beam and beam and shell element models provide reasonable to conservative results for a stress analysis. When deflection and stiffness is critical, it is important to capture the effect of the orthogrid nodes in the model.
Finite Element Modeling of Cracks and Joints
Directory of Open Access Journals (Sweden)
Jozef Čížik
2006-12-01
Full Text Available The application of finite element method to the analysis of discontinuous structural systems has received a considerable interest in recent years. Examples of problems in which discontinuities play a prominent role in the physical behaviour of a system are numerous and include various types of contact problems and layered or jointed systems. This paper gives a state-of-the-art report on the different methods developed to date for the finite element modelling of cracks and joints in discontinuous systems. Particular attention, however, has been given to the use of joint/interface elements, since their application is considered to be most appropriate for modelling of all kinds of discontinuities that may present in a structural system. A chronology of development of the main types of joint elements, including their pertinent characteristics, is also given. Advantages and disadvantages of the individual methods and types of joint elements presented are briefly discussed, together with various applications of interest.
On the reliability of finite element solutions
International Nuclear Information System (INIS)
Prasad, K.S.R.K.
1975-01-01
The extent of reliability of the finite element method for analysis of nuclear reactor structures, and that of reactor vessels in particular and the need for the engineer to guard against the pitfalls that may arise out of both physical and mathematical models have been high-lighted. A systematic way of checking the model to obtain reasonably accurate solutions is presented. Quite often sophisticated elements are suggested for specific design and stress concentration problems. The desirability or otherwise of these elements, their scope and utility vis-a-vis the use of large stack of conventional elements are discussed from the view point of stress analysts. The methods of obtaining a check on the reliability of the finite element solutions either through modelling changes or an extrapolation technique are discussed. (author)
Visualizing higher order finite elements. Final report
Energy Technology Data Exchange (ETDEWEB)
Thompson, David C; Pebay, Philippe Pierre
2005-11-01
This report contains an algorithm for decomposing higher-order finite elements into regions appropriate for isosurfacing and proves the conditions under which the algorithm will terminate. Finite elements are used to create piecewise polynomial approximants to the solution of partial differential equations for which no analytical solution exists. These polynomials represent fields such as pressure, stress, and momentum. In the past, these polynomials have been linear in each parametric coordinate. Each polynomial coefficient must be uniquely determined by a simulation, and these coefficients are called degrees of freedom. When there are not enough degrees of freedom, simulations will typically fail to produce a valid approximation to the solution. Recent work has shown that increasing the number of degrees of freedom by increasing the order of the polynomial approximation (instead of increasing the number of finite elements, each of which has its own set of coefficients) can allow some types of simulations to produce a valid approximation with many fewer degrees of freedom than increasing the number of finite elements alone. However, once the simulation has determined the values of all the coefficients in a higher-order approximant, tools do not exist for visual inspection of the solution. This report focuses on a technique for the visual inspection of higher-order finite element simulation results based on decomposing each finite element into simplicial regions where existing visualization algorithms such as isosurfacing will work. The requirements of the isosurfacing algorithm are enumerated and related to the places where the partial derivatives of the polynomial become zero. The original isosurfacing algorithm is then applied to each of these regions in turn.
Finite elements for analysis and design
Akin, J E; Davenport, J H
1994-01-01
The finite element method (FEM) is an analysis tool for problem-solving used throughout applied mathematics, engineering, and scientific computing. Finite Elements for Analysis and Design provides a thoroughlyrevised and up-to-date account of this important tool and its numerous applications, with added emphasis on basic theory. Numerous worked examples are included to illustrate the material.Key Features* Akin clearly explains the FEM, a numerical analysis tool for problem-solving throughout applied mathematics, engineering and scientific computing* Basic theory has bee
Finite element code development for modeling detonation of HMX composites
Duran, Adam V.; Sundararaghavan, Veera
2017-01-01
In this work, we present a hydrodynamics code for modeling shock and detonation waves in HMX. A stable efficient solution strategy based on a Taylor-Galerkin finite element (FE) discretization was developed to solve the reactive Euler equations. In our code, well calibrated equations of state for the solid unreacted material and gaseous reaction products have been implemented, along with a chemical reaction scheme and a mixing rule to define the properties of partially reacted states. A linear Gruneisen equation of state was employed for the unreacted HMX calibrated from experiments. The JWL form was used to model the EOS of gaseous reaction products. It is assumed that the unreacted explosive and reaction products are in both pressure and temperature equilibrium. The overall specific volume and internal energy was computed using the rule of mixtures. Arrhenius kinetics scheme was integrated to model the chemical reactions. A locally controlled dissipation was introduced that induces a non-oscillatory stabilized scheme for the shock front. The FE model was validated using analytical solutions for SOD shock and ZND strong detonation models. Benchmark problems are presented for geometries in which a single HMX crystal is subjected to a shock condition.
A finite-elements method for turbulent flow analysis
International Nuclear Information System (INIS)
Autret, A.
1986-03-01
The work discussed here covers turbulent flow calculations using GALERKIN's finite-element method. In our specific case, we have to deal with monophasic incompressible flow in Boussinesq approximation in the normal operating conditions of a primary circuit of nuclear power plant. Turbulence effects on the mean field are taken into account by the k-epsilon model with two evolution equations: one for the kinetic energy of the turbulence, and one for the energy dissipation rate. The wall zone is covered by wall laws, and by REICHARDT's law in particular. A Law is advanced for the epsilon input profile, and a numerical solution is proposed for the physically aberrant values of k and epsilon generated by the model. Single-equation models are reviewed comparatively with the k-epsilon model. A comparison between calculated and analytical solutions or calculated and experimental results is presented for decreasing turbulence behind a grid, for the flow between parallel flat plates with three REYNOLDS numbers, and for backward facing step [fr
Quadrilateral/hexahedral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Scott, Michael A; Benzley, Steven E
2012-10-16
A technique for coarsening a finite element mesh ("FEM") is described. This technique includes identifying a coarsening region within the FEM to be coarsened. Perimeter chords running along perimeter boundaries of the coarsening region are identified. The perimeter chords are redirected to create an adaptive chord separating the coarsening region from a remainder of the FEM. The adaptive chord runs through mesh elements residing along the perimeter boundaries of the coarsening region. The adaptive chord is then extracted to coarsen the FEM.
Slave finite elements: The temporal element approach to nonlinear analysis
Gellin, S.
1984-01-01
A formulation method for finite elements in space and time incorporating nonlinear geometric and material behavior is presented. The method uses interpolation polynomials for approximating the behavior of various quantities over the element domain, and only explicit integration over space and time. While applications are general, the plate and shell elements that are currently being programmed are appropriate to model turbine blades, vanes, and combustor liners.
Finite element analysis of photonic crystal fibers
Uranus, H.P.; Hoekstra, Hugo; van Groesen, Embrecht W.C.
2005-01-01
A finite-element-based vectorial optical mode solver, furnished with Bayliss-Gunzburger-Turkel-like transparent boundary conditions, is used to rigorously analyze photonic crystal fibers (PCFs). Both the real and imaginary part of the modal indices can be computed in a relatively small computational
Finite element simulation of asphalt fatigue testing
DEFF Research Database (Denmark)
Ullidtz, Per; Kieler, Thomas Lau; Kargo, Anders
1997-01-01
damage mechanics.The paper describes how continuum damage mechanics may be used with a finite element program to explain the progressive deterioration of asphalt mixes under laboratory fatigue testing. Both constant stress and constant strain testing are simulated, and compared to the actual results from...
Finite element modelling of solidification phenomena
Indian Academy of Sciences (India)
The process of solidification process is complex in nature and the simulation of such process is required in industry before it is actually undertaken. Finite element method is used to simulate the heat transfer process accompanying the solidification process. The metal and the mould along with the air gap formation is ...
Fast finite elements for surgery simulation
DEFF Research Database (Denmark)
Bro-Nielsen, Morten
1997-01-01
This paper discusses volumetric deformable models for modeling human body parts and organs in surgery simulation systems. These models are built using finite element models for linear elastic materials. To achieve real-time response condensation has been applied to the system stiffness matrix...
Equivalent drawbead model in finite element simulations
Carleer, Bart D.; Carleer, B.D.; Meinders, Vincent T.; Huetink, Han; Lee, J.K.; Kinzel, G.L.; Wagoner, R.
1996-01-01
In 3D simulations of the deep drawing process the drawbead geometries are seldom included. Therefore equivalent drawbeads are used. In order to investigate the drawbead behaviour a 2D plane strain finite element model was used. For verification of this model experiments were performed. The analyses
Simplicial Finite Elements in Higher Dimensions
Czech Academy of Sciences Publication Activity Database
Brandts, J.; Korotov, S.; Křížek, Michal
2007-01-01
Roč. 52, č. 3 (2007), s. 251-265 ISSN 0862-7940 R&D Projects: GA ČR GA201/04/1503 Institutional research plan: CEZ:AV0Z10190503 Keywords : n-simplex * finite element method * superconvergence Subject RIV: BA - General Mathematics
Kawahara, Mutsuto
2016-01-01
This book focuses on the finite element method in fluid flows. It is targeted at researchers, from those just starting out up to practitioners with some experience. Part I is devoted to the beginners who are already familiar with elementary calculus. Precise concepts of the finite element method remitted in the field of analysis of fluid flow are stated, starting with spring structures, which are most suitable to show the concepts of superposition/assembling. Pipeline system and potential flow sections show the linear problem. The advection–diffusion section presents the time-dependent problem; mixed interpolation is explained using creeping flows, and elementary computer programs by FORTRAN are included. Part II provides information on recent computational methods and their applications to practical problems. Theories of Streamline-Upwind/Petrov–Galerkin (SUPG) formulation, characteristic formulation, and Arbitrary Lagrangian–Eulerian (ALE) formulation and others are presented with practical results so...
Hall, Eric Joseph
2016-12-08
We derive computable error estimates for finite element approximations of linear elliptic partial differential equations with rough stochastic coefficients. In this setting, the exact solutions contain high frequency content that standard a posteriori error estimates fail to capture. We propose goal-oriented estimates, based on local error indicators, for the pathwise Galerkin and expected quadrature errors committed in standard, continuous, piecewise linear finite element approximations. Derived using easily validated assumptions, these novel estimates can be computed at a relatively low cost and have applications to subsurface flow problems in geophysics where the conductivities are assumed to have lognormal distributions with low regularity. Our theory is supported by numerical experiments on test problems in one and two dimensions.
Variational Multiscale Finite Element Method for Flows in Highly Porous Media
Iliev, O.
2011-10-01
We present a two-scale finite element method (FEM) for solving Brinkman\\'s and Darcy\\'s equations. These systems of equations model fluid flows in highly porous and porous media, respectively. The method uses a recently proposed discontinuous Galerkin FEM for Stokes\\' equations by Wang and Ye and the concept of subgrid approximation developed by Arbogast for Darcy\\'s equations. In order to reduce the "resonance error" and to ensure convergence to the global fine solution, the algorithm is put in the framework of alternating Schwarz iterations using subdomains around the coarse-grid boundaries. The discussed algorithms are implemented using the Deal.II finite element library and are tested on a number of model problems. © 2011 Society for Industrial and Applied Mathematics.
A singular finite element technique for calculating continuum damping of Alfvén eigenmodes
International Nuclear Information System (INIS)
Bowden, G. W.; Hole, M. J.
2015-01-01
Damping due to continuum resonances can be calculated using dissipation-less ideal magnetohydrodynamics provided that the poles due to these resonances are properly treated. We describe a singular finite element technique for calculating the continuum damping of Alfvén waves. A Frobenius expansion is used to determine appropriate finite element basis functions on an inner region surrounding a pole due to the continuum resonance. The location of the pole due to the continuum resonance and mode frequency is calculated iteratively using a Galerkin method. This method is used to find the complex frequency and mode structure of a toroidicity-induced Alfvén eigenmode in a large aspect ratio circular tokamak and is shown to agree closely with a complex contour technique
Error Analysis of a Finite Element Method for the Space-Fractional Parabolic Equation
Jin, Bangti
2014-01-01
© 2014 Society for Industrial and Applied Mathematics We consider an initial boundary value problem for a one-dimensional fractional-order parabolic equation with a space fractional derivative of Riemann-Liouville type and order α ∈ (1, 2). We study a spatial semidiscrete scheme using the standard Galerkin finite element method with piecewise linear finite elements, as well as fully discrete schemes based on the backward Euler method and the Crank-Nicolson method. Error estimates in the L2(D)- and Hα/2 (D)-norm are derived for the semidiscrete scheme and in the L2(D)-norm for the fully discrete schemes. These estimates cover both smooth and nonsmooth initial data and are expressed directly in terms of the smoothness of the initial data. Extensive numerical results are presented to illustrate the theoretical results.
Finite element fluid modeling of axisymmetric magnetized boundary plasma with recycling neutrals
International Nuclear Information System (INIS)
Zanino, R.
1992-01-01
Finite elements should provide a natural and flexible method for fluid modeling of the tokamak SOL, in particular when the SOL geometry is complex, and/or the poloidal magnetic field is very inclined to the limiter/divertor target. Here we present a Galerkin finite element code, FELS, for transport modeling of a 2-fluid magnetized boundary plasma in an axisymmetry domain, in the presence of recycling neutrals. The classical collisional plasma dynamics along magnetic field lines is taken into account, and a simple diffusive Ansatz is used for the fluxes across magnetic surfaces; electric currents and diamagnetic flows are neglected for the time being. An analytical fluid model is used for the recycling neutrals. Results are shown and discussed for the case of a simple geometry. (orig.)
Introduction to finite and spectral element methods using Matlab
Pozrikidis, Constantine
2014-01-01
The Finite Element Method in One Dimension. Further Applications in One Dimension. High-Order and Spectral Elements in One Dimension. The Finite Element Method in Two Dimensions. Quadratic and Spectral Elements in Two Dimensions. Applications in Mechanics. Viscous Flow. Finite and Spectral Element Methods in Three Dimensions. Appendices. References. Index.
Dispersion analysis of the Pn -Pn-1DG mixed finite element pair for atmospheric modelling
Melvin, Thomas
2018-02-01
Mixed finite element methods provide a generalisation of staggered grid finite difference methods with a framework to extend the method to high orders. The ability to generate a high order method is appealing for applications on the kind of quasi-uniform grids that are popular for atmospheric modelling, so that the method retains an acceptable level of accuracy even around special points in the grid. The dispersion properties of such schemes are important to study as they provide insight into the numerical adjustment to imbalance that is an important component in atmospheric modelling. This paper extends the recent analysis of the P2 - P1DG pair, that is a quadratic continuous and linear discontinuous finite element pair, to higher polynomial orders and also spectral element type pairs. In common with the previously studied element pair, and also with other schemes such as the spectral element and discontinuous Galerkin methods, increasing the polynomial order is found to provide a more accurate dispersion relation for the well resolved part of the spectrum but at the cost of a number of unphysical spectral gaps. The effects of these spectral gaps are investigated and shown to have a varying impact depending upon the width of the gap. Finally, the tensor product nature of the finite element spaces is exploited to extend the dispersion analysis into two-dimensions.
Hybridized Multiscale Discontinuous Galerkin Methods for Multiphysics
2015-09-14
Helmholtz equation for a very wide range of wave frequencies, 4) presented an embedded discontinuous Galerkin (EDG) method for numerically ...of freedom (for the same mesh and polynomial degree of the approximation) in comparison to the well-established finite element methods for diffusion... equations , J. Comput. Phys., 228 (2009), pp. 8841–8855. [4] N.-C. NGUYEN, J. PERAIRE, AND B. COCKBURN, A comparison of HDG methods for Stokes flow, J.
Directory of Open Access Journals (Sweden)
Muhammad Azam
2013-01-01
Full Text Available We introduce improved element-free Galerkin method based on block pulse wavelet integration for numerical approximations to the solution of a system of fourth-order boundary-value problems associated with obstacle, unilateral, and contact problems. Moving least squares (MLS approach is used to construct shape functions with optimized weight functions and basis. Numerical results for test problems are presented in this article to elaborate the pertinent features for the proposed technique. Comparison with existing techniques shows that our proposed method based on integration technique provides better approximation at reduced computational cost.
Clearance gap flow: Simulations by discontinuous Galerkin method and experiments
Czech Academy of Sciences Publication Activity Database
Hála, Jindřich; Luxa, Martin; Bublík, O.; Prausová, H.; Vimmr, J.
2016-01-01
Roč. 92, May (2016), 02073-02073 ISSN 2100-014X. [EFM14 – Experimental Fluid Mechanics 2014. Český Krumlov, 18.11.2014-21.11.2014] Institutional support: RVO:61388998 Keywords : compressible fluid flow * narrow channel flow * discontinuous Galerkin finite element method Subject RIV: BK - Fluid Dynamics
Local discontinuous Galerkin methods for phase transition problems
Tian, Lulu
2015-01-01
In this thesis we develop a local discontinuous Galerkin (LDG) finite element method to solve mathematical models for phase transitions in solids and fluids. The first model we study is called a viscosity-capillarity (VC) system associated with phase transitions in elastic bars and Van der Waals
FINITE ELEMENT ANALYSIS OF ELEMENT ANALYSIS OF A FREE ...
African Journals Online (AJOL)
eobe
formulated as functional minimization. Finite Element Method (FEM) is regarde accurate and versatile numerical too differential equations that model phys. The methodology is used in vari engineering in which the problems ar partial differential equations. The met considerable application in structural e related disciplines.
Gao, Kai
2015-04-14
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both boundaries and the interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
International Nuclear Information System (INIS)
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-01-01
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system
Finite element simulations with ANSYS workbench 16
Lee , Huei-Huang
2015-01-01
Finite Element Simulations with ANSYS Workbench 16 is a comprehensive and easy to understand workbook. It utilizes step-by-step instructions to help guide readers to learn finite element simulations. Twenty seven real world case studies are used throughout the book. Many of these cases are industrial or research projects the reader builds from scratch. All the files readers may need if they have trouble are available for download on the publishers website. Companion videos that demonstrate exactly how to preform each tutorial are available to readers by redeeming the access code that comes in the book. Relevant background knowledge is reviewed whenever necessary. To be efficient, the review is conceptual rather than mathematical. Key concepts are inserted whenever appropriate and summarized at the end of each chapter. Additional exercises or extension research problems are provided as homework at the end of each chapter. A learning approach emphasizing hands-on experiences spreads through this entire book. A...
A finite element method for neutron transport
International Nuclear Information System (INIS)
Ackroyd, R.T.
1983-01-01
A completely boundary-free maximum principle for the first-order Boltzmann equation is derived from the completely boundary-free maximum principle for the mixed-parity Boltzmann equation. When continuity is imposed on the trial function for directions crossing interfaces the completely boundary-free principle for the first-order Boltzmann equation reduces to a maximum principle previously established directly from first principles and indirectly by the Euler-Lagrange method. Present finite element methods for the first-order Boltzmann equation are based on a weighted-residual method which permits the use of discontinuous trial functions. The new principle for the first-order equation can be used as a basis for finite-element methods with the same freedom from boundary conditions as those based on the weighted-residual method. The extremum principle as the parent of the variationally-derived weighted-residual equations ensures their good behaviour. (author)
Finite element reliability analysis of fatigue life
International Nuclear Information System (INIS)
Harkness, H.H.; Belytschko, T.; Liu, W.K.
1992-01-01
Fatigue reliability is addressed by the first-order reliability method combined with a finite element method. Two-dimensional finite element models of components with cracks in mode I are considered with crack growth treated by the Paris law. Probability density functions of the variables affecting fatigue are proposed to reflect a setting where nondestructive evaluation is used, and the Rosenblatt transformation is employed to treat non-Gaussian random variables. Comparisons of the first-order reliability results and Monte Carlo simulations suggest that the accuracy of the first-order reliability method is quite good in this setting. Results show that the upper portion of the initial crack length probability density function is crucial to reliability, which suggests that if nondestructive evaluation is used, the probability of detection curve plays a key role in reliability. (orig.)
Finite element analysis of human joints
International Nuclear Information System (INIS)
Bossart, P.L.; Hollerbach, K.
1996-09-01
Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described
Finite element analysis of tibial fractures
DEFF Research Database (Denmark)
Wong, Christian Nai En; Mikkelsen, Mikkel Peter W; Hansen, Leif Berner
2010-01-01
INTRODUCTION: Fractures of the tibial shaft are relatively common injuries. There are indications that tibial shaft fractures share characteristics in terms of site, type and local fracture mechanisms. In this study, we aimed to set up a mathematical, computer-based model using finite element...... analysis of the bones of the lower leg to examine if such a model is adequate for prediction of fracture locations and patterns. In future studies, we aim to use these biomechanical results to examine fracture prevention, among others, and to simulate different types of osteosynthesis and the process...... of bony healing. The biomechanical results are the basis for fracture healing, biomechanical fall analysis and stability analysis of osteosynthesis. MATERIAL AND METHODS: A finite element model of the bony part of the lower leg was generated on the basis of computed tomography data from the Visible Human...
Finite element analysis of human joints
Energy Technology Data Exchange (ETDEWEB)
Bossart, P.L.; Hollerbach, K.
1996-09-01
Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described.
Finite element modelling of TRIP steels
Energy Technology Data Exchange (ETDEWEB)
Papatriantafillou, I.; Aravas, N.; Haidemenopoulos, G.N. [Dept. of Mechanical and Industrial Engineering, Univ. of Thessaly, Volos (Greece)
2004-11-01
A constitutive model that describes the mechanical behaviour of steels exhibiting ''Transformation Induced Plasticity'' (TRIP) during martensitic transformation is presented. Multiphase TRIP steels are considered as composite materials with a ferritic matrix containing bainite and retained austenite, which gradually transforms into martensite. The effective properties and overall behaviour of TRIP steels are determined by using homogenization techniques for non-linear composites. The developed constitutive model considers the different hardening behaviour of the individual phases and estimates the apportionment of plastic strain and stress between the individual phases of the composite. A methodology for the numerical integration of the resulting elastoplastic constitutive equations in the context of the finite element method is developed and the constitutive model is implemented in a general-purpose finite element program. The prediction of the model in uniaxial tension agrees well with the experimental data. The problem of necking of a bar in uniaxial tension is studied in detail. (orig.)
Finite Element Analysis of Honeycomb Impact Attenuator
Yang, Seung-Yong; Choi, Seung-Kyu; Kim, Nohyu
To participate in Student Formula Society of Automotive Engineers (SAE) competitions, it is necessary to build an impact attenuator that would give an average deceleration not to exceed 20g when it runs into a rigid wall. Students can use numerical simulations or experimental test data to show that their car satisfies this safety requirement. A student group to study formula cars at the Korea University of Technology and Education has designed a vehicle to take part in a SAE competition, and a honeycomb structure was adopted as the impact attenuator. In this paper, finite element calculations were carried out to investigate the dynamic behavior of the honeycomb attenuator. Deceleration and deformation behaviors were studied. Effect of the yield strength was checked by comparing the numerical results. ABAQUS/Explicit finite element code was used.
FINITE ELEMENT ANALYSIS FOR PERIFLEX COUPLINGS
Directory of Open Access Journals (Sweden)
URDEA Mihaela
2015-06-01
Full Text Available The Periflex shaft couplings with rubber sleeve have a hig elasticity and link two shafts in diesel-engine and electric drives. They are simple from the point of view of construction, easily mounted and dismounted. The main goal of this paper is to present a finite element analysis for the Periflex coupling using the Generative Structural Analysis from CATIA software package. This paper presents important information about how to prepare an assembly for creating a static analysis case and also the important steps for developing a finite element analysis. It is very important that the analysis model should have the same behavior as the real, also the loading model. The results are images corresponding to Von Mises Stresses and Translational Displacement magnitude.
Introduction to nonlinear finite element analysis
Kim, Nam-Ho
2015-01-01
This book introduces the key concepts of nonlinear finite element analysis procedures. The book explains the fundamental theories of the field and provides instructions on how to apply the concepts to solving practical engineering problems. Instead of covering many nonlinear problems, the book focuses on three representative problems: nonlinear elasticity, elastoplasticity, and contact problems. The book is written independent of any particular software, but tutorials and examples using four commercial programs are included as appendices: ANSYS, NASTRAN, ABAQUS, and MATLAB. In particular, the MATLAB program includes all source codes so that students can develop their own material models, or different algorithms. This book also: · Presents clear explanations of nonlinear finite element analysis for elasticity, elastoplasticity, and contact problems · Includes many informative examples of nonlinear analyses so that students can clearly understand the nonlinear theory · ...
Finite element simulation of heat transfer
Bergheau, Jean-Michel
2010-01-01
This book introduces the finite element method applied to the resolution of industrial heat transfer problems. Starting from steady conduction, the method is gradually extended to transient regimes, to traditional non-linearities, and to convective phenomena. Coupled problems involving heat transfer are then presented. Three types of couplings are discussed: coupling through boundary conditions (such as radiative heat transfer in cavities), addition of state variables (such as metallurgical phase change), and coupling through partial differential equations (such as electrical phenomena).? A re
Finite element methods for incompressible flow problems
John, Volker
2016-01-01
This book explores finite element methods for incompressible flow problems: Stokes equations, stationary Navier-Stokes equations, and time-dependent Navier-Stokes equations. It focuses on numerical analysis, but also discusses the practical use of these methods and includes numerical illustrations. It also provides a comprehensive overview of analytical results for turbulence models. The proofs are presented step by step, allowing readers to more easily understand the analytical techniques.
Coupled finite element modeling of piezothermoelastic materials
Senousy, M. S.; Rajapakse, R. K. N. D.; Gadala, M.
2007-04-01
The governing equations of piezo-thermoelastic materials show full coupling between mechanical, electric, and temperature fields. It is often assumed in the literature that in high-frequency oscillations, the coupling between the temperature and mechanical displacement and electric field is small and, therefore, can be neglected. A solution for the temperature field is then determined from an uncoupled equation. A finite element (FE) model that accounts for full coupling between the mechanical, electric, and thermal fields, nonlinear constitutive behavior and heat generation resulting from dielectric losses under alternating driving fields is under development. This paper presents a linear fully coupled model as an early development of the fully coupled nonlinear FE model. In the linear model, a solution for all field variables is obtained simultaneously and compared with the uncoupled solution. The finite element model is based on the weighted-residual principle and uses 2-D four-node isoparametric finite elements with four degrees of freedom per node. A thin piezoelectric square disk is modeled to obtain some preliminary understanding of the coupled fields in a piezoelectric stack actuator.
Finite element modeling of lipid bilayer membranes
Feng, Feng; Klug, William S.
2006-12-01
A numerical simulation framework is presented for the study of biological membranes composed of lipid bilayers based on the finite element method. The classic model for these membranes employs a two-dimensional-fluid-like elastic constitutive law which is sensitive to curvature, and subjects vesicles to physically imposed constraints on surface area and volume. This model is implemented numerically via the use of C1-conforming triangular Loop subdivision finite elements. The validity of the framework is tested by computing equilibrium shapes from previously-determined axisymmetric shape-phase diagram of lipid bilayer vesicles with homogeneous material properties. Some of the benefits and challenges of finite element modeling of lipid bilayer systems are discussed, and it is indicated how this framework is natural for future investigation of biologically realistic bilayer structures involving nonaxisymmetric geometries, binding and adhesive interactions, heterogeneous mechanical properties, cytoskeletal interactions, and complex loading arrangements. These biologically relevant features have important consequences for the shape mechanics of nonidealized vesicles and cells, and their study requires not simply advances in theory, but also advances in numerical simulation techniques, such as those presented here.
Finite element modelling of contracting skeletal muscle.
Oomens, C W J; Maenhout, M; van Oijen, C H; Drost, M R; Baaijens, F P
2003-09-29
To describe the mechanical behaviour of biological tissues and transport processes in biological tissues, conservation laws such as conservation of mass, momentum and energy play a central role. Mathematically these are cast into the form of partial differential equations. Because of nonlinear material behaviour, inhomogeneous properties and usually a complex geometry, it is impossible to find closed-form analytical solutions for these sets of equations. The objective of the finite element method is to find approximate solutions for these problems. The concepts of the finite element method are explained on a finite element continuum model of skeletal muscle. In this case, the momentum equations have to be solved with an extra constraint, because the material behaves as nearly incompressible. The material behaviour consists of a highly nonlinear passive part and an active part. The latter is described with a two-state Huxley model. This means that an extra nonlinear partial differential equation has to be solved. The problems and solutions involved with this procedure are explained. The model is used to describe the mechanical behaviour of a tibialis anterior of a rat. The results have been compared with experimentally determined strains at the surface of the muscle. Qualitatively there is good agreement between measured and calculated strains, but the measured strains were higher.
BERSAFE: (BERkeley Structural Analysis by Finite Elements)
International Nuclear Information System (INIS)
Anon.
1991-01-01
BERSAFE is a well-known finite element system which has been under continuous use and development for over 20 years. The BERSAFE system comprises an inter-compatible set of program modules covering static stress analysis, linear dynamics and thermal analysis. Data generation and results presentation modules are also available, along with special supporting functions including automatic crack growth through a model with adaptive meshing. The functionality of BERSAFE, is nowadays very advanced, both in engineering scope and finite element technology. It has seen many firsts, including the front solution and Virtual Crack Extension methods (VCE). More recent additions which have developed out of the Power Industry's requirements are a finite element computational fluid dynamics code, FEAT, and engineering design assessment procedures. These procedures include R6 and R5 for the assessment of the integrity of structures containing defects below and within the creep regime. To use all this software in a user-friendly manner, a new computational environment has been developed, called 'The Harness' which takes advantage of modern hardware and software philosophies. This provides the tool-kit to undertake complete problems, covering determination of fluid loads, structural analysis and failure assessment. In the following sections we describe briefly various components of the BERSAFE suite. (author)
Finite element analysis of coupled electromechanical problems
International Nuclear Information System (INIS)
Melgoza-Vazquez, E.
2001-01-01
The modeling of electromechanical problems is discussed. The simultaneous consideration of two distinct phenomena is required, as the evolution of the electromagnetic and the mechanical parts are influenced by each other. In this work the equations of the coupled problem are described and possible methods of solution are considered. Three general approaches with varying degrees of detail are considered. In the first, a lumped parameter model of the device is constructed from the finite element solution of the electromagnetic problem. A second approach links the electromagnetic field directly with the lumped mechanical part. Lastly, both the electromagnetic and the mechanical systems are considered to be distributed, with the individual domains solved by using the finite element method. In the process of solution of transient problems the need to solve differential-algebraic systems of equations arises and some approaches are presented. It is shown that traditional finite difference formulas may be applied as long as the discretization is made at the element level. Higher order methods and step adaptation are discussed. (author)
Stenroos, Matti; Haueisen, Jens
2008-09-01
In electrocardiographic imaging, epicardial potentials are reconstructed computationally from electrocardiographic measurements. The reconstruction is typically done with the help of the boundary element method (BEM), using the point collocation weighting and constant or linear basis functions. In this paper, we evaluated the performance of constant and linear point collocation and Galerkin BEMs in the epicardial potential problem. The integral equations and discretizations were formulated in terms of the single- and double-layer operators. All inner element integrals were calculated analytically. The computational methods were validated against analytical solutions in a simplified geometry. On the basis of the validation, no method was optimal in all testing scenarios. In the forward computation of the epicardial potential, the linear Galerkin (LG) method produced the smallest errors. The LG method also produced the smallest discretization error on the epicardial surface. In the inverse computation of epicardial potential, the electrode-specific transfer matrix performed better than the full transfer matrix. The Tikhonov 2 regularization outperformed the Tikhonov 0. In the optimal modeling conditions, the best BEM technique depended on electrode positions and chosen error measure. When large modeling errors such as omission of the lungs were present, the choice of the basis and weighting functions was not significant.
Transport and dispersion of pollutants in surface impoundments: a finite element model
International Nuclear Information System (INIS)
Yeh, G.T.
1980-07-01
A surface impoundment model in finite element (SIMFE) is presented to enable the simulation of flow circulations and pollutant transport and dispersion in natural or artificial lakes, reservoirs or ponds with any number of islands. This surface impoundment model consists of two sub-models: hydrodynamic and pollutant transport models. Both submodels are simulated by the finite element method. While the hydrodynamic model is solved by the standard Galerkin finite element scheme, the pollutant transport model can be solved by any of the twelve optional finite element schemes built in the program. Theoretical approximations and the numerical algorithm of SIMFE are described. Detail instruction of the application are given and listing of FORTRAN IV source program are provided. Two sample problems are given. One is for an idealized system with a known solution to show the accuracy and partial validation of the models. The other is applied to Prairie Island for a set of hypothetical input data, typifying a class of problems to which SIMFE may be applied
Transport and dispersion of pollutants in surface impoundments: a finite element model
Energy Technology Data Exchange (ETDEWEB)
Yeh, G.T.
1980-07-01
A surface impoundment model in finite element (SIMFE) is presented to enable the simulation of flow circulations and pollutant transport and dispersion in natural or artificial lakes, reservoirs or ponds with any number of islands. This surface impoundment model consists of two sub-models: hydrodynamic and pollutant transport models. Both submodels are simulated by the finite element method. While the hydrodynamic model is solved by the standard Galerkin finite element scheme, the pollutant transport model can be solved by any of the twelve optional finite element schemes built in the program. Theoretical approximations and the numerical algorithm of SIMFE are described. Detail instruction of the application are given and listing of FORTRAN IV source program are provided. Two sample problems are given. One is for an idealized system with a known solution to show the accuracy and partial validation of the models. The other is applied to Prairie Island for a set of hypothetical input data, typifying a class of problems to which SIMFE may be applied.
A Runge-Kutta discontinuous finite element method for high speed flows
Bey, Kim S.; Oden, J. T.
1991-01-01
A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.
Finite element modeling for materials engineers using Matlab
Oluwole, Oluleke
2014-01-01
Finite Element Modeling for Materials Engineers Using MATLAB® combines the finite element method with MATLAB to offer materials engineers a fast and code-free way of modeling for many materials processes.
Finite Element Based Design and Optimization for Piezoelectric Accelerometers
DEFF Research Database (Denmark)
Liu, Bin; Kriegbaum, B.; Yao, Q.
1998-01-01
A systematic Finite Element design and optimisation procedure is implemented for the development of piezoelectric accelerometers. Most of the specifications of accelerometers can be obtained using the Finite Element simulations. The deviations between the simulated and calibrated sensitivities...
International Nuclear Information System (INIS)
Runca, E.; Melli, P.; Sardei, F.
1985-01-01
A finite-difference scheme and a Galerkin scheme are compared with respect to a very accurate solution describing time-dependent advection and diffusion of air pollutants from a line source in an atmosphere vertically stratified and limited by an inversion layer. The accurate solution was achieved by applying the finite-difference scheme on a very refined grid with a very small time step. The grid size and time step were defined according to stability and accuracy criteria discussed in the text. It is found that for the problem considered the two methods can be considered equally accurate. However, the Galerkin method gives a better approximation in the vicinity of the source. This was assumed to be partly due to the different way the source term is taken into account in the two methods. Improvement of the accuracy of the finite-difference scheme was achieved by approximating, at every step, the contribution of the source term by a Gaussian puff moving and diffusing with the velocity and diffusivity of the source location, instead of utilizing a stepwise function for the numerical approximation of the delta function representing the source term
Shock capturing techniques for hphp-adaptive finite elements
Hierro, A.; Kůs, P. (Pavel); Badia, S.
2016-01-01
The aim of this work is to propose an hp-adaptive algorithm for discontinuous Galerkin methods that is capable to detect the discontinuities and sharp layers and avoid the spurious oscillation of the solution around them. In order to control the spurious oscillations, artificial viscosity is used with the particularity that it is only applied around the layers where the solution changes abruptly. To do so, a novel troubled-cell detector has been developed in order to mark the elements around ...
Error-controlled adaptive finite elements in solid mechanics
National Research Council Canada - National Science Library
Stein, Erwin; Ramm, E
2003-01-01
... . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Error-controlled Adaptive Finite-element-methods . . . . . . . . . . . . Missing Features and Properties of Today's General Purpose FE Programs for Structural...
Cojocaru, E.
2009-01-01
The finite element method has become a preeminent simulation technique in electromagnetics. For problems involving anisotropic media and metamaterials, proper algorithms should be developed. It has been proved that discretizing in quadratic triangular elements may lead to an improved accuracy. Here we present a collection of elemental matrices evaluated analytically for quadratic triangular elements. They could be useful for the finite element method in advanced electromagnetics.
The Total Number of Parameters in the Finite Element ...
African Journals Online (AJOL)
Rectangular finite elements are important in Finite Element Method. This paper establishes a general formula for obtaining the total number of parameters associated with any finite element rectangulation of a domain. This number is also the dimension of the trail space as well as the size of the associated linear system.
Parallel direct solver for finite element modeling of manufacturing processes
DEFF Research Database (Denmark)
Nielsen, Chris Valentin; Martins, P.A.F.
2017-01-01
The central processing unit (CPU) time is of paramount importance in finite element modeling of manufacturing processes. Because the most significant part of the CPU time is consumed in solving the main system of equations resulting from finite element assemblies, different approaches have been...... developed to optimize solutions and reduce the overall computational costs of large finite element models....
Finite element modelling of helmeted head impact under frontal ...
Indian Academy of Sciences (India)
Abstract. Finite element models of the head and helmet were used to study contact forces during frontal impact of the head with a rigid surface. The finite element model of the head consists of skin, skull, cerebro-spinal fluid (CSF), brain, tentorium and falx. The finite element model of the helmet consists of shell and foam.
Adaptive Smoothed Finite Elements (ASFEM) for history dependent material models
Quak, W.; van den Boogaard, Antonius H.; Menary, Gary
2011-01-01
A successful simulation of a bulk forming process with finite elements can be difficult due to distortion of the finite elements. Nodal smoothed Finite Elements (NSFEM) are an interesting option for such a process since they show good distortion insensitivity and moreover have locking-free behavior
The finite element method in engineering, 2nd edition
International Nuclear Information System (INIS)
Rao, S.S.
1986-01-01
This work provides a systematic introduction to the various aspects of the finite element method as applied to engineering problems. Contents include: introduction to finite element method; solution of finite element equations; solid and structural mechanics; static analysis; dynamic analysis; heat transfer; fluid mechanics and additional applications
Finite element simulation and testing of ISW CFRP anchorage
DEFF Research Database (Denmark)
Schmidt, Jacob Wittrup; Goltermann, Per; Hertz, Kristian Dahl
2013-01-01
is modelled in the 3D finite Element program ABAQUS, just as digital image correlation (DIC) testing was performed to verify the finite element simulation. Also a new optimized design was produced to ensure that the finite element simulation and anchorage behaviour correlated well. It is seen...
Modelling bucket excavation by finite element
Pecingina, O. M.
2015-11-01
Changes in geological components of the layers from lignite pits have an impact on the sustainability of the cup path elements and under the action of excavation force appear efforts leading to deformation of the entire assembly. Application of finite element method in the optimization of components leads to economic growth, to increase the reliability and durability of the studied machine parts thus the machine. It is obvious usefulness of knowledge the state of mechanical tensions that the designed piece or the assembly not to break under the action of tensions that must cope during operation. In the course of excavation work on all bucket cutting force components, the first coming into contact with the material being excavated cutting edge. Therefore in the study with finite element analysis is retained only cutting edge. To study the field of stress and strain on the cutting edge will be created geometric patterns for each type of cup this will be subject to static analysis. The geometric design retains the cutting edge shape and on this on the tooth cassette location will apply an areal force on the abutment tooth. The cutting edge real pattern is subjected to finite element study for the worst case of rock cutting by symmetrical and asymmetrical cups whose profile is different. The purpose of this paper is to determine the displacement and tensions field for both profiles considering the maximum force applied on the cutting edge and the depth of the cutting is equal with the width of the cutting edge of the tooth. It will consider the worst case when on the structure will act both the tangential force and radial force on the bucket profile. For determination of stress and strain field on the form design of cutting edge profile will apply maximum force assuming uniform distribution and on the edge surface force will apply a radial force. After geometric patterns discretization on the cutting knives and determining stress field, can be seen that at the
Finite rotation shells basic equations and finite elements for Reissner kinematics
Wisniewski, K
2010-01-01
This book covers theoretical and computational aspects of non-linear shells. Several advanced topics of shell equations and finite elements - not included in standard textbooks on finite elements - are addressed, and the book includes an extensive bibliography.
Giraldo, F.; Mueller, A.; Kopera, M. A.; Abdi, D. S.; Wilcox, L.
2015-12-01
In this talk, we shall describe the NUMA atmospheric model, focusing in particular on its unified continuous/discontinuous (CG and DG) Galerkin numerical methods that are used to represent the spatial derivatives. We shall describe how these two methods are formulated in a unified approach and the advantages that this brings. We will also report on the progress in extending NUMA to using adaptive mesh refinement. Lastly, we will report on the scalability and performance of NUMA on the leadership computing facilities (LCF) of the Department of Energy where we have scaled NUMA to over 3 million MPI threads achieving a 90% efficiency.
Energy Technology Data Exchange (ETDEWEB)
Lee, Sang Jin; Seo, Jeong Moon
2000-08-01
The main goal of this research is to establish a methodology of finite element analysis of containment building predicting not only global behaviour but also local failure mode. In this report, we summerize some existing numerical analysis techniques to be improved for containment building. In other words, a complete description of the standard degenerated shell finite element formulation is provided for nonlinear stress analysis of nuclear containment structure. A shell finite element is derived using the degenerated solid concept which does not rely on a specific shell theory. Reissner-Mindlin assumptions are adopted to consider the transverse shear deformation effect. In order to minimize the sensitivity of the constitutive equation to structural types, microscopic material model is adopted. The four solution algorithms based on the standard Newton-Raphson method are discussed. Finally, two numerical examples are carried out to test the performance of the adopted shell medel.
Finite element simulation of piezoelectric transformers.
Tsuchiya, T; Kagawa, Y; Wakatsuki, N; Okamura, H
2001-07-01
Piezoelectric transformers are nothing but ultrasonic resonators with two pairs of electrodes provided on the surface of a piezoelectric substrate in which electrical energy is carried in the mechanical form. The input and output electrodes are arranged to provide the impedance transformation, which results in the voltage transformation. As they are operated at a resonance, the electrical equivalent circuit approach has traditionally been developed in a rather empirical way and has been used for analysis and design. The present paper deals with the analysis of the piezoelectric transformers based on the three-dimensional finite element modelling. The PIEZO3D code that we have developed is modified to include the external loading conditions. The finite element approach is now available for a wide variety of the electrical boundary conditions. The equivalent circuit of lumped parameters can also be derived from the finite element method (FEM) solution if required. The simulation of the present transformers is made for the low intensity operation and compared with the experimental results. Demonstration is made for basic Rosen-type transformers in which the longitudinal mode of a plate plays an important role; in which the equivalent circuit of lumped constants has been used. However, there are many modes of vibration associated with the plate, the effect of which cannot always be ignored. In the experiment, the double resonances are sometimes observed in the vicinity of the operating frequency. The simulation demonstrates that this is due to the coupling of the longitudinal mode with the flexural mode. Thus, the simulation provides an invaluable guideline to the transformer design.
A finite element method for neutron transport
International Nuclear Information System (INIS)
Ackroyd, R.T.
1978-01-01
A variational treatment of the finite element method for neutron transport is given based on a version of the even-parity Boltzmann equation which does not assume that the differential scattering cross-section has a spherical harmonic expansion. The theory of minimum and maximum principles is based on the Cauchy-Schwartz equality and the properties of a leakage operator G and a removal operator C. For systems with extraneous sources, two maximum and one minimum principles are given in boundary free form, to ease finite element computations. The global error of an approximate variational solution is given, the relationship of one the maximum principles to the method of least squares is shown, and the way in which approximate solutions converge locally to the exact solution is established. A method for constructing local error bounds is given, based on the connection between the variational method and the method of the hypercircle. The source iteration technique and a maximum principle for a system with extraneous sources suggests a functional for a variational principle for a self-sustaining system. The principle gives, as a consequence of the properties of G and C, an upper bound to the lowest eigenvalue. A related functional can be used to determine both upper and lower bounds for the lowest eigenvalue from an inspection of any approximate solution for the lowest eigenfunction. The basis for the finite element is presented in a general form so that two modes of exploitation can be undertaken readily. The model can be in phase space, with positional and directional co-ordinates defining points of the model, or it can be restricted to the positional co-ordinates and an expansion in orthogonal functions used for the directional co-ordinates. Suitable sets of functions are spherical harmonics and Walsh functions. The latter set is appropriate if a discrete direction representation of the angular flux is required. (author)
Gelinas, R. J.; Doss, S. K.; Vajk, J. P.; Djomehri, J.; Miller, K.
1983-01-01
The mathematical background regarding the moving finite element (MFE) method of Miller and Miller (1981) is discussed, taking into account a general system of partial differential equations (PDE) and the amenability of the MFE method in two dimensions to code modularization and to semiautomatic user-construction of numerous PDE systems for both Dirichlet and zero-Neumann boundary conditions. A description of test problem results is presented, giving attention to aspects of single square wave propagation, and a solution of the heat equation.
Finite Element analysis of jar connections
DEFF Research Database (Denmark)
Kristensen, A.; Toor, Kashif; Solem, Sigurd
2005-01-01
A new tool joint system is considered. Traditionally these rotary connections have been designed with only one shoulder geometry. However, in order to increase the torque rating of the tool joint, a new design is introduced using two shoulders. This design allow reduced tool joint dimensions wher...... whereby down-hole equipment more easily can be fitted. In order to evaluate the validity of the design, finite element analysis have been performed in ANSYS. The results obtained indicate that the new design is valid and further tests can be performed....
Finite element modeling methods for photonics
Rahman, B M Azizur
2013-01-01
The term photonics can be used loosely to refer to a vast array of components, devices, and technologies that in some way involve manipulation of light. One of the most powerful numerical approaches available to engineers developing photonic components and devices is the Finite Element Method (FEM), which can be used to model and simulate such components/devices and analyze how they will behave in response to various outside influences. This resource provides a comprehensive description of the formulation and applications of FEM in photonics applications ranging from telecommunications, astron
Mixed finite elements for global tide models.
Cotter, Colin J; Kirby, Robert C
2016-01-01
We study mixed finite element methods for the linearized rotating shallow water equations with linear drag and forcing terms. By means of a strong energy estimate for an equivalent second-order formulation for the linearized momentum, we prove long-time stability of the system without energy accumulation-the geotryptic state. A priori error estimates for the linearized momentum and free surface elevation are given in [Formula: see text] as well as for the time derivative and divergence of the linearized momentum. Numerical results confirm the theoretical results regarding both energy damping and convergence rates.
Generalized multiscale finite element methods: Oversampling strategies
Efendiev, Yalchin R.
2014-01-01
In this paper, we propose oversampling strategies in the generalized multiscale finite element method (GMsFEM) framework. The GMsFEM, which has been recently introduced in Efendiev et al. (2013b) [Generalized Multiscale Finite Element Methods, J. Comput. Phys., vol. 251, pp. 116-135, 2013], allows solving multiscale parameter-dependent problems at a reduced computational cost by constructing a reduced-order representation of the solution on a coarse grid. The main idea of the method consists of (1) the construction of snapshot space, (2) the construction of the offline space, and (3) construction of the online space (the latter for parameter-dependent problems). In Efendiev et al. (2013b) [Generalized Multiscale Finite Element Methods, J. Comput. Phys., vol. 251, pp. 116-135, 2013], it was shown that the GMsFEM provides a flexible tool to solve multiscale problems with a complex input space by generating appropriate snapshot, offline, and online spaces. In this paper, we develop oversampling techniques to be used in this context (see Hou and Wu (1997) where oversampling is introduced for multiscale finite element methods). It is known (see Hou and Wu (1997)) that the oversampling can improve the accuracy of multiscale methods. In particular, the oversampling technique uses larger regions (larger than the target coarse block) in constructing local basis functions. Our motivation stems from the analysis presented in this paper, which shows that when using oversampling techniques in the construction of the snapshot space and offline space, GMsFEM will converge independent of small scales and high contrast under certain assumptions. We consider the use of a multiple eigenvalue problems to improve the convergence and discuss their relation to single spectral problems that use oversampled regions. The oversampling procedures proposed in this paper differ from those in Hou and Wu (1997). In particular, the oversampling domains are partially used in constructing local
Approximation of functionally graded plates with non-conforming finite elements
Chinosi, Claudia; Della Croce, Lucia
2007-12-01
In this paper rectangular plates made of functionally graded materials (FGMs) are studied. A two-constituent material distribution through the thickness is considered, varying with a simple power rule of mixture. The equations governing the FGM plates are determined using a variational formulation arising from the Reissner-Mindlin theory. To approximate the problem a simple locking-free Discontinuous Galerkin finite element of non-conforming type is used, choosing a piecewise linear non-conforming approximation for both rotations and transversal displacement. Several numerical simulations are carried out in order to show the capability of the proposed element to capture the properties of plates of various gradings, subjected to thermo-mechanical loads.
Spilker, R L; Suh, J K; Mow, V C
1990-05-01
A finite element analysis is used to study a previously unresolved issue of the effects of platen-specimen friction on the response of the unconfined compression test; effects of platen permeability are also determined. The finite element formulation is based on the linear KLM biphasic model for articular cartilage and other hydrated soft tissues. A Galerkin weighted residual method is applied to both the solid phase and the fluid phase, and the continuity equation for the intrinsically incompressible binary mixture is introduced via a penalty method. The solid phase displacements and fluid phase velocities are interpolated for each element in terms of unknown nodal values, producing a system of first order differential equations which are solved using a standard numerical finite difference technique. An axisymmetric element of quadrilateral cross-section is developed and applied to the mechanical test problem of a cylindrical specimen of soft tissue in unconfined compression. These studies show that interfacial friction plays a major role in the unconfined compression response of articular cartilage specimens with small thickness to diameter ratios.
A multigrid solution method for mixed hybrid finite elements
Energy Technology Data Exchange (ETDEWEB)
Schmid, W. [Universitaet Augsburg (Germany)
1996-12-31
We consider the multigrid solution of linear equations arising within the discretization of elliptic second order boundary value problems of the form by mixed hybrid finite elements. Using the equivalence of mixed hybrid finite elements and non-conforming nodal finite elements, we construct a multigrid scheme for the corresponding non-conforming finite elements, and, by this equivalence, for the mixed hybrid finite elements, following guidelines from Arbogast/Chen. For a rectangular triangulation of the computational domain, this non-conforming schemes are the so-called nodal finite elements. We explicitly construct prolongation and restriction operators for this type of non-conforming finite elements. We discuss the use of plain multigrid and the multilevel-preconditioned cg-method and compare their efficiency in numerical tests.
Finite element modeling of piezoelectric elements with complex electrode configuration
International Nuclear Information System (INIS)
Paradies, R; Schläpfer, B
2009-01-01
It is well known that the material properties of piezoelectric materials strongly depend on the state of polarization of the individual element. While an unpolarized material exhibits mechanically isotropic material properties in the absence of global piezoelectric capabilities, the piezoelectric material properties become transversally isotropic with respect to the polarization direction after polarization. Therefore, for evaluating piezoelectric elements the material properties, including the coupling between the mechanical and the electromechanical behavior, should be addressed correctly. This is of special importance for the micromechanical description of piezoelectric elements with interdigitated electrodes (IDEs). The best known representatives of this group are active fiber composites (AFCs), macro fiber composites (MFCs) and the radial field diaphragm (RFD), respectively. While the material properties are available for a piezoelectric wafer with a homogeneous polarization perpendicular to its plane as postulated in the so-called uniform field model (UFM), the same information is missing for piezoelectric elements with more complex electrode configurations like the above-mentioned ones with IDEs. This is due to the inhomogeneous field distribution which does not automatically allow for the correct assignment of the material, i.e. orientation and property. A variation of the material orientation as well as the material properties can be accomplished by including the polarization process of the piezoelectric transducer in the finite element (FE) simulation prior to the actual load case to be investigated. A corresponding procedure is presented which automatically assigns the piezoelectric material properties, e.g. elasticity matrix, permittivity, and charge vector, for finite element models (FEMs) describing piezoelectric transducers according to the electric field distribution (field orientation and strength) in the structure. A corresponding code has been
Shobeiri, Vahid
2016-03-01
In this article, the bi-directional evolutionary structural optimization (BESO) method based on the element-free Galerkin (EFG) method is presented for topology optimization of continuum structures. The mathematical formulation of the topology optimization is developed considering the nodal strain energy as the design variable and the minimization of compliance as the objective function. The EFG method is used to derive the shape functions using the moving least squares approximation. The essential boundary conditions are enforced by the method of Lagrange multipliers. Several topology optimization problems are presented to show the effectiveness of the proposed method. Many issues related to topology optimization of continuum structures, such as chequerboard patterns and mesh dependency, are studied in the examples.
Modeling shallow water flows using the discontinuous Galerkin method
Khan, Abdul A
2014-01-01
Replacing the Traditional Physical Model Approach Computational models offer promise in improving the modeling of shallow water flows. As new techniques are considered, the process continues to change and evolve. Modeling Shallow Water Flows Using the Discontinuous Galerkin Method examines a technique that focuses on hyperbolic conservation laws and includes one-dimensional and two-dimensional shallow water flows and pollutant transports. Combines the Advantages of Finite Volume and Finite Element Methods This book explores the discontinuous Galerkin (DG) method, also known as the discontinuous finite element method, in depth. It introduces the DG method and its application to shallow water flows, as well as background information for implementing and applying this method for natural rivers. It considers dam-break problems, shock wave problems, and flows in different regimes (subcritical, supercritical, and transcritical). Readily Adaptable to the Real World While the DG method has been widely used in the fie...
Modeling shallow water flows using the discontinuous galerkin method
Khan, Abdul A
2014-01-01
Replacing the Traditional Physical Model Approach Computational models offer promise in improving the modeling of shallow water flows. As new techniques are considered, the process continues to change and evolve. Modeling Shallow Water Flows Using the Discontinuous Galerkin Method examines a technique that focuses on hyperbolic conservation laws and includes one-dimensional and two-dimensional shallow water flows and pollutant transports. Combines the Advantages of Finite Volume and Finite Element Methods This book explores the discontinuous Galerkin (DG) method, also known as the discontinuous finite element method, in depth. It introduces the DG method and its application to shallow water flows, as well as background information for implementing and applying this method for natural rivers. It considers dam-break problems, shock wave problems, and flows in different regimes (subcritical, supercritical, and transcritical). Readily Adaptable to the Real World While the DG method has been widely used in the fie...
Solution of the diffusion equations for several groups by the finite elements method
International Nuclear Information System (INIS)
Arredondo S, C.
1975-01-01
The code DELFIN has been implemented for the solution of the neutrons diffusion equations in two dimensions obtained by applying the approximation of several groups of energy. The code works with any number of groups and regions, and can be applied to thermal reactors as well as fast reactor. Providing it with the diffusion coefficients, the effective sections and the fission spectrum we obtain the results for the systems multiplying constant and the flows of each groups. The code was established using the method of finite elements, which is a form of resolution of the variational formulation of the equations applying the Ritz-Galerkin method with continuous polynomial functions by parts, in one case of the Lagrange type with rectangular geometry and up to the third grade. The obtained results and the comparison with the results in the literature, permit to reach the conclusion that it is convenient, to use the rectangular elements in all the cases where the geometry permits it, and demonstrate also that the finite elements method is better than the finite differences method. (author)
A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners
Watson, W. R.; Jones, M. G.
2009-01-01
A non-modal finite element theory for predicting the attenuation of an extended-reacting liner containing a porous facesheet and located in a no-flow duct is presented. The mathematical approach is to solve separate wave equations in the liner and duct airway and to couple these two solutions by invoking kinematic constraints at the facesheet that are consistent with a continuum theory of fluid motion. Given the liner intrinsic properties, a weak Galerkin finite element formulation with cubic polynomial basis functions is used as the basis for generating a discrete system of acoustic equations that are solved to obtain the coupled acoustic field. A state-of-the-art, asymmetric, parallel, sparse equation solver is implemented that allows tens of thousands of grid points to be analyzed. A grid refinement study is presented to show that the predicted attenuation converges. Excellent comparison of the numerically predicted attenuation to that of a mode theory (using a Haynes 25 metal foam liner) is used to validate the computational approach. Simulations are also presented for fifteen porous plate, extended-reacting liners. The construction of some of the porous plate liners suggest that they should behave as resonant liners while the construction of others suggest that they should behave as broadband attenuators. In each case the finite element theory is observed to predict the proper attenuation trend.
Nodally Integrated Finite Element Formulation for Mindlin-Reissner Plates
Simoes, D. A.; Jadhav, T. A.
2014-01-01
This work describes a nodally integrated finite element formulation for plates under the Mindlin-Reissner theory. The formulation makes use of the weighted residual method and nodal integration to derive the assumed strain relations. An element formulation for four-node quadrilateral elements is implemented in the nonlinear finite element solver Abaqus using the UEL user element subroutine. Numerical tests are carried out on the new element and the results are presented.
A finite element model for protein transport in vivo
Directory of Open Access Journals (Sweden)
Montas Hubert J
2007-06-01
Full Text Available Abstract Background Biological mass transport processes determine the behavior and function of cells, regulate interactions between synthetic agents and recipient targets, and are key elements in the design and use of biosensors. Accurately predicting the outcomes of such processes is crucial to both enhancing our understanding of how these systems function, enabling the design of effective strategies to control their function, and verifying that engineered solutions perform according to plan. Methods A Galerkin-based finite element model was developed and implemented to solve a system of two coupled partial differential equations governing biomolecule transport and reaction in live cells. The simulator was coupled, in the framework of an inverse modeling strategy, with an optimization algorithm and an experimental time series, obtained by the Fluorescence Recovery after Photobleaching (FRAP technique, to estimate biomolecule mass transport and reaction rate parameters. In the inverse algorithm, an adaptive method was implemented to calculate sensitivity matrix. A multi-criteria termination rule was developed to stop the inverse code at the solution. The applicability of the model was illustrated by simulating the mobility and binding of GFP-tagged glucocorticoid receptor in the nucleoplasm of mouse adenocarcinoma. Results The numerical simulator shows excellent agreement with the analytic solutions and experimental FRAP data. Detailed residual analysis indicates that residuals have zero mean and constant variance and are normally distributed and uncorrelated. Therefore, the necessary and sufficient criteria for least square parameter optimization, which was used in this study, were met. Conclusion The developed strategy is an efficient approach to extract as much physiochemical information from the FRAP protocol as possible. Well-posedness analysis of the inverse problem, however, indicates that the FRAP protocol provides insufficient
Computational structural analysis and finite element methods
Kaveh, A
2014-01-01
Graph theory gained initial prominence in science and engineering through its strong links with matrix algebra and computer science. Moreover, the structure of the mathematics is well suited to that of engineering problems in analysis and design. The methods of analysis in this book employ matrix algebra, graph theory and meta-heuristic algorithms, which are ideally suited for modern computational mechanics. Efficient methods are presented that lead to highly sparse and banded structural matrices. The main features of the book include: application of graph theory for efficient analysis; extension of the force method to finite element analysis; application of meta-heuristic algorithms to ordering and decomposition (sparse matrix technology); efficient use of symmetry and regularity in the force method; and simultaneous analysis and design of structures.
Quality management of finite element analysis
Barlow, John
1991-09-01
A quality management system covering the use of finite element analysis is described. The main topics are as follows: acquisition, development and verification of software (including the software suppliers software quality control system), support, documentation, error control, internal software, software acceptance and release; development and qualification of analysis methods, including software evaluation, analysis procedure qualification and documentation, procedure quality checks, control of analysis procedure errors; product design and integrity analysis, including project quality assurance and analysis planning, task specification and allocation, analysis, execution, results checking and analysis records. Other issues include the commercial and business advantages of quality systems, project and technical management and the training and experience of personnel. The items are correlated with the requirements of International Standard Organization 9001.
Adaptive finite element method for shape optimization
Morin, Pedro
2012-01-16
We examine shape optimization problems in the context of inexact sequential quadratic programming. Inexactness is a consequence of using adaptive finite element methods (AFEM) to approximate the state and adjoint equations (via the dual weighted residual method), update the boundary, and compute the geometric functional. We present a novel algorithm that equidistributes the errors due to shape optimization and discretization, thereby leading to coarse resolution in the early stages and fine resolution upon convergence, and thus optimizing the computational effort. We discuss the ability of the algorithm to detect whether or not geometric singularities such as corners are genuine to the problem or simply due to lack of resolution - a new paradigm in adaptivity. © EDP Sciences, SMAI, 2012.
Finite-Element Modelling of Biotransistors
Directory of Open Access Journals (Sweden)
Selvaganapathy PR
2010-01-01
Full Text Available Abstract Current research efforts in biosensor design attempt to integrate biochemical assays with semiconductor substrates and microfluidic assemblies to realize fully integrated lab-on-chip devices. The DNA biotransistor (BioFET is an example of such a device. The process of chemical modification of the FET and attachment of linker and probe molecules is a statistical process that can result in variations in the sensed signal between different BioFET cells in an array. In order to quantify these and other variations and assess their importance in the design, complete physical simulation of the device is necessary. Here, we perform a mean-field finite-element modelling of a short channel, two-dimensional BioFET device. We compare the results of this model with one-dimensional calculation results to show important differences, illustrating the importance of the molecular structure, placement and conformation of DNA in determining the output signal.
Static-kinematic duality in beams, plates, shells and its central role in the finite element method
Directory of Open Access Journals (Sweden)
Carpinteri Alberto
2017-01-01
Full Text Available Static and kinematic matrix operator equations are revisited for one-, two-, and three-dimensional deformable bodies. In particular, the elastic problem is formulated in the details in the case of arches, cylinders, circular plates, thin domes, and, through an induction process, shells of revolution. It is emphasized how the static and kinematic matrix operators are one the adjoint of the other, and then demonstrated through the definition of stiffness matrix and the application of virtual work principle. From the matrix operator formulation it clearly emerges the identity of the usual Finite Element Method definition of elastic stiffness matrix and the classical definition of Ritz-Galerkin matrix.
Static-kinematic duality in beams, plates, shells and its central role in the finite element method
Carpinteri, Alberto
2017-01-01
Static and kinematic matrix operator equations are revisited for one-, two-, and three-dimensional deformable bodies. In particular, the elastic problem is formulated in the details in the case of arches, cylinders, circular plates, thin domes, and, through an induction process, shells of revolution. It is emphasized how the static and kinematic matrix operators are one the adjoint of the other, and then demonstrated through the definition of stiffness matrix and the application of virtual work principle. From the matrix operator formulation it clearly emerges the identity of the usual Finite Element Method definition of elastic stiffness matrix and the classical definition of Ritz-Galerkin matrix.
International Nuclear Information System (INIS)
Cunha Furtado, F. da; Galeao, A.C.N.R.
1984-01-01
A numerical procedure for the integration of the incompressible Navier-Stokes equations, when expressed in terms of a stream function equation and a vorticity transport equation, is presented. This procedure comprises: the variational formulation of the equations, the construction of the approximation spaces by the finite element method and the discretization via the Galerkin method. For the stationary problems, the system of non-linear algebraic equations resulting from the discretization is solved by the Newton-Raphson algorithm. Finally, for the transient problems, the solution of the non-linear ordinary differential equations resulting from the spatial discretization is accomplished through a Crank-Nicolson scheme. (Author) [pt
On the design of general-purpose flux limiters for finite element schemes. I. Scalar convection
Kuzmin, D.
2006-12-01
The algebraic flux correction (AFC) paradigm is extended to finite element discretizations with a consistent mass matrix. It is shown how to render an implicit Galerkin scheme positivity-preserving and remove excessive artificial diffusion in regions where the solution is sufficiently smooth. To this end, the original discrete operators are modified in a mass-conserving fashion so as to enforce the algebraic constraints to be satisfied by the numerical solution. A node-oriented limiting strategy is employed to control the raw antidiffusive fluxes which consist of a convective part and a contribution of the consistent mass matrix. The former offsets the artificial diffusion due to 'upwinding' of the spatial differential operator and lends itself to an upwind-biased flux limiting. The latter eliminates the error induced by mass lumping and calls for the use of a symmetric flux limiter. The concept of a target flux and a new definition of upper/lower bounds make it possible to combine the advantages of algebraic FCT and TVD schemes introduced previously by the author and his coworkers. Unlike other high-resolution schemes for unstructured meshes, the new algorithm reduces to a consistent (high-order) Galerkin scheme in smooth regions and is designed to provide an optimal treatment of both stationary and time-dependent problems. Its performance is illustrated by application to the linear advection equation for a number of 1D and 2D configurations.
Finite Element Solutions for the Space Fractional Diffusion Equation with a Nonlinear Source Term
Directory of Open Access Journals (Sweden)
Y. J. Choi
2012-01-01
Full Text Available We consider finite element Galerkin solutions for the space fractional diffusion equation with a nonlinear source term. Existence, stability, and order of convergence of approximate solutions for the backward Euler fully discrete scheme have been discussed as well as for the semidiscrete scheme. The analytical convergent orders are obtained as O(k+hγ˜, where γ˜ is a constant depending on the order of fractional derivative. Numerical computations are presented, which confirm the theoretical results when the equation has a linear source term. When the equation has a nonlinear source term, numerical results show that the diffusivity depends on the order of fractional derivative as we expect.
Error analysis of semidiscrete finite element methods for inhomogeneous time-fractional diffusion
Jin, B.
2014-05-30
© 2014 Published by Oxford University Press on behalf of the Institute of Mathematics and its Applications. All rights reserved. We consider the initial-boundary value problem for an inhomogeneous time-fractional diffusion equation with a homogeneous Dirichlet boundary condition, a vanishing initial data and a nonsmooth right-hand side in a bounded convex polyhedral domain. We analyse two semidiscrete schemes based on the standard Galerkin and lumped mass finite element methods. Almost optimal error estimates are obtained for right-hand side data f (x, t) ε L∞ (0, T; Hq(ω)), ≤1≥ 1, for both semidiscrete schemes. For the lumped mass method, the optimal L2(ω)-norm error estimate requires symmetric meshes. Finally, twodimensional numerical experiments are presented to verify our theoretical results.
Distributed parameter and finite element models for wave propagation in railway contact lines
Sorrentino, S.; Anastasio, D.; Fasana, A.; Marchesiello, S.
2017-12-01
A distributed parameter model of a railway two-level catenary system is presented for the analysis of the coupled wave dynamics. The wires are modelled as two straight axis parallel beams, with linear equilibrium equations, and the moving load applied by the pantograph is modelled as a constant concentrated travelling force. The general solution is sought by an application of the Ritz-Galerkin method, and then compared with direct time integrations of a finite element model (FEM), achieved by two different integration schemes. The proposed model provides a valid reference for appropriately selecting the FEM parameters, in order to reduce the errors due to spurious modes, affecting the numerical integrations especially at high speeds of the moving pantograph.
Stupishin, L. U.; Nikitin, K. E.; Kolesnikov, A. G.
2018-02-01
The article is concerned with a methodology of optimal design of geometrically nonlinear (flexible) shells of revolution of minimum weight with strength, stability and strain constraints. The problem of optimal design with constraints is reduced to the problem of unconstrained minimization using the penalty functions method. Stress-strain state of shell is determined within the geometrically nonlinear deformation theory. A special feature of the methodology is the use of a mixed finite-element formulation based on the Galerkin method. Test problems for determining the optimal form and thickness distribution of a shell of minimum weight are considered. The validity of the results obtained using the developed methodology is analyzed, and the efficiency of various optimization algorithms is compared to solve the set problem. The developed methodology has demonstrated the possibility and accuracy of finding the optimal solution.
The computation of linear triangular matrices in the finite element ...
African Journals Online (AJOL)
An algorithm is developed for generating the system matrices for the Finite Element Method of solving some classes of second order partial differential equations problems using the linear triangular elements. This algorithm reduces the complexity normally associated with the finite element approximation and makes the ...
Review on finite element method | Erhunmwun | Journal of Applied ...
African Journals Online (AJOL)
... finite elements, so that it is possible to systematically construct the approximation functions needed in a variational or weighted-residual approximation of the solution of a problem over each element. Keywords: Weak Formulation, Discretisation, Numerical methods, Finite element method, Global equations, Nodal solution ...
International Nuclear Information System (INIS)
Uchiyama, Tomomi; Minemura, Kiyoshi
1996-01-01
A two-dimensional finite element method for an incompressible two-fluid model is developed in order to enhance the applicability of the two-fluid model to the analyses in arbitrarily shaped channels. The solution algorithm is based on the SMAC method for incompressible single-phase flows. Quadrilateral element with four nodes is employed. The velocities and volumetric fractions of both phases are defined on the nodes. The pressure is defined at a centroid of each element and assumed to be constant within the element. The Galerkin method is applied to the finite element formulations. Air-water two-phase flow around a circular cylinder is analyzed by the finite element method. The calculated distributions of the volumetric fraction of the gas-phase show good agreement with measurements. (author)
Finite element analysis theory and application with ANSYS
Moaveni, Saeed
2015-01-01
For courses in Finite Element Analysis, offered in departments of Mechanical or Civil and Environmental Engineering. While many good textbooks cover the theory of finite element modeling, Finite Element Analysis: Theory and Application with ANSYS is the only text available that incorporates ANSYS as an integral part of its content. Moaveni presents the theory of finite element analysis, explores its application as a design/modeling tool, and explains in detail how to use ANSYS intelligently and effectively. Teaching and Learning Experience This program will provide a better teaching and learning experience-for you and your students. It will help: *Present the Theory of Finite Element Analysis: The presentation of theoretical aspects of finite element analysis is carefully designed not to overwhelm students. *Explain How to Use ANSYS Effectively: ANSYS is incorporated as an integral part of the content throughout the book. *Explore How to Use FEA as a Design/Modeling Tool: Open-ended design problems help stude...
Symplectic and multisymplectic schemes with the simple finite element method
International Nuclear Information System (INIS)
Zhen Liu; Bai Yongqiang; Li Qisheng; Wu Ke
2003-01-01
We study the numerical scheme of elliptic equations by the finite element method. With the special finite element domain, we can find that the scheme can keep a preserved symplectic structure in one-dimensional case and a preserved multisymplectic structure in two-dimensional case. Then we consider the discrete variational principle with the finite element method in the corresponding Lagrangian formalism for classical mechanics and field theory and get the symplectic or multisymplectic scheme of the Euler-Lagrangian equation
The finite element method its basis and fundamentals
Zienkiewicz, Olek C; Zhu, JZ
2013-01-01
The Finite Element Method: Its Basis and Fundamentals offers a complete introduction to the basis of the finite element method, covering fundamental theory and worked examples in the detail required for readers to apply the knowledge to their own engineering problems and understand more advanced applications. This edition sees a significant rearrangement of the book's content to enable clearer development of the finite element method, with major new chapters and sections added to cover: Weak forms Variational forms Multi-dimensional field prob
Finite Element Modeling of Burr Formation in Metal Cutting
Min, Sangkee; Dornfeld, David; Kim, J.; Shyu, B.
2007-01-01
In order to advance understanding of the burr formation process, a series of finite element models are introduced. First a finite element model of the burr formation of two-dimensional orthogonal cutting is introduced and validated with experimental observations. A detailed and thorough examination of the drilling burr forming process is undertaken. This information is then used in the construction of an analytical model and, leads to development of a three-dimensional finite element mode...
Impact of new computing systems on finite element computations
International Nuclear Information System (INIS)
Noor, A.K.; Fulton, R.E.; Storaasi, O.O.
1983-01-01
Recent advances in computer technology that are likely to impact finite element computations are reviewed. The characteristics of supersystems, highly parallel systems, and small systems (mini and microcomputers) are summarized. The interrelations of numerical algorithms and software with parallel architectures are discussed. A scenario is presented for future hardware/software environment and finite element systems. A number of research areas which have high potential for improving the effectiveness of finite element analysis in the new environment are identified
Taneja, Ankur; Higdon, Jonathan
2018-01-01
A high-order spectral element discontinuous Galerkin method is presented for simulating immiscible two-phase flow in petroleum reservoirs. The governing equations involve a coupled system of strongly nonlinear partial differential equations for the pressure and fluid saturation in the reservoir. A fully implicit method is used with a high-order accurate time integration using an implicit Rosenbrock method. Numerical tests give the first demonstration of high order hp spatial convergence results for multiphase flow in petroleum reservoirs with industry standard relative permeability models. High order convergence is shown formally for spectral elements with up to 8th order polynomials for both homogeneous and heterogeneous permeability fields. Numerical results are presented for multiphase fluid flow in heterogeneous reservoirs with complex geometric or geologic features using up to 11th order polynomials. Robust, stable simulations are presented for heterogeneous geologic features, including globally heterogeneous permeability fields, anisotropic permeability tensors, broad regions of low-permeability, high-permeability channels, thin shale barriers and thin high-permeability fractures. A major result of this paper is the demonstration that the resolution of the high order spectral element method may be exploited to achieve accurate results utilizing a simple cartesian mesh for non-conforming geological features. Eliminating the need to mesh to the boundaries of geological features greatly simplifies the workflow for petroleum engineers testing multiple scenarios in the face of uncertainty in the subsurface geology.
Energy Technology Data Exchange (ETDEWEB)
Kim, S. [Purdue Univ., West Lafayette, IN (United States)
1994-12-31
Parallel iterative procedures based on domain decomposition techniques are defined and analyzed for the numerical solution of wave propagation by finite element and finite difference methods. For finite element methods, in a Lagrangian framework, an efficient way for choosing the algorithm parameter as well as the algorithm convergence are indicated. Some heuristic arguments for finding the algorithm parameter for finite difference schemes are addressed. Numerical results are presented to indicate the effectiveness of the methods.
Solving nonlinear nonstationary problem of heat-conductivity by finite element method
Directory of Open Access Journals (Sweden)
Антон Янович Карвацький
2016-11-01
Full Text Available Methodology and effective solving algorithm of non-linear dynamic problems of thermal and electric conductivity with significant temperature dependence of thermal and physical properties are given on the basis of finite element method (FEM and Newton linearization method. Discrete equations system FEM was obtained with the use of Galerkin method, where the main function is the finite element form function. The methodology based on successive solving problems of thermal and electrical conductivity has been examined in the work in order to minimize the requirements for calculating resources (RAM. in particular. Having used Mathcad software original programming code was developed to solve the given problem. After investigation of the received results, comparative analyses of accurate solution data and results of numerical solutions, obtained with the use of Matlab programming products, was held. The geometry of one fourth part of the finite sized cylinder was used to test the given numerical model. The discretization of the calculation part was fulfilled using the open programming software for automated Gmsh nets with tetrahedral units, while ParaView, which is an open programming code as well, was used to visualize the calculation results. It was found out that the maximum value violation of potential and temperature determination doesn`t exceed 0,2-0,83% in the given work according to the problem conditions
Adaptive Smoothed Finite Elements (ASFEM) for history dependent material models
International Nuclear Information System (INIS)
Quak, W.; Boogaard, A. H. van den
2011-01-01
A successful simulation of a bulk forming process with finite elements can be difficult due to distortion of the finite elements. Nodal smoothed Finite Elements (NSFEM) are an interesting option for such a process since they show good distortion insensitivity and moreover have locking-free behavior and good computational efficiency. In this paper a method is proposed which takes advantage of the nodally smoothed field. This method, named adaptive smoothed finite elements (ASFEM), revises the mesh for every step of a simulation without mapping the history dependent material parameters. In this paper an updated-Lagrangian implementation is presented. Several examples are given to illustrate the method and to show its properties.
Finite element analysis in a minicomputer/mainframe environment
Storaasli, O. O.; Murphy, R. C.
1978-01-01
Design considerations were evaluated for general purpose finite element systems to maximize performance when installed on distributed computer hardware/software systems. It is shown how the features of current minicomputers complement those of a modular implementation of the finite element method for increasing the control, speed, and visibility (interactive graphics) in solving structural problems at reduced cost. The approach used is to implement a finite element system in a distributed computer environment to solve structural problems and to explore alternatives in distributing finite element computations.
Introduction to finite element analysis using MATLAB and Abaqus
Khennane, Amar
2013-01-01
There are some books that target the theory of the finite element, while others focus on the programming side of things. Introduction to Finite Element Analysis Using MATLAB(R) and Abaqus accomplishes both. This book teaches the first principles of the finite element method. It presents the theory of the finite element method while maintaining a balance between its mathematical formulation, programming implementation, and application using commercial software. The computer implementation is carried out using MATLAB, while the practical applications are carried out in both MATLAB and Abaqus. MA
An efficient finite element solution for gear dynamics
Cooley, C. G.; Parker, R. G.; Vijayakar, S. M.
2010-06-01
A finite element formulation for the dynamic response of gear pairs is proposed. Following an established approach in lumped parameter gear dynamic models, the static solution is used as the excitation in a frequency domain solution of the finite element vibration model. The nonlinear finite element/contact mechanics formulation provides accurate calculation of the static solution and average mesh stiffness that are used in the dynamic simulation. The frequency domain finite element calculation of dynamic response compares well with numerically integrated (time domain) finite element dynamic results and previously published experimental results. Simulation time with the proposed formulation is two orders of magnitude lower than numerically integrated dynamic results. This formulation admits system level dynamic gearbox response, which may include multiple gear meshes, flexible shafts, rolling element bearings, housing structures, and other deformable components.
A Finite Element Analysis of Optimal Variable Thickness Sheets
DEFF Research Database (Denmark)
Petersson, Joakim S
1996-01-01
A quasimixed Finite Element (FE) method for maximum stiffness of variablethickness sheets is analysed. The displacement is approximated with ninenode Lagrange quadrilateral elements and the thickness is approximated aselementwise constant. One is guaranteed that the FE displacement solutionswill ...
Control volume finite element method for radiation
International Nuclear Information System (INIS)
Ben Salah, M.; Askri, F.; Rousse, D.; Ben Nasrallah, S.
2005-01-01
In this paper a new methodology is presented by the authors for the numerical treatment of radiative heat transfer in emitting, absorbing and scattering media. This methodology is based on the utilisation of Control Volume Finite Element Method (CVFEM) and the use, for the first time, of matrix formulation of the discretized Radiative Transfer Equation (RTE). The advantages of the proposed methodology is to avoid problems that confronted when previous techniques are used to predict radiative heat transfer, essentially, in complex geometries and when there is scattering and/or non-black boundaries surfaces. Besides, the new formulation of the discretized RTE presented in this paper makes it possible to solve the algebraic system by direct or iterative numerical methods. The theoretical background of CVFEM and matrix formulation is presented in the text. The proposed technique is applied to different test problems, and the results compared favourably against other published works. Moreover this paper discusses in detail the effects of some radiative parameters, such as optical thickness and walls emissivities on the spatial evolution of the radiant heat flux. The numerical simulation of radiative heat transfer for different cases using the algorithm proposed in this work has shown that the developed computer procedure needs an accurate CPU time and is exempt of any numerical oscillations
Finite element modelling of composite castellated beam
Directory of Open Access Journals (Sweden)
Frans Richard
2017-01-01
Full Text Available Nowadays, castellated beam becomes popular in building structural as beam members. This is due to several advantages of castellated beam such as increased depth without any additional mass, passing the underfloor service ducts without changing of story elevation. However, the presence of holes can develop various local effects such as local buckling, lateral torsional buckling caused by compression force at the flange section of the steel beam. Many studies have investigated the failure mechanism of castellated beam and one technique which can prevent the beam fall into local failure is the use of reinforced concrete slab as lateral support on castellated beam, so called composite castellated beam. Besides of preventing the local failure of castellated beam, the concrete slab can increase the plasticity moment of the composite castellated beam section which can deliver into increasing the ultimate load of the beam. The aim of this numerical studies of composite castellated beam on certain loading condition (monotonic quasi-static loading. ABAQUS was used for finite element modelling purpose and compared with the experimental test for checking the reliability of the model. The result shows that the ultimate load of the composite castellated beam reached 6.24 times than the ultimate load of the solid I beam and 1.2 times compared the composite beam.
Finite element modeling of retinal prosthesis mechanics
Basinger, B. C.; Rowley, A. P.; Chen, K.; Humayun, M. S.; Weiland, J. D.
2009-10-01
Epiretinal prostheses used to treat degenerative retina diseases apply stimulus via an electrode array fixed to the ganglion cell side of the retina. Mechanical pressure applied by these arrays to the retina, both during initial insertion and throughout chronic use, could cause sufficient retinal damage to reduce the device's effectiveness. In order to understand and minimize potential mechanical damage, we have used finite element analysis to model mechanical interactions between an electrode array and the retina in both acute and chronic loading configurations. Modeling indicates that an acute tacking force distributes stress primarily underneath the tack site and heel edge of the array, while more moderate chronic stresses are distributed more evenly underneath the array. Retinal damage in a canine model chronically implanted with a similar array occurred in correlating locations, and model predictions correlate well with benchtop eyewall compression tests. This model provides retinal prosthesis researchers with a tool to optimize the mechanical electrode array design, but the techniques used here represent a unique effort to combine a modifiable device and soft biological tissues in the same model and those techniques could be extended to other devices that come into mechanical contact with soft neural tissues.
Probabilistic finite element modeling of waste rollover
International Nuclear Information System (INIS)
Khaleel, M.A.; Cofer, W.F.; Al-fouqaha, A.A.
1995-09-01
Stratification of the wastes in many Hanford storage tanks has resulted in sludge layers which are capable of retaining gases formed by chemical and/or radiolytic reactions. As the gas is produced, the mechanisms of gas storage evolve until the resulting buoyancy in the sludge leads to instability, at which point the sludge ''rolls over'' and a significant volume of gas is suddenly released. Because the releases may contain flammable gases, these episodes of release are potentially hazardous. Mitigation techniques are desirable for more controlled releases at more frequent intervals. To aid the mitigation efforts, a methodology for predicting of sludge rollover at specific times is desired. This methodology would then provide a rational basis for the development of a schedule for the mitigation procedures. In addition, a knowledge of the sensitivity of the sludge rollovers to various physical and chemical properties within the tanks would provide direction for efforts to reduce the frequency and severity of these events. In this report, the use of probabilistic finite element analyses for computing the probability of rollover and the sensitivity of rollover probability to various parameters is described
An Abaqus UEL implementation of the smoothed finite element method
Kumbhar, Pramod Y; Francis, Amrita; Swaminathan, Narasimhan; Annabattula, Ratna Kumar; Natarajan, Sundararajan
2017-01-01
In this paper, we discuss the implementation of a cell based smoothed finite element method (CSFEM) within the commercial finite element software Abaqus. The salient feature of the CSFEM is that it does not require an explicit form of the derivative of the shape functions and there is no isoparametric mapping. This implementation is accomplished by employing the user element subroutine (UEL) feature of the software. The details on the input data format together with the proposed user element ...
Hydrothermal analysis in engineering using control volume finite element method
Sheikholeslami, Mohsen
2015-01-01
Control volume finite element methods (CVFEM) bridge the gap between finite difference and finite element methods, using the advantages of both methods for simulation of multi-physics problems in complex geometries. In Hydrothermal Analysis in Engineering Using Control Volume Finite Element Method, CVFEM is covered in detail and applied to key areas of thermal engineering. Examples, exercises, and extensive references are used to show the use of the technique to model key engineering problems such as heat transfer in nanofluids (to enhance performance and compactness of energy systems),
An adaptive finite element model for hemolysis prediction
International Nuclear Information System (INIS)
Lacasse, D.; Pelletier, D.; Garon, A.
2002-01-01
This paper presents a numerical tool to predict hemolysis caused by mechanical devices that come in contact with blood. Predictions of hemolysis are obtained as a field variable H L , defined as an hemolysis index that takes values between 0 and 1. The field of H L is obtained by solving a hyperbolic equation using the discontinuous Galerkin method. Instead of solving the equation on a element by element basis, as normally done with this method, elemental systems are assembled in a global matrix. The technique readily yields solutions over any computational domain, including domains with recirculating flow fields. An adaptive remeshing algorithm is used to obtain grid independent results of hemolysis for a laminar flow of blood (Re=1600) in a pipe with a sudden contraction. (author)
Geotechnical Ultimate Limit State Design Using Finite Elements
Brinkgreve, R.B.J.; Post, M.
2015-01-01
Displacement-based finite element calculations are primarily used for serviceability limit state (SLS) analysis, but the finite element method also offers possibilities for ultimate limit state (ULS) design in geotechnical engineering. The combined use of SLS and ULS calculations with partial safety
Parallel direct solver for finite element modeling of manufacturing processes
DEFF Research Database (Denmark)
Nielsen, Chris Valentin; Martins, P.A.F.
2017-01-01
The central processing unit (CPU) time is of paramount importance in finite element modeling of manufacturing processes. Because the most significant part of the CPU time is consumed in solving the main system of equations resulting from finite element assemblies, different approaches have been d...
A geometric toolbox for tetrahedral finite element partitions
Brandts, J.; Korotov, S.; Křížek, M.; Axelsson, O.; Karátson, J.
2011-01-01
In this work we present a survey of some geometric results on tetrahedral partitions and their refinements in a unified manner. They can be used for mesh generation and adaptivity in practical calculations by the finite element method (FEM), and also in theoretical finite element (FE) analysis.
An introduction to the UNCLE finite element scheme
International Nuclear Information System (INIS)
Enderby, J.A.
1983-01-01
UNCLE is a completely general finite element scheme which provides common input, output, equation-solving and other facilities for a family of finite element codes for linear and non-linear stress analysis, heat transfer etc. This report describes the concepts on which UNCLE is based and gives a general account of the facilities provided. (author)
About the Finite Element Method Applied to Thick Plates
Directory of Open Access Journals (Sweden)
Mihaela Ibănescu
2006-01-01
Full Text Available The present paper approaches of plates subjected to transverse loads, when the shear force and the actual boundary conditions are considered, by using the Finite Element Method. The isoparametric finite elements create real facilities in formulating the problems and great possibilities in creating adequate computer programs.
Simulation of temperature distribution by finite element analysis on ...
Indian Academy of Sciences (India)
on exposure to the synchrotron beam has been simulated by finite element analysis. Design of the cooling mechanism for each of these components has been carried out and estimation of the temperature rise has also been done incorporating the cooling mechanism. Keywords. Synchrotron; EXAFS; finite element analysis.
THE PRACTICAL ANALYSIS OF FINITE ELEMENTS METHOD ERRORS
Directory of Open Access Journals (Sweden)
Natalia Bakhova
2011-03-01
Full Text Available Abstract. The most important in the practical plan questions of reliable estimations of finite elementsmethod errors are considered. Definition rules of necessary calculations accuracy are developed. Methodsand ways of the calculations allowing receiving at economical expenditures of computing work the best finalresults are offered.Keywords: error, given the accuracy, finite element method, lagrangian and hermitian elements.
Finite Element Method for Analysis of Material Properties
DEFF Research Database (Denmark)
Rauhe, Jens Christian
and the finite element method. The material microstructure of the heterogeneous material is non-destructively determined using X-ray microtomography. A software program has been generated which uses the X-ray tomographic data as an input for the mesh generation of the material microstructure. To obtain a proper...... description of the material microstructure the finite element models must contain a large number of elements and this problem is solved by using the preconditioned conjugated gradient solver with an Element-By-Element preconditioner. Finite element analysis provides the volume averaged stresses and strains...... which are used for the determination of the effective properties of the heterogeneous material. Generally, the properties determined using the finite element method coupled with X-ray microtomography are in good agreement with both experimentally determined properties and properties determined using...
Eversman, W.; Astley, R. J.
1981-01-01
The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.
Wavelet=Galerkin discretization of hyperbolic equations
Energy Technology Data Exchange (ETDEWEB)
Restrepo, J.M.; Leaf, G.K.
1994-12-31
The relative merits of the wavelet-Galerkin solution of hyperbolic partial differential equations, typical of geophysical problems, are quantitatively and qualitatively compared to traditional finite difference and Fourier-pseudo-spectral methods. The wavelet-Galerkin solution presented here is found to be a viable alternative to the two conventional techniques.
Legner, D; Skatulla, S; MBewu, J; Rama, R R; Reddy, B D; Sansour, C; Davies, N H; Franz, T
2014-03-01
Myocardial infarction is an increasing health problem worldwide. Because of an under-supply of blood, the cardiomyocytes in the affected region permanently lose their ability to contract. This in turn gradually weakens the overall heart function. A new therapeutic approach based on the injection of a gel into the infarcted area aims to support the healing and to inhibit adverse remodelling that can lead to heart failure. A computational model is the basis for obtaining a better understanding of the heart mechanics, in particular, how myocardial infarction and gel injections affect its pumping performance. A strain invariant-based stored energy function is proposed to account for the passive mechanical behaviour of the model, which also makes provision for active contraction. To incorporate injections an additive homogenization approach is introduced. The numerical framework is developed using an in-house code based on the element-free Galerkin method. The main focus of this contribution is to investigate the influence of gel injections on the mechanics of the left ventricle during the diastolic filling and systolic isovolumetric (isochoric) contraction phases. It is found that gel injections are able to reduce the elevated fibre stresses caused by an infarct. Copyright © 2013 John Wiley & Sons, Ltd.
Finite Element Analysis of Fluid-Conveying Timoshenko Pipes
Directory of Open Access Journals (Sweden)
Chih-Liang Chu
1995-01-01
Full Text Available A general finite element formulation using cubic Hermitian interpolation for dynamic analysis of pipes conveying fluid is presented. Both the effects of shearing deformations and rotary inertia are considered. The development retains the use of the classical four degrees-of-freedom for a two-node element. The effect of moving fluid is treated as external distributed forces on the support pipe and the fluid finite element matrices are derived from the virtual work done due to the fluid inertia forces. Finite element matrices for both the support pipe and moving fluid are derived and given explicitly. A numerical example is given to demonstrate the validity of the model.
Shear beams in finite element modelling : Software implementation and validation
Schreppers, G.J.; Hendriks, M.A.N.; Boer, A.; Ferreira, D.; Kikstra, W.P.
2015-01-01
Fiber models for beam and shell elements allow for relatively rapid finite element analysis of concrete structures and structural elements. This project aims at the development of the formulation of such elements and a pilot implementation. The reduction of calculation time and degrees of freedom
Space-time discontinuous Galerkin method for parabolic problems in time-dependent domains
Janivita Joto Sudirham, J.J.S.; Sudirham, J.J.; van der Vegt, Jacobus J.W.; van Damme, Rudolf M.J.
2004-01-01
In this report a space-time discontinuous Galerkin (DG) finite element method for the solution of the advection-diffusion-reaction equation in time-dependent domains is presented and analyzed. The variational formulation is based on a combination of the space-time DG method developed by van der Vegt
Discontinuous Galerkin Approximations for Computing Electromagnetic Bloch Modes in Photonic Crystals
Lu, Zhongjie; Cesmelioglu, A.; van der Vegt, Jacobus J.W.; Xu, Yan
We analyze discontinuous Galerkin finite element discretizations of the Maxwell equations with periodic coefficients. These equations are used to model the behavior of light in photonic crystals, which are materials containing a spatially periodic variation of the refractive index commensurate with
Finite element modeling of mass transport in high-Péclet cardiovascular flows
Hansen, Kirk; Arzani, Amirhossein; Shadden, Shawn
2016-11-01
Mass transport plays an important role in many important cardiovascular processes, including thrombus formation and atherosclerosis. These mass transport problems are characterized by Péclet numbers of up to 108, leading to several numerical difficulties. The presence of thin near-wall concentration boundary layers requires very fine mesh resolution in these regions, while large concentration gradients within the flow cause numerical stabilization issues. In this work, we will discuss some guidelines for solving mass transport problems in cardiovascular flows using a stabilized Galerkin finite element method. First, we perform mesh convergence studies in a series of idealized and patient-specific geometries to determine the required near-wall mesh resolution for these types of problems, using both first- and second-order tetrahedral finite elements. Second, we investigate the use of several boundary condition types at outflow boundaries where backflow during some parts of the cardiac cycle can lead to convergence issues. Finally, we evaluate the effect of reducing Péclet number by increasing mass diffusivity as has been proposed by some researchers. This work was supported by the NSF GRFP and NSF Career Award #1354541.
Dual-scale Galerkin methods for Darcy flow
Wang, Guoyin; Scovazzi, Guglielmo; Nouveau, Léo; Kees, Christopher E.; Rossi, Simone; Colomés, Oriol; Main, Alex
2018-02-01
The discontinuous Galerkin (DG) method has found widespread application in elliptic problems with rough coefficients, of which the Darcy flow equations are a prototypical example. One of the long-standing issues of DG approximations is the overall computational cost, and many different strategies have been proposed, such as the variational multiscale DG method, the hybridizable DG method, the multiscale DG method, the embedded DG method, and the Enriched Galerkin method. In this work, we propose a mixed dual-scale Galerkin method, in which the degrees-of-freedom of a less computationally expensive coarse-scale approximation are linked to the degrees-of-freedom of a base DG approximation. We show that the proposed approach has always similar or improved accuracy with respect to the base DG method, with a considerable reduction in computational cost. For the specific definition of the coarse-scale space, we consider Raviart-Thomas finite elements for the mass flux and piecewise-linear continuous finite elements for the pressure. We provide a complete analysis of stability and convergence of the proposed method, in addition to a study on its conservation and consistency properties. We also present a battery of numerical tests to verify the results of the analysis, and evaluate a number of possible variations, such as using piecewise-linear continuous finite elements for the coarse-scale mass fluxes.
A modified finite element procedure for underwater shock analysis
International Nuclear Information System (INIS)
Chan, S.K.
1990-01-01
Using the regular finite element method for analyzing wave propagation problems presents difficulties: (a) The finite element mesh gives spurious reflection of the traveling wave and (b) Since a finite element model has to have a finite boundary, the wave is reflected by the outside boundary. However, for underwater shock problems, only the response of the structure is of major interest, not the behavior of the wave itself, and the shock wave can be assumed to be spherical. By taking advantage of the limited scope of the underwater shock problem, a finite element procedure can be developed that eliminates the above difficulties. This procedure not only can give very accurate solutions but it may also include structural nonlinearities and effect of cavitation
On Finite Element Computations of Contact Problems in Micropolar Elasticity
Eremeyev, Victor A.; Skrzat, Andrzej; Stachowicz, Feliks
2016-01-01
Within the linear micropolar elasticity we discuss the development of new finite element and its implementation in commercial software. Here we implement the developed 8-node hybrid isoparametric element into ABAQUS and perform solutions of contact problems. We consider the contact of polymeric stamp modelled within the micropolar elasticity with an elastic substrate. The peculiarities of modelling of contact problems with a user defined finite element in ABAQUS are discussed. The provided co...
Finite element and boundary element applications in quantum mechanics
International Nuclear Information System (INIS)
Ueta, Tsuyoshi
2003-01-01
Although this book is one of the Oxford Texts in Applied and Engineering Mathematics, we may think of it as a physics book. It explains how to solve the problem of quantum mechanics using the finite element method (FEM) and the boundary element method (BEM). Many examples analysing actual problems are also shown. As for the ratio of the number of pages of FEM and BEM, the former occupies about 80%. This is, however, reasonable reflecting the flexibility of FEM. Although many explanations of FEM and BEM exist, most are written using special mathematical expressions and numerical computation fields. However, this book is written in the 'language of physicists' throughout. I think that it is very readable and easy to understand for physicists. In the derivation of FEM and the argument on calculation accuracy, the action integral and a variation principle are used consistently. In the numerical computation of matrices, such as simultaneous equations and eigen value problems, a description of important points is also fully given. Moreover, the practical problems which become important in the electron device design field and the condensed matter physics field are dealt with as example computations, so that this book is very practical and applicable. It is characteristic and interesting that FEM is applied to solve the Schroedinger and Poisson equations consistently, and to the solution of the Ginzburg--Landau equation in superconductivity. BEM is applied to treat electric field enhancements due to surface plasmon excitations at metallic surfaces. A number of references are cited at the end of all the chapters, and this is very helpful. The description of quantum mechanics is also made appropriately and the actual application of quantum mechanics in condensed matter physics can also be surveyed. In the appendices, the mathematical foundation, such as numerical quadrature formulae and Green's functions, is conveniently described. I recommend this book to those who need to
A Review of Element-Based Galerkin Methods for Numerical Weather Prediction
2015-04-01
is the method of choice for LAMs for several reasons. First, it is simple to implement on a Cartesian grid, especially if the curvature of the earth is... Mills , M., Smith, A.K., Tilmes, S., Vitt, F., Morrison, H., Cameron-Smith, P., Collins, W.D., Iacono, M.J., Easter, R.C., Ghan, S.J., Liu, X., Rasch...simulation and prediction model. part i: Model dynamics and verification. Meteorology and Atmospheric Physics 75, 161–193 (2000) 325. Yang , H.: Finite
Error Analysis for Discontinuous Galerkin Method for Parabolic Problems
Kaneko, Hideaki
2004-01-01
In the proposal, the following three objectives are stated: (1) A p-version of the discontinuous Galerkin method for a one dimensional parabolic problem will be established. It should be recalled that the h-version in space was used for the discontinuous Galerkin method. An a priori error estimate as well as a posteriori estimate of this p-finite element discontinuous Galerkin method will be given. (2) The parameter alpha that describes the behavior double vertical line u(sub t)(t) double vertical line 2 was computed exactly. This was made feasible because of the explicitly specified initial condition. For practical heat transfer problems, the initial condition may have to be approximated. Also, if the parabolic problem is proposed on a multi-dimensional region, the parameter alpha, for most cases, would be difficult to compute exactly even in the case that the initial condition is known exactly. The second objective of this proposed research is to establish a method to estimate this parameter. This will be done by computing two discontinuous Galerkin approximate solutions at two different time steps starting from the initial time and use them to derive alpha. (3) The third objective is to consider the heat transfer problem over a two dimensional thin plate. The technique developed by Vogelius and Babuska will be used to establish a discontinuous Galerkin method in which the p-element will be used for through thickness approximation. This h-p finite element approach, that results in a dimensional reduction method, was used for elliptic problems, but the application appears new for the parabolic problem. The dimension reduction method will be discussed together with the time discretization method.
Gleadall, Andrew; Pan, Jingzhe; Ding, Lifeng; Kruft, Marc-Anton; Curcó, David
2015-11-01
Molecular dynamics (MD) simulations are widely used to analyse materials at the atomic scale. However, MD has high computational demands, which may inhibit its use for simulations of structures involving large numbers of atoms such as amorphous polymer structures. An atomic-scale finite element method (AFEM) is presented in this study with significantly lower computational demands than MD. Due to the reduced computational demands, AFEM is suitable for the analysis of Young's modulus of amorphous polymer structures. This is of particular interest when studying the degradation of bioresorbable polymers, which is the topic of an accompanying paper. AFEM is derived from the inter-atomic potential energy functions of an MD force field. The nonlinear MD functions were adapted to enable static linear analysis. Finite element formulations were derived to represent interatomic potential energy functions between two, three and four atoms. Validation of the AFEM was conducted through its application to atomic structures for crystalline and amorphous poly(lactide). Copyright © 2015 Elsevier Ltd. All rights reserved.
MIB Galerkin method for elliptic interface problems.
Xia, Kelin; Zhan, Meng; Wei, Guo-Wei
2014-12-15
Material interfaces are omnipresent in the real-world structures and devices. Mathematical modeling of material interfaces often leads to elliptic partial differential equations (PDEs) with discontinuous coefficients and singular sources, which are commonly called elliptic interface problems. The development of high-order numerical schemes for elliptic interface problems has become a well defined field in applied and computational mathematics and attracted much attention in the past decades. Despite of significant advances, challenges remain in the construction of high-order schemes for nonsmooth interfaces, i.e., interfaces with geometric singularities, such as tips, cusps and sharp edges. The challenge of geometric singularities is amplified when they are associated with low solution regularities, e.g., tip-geometry effects in many fields. The present work introduces a matched interface and boundary (MIB) Galerkin method for solving two-dimensional (2D) elliptic PDEs with complex interfaces, geometric singularities and low solution regularities. The Cartesian grid based triangular elements are employed to avoid the time consuming mesh generation procedure. Consequently, the interface cuts through elements. To ensure the continuity of classic basis functions across the interface, two sets of overlapping elements, called MIB elements, are defined near the interface. As a result, differentiation can be computed near the interface as if there is no interface. Interpolation functions are constructed on MIB element spaces to smoothly extend function values across the interface. A set of lowest order interface jump conditions is enforced on the interface, which in turn, determines the interpolation functions. The performance of the proposed MIB Galerkin finite element method is validated by numerical experiments with a wide range of interface geometries, geometric singularities, low regularity solutions and grid resolutions. Extensive numerical studies confirm the
Finite Element Simulation of Blanking Process
Directory of Open Access Journals (Sweden)
Afzal Ahmed
2012-10-01
daya penembusan sebanyak 42%. Daya tebukan yang diukur melalui eksperimen dan simulasi kekal pada kira-kira 90kN melepasi penembusan punch sebanyak 62%. Apabila ketebalan keputusan kunci ditambah, ketinggian retak dikurangkan dan ini meningkatkan kualiti pengosongan.KEYWORDS: simulation; finite element simulation; blanking; computer aided manufacturing
A finite element calculation of flux pumping
Campbell, A. M.
2017-12-01
A flux pump is not only a fascinating example of the power of Faraday’s concept of flux lines, but also an attractive way of powering superconducting magnets without large electronic power supplies. However it is not possible to do this in HTS by driving a part of the superconductor normal, it must be done by exceeding the local critical density. The picture of a magnet pulling flux lines through the material is attractive, but as there is no direct contact between flux lines in the magnet and vortices, unless the gap between them is comparable to the coherence length, the process must be explicable in terms of classical electromagnetism and a nonlinear V-I characteristic. In this paper a simple 2D model of a flux pump is used to determine the pumping behaviour from first principles and the geometry. It is analysed with finite element software using the A formulation and FlexPDE. A thin magnet is passed across one or more superconductors connected to a load, which is a large rectangular loop. This means that the self and mutual inductances can be calculated explicitly. A wide strip, a narrow strip and two conductors are considered. Also an analytic circuit model is analysed. In all cases the critical state model is used, so the flux flow resistivity and dynamic resistivity are not directly involved, although an effective resistivity appears when J c is exceeded. In most of the cases considered here is a large gap between the theory and the experiments. In particular the maximum flux transferred to the load area is always less than the flux of the magnet. Also once the threshold needed for pumping is exceeded the flux in the load saturates within a few cycles. However the analytic circuit model allows a simple modification to allow for the large reduction in I c when the magnet is over a conductor. This not only changes the direction of the pumped flux but leads to much more effective pumping.
Precise magnetostatic field using the finite element method
International Nuclear Information System (INIS)
Nascimento, Francisco Rogerio Teixeira do
2013-01-01
The main objective of this work is to simulate electromagnetic fields using the Finite Element Method. Even in the easiest case of electrostatic and magnetostatic numerical simulation some problems appear when the nodal finite element is used. It is difficult to model vector fields with scalar functions mainly in non-homogeneous materials. With the aim to solve these problems two types of techniques are tried: the adaptive remeshing using nodal elements and the edge finite element that ensure the continuity of tangential components. Some numerical analysis of simple electromagnetic problems with homogeneous and non-homogeneous materials are performed using first, the adaptive remeshing based in various error indicators and second, the numerical solution of waveguides using edge finite element. (author)
Ruiz-Baier, Ricardo; Lunati, Ivan
2016-10-01
We present a novel discretization scheme tailored to a class of multiphase models that regard the physical system as consisting of multiple interacting continua. In the framework of mixture theory, we consider a general mathematical model that entails solving a system of mass and momentum equations for both the mixture and one of the phases. The model results in a strongly coupled and nonlinear system of partial differential equations that are written in terms of phase and mixture (barycentric) velocities, phase pressure, and saturation. We construct an accurate, robust and reliable hybrid method that combines a mixed finite element discretization of the momentum equations with a primal discontinuous finite volume-element discretization of the mass (or transport) equations. The scheme is devised for unstructured meshes and relies on mixed Brezzi-Douglas-Marini approximations of phase and total velocities, on piecewise constant elements for the approximation of phase or total pressures, as well as on a primal formulation that employs discontinuous finite volume elements defined on a dual diamond mesh to approximate scalar fields of interest (such as volume fraction, total density, saturation, etc.). As the discretization scheme is derived for a general formulation of multicontinuum physical systems, it can be readily applied to a large class of simplified multiphase models; on the other, the approach can be seen as a generalization of these models that are commonly encountered in the literature and employed when the latter are not sufficiently accurate. An extensive set of numerical test cases involving two- and three-dimensional porous media are presented to demonstrate the accuracy of the method (displaying an optimal convergence rate), the physics-preserving properties of the mixed-primal scheme, as well as the robustness of the method (which is successfully used to simulate diverse physical phenomena such as density fingering, Terzaghi's consolidation
A wave finite element analysis of the passive cochlea
Elliott, Stephen J.; Ni, Guangjian; Mace, Brian R.; Lineton, Ben
2013-01-01
Current models of the cochlea can be characterized as being either based on the assumed propagation of a single slow wave, which provides good insight, or involve the solution of a numerical model, such as in the finite element method, which allows the incorporation of more detailed anatomical features. In this paper it is shown how the wave finite element method can be used to decompose the results of a finite element calculation in terms of wave components, which allows the insight of the w...
High convergence order finite elements with lumped mass matrix
DEFF Research Database (Denmark)
Jensen, Morten skårup
1996-01-01
A method for deriving hexahedral finite elements with lumped mass matrices for three-dimensional problems is presented. These elements meet the theoretical conditions for high order convergence, and two numerical examples based on the three-dimensional scalar wave equation show that this is also...... the case in practice and that their accuracy is comparable to elements with consistent mass matrices....
Parrish, Robert M; Hohenstein, Edward G; Schunck, Nicolas F; Sherrill, C David; Martínez, Todd J
2013-09-27
Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the exact tensor hypercontraction method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly separable tensor factorization of the configuration-space matrix elements. This representation allows for substantial computational savings in chemical, atomic, and nuclear physics simulations, particularly with respect to difficult "exchangelike" contractions.
Accuracy, resolution, and computational complexity of a discontinuous Galerkin finite element method
van der Ven, H.; van der Vegt, Jacobus J.W.; Cockburn, B.; Karniadakis, G.E.; Shu, C.W.
2000-01-01
This series contains monographs of lecture notes type, lecture course material, and high-quality proceedings on topics described by the term "computational science and engineering". This includes theoretical aspects of scientific computing such as mathematical modeling, optimization methods,
A Galerkin Solution for Burgers' Equation Using Cubic B-Spline Finite Elements
Soliman, A. A.
2012-01-01
Numerical solutions for Burgers’ equation based on the Galerkins’ method using cubic B-splines as both weight and interpolation functions are set up. It is shown that this method is capable of solving Burgers’ equation accurately for values of viscosity ranging from very small to large. Three standard problems are used to validate the proposed algorithm. A linear stability analysis shows that a numerical scheme based on a Cranck-Nicolson approximation in time is unconditionally stable.
A Galerkin Solution for Burgers' Equation Using Cubic B-Spline Finite Elements
Directory of Open Access Journals (Sweden)
A. A. Soliman
2012-01-01
Full Text Available Numerical solutions for Burgers’ equation based on the Galerkins’ method using cubic B-splines as both weight and interpolation functions are set up. It is shown that this method is capable of solving Burgers’ equation accurately for values of viscosity ranging from very small to large. Three standard problems are used to validate the proposed algorithm. A linear stability analysis shows that a numerical scheme based on a Cranck-Nicolson approximation in time is unconditionally stable.
2014-05-01
govindah sarva karana-karanam - Brahma Samhita 5.1 “Krsna who is known as Govinda is the Supreme Godhead. He has an eternal blissful spiritual body. He is...density ablators, for meaningful model comparison. Architecture of the material is chosen such that the material can be easily manufactured but is
Pesch, L.; van der Vegt, Jacobus J.W.
2008-01-01
Using the generalized variable formulation of the Euler equations of fluid dynamics, we develop a numerical method that is capable of simulating the flow of fluids with widely differing thermodynamic behavior: ideal and real gases can be treated with the same method as an incompressible fluid. The
Space-time discontinuous Galerkin finite element method for two-fluid flows
Sollie, W.E.H.
2010-01-01
Multifluid and multiphase flows involve combinations of fluids and interfaces which separate these. These flows are of importance in many natural and industrial processes including fluidized beds and bubble columns. Often the interface is not static but moves with the fluid flow velocity. Also,
Pesch, L.
2007-01-01
Many numerical methods for fluid dynamics are suitable only for a single, idealized type of fluid. Most prominently, algorithms for compressible flow are often tailored to ideal gases and another class of schemes is designed for incompressible media. This dissertation targets a numerical method for
Extension of a discontinuous Galerkin finite element method to viscous rotor flow simulations
van der Ven, H.; Boelens, O.J.; Klaij, C.M.; van der Vegt, Jacobus J.W.
2005-01-01
Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, op- erational availability, the pilot’s comfort, and the ef- fectiveness of weapon targeting systems. A precise understanding of the source of
Extension of a discontinuous Galerkin finite element method to viscous rotor flow simulations
van der Ven, H.; Boelens, O.J.; Boelens, O.J.; Klaij, C.M.; van der Vegt, Jacobus J.W.
2005-01-01
Heavy vibratory loading of rotorcraft is relevant for many operational aspects of helicopters, such as the structural life span of (rotating) components, operational availability, the pilot's comfort, and the effectiveness of weapon targeting systems. A precise understanding of the source of these
Finite Element-Galerkin Approximation of the Eigenvalues of Eigenvectors of Selfadjoint Problems
1988-07-01
l’ "k, + 1. Combining (3.20), (3.22), and the fact that I-Eh(Ak ) and Ph are orthogonal projections we have I(I-Eh(Xk,)) PhUB 5 Si (I-Eh(xk)) PhT(Ph-I...Its adjoint are equal. (3.23) implies Hf(I-Eh(1kI )Ph)u{1B - P(IPh)UIBI 5 I(I-Eh(Ak )) PhuB -< d i ii ( Ph- I )T II H B_--H,3 1(P h- I ) u liB , and
Reliable finite element methods for self-adjoint singular perturbation ...
African Journals Online (AJOL)
It is well known that the standard finite element method based on the space Vh of continuous piecewise linear functions is not reliable in solving singular perturbation problems. It is also known that the solution of a two-point boundaryvalue singular perturbation problem admits a decomposition into a regular part and a finite ...
Finite element analyses for RF photoinjector gun cavities
International Nuclear Information System (INIS)
Marhauser, F.
2006-01-01
This paper details electromagnetical, thermal and structural 3D Finite Element Analyses (FEA) for normal conducting RF photoinjector gun cavities. The simulation methods are described extensively. Achieved results are presented. (orig.)
Structural Topology Optimization Based on the Smoothed Finite Element Method
Directory of Open Access Journals (Sweden)
Vahid Shobeiri
Full Text Available Abstract In this paper, the smoothed finite element method, incorporated with the level set method, is employed to carry out the topology optimization of continuum structures. The structural compliance is minimized subject to a constraint on the weight of material used. The cell-based smoothed finite element method is employed to improve the accuracy and stability of the standard finite element method. Several numerical examples are presented to prove the validity and utility of the proposed method. The obtained results are compared with those obtained by several standard finite element-based examples in order to access the applicability and effectiveness of the proposed method. The common numerical instabilities of the structural topology optimization problems such as checkerboard pattern and mesh dependency are studied in the examples.
Comparison of different precondtioners for nonsymmtric finite volume element methods
Energy Technology Data Exchange (ETDEWEB)
Mishev, I.D.
1996-12-31
We consider a few different preconditioners for the linear systems arising from the discretization of 3-D convection-diffusion problems with the finite volume element method. Their theoretical and computational convergence rates are compared and discussed.
Finite element analysis of unnotched charpy impact tests
2008-10-01
This paper describes nonlinear finite element analysis (FEA) to examine the energy to : fracture unnotched Charpy specimens under pendulum impact loading. An oversized, : nonstandard pendulum impactor, called the Bulk Fracture Charpy Machine (BFCM), ...
Finite element analyses of railroad tank car head impacts
2008-09-24
This paper describes engineering analyses of a railroad : tank car impacted at its head by a rigid punch. This type of : collision, referred to as a head impact, is examined using : dynamic, nonlinear finite element analysis (FEA). : Commercial softw...
Finite Element Analysis of the Hierarchical Structure of Human Bone
National Research Council Canada - National Science Library
Dolloff, Katherine
2003-01-01
.... Finally, the effective stiffness of the bone was estimated. In order to determine the stiffness of the collagen fiber, a three-dimensional finite element model was developed and a simple analytical model was derived...
Finite element modelling of helmeted head impact under frontal ...
Indian Academy of Sciences (India)
CSF), brain, tentorium and falx. The finite element model of the helmet consists of shell and foam liner. ... mechanical behaviour of motorcycle helmet. ... the latter authors use a SI (Structural Intensity) approach to study power flow distribution.
Structural analysis with the finite element method linear statics
Oñate, Eugenio
2013-01-01
STRUCTURAL ANALYSIS WITH THE FINITE ELEMENT METHOD Linear Statics Volume 1 : The Basis and Solids Eugenio Oñate The two volumes of this book cover most of the theoretical and computational aspects of the linear static analysis of structures with the Finite Element Method (FEM). The content of the book is based on the lecture notes of a basic course on Structural Analysis with the FEM taught by the author at the Technical University of Catalonia (UPC) in Barcelona, Spain for the last 30 years. Volume1 presents the basis of the FEM for structural analysis and a detailed description of the finite element formulation for axially loaded bars, plane elasticity problems, axisymmetric solids and general three dimensional solids. Each chapter describes the background theory for each structural model considered, details of the finite element formulation and guidelines for the application to structural engineering problems. The book includes a chapter on miscellaneous topics such as treatment of inclined supports, elas...
Finite element analysis of bending performance on polyurethane composite panel
Jia, Minli; Li, Hongqiao; Wang, Xiaoming
2017-09-01
The finite element analysis model of polyurethane composite panel (simply named PCP) is established by using ABAQUS software. In view of the PCPs made of different thickness of surface board, their bending performance is carried out on finite element analysis, and the load-deflection curves which come from it are compared with the experimental results. The results show that the values between finite element analysis and experiment agree well with each other. It can be deduced that the established finite element model is fit to simulate the bending test of PCPs. The simulation not only has certain reference significance to the optimal design for the bending performance of PCPs, but also to the choice of PCPs in the practical project.
Finite-element method for above-core structures
International Nuclear Information System (INIS)
Kennedy, J.M.; Belytschko, T.B.
1979-12-01
Three-dimensional finite-element models for the treatment of the nonlinear, transient response of a fast breeder reactor's above-core structures are described. For purposes of treating arbitrarily large rotations, node orientations are described by unit vectors and the deformable elements are treated by a corotational formulation in which the coordinate system is embedded in the elements. Deformable elements may be connected either to nodes directly or through rigid bodies. The time integration is carried out by the Newmark β method. These features have been incorporated to form the finite-element program SAFE/RAS (Safety Analysis by Finite Elements/Reactor Analysis and Safety Division). Computations are presented for semianalytical comparisons, simple scoping studies, and Stanford Research Institute (SRI) test comparisons
Application of Mass Lumped Higher Order Finite Elements
International Nuclear Information System (INIS)
J. Chen, H.R. Strauss, S.C. Jardin, W. Park, L.E. Sugiyama, G. Fu, J. Breslau
2005-01-01
There are many interesting phenomena in extended-MHD such as anisotropic transport, mhd, 2-fluid effects stellarator and hot particles. Any one of them challenges numerical analysts, and researchers are seeking for higher order methods, such as higher order finite difference, higher order finite elements and hp/spectral elements. It is true that these methods give more accurate solution than their linear counterparts. However, numerically they are prohibitively expensive. Here we give a successful solution of this conflict by applying mass lumped higher order finite elements. This type of elements not only keep second/third order accuracy but also scale closely to linear elements by doing mass lumping. This is especially true for second order lump elements. Full M3D and anisotropic transport models are studied
The finite-difference and finite-element modeling of seismic wave propagation and earthquake motion
International Nuclear Information System (INIS)
Moszo, P.; Kristek, J.; Galis, M.; Pazak, P.; Balazovijech, M.
2006-01-01
Numerical modeling of seismic wave propagation and earthquake motion is an irreplaceable tool in investigation of the Earth's structure, processes in the Earth, and particularly earthquake phenomena. Among various numerical methods, the finite-difference method is the dominant method in the modeling of earthquake motion. Moreover, it is becoming more important in the seismic exploration and structural modeling. At the same time we are convinced that the best time of the finite-difference method in seismology is in the future. This monograph provides tutorial and detailed introduction to the application of the finite-difference, finite-element, and hybrid finite-difference-finite-element methods to the modeling of seismic wave propagation and earthquake motion. The text does not cover all topics and aspects of the methods. We focus on those to which we have contributed. (Author)
Symmetry-preserving finite element schemes: An introductory investigation
Bihlo, Alexander; Valiquette, Francis
2018-01-01
Using the method of equivariant moving frames, we present a procedure for constructing symmetry-preserving finite element methods for second-order ordinary differential equations. Using the method of lines, we then indicate how our constructions can be extended to (1+1)-dimensional evolutionary partial differential equations, using Burgers' equation as an example. Numerical simulations verify that the symmetry-preserving finite element schemes constructed converge at the expected rate and tha...
Examples of finite element mesh generation using SDRC IDEAS
Zapp, John; Volakis, John L.
1990-01-01
IDEAS (Integrated Design Engineering Analysis Software) offers a comprehensive package for mechanical design engineers. Due to its multifaceted capabilities, however, it can be manipulated to serve the needs of electrical engineers, also. IDEAS can be used to perform the following tasks: system modeling, system assembly, kinematics, finite element pre/post processing, finite element solution, system dynamics, drafting, test data analysis, and project relational database.
Non-intrusive finite element reliability analysis methods
Papaioannou, Iason
2014-01-01
This thesis focuses on the modeling of uncertainties in structural systems and on strategies for the reliability assessment of structures analysed by finite element programs. New concepts are introduced for the numerical treatment of spatially varied uncertain quantities through the discretization of the relevant random fields as well as for robust and efficient finite element reliability analysis and updating of the reliability in light of new information. The methods have been implemented i...
Finite element modeling of the filament winding process using ABAQUS
Miltenberger, Louis C.
1992-01-01
A comprehensive stress model of the filament winding fabrication process, previously implemented in the finite element program, WACSAFE, was implemented using the ABAQUS finite element software package. This new implementation, referred to as the ABWACSAFE procedure, consists of the ABAQUS software and a pre/postprocessing routine that was developed to prepare necessary ABAQUS input files and process ABAQUS displacement results for stress and strain computation. The ABWACSAF...
A Comparison of Continuous Mass-lumped Finite Elements and Finite Differences for 3D
Zhebel, E.; Minisini, S.; Kononov, A.; Mulder, W.A.
2012-01-01
The finite-difference method is widely used for time-domain modelling of the wave equation because of its ease of implementation of high-order spatial discretization schemes, parallelization and computational efficiency. However, finite elements on tetrahedral meshes are more accurate in complex
An adaptive discontinuous finite element method for the transport equation
International Nuclear Information System (INIS)
Lang, J.; Walter, A.
1995-01-01
In this paper we introduce a discontinuous finite element method. In our approach, it is possible to combine the advantages of finite element and finite difference methods. The main ingredients are numerical flux approximation and local orthogonal basis functions. The scheme is defined on arbitrary triangulations and two different error indicators are derived. Especially the second one is closely connected to our approach and able to handle arbitrary varying flow directions. Numerical results are given for boundary value problems in two dimensions. They demonstrate the performance of the scheme, combined with the two error indicators
Comparison of 3-D finite elements for incompressible fluid flow
International Nuclear Information System (INIS)
Robichaud, M.; Tanguy, P.A.
1985-01-01
In recent years, the finite element method applied to the solution of incompressible fluid flow has been in constant evolution. In the present state-of-the-art, 2-D problems are solved routinely and reliable results are obtained at a reasonable cost. In 3-D the finite element method is still undergoing active research and many methods have been proposed to solve the Navier-Stokes equations at 'low cost'. These methods have in common the choice of the element which has a trilinear velocity and a discontinuous constant pressure (Q1-PO). The prohibitive cost of 3-D finite element method in fluid flow is the reason for this choice: the Q1-PO is the simplest and the cheapest 3-D element. However, as mentioned in (5) and (6), it generates 'spurious' pressure modes phenomenon called checkerboarding. On regular mesh these spurious modes can be filtered but on distorted mesh the pressure solution is meaningless. (author)
Banks, H T; Birch, Malcolm J; Brewin, Mark P; Greenwald, Stephen E; Hu, Shuhua; Kenz, Zackary R; Kruse, Carola; Maischak, Matthias; Shaw, Simon; Whiteman, John R
2014-04-13
We revisit a method originally introduced by Werder et al. (in Comput. Methods Appl. Mech. Engrg., 190:6685-6708, 2001) for temporally discontinuous Galerkin FEMs applied to a parabolic partial differential equation. In that approach, block systems arise because of the coupling of the spatial systems through inner products of the temporal basis functions. If the spatial finite element space is of dimension D and polynomials of degree r are used in time, the block system has dimension ( r + 1) D and is usually regarded as being too large when r > 1. Werder et al. found that the space-time coupling matrices are diagonalizable over [Formula: see text] for r ⩽ 100, and this means that the time-coupled computations within a time step can actually be decoupled. By using either continuous Galerkin or spectral element methods in space, we apply this DG-in-time methodology, for the first time, to second-order wave equations including elastodynamics with and without Kelvin-Voigt and Maxwell-Zener viscoelasticity. An example set of numerical results is given to demonstrate the favourable effect on error and computational work of the moderately high-order (up to degree 7) temporal and spatio-temporal approximations, and we also touch on an application of this method to an ambitious problem related to the diagnosis of coronary artery disease. Copyright © 2014 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.
Kannan, Kidambi S.; Dasgupta, Abhijit
1998-04-01
Deformation control of smart structures and damage detection in smart composites by magneto-mechanical tagging are just a few of the increasing number of applications of polydomain, polycrystalline magnetostrictive materials that are currently being researched. Robust computational models of bulk magnetostriction will be of great assistance to designers of smart structures for optimization of performance and development of control strategies. This paper discusses the limitations of existing tools, and reports on the work of the authors in developing a 3D nonlinear continuum finite element scheme for magnetostrictive structures, based on an appropriate Galerkin variational principle and incremental constitutive relations. The unique problems posed by the form of the equations governing magneto-mechanical interactions as well as their impact on the proper choice of variational and finite element discretization schemes are discussed. An adaptation of vectorial edge functions for interpolation of magnetic field in hexahedral elements is outlined. The differences between the proposed finite element scheme and available formations are also discussed in this paper. Computational results obtained from the newly proposed scheme will be presented in a future paper.
Numerical experiment on finite element method for matching data
International Nuclear Information System (INIS)
Tokuda, Shinji; Kumakura, Toshimasa; Yoshimura, Koichi.
1993-03-01
Numerical experiments are presented on the finite element method by Pletzer-Dewar for matching data of an ordinary differential equation with regular singular points by using model equation. Matching data play an important role in nonideal MHD stability analysis of a magnetically confined plasma. In the Pletzer-Dewar method, the Frobenius series for the 'big solution', the fundamental solution which is not square-integrable at the regular singular point, is prescribed. The experiments include studies of the convergence rate of the matching data obtained by the finite element method and of the effect on the results of computation by truncating the Frobenius series at finite terms. It is shown from the present study that the finite element method is an effective method for obtaining the matching data with high accuracy. (author)
Vafaeian, B; Le, L H; Tran, T N H T; El-Rich, M; El-Bialy, T; Adeeb, S
2016-05-01
The present study investigated the accuracy of micro-scale finite element modeling for simulating broadband ultrasound propagation in water-saturated trabecular bone-mimicking phantoms. To this end, five commercially manufactured aluminum foam samples as trabecular bone-mimicking phantoms were utilized for ultrasonic immersion through-transmission experiments. Based on micro-computed tomography images of the same physical samples, three-dimensional high-resolution computational samples were generated to be implemented in the micro-scale finite element models. The finite element models employed the standard Galerkin finite element method (FEM) in time domain to simulate the ultrasonic experiments. The numerical simulations did not include energy dissipative mechanisms of ultrasonic attenuation; however, they expectedly simulated reflection, refraction, scattering, and wave mode conversion. The accuracy of the finite element simulations were evaluated by comparing the simulated ultrasonic attenuation and velocity with the experimental data. The maximum and the average relative errors between the experimental and simulated attenuation coefficients in the frequency range of 0.6-1.4 MHz were 17% and 6% respectively. Moreover, the simulations closely predicted the time-of-flight based velocities and the phase velocities of ultrasound with maximum relative errors of 20 m/s and 11 m/s respectively. The results of this study strongly suggest that micro-scale finite element modeling can effectively simulate broadband ultrasound propagation in water-saturated trabecular bone-mimicking structures. Copyright © 2016 Elsevier B.V. All rights reserved.
Finite element analysis of soil-sheet pile interaction
Nyby, D. W.
A finite element model which accurately and economically models soil-sheet pile structures was developed. The model was used to analyze cantilever and anchored sheet pile walls. The finite element model includes transition and interface elements. The transition element has the capability of conforming to the displaced shape of the sheet pile elements on one side (cubic element) and soil elements on the other sides (bilinear element). The interface element models the frictional resistance between the soil and the sheet pile. It behaves elastically below a threshold force level (Coulomb friction) and perfectly plastic above this value. The soil is modeled using nonlinear constitutive relations. These relations are used for both the transition elements and the bilinear elements. The economy of the finite element model was increased in two ways. Closed-form integration was used to reduce the computational effort and an equation solver was used which takes advantage of the banded, symmetric, and positive-definite characteristics of the global stiffness matrix.
Finite Element Aircraft Simulation of Turbulence
1997-02-01
A Simulation of Rotor Blade Element Turbulence (SORBET) model has been : developed for realtime aircraft simulation that accommodates stochastic : turbulence and distributed discrete gusts as a function of the terrain. This : model is applicable to c...
The Cost of Continuity: Performance of Iterative Solvers on Isogeometric Finite Elements
Collier, Nathan
2013-03-19
In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional finite element spaces, and Cp-1 B-splines, which represent maximum continuity We provide theoretical estimates for the increase in cost of the matrix-vector product as well as for the construction and application of black-box preconditioners. We accompany these estimates with numerical results and study their sensitivity to various grid parameters such as element size h and polynomial order of approximation p in addition to the aforementioned continuity of the basis. Finally, we present timing results for a range of preconditioning options for the Laplace problem. We conclude that the matrix-vector product operation is at most 33p2/8 times more expensive for the more continuous space, although for moderately low p, this number is significantly reduced. Moreover, if static condensation is not employed, this number further reduces to at most a value of 8, even for high p. Preconditioning options can be up to p3 times more expensive to set up, although this difference significantly decreases for some popular preconditioners such as incomplete LU factorization. © 2013 Society for Industrial and Applied Mathematics.
A finite element analysis of the distribution velocity in viscous ...
African Journals Online (AJOL)
In this work we use the finite element method to analyze the distribution of velocity in a viscous incompressible fluid flow using Lagrange interpolation function. The results obtained are highly accurate and converge fast to the exact solution as the number of elements increase.
Finite element stress analysis of brick-mortar masonry under ...
African Journals Online (AJOL)
Stress analysis of a brick-mortar couplet as a substitute for brick wall structure has been performed by finite element method, and algorithm for determining the element stiffness matrix for a plane stress problem using the displacement approach was developed. The nodal displacements were derived for the stress in each ...
Finite element solution of the Boussinesq wave equation | Akpobi ...
African Journals Online (AJOL)
In this work, we investigate a Boussinesq-type flow model for nonlinear dispersive waves by developing a computational model based on the finite element discretisation technique. Hermite interpolation functions were used to interpolate approximation elements. The system is modeled using a time dependent equation.
Stress distributions in finite element analysis of concrete gravity dam ...
African Journals Online (AJOL)
Gravity dams are solid structures built of mass concrete material; they maintain their stability against the design loads from the geometric shape, the mass, and the strength of the concrete. The model was meshed with an 8-node biquadratic plane strain quadrilateral (CPE8R) elements, using ABAQUS, a finite element ...
A nonlinear dynamic corotational finite element model for submerged pipes
De Vries, F. H.; Geijselaers, H. J.M.; Van Den Boogaard, A. H.; Huisman, A.
2017-01-01
A three dimensional finite element model is built to compute the motions of a pipe that is being laid on the seabed. This process is geometrically nonlinear, therefore co-rotational beam elements are used. The pipe is subject to static and dynamic forces. Static forces are due to gravity, current
Behaviour of Lagrangian triangular mixed fluid finite elements
Indian Academy of Sciences (India)
The behaviour of mixed fluid finite elements, formulated based on the Lagrangian frame of reference, is investigated to understand the effects of locking due to incompressibility and irrotational constraints. For this purpose, both linear and quadratic mixed triangular fluid elements are formulated. It is found that there exists a ...
Modelling Convergence of Finite Element Analysis of Cantilever Beam
African Journals Online (AJOL)
Convergence studies are carried out by investigating the convergence of numerical results as the number of elements is increased. If convergence is not obtained, the engineer using the finite element method has absolutely no indication whether the results are indicative of a meaningful approximation to the correct solution ...
A set of pathological tests to validate new finite elements
Indian Academy of Sciences (India)
M. Senthilkumar (Newgen Imaging) 1461 1996 Oct 15 13:05:22
End shear. 16. 16. Thick cylinder problem. ∗. Radial pressure. 17. 17. Membrane problem. End shear. 18. 18. Cantilever plate test. Tip moment. 19. End shear. 19 ... Shell finite elements testing. Individual element tests. Patch test. &. FEM convergence. 1989. 1995. 1997. 1998. &. White. & Taylor. Patch test revisited. Zhang.
Behaviour of Lagrangian triangular mixed fluid finite elements
Indian Academy of Sciences (India)
relationship with the penalty finite element approach. Since two constraints are required to be enforced simultaneously, three-field mixed elements involving the displacement and the two Lagrange multipliers, are formulated. The study also includes the effect of bubble functions (or incompatible modes) on the behaviour of ...
Natarajan, Sundararajan; Bordas, Stéphane; Ooi, Ean Tat
2015-01-01
We show both theoretically and numerically a connection between the smoothed finite element method (SFEM) and the virtual element method and use this approach to derive stable, cheap and optimally convergent polyhedral FEM.We show that the stiffness matrix computed with one subcell SFEM is identical to the consistency term of the virtual element method, irrespective of the topology of the element, as long as the shape functions vary linearly on the boundary. Using this connection, we propose ...
Finite element approximation to the even-parity transport equation
International Nuclear Information System (INIS)
Lewis, E.E.
1981-01-01
This paper studies the finite element method, a procedure for reducing partial differential equations to sets of algebraic equations suitable for solution on a digital computer. The differential equation is cast into the form of a variational principle, the resulting domain then subdivided into finite elements. The dependent variable is then approximated by a simple polynomial, and these are linked across inter-element boundaries by continuity conditions. The finite element method is tailored to a variety of transport problems. Angular approximations are formulated, and the extent of ray effect mitigation is examined. Complex trial functions are introduced to enable the inclusion of buckling approximations. The ubiquitous curved interfaces of cell calculations, and coarse mesh methods are also treated. A concluding section discusses limitations of the work to date and suggests possible future directions
Ship Impact Study: Analytical Approaches and Finite Element Modeling
Directory of Open Access Journals (Sweden)
Pawel Woelke
2012-01-01
Full Text Available The current paper presents the results of a ship impact study conducted using various analytical approaches available in the literature with the results obtained from detailed finite element analysis. Considering a typical container vessel impacting a rigid wall with an initial speed of 10 knots, the study investigates the forces imparted on the struck obstacle, the energy dissipated through inelastic deformation, penetration, local deformation patterns, and local failure of the ship elements. The main objective of the paper is to study the accuracy and generality of the predictions of the vessel collision forces, obtained by means of analytical closed-form solutions, in reference to detailed finite element analyses. The results show that significant discrepancies between simplified analytical approaches and detailed finite element analyses can occur, depending on the specific impact scenarios under consideration.
One of the most important considerations in designing a high-speed, gas-lubricated, self acting journal bearing is the instability of the journal under a given operating condition. Intensive research in this direction has led to a number of significant contributions in the past five years. The present paper extends the method of Cheng and Trumpler to study the stability of plain cylindrical gas journal bearings of finite length. Both equilibrium and stability results were obtained.
A wave finite element analysis of the passive cochlea.
Elliott, Stephen J; Ni, Guangjian; Mace, Brian R; Lineton, Ben
2013-03-01
Current models of the cochlea can be characterized as being either based on the assumed propagation of a single slow wave, which provides good insight, or involve the solution of a numerical model, such as in the finite element method, which allows the incorporation of more detailed anatomical features. In this paper it is shown how the wave finite element method can be used to decompose the results of a finite element calculation in terms of wave components, which allows the insight of the wave approach to be brought to bear on more complicated numerical models. In order to illustrate the method, a simple box model is considered, of a passive, locally reacting, basilar membrane interacting via three-dimensional fluid coupling. An analytic formulation of the dispersion equation is used initially to illustrate the types of wave one would expect in such a model. The wave finite element is then used to calculate the wavenumbers of all the waves in the finite element model. It is shown that only a single wave type dominates the response until this peaks at the best place in the cochlea, where an evanescent, higher order fluid wave can make a significant contribution.
Finite element analysis of thrust angle contact ball slewing bearing
Deng, Biao; Guo, Yuan; Zhang, An; Tang, Shengjin
2017-12-01
In view of the large heavy slewing bearing no longer follows the rigid ring hupothesis under the load condition, the entity finite element model of thrust angular contact ball bearing was established by using finite element analysis software ANSYS. The boundary conditions of the model were set according to the actual condition of slewing bearing, the internal stress state of the slewing bearing was obtained by solving and calculation, and the calculated results were compared with the numerical results based on the rigid ring assumption. The results show that more balls are loaded in the result of finite element method, and the maximum contact stresses between the ball and raceway have some reductions. This is because the finite element method considers the ferrule as an elastic body. The ring will produce structure deformation in the radial plane when the heavy load slewing bearings are subjected to external loads. The results of the finite element method are more in line with the actual situation of the slewing bearing in the engineering.
Two-dimensional isostatic meshes in the finite element method
Martínez Marín, Rubén; Samartín, Avelino
2002-01-01
In a Finite Element (FE) analysis of elastic solids several items are usually considered, namely, type and shape of the elements, number of nodes per element, node positions, FE mesh, total number of degrees of freedom (dot) among others. In this paper a method to improve a given FE mesh used for a particular analysis is described. For the improvement criterion different objective functions have been chosen (Total potential energy and Average quadratic error) and the number of nodes and dof's...
Finite element modelling of helmeted head impact under frontal ...
Indian Academy of Sciences (India)
2016-08-26
Aug 26, 2016 ... Finite element models of the head and helmet were used to study contact forces during frontal impact of the head with a rigid surface. The ﬁnite element model of the head consists of skin, skull, cerebro-spinal ﬂuid (CSF), brain, tentorium and falx. The ﬁnite element model of the helmet consists of shell and ...
DEFF Research Database (Denmark)
Marhadi, Kun Saptohartyadi; Evgrafov, Anton; Sørensen, Mads Peter
2011-01-01
We demonstrate the use of a C0 discontinuous Galerkin method for topology optimization of nano-mechanical sensors, namely temperature, surface stress, and mass sensors. The sensors are modeled using classical thin plate theory, which requires C1 basis functions in the standard finite element method...
Error Analysis of p-Version Discontinuous Galerkin Method for Heat Transfer in Built-up Structures
Kaneko, Hideaki; Bey, Kim S.
2004-01-01
The purpose of this paper is to provide an error analysis for the p-version of the discontinuous Galerkin finite element method for heat transfer in built-up structures. As a special case of the results in this paper, a theoretical error estimate for the numerical experiments recently conducted by James Tomey is obtained.
Finite Element Method for Capturing Ultra-relativistic Shocks
Richardson, G. A.; Chung, T. J.
2003-01-01
While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.
Shekar, Balla Chandra; Kishan, Naikoti
2015-12-01
Free convection heat transfer in a square cavity filled with nanofluid-saturated porous medium with the effects of different nanoparticles in the presence of thermal radiation is investigated in this paper. The top and bottom horizontal walls of cavity are considered adiabatic, while the vertical walls are kept at constant temperatures. The governing partial differential equations are solved by finite element method of Galerkin weighted residual scheme. Numerical results are obtained for different values of the Rayleigh number, radiation parameter and nanofluid volume fraction. The overall investigation of variation of streamlines, isotherms and Nusselt numbers is presented graphically. To examine the accuracy, the present results are compared with the available results.
Engineering computation of structures the finite element method
Neto, Maria Augusta; Roseiro, Luis; Cirne, José; Leal, Rogério
2015-01-01
This book presents theories and the main useful techniques of the Finite Element Method (FEM), with an introduction to FEM and many case studies of its use in engineering practice. It supports engineers and students to solve primarily linear problems in mechanical engineering, with a main focus on static and dynamic structural problems. Readers of this text are encouraged to discover the proper relationship between theory and practice, within the finite element method: Practice without theory is blind, but theory without practice is sterile. Beginning with elasticity basic concepts and the classical theories of stressed materials, the work goes on to apply the relationship between forces, displacements, stresses and strains on the process of modeling, simulating and designing engineered technical systems. Chapters discuss the finite element equations for static, eigenvalue analysis, as well as transient analyses. Students and practitioners using commercial FEM software will find this book very helpful. It us...
The Finite Element Numerical Modelling of 3D Magnetotelluric
Directory of Open Access Journals (Sweden)
Ligang Cao
2014-01-01
Full Text Available The ideal numerical simulation of 3D magnetotelluric was restricted by the methodology complexity and the time-consuming calculation. Boundary values, the variation of weighted residual equation, and the hexahedral mesh generation method of finite element are three major causes. A finite element method for 3D magnetotelluric numerical modeling is presented in this paper as a solution for the problem mentioned above. In this algorithm, a hexahedral element coefficient matrix for magnetoelluric finite method is developed, which solves large-scale equations using preconditioned conjugate gradient of the first-type boundary conditions. This algorithm is verified using the homogeneous model, and the positive landform model, as well as the low resistance anomaly model.
Finite Element Modelling of Cold Formed Stainless Steel Columns
Directory of Open Access Journals (Sweden)
M. Macdonald
2005-01-01
Full Text Available This paper describes the results obtained from a finite element investigation into the load capacity of column members of lipped channel cross-section, cold formed from Type 304 stainless steel, subjected to concentric and eccentric compression loading. The main aims of this investigation were to determine the effects which the non-linearity of the stress-strain behaviour of the material would have on the column behaviour under concentric or eccentric loading. Stress-strain curves derived from tests and design codes are incorporated into non-linear finite element analyses of eccentrically loaded columns and the results obtained are compared with those obtained on the basis of experiments on stainless steel channel columns with the same properties and dimensions. Comparisons of the finite element results and the test results are also made with existing design specifications and conclusions are drawn on the basis of the comparisons.
Finite element method for eigenvalue problems in electromagnetics
Reddy, C. J.; Deshpande, Manohar D.; Cockrell, C. R.; Beck, Fred B.
1994-01-01
Finite element method (FEM) has been a very powerful tool to solve many complex problems in electromagnetics. The goal of the current research at the Langley Research Center is to develop a combined FEM/method of moments approach to three-dimensional scattering/radiation problem for objects with arbitrary shape and filled with complex materials. As a first step toward that goal, an exercise is taken to establish the power of FEM, through closed boundary problems. This paper demonstrates the developed of FEM tools for two- and three-dimensional eigenvalue problems in electromagnetics. In section 2, both the scalar and vector finite elements have been used for various waveguide problems to demonstrate the flexibility of FEM. In section 3, vector finite element method has been extended to three-dimensional eigenvalue problems.
Directory of Open Access Journals (Sweden)
P.B. Silva
2013-01-01
Full Text Available Structural spectral elements are formulated using the analytical solution of the applicable elastodynamic equations and, therefore, mesh refinement is not needed to analyze high frequency behavior provided the elastodynamic equations used remain valid. However, for modeling complex structures, standard spectral elements require long and cumbersome analytical formulation. In this work, a method to build spectral finite elements from a finite element model of a slice of a structural waveguide (a structure with one dimension much larger than the other two is proposed. First, the transfer matrix of the structural waveguide is obtained from the finite element model of a thin slice. Then, the wavenumbers and wave propagation modes are obtained from the transfer matrix and used to build the spectral element matrix. These spectral elements can be used to model homogeneous waveguides with constant cross section over long spans without the need of refining the finite element mesh along the waveguide. As an illustrating example, spectral elements are derived for straight uniform rods and beams and used to calculate the forced response in the longitudinal and transverse directions. Results obtained with the spectral element formulation are shown to agree well with results obtained with a finite element model of the whole beam. The proposed approach can be used to generate spectral elements of waveguides of arbitrary cross section and, potentially, of arbitrary order.
A least-squares finite-element Sn method for solving first-order neutron transport equation
International Nuclear Information System (INIS)
Ju Haitao; Wu Hongchun; Zhou Yongqiang; Cao Liangzhi; Yao Dong; Xian, Chun-Yu
2007-01-01
A discrete ordinates finite-element method for solving the two-dimensional first-order neutron transport equation is derived using the least-squares variation. It avoids the singularity in void regions of the method derived from the second-order equation which contains the inversion of the cross-section. Different from using the standard Galerkin variation to the first-order equation, the least-squares variation results in a symmetric matrix, which can be solved easily and effectively. To eliminate the discontinuity of the angular flux on the vacuum boundary in the spherical harmonics method, the angle variable is discretized by the discrete ordinates method. A two-dimensional transport simulation code is developed and applied to some benchmark problems with unstructured geometry. The numerical results verified the validity of this method
A Mixed Finite Element Method to Solve the EEG Forward Problem.
Vorwerk, J; Engwer, C; Pursiainen, S; Wolters, C H
2017-04-01
Finite element methods have been shown to achieve high accuracies in numerically solving the EEG forward problem and they enable the realistic modeling of complex geometries and important conductive features such as anisotropic conductivities. To date, most of the presented approaches rely on the same underlying formulation, the continuous Galerkin (CG)-FEM. In this article, a novel approach to solve the EEG forward problem based on a mixed finite element method (Mixed-FEM) is introduced. To obtain the Mixed-FEM formulation, the electric current is introduced as an additional unknown besides the electric potential. As a consequence of this derivation, the Mixed-FEM is, by construction, current preserving, in contrast to the CG-FEM. Consequently, a higher simulation accuracy can be achieved in certain scenarios, e.g., when the diameter of thin insulating structures, such as the skull, is in the range of the mesh resolution. A theoretical derivation of the Mixed-FEM approach for EEG forward simulations is presented, and the algorithms implemented for solving the resulting equation systems are described. Subsequently, first evaluations in both sphere and realistic head models are presented, and the results are compared to previously introduced CG-FEM approaches. Additional visualizations are shown to illustrate the current preserving property of the Mixed-FEM. Based on these results, it is concluded that the newly presented Mixed-FEM can at least complement and in some scenarios even outperform the established CG-FEM approaches, which motivates a further evaluation of the Mixed-FEM for applications in bioelectromagnetism.
Generalized Multiscale Finite Element Methods for Wave Propagation in Heterogeneous Media
Chung, Eric T.
2014-11-13
Numerical modeling of wave propagation in heterogeneous media is important in many applications. Due to their complex nature, direct numerical simulations on the fine grid are prohibitively expensive. It is therefore important to develop efficient and accurate methods that allow the use of coarse grids. In this paper, we present a multiscale finite element method for wave propagation on a coarse grid. The proposed method is based on the generalized multiscale finite element method (GMsFEM) (see [Y. Efendiev, J. Galvis, and T. Hou, J. Comput. Phys., 251 (2012), pp. 116--135]). To construct multiscale basis functions, we start with two snapshot spaces in each coarse-grid block, where one represents the degrees of freedom on the boundary and the other represents the degrees of freedom in the interior. We use local spectral problems to identify important modes in each snapshot space. These local spectral problems are different from each other and their formulations are based on the analysis. To the best of knowledge, this is the first time that multiple snapshot spaces and multiple spectral problems are used and necessary for efficient computations. Using the dominant modes from local spectral problems, multiscale basis functions are constructed to represent the solution space locally within each coarse block. These multiscale basis functions are coupled via the symmetric interior penalty discontinuous Galerkin method which provides a block diagonal mass matrix and, consequently, results in fast computations in an explicit time discretization. Our methods\\' stability and spectral convergence are rigorously analyzed. Numerical examples are presented to show our methods\\' performance. We also test oversampling strategies. In particular, we discuss how the modes from different snapshot spaces can affect the proposed methods\\' accuracy.
AEROTAXI ground static test and finite element model validation
Directory of Open Access Journals (Sweden)
Radu BISCA
2011-06-01
Full Text Available In this presentation, we will concentrate on typical Ground Static Test (GST and Finite Element (FE software comparisons. It is necessary to note, that standard GST are obligatory for any new aircraft configuration. We can mention here the investigations of the AeroTAXITM, a small aircraft configuration, using PRODERA® equipment. A Finite Element Model (FEM of the AeroTAXITM has been developed in PATRAN/NASTRAN®, partly from a previous ANSYS® model. FEM can be used to investigate potential structural modifications or changes with realistic component corrections. Model validation should be part of every modern engineering analysis and quality assurance procedure.
Stochastic Finite Elements in Reliability-Based Structural Optimization
DEFF Research Database (Denmark)
Sørensen, John Dalsgaard; Engelund, S.
1995-01-01
Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect to optimi......Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect...... to optimization variables can be performed. A computer implementation is described and an illustrative example is given....
Fourier analysis of finite element preconditioned collocation schemes
Deville, Michel O.; Mund, Ernest H.
1990-01-01
The spectrum of the iteration operator of some finite element preconditioned Fourier collocation schemes is investigated. The first part of the paper analyses one-dimensional elliptic and hyperbolic model problems and the advection-diffusion equation. Analytical expressions of the eigenvalues are obtained with use of symbolic computation. The second part of the paper considers the set of one-dimensional differential equations resulting from Fourier analysis (in the tranverse direction) of the 2-D Stokes problem. All results agree with previous conclusions on the numerical efficiency of finite element preconditioning schemes.
Stochastic Finite Elements in Reliability-Based Structural Optimization
DEFF Research Database (Denmark)
Sørensen, John Dalsgaard; Engelund, S.
Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect to optimi......Application of stochastic finite elements in structural optimization is considered. It is shown how stochastic fields modelling e.g. the modulus of elasticity can be discretized in stochastic variables and how a sensitivity analysis of the reliability of a structural system with respect...
Preconditioning for Mixed Finite Element Formulations of Elliptic Problems
Wildey, Tim
2013-01-01
In this paper, we discuss a preconditioning technique for mixed finite element discretizations of elliptic equations. The technique is based on a block-diagonal approximation of the mass matrix which maintains the sparsity and positive definiteness of the corresponding Schur complement. This preconditioner arises from the multipoint flux mixed finite element method and is robust with respect to mesh size and is better conditioned for full permeability tensors than a preconditioner based on a diagonal approximation of the mass matrix. © Springer-Verlag Berlin Heidelberg 2013.
Diffusive mesh relaxation in ALE finite element numerical simulations
Energy Technology Data Exchange (ETDEWEB)
Dube, E.I.
1996-06-01
The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.
FINITE ELEMENT EVALUATION AND OPTIMIZATION OF GEOMETRY WITH DOE
Directory of Open Access Journals (Sweden)
Janko D. Jovanovic
2011-03-01
Full Text Available Since 1960, Taguchi methods have been used for improving the quality of Japanese products with great success. Basic assumption of Taguchi's design for six sigma or robust design is that quality must be designed into a product from the start at both the product and process design stage in order to improve product reliability and manufacturability. This paper deals with case study of product design based on Taguchi's approach that involves parametric optimization of piston rod geometry aiming mass reduction with stress restriction. Finite element analysis software ANSYS Workbench was used to get access to CAD parameters of piston rod within a process of parametric finite element evaluation and optimization.
FINITE ELEMENT MODELING OF THIN CIRCULAR SANDWICH PLATES DEFLECTION
Directory of Open Access Journals (Sweden)
K. S. Kurachka
2014-01-01
Full Text Available A mathematical model of a thin circular sandwich plate being under the vertical load is proposed. The model employs the finite element method and takes advantage of an axisymmetric finite element that leads to the small dimension of the resulting stiffness matrix and sufficient accuracy for practical calculations. The analytical expressions for computing local stiffness matrices are found, which can significantly speed up the process of forming the global stiffness matrix and increase the accuracy of calculations. A software is under development and verification. The discrepancy between the results of the mathematical model and those of analytical formulas for homogeneous thin circularsandwich plates does not exceed 7%.
COMPUTER EXPERIMENTS WITH FINITE ELEMENTS OF HIGHER ORDER
Directory of Open Access Journals (Sweden)
Khomchenko A.
2017-12-01
Full Text Available The paper deals with the problem of constructing the basic functions of a quadrilateral finite element of the fifth order by the means of the computer algebra system Maple. The Lagrangian approximation of such a finite element contains 36 nodes: 20 nodes perimeter and 16 internal nodes. Alternative models with reduced number of internal nodes are considered. Graphs of basic functions and cognitive portraits of lines of zero level are presented. The work is aimed at studying the possibilities of using modern information technologies in the teaching of individual mathematical disciplines.
A multiscale mortar multipoint flux mixed finite element method
Wheeler, Mary Fanett
2012-02-03
In this paper, we develop a multiscale mortar multipoint flux mixed finite element method for second order elliptic problems. The equations in the coarse elements (or subdomains) are discretized on a fine grid scale by a multipoint flux mixed finite element method that reduces to cell-centered finite differences on irregular grids. The subdomain grids do not have to match across the interfaces. Continuity of flux between coarse elements is imposed via a mortar finite element space on a coarse grid scale. With an appropriate choice of polynomial degree of the mortar space, we derive optimal order convergence on the fine scale for both the multiscale pressure and velocity, as well as the coarse scale mortar pressure. Some superconvergence results are also derived. The algebraic system is reduced via a non-overlapping domain decomposition to a coarse scale mortar interface problem that is solved using a multiscale flux basis. Numerical experiments are presented to confirm the theory and illustrate the efficiency and flexibility of the method. © EDP Sciences, SMAI, 2012.
[Finite element analysis on stress change of lumbar spine].
Yan, Jia-zhi; Wu, Zhi-hong; Wang, Xue-song; Xing, Ze-jun; Song, Hai-feng; Zhao, Yu; Zhang, Jian-guo; Wang, Yi-peng; Qiu, Gui-xing
2009-05-05
To build a 3D finite element model of whole lumbar spine and verify its efficiency and analyze the biomechanical change of L3-4 motion segment. L1-L5 segment data were obtained from computed tomography (CT) scans of the lumbar spine of a 40-year-old man with no abnormal findings. A three-dimensional finite element model of the human whole lumbar spine was built in the Mimics and the ABAQUS software. The model was composed of bony vertebrae, articulating facets, intervertebral disc and associated ligaments. The basic stress analysis of L3-4 motion segment was made under the considerations of different material properties of bone, ligaments and facet joints contacting frictional property. The stress on annulus fiber, nucleus pulposus, endplate and facet joints under axial pressure (0.3 MPa, 0.5 MPa, 1.0 MPa, 2.0 MPa & 4.0 MPa) were analyzed. A three-dimensional finite element model of human L3-L4 motion segment has 272, 619 elements, the stresses were higher in the posterior of annulus fiber, the Max pressure stress (S33) distributed in nucleus pulposus and the center of endplate. The stresses increased as axial pressure rose. 3D finite element model of whole lumbar spine and L3-4 motion segment were established successfully and the stress analyses were feasible and reliable.
Comparison study of finite element and basis set methods for finite size scaling
International Nuclear Information System (INIS)
Antillon, Edwin; Moy, Winton; Wei Qi; Kais, Sabre
2009-01-01
We compare two methods of obtaining critical parameters for a quantum Hamiltonian using a finite size scaling approach. A finite element and basis set method were used in conjunction with the finite size scaling to obtain the critical parameters for the Hulthen potential. The critical parameters obtained analytically were the coupling constant λ c =(1/2), the critical exponents for the energy α=2 and for the 'correlation length 'ν=1. The extrapolated results for finite size scaling with the basis set method are λ c =0.499 99, α=1.9960, and ν=0.999 10. The results for the finite element solutions are λ c =0.501 84, α=1.999 93, and ν=1.000 79 for the linear interpolation and λ c =0.500 00, α=2.000 11, and ν=1.000 32 for the Hermite interpolation. The results for each method compare very well with the analytical results obtained for the Hulthen potential. However, the finite element method is easier to implement and may be combined with ab initio and density functional theory to obtain quantum critical parameters for more complex systems.
Chung, Eric
2015-12-11
In this paper, we develop a mass conservative multiscale method for coupled flow and transport in heterogeneous porous media. We consider a coupled system consisting of a convection-dominated transport equation and a flow equation. We construct a coarse grid solver based on the Generalized Multiscale Finite Element Method (GMsFEM) for a coupled system. In particular, multiscale basis functions are constructed based on some snapshot spaces for the pressure and the concentration equations and some local spectral decompositions in the snapshot spaces. The resulting approach uses a few multiscale basis functions in each coarse block (for both the pressure and the concentration) to solve the coupled system. We use the mixed framework, which allows mass conservation. Our main contributions are: (1) the development of a mass conservative GMsFEM for the coupled flow and transport; (2) the development of a robust multiscale method for convection-dominated transport problems by choosing appropriate test and trial spaces within Petrov-Galerkin mixed formulation. We present numerical results and consider several heterogeneous permeability fields. Our numerical results show that with only a few basis functions per coarse block, we can achieve a good approximation.
Analysis of elastic-plastic problems using edge-based smoothed finite element method
International Nuclear Information System (INIS)
Cui, X.Y.; Liu, G.R.; Li, G.Y.; Zhang, G.Y.; Sun, G.Y.
2009-01-01
In this paper, an edge-based smoothed finite element method (ES-FEM) is formulated for stress field determination of elastic-plastic problems using triangular meshes, in which smoothing domains associated with the edges of the triangles are used for smoothing operations to improve the accuracy and the convergence rate of the method. The smoothed Galerkin weak form is adopted to obtain the discretized system equations, and the numerical integration becomes a simple summation over the edge-based smoothing domains. The pseudo-elastic method is employed for the determination of stress field and Hencky's total deformation theory is used to define effective elastic material parameters, which are treated as field variables and considered as functions of the final state of stress fields. The effective elastic material parameters are then obtained in an iterative manner based on the strain controlled projection method from the uniaxial material curve. Some numerical examples are investigated and excellent results have been obtained demonstrating the effectivity of the present method.
A generalized multiscale finite element method for elastic wave propagation in fractured media
Chung, Eric T.
2016-02-26
In this paper, we consider elastic wave propagation in fractured media applying a linear-slip model to represent the effects of fractures on the wavefield. Fractured media, typically, are highly heterogeneous due to multiple length scales. Direct numerical simulations for wave propagation in highly heterogeneous fractured media can be computationally expensive and require some type of model reduction. We develop a multiscale model reduction technique that captures the complex nature of the media (heterogeneities and fractures) in the coarse scale system. The proposed method is based on the generalized multiscale finite element method, where the multiscale basis functions are constructed to capture the fine-scale information of the heterogeneous, fractured media and effectively reduce the degrees of freedom. These multiscale basis functions are coupled via the interior penalty discontinuous Galerkin method, which provides a block-diagonal mass matrix. The latter is needed for fast computation in an explicit time discretization, which is used in our simulations. Numerical results are presented to show the performance of the presented multiscale method for fractured media. We consider several cases where fractured media contain fractures of multiple lengths. Our numerical results show that the proposed reduced-order models can provide accurate approximations for the fine-scale solution.
Divergence of finite element formulations for inverse problems treated as optimization problems
International Nuclear Information System (INIS)
Rivas, Carlos; Barbone, Paul; Oberai, Assad
2008-01-01
Many inverse problems are formulated and solved as optimization problems. In this approach, the data mismatch between a predicted field and a measured field is minimized, subject to a constraint. The constraint represents the 'forward' model of the system under consideration. In this paper, the model considered is plane stress incompressible elasticity. This pde is discretized using several standard Galerkin finite element methods. These are known to yield stable and convergent discrete solutions that converge with mesh refinement to the exact solution of the forward problem. It is usually taken for granted that if the constraint equation is discretized by a stable, convergent numerical method, then the inverse problem will also converge to the exact solution with mesh refinement. We show examples in this paper, however, where this is not the case. These are based on inverse problems with interior data, which have provably unique solutions. Even so, the use of classical discretization techniques for the forward constraint within the optimization formulation leads to ill-posed discrete problems. We analyze the discrete systems of equations and show the source of the instability. We discuss variational properties of the continuous inverse optimization problem, and describe a novel B-spline FEM to solve it. We present computational evidence that suggests the B-spline FEM inverse problem solution converges to the exact inverse problem solution with mesh refinement.
Fluid structure interaction in electrohydraulic servovalve: a finite element approach
Hiremath, Somashekhar S.; Singaperumal, M.
2010-01-01
Electrohydraulic servovalves (EHSV) promise unique application opportunities and high performance, unmatched by other drive technologies. Typical applications include aerospace, robotic manipulators, motion simulators, injection molding, CNC machines and material testing machines. EHSV available are either a flapper/nozzle type or a jet pipe type. In the present paper an attempt has been made to study the dynamics of jet pipe EHSV with built-in mechanical feedback using Finite Element Method (FEM). In jet pipe EHSV, the dynamics of spool greatly depends on pressure recovery and hence the fluid flow at spool ends. The effect of pressure recovery on spool dynamics is studied using FEM by creating the fluid-structure-interaction. The mechanical parts were created using general purpose finite elements like shell, beam, and solid elements while fluid cavities were created using hydrostatic fluid elements. The analysis was carried out using the commercially available FE code ABAQUS. The jet pipe and spool dynamics are presented in the paper.
Magnetoelastic energy calculations for finite element analysis of superconductors
International Nuclear Information System (INIS)
Akin, J.E.; Stoddart, W.C.T.
1977-01-01
It has been shown that the high current density and magnetic flux density associated with superconductors can make the magnetoelastic energy a significant portion of the total energy in a structural system. The present work presents a procedure for evaluating this magnetoelastic energy for use in the finite element analysis of the structural dynamics and stability of the superconductor. A simple, special case of the element matrices is illustrated
A finite element field solver for dipole modes
International Nuclear Information System (INIS)
Nelson, E.M.
1992-01-01
A finite element field solver for dipole modes in axisymmetric structures has been written. The second-order elements used in this formulation yield accurate mode frequencies with no spurious modes. Quasi-periodic boundaries are included to allow travelling waves in periodic structures. The solver is useful in applications requiring precise frequency calculations such as detuned accelerator structures for linear colliders. Comparisons are made with measurements and with the popular but less accurate field solver URMEL. (author). 7 refs., 4 figs
Choice of input fields in stochastic finite elements
DEFF Research Database (Denmark)
Ditlevsen, Ove Dalager; Tarp-Johansen, Niels Jacob
1996-01-01
variables. Several reported discretization methods define these random variables as integrals of the product of the field and some suitable weight functions. In particular, the weight functions can be Dirac delta functions whereby the random variables become the field values at a finite set of given points......The problem of the arbitrary choice of variables for random field modelling in structural mechanics or in soil mechanics is treated. For example, it is relevant to ask the question of whether it is best to choose a stiffness field along a beam element or to choose its reciprocal field......, the flexibility field, as the input to the stochastic finite element model. To answer this question the focus should be on the error of the output of the mechanical model rather than on the input field itself when discretizing the field through replacing it by a field defined in terms of a finite number of random...
A particle finite element method for machining simulations
Sabel, Matthias; Sator, Christian; Müller, Ralf
2014-07-01
The particle finite element method (PFEM) appears to be a convenient technique for machining simulations, since the geometry and topology of the problem can undergo severe changes. In this work, a short outline of the PFEM-algorithm is given, which is followed by a detailed description of the involved operations. The -shape method, which is used to track the topology, is explained and tested by a simple example. Also the kinematics and a suitable finite element formulation are introduced. To validate the method simple settings without topological changes are considered and compared to the standard finite element method for large deformations. To examine the performance of the method, when dealing with separating material, a tensile loading is applied to a notched plate. This investigation includes a numerical analysis of the different meshing parameters, and the numerical convergence is studied. With regard to the cutting simulation it is found that only a sufficiently large number of particles (and thus a rather fine finite element discretisation) leads to converged results of process parameters, such as the cutting force.
Aranha: a 2D mesh generator for triangular finite elements
International Nuclear Information System (INIS)
Fancello, E.A.; Salgado, A.C.; Feijoo, R.A.
1990-01-01
A method for generating unstructured meshes for linear and quadratic triangular finite elements is described in this paper. Some topics on the C language data structure used in the development of the program Aranha are also presented. The applicability for adaptive remeshing is shown and finally several examples are included to illustrate the performance of the method in irregular connected planar domains. (author)
Finite Element Vibration and Dynamic Response Analysis of Engineering Structures
Directory of Open Access Journals (Sweden)
Jaroslav Mackerle
2000-01-01
Full Text Available This bibliography lists references to papers, conference proceedings, and theses/dissertations dealing with finite element vibration and dynamic response analysis of engineering structures that were published from 1994 to 1998. It contains 539 citations. The following types of structures are included: basic structural systems; ground structures; ocean and coastal structures; mobile structures; and containment structures.
A Finite Element Approach to Modeling Abrasive Wear Modes
Woldman, M.; van der Heide, Emile; Tinga, Tiedo; Masen, Marc Arthur
2016-01-01
Machine components operating in sandy environments will wear because of the abrasive interaction with sand particles. In this work, a method is derived to predict the amount of wear caused by such abrasive action, in order to improve the maintenance concept of the components. A finite element model
Finite element concept to derive isostatic residual maps ...
Indian Academy of Sciences (India)
lies are isolated so as to construct the isostatic residual maps. Very accurate geophysical studies have ... Finite element concept; isostatic anomaly; Gorda Plate; Sierra Nevada. Proc. Indian Acad. Sci. (Earth Planet. Sci.), 110 .... The continuous line is obtained by regression analysis. The broken line shows the. FEA regional ...
An Orthogonal Residual Procedure for Nonlinear Finite Element Equations
DEFF Research Database (Denmark)
Krenk, S.
A general and robust solution procedure for nonlinear finite element equations with limit points is developed. At each equilibrium iteration the magnitude of the load is adjusted such that the residual force is orthogonal to the current displacement increment from the last equilibrium state...
Finite Element Analysis of Boron Diffusion in Wooden Poles
DEFF Research Database (Denmark)
Krabbenhøft, Kristian; Hoffmeyer, P.; Bechgaard, C.
2003-01-01
The problem of describing the migration of dissolved boron in wood is treated with special reference to the commonly used remedial treatment of wooden poles. The governing equations are derived and discussed together with some of the material parameters required. The equations are solved by the f...... by the finite element method and, finally, results showing the effect of different treatment strategies are presented....
2-D Finite Element Analysis of Massive RC Structures
DEFF Research Database (Denmark)
Saabye Ottosen, Niels
1982-01-01
Nonlinear analysis of concrete structures using finite elements is discussed. The applications include a thick-walled top-closure for a pressure vessel as well as the delicate problems of beams failing in shear. The top-closure analysis evaluates the effect of two different failure criteria...
A 2-dimensional finite element simulation of cooling in castings ...
African Journals Online (AJOL)
In this work we present a 2 dimensional finite element simulation of the cooling process in castings. A one way coupling +technique was used to predict the behavior of thermal strains and stresses from the temperature history of casting. The temperature distribution across the casting at different times, the cooling pattern of ...
Finite element analysis of boron diffusion in wooden Poles
DEFF Research Database (Denmark)
Krabbenhøft, Kristian; Hoffmeyer, Preben; Bechgaard, Carl
2004-01-01
The problem of describing the migration of dissolved boron in wood is treated with special reference to the commonly used remedial treatment of wooden poles. The governing equations are derived and discussed together with some of the material parameters required. The equations are solved by the f...... by the finite element method and, finally, results showing the effect of different treatment strategies are presented....
The Development of Piezoelectric Accelerometers Using Finite Element Analysis
DEFF Research Database (Denmark)
Liu, Bin
1999-01-01
This paper describes the application of Finite Element (FE) approach for the development of piezoelectric accelerometers. An accelerometer is simulated using the FE approach as an example. Good agreement is achieved between simulated results and calibrated results. It is proved that the FE modeling...
Efficient implicit finite element analysis of sheet forming processes
van den Boogaard, Antonius H.; Meinders, Vincent T.; Huetink, Han
2003-01-01
The computation time for implicit finite element analyses tends to increase disproportionally with increasing problem size. This is due to the repeated solution of linear sets of equations, if direct solvers are used. By using iterative linear equation solvers the total analysis time can be reduced
Bending analysis of laminated composite plates using finite element ...
African Journals Online (AJOL)
In this paper, a number of finite element analyses have been carried out for various side-to-thickness ratios, aspect ratios and modulus ratios to study the effect of transverse shear deformation on deflection and stresses of laminated composite plates subjected to uniformly distributed load. The numerical results showed, ...
Finite Element Analysis of Boron Diffusion in Wooden Poles
DEFF Research Database (Denmark)
Krabbenhøft, Kristian; Hoffmeyer, Preben; Bechgaard, Carl
2004-01-01
The problem of describing the migration of dissolved boron in wood is treated with special reference to the commonly used remedial treatment of wooden poles. The governing equations are derived and discussed together with some of the material parameters required. The equations are solved by the f...... by the finite element method and, finally, results showing the effect of different treatment strategies are presented....
Finite element simulation of laser transmission welding of dissimilar ...
African Journals Online (AJOL)
Now-a-days, metal to plastic micro-welding is of great interest in the field of biomedical and electronics applications. Laser transmission welding (LTW) has emerged as the most suitable technique for such applications. In this paper, a three-dimensional finite element (FE) thermal model is developed to simulate the laser ...
Assessment of Finite Element Approximations for Nonlinear Flexible Multibody Dynamics
1991-05-01
dynamics. Two nonlinear beam finite elements are consistently derived from virtual work principle using Bernoulli Euler and Timoshenko beam...and dynamic buckling. Equations of motion are derived for rigid central body with flexible appendage using virtual work principle. Virtual work principle
Design, development and use of the finite element machine
Adams, L. M.; Voigt, R. C.
1983-01-01
Some of the considerations that went into the design of the Finite Element Machine, a research asynchronous parallel computer are described. The present status of the system is also discussed along with some indication of the type of results that were obtained.
Finite element modeling of intermuscular interactions and myofascial force transmission
Yucesoy, C.A.; Koopman, Hubertus F.J.M.; Huijing, P.A.J.B.M.; Grootenboer, H.J.
2001-01-01
A finite element muscle model to study the principles of intermuscular myofascial force transmission is developed. The results obtained explain force differences at the distal and proximal tendons of muscles that have mechanical interaction. This is in agreement with experimental findings in other
Finite element analysis of thermoelastic instability with intermittent contact
Geijselaers, Hubertus J.M.; Koning, A.J.E.
2000-01-01
The equations that describe the development of corrugations on block braked wheel treads caused by thermoelastic instability are discretized using the finite element method. The perturbations of temperatures and distortions are described by an amplitude function, which is spatially fixed multiplied
Finite element and perturbative study of buffered leaky planar waveguides
Uranus, H.P.; Hoekstra, Hugo; van Groesen, Embrecht W.C.
2005-01-01
The effects of the presence of a high-index medium in the proximity of planar waveguiding structures that makes up buffered leaky waveguides, were studied using a finite element method (FEM) leaky mode solver and a perturbation method. Various phenomena observed in the FEM results were interpreted
Hands on applied finite element analysis application with ANSYS
Arslan, Mehmet Ali
2015-01-01
Hands on Applied Finite Element Analysis Application with Ansys is truly an extraordinary book that offers practical ways of tackling FEA problems in machine design and analysis. In this book, 35 good selection of example problems have been presented, offering students the opportunity to apply their knowledge to real engineering FEA problem solutions by guiding them with real life hands on experience.
The future of the finite element method in geotechnics
Brinkgreve, R.B.J.
2012-01-01
In this presentation a vision is given on tlie fiiture of the finite element method (FEM) for geotechnical engineering and design. In the past 20 years the FEM has proven to be a powerful method for estimating deformation, stability and groundwater flow in geoteclmical stmctures. Much has been
A direct implementation for influence lines in finite element software
DEFF Research Database (Denmark)
Jepsen, Michael S.; Damkilde, Lars
2014-01-01
The use of influence lines is a recognized method for determining the critical design load conditions and this paper shows a direct method for applying influence lines in any structural finite element software. The main idea is to equate displacement or angular discontinuities with nodal forces...
Finite element modelling of fibre-reinforced brittle materials
Kullaa, J.
1997-01-01
The tensile constitutive behaviour of fibre-reinforced brittle materials can be extended to two or three dimensions by using the finite element method with crack models. The three approaches in this study include the smeared and discrete crack concepts and a multi-surface plasticity model. The
Finite element analysis of tubular joints in offshore structures ...
African Journals Online (AJOL)
... representing a 2-D model of the joint between the brace and the chord walls. This was subsequently followed but finite element analysis of six tubular joints. A global analysis was initially undertaken, then the submodel analysis carried in the areas of stress concentration. Journal of Civil Engineering, JKUAT (2001) Vol 6, ...
Piezoelectric Accelerometers Modification Based on the Finite Element Method
DEFF Research Database (Denmark)
Liu, Bin; Kriegbaum, B.
2000-01-01
The paper describes the modification of piezoelectric accelerometers using a Finite Element (FE) method. Brüel & Kjær Accelerometer Type 8325 is chosen as an example to illustrate the advanced accelerometer development procedure. The deviation between the measurement and FE simulation results...
A Monte Carlo adapted finite element method for dislocation ...
Indian Academy of Sciences (India)
P Zakian
2017-10-10
Oct 10, 2017 ... simulations are proposed. Various comparisons are examined to illustrate the capability of both methods for random simulation of faults. Keywords. Monte Carlo simulation; stochastic modeling; split node technique; finite element method; earthquake fault dislocation. 1. Introduction. In material science, a ...
Reliability-Based Shape Optimization using Stochastic Finite Element Methods
DEFF Research Database (Denmark)
Enevoldsen, Ib; Sørensen, John Dalsgaard; Sigurdsson, G.
1991-01-01
stochastic fields (e.g. loads and material parameters such as Young's modulus and the Poisson ratio). In this case stochastic finite element techniques combined with FORM analysis can be used to obtain measures of the reliability of the structural systems, see Der Kiureghian & Ke (6) and Liu & Der Kiureghian...
Finite Element Analysis of a Free-Standing Staircase | Ajagbe ...
African Journals Online (AJOL)
The existing approximate analytical methods of analyzing free-standing stairs fail to predict the distribution of any stress resultant and the actual three dimensional behavior of the stair slab system. A more rationale but simple and accurate method of analysis based on finite element method is presented. Plate flexural ...
Finite element modelling of elastic intraplate stresses due to ...
Indian Academy of Sciences (India)
Finite element modelling of elastic intraplate stresses due to heterogeneities in crustal density and mechanical properties for the Jabalpur earthquake region, central India. A Manglik1,∗. , S Thiagarajan. 1. , A V Mikhailova. 2 and Yu Rebetsky. 2. 1. National Geophysical Research Institute, Uppal Road, Hyderabad 500 007, ...
Finite element concept to derive isostatic residual maps ...
Indian Academy of Sciences (India)
A new space-domain operator based on the shape function concept of finite element analysis has been developed to derive the residual maps of the Gorda Plate of western United States. The technique does not require explicit assumptions on isostatic models. Besides delineating the Gorda Plate boundary, the residual ...
Superconvergence for tetrahedral quadratic finite element methods for elliptic equations
Brandts, J.H.; Krizek, M.
2005-01-01
For a model elliptic boundary value problem we will prove that on strongly regular families of uniform tetrahedral partitions of the domain, the gradient of the quadratic finite element approximation is superclose to the gradient of the quadratic Lagrange interpolant of the exact solution. This
A mixed finite element method for particle simulation in lasertron
International Nuclear Information System (INIS)
Le Meur, G.
1987-03-01
A particle simulation code is being developed with the aim to treat the motion of charged particles in electromagnetic devices, such as Lasertron. The paper describes the use of mixed finite element methods in computing the field components, without derivating them from scalar or vector potentials. Graphical results are shown
Finite element model to study two dimensional unsteady state ...
African Journals Online (AJOL)
Kunal Pathak
2015-10-20
Oct 20, 2015 ... Excess buffer;. Finite element method. Abstract The calcium signaling plays a crucial role in expansion and contraction of cardiac myo- cytes. This calcium signaling is achieved by calcium diffusion, buffering mechanisms and influx in cardiac myocytes. The various calcium distribution patterns required for ...
Finite element simulations of two rock mechanics tests
International Nuclear Information System (INIS)
Dahlke, H.J.; Lott, S.A.
1986-04-01
Rock mechanics tests are performed to determine in situ stress conditions and material properties of an underground rock mass. To design stable underground facilities for the permanent storage of high-level nuclear waste, determination of these properties and conditions is a necessary first step. However, before a test and its associated equipment can be designed, the engineer needs to know the range of expected values to be measured by the instruments. Sensitivity studies by means of finite element simulations are employed in this preliminary design phase to evaluate the pertinent parameters and their effects on the proposed measurements. The simulations, of two typical rock mechanics tests, the plate bearing test and the flat-jack test, by means of the finite element analysis, are described. The plate bearing test is used to determine the rock mass deformation modulus. The flat-jack test is used to determine the in situ stress conditions of the host rock. For the plate bearing test, two finite element models are used to simulate the classic problem of a load on an elastic half space and the actual problem of a plate bearing test in an underground tunnel of circular cross section. For the flat-jack simulation, a single finite element model is used to simulate both horizontal and vertical slots. Results will be compared to closed-form solutions available in the literature
Finite Element Modelling Of Solidification Of Zinc Alloy | Osinkolu ...
African Journals Online (AJOL)
The solidification process of Zinc alloy is modelled by solving heat transfer equations with the aid of finite element method (FEM) using appropriate boundary conditions at the mould walls. The commercial software, Matlab, has been used to model the solidification process. The temperature profiles for each casting condition ...
Deflation in preconditioned conjugate gradient methods for Finite Element Problems
Vermolen, F.J.; Vuik, C.; Segal, A.
2002-01-01
We investigate the influence of the value of deflation vectors at interfaces on the rate of convergence of preconditioned conjugate gradient methods applied to a Finite Element discretization for an elliptic equation. Our set-up is a Poisson problem in two dimensions with continuous or discontinuous
Bending analysis of laminated composite plates using finite element ...
African Journals Online (AJOL)
user
In the past, the structural behavior of plates and shells using the finite element method has been studied by a variety of approaches. Choudhary and Tungikar ... (2011) presented the nonlinear static analysis of a rectangular laminated composite thick plate resting on nonlinear two-parameter elastic foundation with cubic.
Can finite element models detect clinically inferior cemented hip implants?
Stolk, J.; Maher, S.A.; Verdonschot, N.J.J.; Prendergast, P.J.; Huiskes, R.
2003-01-01
Rigorous preclinical testing of cemented hip prostheses against the damage accumulation failure scenario will reduce the incidence of aseptic loosening. For that purpose, a finite element simulation is proposed that predicts damage accumulation in the cement mantle and prosthetic migration. If the
Finite element analysis of one–dimensional hydrodynamic ...
African Journals Online (AJOL)
In this research work, we consider the one dimensional hydrodynamic dispersion of a reactive solute in electroosmotic flow. We present results demonstrating the utility of finite element methods to simulate and visualize hydrodynamic dispersion in the electroosmotic flow. From examination of concentration profile, effective ...
Finite element concept to derive isostatic residual maps-Application ...
Indian Academy of Sciences (India)
A new space-domain operator based on the shape function concept of finite element analysis has been developed to derive the residual maps of the Gorda Plate of western United States. The technique does not require explicit assumptions on isostatic models. Besides delineating the Gorda Plate boundary, the residual ...
Finite element analysis of bone loss around failing implants
Wolff, J.E.H.; Narra, N.; Antalainen, A.K.; Valasek, J.; Kaiser, J.; Sandor, G.K.; Marcian, P.
2014-01-01
Dental implants induce diverse forces on their surrounding bone. However, when excessive unphysiological forces are applied, resorption of the neighbouring bone may occur. The aim of this study was to assess possible causes of bone loss around failing dental implants using finite element analysis. A
Appendix F : finite element analysis of end region.
2013-03-01
FE (finite element) modeling was conducted to 1) provide a better understanding of the : elastic behavior of the end region prior to cracking and 2) to evaluate the effects of bearing pad : stiffness and width on end region elastic stresses. The FEA ...
Finite element investigation of the prestressed jointed concrete ...
African Journals Online (AJOL)
Precast prestressed concrete pavement (PCP) technology is of recent origin, and the information on PCP performance is not available in literature. This research presents a finite-element analysis of the potential benefits of prestressing on the jointed concrete pavements (JCP). With using a 3-dimensional (3D) ...
Material Models for the Human Torso Finite Element Model
2018-04-04
ARL-TR-8338 ● Apr 2018 US Army Research Laboratory Material Models for the Human Torso Finite Element Model by Carolyn E...longer needed. Do not return it to the originator. ARL-TR-8338 ● Apr 2018 US Army Research Laboratory Material Models for the...Weapons and Materials Research Directorate, ARL Approved for public release; distribution is unlimited. ii REPORT
finite element model for predicting residual stresses in shielded
African Journals Online (AJOL)
eobe
Diffractometer (XRD 6000). From the Finite Element Model Simulation, the transverse residual stress in the x ... Keywords: Residual stress, 3D FEM, Shielded manual metal arc welding, Low Carbon Steel (ASTM A36), X-Ray diffraction, degree of ..... I. ''Residual stress effects on fatigue life of welded structures using LEFM'',.
Finite element analyses of wood laminated composite poles
Cheng Piao; Todd F. Shupe; R.C. Tang; Chung Y. Hse
2005-01-01
Finite element analyses using ANSYS were conducted on orthotropic, polygonal, wood laminated composite poles subjected to a body force and a concentrated load at the free end. Deflections and stress distributions of small-scale and full-size composite poles were analyzed and compared to the results obtained in an experimental study. The predicted deflection for both...
GRIZ: Visualization of finite element analysis results on unstructured grids
International Nuclear Information System (INIS)
Dovey, D.; Loomis, M.D.
1994-01-01
GRIZ is a general-purpose post-processing application that supports interactive visualization of finite element analysis results on three-dimensional unstructured grids. GRIZ includes direct-to-videodisc animation capabilities and is being used as a production tool for creating engineering animations
CSIR Research Space (South Africa)
Suliman, Ridhwaan
2015-01-01
Full Text Available A fully-coupled partitioned finite volume–finite volume and hybrid finite volume–finite element fluid-structure interaction scheme is presented. The fluid domain is modelled as a viscous incompressible isothermal region governed by the Navier...
Integral finite element analysis of turntable bearing with flexible rings
Deng, Biao; Liu, Yunfei; Guo, Yuan; Tang, Shengjin; Su, Wenbin; Lei, Zhufeng; Wang, Pengcheng
2018-03-01
This paper suggests a method to calculate the internal load distribution and contact stress of the thrust angular contact ball turntable bearing by FEA. The influence of the stiffness of the bearing structure and the plastic deformation of contact area on the internal load distribution and contact stress of the bearing is considered. In this method, the load-deformation relationship of the rolling elements is determined by the finite element contact analysis of a single rolling element and the raceway. Based on this, the nonlinear contact between the rolling elements and the inner and outer ring raceways is same as a nonlinear compression spring and bearing integral finite element analysis model including support structure was established. The effects of structural deformation and plastic deformation on the built-in stress distribution of slewing bearing are investigated on basis of comparing the consequences of load distribution, inner and outer ring stress, contact stress and other finite element analysis results with the traditional bearing theory, which has guiding function for improving the design of slewing bearing.
Dedicated finite elements for electrode thin films on quartz resonators.
Srivastava, Sonal A; Yong, Yook-Kong; Tanaka, Masako; Imai, Tsutomu
2008-08-01
The accuracy of the finite element analysis for thickness shear quartz resonators is a function of the mesh resolution; the finer the mesh resolution, the more accurate the finite element solution. A certain minimum number of elements are required in each direction for the solution to converge. This places a high demand on memory for computation, and often the available memory is insufficient. Typically the thickness of the electrode films is very small compared with the thickness of the resonator itself; as a result, electrode elements have very poor aspect ratios, and this is detrimental to the accuracy of the result. In this paper, we propose special methods to model the electrodes at the crystal interface of an AT cut crystal. This reduces the overall problem size and eliminates electrode elements having poor aspect ratios. First, experimental data are presented to demonstrate the effects of electrode film boundary conditions on the frequency-temperature curves of an AT cut plate. Finite element analysis is performed on a mesh representing the resonator, and the results are compared for testing the accuracy of the analysis itself and thus validating the results of analysis. Approximations such as lumping and Guyan reduction are then used to model the electrode thin films at the electrode interface and their results are studied. In addition, a new approximation called merging is proposed to model electrodes at the electrode interface.
Finite element analysis of degraded concrete structures - Workshop proceedings
International Nuclear Information System (INIS)
1999-09-01
This workshop is related to the finite element analysis of degraded concrete structures. It is composed of three sessions. The first session (which title is: the use of finite element analysis in safety assessments) comprises six papers which titles are: Historical Development of Concrete Finite Element Modeling for Safety Evaluation of Accident-Challenged and Aging Concrete Structures; Experience with Finite Element Methods for Safety Assessments in Switzerland; Stress State Analysis of the Ignalina NPP Confinement System; Prestressed Containment: Behaviour when Concrete Cracking is Modelled; Application of FEA for Design and Support of NPP Containment in Russia; Verification Problems of Nuclear Installations Safety Software of Strength Analysis (NISS SA). The second session (title: concrete containment structures under accident loads) comprises seven papers which titles are: Two Application Examples of Concrete Containment Structures under Accident Load Conditions Using Finite Element Analysis; What Kind of Prediction for Leak rates for Nuclear Power Plant Containments in Accidental Conditions; Influence of Different Hypotheses Used in Numerical Models for Concrete At Elevated Temperatures on the Predicted Behaviour of NPP Core Catchers Under Severe Accident Conditions; Observations on the Constitutive Modeling of Concrete Under Multi-Axial States at Elevated Temperatures; Analyses of a Reinforced Concrete Containment with Liner Corrosion Damage; Program of Containment Concrete Control During Operation for the Temelin Nuclear Power Plant; Static Limit Load of a Deteriorated Hyperbolic Cooling Tower. The third session (concrete structures under extreme environmental load) comprised five papers which titles are: Shear Transfer Mechanism of RC Plates After Cracking; Seismic Back Calculation of an Auxiliary Building of the Nuclear Power Plant Muehleberg, Switzerland; Seismic Behaviour of Slightly Reinforced Shear Wall Structures; FE Analysis of Degraded Concrete
Plane-wave basis finite elements and boundary elements for three-dimensional wave scattering.
Perrey-Debain, E; Laghrouche, O; Bettess, P; Trevelyan, J
2004-03-15
Classical finite-element and boundary-element formulations for the Helmholtz equation are presented, and their limitations with respect to the number of variables needed to model a wavelength are explained. A new type of approximation for the potential is described in which the usual finite-element and boundary-element shape functions are modified by the inclusion of a set of plane waves, propagating in a range of directions evenly distributed on the unit sphere. Compared with standard piecewise polynomial approximation, the plane-wave basis is shown to give considerable reduction in computational complexity. In practical terms, it is concluded that the frequency for which accurate results can be obtained, using these new techniques, can be up to 60 times higher than that of the conventional finite-element method, and 10 to 15 times higher than that of the conventional boundary-element method.
International Nuclear Information System (INIS)
Correia Filho, A.
1981-04-01
The Neutron Diffusion Equation at two groups of energy is solved with the use of the Finite - Element Method with first order triangular elements. The program EFTDN (Triangular Finite Elements on Neutron Diffusion) was developed using the language FORTRAN IV. The discrete formulation of the Diffusion Equation is obtained with the application of the Galerkin's Method. In order to solve the eigenvalue - problem, the Method of the Power is applied and, with the purpose of the convergence of the results, Chebshev's polynomial expressions are applied. On the solution of the systems of equations Gauss' Method is applied, divided in two different parts: triangularization of the matrix of coeficients and retrosubstitution taking in account the sparsity of the system. Several test - problems are solved, among then two P.W.R. type reactors, the ZION-1 with 1300 MWe and the 2D-IAEA - Benchmark. Comparision of results with standard solutions show the validity of application of the EFM and precision of the results. (Author) [pt
Investigations on Actuator Dynamics through Theoretical and Finite Element Approach
Directory of Open Access Journals (Sweden)
Somashekhar S. Hiremath
2010-01-01
Full Text Available This paper gives a new approach for modeling the fluid-structure interaction of servovalve component-actuator. The analyzed valve is a precision flow control valve-jet pipe electrohydraulic servovalve. The positioning of an actuator depends upon the flow rate from control ports, in turn depends on the spool position. Theoretical investigation is made for No-load condition and Load condition for an actuator. These are used in finite element modeling of an actuator. The fluid-structure-interaction (FSI is established between the piston and the fluid cavities at the piston end. The fluid cavities were modeled with special purpose hydrostatic fluid elements while the piston is modeled with brick elements. The finite element method is used to simulate the variation of cavity pressure, cavity volume, mass flow rate, and the actuator velocity. The finite element analysis is extended to study the system's linearized response to harmonic excitation using direct solution steady-state dynamics. It was observed from the analysis that the natural frequency of the actuator depends upon the position of the piston in the cylinder. This is a close match with theoretical and simulation results. The effect of bulk modulus is also presented in the paper.
FEM and BIEM - A new infinite hybrid finite element
International Nuclear Information System (INIS)
Drakaliev, P.
1993-01-01
The finite element method (MFE) and the boundary integral equation method (BIEM) are general approximation procedures applicable to a wide variety of engineering problems. Each of them has many variants and each possesses certain merits and limitations of its own. The FEM may be easier to apply in domains with anisotropic or nonlinear behaviour. On the other hand the BEM is more attractive for unbounded domains or regions of high stress concentration. Therefore, the idea of combining both numerical techniques is of great interest in many practical problems, especially in solid and fluid mechanics, such as soil-structure and structure-fluid interaction problems. In the developments to follow an energy approach for symmetrizing the indirect BIEM is being used to obtain the stiffness matrix for the infinite or semi-infinite elastic medium. Thus the subdomain is considered as an infinite super element with an arbitrary shaped boundary and can be easily implemented into existing finite element codes
ADAPTIVE FINITE ELEMENT MODELING TECHNIQUES FOR THE POISSON-BOLTZMANN EQUATION
HOLST, MICHAEL; MCCAMMON, JAMES ANDREW; YU, ZEYUN; ZHOU, YOUNGCHENG; ZHU, YUNRONG
2011-01-01
We consider the design of an effective and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the two-term regularization technique for the continuous problem recently proposed by Chen, Holst, and Xu based on the removal of the singular electrostatic potential inside biomolecules; this technique made possible the development of the first complete solution and approximation theory for the Poisson-Boltzmann equation, the first provably convergent discretization, and also allowed for the development of a provably convergent AFEM. However, in practical implementation, this two-term regularization exhibits numerical instability. Therefore, we examine a variation of this regularization technique which can be shown to be less susceptible to such instability. We establish a priori estimates and other basic results for the continuous regularized problem, as well as for Galerkin finite element approximations. We show that the new approach produces regularized continuous and discrete problems with the same mathematical advantages of the original regularization. We then design an AFEM scheme for the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori L∞ estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. All of the algorithms described in the article are implemented in the Finite Element Toolkit (FETK), developed and maintained at UCSD. The stability advantages of the new regularization scheme
Energy Technology Data Exchange (ETDEWEB)
Qi, M.; Wegner, J.; Ganzer, L. [Technische Univ. Clausthal, Clausthal-Zellerfeld (Germany). ITE
2013-08-01
Polymer flooding, as an EOR method, has become one of the most important driving forces after water flooding. The conventional believe is that polymer flooding can only improve sweep efficiency, but it has no contribution to residual oil saturation reduction. However, experimental studies indicated that polymer solution can also improve displacement efficiency and decrease residual oil saturation. To get a better understanding of the mechanism to increase the microscopic sweep efficiency and the displacement efficiency, theoretical studies are required. In this paper, we studied the viscoelasticity effect of polymer by using a numerical simulator, which is based on Finite Element Analysis. Since it is showed experimentally that the first normal stress difference of viscoelastic polymer solution is higher than the second stress difference, the Oldroyd-B model was selected as the constitutive equation in the simulation. Numerical modelling of Oldroyd-B viscoelastic fluids is notoriously difficult. Standard Galerkin finite element methods are prone to numerical oscillations, and there is no convergence as the elasticity of fluid increases. Therefore, we use a stabilised finite element model. In order to verify our model, we first built up a model with the same geometry and fluid properties as presented in literature and compared the results. Then, with the tested model we simulated the effect of viscoelastic polymer fluid on dead pores in three simplified pore structures, which are contraction structure, expansion structure and expansion-contraction structure. Correspondingly, the streamlines and velocity contours of polymer solution, with different Reynolds numbers (Re) and Weissenberg numbers (We), flowing in these three structures are showed. The simulation results indicate that the viscoelasticity of polymer solution is the main contribution to increase the micro-scale sweep efficiency. With higher elasticity, the velocity of polymer solution is getting bigger at
Solution of Fokker–Planck equation by finite element and finite ...
Indian Academy of Sciences (India)
Abstract. The response of a structural system to white noise excitation (delta- correlated) constitutes a Markov vector process whose transitional probability den- sity function (TPDF) is governed by both the forward Fokker–Planck and backward. Kolmogorov equations. Numerical solution of these equations by finite element ...
Solution of Fokker–Planck equation by finite element and finite ...
Indian Academy of Sciences (India)
hindered by the problem of dimensionality. In this paper numerical solution of the stationary and transient form of the Fokker–Planck (FP) equation corresponding to two state nonlinear systems is obtained by standard sequential finite element method. (FEM) using C0 shape function and Crank–Nicholson time integration ...
Finite element simulations of surface effect on Rayleigh waves
He, Jin; Zhao, Jinling
2018-03-01
Rayleigh waves influenced by surface effect are investigated by using finite element methods, in which eigenfrequency analysis are performed on a model composed of a half-space covered by the surface effect dominated domain. For a given wavelength, the frequency of the Rayleigh wave is obtained as the eigenfrequency of the model satisfying Floquet periodic boundary conditions. The thickness of the surface effect can be set to be infinitely small or a finite value in the finite element methods. The curvature-dependent out-of-plane force induced by surface tension as described by the generalized Young-Laplace equation is realized through geometric nonlinear analysis. The finite element simulations show that the assumptions of small curvature and infinitely small thickness of the surface effect widely used in theoretical approaches become invalid when Rayleigh waves are highly influenced by the surface effect. This work gives a more accurate insight into the surface effect on Rayleigh waves and provides a potential method for measuring the thickness of the surface effect from the dispersion curves of surface effect influenced Rayleigh wave velocities.
Finite-element-analysis of fields radiated from ICRF antenna
International Nuclear Information System (INIS)
Yamanaka, Kaoru; Sugihara, Ryo.
1984-01-01
The electromagnetic fields radiated from a loop antenna on which an oscillating current flows across the static magnetic field B 0 are calculated in several simple geometries by the finite element method (FEM) and by analytical methods in a cross section of a plasma cylinder. The wave number along B 0 is assumed to be finite. Good agreement between FEM and the analytical solutions is obtained, demonstrating the accuracy of the FEM solutions. The method is used to calculate the fields from a half-turn antenna, and acceptable results are obtained. (author)
Solution of Exterior Helmholtz Problems Using Finite and Infinite Elements
Shirron, Joseph James
This dissertation discusses methods for the computation of solutions of the Helmholtz equation in unbounded domains. Two classes of methods are considered: one in which the infinite exterior domain is truncated and finite elements are used to discretize the resultant computational domain, and another in which the exterior domain is discretized by infinite elements. For the first class of methods a generalized Robin boundary condition is imposed on the truncating surface to replace the Sommerfeld radiation condition at infinity and to ensure uniqueness of the solution. Several of these approximate radiation conditions are discussed and a comparison is presented to illustrate their efficacy. For the second class of methods finite elements are used to discretize the exterior domain out to an enclosing circle or prolate spheroid, then infinite elements are used to discretize the remaining unbounded domain. Strikingly different approximation and convergence behavior is observed depending on whether a bilinear or sesquilinear form is chosen for the variational formulation of the problem. Convergence analysis for the infinite element methods is presented for both two and three spatial dimensions. A solution method based on the idea of domain decomposition is also discussed, as are various techniques for obtaining the solution in the far field. Numerical experiments for problems of acoustic scattering by bodies of revolution convincingly demonstrate the superiority in terms of computational expense of the infinite element methods over boundary element methods.
Finite Element Framework for Computational Fluid Dynamics in FEBio.
Ateshian, Gerard A; Shim, Jay J; Maas, Steve A; Weiss, Jeffrey A
2018-02-01
The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid-structure interactions (FSIs).
Finite element analysis of structures through unified formulation
Carrera, Erasmo; Petrolo, Marco; Zappino, Enrico
2014-01-01
The finite element method (FEM) is a computational tool widely used to design and analyse complex structures. Currently, there are a number of different approaches to analysis using the FEM that vary according to the type of structure being analysed: beams and plates may use 1D or 2D approaches, shells and solids 2D or 3D approaches, and methods that work for one structure are typically not optimized to work for another. Finite Element Analysis of Structures Through Unified Formulation deals with the FEM used for the analysis of the mechanics of structures in the case of linear elasticity. The novelty of this book is that the finite elements (FEs) are formulated on the basis of a class of theories of structures known as the Carrera Unified Formulation (CUF). It formulates 1D, 2D and 3D FEs on the basis of the same ''fundamental nucleus'' that comes from geometrical relations and Hooke''s law, and presents both 1D and 2D refined FEs that only have displacement variables as in 3D elements. It also covers 1D...
Choice of input fields in stochastic finite elements
DEFF Research Database (Denmark)
Ditlevsen, Ove Dalager; Tarp-Johansen, Niels Jacob
1999-01-01
, the flexibility field, as the input to the stochastic finite element model. To answer this question the focus should be on the error of the output of the mechanical model rather than on the input field itself when discretizing the held through replacing it by a field defined in terms of a finite number of random...... variables. Several reported discretization methods define these random variables as integrals of the product of the held and some suitable weight functions. In particular, the weight functions can be Dirac delta functions whereby the random variables become the field values at a finite set of given points...... the differential equation of the column displacement and the relevant boundary conditions, it can be expected that the discretization of the flexibility field is preferable over the discretization of the stiffness field. Direct mechanical considerations support this expectation. (C) 1998 Published by Elsevier...
Choice of input fields in stochastic finite elements
DEFF Research Database (Denmark)
Ditlevsen, Ove Dalager; Tarp-Johansen, Niels Jacob
1996-01-01
, the flexibility field, as the input to the stochastic finite element model. To answer this question the focus should be on the error of the output of the mechanical model rather than on the input field itself when discretizing the field through replacing it by a field defined in terms of a finite number of random...... variables. Several reported discretization methods define these random variables as integrals of the product of the field and some suitable weight functions. In particular, the weight functions can be Dirac delta functions whereby the random variables become the field values at a finite set of given points...... the differential equation of the column displacement and the relevant boundarv conditions, it can be expected that the discretization of the flexibility field is preferable over the discretization of the stiffness field. Direct mechanical considerations support this expectation.Keywords: Random stiffness...