A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids

Energy Technology Data Exchange (ETDEWEB)

Bailey, T.S.; Adams, M.L. [Texas A M Univ., Dept. of Nuclear Engineering, College Station, TX (United States); Yang, B.; Zika, M.R. [Lawrence Livermore National Lab., Livermore, CA (United States)

2005-07-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)

Fracture Failure of Reinforced Concrete Slabs Subjected to Blast Loading Using the Combined Finite-Discrete Element Method

Directory of Open Access Journals (Sweden)

Z. M. Jaini

Full Text Available Abstract Numerical modeling of fracture failure is challenging due to various issues in the constitutive law and the transition of continuum to discrete bodies. Therefore, this study presents the application of the combined finite-discrete element method to investigate the fracture failure of reinforced concrete slabs subjected to blast loading. In numerical modeling, the interaction of non-uniform blast loading on the concrete slab was modeled using the incorporation of the finite element method with a crack rotating approach and the discrete element method to model crack, fracture onset and its post-failures. A time varying pressure-time history based on the mapping method was adopted to define blast loading. The Mohr-Coulomb with Rankine cut-off and von-Mises criteria were applied for concrete and steel reinforcement respectively. The results of scabbing, spalling and fracture show a reliable prediction of damage and fracture.

Finite elements and approximation

CERN Document Server

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

Finite Volume Element (FVE) discretization and multilevel solution of the axisymmetric heat equation

Science.gov (United States)

Litaker, Eric T.

1994-12-01

The axisymmetric heat equation, resulting from a point-source of heat applied to a metal block, is solved numerically; both iterative and multilevel solutions are computed in order to compare the two processes. The continuum problem is discretized in two stages: finite differences are used to discretize the time derivatives, resulting is a fully implicit backward time-stepping scheme, and the Finite Volume Element (FVE) method is used to discretize the spatial derivatives. The application of the FVE method to a problem in cylindrical coordinates is new, and results in stencils which are analyzed extensively. Several iteration schemes are considered, including both Jacobi and Gauss-Seidel; a thorough analysis of these schemes is done, using both the spectral radii of the iteration matrices and local mode analysis. Using this discretization, a Gauss-Seidel relaxation scheme is used to solve the heat equation iteratively. A multilevel solution process is then constructed, including the development of intergrid transfer and coarse grid operators. Local mode analysis is performed on the components of the amplification matrix, resulting in the two-level convergence factors for various combinations of the operators. A multilevel solution process is implemented by using multigrid V-cycles; the iterative and multilevel results are compared and discussed in detail. The computational savings resulting from the multilevel process are then discussed.

Analysis of Tire Tractive Performance on Deformable Terrain by Finite Element-Discrete Element Method

Science.gov (United States)

Nakashima, Hiroshi; Takatsu, Yuzuru

The goal of this study is to develop a practical and fast simulation tool for soil-tire interaction analysis, where finite element method (FEM) and discrete element method (DEM) are coupled together, and which can be realized on a desktop PC. We have extended our formerly proposed dynamic FE-DE method (FE-DEM) to include practical soil-tire system interaction, where not only the vertical sinkage of a tire, but also the travel of a driven tire was considered. Numerical simulation by FE-DEM is stable, and the relationships between variables, such as load-sinkage and sinkage-travel distance, and the gross tractive effort and running resistance characteristics, are obtained. Moreover, the simulation result is accurate enough to predict the maximum drawbar pull for a given tire, once the appropriate parameter values are provided. Therefore, the developed FE-DEM program can be applied with sufficient accuracy to interaction problems in soil-tire systems.

Finite Discrete Gabor Analysis

DEFF Research Database (Denmark)

Søndergaard, Peter Lempel

2007-01-01

frequency bands at certain times. Gabor theory can be formulated for both functions on the real line and for discrete signals of finite length. The two theories are largely the same because many aspects come from the same underlying theory of locally compact Abelian groups. The two types of Gabor systems...... can also be related by sampling and periodization. This thesis extends on this theory by showing new results for window construction. It also provides a discussion of the problems associated to discrete Gabor bases. The sampling and periodization connection is handy because it allows Gabor systems...... on the real line to be well approximated by finite and discrete Gabor frames. This method of approximation is especially attractive because efficient numerical methods exists for doing computations with finite, discrete Gabor systems. This thesis presents new algorithms for the efficient computation of finite...

Evaluation and optimization of footwear comfort parameters using finite element analysis and a discrete optimization algorithm

Science.gov (United States)

Papagiannis, P.; Azariadis, P.; Papanikos, P.

2017-10-01

Footwear is subject to bending and torsion deformations that affect comfort perception. Following review of Finite Element Analysis studies of sole rigidity and comfort, a three-dimensional, linear multi-material finite element sole model for quasi-static bending and torsion simulation, overcoming boundary and optimisation limitations, is described. Common footwear materials properties and boundary conditions from gait biomechanics are used. The use of normalised strain energy for product benchmarking is demonstrated along with comfort level determination through strain energy density stratification. Sensitivity of strain energy against material thickness is greater for bending than for torsion, with results of both deformations showing positive correlation. Optimization for a targeted performance level and given layer thickness is demonstrated with bending simulations sufficing for overall comfort assessment. An algorithm for comfort optimization w.r.t. bending is presented, based on a discrete approach with thickness values set in line with practical manufacturing accuracy. This work illustrates the potential of the developed finite element analysis applications to offer viable and proven aids to modern footwear sole design assessment and optimization.

Finite element bending behaviour of discretely delaminated ...

African Journals Online (AJOL)

user

due to their light weight, high specific strength and stiffness properties. ... cylindrical shell roofs respectively using finite element method with centrally located .... where { }ε and { }γ are the direct and shear strains in midplane and { }κ denotes ...

On the Robustness and Prospects of Adaptive BDDC Methods for Finite Element Discretizations of Elliptic PDEs with High-Contrast Coefficients

KAUST Repository

Zampini, Stefano; Keyes, David E.

2016-01-01

Balancing Domain Decomposition by Constraints (BDDC) methods have proven to be powerful preconditioners for large and sparse linear systems arising from the finite element discretization of elliptic PDEs. Condition number bounds can be theoretically

High-Order Curvilinear Finite Element Methods for Lagrangian Hydrodynamics [High Order Curvilinear Finite Elements for Lagrangian Hydrodynamics

Energy Technology Data Exchange (ETDEWEB)

Dobrev, Veselin A. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Kolev, Tzanio V. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Rieben, Robert N. [Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)

2012-09-20

The numerical approximation of the Euler equations of gas dynamics in a movingLagrangian frame is at the heart of many multiphysics simulation algorithms. Here, we present a general framework for high-order Lagrangian discretization of these compressible shock hydrodynamics equations using curvilinear finite elements. This method is an extension of the approach outlined in [Dobrev et al., Internat. J. Numer. Methods Fluids, 65 (2010), pp. 1295--1310] and can be formulated for any finite dimensional approximation of the kinematic and thermodynamic fields, including generic finite elements on two- and three-dimensional meshes with triangular, quadrilateral, tetrahedral, or hexahedral zones. We discretize the kinematic variables of position and velocity using a continuous high-order basis function expansion of arbitrary polynomial degree which is obtained via a corresponding high-order parametric mapping from a standard reference element. This enables the use of curvilinear zone geometry, higher-order approximations for fields within a zone, and a pointwise definition of mass conservation which we refer to as strong mass conservation. Moreover, we discretize the internal energy using a piecewise discontinuous high-order basis function expansion which is also of arbitrary polynomial degree. This facilitates multimaterial hydrodynamics by treating material properties, such as equations of state and constitutive models, as piecewise discontinuous functions which vary within a zone. To satisfy the Rankine--Hugoniot jump conditions at a shock boundary and generate the appropriate entropy, we introduce a general tensor artificial viscosity which takes advantage of the high-order kinematic and thermodynamic information available in each zone. Finally, we apply a generic high-order time discretization process to the semidiscrete equations to develop the fully discrete numerical algorithm. Our method can be viewed as the high-order generalization of the so-called staggered

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

Wheeler, Mary Fanett; Xue, Guangri; Yotov, Ivan

2012-01-01

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

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.

SEACAS Theory Manuals: Part III. Finite Element Analysis in Nonlinear Solid Mechanics

Energy Technology Data Exchange (ETDEWEB)

Laursen, T.A.; Attaway, S.W.; Zadoks, R.I.

1999-03-01

This report outlines the application of finite element methodology to large deformation solid mechanics problems, detailing also some of the key technological issues that effective finite element formulations must address. The presentation is organized into three major portions: first, a discussion of finite element discretization from the global point of view, emphasizing the relationship between a virtual work principle and the associated fully discrete system, second, a discussion of finite element technology, emphasizing the important theoretical and practical features associated with an individual finite element; and third, detailed description of specific elements that enjoy widespread use, providing some examples of the theoretical ideas already described. Descriptions of problem formulation in nonlinear solid mechanics, nonlinear continuum mechanics, and constitutive modeling are given in three companion reports.

Finite elements of nonlinear continua

CERN Document Server

Oden, John Tinsley

1972-01-01

Geared toward undergraduate and graduate students, this text extends applications of the finite element method from linear problems in elastic structures to a broad class of practical, nonlinear problems in continuum mechanics. It treats both theory and applications from a general and unifying point of view.The text reviews the thermomechanical principles of continuous media and the properties of the finite element method, and then brings them together to produce discrete physical models of nonlinear continua. The mathematical properties of these models are analyzed, along with the numerical s

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

Variational approach to probabilistic finite elements

Science.gov (United States)

Belytschko, T.; Liu, W. K.; Mani, A.; Besterfield, G.

1991-08-01

Probabilistic finite element methods (PFEM), synthesizing the power of finite element methods with second-moment techniques, are formulated for various classes of problems in structural and solid mechanics. Time-invariant random materials, geometric properties and loads are incorporated in terms of their fundamental statistics viz. second-moments. Analogous to the discretization of the displacement field in finite element methods, the random fields are also discretized. Preserving the conceptual simplicity, the response moments are calculated with minimal computations. By incorporating certain computational techniques, these methods are shown to be capable of handling large systems with many sources of uncertainties. By construction, these methods are applicable when the scale of randomness is not very large and when the probabilistic density functions have decaying tails. The accuracy and efficiency of these methods, along with their limitations, are demonstrated by various applications. Results obtained are compared with those of Monte Carlo simulation and it is shown that good accuracy can be obtained for both linear and nonlinear problems. The methods are amenable to implementation in deterministic FEM based computer codes.

About solution of multipoint boundary problem of static analysis of deep beam with the use of combined application of finite element method and discrete-continual finite element method. part 1: formulation of the problem and general principles of approximation

Directory of Open Access Journals (Sweden)

Lyakhovich Leonid

2017-01-01

Full Text Available This paper is devoted to formulation and general principles of approximation of multipoint boundary problem of static analysis of deep beam with the use of combined application of finite element method (FEM discrete-continual finite element method (DCFEM. The field of application of DCFEM comprises structures with regular physical and geometrical parameters in some dimension (“basic” dimension. DCFEM presupposes finite element approximation for non-basic dimension while in the basic dimension problem remains continual. DCFEM is based on analytical solutions of resulting multipoint boundary problems for systems of ordinary differential equations with piecewise-constant coefficients.

Domain decomposition based iterative methods for nonlinear elliptic finite element problems

Energy Technology Data Exchange (ETDEWEB)

Cai, X.C. [Univ. of Colorado, Boulder, CO (United States)

1994-12-31

The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.

High order curvilinear finite elements for elastic–plastic Lagrangian dynamics

International Nuclear Information System (INIS)

Dobrev, Veselin A.; Kolev, Tzanio V.; Rieben, Robert N.

2014-01-01

This paper presents a high-order finite element method for calculating elastic–plastic flow on moving curvilinear meshes and is an extension of our general high-order curvilinear finite element approach for solving the Euler equations of gas dynamics in a Lagrangian frame [1,2]. In order to handle transition to plastic flow, we formulate the stress–strain relation in rate (or incremental) form and augment our semi-discrete equations for Lagrangian hydrodynamics with an additional evolution equation for the deviatoric stress which is valid for arbitrary order spatial discretizations of the kinematic and thermodynamic variables. The semi-discrete equation for the deviatoric stress rate is developed for 2D planar, 2D axisymmetric and full 3D geometries. For each case, the strain rate is approximated via a collocation method at zone quadrature points while the deviatoric stress is approximated using an L 2 projection onto the thermodynamic basis. We apply high order, energy conserving, explicit time stepping methods to the semi-discrete equations to develop the fully discrete method. We conclude with numerical results from an extensive series of verification tests that demonstrate several practical advantages of using high-order finite elements for elastic–plastic flow

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

Finite element formulation for a digital image correlation method

International Nuclear Information System (INIS)

Sun Yaofeng; Pang, John H. L.; Wong, Chee Khuen; Su Fei

2005-01-01

A finite element formulation for a digital image correlation method is presented that will determine directly the complete, two-dimensional displacement field during the image correlation process on digital images. The entire interested image area is discretized into finite elements that are involved in the common image correlation process by use of our algorithms. This image correlation method with finite element formulation has an advantage over subset-based image correlation methods because it satisfies the requirements of displacement continuity and derivative continuity among elements on images. Numerical studies and a real experiment are used to verify the proposed formulation. Results have shown that the image correlation with the finite element formulation is computationally efficient, accurate, and robust

Analysis of a block Gauss-Seidel iterative method for a finite element discretization of the neutron transport equation

International Nuclear Information System (INIS)

Lorence, L.J. Jr.; Martin, W.R.; Luskin, M.

1985-01-01

We prove the convergence of a finite element discretization of the neutron transport equation. The iterative solution of the resulting linear system by a block Gauss-Seidel method is also analyzed. This procedure is shown to require less storage than the direct solution by Gaussian elimination, and an estimate for the rate of convergence is used to show that fewer arithmetic operations are required

Applications of discrete element method in modeling of grain postharvest operations

Science.gov (United States)

Grain kernels are finite and discrete materials. Although flowing grain can behave like a continuum fluid at times, the discontinuous behavior exhibited by grain kernels cannot be simulated solely with conventional continuum-based computer modeling such as finite-element or finite-difference methods...

Probabilistic finite elements for transient analysis in nonlinear continua

Science.gov (United States)

Liu, W. K.; Belytschko, T.; Mani, A.

1985-01-01

The probabilistic finite element method (PFEM), which is a combination of finite element methods and second-moment analysis, is formulated for linear and nonlinear continua with inhomogeneous random fields. Analogous to the discretization of the displacement field in finite element methods, the random field is also discretized. The formulation is simplified by transforming the correlated variables to a set of uncorrelated variables through an eigenvalue orthogonalization. Furthermore, it is shown that a reduced set of the uncorrelated variables is sufficient for the second-moment analysis. Based on the linear formulation of the PFEM, the method is then extended to transient analysis in nonlinear continua. The accuracy and efficiency of the method is demonstrated by application to a one-dimensional, elastic/plastic wave propagation problem. The moments calculated compare favorably with those obtained by Monte Carlo simulation. Also, the procedure is amenable to implementation in deterministic FEM based computer programs.

Semianalytical analysis of shear walls with the use of discrete-continual finite element method. Part 1: Mathematical foundations

Directory of Open Access Journals (Sweden)

Akimov Pavel

2016-01-01

Full Text Available The distinctive paper is devoted to the two-dimensional semi-analytical solution of boundary problems of analysis of shear walls with the use of discrete-continual finite element method (DCFEM. This approach allows obtaining the exact analytical solution in one direction (so-called “basic” direction, also decrease the size of the problem to one-dimensional common finite element analysis. The resulting multipoint boundary problem for the first-order system of ordinary differential equations with piecewise constant coefficients is solved analytically. The proposed method is rather efficient for evaluation of the boundary effect (such as the stress field near the concentrated force. DCFEM also has a completely computer-oriented algorithm, computational stability, optimal conditionality of resultant system and it is applicable for the various loads at an arbitrary point or a region of the wall.

Shale Fracture Analysis using the Combined Finite-Discrete Element Method

Science.gov (United States)

Carey, J. W.; Lei, Z.; Rougier, E.; Knight, E. E.; Viswanathan, H.

2014-12-01

Hydraulic fracturing (hydrofrac) is a successful method used to extract oil and gas from highly carbonate rocks like shale. However, challenges exist for industry experts estimate that for a single $10 million dollar lateral wellbore fracking operation, only 10% of the hydrocarbons contained in the rock are extracted. To better understand how to improve hydrofrac recovery efficiencies and to lower its costs, LANL recently funded the Laboratory Directed Research and Development (LDRD) project: "Discovery Science of Hydraulic Fracturing: Innovative Working Fluids and Their Interactions with Rocks, Fractures, and Hydrocarbons". Under the support of this project, the LDRD modeling team is working with the experimental team to understand fracture initiation and propagation in shale rocks. LANL's hybrid hydro-mechanical (HM) tool, the Hybrid Optimization Software Suite (HOSS), is being used to simulate the complex fracture and fragment processes under a variety of different boundary conditions. HOSS is based on the combined finite-discrete element method (FDEM) and has been proven to be a superior computational tool for multi-fracturing problems. In this work, the comparison of HOSS simulation results to triaxial core flooding experiments will be presented.

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, \

Simulation of Semi-Solid Material Mechanical Behavior Using a Combined Discrete/Finite Element Method

Science.gov (United States)

Sistaninia, M.; Phillion, A. B.; Drezet, J.-M.; Rappaz, M.

2011-01-01

As a necessary step toward the quantitative prediction of hot tearing defects, a three-dimensional stress-strain simulation based on a combined finite element (FE)/discrete element method (DEM) has been developed that is capable of predicting the mechanical behavior of semisolid metallic alloys during solidification. The solidification model used for generating the initial solid-liquid structure is based on a Voronoi tessellation of randomly distributed nucleation centers and a solute diffusion model for each element of this tessellation. At a given fraction of solid, the deformation is then simulated with the solid grains being modeled using an elastoviscoplastic constitutive law, whereas the remaining liquid layers at grain boundaries are approximated by flexible connectors, each consisting of a spring element and a damper element acting in parallel. The model predictions have been validated against Al-Cu alloy experimental data from the literature. The results show that a combined FE/DEM approach is able to express the overall mechanical behavior of semisolid alloys at the macroscale based on the morphology of the grain structure. For the first time, the localization of strain in the intergranular regions is taken into account. Thus, this approach constitutes an indispensible step towards the development of a comprehensive model of hot tearing.

An element-based finite-volume method approach for naturally fractured compositional reservoir simulation

Energy Technology Data Exchange (ETDEWEB)

Marcondes, Francisco [Federal University of Ceara, Fortaleza (Brazil). Dept. of Metallurgical Engineering and Material Science], e-mail: marcondes@ufc.br; Varavei, Abdoljalil; Sepehrnoori, Kamy [The University of Texas at Austin (United States). Petroleum and Geosystems Engineering Dept.], e-mails: varavei@mail.utexas.edu, kamys@mail.utexas.edu

2010-07-01

An element-based finite-volume approach in conjunction with unstructured grids for naturally fractured compositional reservoir simulation is presented. In this approach, both the discrete fracture and the matrix mass balances are taken into account without any additional models to couple the matrix and discrete fractures. The mesh, for two dimensional domains, can be built of triangles, quadrilaterals, or a mix of these elements. However, due to the available mesh generator to handle both matrix and discrete fractures, only results using triangular elements will be presented. The discrete fractures are located along the edges of each element. To obtain the approximated matrix equation, each element is divided into three sub-elements and then the mass balance equations for each component are integrated along each interface of the sub-elements. The finite-volume conservation equations are assembled from the contribution of all the elements that share a vertex, creating a cell vertex approach. The discrete fracture equations are discretized only along the edges of each element and then summed up with the matrix equations in order to obtain a conservative equation for both matrix and discrete fractures. In order to mimic real field simulations, the capillary pressure is included in both matrix and discrete fracture media. In the implemented model, the saturation field in the matrix and discrete fractures can be different, but the potential of each phase in the matrix and discrete fracture interface needs to be the same. The results for several naturally fractured reservoirs are presented to demonstrate the applicability of the method. (author)

Discrete Maximum Principle for Higher-Order Finite Elements in 1D

Czech Academy of Sciences Publication Activity Database

Vejchodský, Tomáš; Šolín, Pavel

2007-01-01

Roč. 76, č. 260 (2007), s. 1833-1846 ISSN 0025-5718 R&D Projects: GA ČR GP201/04/P021 Institutional research plan: CEZ:AV0Z10190503; CEZ:AV0Z20760514 Keywords : discrete maximum principle * discrete Grren´s function * higher-order elements Subject RIV: BA - General Mathematics Impact factor: 1.230, year: 2007

Image segmentation with a finite element method

DEFF Research Database (Denmark)

Bourdin, Blaise

1999-01-01

regularization results, make possible to imagine a finite element resolution method.In a first time, the Mumford-Shah functional is introduced and some existing results are quoted. Then, a discrete formulation for the Mumford-Shah problem is proposed and its $\\Gamma$-convergence is proved. Finally, some...

A combined finite volume-nonconforming finite element scheme for compressible two phase flow in porous media

KAUST Repository

Saad, Bilal Mohammed; Saad, Mazen Naufal B M

2014-01-01

We propose and analyze a combined finite volume-nonconforming finite element scheme on general meshes to simulate the two compressible phase flow in porous media. The diffusion term, which can be anisotropic and heterogeneous, is discretized by piecewise linear nonconforming triangular finite elements. The other terms are discretized by means of a cell-centered finite volume scheme on a dual mesh, where the dual volumes are constructed around the sides of the original mesh. The relative permeability of each phase is decentred according the sign of the velocity at the dual interface. This technique also ensures the validity of the discrete maximum principle for the saturation under a non restrictive shape regularity of the space mesh and the positiveness of all transmissibilities. Next, a priori estimates on the pressures and a function of the saturation that denote capillary terms are established. These stabilities results lead to some compactness arguments based on the use of the Kolmogorov compactness theorem, and allow us to derive the convergence of a subsequence of the sequence of approximate solutions to a weak solution of the continuous equations, provided the mesh size tends to zero. The proof is given for the complete system when the density of the each phase depends on its own pressure. © 2014 Springer-Verlag Berlin Heidelberg.

A combined finite volume-nonconforming finite element scheme for compressible two phase flow in porous media

KAUST Repository

Saad, Bilal Mohammed

2014-06-28

We propose and analyze a combined finite volume-nonconforming finite element scheme on general meshes to simulate the two compressible phase flow in porous media. The diffusion term, which can be anisotropic and heterogeneous, is discretized by piecewise linear nonconforming triangular finite elements. The other terms are discretized by means of a cell-centered finite volume scheme on a dual mesh, where the dual volumes are constructed around the sides of the original mesh. The relative permeability of each phase is decentred according the sign of the velocity at the dual interface. This technique also ensures the validity of the discrete maximum principle for the saturation under a non restrictive shape regularity of the space mesh and the positiveness of all transmissibilities. Next, a priori estimates on the pressures and a function of the saturation that denote capillary terms are established. These stabilities results lead to some compactness arguments based on the use of the Kolmogorov compactness theorem, and allow us to derive the convergence of a subsequence of the sequence of approximate solutions to a weak solution of the continuous equations, provided the mesh size tends to zero. The proof is given for the complete system when the density of the each phase depends on its own pressure. © 2014 Springer-Verlag Berlin Heidelberg.

Finite element model for heat conduction in jointed rock masses

International Nuclear Information System (INIS)

Gartling, D.K.; Thomas, R.K.

1981-01-01

A computatonal procedure for simulating heat conduction in a fractured rock mass is proposed and illustrated in the present paper. The method makes use of a simple local model for conduction in the vicinity of a single open fracture. The distributions of fractures and fracture properties within the finite element model are based on a statistical representation of geologic field data. Fracture behavior is included in the finite element computation by locating local, discrete fractures at the element integration points

Simplified Qualitative Discrete Numerical Model to Determine Cracking Pattern in Brittle Materials by Means of Finite Element Method

Directory of Open Access Journals (Sweden)

J. Ochoa-Avendaño

2017-01-01

Full Text Available This paper presents the formulation, implementation, and validation of a simplified qualitative model to determine the crack path of solids considering static loads, infinitesimal strain, and plane stress condition. This model is based on finite element method with a special meshing technique, where nonlinear link elements are included between the faces of the linear triangular elements. The stiffness loss of some link elements represents the crack opening. Three experimental tests of bending beams are simulated, where the cracking pattern calculated with the proposed numerical model is similar to experimental result. The advantages of the proposed model compared to discrete crack approaches with interface elements can be the implementation simplicity, the numerical stability, and the very low computational cost. The simulation with greater values of the initial stiffness of the link elements does not affect the discontinuity path and the stability of the numerical solution. The exploded mesh procedure presented in this model avoids a complex nonlinear analysis and regenerative or adaptive meshes.

Linear finite element method for one-dimensional diffusion problems

Energy Technology Data Exchange (ETDEWEB)

Brandao, Michele A.; Dominguez, Dany S.; Iglesias, Susana M., E-mail: micheleabrandao@gmail.com, E-mail: dany@labbi.uesc.br, E-mail: smiglesias@uesc.br [Universidade Estadual de Santa Cruz (LCC/DCET/UESC), Ilheus, BA (Brazil). Departamento de Ciencias Exatas e Tecnologicas. Laboratorio de Computacao Cientifica

2011-07-01

We describe in this paper the fundamentals of Linear Finite Element Method (LFEM) applied to one-speed diffusion problems in slab geometry. We present the mathematical formulation to solve eigenvalue and fixed source problems. First, we discretized a calculus domain using a finite set of elements. At this point, we obtain the spatial balance equations for zero order and first order spatial moments inside each element. Then, we introduce the linear auxiliary equations to approximate neutron flux and current inside the element and architect a numerical scheme to obtain the solution. We offer numerical results for fixed source typical model problems to illustrate the method's accuracy for coarse-mesh calculations in homogeneous and heterogeneous domains. Also, we compare the accuracy and computational performance of LFEM formulation with conventional Finite Difference Method (FDM). (author)

Solving the incompressible surface Navier-Stokes equation by surface finite elements

Science.gov (United States)

Reuther, Sebastian; Voigt, Axel

2018-01-01

We consider a numerical approach for the incompressible surface Navier-Stokes equation on surfaces with arbitrary genus g (S ) . The approach is based on a reformulation of the equation in Cartesian coordinates of the embedding R3, penalization of the normal component, a Chorin projection method, and discretization in space by surface finite elements for each component. The approach thus requires only standard ingredients which most finite element implementations can offer. We compare computational results with discrete exterior calculus simulations on a torus and demonstrate the interplay of the flow field with the topology by showing realizations of the Poincaré-Hopf theorem on n-tori.

Preconditioning for Mixed Finite Element Formulations of Elliptic Problems

KAUST Repository

Wildey, Tim; Xue, Guangri

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.

Hybrid finite difference/finite element immersed boundary method.

Science.gov (United States)

E Griffith, Boyce; Luo, Xiaoyu

2017-12-01

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach uses a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach. © 2017 The Authors International  Journal  for  Numerical  Methods  in  Biomedical  Engineering Published by John Wiley & Sons Ltd.

A finite element-based algorithm for rubbing induced vibration prediction in rotors

Science.gov (United States)

Behzad, Mehdi; Alvandi, Mehdi; Mba, David; Jamali, Jalil

2013-10-01

In this paper, an algorithm is developed for more realistic investigation of rotor-to-stator rubbing vibration, based on finite element theory with unilateral contact and friction conditions. To model the rotor, cross sections are assumed to be radially rigid. A finite element discretization based on traditional beam theories which sufficiently accounts for axial and transversal flexibility of the rotor is used. A general finite element discretization model considering inertial and viscoelastic characteristics of the stator is used for modeling the stator. Therefore, for contact analysis, only the boundary of the stator is discretized. The contact problem is defined as the contact between the circular rigid cross section of the rotor and “nodes” of the stator only. Next, Gap function and contact conditions are described for the contact problem. Two finite element models of the rotor and the stator are coupled via the Lagrange multipliers method in order to obtain the constrained equation of motion. A case study of the partial rubbing is simulated using the algorithm. The synchronous and subsynchronous responses of the partial rubbing are obtained for different rotational speeds. In addition, a sensitivity analysis is carried out with respect to the initial clearance, the stator stiffness, the damping parameter, and the coefficient of friction. There is a good agreement between the result of this research and the experimental result in the literature.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Directory of Open Access Journals (Sweden)

Qi Zhao

2014-12-01

Full Text Available Hydraulic fracturing (HF technique has been extensively used for the exploitation of unconventional oil and gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formations by fluid injection, which creates an interconnected fracture network and increases the hydrocarbon production. Meanwhile, microseismic (MS monitoring is one of the most effective approaches to evaluate such stimulation process. In this paper, the combined finite-discrete element method (FDEM is adopted to numerically simulate HF and associated MS. Several post-processing tools, including frequency-magnitude distribution (b-value, fractal dimension (D-value, and seismic events clustering, are utilized to interpret numerical results. A non-parametric clustering algorithm designed specifically for FDEM is used to reduce the mesh dependency and extract more realistic seismic information. Simulation results indicated that at the local scale, the HF process tends to propagate following the rock mass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to the maximum in-situ stress.

Mathematical aspects of finite element methods for incompressible viscous flows

Science.gov (United States)

Gunzburger, M. D.

1986-01-01

Mathematical aspects of finite element methods are surveyed for incompressible viscous flows, concentrating on the steady primitive variable formulation. The discretization of a weak formulation of the Navier-Stokes equations are addressed, then the stability condition is considered, the satisfaction of which insures the stability of the approximation. Specific choices of finite element spaces for the velocity and pressure are then discussed. Finally, the connection between different weak formulations and a variety of boundary conditions is explored.

A framework for grand scale parallelization of the combined finite discrete element method in 2d

Science.gov (United States)

Lei, Z.; Rougier, E.; Knight, E. E.; Munjiza, A.

2014-09-01

Within the context of rock mechanics, the Combined Finite-Discrete Element Method (FDEM) has been applied to many complex industrial problems such as block caving, deep mining techniques (tunneling, pillar strength, etc.), rock blasting, seismic wave propagation, packing problems, dam stability, rock slope stability, rock mass strength characterization problems, etc. The reality is that most of these were accomplished in a 2D and/or single processor realm. In this work a hardware independent FDEM parallelization framework has been developed using the Virtual Parallel Machine for FDEM, (V-FDEM). With V-FDEM, a parallel FDEM software can be adapted to different parallel architecture systems ranging from just a few to thousands of cores.

Finite element modelling of fibre-reinforced brittle materials

NARCIS (Netherlands)

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

An embedded crack in a constant strain triangle utilizing extended finite element concepts

DEFF Research Database (Denmark)

Olesen, J.F.; Poulsen, P.N.

2013-01-01

This paper revisits the formulation of the CST element with an embedded discrete crack taking advantage of the direct formulations developed within the framework of the extended finite element method, XFEM. The result is a simple element for modeling cohesive fracture processes in quasi-brittle m......This paper revisits the formulation of the CST element with an embedded discrete crack taking advantage of the direct formulations developed within the framework of the extended finite element method, XFEM. The result is a simple element for modeling cohesive fracture processes in quasi......-element discontinuity of displacements. The formulation is based on a variational principle of virtual work involving only the interpolation of displacements. The good performance of the element is demonstrated through the comparison with three benchmark tests in which a single crack is propagated: The center cracked...

A multiscale mortar multipoint flux mixed finite element method

KAUST Repository

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 of three dimensional crack growth by the use of a boundary element sub model

DEFF Research Database (Denmark)

Lucht, Tore

2009-01-01

A new automated method to model non-planar three dimensional crack growth is proposed which combines the advantages of both the boundary element method and the finite element method. The proposed method links the two methods by a submodelling strategy in which the solution of a global finite...... element model containing an approximation of the crack is interpolated to a much smaller boundary element model containing a fine discretization of the real crack. The method is validated through several numerical comparisons and by comparison to crack growth measured in a test specimen for an engineering...

Discrete phase space based on finite fields

International Nuclear Information System (INIS)

Gibbons, Kathleen S.; Hoffman, Matthew J.; Wootters, William K.

2004-01-01

The original Wigner function provides a way of representing in phase space the quantum states of systems with continuous degrees of freedom. Wigner functions have also been developed for discrete quantum systems, one popular version being defined on a 2Nx2N discrete phase space for a system with N orthogonal states. Here we investigate an alternative class of discrete Wigner functions, in which the field of real numbers that labels the axes of continuous phase space is replaced by a finite field having N elements. There exists such a field if and only if N is a power of a prime; so our formulation can be applied directly only to systems for which the state-space dimension takes such a value. Though this condition may seem limiting, we note that any quantum computer based on qubits meets the condition and can thus be accommodated within our scheme. The geometry of our NxN phase space also leads naturally to a method of constructing a complete set of N+1 mutually unbiased bases for the state space

Finite Mathematics and Discrete Mathematics: Is There a Difference?

Science.gov (United States)

Johnson, Marvin L.

Discrete mathematics and finite mathematics differ in a number of ways. First, finite mathematics has a longer history and is therefore more stable in terms of course content. Finite mathematics courses emphasize certain particular mathematical tools which are useful in solving the problems of business and the social sciences. Discrete mathematics…

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.

Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner

KAUST Repository

Gao, Longfei

2018-02-22

We consider numerical simulation of the isotropic elastic wave equations arising from seismic applications with non-trivial land topography. The more flexible finite element method is applied to the shallow region of the simulation domain to account for the topography, and combined with the more efficient finite difference method that is applied to the deep region of the simulation domain. We demonstrate that these two discretization methods, albeit starting from different formulations of the elastic wave equation, can be joined together smoothly via weakly imposed interface conditions. Discrete energy analysis is employed to derive the proper interface treatment, leading to an overall discretization that is energy-conserving. Numerical examples are presented to demonstrate the efficacy of the proposed interface treatment.

Combining finite element and finite difference methods for isotropic elastic wave simulations in an energy-conserving manner

KAUST Repository

Gao, Longfei; Keyes, David E.

2018-01-01

We consider numerical simulation of the isotropic elastic wave equations arising from seismic applications with non-trivial land topography. The more flexible finite element method is applied to the shallow region of the simulation domain to account for the topography, and combined with the more efficient finite difference method that is applied to the deep region of the simulation domain. We demonstrate that these two discretization methods, albeit starting from different formulations of the elastic wave equation, can be joined together smoothly via weakly imposed interface conditions. Discrete energy analysis is employed to derive the proper interface treatment, leading to an overall discretization that is energy-conserving. Numerical examples are presented to demonstrate the efficacy of the proposed interface treatment.

Finite element method for solving neutron transport problems

International Nuclear Information System (INIS)

Ferguson, J.M.; Greenbaum, A.

1984-01-01

A finite element method is introduced for solving the neutron transport equations. Our method falls into the category of Petrov-Galerkin solution, since the trial space differs from the test space. The close relationship between this method and the discrete ordinate method is discussed, and the methods are compared for simple test problems

3D unstructured mesh discontinuous finite element hydro

International Nuclear Information System (INIS)

Prasad, M.K.; Kershaw, D.S.; Shaw, M.J.

1995-01-01

The authors present detailed features of the ICF3D hydrodynamics code used for inertial fusion simulations. This code is intended to be a state-of-the-art upgrade of the well-known fluid code, LASNEX. ICF3D employs discontinuous finite elements on a discrete unstructured mesh consisting of a variety of 3D polyhedra including tetrahedra, prisms, and hexahedra. The authors discussed details of how the ROE-averaged second-order convection was applied on the discrete elements, and how the C++ coding interface has helped to simplify implementing the many physics and numerics modules within the code package. The author emphasized the virtues of object-oriented design in large scale projects such as ICF3D

Quadratic Finite Element Method for 1D Deterministic Transport

International Nuclear Information System (INIS)

Tolar, D R Jr.; Ferguson, J M

2004-01-01

In the discrete ordinates, or SN, numerical solution of the transport equation, both the spatial ((und r)) and angular ((und (Omega))) dependences on the angular flux ψ(und r),(und (Omega))are modeled discretely. While significant effort has been devoted toward improving the spatial discretization of the angular flux, we focus on improving the angular discretization of ψ(und r),(und (Omega)). Specifically, we employ a Petrov-Galerkin quadratic finite element approximation for the differencing of the angular variable (μ) in developing the one-dimensional (1D) spherical geometry S N equations. We develop an algorithm that shows faster convergence with angular resolution than conventional S N algorithms

A Kohn–Sham equation solver based on hexahedral finite elements

International Nuclear Information System (INIS)

Fang Jun; Gao Xingyu; Zhou Aihui

2012-01-01

We design a Kohn–Sham equation solver based on hexahedral finite element discretizations. The solver integrates three schemes proposed in this paper. The first scheme arranges one a priori locally-refined hexahedral mesh with appropriate multiresolution. The second one is a modified mass-lumping procedure which accelerates the diagonalization in the self-consistent field iteration. The third one is a finite element recovery method which enhances the eigenpair approximations with small extra work. We carry out numerical tests on each scheme to investigate the validity and efficiency, and then apply them to calculate the ground state total energies of nanosystems C 60 , C 120 , and C 275 H 172 . It is shown that our solver appears to be computationally attractive for finite element applications in electronic structure study.

Application of the finite element method to the neutron transport equation

International Nuclear Information System (INIS)

Martin, W.R.

1976-01-01

This paper examines the theoretical and practical application of the finite element method to the neutron transport equation. It is shown that in principle the system of equations obtained by application of the finite element method can be solved with certain physical restrictions concerning the criticality of the medium. The convergence of this approximate solution to the exact solution with mesh refinement is examined, and a non-optical estimate of the convergence rate is obtained analytically. It is noted that the numerical results indicate a faster convergence rate and several approaches to obtain this result analytically are outlined. The practical application of the finite element method involved the development of a computer code capable of solving the neutron transport equation in 1-D plane geometry. Vacuum, reflecting, or specified incoming boundary conditions may be analyzed, and all are treated as natural boundary conditions. The time-dependent transport equation is also examined and it is shown that the application of the finite element method in conjunction with the Crank-Nicholson time discretization method results in a system of algebraic equations which is readily solved. Numerical results are given for several critical slab eigenvalue problems, including anisotropic scattering, and the results compare extremely well with benchmark results. It is seen that the finite element code is more efficient than a standard discrete ordinates code for certain problems. A problem with severe heterogeneities is considered and it is shown that the use of discontinuous spatial and angular elements results in a marked improvement in the results. Finally, time-dependent problems are examined and it is seen that the phenomenon of angular mode separation makes the numerical treatment of the transport equation in slab geometry a considerable challenge, with the result that the angular mesh has a dominant effect on obtaining acceptable solutions

Finite element analysis of multilayer coextrusion.

Energy Technology Data Exchange (ETDEWEB)

Hopkins, Matthew Morgan; Schunk, Peter Randall; Baer, Thomas A. (Proctor & Gamble Company, West Chester, OH); Mrozek, Randy A. (Army Research Laboratory, Adelphi, MD); Lenhart, Joseph Ludlow (Army Research Laboratory, Adelphi, MD); Rao, Rekha Ranjana; Collins, Robert (Oak Ridge National Laboratory); Mondy, Lisa Ann

2011-09-01

Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.

A generalized nodal finite element formalism for discrete ordinates equations in slab geometry Part I: Theory in the continuous moment case

International Nuclear Information System (INIS)

Hennart, J.P.; Valle, E. del.

1995-01-01

A generalized nodal finite element formalism is presented, which covers virtually all known finit difference approximation to the discrete ordinates equations in slab geometry. This paper (Part 1) presents the theory of the so called open-quotes continuous moment methodsclose quotes, which include such well-known methods as the open-quotes diamond differenceclose quotes and the open-quotes characteristicclose quotes schemes. In a second paper (hereafter referred to as Part II), the authors will present the theory of the open-quotes discontinuous moment methodsclose quotes, consisting in particular of the open-quotes linear discontinuousclose quotes scheme as well as of an entire new class of schemes. Corresponding numerical results are available for all these schemes and will be presented in a third paper (Part III). 12 refs

Lagrangian Finite-Element Method for the Simulation of K-BKZ Fluids with Third Order Accuracy

DEFF Research Database (Denmark)

Marin, José Manuel Román; Rasmussen, Henrik K.

2009-01-01

system attached to the particles is discretized by ten-node quadratic tetrahedral elements using Cartesian coordinates and the pressure by linear interpolation inside these elements. The spatial discretization of the governing equations follows the mixed Galerkin finite element method. The time integral...... is discretized by a quadratic interpolation in time. The convergence of the method in time and space was demonstrated on the free surface problem of a filament stretched between two plates, considering the axisymmetric case as well as the growth of non-axisymmetric disturbances on the free surface. The scheme...

Deflation in preconditioned conjugate gradient methods for Finite Element Problems

NARCIS (Netherlands)

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

The finite volume element (FVE) and multigrid method for the incompressible Navier-Stokes equations

International Nuclear Information System (INIS)

Gu Lizhen; Bao Weizhu

1992-01-01

The authors apply FVE method to discrete INS equations with the original variable, in which the bilinear square finite element and the square finite volume are chosen. The discrete schemes of INS equations are presented. The FMV multigrid algorithm is applied to solve that discrete system, where DGS iteration is used as smoother, DGS distributive mode for the INS discrete system is also presented. The sample problems for the square cavity flow with Reynolds number Re≤100 are successfully calculated. The numerical solutions show that the results with 1 FMV is satisfactory and when Re is not large, The FVE discrete scheme of the conservative INS equations and that of non-conservative INS equations with linearization both can provide almost same accuracy

Finite element modelling of the mechanics of discrete carbon nanotubes filled with ZnS and comparison with experimental observations

KAUST Repository

Monteiro, André O.; Da Costa, Pedro M. F. J.; Cachim, Paulo Barreto

2013-01-01

The mechanical response to a uniaxial compressive force of a single carbon nanotube (CNT) filled (or partially-filled) with ZnS has been modelled. A semi-empirical approach based on the finite element method was used whereby modelling outcomes were closely matched to experimental observations. This is the first example of the use of the continuum approach to model the mechanical behaviour of discrete filled CNTs. In contrast to more computationally demanding methods such as density functional theory or molecular dynamics, our approach provides a viable and expedite alternative to model the mechanics of filled multi-walled CNTs. © 2013 Springer Science+Business Media New York.

Finite element modelling of the mechanics of discrete carbon nanotubes filled with ZnS and comparison with experimental observations

KAUST Repository

Monteiro, André O.

2013-09-25

The mechanical response to a uniaxial compressive force of a single carbon nanotube (CNT) filled (or partially-filled) with ZnS has been modelled. A semi-empirical approach based on the finite element method was used whereby modelling outcomes were closely matched to experimental observations. This is the first example of the use of the continuum approach to model the mechanical behaviour of discrete filled CNTs. In contrast to more computationally demanding methods such as density functional theory or molecular dynamics, our approach provides a viable and expedite alternative to model the mechanics of filled multi-walled CNTs. © 2013 Springer Science+Business Media New York.

The discontinuous finite element method for solving Eigenvalue problems of transport equations

International Nuclear Information System (INIS)

Yang, Shulin; Wang, Ruihong

2011-01-01

In this paper, the multigroup transport equations for solving the eigenvalues λ and K_e_f_f under two dimensional cylindrical coordinate are discussed. Aimed at the equations, the discretizing way combining discontinuous finite element method (DFE) with discrete ordinate method (SN) is developed, and the iterative algorithms and steps are studied. The numerical results show that the algorithms are efficient. (author)

A study of unstable rock failures using finite difference and discrete element methods

Science.gov (United States)

Garvey, Ryan J.

Case histories in mining have long described pillars or faces of rock failing violently with an accompanying rapid ejection of debris and broken material into the working areas of the mine. These unstable failures have resulted in large losses of life and collapses of entire mine panels. Modern mining operations take significant steps to reduce the likelihood of unstable failure, however eliminating their occurrence is difficult in practice. Researchers over several decades have supplemented studies of unstable failures through the application of various numerical methods. The direction of the current research is to extend these methods and to develop improved numerical tools with which to study unstable failures in underground mining layouts. An extensive study is first conducted on the expression of unstable failure in discrete element and finite difference methods. Simulated uniaxial compressive strength tests are run on brittle rock specimens. Stable or unstable loading conditions are applied onto the brittle specimens by a pair of elastic platens with ranging stiffnesses. Determinations of instability are established through stress and strain histories taken for the specimen and the system. Additional numerical tools are then developed for the finite difference method to analyze unstable failure in larger mine models. Instability identifiers are established for assessing the locations and relative magnitudes of unstable failure through measures of rapid dynamic motion. An energy balance is developed which calculates the excess energy released as a result of unstable equilibria in rock systems. These tools are validated through uniaxial and triaxial compressive strength tests and are extended to models of coal pillars and a simplified mining layout. The results of the finite difference simulations reveal that the instability identifiers and excess energy calculations provide a generalized methodology for assessing unstable failures within potentially complex

On the Robustness and Prospects of Adaptive BDDC Methods for Finite Element Discretizations of Elliptic PDEs with High-Contrast Coefficients

KAUST Repository

Zampini, Stefano

2016-06-02

Balancing Domain Decomposition by Constraints (BDDC) methods have proven to be powerful preconditioners for large and sparse linear systems arising from the finite element discretization of elliptic PDEs. Condition number bounds can be theoretically established that are independent of the number of subdomains of the decomposition. The core of the methods resides in the design of a larger and partially discontinuous finite element space that allows for fast application of the preconditioner, where Cholesky factorizations of the subdomain finite element problems are additively combined with a coarse, global solver. Multilevel and highly-scalable algorithms can be obtained by replacing the coarse Cholesky solver with a coarse BDDC preconditioner. BDDC methods have the remarkable ability to control the condition number, since the coarse space of the preconditioner can be adaptively enriched at the cost of solving local eigenproblems. The proper identification of these eigenproblems extends the robustness of the methods to any heterogeneity in the distribution of the coefficients of the PDEs, not only when the coefficients jumps align with the subdomain boundaries or when the high contrast regions are confined to lie in the interior of the subdomains. The specific adaptive technique considered in this paper does not depend upon any interaction of discretization and partition; it relies purely on algebraic operations. Coarse space adaptation in BDDC methods has attractive algorithmic properties, since the technique enhances the concurrency and the arithmetic intensity of the preconditioning step of the sparse implicit solver with the aim of controlling the number of iterations of the Krylov method in a black-box fashion, thus reducing the number of global synchronization steps and matrix vector multiplications needed by the iterative solver; data movement and memory bound kernels in the solve phase can be thus limited at the expense of extra local ops during the setup of

Finite-element discretization of 3D energy-transport equations for semiconductors

Energy Technology Data Exchange (ETDEWEB)

Gadau, Stephan

2007-07-01

In this thesis a mathematical model was derived that describes the charge and energy transport in semiconductor devices like transistors. Moreover, numerical simulations of these physical processes are performed. In order to accomplish this, methods of theoretical physics, functional analysis, numerical mathematics and computer programming are applied. After an introduction to the status quo of semiconductor device simulation methods and a brief review of historical facts up to now, the attention is shifted to the construction of a model, which serves as the basis of the subsequent derivations in the thesis. Thereby the starting point is an important equation of the theory of dilute gases. From this equation the model equations are derived and specified by means of a series expansion method. This is done in a multi-stage derivation process, which is mainly taken from a scientific paper and which does not constitute the focus of this thesis. In the following phase we specify the mathematical setting and make precise the model assumptions. Thereby we make use of methods of functional analysis. Since the equations we deal with are coupled, we are concerned with a nonstandard problem. In contrary, the theory of scalar elliptic equations is established meanwhile. Subsequently, we are preoccupied with the numerical discretization of the equations. A special finite-element method is used for the discretization. This special approach has to be done in order to make the numerical results appropriate for practical application. By a series of transformations from the discrete model we derive a system of algebraic equations that are eligible for numerical evaluation. Using self-made computer programs we solve the equations to get approximate solutions. These programs are based on new and specialized iteration procedures that are developed and thoroughly tested within the frame of this research work. Due to their importance and their novel status, they are explained and

Coupling multipoint flux mixed finite element methodswith continuous Galerkin methods for poroelasticity

KAUST Repository

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.

Dual-mixed finite elements for the three-field Stokes model as a finite volume method on staggered grids

KAUST Repository

Kou, Jisheng

2017-06-09

In this paper, a new three-field weak formulation for Stokes problems is developed, and from this, a dual-mixed finite element method is proposed on a rectangular mesh. In the proposed mixed methods, the components of stress tensor are approximated by piecewise constant functions or Q1 functions, while the velocity and pressure are discretized by the lowest-order Raviart-Thomas element and the piecewise constant functions, respectively. Using quadrature rules, we demonstrate that this scheme can be reduced into a finite volume method on staggered grid, which is extensively used in computational fluid mechanics and engineering.

Discrete elements method of neutron transport

International Nuclear Information System (INIS)

Mathews, K.A.

1988-01-01

In this paper a new neutron transport method, called discrete elements (L N ) is derived and compared to discrete ordinates methods, theoretically and by numerical experimentation. The discrete elements method is based on discretizing the Boltzmann equation over a set of elements of angle. The discrete elements method is shown to be more cost-effective than discrete ordinates, in terms of accuracy versus execution time and storage, for the cases tested. In a two-dimensional test case, a vacuum duct in a shield, the L N method is more consistently convergent toward a Monte Carlo benchmark solution

Coupling multipoint flux mixed finite element methodswith continuous Galerkin methods for poroelasticity

KAUST Repository

Wheeler, Mary; Xue, Guangri; Yotov, Ivan

2013-01-01

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

A (Dis)continuous finite element model for generalized 2D vorticity dynamics

NARCIS (Netherlands)

Bernsen, E.; Bokhove, Onno; van der Vegt, Jacobus J.W.

2005-01-01

A mixed continuous and discontinuous Galerkin finite element discretization is constructed for a generalized vorticity streamfunction formulation in two spatial dimensions. This formulation consists of a hyperbolic (potential) vorticity equation and a linear elliptic equation for a (transport)

Basic Finite Element Method

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.

Finite element method for radiation heat transfer in multi-dimensional graded index medium

International Nuclear Information System (INIS)

Liu, L.H.; Zhang, L.; Tan, H.P.

2006-01-01

In graded index medium, ray goes along a curved path determined by Fermat principle, and curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectories, a finite element method based on discrete ordinate equation is developed to solve the radiative transfer problem in a multi-dimensional semitransparent graded index medium. Two particular test problems of radiative transfer are taken as examples to verify this finite element method. The predicted dimensionless net radiative heat fluxes are determined by the proposed method and compared with the results obtained by finite volume method. The results show that the finite element method presented in this paper has a good accuracy in solving the multi-dimensional radiative transfer problem in semitransparent graded index medium

Assembly of finite element methods on graphics processors

KAUST Repository

Cecka, Cris

2010-08-23

Recently, graphics processing units (GPUs) have had great success in accelerating many numerical computations. We present their application to computations on unstructured meshes such as those in finite element methods. Multiple approaches in assembling and solving sparse linear systems with NVIDIA GPUs and the Compute Unified Device Architecture (CUDA) are created and analyzed. Multiple strategies for efficient use of global, shared, and local memory, methods to achieve memory coalescing, and optimal choice of parameters are introduced. We find that with appropriate preprocessing and arrangement of support data, the GPU coprocessor using single-precision arithmetic achieves speedups of 30 or more in comparison to a well optimized double-precision single core implementation. We also find that the optimal assembly strategy depends on the order of polynomials used in the finite element discretization. © 2010 John Wiley & Sons, Ltd.

Computation of 2D compressible flows with a finite element method

International Nuclear Information System (INIS)

Montagne, J.L.

1981-04-01

When the homogeneous modelisation of the two phase flow is used the set of equations describing the flow is similar to an Euler system. Mixed finite elements are appropriate to discretize the equations. First, main properties of this kind of elements are reminded. Then, some properties of semi-implicite schemes on stability and entropy are given. Numerical tests have been performed, and the scheme gave satisfactory results

Nonparametric Identification and Estimation of Finite Mixture Models of Dynamic Discrete Choices

OpenAIRE

Hiroyuki Kasahara; Katsumi Shimotsu

2006-01-01

In dynamic discrete choice analysis, controlling for unobserved heterogeneity is an important issue, and finite mixture models provide flexible ways to account for unobserved heterogeneity. This paper studies nonparametric identifiability of type probabilities and type-specific component distributions in finite mixture models of dynamic discrete choices. We derive sufficient conditions for nonparametric identification for various finite mixture models of dynamic discrete choices used in appli...

A mixed finite element domain decomposition method for nearly elastic wave equations in the frequency domain

Energy Technology Data Exchange (ETDEWEB)

Feng, Xiaobing [Univ. of Tennessee, Knoxville, TN (United States)

1996-12-31

A non-overlapping domain decomposition iterative method is proposed and analyzed for mixed finite element methods for a sequence of noncoercive elliptic systems with radiation boundary conditions. These differential systems describe the motion of a nearly elastic solid in the frequency domain. The convergence of the iterative procedure is demonstrated and the rate of convergence is derived for the case when the domain is decomposed into subdomains in which each subdomain consists of an individual element associated with the mixed finite elements. The hybridization of mixed finite element methods plays a important role in the construction of the discrete procedure.

Error Analysis of a Finite Element Method for the Space-Fractional Parabolic Equation

KAUST Repository

Jin, Bangti; Lazarov, Raytcho; Pasciak, Joseph; Zhou, Zhi

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.

An implicit finite element method for discrete dynamic fracture

Energy Technology Data Exchange (ETDEWEB)

Gerken, Jobie M. [Colorado State Univ., Fort Collins, CO (United States)

1999-12-01

A method for modeling the discrete fracture of two-dimensional linear elastic structures with a distribution of small cracks subject to dynamic conditions has been developed. The foundation for this numerical model is a plane element formulated from the Hu-Washizu energy principle. The distribution of small cracks is incorporated into the numerical model by including a small crack at each element interface. The additional strain field in an element adjacent to this crack is treated as an externally applied strain field in the Hu-Washizu energy principle. The resulting stiffness matrix is that of a standard plane element. The resulting load vector is that of a standard plane element with an additional term that includes the externally applied strain field. Except for the crack strain field equations, all terms of the stiffness matrix and load vector are integrated symbolically in Maple V so that fully integrated plane stress and plane strain elements are constructed. The crack strain field equations are integrated numerically. The modeling of dynamic behavior of simple structures was demonstrated within acceptable engineering accuracy. In the model of axial and transverse vibration of a beam and the breathing mode of vibration of a thin ring, the dynamic characteristics were shown to be within expected limits. The models dominated by tensile forces (the axially loaded beam and the pressurized ring) were within 0.5% of the theoretical values while the shear dominated model (the transversely loaded beam) is within 5% of the calculated theoretical value. The constant strain field of the tensile problems can be modeled exactly by the numerical model. The numerical results should therefore, be exact. The discrepancies can be accounted for by errors in the calculation of frequency from the numerical results. The linear strain field of the transverse model must be modeled by a series of constant strain elements. This is an approximation to the true strain field, so some

A spatial discretization of the MHD equations based on the finite volume - spectral method

International Nuclear Information System (INIS)

Miyoshi, Takahiro

2000-05-01

Based on the finite volume - spectral method, we present new discretization formulae for the spatial differential operators in the full system of the compressible MHD equations. In this approach, the cell-centered finite volume method is adopted in a bounded plane (poloidal plane), while the spectral method is applied to the differential with respect to the periodic direction perpendicular to the poloidal plane (toroidal direction). Here, an unstructured grid system composed of the arbitrary triangular elements is utilized for constructing the cell-centered finite volume method. In order to maintain the divergence free constraint of the magnetic field numerically, only the poloidal component of the rotation is defined at three edges of the triangular element. This poloidal component is evaluated under the assumption that the toroidal component of the operated vector times the radius, RA φ , is linearly distributed in the element. The present method will be applied to the nonlinear MHD dynamics in an realistic torus geometry without the numerical singularities. (author)

Least-squares finite element discretizations of neutron transport equations in 3 dimensions

Energy Technology Data Exchange (ETDEWEB)

Manteuffel, T.A [Univ. of Colorado, Boulder, CO (United States); Ressel, K.J. [Interdisciplinary Project Center for Supercomputing, Zurich (Switzerland); Starkes, G. [Universtaet Karlsruhe (Germany)

1996-12-31

The least-squares finite element framework to the neutron transport equation introduced in is based on the minimization of a least-squares functional applied to the properly scaled neutron transport equation. Here we report on some practical aspects of this approach for neutron transport calculations in three space dimensions. The systems of partial differential equations resulting from a P{sub 1} and P{sub 2} approximation of the angular dependence are derived. In the diffusive limit, the system is essentially a Poisson equation for zeroth moment and has a divergence structure for the set of moments of order 1. One of the key features of the least-squares approach is that it produces a posteriori error bounds. We report on the numerical results obtained for the minimum of the least-squares functional augmented by an additional boundary term using trilinear finite elements on a uniform tesselation into cubes.

Robust mixed finite element methods to deal with incompressibility in finite strain in an industrial framework

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)

Navier-Stokes equations by the finite element method

International Nuclear Information System (INIS)

Portella, P.E.

1984-01-01

A computer program to solve the Navier-Stokes equations by using the Finite Element Method is implemented. The solutions variables investigated are stream-function/vorticity in the steady case and velocity/pressure in the steady state and transient cases. For steady state flow the equations are solved simultaneously by the Newton-Raphson method. For the time dependent formulation, a fractional step method is employed to discretize in time and artificial viscosity is used to preclude spurious oscilations in the solution. The element used is the three node triangle. Some numerical examples are presented and comparisons are made with applications already existent. (Author) [pt

2D deterministic radiation transport with the discontinuous finite element method

International Nuclear Information System (INIS)

Kershaw, D.; Harte, J.

1993-01-01

This report provides a complete description of the analytic and discretized equations for 2D deterministic radiation transport. This computational model has been checked against a wide variety of analytic test problems and found to give excellent results. We make extensive use of the discontinuous finite element method

Hybrid finite volume/ finite element method for radiative heat transfer in graded index media

Science.gov (United States)

Zhang, L.; Zhao, J. M.; Liu, L. H.; Wang, S. Y.

2012-09-01

The rays propagate along curved path determined by the Fermat principle in the graded index medium. The radiative transfer equation in graded index medium (GRTE) contains two specific redistribution terms (with partial derivatives to the angular coordinates) accounting for the effect of the curved ray path. In this paper, the hybrid finite volume with finite element method (hybrid FVM/FEM) (P.J. Coelho, J. Quant. Spectrosc. Radiat. Transf., vol. 93, pp. 89-101, 2005) is extended to solve the radiative heat transfer in two-dimensional absorbing-emitting-scattering graded index media, in which the spatial discretization is carried out using a FVM, while the angular discretization is by a FEM. The FEM angular discretization is demonstrated to be preferable in dealing with the redistribution terms in the GRTE. Two stiff matrix assembly schemes of the angular FEM discretization, namely, the traditional assembly approach and a new spherical assembly approach (assembly on the unit sphere of the solid angular space), are discussed. The spherical assembly scheme is demonstrated to give better results than the traditional assembly approach. The predicted heat flux distributions and temperature distributions in radiative equilibrium are determined by the proposed method and compared with the results available in other references. The proposed hybrid FVM/FEM method can predict the radiative heat transfer in absorbing-emitting-scattering graded index medium with good accuracy.

Hybrid finite volume/ finite element method for radiative heat transfer in graded index media

International Nuclear Information System (INIS)

Zhang, L.; Zhao, J.M.; Liu, L.H.; Wang, S.Y.

2012-01-01

The rays propagate along curved path determined by the Fermat principle in the graded index medium. The radiative transfer equation in graded index medium (GRTE) contains two specific redistribution terms (with partial derivatives to the angular coordinates) accounting for the effect of the curved ray path. In this paper, the hybrid finite volume with finite element method (hybrid FVM/FEM) (P.J. Coelho, J. Quant. Spectrosc. Radiat. Transf., vol. 93, pp. 89-101, 2005) is extended to solve the radiative heat transfer in two-dimensional absorbing-emitting-scattering graded index media, in which the spatial discretization is carried out using a FVM, while the angular discretization is by a FEM. The FEM angular discretization is demonstrated to be preferable in dealing with the redistribution terms in the GRTE. Two stiff matrix assembly schemes of the angular FEM discretization, namely, the traditional assembly approach and a new spherical assembly approach (assembly on the unit sphere of the solid angular space), are discussed. The spherical assembly scheme is demonstrated to give better results than the traditional assembly approach. The predicted heat flux distributions and temperature distributions in radiative equilibrium are determined by the proposed method and compared with the results available in other references. The proposed hybrid FVM/FEM method can predict the radiative heat transfer in absorbing-emitting-scattering graded index medium with good accuracy.

Explicit Dynamic Finite Element Method for Predicting Implosion/Explosion Induced Failure of Shell Structures

Directory of Open Access Journals (Sweden)

Jeong-Hoon Song

2013-01-01

Full Text Available A simplified implementation of the conventional extended finite element method (XFEM for dynamic fracture in thin shells is presented. Though this implementation uses the same linear combination of the conventional XFEM, it allows for considerable simplifications of the discontinuous displacement and velocity fields in shell finite elements. The proposed method is implemented for the discrete Kirchhoff triangular (DKT shell element, which is one of the most popular shell elements in engineering analysis. Numerical examples for dynamic failure of shells under impulsive loads including implosion and explosion are presented to demonstrate the effectiveness and robustness of the method.

Finite Element Method Based Modeling of Resistance Spot-Welded Mild Steel

Directory of Open Access Journals (Sweden)

Miloud Zaoui

Full Text Available Abstract This paper deals with Finite Element refined and simplified models of a mild steel spot-welded specimen, developed and validated based on quasi-static cross-tensile experimental tests. The first model was constructed with a fine discretization of the metal sheet and the spot weld was defined as a special geometric zone of the specimen. This model provided, in combination with experimental tests, the input data for the development of the second model, which was constructed with respect to the mesh size used in the complete car finite element model. This simplified model was developed with coarse shell elements and a spring-type beam element was used to model the spot weld behavior. The global accuracy of the two models was checked by comparing simulated and experimental load-displacement curves and by studying the specimen deformed shapes and the plastic deformation growth in the metal sheets. The obtained results show that both fine and coarse finite element models permit a good prediction of the experimental tests.

Multigrid Finite Element Method in Calculation of 3D Homogeneous and Composite Solids

Directory of Open Access Journals (Sweden)

A.D. Matveev

2016-12-01

Full Text Available In the present paper, a method of multigrid finite elements to calculate elastic three-dimensional homogeneous and composite solids under static loading has been suggested. The method has been developed based on the finite element method algorithms using homogeneous and composite three-dimensional multigrid finite elements (MFE. The procedures for construction of MFE of both rectangular parallelepiped and complex shapes have been shown. The advantages of MFE are that they take into account, following the rules of the microapproach, heterogeneous and microhomogeneous structures of the bodies, describe the three-dimensional stress-strain state (without any simplifying hypotheses in homogeneous and composite solids, as well as generate small dimensional discrete models and numerical solutions with a high accuracy.

Choice of input fields in stochastic finite elements

DEFF Research Database (Denmark)

Ditlevsen, Ove Dalager; Tarp-Johansen, Niels Jacob

1996-01-01

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

Choice of input fields in stochastic finite elements

DEFF Research Database (Denmark)

Ditlevsen, Ove Dalager; Tarp-Johansen, Niels Jacob

1999-01-01

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

Sensitivity analysis of the Galerkin finite element method neutron diffusion solver to the shape of the elements

Energy Technology Data Exchange (ETDEWEB)

Hosseini, Seyed Abolfaz [Dept. of Energy Engineering, Sharif University of Technology, Tehran (Iran, Islamic Republic of)

2017-02-15

The purpose of the present study is the presentation of the appropriate element and shape function in the solution of the neutron diffusion equation in two-dimensional (2D) geometries. To this end, the multigroup neutron diffusion equation is solved using the Galerkin finite element method in both rectangular and hexagonal reactor cores. The spatial discretization of the equation is performed using unstructured triangular and quadrilateral finite elements. Calculations are performed using both linear and quadratic approximations of shape function in the Galerkin finite element method, based on which results are compared. Using the power iteration method, the neutron flux distributions with the corresponding eigenvalue are obtained. The results are then validated against the valid results for IAEA-2D and BIBLIS-2D benchmark problems. To investigate the dependency of the results to the type and number of the elements, and shape function order, a sensitivity analysis of the calculations to the mentioned parameters is performed. It is shown that the triangular elements and second order of the shape function in each element give the best results in comparison to the other states.

Adaptive mixed finite element methods for Darcy flow in fractured porous media

KAUST Repository

Chen, Huangxin; Salama, Amgad; Sun, Shuyu

2016-01-01

In this paper, we propose adaptive mixed finite element methods for simulating the single-phase Darcy flow in two-dimensional fractured porous media. The reduced model that we use for the simulation is a discrete fracture model coupling Darcy flows in the matrix and the fractures, and the fractures are modeled by one-dimensional entities. The Raviart-Thomas mixed finite element methods are utilized for the solution of the coupled Darcy flows in the matrix and the fractures. In order to improve the efficiency of the simulation, we use adaptive mixed finite element methods based on novel residual-based a posteriori error estimators. In addition, we develop an efficient upscaling algorithm to compute the effective permeability of the fractured porous media. Several interesting examples of Darcy flow in the fractured porous media are presented to demonstrate the robustness of the algorithm.

Adaptive mixed finite element methods for Darcy flow in fractured porous media

KAUST Repository

Chen, Huangxin

2016-09-21

In this paper, we propose adaptive mixed finite element methods for simulating the single-phase Darcy flow in two-dimensional fractured porous media. The reduced model that we use for the simulation is a discrete fracture model coupling Darcy flows in the matrix and the fractures, and the fractures are modeled by one-dimensional entities. The Raviart-Thomas mixed finite element methods are utilized for the solution of the coupled Darcy flows in the matrix and the fractures. In order to improve the efficiency of the simulation, we use adaptive mixed finite element methods based on novel residual-based a posteriori error estimators. In addition, we develop an efficient upscaling algorithm to compute the effective permeability of the fractured porous media. Several interesting examples of Darcy flow in the fractured porous media are presented to demonstrate the robustness of the algorithm.

Generalized finite elements

International Nuclear Information System (INIS)

Wachspress, E.

2009-01-01

Triangles and rectangles are the ubiquitous elements in finite element studies. Only these elements admit polynomial basis functions. Rational functions provide a basis for elements having any number of straight and curved sides. Numerical complexities initially associated with rational bases precluded extensive use. Recent analysis has reduced these difficulties and programs have been written to illustrate effectiveness. Although incorporation in major finite element software requires considerable effort, there are advantages in some applications which warrant implementation. An outline of the basic theory and of recent innovations is presented here. (authors)

Superconvergence of mixed finite element approximations to 3-D Maxwell's equations in metamaterials

KAUST Repository

Huang, Yunqing

2011-09-01

Numerical simulation of metamaterials has attracted more and more attention since 2000, after the first metamaterial with negative refraction index was successfully constructed. In this paper we construct a fully-discrete leap-frog type finite element scheme to solve the three-dimensional time-dependent Maxwell\\'s equations when metamaterials are involved. First, we obtain some superclose results between the interpolations of the analytical solutions and finite element solutions obtained using arbitrary orders of Raviart-Thomas-Nédélec mixed spaces on regular cubic meshes. Then we prove the superconvergence result in the discrete l2 norm achieved for the lowest-order Raviart-Thomas-Nédélec space. To our best knowledge, such superconvergence results have never been obtained elsewhere. Finally, we implement the leap-frog scheme and present numerical results justifying our theoretical analysis. © 2011 Elsevier Inc.

Superconvergence of mixed finite element approximations to 3-D Maxwell's equations in metamaterials

KAUST Repository

Huang, Yunqing; Li, Jichun; Yang, Wei; Sun, Shuyu

2011-01-01

Numerical simulation of metamaterials has attracted more and more attention since 2000, after the first metamaterial with negative refraction index was successfully constructed. In this paper we construct a fully-discrete leap-frog type finite element scheme to solve the three-dimensional time-dependent Maxwell's equations when metamaterials are involved. First, we obtain some superclose results between the interpolations of the analytical solutions and finite element solutions obtained using arbitrary orders of Raviart-Thomas-Nédélec mixed spaces on regular cubic meshes. Then we prove the superconvergence result in the discrete l2 norm achieved for the lowest-order Raviart-Thomas-Nédélec space. To our best knowledge, such superconvergence results have never been obtained elsewhere. Finally, we implement the leap-frog scheme and present numerical results justifying our theoretical analysis. © 2011 Elsevier Inc.

Adaptive Finite Element-Discrete Element Analysis for Microseismic Modelling of Hydraulic Fracture Propagation of Perforation in Horizontal Well considering Pre-Existing Fractures

Directory of Open Access Journals (Sweden)

Yongliang Wang

2018-01-01

Full Text Available Hydrofracturing technology of perforated horizontal well has been widely used to stimulate the tight hydrocarbon reservoirs for gas production. To predict the hydraulic fracture propagation, the microseismicity can be used to infer hydraulic fractures state; by the effective numerical methods, microseismic events can be addressed from changes of the computed stresses. In numerical models, due to the challenges in accurately representing the complex structure of naturally fractured reservoir, the interaction between hydraulic and pre-existing fractures has not yet been considered and handled satisfactorily. To overcome these challenges, the adaptive finite element-discrete element method is used to refine mesh, effectively identify the fractures propagation, and investigate microseismic modelling. Numerical models are composed of hydraulic fractures, pre-existing fractures, and microscale pores, and the seepage analysis based on the Darcy’s law is used to determine fluid flow; then moment tensors in microseismicity are computed based on the computed stresses. Unfractured and naturally fractured models are compared to assess the influences of pre-existing fractures on hydrofracturing. The damaged and contact slip events were detected by the magnitudes, B-values, Hudson source type plots, and focal spheres.

Introduction to assembly of finite element methods on graphics processors

International Nuclear Information System (INIS)

Cecka, Cristopher; Lew, Adrian; Darve, Eric

2010-01-01

Recently, graphics processing units (GPUs) have had great success in accelerating numerical computations. We present their application to computations on unstructured meshes such as those in finite element methods. Multiple approaches in assembling and solving sparse linear systems with NVIDIA GPUs and the Compute Unified Device Architecture (CUDA) are presented and discussed. Multiple strategies for efficient use of global, shared, and local memory, methods to achieve memory coalescing, and optimal choice of parameters are introduced. We find that with appropriate preprocessing and arrangement of support data, the GPU coprocessor achieves speedups of 30x or more in comparison to a well optimized serial implementation on the CPU. We also find that the optimal assembly strategy depends on the order of polynomials used in the finite-element discretization.

Straightened cervical lordosis causes stress concentration: a finite element model study

Energy Technology Data Exchange (ETDEWEB)

Wei, Wei; Shi, Shiyuan; Fei, Jun; Wang, Yifan; Chen, Chunyue [Hangzhou Red Cross Hospital, Hangzhou, Zhejiang, (China); Liao, Shenhui [School of Information Science and Engineering, Central South University, Changsha, Hunan (China)

2013-03-15

In this study, we propose a finite element analysis of the complete cervical spine with straightened and normal physiological curvature by using a specially designed modelling system. An accurate finite element model is established to recommend plausible approaches to treatment of cervical spondylosis through the finite element analysis results. There are few reports of biomechanics influence of the straightened cervical curve. It is difficult to measure internal responses of cervical spine directly. However, the finite element method has been reported to have the capability to quantify both external and internal responses to mechanical loading, such as the strain and stress distribution of spinal components. We choose a subject with a straightened cervical spine from whom to collect the CT scan data, which formed the basis of the finite element analysis. By using a specially designed modelling system, a high quality finite element model of the complete cervical spine with straightened curvature was generated, which was then mapped to reconstruct a normal physiological curvature model by a volumetric mesh deformation method based on discrete differential properties. Then, the same boundary conditions were applied to do a comparison. The result demonstrated that the active movement range of straightened cervical spine decreased by 24–33 %, but the stress increased by 5–95 %. The stress was concentrated at the facet joint cartilage, uncovertebral joint and the disk. The results suggest that cervical lordosis may have a direct impact on cervical spondylosis treatment. These results may be useful for clinical treatment of cervical spondylosis with straightened curvature.

Straightened cervical lordosis causes stress concentration: a finite element model study

International Nuclear Information System (INIS)

Wei, Wei; Shi, Shiyuan; Fei, Jun; Wang, Yifan; Chen, Chunyue; Liao, Shenhui

2013-01-01

In this study, we propose a finite element analysis of the complete cervical spine with straightened and normal physiological curvature by using a specially designed modelling system. An accurate finite element model is established to recommend plausible approaches to treatment of cervical spondylosis through the finite element analysis results. There are few reports of biomechanics influence of the straightened cervical curve. It is difficult to measure internal responses of cervical spine directly. However, the finite element method has been reported to have the capability to quantify both external and internal responses to mechanical loading, such as the strain and stress distribution of spinal components. We choose a subject with a straightened cervical spine from whom to collect the CT scan data, which formed the basis of the finite element analysis. By using a specially designed modelling system, a high quality finite element model of the complete cervical spine with straightened curvature was generated, which was then mapped to reconstruct a normal physiological curvature model by a volumetric mesh deformation method based on discrete differential properties. Then, the same boundary conditions were applied to do a comparison. The result demonstrated that the active movement range of straightened cervical spine decreased by 24–33 %, but the stress increased by 5–95 %. The stress was concentrated at the facet joint cartilage, uncovertebral joint and the disk. The results suggest that cervical lordosis may have a direct impact on cervical spondylosis treatment. These results may be useful for clinical treatment of cervical spondylosis with straightened curvature.

On using moving windows in finite element time domain simulation for long accelerator structures

International Nuclear Information System (INIS)

Lee, L.-Q.; Candel, Arno; Ng, Cho; Ko, Kwok

2010-01-01

A finite element moving window technique is developed to simulate the propagation of electromagnetic waves induced by the transit of a charged particle beam inside large and long structures. The window moving along with the beam in the computational domain adopts high-order finite element basis functions through p refinement and/or a high-resolution mesh through h refinement so that a sufficient accuracy is attained with substantially reduced computational costs. Algorithms to transfer discretized fields from one mesh to another, which are the keys to implementing a moving window in a finite element unstructured mesh, are presented. Numerical experiments are carried out using the moving window technique to compute short-range wakefields in long accelerator structures. The results are compared with those obtained from the normal finite element time domain (FETD) method and the advantages of using the moving window technique are discussed.

SURF: a subroutine code to draw the axonometric projection of a surface generated by a scalar function over a discretized plane domain using finite element computations

International Nuclear Information System (INIS)

Giuliani, Giovanni; Giuliani, Silvano.

1980-01-01

The FORTRAN IV subroutine SURF has been designed to help visualising the results of Finite Element computations. It drawns the axonometric projection of a surface generated in 3-dimensional space by a scalar function over a discretized plane domain. The most important characteristic of the routine is to remove the hidden lines and in this way it enables a clear vision of the details of the generated surface

The Galerkin finite element method for a multi-term time-fractional diffusion equation

KAUST Repository

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.

Finite cover method with mortar elements for elastoplasticity problems

Science.gov (United States)

Kurumatani, M.; Terada, K.

2005-06-01

Finite cover method (FCM) is extended to elastoplasticity problems. The FCM, which was originally developed under the name of manifold method, has recently been recognized as one of the generalized versions of finite element methods (FEM). Since the mesh for the FCM can be regular and squared regardless of the geometry of structures to be analyzed, structural analysts are released from a burdensome task of generating meshes conforming to physical boundaries. Numerical experiments are carried out to assess the performance of the FCM with such discretization in elastoplasticity problems. Particularly to achieve this accurately, the so-called mortar elements are introduced to impose displacement boundary conditions on the essential boundaries, and displacement compatibility conditions on material interfaces of two-phase materials or on joint surfaces between mutually incompatible meshes. The validity of the mortar approximation is also demonstrated in the elastic-plastic FCM.

A strongly conservative finite element method for the coupling of Stokes and Darcy flow

KAUST Repository

Kanschat, G.; Riviè re, B.

2010-01-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

TRIP: a finite element computer program for the solution of convection heat transfer problems

International Nuclear Information System (INIS)

Slagter, W.; Roodbergen, H.A.

1976-01-01

The theory and use of the finite element code TRIP are described. The code calculates temperature distributions in three-dimensional continua subjected to convection heat transfer. A variational principle for transport phenomena is applied to solve the convection heat transfer problem with temperature and heat flux boundary conditions. The finite element discretization technique is used to reduce the continuous spatial solution into a finite number of unknowns. The method is developed in detail to determine temperature distributions in coolant passages of fuel rod bundles which are idealized by hexahedral elements. The development of the TRIP code is discussed and the listing of the program is given in FORTRAN IV. An example is given to illustrate the validity and practicality of the method

Proposal for element size and time increment selection guideline by 3-D finite element method for elastic waves propagation analysis

International Nuclear Information System (INIS)

Ishida, Hitoshi; Meshii, Toshiyuki

2008-01-01

This paper proposes a guideline for selection of element size and time increment by 3-D finite element method, which is applied to elastic wave propagation analysis for a long distance of a large structure. An element size and a time increment are determined by quantitative evaluation of strain, which must be 0 on the analysis model with a uniform motion, caused by spatial and time discretization. (author)

A stabilized second-order time accurate finite element formulation for incompressible viscous flow with heat transfer

International Nuclear Information System (INIS)

Curi, Marcos Filardy

2011-01-01

In view of the problem of global warming and the search for clean energy sources, a worldwide expansion on the use of nuclear energy is foreseen. Thus, the development of science and technology regarding nuclear power plants is essential, in particular in the field of reactor engineering. Fluid mechanics and heat transfer play an important role in the development of nuclear reactors. Computational Fluid Mechanics (CFD) is becoming ever more important in the optimization of cost and safety of the designs. This work presents a stabilized second-order time accurate finite element formulation for incompressible flows with heat transfer. A second order time discretization precedes a spatial discretization using finite elements. The terms that stabilize the finite element method arise naturally from the discretization process, rather than being introduced a priori in the variational formulation. The method was implemented in the program 'ns n ew s olvec2d av 2 M PI' written in FORTRAN90, developed in the Parallel Computing Laboratory at the Institute of Nuclear Engineering (LCP/IEN). Numerical solutions of some representative examples, including free, mixed and forced convection, demonstrate that the proposed stabilized formulation attains very good agreement with experimental and computational results available in the literature. (author)

Computational contact and impact mechanics fundamentals of modeling interfacial phenomena in nonlinear finite element analysis

CERN Document Server

Laursen, Tod A

2003-01-01

This book comprehensively treats the formulation and finite element approximation of contact and impact problems in nonlinear mechanics. Intended for students, researchers and practitioners interested in numerical solid and structural analysis, as well as for engineers and scientists dealing with technologies in which tribological response must be characterized, the book includes an introductory but detailed overview of nonlinear finite element formulations before dealing with contact and impact specifically. Topics encompassed include the continuum mechanics, mathematical structure, variational framework, and finite element implementations associated with contact/impact interaction. Additionally, important and currently emerging research topics in computational contact mechanics are introduced, encompassing such topics as tribological complexity, conservative treatment of inelastic impact interaction, and novel spatial discretization strategies.

Semianalytic Design Sensitivity Analysis of Nonlinear Structures With a Commercial Finite Element Package

International Nuclear Information System (INIS)

Lee, Tae Hee; Yoo, Jung Hun; Choi, Hyeong Cheol

2002-01-01

A finite element package is often used as a daily design tool for engineering designers in order to analyze and improve the design. The finite element analysis can provide the responses of a system for given design variables. Although finite element analysis can quite well provide the structural behaviors for given design variables, it cannot provide enough information to improve the design such as design sensitivity coefficients. Design sensitivity analysis is an essential step to predict the change in responses due to a change in design variables and to optimize a system with the aid of the gradient-based optimization techniques. To develop a numerical method of design sensitivity analysis, analytical derivatives that are based on analytical differentiation of the continuous or discrete finite element equations are effective but analytical derivatives are difficult because of the lack of internal information of the commercial finite element package such as shape functions. Therefore, design sensitivity analysis outside of the finite element package is necessary for practical application in an industrial setting. In this paper, the semi-analytic method for design sensitivity analysis is used for the development of the design sensitivity module outside of a commercial finite element package of ANSYS. The direct differentiation method is employed to compute the design derivatives of the response and the pseudo-load for design sensitivity analysis is effectively evaluated by using the design variation of the related internal nodal forces. Especially, we suggest an effective method for stress and nonlinear design sensitivity analyses that is independent of the commercial finite element package is also discussed. Numerical examples are illustrated to show the accuracy and efficiency of the developed method and to provide insights for implementation of the suggested method into other commercial finite element packages

Finite element solution of nonlinear eddy current problems with periodic excitation and its industrial applications.

Science.gov (United States)

Bíró, Oszkár; Koczka, Gergely; Preis, Kurt

2014-05-01

An efficient finite element method to take account of the nonlinearity of the magnetic materials when analyzing three-dimensional eddy current problems is presented in this paper. The problem is formulated in terms of vector and scalar potentials approximated by edge and node based finite element basis functions. The application of Galerkin techniques leads to a large, nonlinear system of ordinary differential equations in the time domain. The excitations are assumed to be time-periodic and the steady-state periodic solution is of interest only. This is represented either in the frequency domain as a finite Fourier series or in the time domain as a set of discrete time values within one period for each finite element degree of freedom. The former approach is the (continuous) harmonic balance method and, in the latter one, discrete Fourier transformation will be shown to lead to a discrete harmonic balance method. Due to the nonlinearity, all harmonics, both continuous and discrete, are coupled to each other. The harmonics would be decoupled if the problem were linear, therefore, a special nonlinear iteration technique, the fixed-point method is used to linearize the equations by selecting a time-independent permeability distribution, the so-called fixed-point permeability in each nonlinear iteration step. This leads to uncoupled harmonics within these steps. As industrial applications, analyses of large power transformers are presented. The first example is the computation of the electromagnetic field of a single-phase transformer in the time domain with the results compared to those obtained by traditional time-stepping techniques. In the second application, an advanced model of the same transformer is analyzed in the frequency domain by the harmonic balance method with the effect of the presence of higher harmonics on the losses investigated. Finally a third example tackles the case of direct current (DC) bias in the coils of a single-phase transformer.

Three dimensional finite element linear analysis of reinforced concrete structures

International Nuclear Information System (INIS)

Inbasakaran, M.; Pandarinathan, V.G.; Krishnamoorthy, C.S.

1979-01-01

A twenty noded isoparametric reinforced concrete solid element for the three dimensional linear elastic stress analysis of reinforced concrete structures is presented. The reinforcement is directly included as an integral part of the element thus facilitating discretization of the structure independent of the orientation of reinforcement. Concrete stiffness is evaluated by taking 3 x 3 x 3 Gauss integration rule and steel stiffness is evaluated numerically by considering three Gaussian points along the length of reinforcement. The numerical integration for steel stiffness necessiates the conversion of global coordiantes of the Gaussian points to nondimensional local coordinates and this is done by Newton Raphson iterative method. Subroutines for the above formulation have been developed and added to SAP and STAP routines for solving the examples. The validity of the reinforced concrete element is verified by comparison of results from finite element analysis and analytical results. It is concluded that this finite element model provides a valuable analytical tool for the three dimensional elastic stress analysis of concrete structures like beams curved in plan and nuclear containment vessels. (orig.)

Simulation of three-dimensional, time-dependent, incompressible flows by a finite element method

International Nuclear Information System (INIS)

Chan, S.T.; Gresho, P.M.; Lee, R.L.; Upson, C.D.

1981-01-01

A finite element model has been developed for simulating the dynamics of problems encountered in atmospheric pollution and safety assessment studies. The model is based on solving the set of three-dimensional, time-dependent, conservation equations governing incompressible flows. Spatial discretization is performed via a modified Galerkin finite element method, and time integration is carried out via the forward Euler method (pressure is computed implicitly, however). Several cost-effective techniques (including subcycling, mass lumping, and reduced Gauss-Legendre quadrature) which have been implemented are discussed. Numerical results are presented to demonstrate the applicability of the model

An h-adaptive finite element solver for the calculations of the electronic structures

International Nuclear Information System (INIS)

Bao Gang; Hu Guanghui; Liu Di

2012-01-01

In this paper, a framework of using h-adaptive finite element method for the Kohn–Sham equation on the tetrahedron mesh is presented. The Kohn–Sham equation is discretized by the finite element method, and the h-adaptive technique is adopted to optimize the accuracy and the efficiency of the algorithm. The locally optimal block preconditioned conjugate gradient method is employed for solving the generalized eigenvalue problem, and an algebraic multigrid preconditioner is used to accelerate the solver. A variety of numerical experiments demonstrate the effectiveness of our algorithm for both the all-electron and the pseudo-potential calculations.

A non-conformal finite element/finite volume scheme for the non-structured grid-based approximation of low Mach number flows

International Nuclear Information System (INIS)

Ansanay-Alex, G.

2009-01-01

The development of simulation codes aimed at a precise simulation of fires requires a precise approach of flame front phenomena by using very fine grids. The need to take different spatial scale into consideration leads to a local grid refinement and to a discretization with homogeneous grid for computing time and memory purposes. The author reports the approximation of the non-linear convection term, the scalar advection-diffusion in finite volumes, numerical simulations of a flow in a bent tube, of a three-dimensional laminar flame and of a low Mach number an-isotherm flow. Non conformal finite elements are also presented (Rannacher-Turek and Crouzeix-Raviart elements)

A self-adaptive finite element approach for simulation of mixed-mode delamination using cohesive zone models

NARCIS (Netherlands)

Samimi, M.; Dommelen, van J.A.W.; Geers, M.G.D.

2011-01-01

Oscillations observed in the load–displacement response of brittle interfaces modeled by cohesive zone elements in a quasi-static finite element framework are artifacts of the discretization. The typical limit points in this oscillatory path can be traced by application of path-following techniques,

Algorithms and data structures for massively parallel generic adaptive finite element codes

KAUST Repository

Bangerth, Wolfgang

2011-12-01

Today\\'s largest supercomputers have 100,000s of processor cores and offer the potential to solve partial differential equations discretized by billions of unknowns. However, the complexity of scaling to such large machines and problem sizes has so far prevented the emergence of generic software libraries that support such computations, although these would lower the threshold of entry and enable many more applications to benefit from large-scale computing. We are concerned with providing this functionality for mesh-adaptive finite element computations. We assume the existence of an "oracle" that implements the generation and modification of an adaptive mesh distributed across many processors, and that responds to queries about its structure. Based on querying the oracle, we develop scalable algorithms and data structures for generic finite element methods. Specifically, we consider the parallel distribution of mesh data, global enumeration of degrees of freedom, constraints, and postprocessing. Our algorithms remove the bottlenecks that typically limit large-scale adaptive finite element analyses. We demonstrate scalability of complete finite element workflows on up to 16,384 processors. An implementation of the proposed algorithms, based on the open source software p4est as mesh oracle, is provided under an open source license through the widely used deal.II finite element software library. © 2011 ACM 0098-3500/2011/12-ART10 $10.00.

Application of the finite element method to neutronics problems with inhomogeneous boundray conditions

International Nuclear Information System (INIS)

Yoo, K.J.

1982-01-01

The albedo boundary conditions are incorporated into the finite element method using bicubic Hermite element functions in order to reduce the computer memory and computation time in two-group diffusion calculations by excluding the reflector regions in computation space. The basis functions at the core-reflector interfaces are newly established to satisfy the albedo boundary conditions, and then the ''weak'' form of two-group diffusion equations is discretized using the principle of the weighted residual method in combination with the Galerkin approximation. The discretized two-group diffusion equation is then solved by the Gaussian elimination method with the scaled column pivoting algorithm in one-dimensional problem and Gauss-Seidel method in two-dimensional problem. Prior to the application of the method to two-group diffusion problems, the same method is applied to the one-speed neutron transport equation in a bare slab reactor with the vacuum boundary condition to confirm its usefulness in the diffusion calculations. To investigate the applicability of our diffusion method, several numerical calculations are performed: two-dimensional IAEA benchmark problem and two-dimensional ZION problem. The results are compared with the available results from the conventional finite difference and other finite element methods. If the albedo values are appropriately adjusted, our results of the two-dimensional IAEA benchmark problem are agreed within 0.002% of ksub(eff) with the fine mesh PDQ results. Comparing with CITATION results, one-eighth of core memory and one-fifteenth of computing time are required to obtain the same accuracy even though no acceleration technique is used in the present case. Also, it is found that the results are comparable with the other finite element results. However, no significant saving is obtained in computation time comparing with the other finite element results, where the reflector regions are explicity included. This mainly comes from

Lagrangian finite element method for 3D time-dependent non-isothermal flow of K-BKZ fluids

DEFF Research Database (Denmark)

Román Marín, José Manuel; Rasmussen, Henrik K.

2009-01-01

equation is replaced with a temperature dependent pseudo time. The spatial coordinate system attached to the particles is discretized by 10-node quadratic tetrahedral elements using Cartesian coordinates. The temperature and the pressure are discretized by 10-node quadratic and linear interpolation...... utilizing an implicit variable step backwards differencing (BDF2) scheme, obtaining second order convergence of the temperature in time. A quadratic interpolation in time is applied to approximate the time integral in the K-BKZ equation. This type of scheme ensures third order accuracy with respect......, respectively, in the tetrahedral particle elements. The spatial discretization of the governing equations follows a mixed Galerkin finite element method. This type of scheme ensures third order accuracy with respect to the discretization of spatial dimension. The temperature equation is solved in time...

The Galerkin Finite Element Method for A Multi-term Time-Fractional Diffusion equation

OpenAIRE

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

Numerical analysis of creep brittle rupture by the finite element method

International Nuclear Information System (INIS)

Goncalves, O.J.A.; Owen, D.R.J.

1983-01-01

In this work an implicit algorithm is proposed for the numerical analysis of creep brittle rupture problems by the finite element method. This kind of structural failure, typical in components operating at high temperatures for long periods of time, is modelled using either a three dimensional generalization of the Kachanov-Rabotnov equations due to Leckie and Hayhurst or the Monkman-Grant fracture criterion together with the Linear Life Fraction Rule. The finite element equations are derived by the displacement method and isoparametric elements are used for the spatial discretization. Geometric nonlinear effects (large displacements) are accounted for by an updated Lagrangian formulation. Attention is also focussed on the solution of the highly stiff differential equations that govern damage growth. Finally the numerical results of a three-dimensional analysis of a pressurized thin cylinder containing oxidised pits in its external wall are discussed. (orig.)

Massively Parallel Finite Element Programming

KAUST Repository

Heister, Timo; Kronbichler, Martin; Bangerth, Wolfgang

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.

Massively Parallel Finite Element Programming

KAUST Repository

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.

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

Finite element formulation for fluid-structure interaction in three-dimensional space

International Nuclear Information System (INIS)

Kulak, R.F.

1979-01-01

A development is presented for a three-dimension hexahedral hydrodynamic finite-element. Using trilinear shape functions and assuming a constant pressure field in each element, simple relations were obtained for internal nodal forces. Because the formulation was based upon a rate approach it was applicable to problems involving large displacements. This element was incorporated into an existing plate-shell finite element code. Diagonal mass matrices were used and the resulting discrete equations of motion were solved using explicit temporal integrator. Results for several problems were presented which compare numerical predictions to closed form analytical solutions. In addition, the fluid-structure interaction problem of a fluid-filled, cylindrical vessel containing internal cylinders was studied. The internal cylinders were cantilever supported from the top cover of the vessel and were periodically located circumferentially at a fixed radius. A pressurized cylindrical cavity located at the bottom of the vessel at its centerline provided the loading

Comparison between a finite difference model (PUMA) and a finite element model (DELFIN) for simulation of the reactor of the atomic power plant of Atucha I

International Nuclear Information System (INIS)

Grant, C.R.

1996-01-01

The reactor code PUMA, developed in CNEA, simulates nuclear reactors discretizing space in finite difference elements. Core representation is performed by means a cylindrical mesh, but the reactor channels are arranged in an hexagonal lattice. That is why a mapping using volume intersections must be used. This spatial treatment is the reason of an overestimation of the control rod reactivity values, which must be adjusted modifying the incremental cross sections. Also, a not very good treatment of the continuity conditions between core and reflector leads to an overestimation of channel power of the peripherical fuel elements between 5 to 8 per cent. Another code, DELFIN, developed also in CNEA, treats the spatial discretization using heterogeneous finite elements, allowing a correct treatment of the continuity of fluxes and current among elements and a more realistic representation of the hexagonal lattice of the reactor. A comparison between results obtained using both methods in done in this paper. (author). 4 refs., 3 figs

Vibrations And Deformations Of Moderately Thick Plates In Stochastic Finite Element Method

Directory of Open Access Journals (Sweden)

Grzywiński Maksym

2015-12-01

Full Text Available The paper deals with some chosen aspects of stochastic dynamical analysis of moderately thick plates. The discretization of the governing equations is described by the finite element method. The main aim of the study is to provide the generalized stochastic perturbation technique based on classical Taylor expansion with a single random variable.

L∞-error estimates of a finite element method for the Hamilton-Jacobi-Bellman equations

International Nuclear Information System (INIS)

Bouldbrachene, M.

1994-11-01

We study the finite element approximation for the solution of the Hamilton-Jacobi-Bellman equations involving a system of quasi-variational inequalities (QVI). We also give the optimal L ∞ -error estimates, using the concepts of subsolutions and discrete regularity. (author). 7 refs

Wavelet-based spectral finite element dynamic analysis for an axially moving Timoshenko beam

Science.gov (United States)

Mokhtari, Ali; Mirdamadi, Hamid Reza; Ghayour, Mostafa

2017-08-01

In this article, wavelet-based spectral finite element (WSFE) model is formulated for time domain and wave domain dynamic analysis of an axially moving Timoshenko beam subjected to axial pretension. The formulation is similar to conventional FFT-based spectral finite element (SFE) model except that Daubechies wavelet basis functions are used for temporal discretization of the governing partial differential equations into a set of ordinary differential equations. The localized nature of Daubechies wavelet basis functions helps to rule out problems of SFE model due to periodicity assumption, especially during inverse Fourier transformation and back to time domain. The high accuracy of WSFE model is then evaluated by comparing its results with those of conventional finite element and SFE results. The effects of moving beam speed and axial tensile force on vibration and wave characteristics, and static and dynamic stabilities of moving beam are investigated.

A collocation--Galerkin finite element model of cardiac action potential propagation.

Science.gov (United States)

Rogers, J M; McCulloch, A D

1994-08-01

A new computational method was developed for modeling the effects of the geometric complexity, nonuniform muscle fiber orientation, and material inhomogeneity of the ventricular wall on cardiac impulse propagation. The method was used to solve a modification to the FitzHugh-Nagumo system of equations. The geometry, local muscle fiber orientation, and material parameters of the domain were defined using linear Lagrange or cubic Hermite finite element interpolation. Spatial variations of time-dependent excitation and recovery variables were approximated using cubic Hermite finite element interpolation, and the governing finite element equations were assembled using the collocation method. To overcome the deficiencies of conventional collocation methods on irregular domains, Galerkin equations for the no-flux boundary conditions were used instead of collocation equations for the boundary degrees-of-freedom. The resulting system was evolved using an adaptive Runge-Kutta method. Converged two-dimensional simulations of normal propagation showed that this method requires less CPU time than a traditional finite difference discretization. The model also reproduced several other physiologic phenomena known to be important in arrhythmogenesis including: Wenckebach periodicity, slowed propagation and unidirectional block due to wavefront curvature, reentry around a fixed obstacle, and spiral wave reentry. In a new result, we observed wavespeed variations and block due to nonuniform muscle fiber orientation. The findings suggest that the finite element method is suitable for studying normal and pathological cardiac activation and has significant advantages over existing techniques.

Topology optimization for three-dimensional electromagnetic waves using an edge element-based finite-element method.

Science.gov (United States)

Deng, Yongbo; Korvink, Jan G

2016-05-01

This paper develops a topology optimization procedure for three-dimensional electromagnetic waves with an edge element-based finite-element method. In contrast to the two-dimensional case, three-dimensional electromagnetic waves must include an additional divergence-free condition for the field variables. The edge element-based finite-element method is used to both discretize the wave equations and enforce the divergence-free condition. For wave propagation described in terms of the magnetic field in the widely used class of non-magnetic materials, the divergence-free condition is imposed on the magnetic field. This naturally leads to a nodal topology optimization method. When wave propagation is described using the electric field, the divergence-free condition must be imposed on the electric displacement. In this case, the material in the design domain is assumed to be piecewise homogeneous to impose the divergence-free condition on the electric field. This results in an element-wise topology optimization algorithm. The topology optimization problems are regularized using a Helmholtz filter and a threshold projection method and are analysed using a continuous adjoint method. In order to ensure the applicability of the filter in the element-wise topology optimization version, a regularization method is presented to project the nodal into an element-wise physical density variable.

Finite element computational fluid mechanics

International Nuclear Information System (INIS)

Baker, A.J.

1983-01-01

This book analyzes finite element theory as applied to computational fluid mechanics. It includes a chapter on using the heat conduction equation to expose the essence of finite element theory, including higher-order accuracy and convergence in a common knowledge framework. Another chapter generalizes the algorithm to extend application to the nonlinearity of the Navier-Stokes equations. Other chapters are concerned with the analysis of a specific fluids mechanics problem class, including theory and applications. Some of the topics covered include finite element theory for linear mechanics; potential flow; weighted residuals/galerkin finite element theory; inviscid and convection dominated flows; boundary layers; parabolic three-dimensional flows; and viscous and rotational flows

Discrete finite nilpotent Lie analogs: New models for unified gauge field theory

International Nuclear Information System (INIS)

Kornacker, K.

1978-01-01

To each finite dimensional real Lie algebra with integer structure constants there corresponds a countable family of discrete finite nilpotent Lie analogs. Each finite Lie analog maps exponentially onto a finite unipotent group G, and is isomorphic to the Lie algebra of G. Reformulation of quantum field theory in discrete finite form, utilizing nilpotent Lie analogs, should elminate all divergence problems even though some non-Abelian gauge symmetry may not be spontaneously broken. Preliminary results in the new finite representation theory indicate that a natural hierarchy of spontaneously broken symmetries can arise from a single unbroken non-Abelian gauge symmetry, and suggest the possibility of a new unified group theoretic interpretation for hadron colors and flavors

Assessment of Structural Behavior of Non-corroded and Corroded RCC Beams Using Finite Element Method

Directory of Open Access Journals (Sweden)

Anand Parande

2008-09-01

Full Text Available A three dimensional finite element model is developed to examine the structural behaviour of corroded reinforced concrete beam and non corroded reinforced concrete beam. Non linear finite element analysis is performed using the ANSYS program. SOLID 65, LINK 8 element represent concrete and discrete reinforcing steel bars, based on each component actual characteristics, non linear material properties are defined for both elements. The effect of corrosion in reinforced concrete is studied by finite element analysis; an approach is developed to model the corrosion product expansion causing concrete cover cracking for this, beam has been modeled using ANSYS and using this data the beam has been casted with M20 concrete after 28 days the beam will be tested for flexural strength. The comparison between ANSYS prediction and field data are made in terms of deflection, stress, strain, bond strength and crack pattern of concrete beam.

A finite strain Eulerian formulation for compressible and nearly incompressible hyperelasticity using high-order B-spline finite elements

KAUST Repository

Duddu, Ravindra

2011-10-05

We present a numerical formulation aimed at modeling the nonlinear response of elastic materials using large deformation continuum mechanics in three dimensions. This finite element formulation is based on the Eulerian description of motion and the transport of the deformation gradient. When modeling a nearly incompressible solid, the transport of the deformation gradient is decomposed into its isochoric part and the Jacobian determinant as independent fields. A homogeneous isotropic hyperelastic solid is assumed and B-splines-based finite elements are used for the spatial discretization. A variational multiscale residual-based approach is employed to stabilize the transport equations. The performance of the scheme is explored for both compressible and nearly incompressible applications. The numerical results are in good agreement with theory illustrating the viability of the computational scheme. © 2011 John Wiley & Sons, Ltd.

Discontinuous Galerkin finite element method with anisotropic local grid refinement for inviscid compressible flows

NARCIS (Netherlands)

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

Modelling bucket excavation by finite element

Science.gov (United States)

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

Spectral/hp least-squares finite element formulation for the Navier-Stokes equations

International Nuclear Information System (INIS)

Pontaza, J.P.; Reddy, J.N.

2003-01-01

We consider the application of least-squares finite element models combined with spectral/hp methods for the numerical solution of viscous flow problems. The paper presents the formulation, validation, and application of a spectral/hp algorithm to the numerical solution of the Navier-Stokes equations governing two- and three-dimensional stationary incompressible and low-speed compressible flows. The Navier-Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity or velocity gradients as additional independent variables and the least-squares method is used to develop the finite element model. High-order element expansions are used to construct the discrete model. The discrete model thus obtained is linearized by Newton's method, resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method. Spectral convergence of the L 2 least-squares functional and L 2 error norms is verified using smooth solutions to the two-dimensional stationary Poisson and incompressible Navier-Stokes equations. Numerical results for flow over a backward-facing step, steady flow past a circular cylinder, three-dimensional lid-driven cavity flow, and compressible buoyant flow inside a square enclosure are presented to demonstrate the predictive capability and robustness of the proposed formulation

A Lagrangian finite element method for the simulation of flow of non-newtonian liquids

DEFF Research Database (Denmark)

Hassager, Ole; Bisgaard, C

1983-01-01

A Lagrangian method for the simulation of flow of non-Newtonian liquids is implemented. The fluid mechanical equations are formulated in the form of a variational principle, and a discretization is performed by finite elements. The method is applied to the slow of a contravariant convected Maxwell...

SIMULATION OF THE BEHAVIOR OF THE WATER TABLE IN A COASTAL AQUIFER SYSTEM FINITE ELEMENT

Directory of Open Access Journals (Sweden)

Luis Lara Romero

2016-06-01

Full Text Available This paper presents the application of Galerkin method to discretize the model equation of groundwater ow in a conned aquifer semipermeable with tidal boundary conditions on one of its borders, the other borders remain constant. For the simulations was generated a numerical program, Ground Water Finite Element Method, which implements the method of nite elements with triangular elements with three nodes and a degree of freedom per node.

Additive Schwarz preconditioner for the finite volume element discretization of symmetric elliptic problems

Czech Academy of Sciences Publication Activity Database

Marcinkowski, L.; Rahman, T.; Loneland, A.; Valdman, Jan

2016-01-01

Roč. 56, č. 3 (2016), s. 967-993 ISSN 0006-3835 R&D Projects: GA ČR GA13-18652S Institutional support: RVO:67985556 Keywords : Domain decomposition * Additive Schwarz method * Finite volume element * GMRES Subject RIV: BA - General Mathematics Impact factor: 1.670, year: 2016 http://library.utia.cas.cz/separaty/2015/MTR/valdman-0447835.pdf

The finite element method and applications in engineering using ANSYS

CERN Document Server

Madenci, Erdogan

2015-01-01

This textbook offers theoretical and practical knowledge of the finite element method. The book equips readers with the skills required to analyze engineering problems using ANSYS®, a commercially available FEA program. Revised and updated, this new edition presents the most current ANSYS® commands and ANSYS® screen shots, as well as modeling steps for each example problem. This self-contained, introductory text minimizes the need for additional reference material by covering both the fundamental topics in finite element methods and advanced topics concerning modeling and analysis. It focuses on the use of ANSYS® through both the Graphics User Interface (GUI) and the ANSYS® Parametric Design Language (APDL). Extensive examples from a range of engineering disciplines are presented in a straightforward, step-by-step fashion. Key topics include: • An introduction to FEM • Fundamentals and analysis capabilities of ANSYS® • Fundamentals of discretization and approximation functions • Modeling techniq...

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

Implicit three-dimensional finite-element formulation for the nonlinear structural response of reactor components

International Nuclear Information System (INIS)

Kulak, R.F.; Belytschko, T.B.

1975-09-01

The formulation of a finite-element procedure for the implicit transient and static analysis of plate/shell type structures in three-dimensional space is described. The triangular plate/shell element can sustain both membrane and bending stresses. Both geometric and material nonlinearities can be treated, and an elastic-plastic material law has been incorporated. The formulation permits the element to undergo arbitrarily large rotations and translations; but, in its present form it is restricted to small strains. The discretized equations of motion are obtained by a stiffness method. An implicit integration algorithm based on trapezoidal integration formulas is used to integrate the discretized equations of motion in time. To ensure numerical stability, an iterative solution procedure with equilibrium checks is used

A first course in finite elements

CERN Document Server

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

On the time-stepping stability of continuous mass-lumped and discontinuous Galerkin finite elements for the 3D acoustic wave equation

NARCIS (Netherlands)

Zhebel, E.; Minisini, S.; Mulder, W.A.

2012-01-01

We solve the three-dimensional acoustic wave equation, discretized on tetrahedral meshes. Two methods are considered: mass-lumped continuous finite elements and the symmetric interior-penalty discontinuous Galerkin method (SIP-DG). Combining the spatial discretization with the leap-frog

A mixed finite element method for nonlinear diffusion equations

KAUST Repository

Burger, Martin; Carrillo, José ; Wolfram, Marie-Therese

2010-01-01

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model. © American Institute of Mathematical Sciences.

On Stochastic Finite-Time Control of Discrete-Time Fuzzy Systems with Packet Dropout

Directory of Open Access Journals (Sweden)

Yingqi Zhang

2012-01-01

Full Text Available This paper is concerned with the stochastic finite-time stability and stochastic finite-time boundedness problems for one family of fuzzy discrete-time systems over networks with packet dropout, parametric uncertainties, and time-varying norm-bounded disturbance. Firstly, we present the dynamic model description studied, in which the discrete-time fuzzy T-S systems with packet loss can be described by one class of fuzzy Markovian jump systems. Then, the concepts of stochastic finite-time stability and stochastic finite-time boundedness and problem formulation are given. Based on Lyapunov function approach, sufficient conditions on stochastic finite-time stability and stochastic finite-time boundedness are established for the resulting closed-loop fuzzy discrete-time system with Markovian jumps, and state-feedback controllers are designed to ensure stochastic finite-time stability and stochastic finite-time boundedness of the class of fuzzy systems. The stochastic finite-time stability and stochastic finite-time boundedness criteria can be tackled in the form of linear matrix inequalities with a fixed parameter. As an auxiliary result, we also give sufficient conditions on the stochastic stability of the class of fuzzy T-S systems with packet loss. Finally, two illustrative examples are presented to show the validity of the developed methodology.

Finite Element Simulation of Sheet Metal Forming Process Using Local Interpolation for Tool Surfaces

International Nuclear Information System (INIS)

Hama, Takayuki; Takuda, Hirohiko; Takamura, Masato; Makinouchi, Akitake; Teodosiu, Cristian

2005-01-01

Treatment of contact between a sheet and tools is one of the most difficult problems to deal with in finite-element simulations of sheet forming processes. In order to obtain more accurate tool models without increasing the number of elements, this paper describes a new formulation for contact problems using interpolation proposed by Nagata for tool surfaces. A contact search algorithm between sheet nodes and the interpolated tool surfaces was developed and was introduced into the static-explicit elastoplastic finite-element method code STAMP3D. Simulations of a square cup deep drawing process with a very coarsely discretized punch model were carried out. The simulated results showed that the proposed algorithm gave the proper drawn shape, demonstrating the validity of the proposed algorithm

Nonlinear magnetohydrodynamics simulation using high-order finite elements

International Nuclear Information System (INIS)

Plimpton, Steven James; Schnack, D.D.; Tarditi, A.; Chu, M.S.; Gianakon, T.A.; Kruger, S.E.; Nebel, R.A.; Barnes, D.C.; Sovinec, C.R.; Glasser, A.H.

2005-01-01

A conforming representation composed of 2D finite elements and finite Fourier series is applied to 3D nonlinear non-ideal magnetohydrodynamics using a semi-implicit time-advance. The self-adjoint semi-implicit operator and variational approach to spatial discretization are synergistic and enable simulation in the extremely stiff conditions found in high temperature plasmas without sacrificing the geometric flexibility needed for modeling laboratory experiments. Growth rates for resistive tearing modes with experimentally relevant Lundquist number are computed accurately with time-steps that are large with respect to the global Alfven time and moderate spatial resolution when the finite elements have basis functions of polynomial degree (p) two or larger. An error diffusion method controls the generation of magnetic divergence error. Convergence studies show that this approach is effective for continuous basis functions with p (ge) 2, where the number of test functions for the divergence control terms is less than the number of degrees of freedom in the expansion for vector fields. Anisotropic thermal conduction at realistic ratios of parallel to perpendicular conductivity (x(parallel)/x(perpendicular)) is computed accurately with p (ge) 3 without mesh alignment. A simulation of tearing-mode evolution for a shaped toroidal tokamak equilibrium demonstrates the effectiveness of the algorithm in nonlinear conditions, and its results are used to verify the accuracy of the numerical anisotropic thermal conduction in 3D magnetic topologies.

Modelling optimization involving different types of elements in finite element analysis

International Nuclear Information System (INIS)

Wai, C M; Rivai, Ahmad; Bapokutty, Omar

2013-01-01

Finite elements are used to express the mechanical behaviour of a structure in finite element analysis. Therefore, the selection of the elements determines the quality of the analysis. The aim of this paper is to compare and contrast 1D element, 2D element, and 3D element used in finite element analysis. A simple case study was carried out on a standard W460x74 I-beam. The I-beam was modelled and analyzed statically with 1D elements, 2D elements and 3D elements. The results for the three separate finite element models were compared in terms of stresses, deformation and displacement of the I-beam. All three finite element models yield satisfactory results with acceptable errors. The advantages and limitations of these elements are discussed. 1D elements offer simplicity although lacking in their ability to model complicated geometry. 2D elements and 3D elements provide more detail yet sophisticated results which require more time and computer memory in the modelling process. It is also found that the choice of element in finite element analysis is influence by a few factors such as the geometry of the structure, desired analysis results, and the capability of the computer

A strongly conservative finite element method for the coupling of Stokes and Darcy flow

KAUST Repository

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.

Adaptive finite element method for shape optimization

KAUST Repository

Morin, Pedro; Nochetto, Ricardo H.; Pauletti, Miguel S.; Verani, Marco

2012-01-01

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.

Adaptive finite element method for shape optimization

KAUST Repository

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.

Pellet Cladding Mechanical Interaction Modeling Using the Extended Finite Element Method

Energy Technology Data Exchange (ETDEWEB)

Spencer, Benjamin W.; Jiang, Wen; Dolbow, John E.; Peco, Christian

2016-09-01

As a brittle material, the ceramic UO2 used as light water reactor fuel experiences significant fracturing throughout its life, beginning with the first rise to power of fresh fuel. This has multiple effects on the thermal and mechanical response of the fuel/cladding system. One such effect that is particularly important is that when there is mechanical contact between the fuel and cladding, cracks that extending from the outer surface of the fuel into the volume of the fuel cause elevated stresses in the adjacent cladding, which can potentially lead to cladding failure. Modeling the thermal and mechanical response of the cladding in the vicinity of these surface-breaking cracks in the fuel can provide important insights into this behavior to help avoid operating conditions that could lead to cladding failure. Such modeling has traditionally been done in the context of finite-element-based fuel performance analysis by modifying the fuel mesh to introduce discrete cracks. While this approach is effective in capturing the important behavior at the fuel/cladding interface, there are multiple drawbacks to explicitly incorporating the cracks in the finite element mesh. Because the cracks are incorporated in the original mesh, the mesh must be modified for cracks of specified location and depth, so it is difficult to account for crack propagation and the formation of new cracks at other locations. The extended finite element method (XFEM) has emerged in recent years as a powerful method to represent arbitrary, evolving, discrete discontinuities within the context of the finite element method. Development work is underway by the authors to implement XFEM in the BISON fuel performance code, and this capability has previously been demonstrated in simulations of fracture propagation in ceramic nuclear fuel. These preliminary demonstrations have included only the fuel, and excluded the cladding for simplicity. This paper presents initial results of efforts to apply XFEM to

Adaptive Multilevel Methods with Local Smoothing for $H^1$- and $H^{\\mathrm{curl}}$-Conforming High Order Finite Element Methods

KAUST Repository

Janssen, Bä rbel; Kanschat, Guido

2011-01-01

A multilevel method on adaptive meshes with hanging nodes is presented, and the additional matrices appearing in the implementation are derived. Smoothers of overlapping Schwarz type are discussed; smoothing is restricted to the interior of the subdomains refined to the current level; thus it has optimal computational complexity. When applied to conforming finite element discretizations of elliptic problems and Maxwell equations, the method's convergence rates are very close to those for the nonadaptive version. Furthermore, the smoothers remain efficient for high order finite elements. We discuss the implementation in a general finite element code using the example of the deal.II library. © 2011 Societ y for Industrial and Applied Mathematics.

Mechanical strength calculation of the disk type windings with elastic couplings by the finite element method

International Nuclear Information System (INIS)

Sivkova, G.N.; Spirchenko, Yu.V.; Chvartatskij, P.V.

1981-01-01

Stressed-deformed state of toroidal field coils of the disc type with elastic couplings of the tokamaks has been investigated with provision for the effect of the central core pliability by means of the two-dimensional version of the finite element method. Numerical solution of the finite element method is performed by means of the ES 1040 computer according to the computer code permitting taking account of boundary conditions of elastic support. The calculation has been performed using as the example the project of T-20 facility coil of the disc type. Consideration of pliability of the central core of the facility inductor is accomplished by the introduction of additional rigidities to the complete matrix of rigidity. Scheme of the structure distretization includes 141 units, 211 elements. The accuracy of solution depends on the reduction accuracy of the volume load to unit forces and on the number of finite elements. Analysis of the solution convergence is performed by the comparison of solutions obtained for three different schemes of the disk discretization without regard for the inductor pliability. The comparative analysis of the results shows that transfer epures for all the three discretization versions practically coincide and stresses differ not more than by 10%. On the whole the above investigation has demonstrated good convergence of the problem solution [ru

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

Finite discrete field theory

International Nuclear Information System (INIS)

Souza, Manoelito M. de

1997-01-01

We discuss the physical meaning and the geometric interpretation of implementation in classical field theories. The origin of infinities and other inconsistencies in field theories is traced to fields defined with support on the light cone; a finite and consistent field theory requires a light-cone generator as the field support. Then, we introduce a classical field theory with support on the light cone generators. It results on a description of discrete (point-like) interactions in terms of localized particle-like fields. We find the propagators of these particle-like fields and discuss their physical meaning, properties and consequences. They are conformally invariant, singularity-free, and describing a manifestly covariant (1 + 1)-dimensional dynamics in a (3 = 1) spacetime. Remarkably this conformal symmetry remains even for the propagation of a massive field in four spacetime dimensions. We apply this formalism to Classical electrodynamics and to the General Relativity Theory. The standard formalism with its distributed fields is retrieved in terms of spacetime average of the discrete field. Singularities are the by-products of the averaging process. This new formalism enlighten the meaning and the problem of field theory, and may allow a softer transition to a quantum theory. (author)

Structural modeling techniques by finite element method

International Nuclear Information System (INIS)

Kang, Yeong Jin; Kim, Geung Hwan; Ju, Gwan Jeong

1991-01-01

This book includes introduction table of contents chapter 1 finite element idealization introduction summary of the finite element method equilibrium and compatibility in the finite element solution degrees of freedom symmetry and anti symmetry modeling guidelines local analysis example references chapter 2 static analysis structural geometry finite element models analysis procedure modeling guidelines references chapter 3 dynamic analysis models for dynamic analysis dynamic analysis procedures modeling guidelines and modeling guidelines.

Discrete elements method of neutral particle transport

International Nuclear Information System (INIS)

Mathews, K.A.

1983-01-01

A new discrete elements (L/sub N/) transport method is derived and compared to the discrete ordinates S/sub N/ method, theoretically and by numerical experimentation. The discrete elements method is more accurate than discrete ordinates and strongly ameliorates ray effects for the practical problems studied. The discrete elements method is shown to be more cost effective, in terms of execution time with comparable storage to attain the same accuracy, for a one-dimensional test case using linear characteristic spatial quadrature. In a two-dimensional test case, a vacuum duct in a shield, L/sub N/ is more consistently convergent toward a Monte Carlo benchmark solution than S/sub N/, using step characteristic spatial quadrature. An analysis of the interaction of angular and spatial quadrature in xy-geometry indicates the desirability of using linear characteristic spatial quadrature with the L/sub N/ method

Compatible-strain mixed finite element methods for incompressible nonlinear elasticity

Science.gov (United States)

Faghih Shojaei, Mostafa; Yavari, Arash

2018-05-01

We introduce a new family of mixed finite elements for incompressible nonlinear elasticity - compatible-strain mixed finite element methods (CSFEMs). Based on a Hu-Washizu-type functional, we write a four-field mixed formulation with the displacement, the displacement gradient, the first Piola-Kirchhoff stress, and a pressure-like field as the four independent unknowns. Using the Hilbert complexes of nonlinear elasticity, which describe the kinematics and the kinetics of motion, we identify the solution spaces of the independent unknown fields. In particular, we define the displacement in H1, the displacement gradient in H (curl), the stress in H (div), and the pressure field in L2. The test spaces of the mixed formulations are chosen to be the same as the corresponding solution spaces. Next, in a conforming setting, we approximate the solution and the test spaces with some piecewise polynomial subspaces of them. Among these approximation spaces are the tensorial analogues of the Nédélec and Raviart-Thomas finite element spaces of vector fields. This approach results in compatible-strain mixed finite element methods that satisfy both the Hadamard compatibility condition and the continuity of traction at the discrete level independently of the refinement level of the mesh. By considering several numerical examples, we demonstrate that CSFEMs have a good performance for bending problems and for bodies with complex geometries. CSFEMs are capable of capturing very large strains and accurately approximating stress and pressure fields. Using CSFEMs, we do not observe any numerical artifacts, e.g., checkerboarding of pressure, hourglass instability, or locking in our numerical examples. Moreover, CSFEMs provide an efficient framework for modeling heterogeneous solids.

Solution of finite element problems using hybrid parallelization with MPI and OpenMP Solution of finite element problems using hybrid parallelization with MPI and OpenMP

Directory of Open Access Journals (Sweden)

José Miguel Vargas-Félix

2012-11-01

Full Text Available The Finite Element Method (FEM is used to solve problems like solid deformation and heat diffusion in domains with complex geometries. This kind of geometries requires discretization with millions of elements; this is equivalent to solve systems of equations with sparse matrices and tens or hundreds of millions of variables. The aim is to use computer clusters to solve these systems. The solution method used is Schur substructuration. Using it is possible to divide a large system of equations into many small ones to solve them more efficiently. This method allows parallelization. MPI (Message Passing Interface is used to distribute the systems of equations to solve each one in a computer of a cluster. Each system of equations is solved using a solver implemented to use OpenMP as a local parallelization method.The Finite Element Method (FEM is used to solve problems like solid deformation and heat diffusion in domains with complex geometries. This kind of geometries requires discretization with millions of elements; this is equivalent to solve systems of equations with sparse matrices and tens or hundreds of millions of variables. The aim is to use computer clusters to solve these systems. The solution method used is Schur substructuration. Using it is possible to divide a large system of equations into many small ones to solve them more efficiently. This method allows parallelization. MPI (Message Passing Interface is used to distribute the systems of equations to solve each one in a computer of a cluster. Each system of equations is solved using a solver implemented to use OpenMP as a local parallelization method.

Element-topology-independent preconditioners for parallel finite element computations

Science.gov (United States)

Park, K. C.; Alexander, Scott

1992-01-01

A family of preconditioners for the solution of finite element equations are presented, which are element-topology independent and thus can be applicable to element order-free parallel computations. A key feature of the present preconditioners is the repeated use of element connectivity matrices and their left and right inverses. The properties and performance of the present preconditioners are demonstrated via beam and two-dimensional finite element matrices for implicit time integration computations.

A quasi-Lagrangian finite element method for the Navier-Stokes equations in a time-dependent domain

Science.gov (United States)

Lozovskiy, Alexander; Olshanskii, Maxim A.; Vassilevski, Yuri V.

2018-05-01

The paper develops a finite element method for the Navier-Stokes equations of incompressible viscous fluid in a time-dependent domain. The method builds on a quasi-Lagrangian formulation of the problem. The paper provides stability and convergence analysis of the fully discrete (finite-difference in time and finite-element in space) method. The analysis does not assume any CFL time-step restriction, it rather needs mild conditions of the form $\\Delta t\\le C$, where $C$ depends only on problem data, and $h^{2m_u+2}\\le c\\,\\Delta t$, $m_u$ is polynomial degree of velocity finite element space. Both conditions result from a numerical treatment of practically important non-homogeneous boundary conditions. The theoretically predicted convergence rate is confirmed by a set of numerical experiments. Further we apply the method to simulate a flow in a simplified model of the left ventricle of a human heart, where the ventricle wall dynamics is reconstructed from a sequence of contrast enhanced Computed Tomography images.

Finite element modeling of tornado missile impact on reinforced concrete wall panels

International Nuclear Information System (INIS)

Zhang, Y.; Vallabhan, C.V.G.; McDonald, J.R.

1993-01-01

This paper describes a finite element model for the impact of large tornado-generated missiles with reinforced concrete wall panels. The analysis predicts the dynamic response of a wall panel when impacted by a missile with a large contact area such as an automobile. Quadratic finite elements are used to discretize the domain of the wall panel. Fundamental assumptions are based on the Mindlin and the related Reinsser plate theories. An 'embedded' model is employed to account for the reinforcing bars. The nonlinear behavior of concrete and steel bars are analyzed by means of time-dependent constitutive relationships. A model is proposed to describe the initial and subsequent yield surfaces of concrete material, which avoids underestimation of the effect of high hydrostatic stresses on the yielding behavior of concrete. Ottosen's four-parameter failure criterion is used to define the failure surface of concrete. A crack monitoring algorithm accounts for post-cracking and post-crushing behavior of concrete. Explicit time step integration of nonlinear dynamic equations are carried out using the finite element discretization of a concrete wall panel. As a practical application of the analysis technique, the contact failure pressure for a particular panel geometry can be calculated. The contact failure pressure and the elapsed time to failure after missile contact define a rectangular or triangular impulse loading to produce failure of the panel. Since automobile crashes are known to produce triangular impulse loads, the two pulses (failure and impact) can be compared to determine if a particular impact will fail the panels. Thus, a particular concrete panel can be analyzed to determine if it will fail under a postulated missile impact

Theoretical formulation of finite-dimensional discrete phase spaces: I. Algebraic structures and uncertainty principles

International Nuclear Information System (INIS)

Marchiolli, M.A.; Ruzzi, M.

2012-01-01

We propose a self-consistent theoretical framework for a wide class of physical systems characterized by a finite space of states which allows us, within several mathematical virtues, to construct a discrete version of the Weyl–Wigner–Moyal (WWM) formalism for finite-dimensional discrete phase spaces with toroidal topology. As a first and important application from this ab initio approach, we initially investigate the Robertson–Schrödinger (RS) uncertainty principle related to the discrete coordinate and momentum operators, as well as its implications for physical systems with periodic boundary conditions. The second interesting application is associated with a particular uncertainty principle inherent to the unitary operators, which is based on the Wiener–Khinchin theorem for signal processing. Furthermore, we also establish a modified discrete version for the well-known Heisenberg–Kennard–Robertson (HKR) uncertainty principle, which exhibits additional terms (or corrections) that resemble the generalized uncertainty principle (GUP) into the context of quantum gravity. The results obtained from this new algebraic approach touch on some fundamental questions inherent to quantum mechanics and certainly represent an object of future investigations in physics. - Highlights: ► We construct a discrete version of the Weyl–Wigner–Moyal formalism. ► Coherent states for finite-dimensional discrete phase spaces are established. ► Discrete coordinate and momentum operators are properly defined. ► Uncertainty principles depend on the topology of finite physical systems. ► Corrections for the discrete Heisenberg uncertainty relation are also obtained.

Hierarchical multiscale modeling for flows in fractured media using generalized multiscale finite element method

KAUST Repository

Efendiev, Yalchin R.

2015-06-05

In this paper, we develop a multiscale finite element method for solving flows in fractured media. Our approach is based on generalized multiscale finite element method (GMsFEM), where we represent the fracture effects on a coarse grid via multiscale basis functions. These multiscale basis functions are constructed in the offline stage via local spectral problems following GMsFEM. To represent the fractures on the fine grid, we consider two approaches (1) discrete fracture model (DFM) (2) embedded fracture model (EFM) and their combination. In DFM, the fractures are resolved via the fine grid, while in EFM the fracture and the fine grid block interaction is represented as a source term. In the proposed multiscale method, additional multiscale basis functions are used to represent the long fractures, while short-size fractures are collectively represented by a single basis functions. The procedure is automatically done via local spectral problems. In this regard, our approach shares common concepts with several approaches proposed in the literature as we discuss. We would like to emphasize that our goal is not to compare DFM with EFM, but rather to develop GMsFEM framework which uses these (DFM or EFM) fine-grid discretization techniques. Numerical results are presented, where we demonstrate how one can adaptively add basis functions in the regions of interest based on error indicators. We also discuss the use of randomized snapshots (Calo et al. Randomized oversampling for generalized multiscale finite element methods, 2014), which reduces the offline computational cost.

Probabilistic finite elements

Science.gov (United States)

Belytschko, Ted; Wing, Kam Liu

1987-01-01

In the Probabilistic Finite Element Method (PFEM), finite element methods have been efficiently combined with second-order perturbation techniques to provide an effective method for informing the designer of the range of response which is likely in a given problem. The designer must provide as input the statistical character of the input variables, such as yield strength, load magnitude, and Young's modulus, by specifying their mean values and their variances. The output then consists of the mean response and the variance in the response. Thus the designer is given a much broader picture of the predicted performance than with simply a single response curve. These methods are applicable to a wide class of problems, provided that the scale of randomness is not too large and the probabilistic density functions possess decaying tails. By incorporating the computational techniques we have developed in the past 3 years for efficiency, the probabilistic finite element methods are capable of handling large systems with many sources of uncertainties. Sample results for an elastic-plastic ten-bar structure and an elastic-plastic plane continuum with a circular hole subject to cyclic loadings with the yield stress on the random field are given.

Finite Element in Angle Unit Sphere Meshing for Charged Particle Transport.

Energy Technology Data Exchange (ETDEWEB)

Ortega, Mario Ivan [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States); Drumm, Clifton R. [Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

2017-10-01

Finite element in angle formulations of the charged particle transport equation require the discretization of the unit sphere. In Sceptre, a three-dimensional surface mesh of a sphere is transformed into a two-dimensional mesh. Projection of a sphere onto a two-dimensional surface is well studied with map makers spending the last few centuries attempting to create maps that preserve proportion and area. Using these techniques, various meshing schemes for the unit sphere were investigated.

Finite element and finite difference methods in electromagnetic scattering

CERN Document Server

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

A finite element method for solving the shallow water equations on the sphere

Science.gov (United States)

Comblen, Richard; Legrand, Sébastien; Deleersnijder, Eric; Legat, Vincent

Within the framework of ocean general circulation modeling, the present paper describes an efficient way to discretize partial differential equations on curved surfaces by means of the finite element method on triangular meshes. Our approach benefits from the inherent flexibility of the finite element method. The key idea consists in a dialog between a local coordinate system defined for each element in which integration takes place, and a nodal coordinate system in which all local contributions related to a vectorial degree of freedom are assembled. Since each element of the mesh and each degree of freedom are treated in the same way, the so-called pole singularity issue is fully circumvented. Applied to the shallow water equations expressed in primitive variables, this new approach has been validated against the standard test set defined by [Williamson, D.L., Drake, J.B., Hack, J.J., Jakob, R., Swarztrauber, P.N., 1992. A standard test set for numerical approximations to the shallow water equations in spherical geometry. Journal of Computational Physics 102, 211-224]. Optimal rates of convergence for the P1NC-P1 finite element pair are obtained, for both global and local quantities of interest. Finally, the approach can be extended to three-dimensional thin-layer flows in a straightforward manner.

Electrical machine analysis using finite elements

CERN Document Server

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

Adaptive Multilevel Methods with Local Smoothing for $H^1$- and $H^{\\mathrm{curl}}$-Conforming High Order Finite Element Methods

KAUST Repository

Janssen, Bärbel

2011-01-01

A multilevel method on adaptive meshes with hanging nodes is presented, and the additional matrices appearing in the implementation are derived. Smoothers of overlapping Schwarz type are discussed; smoothing is restricted to the interior of the subdomains refined to the current level; thus it has optimal computational complexity. When applied to conforming finite element discretizations of elliptic problems and Maxwell equations, the method\\'s convergence rates are very close to those for the nonadaptive version. Furthermore, the smoothers remain efficient for high order finite elements. We discuss the implementation in a general finite element code using the example of the deal.II library. © 2011 Societ y for Industrial and Applied Mathematics.

Stochastic ℋ∞ Finite-Time Control of Discrete-Time Systems with Packet Loss

Directory of Open Access Journals (Sweden)

Yingqi Zhang

2012-01-01

Full Text Available This paper investigates the stochastic finite-time stabilization and ℋ∞ control problem for one family of linear discrete-time systems over networks with packet loss, parametric uncertainties, and time-varying norm-bounded disturbance. Firstly, the dynamic model description studied is given, which, if the packet dropout is assumed to be a discrete-time homogenous Markov process, the class of discrete-time linear systems with packet loss can be regarded as Markovian jump systems. Based on Lyapunov function approach, sufficient conditions are established for the resulting closed-loop discrete-time system with Markovian jumps to be stochastic ℋ∞ finite-time boundedness and then state feedback controllers are designed to guarantee stochastic ℋ∞ finite-time stabilization of the class of stochastic systems. The stochastic ℋ∞ finite-time boundedness criteria can be tackled in the form of linear matrix inequalities with a fixed parameter. As an auxiliary result, we also give sufficient conditions on the robust stochastic stabilization of the class of linear systems with packet loss. Finally, simulation examples are presented to illustrate the validity of the developed scheme.

Development of polygon elements based on the scaled boundary finite element method

International Nuclear Information System (INIS)

Chiong, Irene; Song Chongmin

2010-01-01

We aim to extend the scaled boundary finite element method to construct conforming polygon elements. The development of the polygonal finite element is highly anticipated in computational mechanics as greater flexibility and accuracy can be achieved using these elements. The scaled boundary polygonal finite element will enable new developments in mesh generation, better accuracy from a higher order approximation and better transition elements in finite element meshes. Polygon elements of arbitrary number of edges and order have been developed successfully. The edges of an element are discretised with line elements. The displacement solution of the scaled boundary finite element method is used in the development of shape functions. They are shown to be smooth and continuous within the element, and satisfy compatibility and completeness requirements. Furthermore, eigenvalue decomposition has been used to depict element modes and outcomes indicate the ability of the scaled boundary polygonal element to express rigid body and constant strain modes. Numerical tests are presented; the patch test is passed and constant strain modes verified. Accuracy and convergence of the method are also presented and the performance of the scaled boundary polygonal finite element is verified on Cook's swept panel problem. Results show that the scaled boundary polygonal finite element method outperforms a traditional mesh and accuracy and convergence are achieved from fewer nodes. The proposed method is also shown to be truly flexible, and applies to arbitrary n-gons formed of irregular and non-convex polygons.

New non-structured discretizations for fluid flows with reinforced incompressibility

International Nuclear Information System (INIS)

Heib, S.

2003-01-01

This work deals with the discretization of Stokes and Navier-Stokes equations modeling the flow of incompressible fluids on 2-D or 3-D non-structured meshes. Triangles and tetrahedrons are used for 2-D and 3-D meshes, respectively. The developments and calculations are performed with the code Priceles (fast CEA-EdF industrial platform for large Eddy simulation). This code allows to perform simulations both on structured and non-structured meshes. A finite-volume resolution method is used: a finite difference volume (FDV) method is used for the structured meshes and a finite element volume (FEV) method is used for the non-structured meshes. The finite element used in the beginning of this work has several defects. Starting from this situation, the discretization is improved by adding modifications to this element and the new elements introduced are analyzed theoretically. In parallel to these analyses, the new discretizations are implemented in order to test them numerically and to confirm the theoretical analyses. The first chapter presents the physical and mathematical modeling used in this work. The second chapter treats of the discretization of Stokes equations and presents the FEV resolution method. Chapter 3 presents a first attempt of improvement of this finite element and leads to the proposal of a new element which is presented in details. The problem encountered with the new discretization leads to a modification presented in chapter 4. This new discretization gives all the expected convergence results and sometimes shows super-convergence properties. Chapter 5 deals with the study and discretization of the Navier-Stokes equations. The study of the filtered Navier-Stokes equations, used for large Eddy simulations, requires to give a particular attention to the discretization of the diffusive terms. Then, the convective terms are considered. The effects of the convective terms in the initial discretization and in the improved method are compared. The use of

An angularly refineable phase space finite element method with approximate sweeping procedure

International Nuclear Information System (INIS)

Kophazi, J.; Lathouwers, D.

2013-01-01

An angularly refineable phase space finite element method is proposed to solve the neutron transport equation. The method combines the advantages of two recently published schemes. The angular domain is discretized into small patches and patch-wise discontinuous angular basis functions are restricted to these patches, i.e. there is no overlap between basis functions corresponding to different patches. This approach yields block diagonal Jacobians with small block size and retains the possibility for S n -like approximate sweeping of the spatially discontinuous elements in order to provide efficient preconditioners for the solution procedure. On the other hand, the preservation of the full FEM framework (as opposed to collocation into a high-order S n scheme) retains the possibility of the Galerkin interpolated connection between phase space elements at arbitrary levels of discretization. Since the basis vectors are not orthonormal, a generalization of the Riemann procedure is introduced to separate the incoming and outgoing contributions in case of unstructured meshes. However, due to the properties of the angular discretization, the Riemann procedure can be avoided at a large fraction of the faces and this fraction rapidly increases as the level of refinement increases, contributing to the computational efficiency. In this paper the properties of the discretization scheme are studied with uniform refinement using an iterative solver based on the S 2 sweep order of the spatial elements. The fourth order convergence of the scalar flux is shown as anticipated from earlier schemes and the rapidly decreasing fraction of required Riemann faces is illustrated. (authors)

SAFE-PLANE, Stress Analysis of Planar Structure by Finite Elements Method

International Nuclear Information System (INIS)

Cornell, D.C.; Reich, Morris

1967-01-01

1 - Description of problem or function: SAFE-PLANE is applied to two- dimensional structures of arbitrary geometry under in-plane loads. Either plane stress or plane strain conditions may be imposed. Mechanical and thermal loads are permitted. 2 - Method of solution: The finite-element method is used to construct a mathematical model by assembling discrete elements. The total potential energy of the structure is determined and subsequently minimized by iteration on components of the displacement field until static equilibrium of the structure is attained. Strains and stresses are computed from the resulting displacements. 3 - Restrictions on the complexity of the problem: Multi-material structures with varying rigidities converge very slowly. Not valid for incompressible materials. Maximum number of nodal points = 675. Maximum number of elements = 1350

Finite Element Reliability Analysis of Chloride Ingress into Reinforced Concrete Structures

DEFF Research Database (Denmark)

Frier, Christian; Sørensen, John Dalsgaard

2007-01-01

For many reinforced concrete structures corrosion of the reinforcement is an important problem since it can result in maintenance and repair actions. Further, a reduction of the load-bearing capacity can occur. In the present paper the Finite Element Reliability Method (FERM) is employed for obta......For many reinforced concrete structures corrosion of the reinforcement is an important problem since it can result in maintenance and repair actions. Further, a reduction of the load-bearing capacity can occur. In the present paper the Finite Element Reliability Method (FERM) is employed...... concentration and reinforcement cover depth are modelled by stochastic fields, which are discretized using the Expansion Optimum Linear Estimation (EOLE) approach. The response gradients needed for FORM analysis are derived analytically using the Direct Differentiation Method (DDM). As an example, a bridge pier...... in a marine environment is considered and the results are given in terms of distributions of time for initiation of corrosion....

On the solvability of asymmetric quasilinear finite element approximate problems in nonlinear incompressible elasticity

International Nuclear Information System (INIS)

Ruas, V.

1982-09-01

A class of simplicial finite elements for solving incompressible elasticity problems in n-dimensional space, n=2 or 3, is presented. An asymmetric structure of the shape functions with respect to the centroid of the simplex, renders them particularly stable in the large strain case, in which the incompressibility condition is nonlinear. It is proved that under certain assembling conditions of the elements, there exists a solution to the corresponding discrete problems. Numerical examples illustrate the efficiency of the method. (Author) [pt

Finite difference discretization of semiconductor drift-diffusion equations for nanowire solar cells

Science.gov (United States)

Deinega, Alexei; John, Sajeev

2012-10-01

We introduce a finite difference discretization of semiconductor drift-diffusion equations using cylindrical partial waves. It can be applied to describe the photo-generated current in radial pn-junction nanowire solar cells. We demonstrate that the cylindrically symmetric (l=0) partial wave accurately describes the electronic response of a square lattice of silicon nanowires at normal incidence. We investigate the accuracy of our discretization scheme by using different mesh resolution along the radial direction r and compare with 3D (x, y, z) discretization. We consider both straight nanowires and nanowires with radius modulation along the vertical axis. The charge carrier generation profile inside each nanowire is calculated using an independent finite-difference time-domain simulation.

Application of Dynamic Analysis in Semi-Analytical Finite Element Method.

Science.gov (United States)

Liu, Pengfei; Xing, Qinyan; Wang, Dawei; Oeser, Markus

2017-08-30

Analyses of dynamic responses are significantly important for the design, maintenance and rehabilitation of asphalt pavement. In order to evaluate the dynamic responses of asphalt pavement under moving loads, a specific computational program, SAFEM, was developed based on a semi-analytical finite element method. This method is three-dimensional and only requires a two-dimensional FE discretization by incorporating Fourier series in the third dimension. In this paper, the algorithm to apply the dynamic analysis to SAFEM was introduced in detail. Asphalt pavement models under moving loads were built in the SAFEM and commercial finite element software ABAQUS to verify the accuracy and efficiency of the SAFEM. The verification shows that the computational accuracy of SAFEM is high enough and its computational time is much shorter than ABAQUS. Moreover, experimental verification was carried out and the prediction derived from SAFEM is consistent with the measurement. Therefore, the SAFEM is feasible to reliably predict the dynamic response of asphalt pavement under moving loads, thus proving beneficial to road administration in assessing the pavement's state.

The Mixed Finite Element Multigrid Method for Stokes Equations

Science.gov (United States)

Muzhinji, K.; Shateyi, S.; Motsa, S. S.

2015-01-01

The stable finite element discretization of the Stokes problem produces a symmetric indefinite system of linear algebraic equations. A variety of iterative solvers have been proposed for such systems in an attempt to construct efficient, fast, and robust solution techniques. This paper investigates one of such iterative solvers, the geometric multigrid solver, to find the approximate solution of the indefinite systems. The main ingredient of the multigrid method is the choice of an appropriate smoothing strategy. This study considers the application of different smoothers and compares their effects in the overall performance of the multigrid solver. We study the multigrid method with the following smoothers: distributed Gauss Seidel, inexact Uzawa, preconditioned MINRES, and Braess-Sarazin type smoothers. A comparative study of the smoothers shows that the Braess-Sarazin smoothers enhance good performance of the multigrid method. We study the problem in a two-dimensional domain using stable Hood-Taylor Q 2-Q 1 pair of finite rectangular elements. We also give the main theoretical convergence results. We present the numerical results to demonstrate the efficiency and robustness of the multigrid method and confirm the theoretical results. PMID:25945361

Probabilistic finite elements for fatigue and fracture analysis

Science.gov (United States)

Belytschko, Ted; Liu, Wing Kam

1993-04-01

An overview of the probabilistic finite element method (PFEM) developed by the authors and their colleagues in recent years is presented. The primary focus is placed on the development of PFEM for both structural mechanics problems and fracture mechanics problems. The perturbation techniques are used as major tools for the analytical derivation. The following topics are covered: (1) representation and discretization of random fields; (2) development of PFEM for the general linear transient problem and nonlinear elasticity using Hu-Washizu variational principle; (3) computational aspects; (4) discussions of the application of PFEM to the reliability analysis of both brittle fracture and fatigue; and (5) a stochastic computational tool based on stochastic boundary element (SBEM). Results are obtained for the reliability index and corresponding probability of failure for: (1) fatigue crack growth; (2) defect geometry; (3) fatigue parameters; and (4) applied loads. These results show that initial defect is a critical parameter.

Finite element procedures for coupled linear analysis of heat transfer, fluid and solid mechanics

Science.gov (United States)

Sutjahjo, Edhi; Chamis, Christos C.

1993-01-01

Coupled finite element formulations for fluid mechanics, heat transfer, and solid mechanics are derived from the conservation laws for energy, mass, and momentum. To model the physics of interactions among the participating disciplines, the linearized equations are coupled by combining domain and boundary coupling procedures. Iterative numerical solution strategy is presented to solve the equations, with the partitioning of temporal discretization implemented.

FINITE ELEMENT MODEL FOR PREDICTING RESIDUAL ...

African Journals Online (AJOL)

FINITE ELEMENT MODEL FOR PREDICTING RESIDUAL STRESSES IN ... the transverse residual stress in the x-direction (σx) had a maximum value of 375MPa ... the finite element method are in fair agreement with the experimental results.

A high-order multiscale finite-element method for time-domain acoustic-wave modeling

Science.gov (United States)

Gao, Kai; Fu, Shubin; Chung, Eric T.

2018-05-01

Accurate and efficient wave equation modeling is vital for many applications in such as acoustics, electromagnetics, and seismology. However, solving the wave equation in large-scale and highly heterogeneous models is usually computationally expensive because the computational cost is directly proportional to the number of grids in the model. We develop a novel high-order multiscale finite-element method to reduce the computational cost of time-domain acoustic-wave equation numerical modeling by solving the wave equation on a coarse mesh based on the multiscale finite-element theory. In contrast to existing multiscale finite-element methods that use only first-order multiscale basis functions, our new method constructs high-order multiscale basis functions from local elliptic problems which are closely related to the Gauss-Lobatto-Legendre quadrature points in a coarse element. Essentially, these basis functions are not only determined by the order of Legendre polynomials, but also by local medium properties, and therefore can effectively convey the fine-scale information to the coarse-scale solution with high-order accuracy. Numerical tests show that our method can significantly reduce the computation time while maintain high accuracy for wave equation modeling in highly heterogeneous media by solving the corresponding discrete system only on the coarse mesh with the new high-order multiscale basis functions.

Bridges for Pedestrians with Random Parameters using the Stochastic Finite Elements Analysis

Science.gov (United States)

Szafran, J.; Kamiński, M.

2017-02-01

The main aim of this paper is to present a Stochastic Finite Element Method analysis with reference to principal design parameters of bridges for pedestrians: eigenfrequency and deflection of bridge span. They are considered with respect to random thickness of plates in boxed-section bridge platform, Young modulus of structural steel and static load resulting from crowd of pedestrians. The influence of the quality of the numerical model in the context of traditional FEM is shown also on the example of a simple steel shield. Steel structures with random parameters are discretized in exactly the same way as for the needs of traditional Finite Element Method. Its probabilistic version is provided thanks to the Response Function Method, where several numerical tests with random parameter values varying around its mean value enable the determination of the structural response and, thanks to the Least Squares Method, its final probabilistic moments.

Bridges for Pedestrians with Random Parameters using the Stochastic Finite Elements Analysis

Directory of Open Access Journals (Sweden)

Szafran J.

2017-02-01

Full Text Available The main aim of this paper is to present a Stochastic Finite Element Method analysis with reference to principal design parameters of bridges for pedestrians: eigenfrequency and deflection of bridge span. They are considered with respect to random thickness of plates in boxed-section bridge platform, Young modulus of structural steel and static load resulting from crowd of pedestrians. The influence of the quality of the numerical model in the context of traditional FEM is shown also on the example of a simple steel shield. Steel structures with random parameters are discretized in exactly the same way as for the needs of traditional Finite Element Method. Its probabilistic version is provided thanks to the Response Function Method, where several numerical tests with random parameter values varying around its mean value enable the determination of the structural response and, thanks to the Least Squares Method, its final probabilistic moments.

Finite element methods a practical guide

CERN Document Server

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.

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)

Finite rotation shells basic equations and finite elements for Reissner kinematics

CERN Document Server

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.

Dynamic analysis of suspension cable based on vector form intrinsic finite element method

Science.gov (United States)

Qin, Jian; Qiao, Liang; Wan, Jiancheng; Jiang, Ming; Xia, Yongjun

2017-10-01

A vector finite element method is presented for the dynamic analysis of cable structures based on the vector form intrinsic finite element (VFIFE) and mechanical properties of suspension cable. Firstly, the suspension cable is discretized into different elements by space points, the mass and external forces of suspension cable are transformed into space points. The structural form of cable is described by the space points at different time. The equations of motion for the space points are established according to the Newton’s second law. Then, the element internal forces between the space points are derived from the flexible truss structure. Finally, the motion equations of space points are solved by the central difference method with reasonable time integration step. The tangential tension of the bearing rope in a test ropeway with the moving concentrated loads is calculated and compared with the experimental data. The results show that the tangential tension of suspension cable with moving loads is consistent with the experimental data. This method has high calculated precision and meets the requirements of engineering application.

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)

Finite element method for time-space-fractional Schrodinger equation

Directory of Open Access Journals (Sweden)

Xiaogang Zhu

2017-07-01

Full Text Available In this article, we develop a fully discrete finite element method for the nonlinear Schrodinger equation (NLS with time- and space-fractional derivatives. The time-fractional derivative is described in Caputo's sense and the space-fractional derivative in Riesz's sense. Its stability is well derived; the convergent estimate is discussed by an orthogonal operator. We also extend the method to the two-dimensional time-space-fractional NLS and to avoid the iterative solvers at each time step, a linearized scheme is further conducted. Several numerical examples are implemented finally, which confirm the theoretical results as well as illustrate the accuracy of our methods.

Finite approximations in discrete-time stochastic control quantized models and asymptotic optimality

CERN Document Server

Saldi, Naci; Yüksel, Serdar

2018-01-01

In a unified form, this monograph presents fundamental results on the approximation of centralized and decentralized stochastic control problems, with uncountable state, measurement, and action spaces. It demonstrates how quantization provides a system-independent and constructive method for the reduction of a system with Borel spaces to one with finite state, measurement, and action spaces. In addition to this constructive view, the book considers both the information transmission approach for discretization of actions, and the computational approach for discretization of states and actions. Part I of the text discusses Markov decision processes and their finite-state or finite-action approximations, while Part II builds from there to finite approximations in decentralized stochastic control problems. This volume is perfect for researchers and graduate students interested in stochastic controls. With the tools presented, readers will be able to establish the convergence of approximation models to original mo...

A Wavelet-Based Finite Element Method for the Self-Shielding Issue in Neutron Transport

International Nuclear Information System (INIS)

Le Tellier, R.; Fournier, D.; Ruggieri, J. M.

2009-01-01

This paper describes a new approach for treating the energy variable of the neutron transport equation in the resolved resonance energy range. The aim is to avoid recourse to a case-specific spatially dependent self-shielding calculation when considering a broad group structure. This method consists of a discontinuous Galerkin discretization of the energy using wavelet-based elements. A Σ t -orthogonalization of the element basis is presented in order to make the approach tractable for spatially dependent problems. First numerical tests of this method are carried out in a limited framework under the Livolant-Jeanpierre hypotheses in an infinite homogeneous medium. They are mainly focused on the way to construct the wavelet-based element basis. Indeed, the prior selection of these wavelet functions by a thresholding strategy applied to the discrete wavelet transform of a given quantity is a key issue for the convergence rate of the method. The Canuto thresholding approach applied to an approximate flux is found to yield a nearly optimal convergence in many cases. In these tests, the capability of such a finite element discretization to represent the flux depression in a resonant region is demonstrated; a relative accuracy of 10 -3 on the flux (in L 2 -norm) is reached with less than 100 wavelet coefficients per group. (authors)

From Laser Scanning to Finite Element Analysis of Complex Buildings by Using a Semi-Automatic Procedure.

Science.gov (United States)

Castellazzi, Giovanni; D'Altri, Antonio Maria; Bitelli, Gabriele; Selvaggi, Ilenia; Lambertini, Alessandro

2015-07-28

In this paper, a new semi-automatic procedure to transform three-dimensional point clouds of complex objects to three-dimensional finite element models is presented and validated. The procedure conceives of the point cloud as a stacking of point sections. The complexity of the clouds is arbitrary, since the procedure is designed for terrestrial laser scanner surveys applied to buildings with irregular geometry, such as historical buildings. The procedure aims at solving the problems connected to the generation of finite element models of these complex structures by constructing a fine discretized geometry with a reduced amount of time and ready to be used with structural analysis. If the starting clouds represent the inner and outer surfaces of the structure, the resulting finite element model will accurately capture the whole three-dimensional structure, producing a complex solid made by voxel elements. A comparison analysis with a CAD-based model is carried out on a historical building damaged by a seismic event. The results indicate that the proposed procedure is effective and obtains comparable models in a shorter time, with an increased level of automation.

Summary compilation of shell element performance versus formulation.

Energy Technology Data Exchange (ETDEWEB)

Heinstein, Martin Wilhelm; Hales, Jason Dean (Idaho National Laboratory, Idaho Falls, ID); Breivik, Nicole L.; Key, Samuel W. (FMA Development, LLC, Great Falls, MT)

2011-07-01

This document compares the finite element shell formulations in the Sierra Solid Mechanics code. These are finite elements either currently in the Sierra simulation codes Presto and Adagio, or expected to be added to them in time. The list of elements are divided into traditional two-dimensional, plane stress shell finite elements, and three-dimensional solid finite elements that contain either modifications or additional terms designed to represent the bending stiffness expected to be found in shell formulations. These particular finite elements are formulated for finite deformation and inelastic material response, and, as such, are not based on some of the elegant formulations that can be found in an elastic, infinitesimal finite element setting. Each shell element is subjected to a series of 12 verification and validation test problems. The underlying purpose of the tests here is to identify the quality of both the spatially discrete finite element gradient operator and the spatially discrete finite element divergence operator. If the derivation of the finite element is proper, the discrete divergence operator is the transpose of the discrete gradient operator. An overall summary is provided from which one can rank, at least in an average sense, how well the individual formulations can be expected to perform in applications encountered year in and year out. A letter grade has been assigned albeit sometimes subjectively for each shell element and each test problem result. The number of A's, B's, C's, et cetera assigned have been totaled, and a grade point average (GPA) has been computed, based on a 4.0-system. These grades, combined with a comparison between the test problems and the application problem, can be used to guide an analyst to select the element with the best shell formulation.

A finite element model for the quench front evolution problem

International Nuclear Information System (INIS)

Folescu, J.; Galeao, A.C.N.R.; Carmo, E.G.D. do.

1985-01-01

A model for the rewetting problem associated with the loss of coolant accident in a PWR reactor is proposed. A variational formulation for the time-dependent heat conduction problem on fuel rod cladding is used, and appropriate boundary conditions are assumed in order to simulate the thermal interaction between the fuel rod cladding and the fluid. A numerical procedure which uses the finite element method for the spatial discretization and a Crank-Nicolson-like method for the step-by-step integration is developed. Some numerical results are presented showing the quench front evolution and its stationary profile. (Author) [pt

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.

Why do probabilistic finite element analysis ?

CERN Document Server

Thacker, Ben H

2008-01-01

The intention of this book is to provide an introduction to performing probabilistic finite element analysis. As a short guideline, the objective is to inform the reader of the use, benefits and issues associated with performing probabilistic finite element analysis without excessive theory or mathematical detail.

An efficient finite element solution for gear dynamics

International Nuclear Information System (INIS)

Cooley, C G; Parker, R G; Vijayakar, S M

2010-01-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.

Finite element analysis theory and application with ANSYS

CERN Document Server

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

Finite element application to global reactor analysis

International Nuclear Information System (INIS)

Schmidt, F.A.R.

1981-01-01

The Finite Element Method is described as a Coarse Mesh Method with general basis and trial functions. Various consequences concerning programming and application of Finite Element Methods in reactor physics are drawn. One of the conclusions is that the Finite Element Method is a valuable tool in solving global reactor analysis problems. However, problems which can be described by rectangular boxes still can be solved with special coarse mesh programs more efficiently. (orig.) [de

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

Books and monographs on finite element technology

Science.gov (United States)

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.

Probabilistic finite elements for fracture mechanics

Science.gov (United States)

Besterfield, Glen

1988-01-01

The probabilistic finite element method (PFEM) is developed for probabilistic fracture mechanics (PFM). A finite element which has the near crack-tip singular strain embedded in the element is used. Probabilistic distributions, such as expectation, covariance and correlation stress intensity factors, are calculated for random load, random material and random crack length. The method is computationally quite efficient and can be expected to determine the probability of fracture or reliability.

Simplified Qualitative Discrete Numerical Model to Determine Cracking Pattern in Brittle Materials by Means of Finite Element Method

OpenAIRE

Ochoa-Avendaño, J.; Garzon-Alvarado, D. A.; Linero, Dorian L.; Cerrolaza, M.

2017-01-01

This paper presents the formulation, implementation, and validation of a simplified qualitative model to determine the crack path of solids considering static loads, infinitesimal strain, and plane stress condition. This model is based on finite element method with a special meshing technique, where nonlinear link elements are included between the faces of the linear triangular elements. The stiffness loss of some link elements represents the crack opening. Three experimental tests of bending...

Discrete differential geometry. Consistency as integrability

OpenAIRE

Bobenko, Alexander I.; Suris, Yuri B.

2005-01-01

A new field of discrete differential geometry is presently emerging on the border between differential and discrete geometry. Whereas classical differential geometry investigates smooth geometric shapes (such as surfaces), and discrete geometry studies geometric shapes with finite number of elements (such as polyhedra), the discrete differential geometry aims at the development of discrete equivalents of notions and methods of smooth surface theory. Current interest in this field derives not ...

Efficient CUDA Polynomial Preconditioned Conjugate Gradient Solver for Finite Element Computation of Elasticity Problems

Directory of Open Access Journals (Sweden)

Jianfei Zhang

2013-01-01

Full Text Available Graphics processing unit (GPU has obtained great success in scientific computations for its tremendous computational horsepower and very high memory bandwidth. This paper discusses the efficient way to implement polynomial preconditioned conjugate gradient solver for the finite element computation of elasticity on NVIDIA GPUs using compute unified device architecture (CUDA. Sliced block ELLPACK (SBELL format is introduced to store sparse matrix arising from finite element discretization of elasticity with fewer padding zeros than traditional ELLPACK-based formats. Polynomial preconditioning methods have been investigated both in convergence and running time. From the overall performance, the least-squares (L-S polynomial method is chosen as a preconditioner in PCG solver to finite element equations derived from elasticity for its best results on different example meshes. In the PCG solver, mixed precision algorithm is used not only to reduce the overall computational, storage requirements and bandwidth but to make full use of the capacity of the GPU devices. With SBELL format and mixed precision algorithm, the GPU-based L-S preconditioned CG can get a speedup of about 7–9 to CPU-implementation.

Discrete and mesoscopic regimes of finite-size wave turbulence

International Nuclear Information System (INIS)

L'vov, V. S.; Nazarenko, S.

2010-01-01

Bounding volume results in discreteness of eigenmodes in wave systems. This leads to a depletion or complete loss of wave resonances (three-wave, four-wave, etc.), which has a strong effect on wave turbulence (WT) i.e., on the statistical behavior of broadband sets of weakly nonlinear waves. This paper describes three different regimes of WT realizable for different levels of the wave excitations: discrete, mesoscopic and kinetic WT. Discrete WT comprises chaotic dynamics of interacting wave 'clusters' consisting of discrete (often finite) number of connected resonant wave triads (or quarters). Kinetic WT refers to the infinite-box theory, described by well-known wave-kinetic equations. Mesoscopic WT is a regime in which either the discrete and the kinetic evolutions alternate or when none of these two types is purely realized. We argue that in mesoscopic systems the wave spectrum experiences a sandpile behavior. Importantly, the mesoscopic regime is realized for a broad range of wave amplitudes which typically spans over several orders on magnitude, and not just for a particular intermediate level.

A Finite Element Theory for Predicting the Attenuation of Extended-Reacting Liners

Science.gov (United States)

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.

Introduction to finite and spectral element methods using Matlab

CERN Document Server

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.

A novel upwind stabilized discontinuous finite element angular framework for deterministic dose calculations in magnetic fields.

Science.gov (United States)

Yang, R; Zelyak, O; Fallone, B G; St-Aubin, J

2018-01-30

Angular discretization impacts nearly every aspect of a deterministic solution to the linear Boltzmann transport equation, especially in the presence of magnetic fields, as modeled by a streaming operator in angle. In this work a novel stabilization treatment of the magnetic field term is developed for an angular finite element discretization on the unit sphere, specifically involving piecewise partitioning of path integrals along curved element edges into uninterrupted segments of incoming and outgoing flux, with outgoing components updated iteratively. Correct order-of-accuracy for this angular framework is verified using the method of manufactured solutions for linear, quadratic, and cubic basis functions in angle. Higher order basis functions were found to reduce the error especially in strong magnetic fields and low density media. We combine an angular finite element mesh respecting octant boundaries on the unit sphere to spatial Cartesian voxel elements to guarantee an unambiguous transport sweep ordering in space. Accuracy for a dosimetrically challenging scenario involving bone and air in the presence of a 1.5 T parallel magnetic field is validated against the Monte Carlo package GEANT4. Accuracy and relative computational efficiency were investigated for various angular discretization parameters. 32 angular elements with quadratic basis functions yielded a reasonable compromise, with gamma passing rates of 99.96% (96.22%) for a 2%/2 mm (1%/1 mm) criterion. A rotational transformation of the spatial calculation geometry is performed to orient an arbitrary magnetic field vector to be along the z-axis, a requirement for a constant azimuthal angular sweep ordering. Working on the unit sphere, we apply the same rotational transformation to the angular domain to align its octants with the rotated Cartesian mesh. Simulating an oblique 1.5 T magnetic field against GEANT4 yielded gamma passing rates of 99.42% (95.45%) for a 2%/2 mm (1%/1 mm) criterion.

A novel upwind stabilized discontinuous finite element angular framework for deterministic dose calculations in magnetic fields

Science.gov (United States)

Yang, R.; Zelyak, O.; Fallone, B. G.; St-Aubin, J.

2018-02-01

Angular discretization impacts nearly every aspect of a deterministic solution to the linear Boltzmann transport equation, especially in the presence of magnetic fields, as modeled by a streaming operator in angle. In this work a novel stabilization treatment of the magnetic field term is developed for an angular finite element discretization on the unit sphere, specifically involving piecewise partitioning of path integrals along curved element edges into uninterrupted segments of incoming and outgoing flux, with outgoing components updated iteratively. Correct order-of-accuracy for this angular framework is verified using the method of manufactured solutions for linear, quadratic, and cubic basis functions in angle. Higher order basis functions were found to reduce the error especially in strong magnetic fields and low density media. We combine an angular finite element mesh respecting octant boundaries on the unit sphere to spatial Cartesian voxel elements to guarantee an unambiguous transport sweep ordering in space. Accuracy for a dosimetrically challenging scenario involving bone and air in the presence of a 1.5 T parallel magnetic field is validated against the Monte Carlo package GEANT4. Accuracy and relative computational efficiency were investigated for various angular discretization parameters. 32 angular elements with quadratic basis functions yielded a reasonable compromise, with gamma passing rates of 99.96% (96.22%) for a 2%/2 mm (1%/1 mm) criterion. A rotational transformation of the spatial calculation geometry is performed to orient an arbitrary magnetic field vector to be along the z-axis, a requirement for a constant azimuthal angular sweep ordering. Working on the unit sphere, we apply the same rotational transformation to the angular domain to align its octants with the rotated Cartesian mesh. Simulating an oblique 1.5 T magnetic field against GEANT4 yielded gamma passing rates of 99.42% (95.45%) for a 2%/2 mm (1%/1 mm) criterion.

Finite-time stability of discrete fractional delay systems: Gronwall inequality and stability criterion

Science.gov (United States)

Wu, Guo-Cheng; Baleanu, Dumitru; Zeng, Sheng-Da

2018-04-01

This study investigates finite-time stability of Caputo delta fractional difference equations. A generalized Gronwall inequality is given on a finite time domain. A finite-time stability criterion is proposed for fractional differential equations. Then the idea is extended to the discrete fractional case. A linear fractional difference equation with constant delays is considered and finite-time stable conditions are provided. One example is numerically illustrated to support the theoretical result.

Finite Element Methods and Their Applications

CERN Document Server

Chen, Zhangxin

2005-01-01

This book serves as a text for one- or two-semester courses for upper-level undergraduates and beginning graduate students and as a professional reference for people who want to solve partial differential equations (PDEs) using finite element methods. The author has attempted to introduce every concept in the simplest possible setting and maintain a level of treatment that is as rigorous as possible without being unnecessarily abstract. Quite a lot of attention is given to discontinuous finite elements, characteristic finite elements, and to the applications in fluid and solid mechanics including applications to porous media flow, and applications to semiconductor modeling. An extensive set of exercises and references in each chapter are provided.

NONLINEAR FINITE ELEMENT ANALYSIS OF NONSEISMICALLY DETAILED INTERIOR RC BEAM-COLUMN CONNECTION UNDER REVERSED CYCLIC LOAD

Directory of Open Access Journals (Sweden)

Teeraphot Supaviriyakit

2017-11-01

Full Text Available This paper presents a nonlinear finite element analysis of non-seismically detailed RC beam column connections under reversed cyclic load. The test of half-scale nonductile reinforced concrete beam-column joints was conducted. The tested specimens represented those of the actual mid-rise reinforced concrete frame buildings designed according to the non-seismic provisions of the ACI building code.  The test results show that specimens representing small and medium column tributary area failed in brittle joint shear while specimen representing large column tributary area failed by ductile flexure though no ductile reinforcement details were provided. The nonlinear finite element analysis was applied to simulate the behavior of the specimens. The finite element analysis employs the smeared crack approach for modeling beam, column and joint, and employs the discrete crack approach for modeling the interface between beam and joint face. The nonlinear constitutive models of reinforced concrete elements consist of coupled tension-compression model to model normal force orthogonal and parallel to the crack and shear transfer model to capture the shear sliding mechanism. The FEM shows good comparison with test results in terms of load-displacement relations, hysteretic loops, cracking process and the failure mode of the tested specimens. The finite element analysis clarifies that the joint shear failure was caused by the collapse of principal diagonal concrete strut.

Programming the finite element method

CERN Document Server

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

Lumped impulses, discrete displacements and a moving load analysis

NARCIS (Netherlands)

Kok, A.W.M.

1997-01-01

Finite element models are usually presented as relations between lumped forces and discrete displacements. Mostly finite element models are found by the elaboration of the method of the virtual work - which is a special case of the Galerkin's variational principle -. By application of Galerkin's

Mixed Element Formulation for the Finite Element-Boundary Integral Method

National Research Council Canada - National Science Library

Meese, J; Kempel, L. C; Schneider, S. W

2006-01-01

A mixed element approach using right hexahedral elements and right prism elements for the finite element-boundary integral method is presented and discussed for the study of planar cavity-backed antennas...

A new Eulerian-Lagrangian finite element simulator for solute transport in discrete fracture-matrix systems

Energy Technology Data Exchange (ETDEWEB)

Birkholzer, J.; Karasaki, K. [Lawrence Berkeley National Lab., CA (United States). Earth Sciences Div.

1996-07-01

Fracture network simulators have extensively been used in the past for obtaining a better understanding of flow and transport processes in fractured rock. However, most of these models do not account for fluid or solute exchange between the fractures and the porous matrix, although diffusion into the matrix pores can have a major impact on the spreading of contaminants. In the present paper a new finite element code TRIPOLY is introduced which combines a powerful fracture network simulator with an efficient method to account for the diffusive interaction between the fractures and the adjacent matrix blocks. The fracture network simulator used in TRIPOLY features a mixed Lagrangian-Eulerian solution scheme for the transport in fractures, combined with an adaptive gridding technique to account for sharp concentration fronts. The fracture-matrix interaction is calculated with an efficient method which has been successfully used in the past for dual-porosity models. Discrete fractures and matrix blocks are treated as two different systems, and the interaction is modeled by introducing sink/source terms in both systems. It is assumed that diffusive transport in the matrix can be approximated as a one-dimensional process, perpendicular to the adjacent fracture surfaces. A direct solution scheme is employed to solve the coupled fracture and matrix equations. The newly developed combination of the fracture network simulator and the fracture-matrix interaction module allows for detailed studies of spreading processes in fractured porous rock. The authors present a sample application which demonstrate the codes ability of handling large-scale fracture-matrix systems comprising individual fractures and matrix blocks of arbitrary size and shape.

An A Posteriori Error Analysis of Mixed Finite Element Galerkin Approximations to Second Order Linear Parabolic Problems

KAUST Repository

Memon, Sajid; Nataraj, Neela; Pani, Amiya Kumar

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.

A discrete element and ray framework for rapid simulation of acoustical dispersion of microscale particulate agglomerations

Science.gov (United States)

Zohdi, T. I.

2016-03-01

In industry, particle-laden fluids, such as particle-functionalized inks, are constructed by adding fine-scale particles to a liquid solution, in order to achieve desired overall properties in both liquid and (cured) solid states. However, oftentimes undesirable particulate agglomerations arise due to some form of mutual-attraction stemming from near-field forces, stray electrostatic charges, process ionization and mechanical adhesion. For proper operation of industrial processes involving particle-laden fluids, it is important to carefully breakup and disperse these agglomerations. One approach is to target high-frequency acoustical pressure-pulses to breakup such agglomerations. The objective of this paper is to develop a computational model and corresponding solution algorithm to enable rapid simulation of the effect of acoustical pulses on an agglomeration composed of a collection of discrete particles. Because of the complex agglomeration microstructure, containing gaps and interfaces, this type of system is extremely difficult to mesh and simulate using continuum-based methods, such as the finite difference time domain or the finite element method. Accordingly, a computationally-amenable discrete element/discrete ray model is developed which captures the primary physical events in this process, such as the reflection and absorption of acoustical energy, and the induced forces on the particulate microstructure. The approach utilizes a staggered, iterative solution scheme to calculate the power transfer from the acoustical pulse to the particles and the subsequent changes (breakup) of the pulse due to the particles. Three-dimensional examples are provided to illustrate the approach.

Quasistatic field simulations based on finite elements and spectral methods applied to superconducting magnets

International Nuclear Information System (INIS)

Koch, Stephan

2009-01-01

This thesis is concerned with the numerical simulation of electromagnetic fields in the quasi-static approximation which is applicable in many practical cases. Main emphasis is put on higher-order finite element methods. Quasi-static applications can be found, e.g., in accelerator physics in terms of the design of magnets required for beam guidance, in power engineering as well as in high-voltage engineering. Especially during the first design and optimization phase of respective devices, numerical models offer a cheap alternative to the often costly assembly of prototypes. However, large differences in the magnitude of the material parameters and the geometric dimensions as well as in the time-scales of the electromagnetic phenomena involved lead to an unacceptably long simulation time or to an inadequately large memory requirement. Under certain circumstances, the simulation itself and, in turn, the desired design improvement becomes even impossible. In the context of this thesis, two strategies aiming at the extension of the range of application for numerical simulations based on the finite element method are pursued. The first strategy consists in parallelizing existing methods such that the computation can be distributed over several computers or cores of a processor. As a consequence, it becomes feasible to simulate a larger range of devices featuring more degrees of freedom in the numerical model than before. This is illustrated for the calculation of the electromagnetic fields, in particular of the eddy-current losses, inside a superconducting dipole magnet developed at the GSI Helmholtzzentrum fuer Schwerionenforschung as a part of the FAIR project. As the second strategy to improve the efficiency of numerical simulations, a hybrid discretization scheme exploiting certain geometrical symmetries is established. Using this method, a significant reduction of the numerical effort in terms of required degrees of freedom for a given accuracy is achieved. The

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.

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

Using Finite Element Method

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.

Incorporation of exact boundary conditions into a discontinuous galerkin finite element method for accurately solving 2d time-dependent maxwell equations

KAUST Repository

Sirenko, Kostyantyn; Liu, Meilin; Bagci, Hakan

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

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)

A discrete finite element modelling and measurements for powder compaction

International Nuclear Information System (INIS)

Choi, J L; Gethin, D T

2009-01-01

An experimental investigation into friction between powder and a target surface together with numerical modelling of compaction and friction processes at a micro-scale are presented in this paper. The experimental work explores friction mechanisms by using an extended sliding plate apparatus operating at low load while sliding over a long distance. Tests were conducted for copper and 316 steel with variation in loads, surface finish and its orientation. The behaviours of the static and dynamic friction were identified highlighting the important influence of particle size, particle shape, material response and surface topography. The results also highlighted that under light loading the friction coefficient remains at a level lower than that derived from experiments on equipment having a wider dynamic range and this is attributed to the enhanced sensitivity of the measurement equipment. The results also suggest that friction variation with sliding distance is a consequence of damage, rather than presentation of an uncontaminated target sliding surface. The complete experimental cycle was modelled numerically using a combined discrete and finite element scheme enabling exploration of mechanisms that are defined at the particle level. Using compaction as the starting point, a number of simulation factors and process parameters were investigated. Comparisons were made with previously published work, showing reasonable agreement and the simulations were then used to explore the process response to the range of particle scale factors. Models comprising regular packing of round particles exhibited stiff response with high initial density. Models with random packing were explored and were found to reflect trends that are more closely aligned with experimental observation, including rearrangement, followed by compaction under a regime of elastic then plastic deformation. Numerical modelling of the compaction stage was extended to account for the shearing stage of the

Explicit dynamics for numerical simulation of crack propagation by the extended finite element method

International Nuclear Information System (INIS)

Menouillard, T.

2007-09-01

Computerized simulation is nowadays an integrating part of design and validation processes of mechanical structures. Simulation tools are more and more performing allowing a very acute description of the phenomena. Moreover, these tools are not limited to linear mechanics but are developed to describe more difficult behaviours as for instance structures damage which interests the safety domain. A dynamic or static load can thus lead to a damage, a crack and then a rupture of the structure. The fast dynamics allows to simulate 'fast' phenomena such as explosions, shocks and impacts on structure. The application domain is various. It concerns for instance the study of the lifetime and the accidents scenario of the nuclear reactor vessel. It is then very interesting, for fast dynamics codes, to be able to anticipate in a robust and stable way such phenomena: the assessment of damage in the structure and the simulation of crack propagation form an essential stake. The extended finite element method has the advantage to break away from mesh generation and from fields projection during the crack propagation. Effectively, crack is described kinematically by an appropriate strategy of enrichment of supplementary freedom degrees. Difficulties connecting the spatial discretization of this method with the temporal discretization of an explicit calculation scheme has then been revealed; these difficulties are the diagonal writing of the mass matrix and the associated stability time step. Here are presented two methods of mass matrix diagonalization based on the kinetic energy conservation, and studies of critical time steps for various enriched finite elements. The interest revealed here is that the time step is not more penalizing than those of the standard finite elements problem. Comparisons with numerical simulations on another code allow to validate the theoretical works. A crack propagation test in mixed mode has been exploited in order to verify the simulation

ANSYS mechanical APDL for finite element analysis

CERN Document Server

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

Parallel algorithms for solving the diffusion equation by finite elements methods and by nodal methods

International Nuclear Information System (INIS)

Coulomb, F.

1989-06-01

The aim of this work is to study methods for solving the diffusion equation, based on a primal or mixed-dual finite elements discretization and well suited for use on multiprocessors computers; domain decomposition methods are the subject of the main part of this study, the linear systems being solved by the block-Jacobi method. The origin of the diffusion equation is explained in short, and various variational formulations are reminded. A survey of iterative methods is given. The elemination of the flux or current is treated in the case of a mixed method. Numerical tests are performed on two examples of reactors, in order to compare mixed elements and Lagrange elements. A theoretical study of domain decomposition is led in the case of Lagrange finite elements, and convergence conditions for the block-Jacobi method are derived; the dissection decomposition is previously the purpose of a particular numerical analysis. In the case of mixed-dual finite elements, a study is led on examples and is confirmed by numerical tests performed for the dissection decomposition; furthermore, after being justified, decompositions along axes of symmetry are numerically tested. In the case of a decomposition into two subdomains, the dissection decomposition and the decomposition with an integrated interface are compared. Alternative directions methods are defined; the convergence of those relative to Lagrange elements is shown; in the case of mixed elements, convergence conditions are found [fr

Space-time coupled spectral/hp least-squares finite element formulation for the incompressible Navier-Stokes equations

International Nuclear Information System (INIS)

Pontaza, J.P.; Reddy, J.N.

2004-01-01

We consider least-squares finite element models for the numerical solution of the non-stationary Navier-Stokes equations governing viscous incompressible fluid flows. The paper presents a formulation where the effects of space and time are coupled, resulting in a true space-time least-squares minimization procedure, as opposed to a space-time decoupled formulation where a least-squares minimization procedure is performed in space at each time step. The formulation is first presented for the linear advection-diffusion equation and then extended to the Navier-Stokes equations. The formulation has no time step stability restrictions and is spectrally accurate in both space and time. To allow the use of practical C 0 element expansions in the resulting finite element model, the Navier-Stokes equations are expressed as an equivalent set of first-order equations by introducing vorticity as an additional independent variable and the least-squares method is used to develop the finite element model of the governing equations. High-order element expansions are used to construct the discrete model. The discrete model thus obtained is linearized by Newton's method, resulting in a linear system of equations with a symmetric positive definite coefficient matrix that is solved in a fully coupled manner by a preconditioned conjugate gradient method in matrix-free form. Spectral convergence of the L 2 least-squares functional and L 2 error norms in space-time is verified using a smooth solution to the two-dimensional non-stationary incompressible Navier-Stokes equations. Numerical results are presented for impulsively started lid-driven cavity flow, oscillatory lid-driven cavity flow, transient flow over a backward-facing step, and flow around a circular cylinder; the results demonstrate the predictive capability and robustness of the proposed formulation. Even though the space-time coupled formulation is emphasized, we also present the formulation and numerical results for least

A non-conformal finite element/finite volume scheme for the non-structured grid-based approximation of low Mach number flows; Un schema elements finis non-conformes/volumes finis pour l'approximation en maillages non-structures des ecoulements a faible nombre de Mach

Energy Technology Data Exchange (ETDEWEB)

Ansanay-Alex, G.

2009-06-17

The development of simulation codes aimed at a precise simulation of fires requires a precise approach of flame front phenomena by using very fine grids. The need to take different spatial scale into consideration leads to a local grid refinement and to a discretization with homogeneous grid for computing time and memory purposes. The author reports the approximation of the non-linear convection term, the scalar advection-diffusion in finite volumes, numerical simulations of a flow in a bent tube, of a three-dimensional laminar flame and of a low Mach number an-isotherm flow. Non conformal finite elements are also presented (Rannacher-Turek and Crouzeix-Raviart elements)

A discontinous Galerkin finite element method with an efficient time integration scheme for accurate simulations

KAUST Repository

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.

Hierarchical Material Properties in Finite Element Analysis: The Oilfield Infrastructure Problem.

Science.gov (United States)

Weiss, C. J.; Wilson, G. A.

2017-12-01

Geophysical simulation of low-frequency electromagnetic signals within built environments such as urban centers and industrial landscapes facilities is a challenging computational problem because strong conductors (e.g., pipes, fences, rail lines, rebar, etc.) are not only highly conductive and/or magnetic relative to the surrounding geology, but they are very small in one or more of their physical length coordinates. Realistic modeling of such structures as idealized conductors has long been the standard approach; however this strategy carries with it computational burdens such as cumbersome implementation of internal boundary conditions, and limited flexibility for accommodating realistic geometries. Another standard approach is "brute force" discretization (often coupled with an equivalent medium model) whereby 100's of millions of voxels are used to represent these strong conductors, but at the cost of extreme computation times (and mesh design) for a simulation result when possible. To minimize these burdens, a new finite element scheme (Weiss, Geophysics, 2017) has been developed in which the material properties reside on a hierarchy of geometric simplicies (i.e., edges, facets and volumes) within an unstructured tetrahedral mesh. This allows thin sheet—like structures, such as subsurface fractures, to be economically represented by a connected set of triangular facets, for example, that freely conform to arbitrary "real world" geometries. The same holds thin pipe/wire-like structures, such as casings or pipelines. The hierarchical finite element scheme has been applied to problems in electro- and magnetostatics for oilfield problems where the elevated, but finite, conductivity and permeability of the steel-cased oil wells must be properly accounted for, yielding results that are otherwise unobtainable, with run times as low as a few 10s of seconds. Extension of the hierarchical finite element concept to broadband electromagnetics is presently underway, as

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

Low Reynolds number steady state flow through a branching network of rigid vessels: II. A finite element mixture model

NARCIS (Netherlands)

Huyghe, J.M.R.J.; Oomens, C.W.J.; Campen, van D.H.; Heethaar, R.M.

1989-01-01

This research aims at formulating and verifying a finite element mixture formulation for blood perfusion. The equations derived in a companion paper [3] are discretized according to the Galerkin method. A flow experiment in a rigid model of a vascular tree of about 500 vessels is performed in order

Mixed Finite Element Simulation with Stability Analysis for Gas Transport in Low-Permeability Reservoirs

Directory of Open Access Journals (Sweden)

Mohamed F. El-Amin

2018-01-01

Full Text Available Natural gas exists in considerable quantities in tight reservoirs. Tight formations are rocks with very tiny or poorly connected pors that make flow through them very difficult, i.e., the permeability is very low. The mixed finite element method (MFEM, which is locally conservative, is suitable to simulate the flow in porous media. This paper is devoted to developing a mixed finite element (MFE technique to simulate the gas transport in low permeability reservoirs. The mathematical model, which describes gas transport in low permeability formations, contains slippage effect, as well as adsorption and diffusion mechanisms. The apparent permeability is employed to represent the slippage effect in low-permeability formations. The gas adsorption on the pore surface has been described by Langmuir isotherm model, while the Peng-Robinson equation of state is used in the thermodynamic calculations. Important compatibility conditions must hold to guarantee the stability of the mixed method by adding additional constraints to the numerical discretization. The stability conditions of the MFE scheme has been provided. A theorem and three lemmas on the stability analysis of the mixed finite element method (MFEM have been established and proven. A semi-implicit scheme is developed to solve the governing equations. Numerical experiments are carried out under various values of the physical parameters.

Analyses of large quasistatic deformations of inelastic bodies by a new hybrid-stress finite element algorithm

Science.gov (United States)

Reed, K. W.; Atluri, S. N.

1983-01-01

A new hybrid-stress finite element algorithm, suitable for analyses of large, quasistatic, inelastic deformations, is presented. The algorithm is base upon a generalization of de Veubeke's complementary energy principle. The principal variables in the formulation are the nominal stress rate and spin, and thg resulting finite element equations are discrete versions of the equations of compatibility and angular momentum balance. The algorithm produces true rates, time derivatives, as opposed to 'increments'. There results a complete separation of the boundary value problem (for stress rate and velocity) and the initial value problem (for total stress and deformation); hence, their numerical treatments are essentially independent. After a fairly comprehensive discussion of the numerical treatment of the boundary value problem, we launch into a detailed examination of the numerical treatment of the initial value problem, covering the topics of efficiency, stability and objectivity. The paper is closed with a set of examples, finite homogeneous deformation problems, which serve to bring out important aspects of the algorithm.

A finite-element model in vorticity and current function for the numerical solution of the Navier-Stokes equations

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

A study on the nonlinear finite element analysis of reinforced concrete structures: shell finite element formulation

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.

A study on the nonlinear finite element analysis of reinforced concrete structures: shell finite element formulation

International Nuclear Information System (INIS)

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

New Multigrid Method Including Elimination Algolithm Based on High-Order Vector Finite Elements in Three Dimensional Magnetostatic Field Analysis

Science.gov (United States)

Hano, Mitsuo; Hotta, Masashi

A new multigrid method based on high-order vector finite elements is proposed in this paper. Low level discretizations in this method are obtained by using low-order vector finite elements for the same mesh. Gauss-Seidel method is used as a smoother, and a linear equation of lowest level is solved by ICCG method. But it is often found that multigrid solutions do not converge into ICCG solutions. An elimination algolithm of constant term using a null space of the coefficient matrix is also described. In three dimensional magnetostatic field analysis, convergence time and number of iteration of this multigrid method are discussed with the convectional ICCG method.

Finite Element Analysis of Pipe T-Joint

OpenAIRE

P.M.Gedkar; Dr. D.V. Bhope

2012-01-01

This paper reports stress analysis of two pressurized cylindrical intersection using finite element method. The different combinations of dimensions of run pipe and the branch pipe are used to investigate thestresses in pipe at the intersection. In this study the stress analysis is accomplished by finite element package ANSYS.

Introduction to finite element analysis using MATLAB and Abaqus

CERN Document Server

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

Final Report for DOE grant DE-FG02-07ER64432 "New Grid and Discretization Technologies for Ocean and Ice Simulations"

Energy Technology Data Exchange (ETDEWEB)

Gunzburger, Max

2013-03-12

The work reported is in pursuit of these goals: high-quality unstructured, non-uniform Voronoi and Delaunay grids; improved finite element and finite volume discretization schemes; and improved finite element and finite volume discretization schemes. These are sought for application to spherical and three-dimensional applications suitable for ocean, atmosphere, ice-sheet, and other climate modeling applications.

Nonlinear finite element modeling of corrugated board

Science.gov (United States)

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

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

The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the

Bibliography for finite elements. [2200 references

Energy Technology Data Exchange (ETDEWEB)

Whiteman, J R [comp.

1975-01-01

This bibliography cites almost all of the significant papers on advances in the mathematical theory of finite elements. Reported are applications in aeronautical, civil, mechanical, nautical and nuclear engineering. Such topics as classical analysis, functional analysis, approximation theory, fluids, and diffusion are covered. Over 2200 references to publications up to the end of 1974 are included. Publications are listed alphabetically by author and also by keywords. In addition, finite element packages are listed.

Complex finite element sensitivity method for creep analysis

International Nuclear Information System (INIS)

Gomez-Farias, Armando; Montoya, Arturo; Millwater, Harry

2015-01-01

The complex finite element method (ZFEM) has been extended to perform sensitivity analysis for mechanical and structural systems undergoing creep deformation. ZFEM uses a complex finite element formulation to provide shape, material, and loading derivatives of the system response, providing an insight into the essential factors which control the behavior of the system as a function of time. A complex variable-based quadrilateral user element (UEL) subroutine implementing the power law creep constitutive formulation was incorporated within the Abaqus commercial finite element software. The results of the complex finite element computations were verified by comparing them to the reference solution for the steady-state creep problem of a thick-walled cylinder in the power law creep range. A practical application of the ZFEM implementation to creep deformation analysis is the calculation of the skeletal point of a notched bar test from a single ZFEM run. In contrast, the standard finite element procedure requires multiple runs. The value of the skeletal point is that it identifies the location where the stress state is accurate, regardless of the certainty of the creep material properties. - Highlights: • A novel finite element sensitivity method (ZFEM) for creep was introduced. • ZFEM has the capability to calculate accurate partial derivatives. • ZFEM can be used for identification of the skeletal point of creep structures. • ZFEM can be easily implemented in a commercial software, e.g. Abaqus. • ZFEM results were shown to be in excellent agreement with analytical solutions

Dynamic mortar finite element method for modeling of shear rupture on frictional rough surfaces

Science.gov (United States)

Tal, Yuval; Hager, Bradford H.

2017-09-01

This paper presents a mortar-based finite element formulation for modeling the dynamics of shear rupture on rough interfaces governed by slip-weakening and rate and state (RS) friction laws, focusing on the dynamics of earthquakes. The method utilizes the dual Lagrange multipliers and the primal-dual active set strategy concepts, together with a consistent discretization and linearization of the contact forces and constraints, and the friction laws to obtain a semi-smooth Newton method. The discretization of the RS friction law involves a procedure to condense out the state variables, thus eliminating the addition of another set of unknowns into the system. Several numerical examples of shear rupture on frictional rough interfaces demonstrate the efficiency of the method and examine the effects of the different time discretization schemes on the convergence, energy conservation, and the time evolution of shear traction and slip rate.

A suitable low-order, eight-node tetrahedral finite element for solids

Energy Technology Data Exchange (ETDEWEB)

Key, S.W.; Heinstein, M.S.; Stone, C.M.; Mello, F.J.; Blanford, M.L.; Budge, K.G.

1998-03-01

To use the all-tetrahedral mesh generation existing today, the authors have explored the creation of a computationally efficient eight-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four mid-face nodal points). The derivation of the element`s gradient operator, studies in obtaining a suitable mass lumping, and the element`s performance in applications are presented. In particular they examine the eight-node tetrahedral finite element`s behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hour-glass modes. In this regard it bears similarities to the eight-node, mean-quadrature hexahedral finite element. Comparisons with the results obtained from the mean-quadrature eight-node hexahedral finite element and the four-node tetrahedral finite element are included. Given automatic all-tetrahedral meshing, the eight-node, constant-strain tetrahedral finite element is a suitable replacement for the eight-node hexahedral finite element in those cases where mesh generation requires an inordinate amount of user intervention and direction to obtain acceptable mesh properties.

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

Linear triangle finite element formulation for multigroup neutron transport analysis with anisotropic scattering

Energy Technology Data Exchange (ETDEWEB)

Lillie, R.A.; Robinson, J.C.

1976-05-01

The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. A finite element formulation, utilizing a canonical form of the transport equation, is here developed to obtain both integral and pointwise solutions to neutron transport problems. The formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included by employing discrete ordinates-like approximations. In addition, multigroup source outer iteration techniques are employed to perform group-dependent calculations. The ability of the formulation to reduce substantially ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal-type lattice. A small, high-leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation. 8 tables, 29 figures.

Linear triangle finite element formulation for multigroup neutron transport analysis with anisotropic scattering

International Nuclear Information System (INIS)

Lillie, R.A.; Robinson, J.C.

1976-05-01

The discrete ordinates method is the most powerful and generally used deterministic method to obtain approximate solutions of the Boltzmann transport equation. A finite element formulation, utilizing a canonical form of the transport equation, is here developed to obtain both integral and pointwise solutions to neutron transport problems. The formulation is based on the use of linear triangles. A general treatment of anisotropic scattering is included by employing discrete ordinates-like approximations. In addition, multigroup source outer iteration techniques are employed to perform group-dependent calculations. The ability of the formulation to reduce substantially ray effects and its ability to perform streaming calculations are demonstrated by analyzing a series of test problems. The anisotropic scattering and multigroup treatments used in the development of the formulation are verified by a number of one-dimensional comparisons. These comparisons also demonstrate the relative accuracy of the formulation in predicting integral parameters. The applicability of the formulation to nonorthogonal planar geometries is demonstrated by analyzing a hexagonal-type lattice. A small, high-leakage reactor model is analyzed to investigate the effects of varying both the spatial mesh and order of angular quadrature. This analysis reveals that these effects are more pronounced in the present formulation than in other conventional formulations. However, the insignificance of these effects is demonstrated by analyzing a realistic reactor configuration. In addition, this final analysis illustrates the importance of incorporating anisotropic scattering into the finite element formulation. 8 tables, 29 figures

The finite element method its basis and fundamentals

CERN Document Server

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

Analysis of a discrete element method and coupling with a compressible fluid flow method

International Nuclear Information System (INIS)

Monasse, L.

2011-01-01

This work aims at the numerical simulation of compressible fluid/deformable structure interactions. In particular, we have developed a partitioned coupling algorithm between a Finite Volume method for the compressible fluid and a Discrete Element method capable of taking into account fractures in the solid. A survey of existing fictitious domain methods and partitioned algorithms has led to choose an Embedded Boundary method and an explicit coupling scheme. We first showed that the Discrete Element method used for the solid yielded the correct macroscopic behaviour and that the symplectic time-integration scheme ensured the preservation of energy. We then developed an explicit coupling algorithm between a compressible inviscid fluid and an un-deformable solid. Mass, momentum and energy conservation and consistency properties were proved for the coupling scheme. The algorithm was then extended to the coupling with a deformable solid, in the form of a semi implicit scheme. Finally, we applied this method to unsteady inviscid flows around moving structures: comparisons with existing numerical and experimental results demonstrate the excellent accuracy of our method. (author) [fr

A suitable low-order, eight-node tetrahedral finite element for solids

International Nuclear Information System (INIS)

Key, S.W.; Heinstein, M.S.; Stone, C.M.; Mello, F.J.; Blanford, M.L.; Budge, K.G.

1998-03-01

To use the all-tetrahedral mesh generation existing today, the authors have explored the creation of a computationally efficient eight-node tetrahedral finite element (a four-node tetrahedral finite element enriched with four mid-face nodal points). The derivation of the element's gradient operator, studies in obtaining a suitable mass lumping, and the element's performance in applications are presented. In particular they examine the eight-node tetrahedral finite element's behavior in longitudinal plane wave propagation, in transverse cylindrical wave propagation, and in simulating Taylor bar impacts. The element samples only constant strain states and, therefore, has 12 hour-glass modes. In this regard it bears similarities to the eight-node, mean-quadrature hexahedral finite element. Comparisons with the results obtained from the mean-quadrature eight-node hexahedral finite element and the four-node tetrahedral finite element are included. Given automatic all-tetrahedral meshing, the eight-node, constant-strain tetrahedral finite element is a suitable replacement for the eight-node hexahedral finite element in those cases where mesh generation requires an inordinate amount of user intervention and direction to obtain acceptable mesh properties

A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems

KAUST Repository

Bao, Kai

2012-10-01

In this paper, a semi-implicit finite element method is presented for the coupled Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn-Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier-Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples. © 2012 Elsevier Inc..

A finite element method for the numerical solution of the coupled Cahn-Hilliard and Navier-Stokes system for moving contact line problems

KAUST Repository

Bao, Kai; Shi, Yi; Sun, Shuyu; Wang, Xiaoping

2012-01-01

In this paper, a semi-implicit finite element method is presented for the coupled Cahn-Hilliard and Navier-Stokes equations with the generalized Navier boundary condition for the moving contact line problems. In our method, the system is solved in a decoupled way. For the Cahn-Hilliard equations, a convex splitting scheme is used along with a P1-P1 finite element discretization. The scheme is unconditionally stable. A linearized semi-implicit P2-P0 mixed finite element method is employed to solve the Navier-Stokes equations. With our method, the generalized Navier boundary condition is extended to handle the moving contact line problems with complex boundary in a very natural way. The efficiency and capacity of the present method are well demonstrated with several numerical examples. © 2012 Elsevier Inc..

A locally conservative stabilized continuous Galerkin finite element method for two-phase flow in poroelastic subsurfaces

Science.gov (United States)

Deng, Q.; Ginting, V.; McCaskill, B.; Torsu, P.

2017-10-01

We study the application of a stabilized continuous Galerkin finite element method (CGFEM) in the simulation of multiphase flow in poroelastic subsurfaces. The system involves a nonlinear coupling between the fluid pressure, subsurface's deformation, and the fluid phase saturation, and as such, we represent this coupling through an iterative procedure. Spatial discretization of the poroelastic system employs the standard linear finite element in combination with a numerical diffusion term to maintain stability of the algebraic system. Furthermore, direct calculation of the normal velocities from pressure and deformation does not entail a locally conservative field. To alleviate this drawback, we propose an element based post-processing technique through which local conservation can be established. The performance of the method is validated through several examples illustrating the convergence of the method, the effectivity of the stabilization term, and the ability to achieve locally conservative normal velocities. Finally, the efficacy of the method is demonstrated through simulations of realistic multiphase flow in poroelastic subsurfaces.

Application of finite-element-methods in food processing

DEFF Research Database (Denmark)

Risum, Jørgen

2004-01-01

Presentation of the possible use of finite-element-methods in food processing. Examples from diffusion studies are given.......Presentation of the possible use of finite-element-methods in food processing. Examples from diffusion studies are given....

Interpolation of the discrete logarithm in a finite field of characteristic two by Boolean functions

DEFF Research Database (Denmark)

Brandstaetter, Nina; Lange, Tanja; Winterhof, Arne

2005-01-01

We obtain bounds on degree, weight, and the maximal Fourier coefficient of Boolean functions interpolating the discrete logarithm in finite fields of characteristic two. These bounds complement earlier results for finite fields of odd characteristic....

Parametric Design Optimization Of A Novel Permanent Magnet Coupling Using Finite Element Analysis

DEFF Research Database (Denmark)

Högberg, Stig; Mijatovic, Nenad; Holbøll, Joachim

2014-01-01

A parametric design optimization routine has been applied to a novel magnetic coupling with improved recyclability. Coupling designs are modeled in a 3-D finite element environ- ment, and evaluated by three design objectives: pull-out torque, torque density by magnet mass, and torque density...... by total mass. Magnet and outer core thicknesses are varied discretely, whereas outer dimensions and air-gap length are kept constant. Comparative trends as a function of pole number and dimensions are depicted. A compromise exist between the design objectives, in which favoring one might reduce the other...

Residual-based a posteriori error estimation for multipoint flux mixed finite element methods

KAUST Repository

Du, Shaohong; Sun, Shuyu; Xie, Xiaoping

2015-01-01

A novel residual-type a posteriori error analysis technique is developed for multipoint flux mixed finite element methods for flow in porous media in two or three space dimensions. The derived a posteriori error estimator for the velocity and pressure error in L-norm consists of discretization and quadrature indicators, and is shown to be reliable and efficient. The main tools of analysis are a locally postprocessed approximation to the pressure solution of an auxiliary problem and a quadrature error estimate. Numerical experiments are presented to illustrate the competitive behavior of the estimator.

Residual-based a posteriori error estimation for multipoint flux mixed finite element methods

KAUST Repository

Du, Shaohong

2015-10-26

A novel residual-type a posteriori error analysis technique is developed for multipoint flux mixed finite element methods for flow in porous media in two or three space dimensions. The derived a posteriori error estimator for the velocity and pressure error in L-norm consists of discretization and quadrature indicators, and is shown to be reliable and efficient. The main tools of analysis are a locally postprocessed approximation to the pressure solution of an auxiliary problem and a quadrature error estimate. Numerical experiments are presented to illustrate the competitive behavior of the estimator.

Guaranteed Cost Finite-Time Control of Discrete-Time Positive Impulsive Switched Systems

Directory of Open Access Journals (Sweden)

Leipo Liu

2018-01-01

Full Text Available This paper considers the guaranteed cost finite-time boundedness of discrete-time positive impulsive switched systems. Firstly, the definition of guaranteed cost finite-time boundedness is introduced. By using the multiple linear copositive Lyapunov function (MLCLF and average dwell time (ADT approach, a state feedback controller is designed and sufficient conditions are obtained to guarantee that the corresponding closed-loop system is guaranteed cost finite-time boundedness (GCFTB. Such conditions can be solved by linear programming. Finally, a numerical example is provided to show the effectiveness of the proposed method.

Benders decomposition for discrete material optimization in laminate design with local failure criteria

DEFF Research Database (Denmark)

Munoz, Eduardo; Stolpe, Mathias; Bendsøe, Martin P.

2009-01-01

in any discrete angle optimization design, or material selection problems. The mathematical modeling of this problem is more general than the one of standard topology optimization. When considering only two material candidates with a considerable difference in stiffness, it corresponds exactly...... to a topology optimization problem. The problem is modeled as a discrete design problem coming from a finite element discretization of the continuum problem. This discretization is made of shell or plate elements. For each element (selection domain), only one of the material candidates must be selected...... of the relaxed master problem and the current best compliance (weight) found get close enough with respect to certain tolerance. The method is investigated by computational means, using the finite element method to solve the analysis problems, and a commercial branch and cut method for solving the relaxed master...

Advances in 3D electromagnetic finite element modeling

International Nuclear Information System (INIS)

Nelson, E.M.

1997-01-01

Numerous advances in electromagnetic finite element analysis (FEA) have been made in recent years. The maturity of frequency domain and eigenmode calculations, and the growth of time domain applications is briefly reviewed. A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will also be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis is also discussed

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

Ballistic calculation of nonequilibrium Green's function in nanoscale devices using finite element method

International Nuclear Information System (INIS)

Kurniawan, O; Bai, P; Li, E

2009-01-01

A ballistic calculation of a full quantum mechanical system is presented to study 2D nanoscale devices. The simulation uses the nonequilibrium Green's function (NEGF) approach to calculate the transport properties of the devices. While most available software uses the finite difference discretization technique, our work opts to formulate the NEGF calculation using the finite element method (FEM). In calculating a ballistic device, the FEM gives some advantages. In the FEM, the floating boundary condition for ballistic devices is satisfied naturally. This paper gives a detailed finite element formulation of the NEGF calculation applied to a double-gate MOSFET device with a channel length of 10 nm and a body thickness of 3 nm. The potential, electron density, Fermi functions integrated over the transverse energy, local density of states and the transmission coefficient of the device have been studied. We found that the transmission coefficient is significantly affected by the top of the barrier between the source and the channel, which in turn depends on the gate control. This supports the claim that ballistic devices can be modelled by the transport properties at the top of the barrier. Hence, the full quantum mechanical calculation presented here confirms the theory of ballistic transport in nanoscale devices.

Finite-element analysis of dynamic fracture

Science.gov (United States)

Aberson, J. A.; Anderson, J. M.; King, W. W.

1976-01-01

Applications of the finite element method to the two dimensional elastodynamics of cracked structures are presented. Stress intensity factors are computed for two problems involving stationary cracks. The first serves as a vehicle for discussing lumped-mass and consistent-mass characterizations of inertia. In the second problem, the behavior of a photoelastic dynamic tear test specimen is determined for the time prior to crack propagation. Some results of a finite element simulation of rapid crack propagation in an infinite body are discussed.

The finite element method in electromagnetics

CERN Document Server

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

Finite elements for analysis and design

CERN Document Server

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

A locally conservative non-negative finite element formulation for anisotropic advective-diffusive-reactive systems

Science.gov (United States)

Mudunuru, M. K.; Shabouei, M.; Nakshatrala, K.

2015-12-01

Advection-diffusion-reaction (ADR) equations appear in various areas of life sciences, hydrogeological systems, and contaminant transport. Obtaining stable and accurate numerical solutions can be challenging as the underlying equations are coupled, nonlinear, and non-self-adjoint. Currently, there is neither a robust computational framework available nor a reliable commercial package known that can handle various complex situations. Herein, the objective of this poster presentation is to present a novel locally conservative non-negative finite element formulation that preserves the underlying physical and mathematical properties of a general linear transient anisotropic ADR equation. In continuous setting, governing equations for ADR systems possess various important properties. In general, all these properties are not inherited during finite difference, finite volume, and finite element discretizations. The objective of this poster presentation is two fold: First, we analyze whether the existing numerical formulations (such as SUPG and GLS) and commercial packages provide physically meaningful values for the concentration of the chemical species for various realistic benchmark problems. Furthermore, we also quantify the errors incurred in satisfying the local and global species balance for two popular chemical kinetics schemes: CDIMA (chlorine dioxide-iodine-malonic acid) and BZ (Belousov--Zhabotinsky). Based on these numerical simulations, we show that SUPG and GLS produce unphysical values for concentration of chemical species due to the violation of the non-negative constraint, contain spurious node-to-node oscillations, and have large errors in local and global species balance. Second, we proposed a novel finite element formulation to overcome the above difficulties. The proposed locally conservative non-negative computational framework based on low-order least-squares finite elements is able to preserve these underlying physical and mathematical properties

The finite element response matrix method

International Nuclear Information System (INIS)

Nakata, H.; Martin, W.R.

1983-02-01

A new technique is developed with an alternative formulation of the response matrix method implemented with the finite element scheme. Two types of response matrices are generated from the Galerkin solution to the weak form of the diffusion equation subject to an arbitrary current and source. The piecewise polynomials are defined in two levels, the first for the local (assembly) calculations and the second for the global (core) response matrix calculations. This finite element response matrix technique was tested in two 2-dimensional test problems, 2D-IAEA benchmark problem and Biblis benchmark problem, with satisfatory results. The computational time, whereas the current code is not extensively optimized, is of the same order of the well estabilished coarse mesh codes. Furthermore, the application of the finite element technique in an alternative formulation of response matrix method permits the method to easily incorporate additional capabilities such as treatment of spatially dependent cross-sections, arbitrary geometrical configurations, and high heterogeneous assemblies. (Author) [pt

Finite Element Method in Machining Processes

CERN Document Server

Markopoulos, Angelos P

2013-01-01

Finite Element Method in Machining Processes provides a concise study on the way the Finite Element Method (FEM) is used in the case of manufacturing processes, primarily in machining. The basics of this kind of modeling are detailed to create a reference that will provide guidelines for those who start to study this method now, but also for scientists already involved in FEM and want to expand their research. A discussion on FEM, formulations and techniques currently in use is followed up by machining case studies. Orthogonal cutting, oblique cutting, 3D simulations for turning and milling, grinding, and state-of-the-art topics such as high speed machining and micromachining are explained with relevant examples. This is all supported by a literature review and a reference list for further study. As FEM is a key method for researchers in the manufacturing and especially in the machining sector, Finite Element Method in Machining Processes is a key reference for students studying manufacturing processes but al...

Inside finite elements

CERN Document Server

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.

Numerical simulations of granular dynamics: I. Hard-sphere discrete element method and tests

Science.gov (United States)

Richardson, Derek C.; Walsh, Kevin J.; Murdoch, Naomi; Michel, Patrick

2011-03-01

We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body tree code pkdgrav to search for collisions and compute particle trajectories. Collisions are treated as instantaneous point-contact events between rigid spheres. Particle confinement is achieved by combining arbitrary combinations of four provided wall primitives, namely infinite plane, finite disk, infinite cylinder, and finite cylinder, and degenerate cases of these. Various wall movements, including translation, oscillation, and rotation, are supported. We provide full derivations of collision prediction and resolution equations for all geometries and motions. Several tests of the method are described, including a model granular “atmosphere” that achieves correct energy equipartition, and a series of tumbler simulations that show the expected transition from tumbling to centrifuging as a function of rotation rate.

Finite Element Modelling of Seismic Liquefaction in Soils

NARCIS (Netherlands)

Galavi, V.; Petalas, A.; Brinkgreve, R.B.J.

2013-01-01

Numerical aspects of seismic liquefaction in soils as implemented in the finite element code, PLAXIS, is described in this paper. After description of finite element equations of dynamic problems, three practical dynamic boundary conditions, namely viscous boundary tractions, tied degrees of freedom

A novel method for computation of the discrete Fourier transform over characteristic two finite field of even extension degree

OpenAIRE

Fedorenko, Sergei V.

2011-01-01

A novel method for computation of the discrete Fourier transform over a finite field with reduced multiplicative complexity is described. If the number of multiplications is to be minimized, then the novel method for the finite field of even extension degree is the best known method of the discrete Fourier transform computation. A constructive method of constructing for a cyclic convolution over a finite field is introduced.

Solving ill-posed control problems by stabilized finite element methods: an alternative to Tikhonov regularization

Science.gov (United States)

Burman, Erik; Hansbo, Peter; Larson, Mats G.

2018-03-01

Tikhonov regularization is one of the most commonly used methods for the regularization of ill-posed problems. In the setting of finite element solutions of elliptic partial differential control problems, Tikhonov regularization amounts to adding suitably weighted least squares terms of the control variable, or derivatives thereof, to the Lagrangian determining the optimality system. In this note we show that the stabilization methods for discretely ill-posed problems developed in the setting of convection-dominated convection-diffusion problems, can be highly suitable for stabilizing optimal control problems, and that Tikhonov regularization will lead to less accurate discrete solutions. We consider some inverse problems for Poisson’s equation as an illustration and derive new error estimates both for the reconstruction of the solution from the measured data and reconstruction of the source term from the measured data. These estimates include both the effect of the discretization error and error in the measurements.

Guided wave mode selection for inhomogeneous elastic waveguides using frequency domain finite element approach.

Science.gov (United States)

Chillara, Vamshi Krishna; Ren, Baiyang; Lissenden, Cliff J

2016-04-01

This article describes the use of the frequency domain finite element (FDFE) technique for guided wave mode selection in inhomogeneous waveguides. Problems with Rayleigh-Lamb and Shear-Horizontal mode excitation in isotropic homogeneous plates are first studied to demonstrate the application of the approach. Then, two specific cases of inhomogeneous waveguides are studied using FDFE. Finally, an example of guided wave mode selection for inspecting disbonds in composites is presented. Identification of sensitive and insensitive modes for defect inspection is demonstrated. As the discretization parameters affect the accuracy of the results obtained from FDFE, effect of spatial discretization and the length of the domain used for the spatial fast Fourier transform are studied. Some recommendations with regard to the choice of the above parameters are provided. Copyright © 2015 Elsevier B.V. All rights reserved.

Solution of the neutron diffusion equation at two groups of energy by method of triangular finite elements

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

Mesh-size errors in diffusion-theory calculations using finite-difference and finite-element methods

International Nuclear Information System (INIS)

Baker, A.R.

1982-07-01

A study has been performed of mesh-size errors in diffusion-theory calculations using finite-difference and finite-element methods. As the objective was to illuminate the issues, the study was performed for a 1D slab model of a reactor with one neutron-energy group for which analytical solutions were possible. A computer code SLAB was specially written to perform the finite-difference and finite-element calculations and also to obtain the analytical solutions. The standard finite-difference equations were obtained by starting with an expansion of the neutron current in powers of the mesh size, h, and keeping terms as far as h 2 . It was confirmed that these equations led to the well-known result that the criticality parameter varied with the square of the mesh size. An improved form of the finite-difference equations was obtained by continuing the expansion for the neutron current as far as the term in h 4 . In this case, the critical parameter varied as the fourth power of the mesh size. The finite-element solutions for 2 and 3 nodes per element revealed that the criticality parameter varied as the square and fourth power of the mesh size, respectively. Numerical results are presented for a bare reactive core of uniform composition with 2 zones of different uniform mesh and for a reactive core with an absorptive reflector. (author)

A geometric toolbox for tetrahedral finite element partitions

NARCIS (Netherlands)

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.

Coarse mesh finite element method for boiling water reactor physics analysis

International Nuclear Information System (INIS)

Ellison, P.G.

1983-01-01

A coarse mesh method is formulated for the solution of Boiling Water Reactor physics problems using two group diffusion theory. No fuel assembly cross-section homogenization is required; water gaps, control blades and fuel pins of varying enrichments are treated explicitly. The method combines constrained finite element discretization with infinite lattice super cell trial functions to obtain coarse mesh solutions for which the only approximations are along the boundaries between fuel assemblies. The method is applied to bench mark Boiling Water Reactor problems to obtain both the eigenvalue and detailed flux distributions. The solutions to these problems indicate the method is useful in predicting detailed power distributions and eigenvalues for Boiling Water Reactor physics problems

Nonlinear finite element analysis of nuclear reinforced prestressed concrete containments up to ultimate load capacity

International Nuclear Information System (INIS)

Gupta, A.; Singh, R.K.; Kushwaha, H.S.; Mahajan, S.C.; Kakodkar, A.

1996-01-01

For safety evaluation of nuclear structures a finite element code ULCA (Ultimate Load Capacity Assessment) has been developed. Eight/nine noded isoparametric quadrilateral plate/shell element with reinforcement as a through thickness discrete but integral smeared layer of the element is presented to analyze reinforced and prestressed concrete structures. Various constitutive models such as crushing, cracking in tension, tension stiffening and rebar yielding are studied and effect of these parameters on the reserve strength of structures is brought out through a number of benchmark tests. A global model is used to analyze the prestressed concrete containment wall of a typical 220 MWe Pressurized Heavy Water Reactor (PHWR) up to its ultimate capacity. This demonstrates the adequacy of Indian PHWR containment design to withstand severe accident loads

Mixed finite element simulations in two-dimensional groundwater flow problems

International Nuclear Information System (INIS)

Kimura, Hideo

1989-01-01

A computer code of groundwater flow in two-dimensional porous media based on the mixed finite element method was developed for accurate approximations of Darcy velocities in safety evaluation of radioactive waste disposal. The mixed finite element procedure solves for both the Darcy velocities and pressure heads simultaneously in the Darcy equation and continuity equation. Numerical results of a single well pumping at a constant rate in a uniform flow field showed that the mixed finite element method gives more accurate Darcy velocities nearly 50 % on average error than standard finite element method. (author)

New formulation of the discrete element method

Science.gov (United States)

Rojek, Jerzy; Zubelewicz, Aleksander; Madan, Nikhil; Nosewicz, Szymon

2018-01-01

A new original formulation of the discrete element method based on the soft contact approach is presented in this work. The standard DEM has heen enhanced by the introduction of the additional (global) deformation mode caused by the stresses in the particles induced by the contact forces. Uniform stresses and strains are assumed for each particle. The stresses are calculated from the contact forces. The strains are obtained using an inverse constitutive relationship. The strains allow us to obtain deformed particle shapes. The deformed shapes (ellipses) are taken into account in contact detection and evaluation of the contact forces. A simple example of a uniaxial compression of a rectangular specimen, discreti.zed with equal sized particles is simulated to verify the DDEM algorithm. The numerical example shows that a particle deformation changes the particle interaction and the distribution of forces in the discrete element assembly. A quantitative study of micro-macro elastic properties proves the enhanced capabilities of the DDEM as compared to standard DEM.

Finite element analysis of rotating beams physics based interpolation

CERN Document Server

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.

Discontinuous finite element solution of the radiation diffusion equation on arbitrary polygonal meshes and locally adapted quadrilateral grids

International Nuclear Information System (INIS)

Ragusa, Jean C.

2015-01-01

In this paper, we propose a piece-wise linear discontinuous (PWLD) finite element discretization of the diffusion equation for arbitrary polygonal meshes. It is based on the standard diffusion form and uses the symmetric interior penalty technique, which yields a symmetric positive definite linear system matrix. A preconditioned conjugate gradient algorithm is employed to solve the linear system. Piece-wise linear approximations also allow a straightforward implementation of local mesh adaptation by allowing unrefined cells to be interpreted as polygons with an increased number of vertices. Several test cases, taken from the literature on the discretization of the radiation diffusion equation, are presented: random, sinusoidal, Shestakov, and Z meshes are used. The last numerical example demonstrates the application of the PWLD discretization to adaptive mesh refinement

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

A two-dimensional discontinuous heterogeneous finite element method for neutron transport calculations

International Nuclear Information System (INIS)

Masiello, E.; Sanchez, R.

2007-01-01

A discontinuous heterogeneous finite element method is presented and discussed. The method is intended for realistic numerical pin-by-pin lattice calculations when an exact representation of the geometric shape of the pins is made without need for homogenization. The method keeps the advantages of conventional discrete ordinate methods, such as fast execution together with the possibility to deal with a large number of spatial meshes, while minimizing the need for geometric modeling. It also provides a complete factorization in space, angle, and energy for the discretized matrices and minimizes, thus, storage requirements. An angular multigrid acceleration technique has also been developed to speed up the rate of convergence of the inner iterations. A particular aspect of this acceleration is the introduction of boundary restriction and prolongation operators that minimize oscillatory behavior and enhance positivity. Numerical tests are presented that show the high precision of the method and the efficiency of the angular multigrid acceleration. (authors)

Implementation of thermo-viscoplastic constitutive equations into the finite element code ABAQUS

International Nuclear Information System (INIS)

Youn, Sam Son; Lee, Soon Bok; Kim, Jong Bum; Lee, Hyeong Yeon; Yoo, Bong

1998-01-01

Sophisticated viscoplatic constitutive laws describing material behavior at high temperature have been implemented in the general-purpose finite element code ABAQUS to predict the viscoplastic response of structures to cyclic loading. Because of the complexity of viscoplastic constitutive equation, the general implementation methods are developed. The solution of the non-linear system of algebraic equations arising from time discretization is determined using line-search and back-tracking in combination with Newton method. The time integration method of the constitutive equations is based on semi-implicit method with efficient time step control. For numerical examples, the viscoplastic model proposed by Chaboche is implemented and several applications are illustrated

Face-based smoothed finite element method for real-time simulation of soft tissue

Science.gov (United States)

Mendizabal, Andrea; Bessard Duparc, Rémi; Bui, Huu Phuoc; Paulus, Christoph J.; Peterlik, Igor; Cotin, Stéphane

2017-03-01

In soft tissue surgery, a tumor and other anatomical structures are usually located using the preoperative CT or MR images. However, due to the deformation of the concerned tissues, this information suffers from inaccuracy when employed directly during the surgery. In order to account for these deformations in the planning process, the use of a bio-mechanical model of the tissues is needed. Such models are often designed using the finite element method (FEM), which is, however, computationally expensive, in particular when a high accuracy of the simulation is required. In our work, we propose to use a smoothed finite element method (S-FEM) in the context of modeling of the soft tissue deformation. This numerical technique has been introduced recently to overcome the overly stiff behavior of the standard FEM and to improve the solution accuracy and the convergence rate in solid mechanics problems. In this paper, a face-based smoothed finite element method (FS-FEM) using 4-node tetrahedral elements is presented. We show that in some cases, the method allows for reducing the number of degrees of freedom, while preserving the accuracy of the discretization. The method is evaluated on a simulation of a cantilever beam loaded at the free end and on a simulation of a 3D cube under traction and compression forces. Further, it is applied to the simulation of the brain shift and of the kidney's deformation. The results demonstrate that the method outperforms the standard FEM in a bending scenario and that has similar accuracy as the standard FEM in the simulations of the brain-shift and of the kidney's deformation.

Ultimate limit state design of sheet pile walls by finite elements and nonlinear programming

DEFF Research Database (Denmark)

Krabbenhøft, Kristian; Damkilde, Lars; Krabbenhøft, Sven

2005-01-01

The design of sheet pile walls by lower bound limit analysis is considered. The design problem involves the determination of the necessary yield moment of the wall, the wall depth and the anchor force such that the structure is able to sustain the given loads. This problem is formulated...... as a nonlinear programming problem where the yield moment of the wall is minimized subject to equilibrium and yield conditions. The finite element discretization used enables exact fulfillment of these conditions and thus, according to the lower bound theorem, the solutions are safe....

An efficient discontinuous Galerkin finite element method for highly accurate solution of maxwell equations

KAUST Repository

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.

An efficient discontinuous Galerkin finite element method for highly accurate solution of maxwell equations

KAUST Repository

Liu, Meilin; Sirenko, Kostyantyn; Bagci, Hakan

2012-01-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.

Finite size effects of a pion matrix element

International Nuclear Information System (INIS)

Guagnelli, M.; Jansen, K.; Palombi, F.; Petronzio, R.; Shindler, A.; Wetzorke, I.

2004-01-01

We investigate finite size effects of the pion matrix element of the non-singlet, twist-2 operator corresponding to the average momentum of non-singlet quark densities. Using the quenched approximation, they come out to be surprisingly large when compared to the finite size effects of the pion mass. As a consequence, simulations of corresponding nucleon matrix elements could be affected by finite size effects even stronger which could lead to serious systematic uncertainties in their evaluation

Verification of Orthogrid Finite Element Modeling Techniques

Science.gov (United States)

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.

Review of Tomographic Imaging using Finite Element Method

Directory of Open Access Journals (Sweden)

Mohd Fua’ad RAHMAT

2011-12-01

Full Text Available Many types of techniques for process tomography were proposed and developed during the past 20 years. This paper review the techniques and the current state of knowledge and experience on the subject, aimed at highlighting the problems associated with the non finite element methods, such as the ill posed, ill conditioned which relates to the accuracy and sensitivity of measurements. In this paper, considerations for choice of sensors and its applications were outlined and descriptions of non finite element tomography systems were presented. The finite element method tomography system as obtained from recent works, suitable for process control and measurement were also presented.

Magnetic materials and 3D finite element modeling

CERN Document Server

Bastos, Joao Pedro A

2014-01-01

Magnetic Materials and 3D Finite Element Modeling explores material characterization and finite element modeling (FEM) applications. This book relates to electromagnetic analysis based on Maxwell’s equations and application of the finite element (FE) method to low frequency devices. A great source for senior undergraduate and graduate students in electromagnetics, it also supports industry professionals working in magnetics, electromagnetics, ferromagnetic materials science and electrical engineering. The authors present current concepts on ferromagnetic material characterizations and losses. They provide introductory material; highlight basic electromagnetics, present experimental and numerical modeling related to losses and focus on FEM applied to 3D applications. They also explain various formulations, and discuss numerical codes.

New mixed finite-element methods

International Nuclear Information System (INIS)

Franca, L.P.

1987-01-01

New finite-element methods are proposed for mixed variational formulations. The methods are constructed by adding to the classical Galerkin method various least-squares like terms. The additional terms involve integrals over element interiors, and include mesh-parameter dependent coefficients. The methods are designed to enhance stability. Consistency is achieved in the sense that exact solutions identically satisfy the variational equations.Applied to several problems, simple finite-element interpolations are rendered convergent, including convenient equal-order interpolations generally unstable within the Galerkin approach. The methods are subdivided into two classes according to the manner in which stability is attained: (1) circumventing Babuska-Brezzi condition methods; (2) satisfying Babuska-Brezzi condition methods. Convergence is established for each class of methods. Applications of the first class of methods to Stokes flow and compressible linear elasticity are presented. The second class of methods is applied to the Poisson, Timoshenko beam and incompressible elasticity problems. Numerical results demonstrate the good stability and accuracy of the methods, and confirm the error estimates

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

Crack nucleation in solid materials under external load - simulations with the Discrete Element Method

Directory of Open Access Journals (Sweden)

Klejment Piotr

2018-01-01

Full Text Available Numerical analysis of cracking processes require an appropriate numerical technique. Classical engineering approach to the problem has its roots in the continuum mechanics and is based mainly on the Finite Element Method. This technique allows simulations of both elastic and large deformation processes, so it is very popular in the engineering applications. However, a final effect of cracking - fragmentation of an object at hand can hardly be described by this approach in a numerically efficient way since it requires a solution of a problem of nontrivial evolving in time boundary conditions. We focused our attention on the Discrete Element Method (DEM, which by definition implies “molecular” construction of the matter. The basic idea behind DEM is to represent an investigated body as an assemblage of discrete particles interacting with each other. Breaking interaction bonds between particles induced by external forces imeditelly implies creation/evolution of boundary conditions. In this study we used the DEM approach to simulate cracking process in the three dimensional solid material under external tension. The used numerical model, although higly simplified, can be used to describe behaviour of such materials like thin films, biological tissues, metal coatings, to name a few.

On the development of a three-dimensional finite-element groundwater flow model of the saturated zone, Yucca Mountain, Nevada

International Nuclear Information System (INIS)

Czarnecki, J.B.; Faunt, C.C.; Gable, C.W.; Zyvoloski, G.A.

1996-01-01

Development of a preliminary three-dimensional model of the saturated zone at Yucca Mountain, the potential location for a high-level nuclear waste repository, is presented. The development of the model advances the technology of interfacing: (1)complex three-dimensional hydrogeologic framework modeling; (2) fully three-dimensional, unstructured, finite-element mesh generation; and (3) groundwater flow, heat, and transport simulation. The three-dimensional hydrogeologic framework model is developed using maps, cross sections, and well data. The framework model data are used to feed an automated mesh generator, designed to discretize irregular three-dimensional solids,a nd to assign materials properties from the hydrogeologic framework model to the tetrahedral elements. The mesh generator facilitated the addition of nodes to the finite-element mesh which correspond to the exact three-dimensional position of the potentiometric surface based on water-levels from wells. A ground water flow and heat simulator is run with the resulting finite- element mesh, within a parameter-estimation program. The application of the parameter-estimation program is designed to provide optimal values of permeability and specified fluxes over the model domain to minimize the residual between observed and simulated water levels

Validation of High Displacement Piezoelectric Actuator Finite Element Models

Science.gov (United States)

Taleghani, B. K.

2000-01-01

The paper presents the results obtained by using NASTRAN(Registered Trademark) and ANSYS(Regitered Trademark) finite element codes to predict doming of the THUNDER piezoelectric actuators during the manufacturing process and subsequent straining due to an applied input voltage. To effectively use such devices in engineering applications, modeling and characterization are essential. Length, width, dome height, and thickness are important parameters for users of such devices. Therefore, finite element models were used to assess the effects of these parameters. NASTRAN(Registered Trademark) and ANSYS(Registered Trademark) used different methods for modeling piezoelectric effects. In NASTRAN(Registered Trademark), a thermal analogy was used to represent voltage at nodes as equivalent temperatures, while ANSYS(Registered Trademark) processed the voltage directly using piezoelectric finite elements. The results of finite element models were validated by using the experimental results.

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)

A moving finite element model of the high density z-pinch

International Nuclear Information System (INIS)

Glasser, A.H.

1989-01-01

This paper presents the results of computations of the behavior of the fiber-initiated high density Z-pinch (HDZP). It purpose is twofold. One is to study the behavior of the physical system itself as an interesting controlled fusion experiment. The main result of this study is a demonstration of the relaxation of the full inertial behavior of the pinch to simplified self-similar behavior in which the forces on the system are in near balance. The moving free boundary and violent initial behavior of this configuration require carful treatment. This leads to the other purpose of the work, to use this realistic physical system as a test-bed for a general-purpose 1-dimensional code based on moving finite elements. A key step in accomplishing this goal has been the recognition that numerical stability of the discretized equations has required the use of nonconservative quantities as the fundamental dependent variables to be discretized. The main result of this work is a code which is capable of treating a very general class of nonlinear, time-dependent fluid equations. copyright 1989 Academic Press, Inc

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

The finite element method entails several approximations. Hence it ... researchers have designed several pathological tests to validate any new finite element. The .... Three dimensional thick shell elements using a hybrid/mixed formu- lation.

Crack Propagation by Finite Element Method

Directory of Open Access Journals (Sweden)

Luiz Carlos H. Ricardo

2018-01-01

Full Text Available Crack propagation simulation began with the development of the finite element method; the analyses were conducted to obtain a basic understanding of the crack growth. Today structural and materials engineers develop structures and materials properties using this technique. The aim of this paper is to verify the effect of different crack propagation rates in determination of crack opening and closing stress of an ASTM specimen under a standard suspension spectrum loading from FDandE SAE Keyhole Specimen Test Load Histories by finite element analysis. To understand the crack propagation processes under variable amplitude loading, retardation effects are observed

Finite element analysis of plastic recycling machine designed for ...

African Journals Online (AJOL)

... design was evaluated using finite element analysis (FEA) tool in Solid Works Computer ... Also, a minimum factor of safety value of 5.3 was obtained for shredder shaft ... Machine; Design; Recycling; Sustainability; Finite Element; Simulation ...

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)

An introduction to the UNCLE finite element scheme

Energy Technology Data Exchange (ETDEWEB)

Enderby, J A [UK Atomic Energy Authority, Northern Division, Risley Nuclear Power Development Establishment, Risley, Warrington (United Kingdom)

1983-05-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)

A particle finite element method for machining simulations

Science.gov (United States)

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.

A finite element primer for beginners the basics

CERN Document Server

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

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)

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.

Automation of finite element methods

CERN Document Server

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.

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

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.

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, Eric M.

1997-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1996-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed

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....... This paper presents a novel mechanical integrated sleeve wedge anchorage which seem very promising when perusing the scope of ultimate utilization of CFRP 8mm rods (with a tension capacity of approximately 140kN). Compression transverse to the CFRP is evaluated to prevent premature failure. The anchorage...

Numerical Integration Techniques for Curved-Element Discretizations of Molecule–Solvent Interfaces

Science.gov (United States)

Bardhan, Jaydeep P.; Altman, Michael D.; Willis, David J.; Lippow, Shaun M.; Tidor, Bruce; White, Jacob K.

2012-01-01

Surface formulations of biophysical modeling problems offer attractive theoretical and computational properties. Numerical simulations based on these formulations usually begin with discretization of the surface under consideration; often, the surface is curved, possessing complicated structure and possibly singularities. Numerical simulations commonly are based on approximate, rather than exact, discretizations of these surfaces. To assess the strength of the dependence of simulation accuracy on the fidelity of surface representation, we have developed methods to model several important surface formulations using exact surface discretizations. Following and refining Zauhar’s work (J. Comp.-Aid. Mol. Des. 9:149-159, 1995), we define two classes of curved elements that can exactly discretize the van der Waals, solvent-accessible, and solvent-excluded (molecular) surfaces. We then present numerical integration techniques that can accurately evaluate nonsingular and singular integrals over these curved surfaces. After validating the exactness of the surface discretizations and demonstrating the correctness of the presented integration methods, we present a set of calculations that compare the accuracy of approximate, planar-triangle-based discretizations and exact, curved-element-based simulations of surface-generalized-Born (sGB), surface-continuum van der Waals (scvdW), and boundary-element method (BEM) electrostatics problems. Results demonstrate that continuum electrostatic calculations with BEM using curved elements, piecewise-constant basis functions, and centroid collocation are nearly ten times more accurate than planartriangle BEM for basis sets of comparable size. The sGB and scvdW calculations give exceptional accuracy even for the coarsest obtainable discretized surfaces. The extra accuracy is attributed to the exact representation of the solute–solvent interface; in contrast, commonly used planar-triangle discretizations can only offer improved

Accurate and efficient implementation of the von Neumann representation for laser pulses with discrete and finite spectra

International Nuclear Information System (INIS)

Dimler, Frank; Fechner, Susanne; Rodenberg, Alexander; Brixner, Tobias; Tannor, David J

2009-01-01

We recently introduced the von Neumann picture, a joint time-frequency representation, for describing ultrashort laser pulses. The method exploits a discrete phase-space lattice of nonorthogonal Gaussians to represent the pulses; an arbitrary pulse shape can be represented on this lattice in a one-to-one manner. Although the representation was originally defined for signals with an infinite continuous spectrum, it can be adapted to signals with discrete and finite spectrum with great computational savings, provided that discretization and truncation effects are handled with care. In this paper, we present three methods that avoid loss of accuracy due to these effects. The approach has immediate application to the representation and manipulation of femtosecond laser pulses produced by a liquid-crystal mask with a discrete and finite number of pixels.

Hydrothermal analysis in engineering using control volume finite element method

CERN Document Server

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),

A comparison of the lattice discrete particle method to the finite-element method and the K&C material model for simulating the static and dynamic response of concrete.

Energy Technology Data Exchange (ETDEWEB)

Smith, Jovanca J.; Bishop, Joseph E.

2013-11-01

This report summarizes the work performed by the graduate student Jovanca Smith during a summer internship in the summer of 2012 with the aid of mentor Joe Bishop. The projects were a two-part endeavor that focused on the use of the numerical model called the Lattice Discrete Particle Model (LDPM). The LDPM is a discrete meso-scale model currently used at Northwestern University and the ERDC to model the heterogeneous quasi-brittle material, concrete. In the first part of the project, LDPM was compared to the Karagozian and Case Concrete Model (K&C) used in Presto, an explicit dynamics finite-element code, developed at Sandia National Laboratories. In order to make this comparison, a series of quasi-static numerical experiments were performed, namely unconfined uniaxial compression tests on four varied cube specimen sizes, three-point bending notched experiments on three proportional specimen sizes, and six triaxial compression tests on a cylindrical specimen. The second part of this project focused on the application of LDPM to simulate projectile perforation on an ultra high performance concrete called CORTUF. This application illustrates the strengths of LDPM over traditional continuum models.

Parallel iterative procedures for approximate solutions of wave propagation by finite element and finite difference methods

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.

Soft tissue deformation estimation by spatio-temporal Kalman filter finite element method.

Science.gov (United States)

Yarahmadian, Mehran; Zhong, Yongmin; Gu, Chengfan; Shin, Jaehyun

2018-01-01

Soft tissue modeling plays an important role in the development of surgical training simulators as well as in robot-assisted minimally invasive surgeries. It has been known that while the traditional Finite Element Method (FEM) promises the accurate modeling of soft tissue deformation, it still suffers from a slow computational process. This paper presents a Kalman filter finite element method to model soft tissue deformation in real time without sacrificing the traditional FEM accuracy. The proposed method employs the FEM equilibrium equation and formulates it as a filtering process to estimate soft tissue behavior using real-time measurement data. The model is temporally discretized using the Newmark method and further formulated as the system state equation. Simulation results demonstrate that the computational time of KF-FEM is approximately 10 times shorter than the traditional FEM and it is still as accurate as the traditional FEM. The normalized root-mean-square error of the proposed KF-FEM in reference to the traditional FEM is computed as 0.0116. It is concluded that the proposed method significantly improves the computational performance of the traditional FEM without sacrificing FEM accuracy. The proposed method also filters noises involved in system state and measurement data.

High-order solution methods for grey discrete ordinates thermal radiative transfer

Energy Technology Data Exchange (ETDEWEB)

Maginot, Peter G., E-mail: maginot1@llnl.gov [Lawrence Livermore National Laboratory, Livermore, CA 94551 (United States); Ragusa, Jean C., E-mail: jean.ragusa@tamu.edu [Department of Nuclear Engineering, Texas A& M University, College Station, TX 77843 (United States); Morel, Jim E., E-mail: morel@tamu.edu [Department of Nuclear Engineering, Texas A& M University, College Station, TX 77843 (United States)

2016-12-15

This work presents a solution methodology for solving the grey radiative transfer equations that is both spatially and temporally more accurate than the canonical radiative transfer solution technique of linear discontinuous finite element discretization in space with implicit Euler integration in time. We solve the grey radiative transfer equations by fully converging the nonlinear temperature dependence of the material specific heat, material opacities, and Planck function. The grey radiative transfer equations are discretized in space using arbitrary-order self-lumping discontinuous finite elements and integrated in time with arbitrary-order diagonally implicit Runge–Kutta time integration techniques. Iterative convergence of the radiation equation is accelerated using a modified interior penalty diffusion operator to precondition the full discrete ordinates transport operator.

Analysis of 2D reactor core using linear perturbation theory and nodal finite element methods

International Nuclear Information System (INIS)

Adrian Mugica; Edmundo del Valle

2005-01-01

In this work the multigroup steady state neutron diffusion equations are solved using the nodal finite element method (NFEM) and the Linear Perturbation Theory (LPT) for XY geometry. The NFEM used corresponds to the Raviart-Thomas schemes RT0 and RT1, interpolating 5 and 12 parameters respectively in each node of the space discretization. The accuracy of these methods is related with the dimension of the space approximation and the mesh size. Therefore, using fine meshes and the RT0 or RT1 nodal methods leads to a large an interesting eigenvalue problem. The finite element method used to discretize the weak formulation of the diffusion equations is the Galerkin one. The algebraic structure of the discrete eigenvalue problem is obtained and solved using the Wielandt technique and the BGSTAB iterative method using the SPARSKIT package developed by Yousef Saad. The results obtained with LPT show good agreement with the results obtained directly for the perturbed problem. In fact, the cpu time to solve a single problem, the unperturbed and the perturbed one, is practically the same but when one is focused in shuffling many times two different assemblies in the core then the LPT technique becomes quite useful to get good approximations in a short time. This particular problem was solved for one quarter-core with NFEM. Thus, the computer program based on LPT can be used to perform like an analysis tool in the fuel reload optimization or combinatory analysis to get reload patterns in nuclear power plants once that it had been incorporated with the thermohydraulic aspects needed to simulate accurately a real problem. The maximum differences between the NFEM and LPT for the three LWR reactor cores are about 250 pcm. This quantity is considered an acceptable value for this kind of analysis. (authors)

Generalized multiscale finite element methods for problems in perforated heterogeneous domains

KAUST Repository

Chung, Eric T.

2015-06-08

Complex processes in perforated domains occur in many real-world applications. These problems are typically characterized by physical processes in domains with multiple scales. Moreover, these problems are intrinsically multiscale and their discretizations can yield very large linear or nonlinear systems. In this paper, we investigate multiscale approaches that attempt to solve such problems on a coarse grid by constructing multiscale basis functions in each coarse grid, where the coarse grid can contain many perforations. In particular, we are interested in cases when there is no scale separation and the perforations can have different sizes. In this regard, we mention some earlier pioneering works, where the authors develop multiscale finite element methods. In our paper, we follow Generalized Multiscale Finite Element Method (GMsFEM) and develop a multiscale procedure where we identify multiscale basis functions in each coarse block using snapshot space and local spectral problems. We show that with a few basis functions in each coarse block, one can approximate the solution, where each coarse block can contain many small inclusions. We apply our general concept to (1) Laplace equation in perforated domains; (2) elasticity equation in perforated domains; and (3) Stokes equations in perforated domains. Numerical results are presented for these problems using two types of heterogeneous perforated domains. The analysis of the proposed methods will be presented elsewhere. © 2015 Taylor & Francis

Global-Local Finite Element Analysis of Bonded Single-Lap Joints

Science.gov (United States)

Kilic, Bahattin; Madenci, Erdogan; Ambur, Damodar R.

2004-01-01

Adhesively bonded lap joints involve dissimilar material junctions and sharp changes in geometry, possibly leading to premature failure. Although the finite element method is well suited to model the bonded lap joints, traditional finite elements are incapable of correctly resolving the stress state at junctions of dissimilar materials because of the unbounded nature of the stresses. In order to facilitate the use of bonded lap joints in future structures, this study presents a finite element technique utilizing a global (special) element coupled with traditional elements. The global element includes the singular behavior at the junction of dissimilar materials with or without traction-free surfaces.

In-plane material continuity for the discrete material optimization method

DEFF Research Database (Denmark)

Sørensen, Rene; Lund, Erik

2015-01-01

When performing discrete material optimization of laminated composite structures, the variation of the in-plane material continuity is typically governed by the size of the finite element discretization. For a fine mesh, this can lead to designs that cannot be manufactured due to the complexity...

Limit sets for the discrete spectrum of complex Jacobi matrices

International Nuclear Information System (INIS)

Golinskii, L B; Egorova, I E

2005-01-01

The discrete spectrum of complex Jacobi matrices that are compact perturbations of the discrete Laplacian is studied. The precise stabilization rate (in the sense of order) of the matrix elements ensuring the finiteness of the discrete spectrum is found. An example of a Jacobi matrix with discrete spectrum having a unique limit point is constructed. These results are discrete analogues of Pavlov's well-known results on Schroedinger operators with complex potential on a half-axis.

Finite element simulation of laser transmission welding of dissimilar ...

African Journals Online (AJOL)

user

materials between polyvinylidene fluoride and titanium ... finite element (FE) thermal model is developed to simulate the laser ... Keywords: Laser transmission welding, Temperature field, Weld dimension, Finite element analysis, Thermal modeling. 1. .... 4) The heating phenomena due to the phase changes are neglected.

Finite element analysis of tibial fractures

DEFF Research Database (Denmark)

Wong, Christian Nai En; Mikkelsen, Mikkel Peter W; Hansen, Leif Berner

2010-01-01

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

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

Efendiev, Yalchin R.; Galvis, Juan; Lazarov, Raytcho D.; Moon, M.; Sarkis, Marcus V.

2013-01-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.

Generalized multiscale finite element method. Symmetric interior penalty coupling

KAUST Repository

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 Volumes Discretization of Topology Optimization Problems

DEFF Research Database (Denmark)

Evgrafov, Anton; Gregersen, Misha Marie; Sørensen, Mads Peter

, FVMs represent a standard method of discretization within engineering communities dealing with computational uid dy- namics, transport, and convection-reaction problems. Among various avours of FVMs, cell based approaches, where all variables are associated only with cell centers, are particularly...... computations is done using nite element methods (FEMs). Despite some limited recent eorts [1, 2], we have only started to develop our understanding of the interplay between the control in the coecients and FVMs. Recent advances in discrete functional analysis allow us to analyze convergence of FVM...... of the induced parametrization of the design space that allows optimization algorithms to eciently explore it, and the ease of integration with existing computational codes in a variety of application areas, the simplicity and eciency of sensitivity analyses|all stemming from the use of the same grid throughout...

Analysis of anisotropic crack problems using coupled meshless and fractal finite element method

International Nuclear Information System (INIS)

Rao, B N; Rajesh, K N

2010-01-01

This paper presents a coupling technique for integrating the element-free Galerkin method (EFGM) with fractal the finite element method (FFEM) for analyzing homogeneous, anisotropic, and two dimensional linear-elastic cracked structures subjected to mixed-mode (modes I and II) loading conditions. FFEM is adopted for discretization of domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The proposed method combines the best features of EFGM and FFEM, in the sense that no structured mesh or special enriched basis functions are necessary and no post-processing (employing any path independent integrals) is needed to determine fracture parameters such as stress intensity factors (SIFs) and T-stress. The numerical results based on all four orthotropic cases show that SIFs and T-stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study.

Multiscale Modeling of Blood Flow: Coupling Finite Elements with Smoothed Dissipative Particle Dynamics

KAUST Repository

Moreno Chaparro, Nicolas; Vignal, Philippe; Li, Jun; Calo, Victor M.

2013-01-01

A variational multi scale approach to model blood flow through arteries is proposed. A finite element discretization to represent the coarse scales (macro size), is coupled to smoothed dissipative particle dynamics that captures the fine scale features (micro scale). Blood is assumed to be incompressible, and flow is described through the Navier Stokes equation. The proposed cou- pling is tested with two benchmark problems, in fully coupled systems. Further refinements of the model can be incorporated in order to explicitly include blood constituents and non-Newtonian behavior. The suggested algorithm can be used with any particle-based method able to solve the Navier-Stokes equation.

Multiscale Modeling of Blood Flow: Coupling Finite Elements with Smoothed Dissipative Particle Dynamics

KAUST Repository

Moreno Chaparro, Nicolas

2013-06-01

A variational multi scale approach to model blood flow through arteries is proposed. A finite element discretization to represent the coarse scales (macro size), is coupled to smoothed dissipative particle dynamics that captures the fine scale features (micro scale). Blood is assumed to be incompressible, and flow is described through the Navier Stokes equation. The proposed cou- pling is tested with two benchmark problems, in fully coupled systems. Further refinements of the model can be incorporated in order to explicitly include blood constituents and non-Newtonian behavior. The suggested algorithm can be used with any particle-based method able to solve the Navier-Stokes equation.

Finite Element Analysis of Circular Plate using SolidWorks

International Nuclear Information System (INIS)

Kang, Yeo Jin; Jhung, Myung Jo

2011-01-01

Circular plates are used extensively in mechanical engineering for nuclear reactor internal components. The examples in the reactor vessel internals are upper guide structure support plate, fuel alignment plate, lower support plate etc. To verify the structural integrity of these plates, the finite element analyses are performed, which require the development of the finite element model. Sometimes it is very costly and time consuming to make the model especially for the beginners who start their engineering job for the structural analysis, necessitating a simple method to develop the finite element model for the pursuing structural analysis. Therefore in this study, the input decks are generated for the finite element analysis of a circular plate as shown in Fig. 1, which can be used for the structural analysis such as modal analysis, response spectrum analysis, stress analysis, etc using the commercial program Solid Works. The example problems are solved and the results are included for analysts to perform easily the finite element analysis of the mechanical plate components due to various loadings. The various results presented in this study would be helpful not only for the benchmark calculations and results comparisons but also as a part of the knowledge management for the future generation of young designers, scientists and computer analysts

High accuracy 3D electromagnetic finite element analysis

International Nuclear Information System (INIS)

Nelson, E.M.

1997-01-01

A high accuracy 3D electromagnetic finite element field solver employing quadratic hexahedral elements and quadratic mixed-order one-form basis functions will be described. The solver is based on an object-oriented C++ class library. Test cases demonstrate that frequency errors less than 10 ppm can be achieved using modest workstations, and that the solutions have no contamination from spurious modes. The role of differential geometry and geometrical physics in finite element analysis will also be discussed. copyright 1997 American Institute of Physics

Discrete element modeling of deformable particles in YADE

Directory of Open Access Journals (Sweden)

Martin Haustein

2017-01-01

Full Text Available In this paper we describe the open-source discrete element framework YADE and the implementation of a new deformation engine. YADE is a highly expandable software package that allows the simulation of current industrial problems in the field of granular materials using particle-based numerical methods. The description of the compaction of powders and granular material like metal pellets is now possible with a pure and simple discrete element approach in a modern DEM-framework. The deformation is realized by expanding the radius of the spherical particles, depending on their overlap, so that the volume of the material is kept constant.

Parallelized Three-Dimensional Resistivity Inversion Using Finite Elements And Adjoint State Methods

Science.gov (United States)

Schaa, Ralf; Gross, Lutz; Du Plessis, Jaco

2015-04-01

The resistivity method is one of the oldest geophysical exploration methods, which employs one pair of electrodes to inject current into the ground and one or more pairs of electrodes to measure the electrical potential difference. The potential difference is a non-linear function of the subsurface resistivity distribution described by an elliptic partial differential equation (PDE) of the Poisson type. Inversion of measured potentials solves for the subsurface resistivity represented by PDE coefficients. With increasing advances in multichannel resistivity acquisition systems (systems with more than 60 channels and full waveform recording are now emerging), inversion software require efficient storage and solver algorithms. We developed the finite element solver Escript, which provides a user-friendly programming environment in Python to solve large-scale PDE-based problems (see https://launchpad.net/escript-finley). Using finite elements, highly irregular shaped geology and topography can readily be taken into account. For the 3D resistivity problem, we have implemented the secondary potential approach, where the PDE is decomposed into a primary potential caused by the source current and the secondary potential caused by changes in subsurface resistivity. The primary potential is calculated analytically, and the boundary value problem for the secondary potential is solved using nodal finite elements. This approach removes the singularity caused by the source currents and provides more accurate 3D resistivity models. To solve the inversion problem we apply a 'first optimize then discretize' approach using the quasi-Newton scheme in form of the limited-memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) method (see Gross & Kemp 2013). The evaluation of the cost function requires the solution of the secondary potential PDE for each source current and the solution of the corresponding adjoint-state PDE for the cost function gradients with respect to the subsurface

finite element model for predicting residual stresses in shielded

African Journals Online (AJOL)

eobe

This paper investigates the prediction of residual stresses developed ... steel plates through Finite Element Model simulation and experiments. ... The experimental values as measured by the X-Ray diffractometer were of ... Based on this, it can be concluded that Finite Element .... Comparison of Residual Stresses from X.

Intercomparison of the finite difference and nodal discrete ordinates and surface flux transport methods for a LWR pool-reactor benchmark problem in X-Y geometry

International Nuclear Information System (INIS)

O'Dell, R.D.; Stepanek, J.; Wagner, M.R.

1983-01-01

The aim of the present work is to compare and discuss the three of the most advanced two dimensional transport methods, the finite difference and nodal discrete ordinates and surface flux method, incorporated into the transport codes TWODANT, TWOTRAN-NODAL, MULTIMEDIUM and SURCU. For intercomparison the eigenvalue and the neutron flux distribution are calculated using these codes in the LWR pool reactor benchmark problem. Additionally the results are compared with some results obtained by French collision probability transport codes MARSYAS and TRIDENT. Because the transport solution of this benchmark problem is close to its diffusion solution some results obtained by the finite element diffusion code FINELM and the finite difference diffusion code DIFF-2D are included

Finite-Element Modeling of Timber Joints with Punched Metal Plate Fasteners

DEFF Research Database (Denmark)

Ellegaard, Peter

2006-01-01

The focus of this paper is to describe the idea and the theory behind a finite-element model developed for analysis of timber trusses with punched metal plate fasteners (nail plates). The finite-element model includes the semirigid and nonlinear behavior of the joints (nonlinear nail and plate...... elements) and contact between timber beams, if any (bilinear contact elements). The timber beams have linear-elastic properties. The section forces needed for design of the joints are given directly by the finite-element model, since special elements are used to model the nail groups and the nail plate...... the behavior of the joints very well at lower load levels. At higher load levels the stiffness is overestimated due to development of cracks in the timber and the linear-elastic timber properties in the finite-element model....

Finite element analysis of inelastic structural behavior

International Nuclear Information System (INIS)

Argyris, J.H.; Szimmat, J.; Willam, K.J.

1977-01-01

The paper describes recent achievements in the finite element analysis of inelastic material behavior. The main purpose is to examine the interaction of three disciplines; (i) the finite element formulation of large deformation problems in the light of a systematic linearization, (ii) the constitutive modelling of inelastic processes in the rate-dependent and rate-independent response regime and (iii) the numerical solution of nonlinear rate problems via incremental iteration techniques. In the first part, alternative finite element models are developed for the idealization of large deformation problems. A systematic approach is presented to linearize the field equations locally by an incremental procedure. The finite element formulation is then examined for the description of inelastic material processes. In the second part, nonlinear and inelastic material phenomena are classified and illustrated with representative examples of concrete and metal components. In particular, rate-dependent and rate-independent material behavior is examined and representative constitutive models are assessed for their mathematical characterization. Hypoelastic, elastoplastic and endochronic models are compared for the description rate-independent material phenomena. In the third part, the numerial solution of inelastic structural behavior is discussed. In this context, several incremental techniques are developed and compared for tracing the evolution of the inelastic process. The numerical procedures are examined with regard to stability and accuracy to assess the overall efficiency. The 'optimal' incremental technique is then contrasted with the computer storage requirements to retain the data for the 'memory-characteristics' of the constitutive model

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)

Analytical and finite element modeling of grounding systems

Energy Technology Data Exchange (ETDEWEB)

Luz, Mauricio Valencia Ferreira da [University of Santa Catarina (UFSC), Florianopolis, SC (Brazil)], E-mail: mauricio@grucad.ufsc.br; Dular, Patrick [University of Liege (Belgium). Institut Montefiore], E-mail: Patrick.Dular@ulg.ac.be

2007-07-01

Grounding is the art of making an electrical connection to the earth. This paper deals with the analytical and finite element modeling of grounding systems. An electrokinetic formulation using a scalar potential can benefit from floating potentials to define global quantities such as electric voltages and currents. The application concerns a single vertical grounding with one, two and three-layer soil, where the superior extremity stays in the surface of the soil. This problem has been modeled using a 2D axi-symmetric electrokinetic formulation. The grounding resistance obtained by finite element method is compared with the analytical one for one-layer soil. With the results of this paper it is possible to show that finite element method is a powerful tool in the analysis of the grounding systems in low frequencies. (author)

Non-Linear Three Dimensional Finite Elements for Composite Concrete Structures

Directory of Open Access Journals (Sweden)

O. Kohnehpooshi

Full Text Available Abstract The current investigation focused on the development of effective and suitable modelling of reinforced concrete component with and without strengthening. The modelling includes physical and constitutive models. New interface elements have been developed, while modified constitutive law have been applied and new computational algorithm is utilised. The new elements are the Truss-link element to model the interaction between concrete and reinforcement bars, the interface element between two plate bending elements and the interface element to represent the interfacial behaviour between FRP, steel plates and concrete. Nonlinear finite-element (FE codes were developed with pre-processing. The programme was written using FORTRAN language. The accuracy and efficiency of the finite element programme were achieved by analyzing several examples from the literature. The application of the 3D FE code was further enhanced by carrying out the numerical analysis of the three dimensional finite element analysis of FRP strengthened RC beams, as well as the 3D non-linear finite element analysis of girder bridge. Acceptable distributions of slip, deflection, stresses in the concrete and FRP plate have also been found. These results show that the new elements are effective and appropriate to be used for structural component modelling.

Discrete breathers dynamic in a model for DNA chain with a finite stacking enthalpy

Science.gov (United States)

Gninzanlong, Carlos Lawrence; Ndjomatchoua, Frank Thomas; Tchawoua, Clément

2018-04-01

The nonlinear dynamics of a homogeneous DNA chain based on site-dependent finite stacking and pairing enthalpies is studied. A new variant of extended discrete nonlinear Schrödinger equation describing the dynamics of modulated wave is derived. The regions of discrete modulational instability of plane carrier waves are studied, and it appears that these zones depend strongly on the phonon frequency of Fourier's mode. The staggered/unstaggered discrete breather (SDB/USDB) is obtained straightforwardly without the staggering transformation, and it is demonstrated that SDBs are less unstable than USDB. The instability of discrete multi-humped SDB/USDB solution does not depend on the number of peaks of the discrete breather (DB). By using the concept of Peierls-Nabarro energy barrier, it appears that the low-frequency DBs are more mobile.

A fast linearized conservative finite element method for the strongly coupled nonlinear fractional Schrödinger equations

Science.gov (United States)

Li, Meng; Gu, Xian-Ming; Huang, Chengming; Fei, Mingfa; Zhang, Guoyu

2018-04-01

In this paper, a fast linearized conservative finite element method is studied for solving the strongly coupled nonlinear fractional Schrödinger equations. We prove that the scheme preserves both the mass and energy, which are defined by virtue of some recursion relationships. Using the Sobolev inequalities and then employing the mathematical induction, the discrete scheme is proved to be unconditionally convergent in the sense of L2-norm and H α / 2-norm, which means that there are no any constraints on the grid ratios. Then, the prior bound of the discrete solution in L2-norm and L∞-norm are also obtained. Moreover, we propose an iterative algorithm, by which the coefficient matrix is independent of the time level, and thus it leads to Toeplitz-like linear systems that can be efficiently solved by Krylov subspace solvers with circulant preconditioners. This method can reduce the memory requirement of the proposed linearized finite element scheme from O (M2) to O (M) and the computational complexity from O (M3) to O (Mlog ⁡ M) in each iterative step, where M is the number of grid nodes. Finally, numerical results are carried out to verify the correction of the theoretical analysis, simulate the collision of two solitary waves, and show the utility of the fast numerical solution techniques.

ZONE: a finite element mesh generator

International Nuclear Information System (INIS)

Burger, M.J.

1976-05-01

The ZONE computer program is a finite-element mesh generator which produces the nodes and element description of any two-dimensional geometry. The geometry is subdivided into a mesh of quadrilateral and triangular zones arranged sequentially in an ordered march through the geometry. The order of march can be chosen so that the minimum bandwidth is obtained. The node points are defined in terms of the x and y coordinates in a global rectangular coordinate system. The zones generated are quadrilaterals or triangles defined by four node points in a counterclockwise sequence. Node points defining the outside boundary are generated to describe pressure boundary conditions. The mesh that is generated can be used as input to any two-dimensional as well as any axisymmetrical structure program. The output from ZONE is essentially the input file to NAOS, HONDO, and other axisymmetric finite element programs. 14 figures

Finite elements methods in mechanics

CERN Document Server

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

A finite element conjugate gradient FFT method for scattering

Science.gov (United States)

Collins, Jeffery D.; Ross, Dan; Jin, J.-M.; Chatterjee, A.; Volakis, John L.

1991-01-01

Validated results are presented for the new 3D body of revolution finite element boundary integral code. A Fourier series expansion of the vector electric and mangnetic fields is employed to reduce the dimensionality of the system, and the exact boundary condition is employed to terminate the finite element mesh. The mesh termination boundary is chosen such that is leads to convolutional boundary operatores of low O(n) memory demand. Improvements of this code are discussed along with the proposed formulation for a full 3D implementation of the finite element boundary integral method in conjunction with a conjugate gradiant fast Fourier transformation (CGFFT) solution.

Analysis and Calculation of the Fluid Flow and the Temperature Field by Finite Element Modeling

Science.gov (United States)

Dhamodaran, M.; Jegadeesan, S.; Kumar, R. Praveen

2018-04-01

This paper presents a fundamental and accurate approach to study numerical analysis of fluid flow and heat transfer inside a channel. In this study, the Finite Element Method is used to analyze the channel, which is divided into small subsections. The small subsections are discretized using higher number of domain elements and the corresponding number of nodes. MATLAB codes are developed to be used in the analysis. Simulation results showed that the analyses of fluid flow and temperature are influenced significantly by the changing entrance velocity. Also, there is an apparent effect on the temperature fields due to the presence of an energy source in the middle of the domain. In this paper, the characteristics of flow analysis and heat analysis in a channel have been investigated.

A finite element method with overlapping meshes for free-boundary axisymmetric plasma equilibria in realistic geometries

Science.gov (United States)

Heumann, Holger; Rapetti, Francesca

2017-04-01

Existing finite element implementations for the computation of free-boundary axisymmetric plasma equilibria approximate the unknown poloidal flux function by standard lowest order continuous finite elements with discontinuous gradients. As a consequence, the location of critical points of the poloidal flux, that are of paramount importance in tokamak engineering, is constrained to nodes of the mesh leading to undesired jumps in transient problems. Moreover, recent numerical results for the self-consistent coupling of equilibrium with resistive diffusion and transport suggest the necessity of higher regularity when approximating the flux map. In this work we propose a mortar element method that employs two overlapping meshes. One mesh with Cartesian quadrilaterals covers the vacuum chamber domain accessible by the plasma and one mesh with triangles discretizes the region outside. The two meshes overlap in a narrow region. This approach gives the flexibility to achieve easily and at low cost higher order regularity for the approximation of the flux function in the domain covered by the plasma, while preserving accurate meshing of the geometric details outside this region. The continuity of the numerical solution in the region of overlap is weakly enforced by a mortar-like mapping.

The Cost of Continuity: Performance of Iterative Solvers on Isogeometric Finite Elements

KAUST Repository

Collier, Nathan; Dalcin, Lisandro; Pardo, David; Calo, Victor M.

2013-01-01

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.

The Cost of Continuity: Performance of Iterative Solvers on Isogeometric Finite Elements

KAUST Repository

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.

Efficient evaluation of the Coulomb force in the Gaussian and finite-element Coulomb method.

Science.gov (United States)

Kurashige, Yuki; Nakajima, Takahito; Sato, Takeshi; Hirao, Kimihiko

2010-06-28

We propose an efficient method for evaluating the Coulomb force in the Gaussian and finite-element Coulomb (GFC) method, which is a linear-scaling approach for evaluating the Coulomb matrix and energy in large molecular systems. The efficient evaluation of the analytical gradient in the GFC is not straightforward as well as the evaluation of the energy because the SCF procedure with the Coulomb matrix does not give a variational solution for the Coulomb energy. Thus, an efficient approximate method is alternatively proposed, in which the Coulomb potential is expanded in the Gaussian and finite-element auxiliary functions as done in the GFC. To minimize the error in the gradient not just in the energy, the derived functions of the original auxiliary functions of the GFC are used additionally for the evaluation of the Coulomb gradient. In fact, the use of the derived functions significantly improves the accuracy of this approach. Although these additional auxiliary functions enlarge the size of the discretized Poisson equation and thereby increase the computational cost, it maintains the near linear scaling as the GFC and does not affects the overall efficiency of the GFC approach.

Variable order spherical harmonic expansion scheme for the radiative transport equation using finite elements

International Nuclear Information System (INIS)

Surya Mohan, P.; Tarvainen, Tanja; Schweiger, Martin; Pulkkinen, Aki; Arridge, Simon R.

2011-01-01

Highlights: → We developed a variable order global basis scheme to solve light transport in 3D. → Based on finite elements, the method can be applied to a wide class of geometries. → It is computationally cheap when compared to the fixed order scheme. → Comparisons with local basis method and other models demonstrate its accuracy. → Addresses problems encountered n modeling of light transport in human brain. - Abstract: We propose the P N approximation based on a finite element framework for solving the radiative transport equation with optical tomography as the primary application area. The key idea is to employ a variable order spherical harmonic expansion for angular discretization based on the proximity to the source and the local scattering coefficient. The proposed scheme is shown to be computationally efficient compared to employing homogeneously high orders of expansion everywhere in the domain. In addition the numerical method is shown to accurately describe the void regions encountered in the forward modeling of real-life specimens such as infant brains. The accuracy of the method is demonstrated over three model problems where the P N approximation is compared against Monte Carlo simulations and other state-of-the-art methods.

Application of hexagonal element scheme in finite element method to three-dimensional diffusion problem of fast reactors

International Nuclear Information System (INIS)

Ishiguro, Misako; Higuchi, Kenji

1983-01-01

The finite element method is applied in Galerkin-type approximation to three-dimensional neutron diffusion equations of fast reactors. A hexagonal element scheme is adopted for treating the hexagonal lattice which is typical for fast reactors. The validity of the scheme is verified by applying the scheme as well as alternative schemes to the neutron diffusion calculation of a gas-cooled fast reactor of actual scale. The computed results are compared with corresponding values obtained using the currently applied triangular-element and also with conventional finite difference schemes. The hexagonal finite element scheme is found to yield a reasonable solution to the problem taken up here, with some merit in terms of saving in computing time, but the resulting multiplication factor differs by 1% and the flux by 9% compared with the triangular mesh finite difference scheme. The finite element method, even in triangular element scheme, would appear to incur error in inadmissible amount and which could not be easily eliminated by refining the nodes. (author)

WorkStream-- A Design Pattern for Multicore-Enabled Finite Element Computations

KAUST Repository

Turcksin, Bruno

2016-08-31

Many operations that need to be performed in modern finite element codes can be described as an operation that needs to be done independently on every cell, followed by a reduction of these local results into a global data structure. For example, matrix assembly, estimating discretization errors, or converting nodal values into data structures that can be output in visualization file formats all fall into this class of operations. Using this realization, we identify a software design pattern that we callWorkStream and that can be used to model such operations and enables the use of multicore shared memory parallel processing. We also describe in detail how this design pattern can be efficiently implemented, and we provide numerical scalability results from its use in the DEAL.II software library.

Variational integrators for the dynamics of thermo-elastic solids with finite speed thermal waves

International Nuclear Information System (INIS)

Mata, Pablo; Lew, Adrian J.

2014-01-01

This paper formulates variational integrators for finite element discretizations of deformable bodies with heat conduction in the form of finite speed thermal waves. The cornerstone of the construction consists in taking advantage of the fact that the Green–Naghdi theory of type II for thermo-elastic solids has a Hamiltonian structure. Thus, standard techniques to construct variational integrators can be applied to finite element discretizations of the problem. The resulting discrete-in-time trajectories are then consistent with the laws of thermodynamics for these systems: for an isolated system, they exactly conserve the total entropy, and nearly exactly conserve the total energy over exponentially long periods of time. Moreover, linear and angular momenta are also exactly conserved whenever the exact system does. For definiteness, we construct an explicit second-order accurate algorithm for affine tetrahedral elements in two and three dimensions, and demonstrate its performance with numerical examples

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

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

Modelling Dowel Action of Discrete Reinforcing Bars in Cracked Concrete Structures

International Nuclear Information System (INIS)

Kwan, A. K. H.; Ng, P. L.; Lam, J. Y. K.

2010-01-01

Dowel action is one of the component actions for shear force transfer in cracked reinforced concrete. In finite element analysis of concrete structures, the use of discrete representation of reinforcing bars is considered advantageous over the smeared representation due to the relative ease of modelling the bond-slip behaviour. However, there is very limited research on how to simulate the dowel action of discrete reinforcing bars. Herein, a numerical model for dowel action of discrete reinforcing bars crossing cracks in concrete is developed. The model features the derivation of dowel stiffness matrix based on beam-on-elastic-foundation theory and the direct assemblage of dowel stiffness into the concrete element stiffness matrices. The dowel action model is incorporated in a nonlinear finite element programme with secant stiffness formulation. Deep beams tested in the literature are analysed and it is found that the incorporation of dowel action model improves the accuracy of analysis.

Development of a partitioned finite volume-finite element fluid-structure interaction scheme for strongly-coupled problems

CSIR Research Space (South Africa)

Suliman, Ridhwaan

2012-07-01

Full Text Available -linear deformations are accounted for. As will be demonstrated, the finite volume approach exhibits similar disad- vantages to the linear Q4 finite element formulation when undergoing bending. An enhanced finite volume approach is discussed and compared with finite...

The finite element analysis program MSC Marc/Mentat a first introduction

CERN Document Server

Öchsner, Andreas

2016-01-01

Based on simple examples, this book offers a short introduction to the general-purpose finite element program MSC Marc, a specialized program for non-linear problems (implicit solver) distributed by the MSC Software Corporation, which is commonly used in academia and industry. Today the documentation of all finite element programs includes a variety of step-by-step examples of differing complexity, and in addition, all software companies offer professional workshops on different topics. As such, rather than competing with these, the book focuses on providing simple examples, often single-element problems, which can easily be related to the theory that is discussed in finite element lectures. This makes it an ideal companion book to classical introductory courses on the finite element method.

GPU-based discrete element rigid body transport

CSIR Research Space (South Africa)

Govender, Nicolin

2013-08-01

Full Text Available . For applications in coastal engineering and also in pavement engineering, the capture of particle shapes as polyhedra rather than clumped spheres is particularly important. The development of a Discrete Element Model applicable to both fields, and to industrial...

A finite element thermohydrodynamic analyis of profile bore bearing

International Nuclear Information System (INIS)

Shah Nor bin Basri

1994-01-01

A finite element-based method is presented for analysing the thermohydrodynamic (THD) behaviour of profile bore bearing. A variational statement for the governing equation is derived and used to formulate a non-linear quadrilateral finite element of serendipity family. The predicted behaviour is compared with experimental evidence where possible and favorable correlation is obtained

Comparison between results of solution of Burgers' equation and Laplace's equation by Galerkin and least-square finite element methods

Science.gov (United States)

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.

On approximation of non-Newtonian fluid flow by the finite element method

Science.gov (United States)

Svácek, Petr

2008-08-01

In this paper the problem of numerical approximation of non-Newtonian fluid flow with free surface is considered. Namely, the flow of fresh concrete is addressed. Industrial mixtures often behaves like non-Newtonian fluids exhibiting a yield stress that needs to be overcome for the flow to take place, cf. [R.B. Bird, R.C. Armstrong, O. Hassager, Dynamics of Polymeric Liquids, vol. 1, Fluid Mechanics, Wiley, New York, 1987; R.P. Chhabra, J.F. Richardson, Non-Newtonian Flow in the Process Industries, Butterworth-Heinemann, London, 1999]. The main interest is paid to the mathematical formulation of the problem and to discretization with the aid of finite element method. The described numerical procedure is applied onto the solution of several problems.

Hualien forced vibration calculation with a finite element model

International Nuclear Information System (INIS)

Wang, F.; Gantenbein, F.; Nedelec, M.; Duretz, Ch.

1995-01-01

The forced vibration tests of the Hualien mock-up were useful to validate finite element models developed for soil-structure interaction. In this paper the two sets of tests with and without backfill were analysed. the methods used are based on finite element modeling for the soil. Two approaches were considered: calculation of soil impedance followed by the calculation of the transfer functions with a model taking into account the superstructure and the impedance; direct calculation of the soil-structure transfer functions, with the soil and the structure being represented in the same model by finite elements. Blind predictions and post-test calculations are presented and compared with the test results. (author). 4 refs., 8 figs., 2 tabs

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.