Sample records for finite element method

  1. Finite Volume Element Method

    E-Print Network [OSTI]

    2003-12-06T23:59:59.000Z

    FVE is closely related to the control volume finite element method ... simple stencils, to apply to a fairly wide range of fluid flow equations, to effectively treat.

  2. ME 872 -Finite Element Methods Spring 2014

    E-Print Network [OSTI]

    Diaz, Alejandro

    Element Method: Linear Static and Dynamic Finite Element Analysis (Dover Civil and Mechanical Engineering problems Special topics: Lagrange multipliers, adaptive finite elements, sensitivity analysis, nonlinearME 872 - Finite Element Methods Spring 2014 Catalog Description: Theory and application

  3. Combined Finite Element --Finite Volume Method ( Convergence Analysis )

    E-Print Network [OSTI]

    Magdeburg, Universität

    Combined Finite Element -- Finite Volume Method ( Convergence Analysis ) M'aria Luk idea is to combine finite volume and finite element methods in an appropriate way. Thus nonlinear grid. Diffusion terms are discretized by the conforming piecewise linear finite element method

  4. INTRODUCTION TO THE FINITE ELEMENT METHOD

    E-Print Network [OSTI]

    California at Berkeley, University of

    in popularity with the finite difference method in the areas of heat transfer and fluid mechanics. 1DRAFT Chapter 1 INTRODUCTION TO THE FINITE ELEMENT METHOD 1.1 Historical perspective: the origins of the finite el- ement method The finite element method constitutes a general tool for the numerical solution

  5. Generalized finite element method for Helmholtz equation 

    E-Print Network [OSTI]

    Hidajat, Realino Lulie

    2009-05-15T23:59:59.000Z

    reduces the pollution effect due to wave number and we are able to obtain a highly accurate solution with a much smaller number of degrees-of-freedom compared with the classical Finite Element Method. The q-convergence of the handbook functions...

  6. Math/AMath 595: Finite Element Methods (Winter 2001)

    E-Print Network [OSTI]

    Stein, William

    for Elliptic Problems by Philippe G. Ciarlet (1978). 3. An Analysis of the Finite Element Method by GilbertMath/AMath 595: Finite Element Methods (Winter 2001) Lectures: MW 3:30{4:50, room 121 RAI Professor Element Methods by Susanne C. Brenner and L. Ridgway Scott (1994). 2. The Finite Element Method

  7. Finite Element Methods for Radiosity \\Lambda Paul S. Heckbert

    E-Print Network [OSTI]

    Treuille, Adrien

    Finite Element Methods for Radiosity \\Lambda Paul S. Heckbert Carnegie Mellon University Global global illu­ mination: finite element methods and Monte Carlo methods. The former approach yields radiosity methods and the latter approach yields ray tracing methods, generally speaking. Finite element

  8. Interval Finite Element Methods: New Directions Rafi Muhanna1

    E-Print Network [OSTI]

    Kreinovich, Vladik

    Interval Finite Element Methods: New Directions Rafi Muhanna1 , Vladik Kreinovich2 , Pavel Sol´in2 advanced numerical methods to tackle these problems, such as finite element methods (FEM). Lately, new advanced version of these methods have appeared, such as hierarchic higher-order finite element methods (hp

  9. ANALYSIS OF SOME MOVING SPACETIME FINITE ELEMENT METHODS \\Lambda

    E-Print Network [OSTI]

    Bank, Randolph E.

    ANALYSIS OF SOME MOVING SPACE­TIME FINITE ELEMENT METHODS \\Lambda RANDOLPH E. BANK y AND RAFAEL F. SANTOS z Abstract. Two space­time finite element methods for solving time­dependent partial differential equations are defined and analyzed. The methods are based on the use of isoparametric finite elements

  10. Introduction to Finite Element Methods \\Lambda Paul S. Heckbert

    E-Print Network [OSTI]

    Treuille, Adrien

    Introduction to Finite Element Methods \\Lambda Paul S. Heckbert Carnegie Mellon University Global be solved numerically. In this chapter we explain some of the basics of the finite element method Illumination course notes, SIGGRAPH '93, Anaheim, Aug. 1993 4­1 #12; 4­2 Introduction to Finite Element Methods

  11. AN ADAPTIVE SURFACE FINITE ELEMENT METHOD BASED ON VOLUME MESHES

    E-Print Network [OSTI]

    Demlow, Alan

    studied in this paper. The paper of Dziuk [12] contains the first analysis of a finite element methodAN ADAPTIVE SURFACE FINITE ELEMENT METHOD BASED ON VOLUME MESHES ALAN DEMLOW AND MAXIM A. OLSHANSKII Abstract. In this paper we define an adaptive version of a recently introduced finite element

  12. A finite element method and the method of finite spheres enriched for analysis of wave propagations

    E-Print Network [OSTI]

    Ham, Seounghyun, 1982-

    2014-01-01T23:59:59.000Z

    The objective of this thesis is to present a finite element method and the method of finite spheres enriched for the solution of various wave propagation problems. The first part of this thesis is to present an enriched ...

  13. Complex Wavenumber Fourier Analysis of the PVersion Finite Element Method

    E-Print Network [OSTI]

    Thompson, Lonny L.

    Complex Wavenumber Fourier Analysis of the P­Version Finite Element Method Lonny L. Thompson wave solutions. This type of analysis amounts to a discrete Fourier synthesis of the finite element bands in the context of finite element analysis were not fully understood, although there is some

  14. Analysis of a Combined Barycentric Finite Volume --Nonconforming Finite Element Method for Nonlinear

    E-Print Network [OSTI]

    Dolejsi, Vit

    Analysis of a Combined Barycentric Finite Volume -- Nonconforming Finite Element Method, the theoretical analysis of the combined finite volume -- finite element method has been carried out for the case term is discretized by piecewise linear nonconforming triangular finite elements. Under the assumption

  15. Adaptive finite element method for eddy current problems

    E-Print Network [OSTI]

    Adaptive finite element method for eddy current problems. In this chapter we shall discuss Maxwell's equations at very low frequencies. We also refer to [bo98 ...

  16. LEAST-SQUARES FINITE ELEMENT METHODS FOR OPTIMALITY SYSTEMS ARISING

    E-Print Network [OSTI]

    Bochev, Pavel

    LEAST-SQUARES FINITE ELEMENT METHODS FOR OPTIMALITY SYSTEMS ARISING IN OPTIMIZATION AND CONTROL control problems for systems governed by linear, elliptic partial differential equations is considered modern least-squares finite element methods for the solution of the optimality systems. The matrix equa

  17. ADAPTIVE FINITE ELEMENT FREQUENCY DOMAIN METHOD FOR EDDY CURRENT PROBLEMS

    E-Print Network [OSTI]

    Zheng, Weiying

    ADAPTIVE FINITE ELEMENT FREQUENCY DOMAIN METHOD FOR EDDY CURRENT PROBLEMS WEIYING ZHENG-harmonic eddy current problems in the case of three-dimensional isotropic and linear materials. We adopt. Time-harmonic Maxwell's equations, eddy current, adaptive finite element method, multiply connected

  18. Direct Forcing Immersed Boundary Methods: Finite Element Versus Finite Volume Approach

    E-Print Network [OSTI]

    Frisani, Angelo 1980-

    2012-12-07T23:59:59.000Z

    Two immersed boundary methods (IBM) for the simulation of conjugate heat transfer problems with complex geometries are introduced: a finite element (IFEM) and a finite volume (IFVM) immersed boundary methods are discussed. In the IFEM a projection...

  19. A Note on the Finite Element Method with Singular Basis Functions

    E-Print Network [OSTI]

    Kaneko, Hideaki

    finite element analysis that incorporates singular element functions. A need for introducing * *some singular elements as part of basis functions in certain finite element analysis arises o* *ut A Note on the Finite Element Method with Singular Basis

  20. PARTIALLY PENALIZED IMMERSED FINITE ELEMENT METHODS ...

    E-Print Network [OSTI]

    2015-01-05T23:59:59.000Z

    each interface element, it uses IFE functions constructed with piecewise ...... strate this, we plot errors of a classic bilinear IFE solution and a NPP IFE solution in.

  1. Coupled finite element and meshfree method Application to wide-flange steel sections

    E-Print Network [OSTI]

    Yaw, Louis L.

    , 2007 Yaw, Kunnath & Sukumar Meshfree Inelastic Frame Analysis #12;Motivation Coupled finite element Inelastic Frame Analysis #12;Motivation Coupled finite element and meshfree method Application to wide & Sukumar Meshfree Inelastic Frame Analysis #12;Motivation Coupled finite element and meshfree method

  2. Generalized finite element method for multiscale analysis

    E-Print Network [OSTI]

    Zhang, Lin

    2004-11-15T23:59:59.000Z

    by the Partition of Unity Method (PUM). It is shown that the p-version of the Generalized FEM using mesh-based handbook functions is capable of achieving very high accuracy. It is also analyzed that the effect of the main factors affecting the accuracy...

  3. Discontinuous Galerkin finite element methods for gradient plasticity.

    SciTech Connect (OSTI)

    Garikipati, Krishna. (University of Michigan, Ann Arbor, MI); Ostien, Jakob T.

    2010-10-01T23:59:59.000Z

    In this report we apply discontinuous Galerkin finite element methods to the equations of an incompatibility based formulation of gradient plasticity. The presentation is motivated with a brief overview of the description of dislocations within a crystal lattice. A tensor representing a measure of the incompatibility with the lattice is used in the formulation of a gradient plasticity model. This model is cast in a variational formulation, and discontinuous Galerkin machinery is employed to implement the formulation into a finite element code. Finally numerical examples of the model are shown.

  4. A ROBUST FINITE ELEMENT METHOD FOR DARCY--STOKES FLOW

    E-Print Network [OSTI]

    Winther, Ragnar

    . When # = 0 the first equation in (1.1) has the form of Darcy's law for flow in a homogeneous porousA ROBUST FINITE ELEMENT METHOD FOR DARCY--STOKES FLOW KENT ANDRE MARDAL, XUE--CHENG TAI, AND RAGNAR the perturbation parameter is large, while it degenerates to a mixed formulation of Poisson's equation

  5. An Efficient Probabilistic Finite Element Method for Stochastic Groundwater Flow

    E-Print Network [OSTI]

    Osnes, Harald

    of groundwater flow and transport from mathematical models involves considerable uncertainty due to the presenceAn Efficient Probabilistic Finite Element Method for Stochastic Groundwater Flow Harald Osnes Hans by nonlinear effects as oe Y approaches unity in more general domains. 1 #12; 1 Introduction The prediction

  6. FINITE ELEMENT METHOD IN FLUID MECHANICS & HEAT TRANSFER

    E-Print Network [OSTI]

    Camci, Cengiz

    completed this course should be able to perform quick analysis of small problems using the finite element of Fluid Mechanics and Heat Transfer An Introduction to Finite Element Analysis Using "Galerkin Weak of Euler's Equation in Finite Element Analysis Generalized Form of Euler's Equation in Three Dimensional

  7. The Mimetic Finite Element Method and the Virtual Element Method for elliptic problems with arbitrary regularity.

    SciTech Connect (OSTI)

    Manzini, Gianmarco [Los Alamos National Laboratory

    2012-07-13T23:59:59.000Z

    We develop and analyze a new family of virtual element methods on unstructured polygonal meshes for the diffusion problem in primal form, that use arbitrarily regular discrete spaces V{sub h} {contained_in} C{sup {alpha}} {element_of} N. The degrees of freedom are (a) solution and derivative values of various degree at suitable nodes and (b) solution moments inside polygons. The convergence of the method is proven theoretically and an optimal error estimate is derived. The connection with the Mimetic Finite Difference method is also discussed. Numerical experiments confirm the convergence rate that is expected from the theory.

  8. THE BUBBLE TRANSFORM: A NEW TOOL FOR ANALYSIS OF FINITE ELEMENT METHODS

    E-Print Network [OSTI]

    THE BUBBLE TRANSFORM: A NEW TOOL FOR ANALYSIS OF FINITE ELEMENT METHODS RICHARD S. FALK AND RAGNAR estimates for various finite element methods. In particular, for the analysis of mixed finite element be very useful, but hard to construct, is the analysis of the so-called p-version of the finite element

  9. An application of stochastic finite element method in structural reliabiity

    E-Print Network [OSTI]

    Agrawal, Sanjay

    1993-01-01T23:59:59.000Z

    OF SCIENCE August 1993 Major Subject: Civil Engineering AN APPLICATION OF STOCHASTIC FINITE ELEMENT METHOD IN STRUCTURAL RELIABILITY A Thesis by SANJAY AGRAWAL Approved as to style and content by: Kostas apadimitriou (Chair of Committee) John M... of the beam on the buckling load. The statistical methodology can be efficiently used with available FORM/SORM reliability techniques to compute the probability of failure of the beam due to buckling. Random fields are used to quantify the modeling...

  10. Finite Element Methods for Nonlinear Acoustics in Fluids

    E-Print Network [OSTI]

    The non- linear terms in these wave equations imply that the sound speed ... ics, the finite element formulation of nonlinear acoustic equations for fluids has.

  11. Adaptive Finite Element Method for Time-domain Maxwell's ...

    E-Print Network [OSTI]

    Shuhao Cao supervised under Professor Cai

    2013-10-31T23:59:59.000Z

    Oct 11, 2013 ... Avoiding the pollution of the finite element approximation, a.k.a. the spurious modes. Lowest-order element's DoFs associated with edges, and.

  12. THE STEFAN PROBLEM: A STOCHASTIC ANALYSIS USING THE EXTENDED FINITE ELEMENT METHOD

    E-Print Network [OSTI]

    Zabaras, Nicholas J.

    THE STEFAN PROBLEM: A STOCHASTIC ANALYSIS USING THE EXTENDED FINITE ELEMENT METHOD Baskar these uncertainties in a stochastic sense using the extended finite element method as the deterministic solution convection using the extended finite element method (XFEM) and level set methods," Journal of Computational

  13. Adaptive Finite Element Method for Solving the Exact Kohn-Sham...

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Finite Element Method for Solving the Exact Kohn-Sham Equation of Density Functional Theory. Adaptive Finite Element Method for Solving the Exact Kohn-Sham Equation of Density...

  14. Stabilized finite element methods with fast iterative solution ...

    E-Print Network [OSTI]

    2003-04-11T23:59:59.000Z

    Computer methods in applied mechanics and engineering. Stabilized finite ... Department of Mathemutics, Purdue Universi@, 1395 Mathematical Science ...

  15. ANALYSIS OF A FINITE VOLUME ELEMENT METHOD FOR THE STOKES PROBLEM

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    ANALYSIS OF A FINITE VOLUME ELEMENT METHOD FOR THE STOKES PROBLEM ALFIO QUARTERONI, AND RICARDO stabilization, finite volume element method, a priori error estimates, superconvergence analysis. CMCS RUIZ BAIER Abstract. In this paper we propose a stabilized conforming finite volume element method

  16. A multiresolution finite element method based on a new locking-free rectangular Mindlin plate element

    E-Print Network [OSTI]

    Xia, Yi-Ming

    2015-01-01T23:59:59.000Z

    A locking-free rectangular Mindlin plate element with a new multi-resolution analysis (MRA) is proposed and a multireolution finite element method is hence presented. The MRA framework is formulated out of a mutually nesting displacement subspace sequence. The MRA endows the proposed element with the resolution level (RL) to adjust the element node number, thus modulating structural analysis accuracy accordingly. As a result, the traditional 4-node rectangular Mindlin plate element and method is a mono-resolution one and also a special case of the proposed element and method. The meshing for the monoresolution plate element model is based on the empiricism while the RL adjusting for the multiresolution is laid on the rigorous mathematical basis. The accuracy of a structural analysis is actually determined by the RL, not by the mesh. The rational MRA enables the implementation of the multiresolution Mindlin plate element method to be more rational and efficient than that of the conventional monoresolution or o...

  17. Analysis of Conforming and Nonconforming Quadrilateral Finite Element Methods for the

    E-Print Network [OSTI]

    Sheen, Dongwoo

    Analysis of Conforming and Nonconforming Quadrilateral Finite Element Methods for the Helmholtz in this analysis are the standard Q1 conforming finite element, the DSSY nonconforming element [5] and the P1-Santos-Sheen-Bennethum in [4]; also in the same paper a first rigorous finite element error analysis was carried out

  18. EXPANDED MIXED FINITE ELEMENT METHODS FOR QUASILINEAR SECOND ORDER ELLIPTIC PROBLEMS, II

    E-Print Network [OSTI]

    EXPANDED MIXED FINITE ELEMENT METHODS FOR QUASILINEAR SECOND ORDER ELLIPTIC PROBLEMS, II Zhangxin coefficient times the gradient). Based on this formulation, mixed finite element approximations with low permeability and to the transport equation with small dispersion in groundwater modeling

  19. A SpaceTime Finite Element Method for the Exterior Structural Acoustics Problem

    E-Print Network [OSTI]

    Thompson, Lonny L.

    A Space­Time Finite Element Method for the Exterior Structural Acoustics Problem: Time­discontinuous Galerkin space­time finite element method is formu­ lated for the exterior structural acoustics problem Introduction A space­time finite element formulation is presented for solution of the exterior struc­ tural

  20. A PARALLEL GEOMETRIC MULTIGRID METHOD FOR FINITE ELEMENTS ON OCTREE MESHES

    E-Print Network [OSTI]

    Biros, George

    A PARALLEL GEOMETRIC MULTIGRID METHOD FOR FINITE ELEMENTS ON OCTREE MESHES RAHUL S. SAMPATH for the discretized finite element operators and the intergrid transfer operations. The key component of our scheme it has a low setup cost. Key words. Geometric Multigrid, Finite Element Method, Linear Octrees, Adaptive

  1. L.L.Thompson: Finite element methods for acoustics, Preprint: J.Acoust.Soc.Am. A review of finite element methods for time-harmonic acoustics

    E-Print Network [OSTI]

    Thompson, Lonny L.

    L.L.Thompson: Finite element methods for acoustics, Preprint: J.Acoust.Soc.Am. A review of finite element methods for time-harmonic acoustics Lonny L. Thompson Department of Mechanical Engineering, Clemson University Clemson, South Carolina, 29634-0921, USA Email: lonny.thompson@ces.clemson.edu (Dated

  2. Finite size scaling for quantum criticality using the finite-element method

    E-Print Network [OSTI]

    Edwin Antillon; Birgit Wehefritz-Kaufmann; Sabre Kais

    2012-03-08T23:59:59.000Z

    Finite size scaling for the Schr\\"{o}dinger equation is a systematic approach to calculate the quantum critical parameters for a given Hamiltonian. This approach has been shown to give very accurate results for critical parameters by using a systematic expansion with global basis-type functions. Recently, the finite element method was shown to be a powerful numerical method for ab initio electronic structure calculations with a variable real-space resolution. In this work, we demonstrate how to obtain quantum critical parameters by combining the finite element method (FEM) with finite size scaling (FSS) using different ab initio approximations and exact formulations. The critical parameters could be atomic nuclear charges, internuclear distances, electron density, disorder, lattice structure, and external fields for stability of atomic, molecular systems and quantum phase transitions of extended systems. To illustrate the effectiveness of this approach we provide detailed calculations of applying FEM to approximate solutions for the two-electron atom with varying nuclear charge; these include Hartree-Fock, density functional theory under the local density approximation, and an "exact"' formulation using FEM. We then use the FSS approach to determine its critical nuclear charge for stability; here, the size of the system is related to the number of elements used in the calculations. Results prove to be in good agreement with previous Slater-basis set calculations and demonstrate that it is possible to combine finite size scaling with the finite-element method by using ab initio calculations to obtain quantum critical parameters. The combined approach provides a promising first-principles approach to describe quantum phase transitions for materials and extended systems.

  3. Control volume finite element method with multidimensional edge element Scharfetter-Gummel upwinding. Part 1, formulation.

    SciTech Connect (OSTI)

    Bochev, Pavel Blagoveston

    2011-06-01T23:59:59.000Z

    We develop a new formulation of the Control Volume Finite Element Method (CVFEM) with a multidimensional Scharfetter-Gummel (SG) upwinding for the drift-diffusion equations. The formulation uses standard nodal elements for the concentrations and expands the flux in terms of the lowest-order Nedelec H(curl; {Omega})-compatible finite element basis. The SG formula is applied to the edges of the elements to express the Nedelec element degree of freedom on this edge in terms of the nodal degrees of freedom associated with the endpoints of the edge. The resulting upwind flux incorporates the upwind effects from all edges and is defined at the interior of the element. This allows for accurate evaluation of integrals on the boundaries of the control volumes for arbitrary quadrilateral elements. The new formulation admits efficient implementation through a standard loop over the elements in the mesh followed by loops over the element nodes (associated with control volume fractions in the element) and element edges (associated with flux degrees of freedom). The quantities required for the SG formula can be precomputed and stored for each edge in the mesh for additional efficiency gains. For clarity the details are presented for two-dimensional quadrilateral grids. Extension to other element shapes and three dimensions is straightforward.

  4. FINITE ELEMENT METHOD FOR SUBSURFACE HYDROLOGY USING A MIXED EXPLICIT-IMPLICIT SCHEME

    E-Print Network [OSTI]

    Narasimhan, T.N.

    2011-01-01T23:59:59.000Z

    METHOD FOR SUBSURFACE HYDROLOGY USING A MIXED EXPLICIT-arising in subsurface hydrology. These problems includeFinite Element Method in Hydrology," Int. Jour. Num. Meth.

  5. Developments in the method of finite spheres : efficiency and coupling to the traditional finite element method

    E-Print Network [OSTI]

    Hong, Jung-Wuk, 1970-

    2004-01-01T23:59:59.000Z

    In this thesis we develop some advances in the method of finite spheres which is a truly meshless numerical technique for the solution of boundary value problems on geometrically complex domains. We present the development ...

  6. Analysis of the finite element heterogeneous multiscale method for nonlinear elliptic homogenization problems.

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Analysis of the finite element heterogeneous multiscale method for nonlinear elliptic homogenization problems. Assyr Abdulle and Gilles Vilmart September 28, 2012 Abstract An analysis of the finite finite elements. Op- timal a-priori error estimates are obtained for the H1 and L2 norms, error bounds

  7. AN ADAPTIVE FINITE ELEMENT METHOD FOR THE H -FORMULATION OF TIME-DEPENDENT EDDY CURRENT

    E-Print Network [OSTI]

    AN ADAPTIVE FINITE ELEMENT METHOD FOR THE H - FORMULATION OF TIME-DEPENDENT EDDY CURRENT PROBLEMS, eddy current problem, adaptive finite element method, mul- tiply connected conductor, Team Workshop Problem 7 AMS subject classifications. 65M60, 65M50, 78A25 1. Introduction. Eddy currents appear in almost

  8. A finite element method for transient analysis of concurrent large deformation and mass transport in gels

    E-Print Network [OSTI]

    Suo, Zhigang

    A finite element method for transient analysis of concurrent large deformation and mass transport and Aerospace Engineering, Arizona State University, Tempe, Arizona 85287, USA 2 School of Engineering-dependent concurrent process of large deformation and mass transport is studied by developing a finite element method

  9. A Space-Time Finite Element Method for the Exterior Acoustics Problem

    E-Print Network [OSTI]

    Thompson, Lonny L.

    A Space-Time Finite Element Method for the Exterior Acoustics Problem Lonny L. Thompson Department in exterior domains is discussed. The space-time formulation for the exterior acoustics problem is obtained, the development of a space-time finite element method for so- lution of the transient acoustics problem

  10. A Space-Time Finite Element Method for the Exterior Structural Acoustics Problem

    E-Print Network [OSTI]

    Thompson, Lonny L.

    A Space-Time Finite Element Method for the Exterior Structural Acoustics Problem: Time-time finite element method is formu- lated for the exterior structural acoustics problem in two space formulation is presented for solution of the exterior struc- tural acoustics problem in two space dimensions

  11. Adaptive finite element methods for PDE-constrained optimal control problems

    E-Print Network [OSTI]

    Adaptive finite element methods for PDE-constrained optimal control problems R. Becker1 , M. Braack is discretized by a finite element Galerkin method. The accuracy of the discretization is controlled by residual of heat transfer, fluid flow and parameter estimation. The contents of this article is as follows

  12. Numerical Analysis of a Mixed Finite Element Method for a FlowTransport Problem

    E-Print Network [OSTI]

    Chou, So-Hsiang

    element method for the brine, radionuclides, and heat combined with a mixed finite element method dissolves to create a brine, radioactive elements decay to generate heat, and radionuclides are transported), the viscosity of the fluid is dependent on â?? c, the concentration of the brine in the fluid. Here and below

  13. Interval Finite Element Methods: New Directions Rafi Muhanna 1 , Vladik Kreinovich 2 , Pavel Sol'in 2 ,

    E-Print Network [OSTI]

    Ward, Karen

    Interval Finite Element Methods: New Directions Rafi Muhanna 1 , Vladik Kreinovich 2 , Pavel Ÿ Sol advanced numerical methods to tackle these problems, such as finite element methods (FEM). Lately, new advanced version of these methods have appeared, such as hierarchic higher­order finite element methods (hp

  14. An implicit finite element method for discrete dynamic fracture

    SciTech Connect (OSTI)

    Jobie M. Gerken

    1999-12-01T23:59:59.000Z

    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 error is expected.

  15. DIRECT AND INVERSE APPROXIMATION THEOREMS FOR THE p-VERSION OF THE FINITE ELEMENT METHOD IN THE

    E-Print Network [OSTI]

    Guo, Benqi

    of the finite element method, for example, MSC/PROBE, FIESTA, MECHANICA, PHLEX, STRESSCHECK, and STRIPE. In 1980

  16. IMA Journal of Numerical Analysis (1988) 8, 415-433 Finite Element Methods for a Model for Full Waveform Acoustic

    E-Print Network [OSTI]

    Harrison, Mark

    1988-01-01T23:59:59.000Z

    IMA Journal of Numerical Analysis (1988) 8, 415-433 Finite Element Methods for a Model for Full are given and then a discrete-time, explicit finite element procedure is defined and analysed, with finite on the existence, uniqueness and finite element approximation of the solution of Biot's equations were given in [15

  17. A Family of Uniform Strain Tetrahedral Elements and a Method for Connecting Dissimilar Finite Element Meshes

    SciTech Connect (OSTI)

    Dohrmann, C.R.; Heinstein, M.W.; Jung, J.; Key, S.W.

    1999-01-01T23:59:59.000Z

    This report documents a collection of papers on a family of uniform strain tetrahedral finite elements and their connection to different element types. Also included in the report are two papers which address the general problem of connecting dissimilar meshes in two and three dimensions. Much of the work presented here was motivated by the development of the tetrahedral element described in the report "A Suitable Low-Order, Eight-Node Tetrahedral Finite Element For Solids," by S. W. Key {ital et al.}, SAND98-0756, March 1998. Two basic issues addressed by the papers are: (1) the performance of alternative tetrahedral elements with uniform strain and enhanced uniform strain formulations, and (2) the proper connection of tetrahedral and other element types when two meshes are "tied" together to represent a single continuous domain.

  18. A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION

    E-Print Network [OSTI]

    Mottram, Nigel

    A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA results. 1. Introduction The Darcy equation arising in a porous media field belongs to the family of mixed element methods [23, 21, 13, 9, 10] the range of possibilities to tackle the Darcy equation has increased

  19. A SpaceTime Finite Element Method for the Exterior Acoustics Problem

    E-Print Network [OSTI]

    Thompson, Lonny L.

    A Space­Time Finite Element Method for the Exterior Acoustics Problem by Lonny L. Thompson problem in exterior domains is discussed. The space­time formulation for the exterior acoustics problem acoustics problem. i #12; Contents 1 Introduction 1 2 The Exterior Acoustics Problem 3 3 Space­time finite

  20. Decoupled Modeling of Chilled Water Cooling Coils Using a Finite Element Method

    E-Print Network [OSTI]

    Wang, G.; Liu, M.

    2005-01-01T23:59:59.000Z

    Decoupled Modeling of Chilled Water Cooling Coils Using a Finite Element Method Gang Wang Research Associate University of Nebraska – Lincoln Mingsheng Liu Professor University of Nebraska – Lincoln David E. Claridge Professor Texas A... be decoupled using a constant sensible heat ratio (SHR) and the saturation humidity ratio vs. temperature curve can be treated as linear in a small area corresponding to a finite element of the coil. This paper presents the decoupled cooling coil model...

  1. Mixed finite element methods for incompressible flow: Stationary ...

    E-Print Network [OSTI]

    2010-04-14T23:59:59.000Z

    May 14, 2009 ... Received 2 October 2008; accepted 27 January 2009 ... Thomas (RT) element of index k ? 0 and the velocity by piecewise .... We use the standard notations and definitions for the Sobolev spaces Hs( )d and Hs(? )d for.

  2. Numerical analysis of a finite element method for the axisymmetric eddy current model of an induction furnace

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    Numerical analysis of a finite element method for the axisymmetric eddy current model, 27002, Lugo, Spain The aim of this paper is to analyze a finite element method to solve an eddy current of the method are reported. Keywords: low-frequency harmonic Maxwell equations, eddy current problems, finite

  3. 1 INTRODUCTION An application of the finite element method (FEM) for non-linear elastoplastic analysis of rein-

    E-Print Network [OSTI]

    Zornberg, Jorge G.

    1391 1 INTRODUCTION An application of the finite element method (FEM) for non-linear elastoplastic, 2003; Oliveira, 2006). 2 FINITE ELEMENT REPRESENTATION OF REINFORCED SOIL A discrete representation, the reinforcement and the soil-reinforcement interface-- can be represented using a specific finite element with its

  4. Uncoupled thermoelastic structural response by the finite element method

    E-Print Network [OSTI]

    Linn, William Ray

    1982-01-01T23:59:59.000Z

    constant strain triangle (csT) element, when the solid mechanics problem is examined. Nonisothermal effects are obtained from the discretized form of equation (54), written as U . { {' [D] , {, { av - ? z f {*{ [Dj[BI {Q) dr (63) 2 J V V with 31 2 P... partial spacial derivatives and repeated indices imply summation), (k T ) P T[1+ 5( 1 & )( ? )] QT where (3X+2u) o T 6 6 v e (3) The style and citations of this thesis follow that of the 1 t'16166'1MttdJ~Eg where km is the thermal conductivity...

  5. AN ADAPTIVE FINITE ELEMENT METHOD FOR THE EDDY CURRENT MODEL WITH CIRCUIT/FIELD COUPLINGS

    E-Print Network [OSTI]

    AN ADAPTIVE FINITE ELEMENT METHOD FOR THE EDDY CURRENT MODEL WITH CIRCUIT/FIELD COUPLINGS JUNQING for solving the eddy current model with voltage excitations for complicated three dimensional structures of the proposed method. Key words. Eddy current, circuit/field coupling, adaptivity, a posteriori error analysis

  6. A VECTOR FINITE ELEMENT TIME-DOMAIN METHOD FOR ...

    E-Print Network [OSTI]

    SIAM (#1) 1035 2001 Apr 10 12:32:38

    2001-08-14T23:59:59.000Z

    and to conserve energy and charge for arbitrary hexahedral grids. ... experiments are performed to determine the accuracy and efficiency of the method.

  7. spectral stochastic finite-element methods for parametric ...

    E-Print Network [OSTI]

    Mihai Anitescu

    2006-11-30T23:59:59.000Z

    Nov 30, 2006 ... Abstract: We present a method to approximate the solution mapping of ... Category 3: Applications -- Science and Engineering (Mechanical Engineering ) ... Mathematics and Computer Science Division, 9700 S Cass Avenue, ...

  8. Advanced finite-element methods for design and analysis of nanooptical structures: Applications

    E-Print Network [OSTI]

    Burger, S; Pomplun, J; Blome, M; Schmidt, F; 10.1117/12.2001094

    2013-01-01T23:59:59.000Z

    An overview on recent applications of the finite-element method Maxwell-solver JCMsuite to simulation tasks in nanooptics is given. Numerical achievements in the fields of optical metamaterials, plasmonics, photonic crystal fibers, light emitting devices, solar cells, optical lithography, optical metrology, integrated optics, and photonic crystals are summarized.

  9. Modeling Respiratory Lung Motion a Biophysical Approach using Finite Element Methods

    E-Print Network [OSTI]

    Lübeck, Universität zu

    Modeling Respiratory Lung Motion ­ a Biophysical Approach using Finite Element Methods Rene Wernera motion gains in importance. In this paper a biophysical approach for modeling lung motion is described. Main aspects of the process of lung ventilation are identified and outlined as the starting point

  10. Deforming fluid domains within the finite element method: Five mesh-based tracking methods in comparison

    E-Print Network [OSTI]

    Elgeti, Stefanie

    2015-01-01T23:59:59.000Z

    Fluid flow applications can involve a number of coupled problems. One is the simulation of free-surface flows, which require the solution of a free-boundary problem. Within this problem, the governing equations of fluid flow are coupled with a domain deformation approach. This work reviews five of those approaches: interface tracking using a boundary-conforming mesh and, in the interface capturing context, the level-set method, the volume-of-fluid method, particle methods, as well as the phase-field method. The history of each method is presented in combination with the most recent developments in the field. Particularly, the topics of extended finite elements (XFEM) and NURBS-based methods, such as Isogeometric Analysis (IGA), are addressed. For illustration purposes, two applications have been chosen: two-phase flow involving drops or bubbles and sloshing tanks. The challenges of these applications, such as the geometrically correct representation of the free surface or the incorporation of surface tension ...

  11. Prediction of fluid flow in curved pipe using the finite element method

    E-Print Network [OSTI]

    Maitin, Christopher Benjamin

    1987-01-01T23:59:59.000Z

    . Therefore mathematical models have been developed to simulate the elfect of these stresses on the flow field. Until recently these models have only been used for simple geometries. With the advancement of the computer, numerical methods have been de... OF SCIENCE May 1987 Major Subject: Mechanical Engineering PREDICTION OF FLUID FLOW IN CURVED PIPE USING THE FINITE ELEMENT METHOD A Thesis CHRISTOPHER B. MAITIN Approved as to style and content by: Dennis L. O'Neal (Chairman of Committee) Warren...

  12. A study of arbitrary rotations of shells of revolution by the finite element method

    E-Print Network [OSTI]

    Rosenbaum, James Byron

    1974-01-01T23:59:59.000Z

    . Haisler, Jr. Various nonlinear shell theories are reviewed, discussed, and one is selected for implementation into the solution code SNASOR III. The theory selected, Novozhilov's, is discussed in greater detail, comparing his nonlinear equations... to the ones in SNASOR III. The changes required to employ Novozhilov's nonlinear theory in SNASOR III are illustrated. The finite element method is then briefly discussed, nonlinear solution techniques are reviewed and the mod- ified incremental stiffness...

  13. Computation of mode II linear elastic stress intensity factor using the p-version of the finite element method

    SciTech Connect (OSTI)

    Orth, N.J.; Sansguiri, A.M. [Lamar Univ., Beaumont, TX (United States). Dept. of Mechanical Engineering; Orth, F.J. [Edison Welding Inst., Columbus, OH (United States). Structural Integrity Section

    1995-08-01T23:59:59.000Z

    It has been proven that the p-version of the finite element method results in a higher accuracy for a given mesh compared to the h-version of the finite element method. This investigation demonstrates the accuracy of the displacement and stress field methods in the vicinity of the crack tip using the p-version of the finite element method to compute the linear elastic stress intensity factor. The stress intensity factor for mode II is computed for an edge crack panel using stress and displacement methods and is compared to the analytical solution.

  14. Kron`s reduction method applied to the time stepping finite element analysis of induction machines

    SciTech Connect (OSTI)

    Degeneff, R.C.; Gutierrez, M.R.; Salon, S.J.; Burow, D.W. [Rensselaer Polytechnic Inst., Troy, NY (United States). Electric Power Engineering Dept.] [Rensselaer Polytechnic Inst., Troy, NY (United States). Electric Power Engineering Dept.; Nevins, R.J. [Martin Marietta Corp., Schenectady, NY (United States)] [Martin Marietta Corp., Schenectady, NY (United States)

    1995-12-01T23:59:59.000Z

    The behavior of large induction motors during transient as well as steady state running conditions is of significant interest to the power industry. A variety of analytical predictive tools are employed to aid the design and predict their operation under transient and steady state conditions. One of the most powerful method for investigating the transient behavior of induction machines is a coupled time stepping finite element analysis which can combine electromagnetic fields, circuits and mechanical systems. Due to the complexity of the finite element induction machine model and the resulting large number of describing equations, the computation time required for such programs to solve practical problems becomes a major limitation. This becomes even more of a concern when different design options or operating scenarios are evaluated. This paper presents a strategy to reduce the required running time in order to make a parametric study of induction machines such as the assessment of different design options feasible. This is accomplished by reducing the number of finite element equations that must be solved while maintaining the same level of accuracy of solutions. This method is based on Kron`s network reduction work for linear systems and has successfully been applied to large lumped parameter model of transformers. This paper illustrates the reduction method by comparing the flux density in the air gap for a complete FEM model of an induction machine to that of the reduced model. The results are essentially identical with a reduction in computational time of approximately 71%.

  15. HYPERBOLIC DIFFERENTIAL The classical Bubnov-Galerkin finite element method is optimal in the sense of the best ap-

    E-Print Network [OSTI]

    California at Berkeley, University of

    -Galerkin method and subsequently discretize the resulting equations by N + 1 equally-sized finite elementsDRAFT Chapter 9 HYPERBOLIC DIFFERENTIAL EQUATIONS The classical Bubnov-Galerkin finite element equations. In many problems of mechan- ics and convective heat transfer where convection dominates diffusion

  16. Permeability computation on a REV with an immersed finite element method

    SciTech Connect (OSTI)

    Laure, P. [Laboratoire J.-A. Dieudonne, CNRS UMR 6621, Universite de Nice-Sophia Antipolis, Parc Valrose, 06108 Nice, Cedex 02 (France); Puaux, G.; Silva, L.; Vincent, M. [MINES ParisTech, CEMEF-Centre de Mise en Forme des Materiaux, CNRS UMR 7635, BP 207 1 rue Claude, Daunesse 06904 Sophia Antipolis cedex (France)

    2011-05-04T23:59:59.000Z

    An efficient method to compute permeability of fibrous media is presented. An immersed domain approach is used to represent the porous material at its microscopic scale and the flow motion is computed with a stabilized mixed finite element method. Therefore the Stokes equation is solved on the whole domain (including solid part) using a penalty method. The accuracy is controlled by refining the mesh around the solid-fluid interface defined by a level set function. Using homogenisation techniques, the permeability of a representative elementary volume (REV) is computed. The computed permeabilities of regular fibre packings are compared to classical analytical relations found in the bibliography.

  17. Use of the iterative solution method for coupled finite element and boundary element modeling; Yucca Mountain Site Characterization Project

    SciTech Connect (OSTI)

    Koteras, J.R.

    1993-07-01T23:59:59.000Z

    Tunnels buried deep within the earth constitute an important class geomechanics problems. Two numerical techniques used for the analysis of geomechanics problems, the finite element method and the boundary element method, have complementary characteristics for applications to problems of this type. The usefulness of combining these two methods for use as a geomechanics analysis tool has been recognized for some time, and a number of coupling techniques have been proposed. However, not all of them lend themselves to efficient computational implementations for large-scale problems. This report examines a coupling technique that can form the basis for an efficient analysis tool for large scale geomechanics problems through the use of an iterative equation solver.

  18. Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media

    E-Print Network [OSTI]

    Gao, Kai; Gibson, Richard L; Chung, Eric T; Efendiev, Yalchin

    2014-01-01T23:59:59.000Z

    It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both boundaries and the interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and di...

  19. Simulation of borehole-eccentered triaxial induction measurements using a Fourier hp finite-element method

    E-Print Network [OSTI]

    Torres-Verdín, Carlos

    Simulation of borehole-eccentered triaxial induction measurements using a Fourier hp finite-element with a 2D, goal-oriented, high-order, and self- adaptive hp finite-element refinement strategy three mutually orthogonal transmitter coils located at the same vertical position and three collocated

  20. Probabilistic finite element analysis of marine risers

    E-Print Network [OSTI]

    Leder, H. Vern

    1990-01-01T23:59:59.000Z

    Review 1. 2 Research Study 10 2 FORMULATION OF THE SECOND ? MOMENT ANALYSIS METHOD 13 2. 1 Finite Element Equations 2. 2 Random Vector Formulation 2. 3 The Correlation Function 2. 4 Random Field Discretization 2. 5 Taylor Series Expansion 2. 6... ILLUSTRATIVE EXAMPLE. . . . . . . . . . . . . . . . . . 4 APPLICATION OF PROBABILISTIC FINITE ELEMENT METHODS TO MARINE RISER ANALYSES 4. 1 Finite Element Model . 4. 1. 1 Formulation of the Equation of Motion 4. 1. 2 Finite Element Discretization 4. 1. 3...

  1. Least-squares variational principles and the finite element method: theory, formulations, and models for solid and fluid mechanics 

    E-Print Network [OSTI]

    Pontaza, Juan Pablo

    2004-09-30T23:59:59.000Z

    We consider the application of least-squares variational principles and the finite element method to the numerical solution of boundary value problems arising in the fields of solidand fluidmechanics.For manyof these ...

  2. Simulation of the ultrasonic array response from real branched cracks using an efficient finite element method

    SciTech Connect (OSTI)

    Felice, Maria V. [Department of Mechanical Engineering, University of Bristol, Bristol BS8 1TR, United Kingdom and Rolls-Royce plc., Bristol BS34 7QE (United Kingdom); Velichko, Alexander; Wilcox, Paul D. [Department of Mechanical Engineering, University of Bristol, Bristol BS8 1TR (United Kingdom); Barden, Tim J.; Dunhill, Tony K. [Rolls-Royce plc., Bristol BS34 7QE (United Kingdom)

    2014-02-18T23:59:59.000Z

    A hybrid model to simulate the ultrasonic array response from stress corrosion cracks is presented. These cracks are branched and difficult to detect so the model is required to enable optimization of an array design. An efficient frequency-domain finite element method is described and selected to simulate the ultrasonic scattering. Experimental validation results are presented, followed by an example of the simulated ultrasonic array response from a real stress corrosion crack whose geometry is obtained from an X-ray Computed Tomography image. A simulation-assisted array design methodology, which includes the model and use of real crack geometries, is proposed.

  3. Higher-order adaptive finite-element methods for Kohn–Sham density functional theory

    SciTech Connect (OSTI)

    Motamarri, P. [Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109 (United States)] [Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109 (United States); Nowak, M.R. [Department of Electrical Engineering, University of Michigan, Ann Arbor, MI 48109 (United States)] [Department of Electrical Engineering, University of Michigan, Ann Arbor, MI 48109 (United States); Leiter, K.; Knap, J. [U.S. Army Research Labs, Aberdeen Proving Ground, Aberdeen, MD 21001 (United States)] [U.S. Army Research Labs, Aberdeen Proving Ground, Aberdeen, MD 21001 (United States); Gavini, V., E-mail: vikramg@umich.edu [Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI 48109 (United States)

    2013-11-15T23:59:59.000Z

    We present an efficient computational approach to perform real-space electronic structure calculations using an adaptive higher-order finite-element discretization of Kohn–Sham density-functional theory (DFT). To this end, we develop an a priori mesh-adaption technique to construct a close to optimal finite-element discretization of the problem. We further propose an efficient solution strategy for solving the discrete eigenvalue problem by using spectral finite-elements in conjunction with Gauss–Lobatto quadrature, and a Chebyshev acceleration technique for computing the occupied eigenspace. The proposed approach has been observed to provide a staggering 100–200-fold computational advantage over the solution of a generalized eigenvalue problem. Using the proposed solution procedure, we investigate the computational efficiency afforded by higher-order finite-element discretizations of the Kohn–Sham DFT problem. Our studies suggest that staggering computational savings—of the order of 1000-fold—relative to linear finite-elements can be realized, for both all-electron and local pseudopotential calculations, by using higher-order finite-element discretizations. On all the benchmark systems studied, we observe diminishing returns in computational savings beyond the sixth-order for accuracies commensurate with chemical accuracy, suggesting that the hexic spectral-element may be an optimal choice for the finite-element discretization of the Kohn–Sham DFT problem. A comparative study of the computational efficiency of the proposed higher-order finite-element discretizations suggests that the performance of finite-element basis is competing with the plane-wave discretization for non-periodic local pseudopotential calculations, and compares to the Gaussian basis for all-electron calculations to within an order of magnitude. Further, we demonstrate the capability of the proposed approach to compute the electronic structure of a metallic system containing 1688 atoms using modest computational resources, and good scalability of the present implementation up to 192 processors.

  4. Medial Techniques for Automating Finite Element Analysis

    E-Print Network [OSTI]

    Whitton, Mary C.

    Medial Techniques for Automating Finite Element Analysis by Jessica Renee Crawford Crouch Analysis. (Under the direction of Stephen M. Pizer.) Finite element analysis provides a principled method the simulation of tissue deformation. The drawback to using finite element analysis for imaging problems

  5. Finite Element Analysis in Vertebrate Biomechanics

    E-Print Network [OSTI]

    Finite Element Analysis in Vertebrate Biomechanics CALLUM F. ROSS* Department of Organismal Biology presents a series of papers that apply the method of finite element analysis (FEA) to questions ontogenetic or phylogenetic transformations. © 2005 Wiley-Liss, Inc. Key words: finite-element analysis

  6. A stabilized mixed finite element method for Darcy flow Arif Masud a,*, Thomas J.R. Hughes b

    E-Print Network [OSTI]

    Masud, Arif

    , advective­diffusive equations for concentrations are coupled to the Darcy flow equation. For examplesA stabilized mixed finite element method for Darcy flow Arif Masud a,*, Thomas J.R. Hughes b element methods for Darcy flow. Stability and an a priori error estimate in the ``stability norm

  7. Implementation of a python version of a scaled boundary finite element method for plate bending analysis

    E-Print Network [OSTI]

    Chen, Lingfeng, M. Eng. Massachusetts Institute of Technology

    2014-01-01T23:59:59.000Z

    Common finite element programs for plate bending analysis are complicated and limited by the common plate theories. Such programs are usually not user-friendly for designers to implement. Lately, Hou Man et al. from the ...

  8. Discontinuous Galerkin Finite Element Method for the Nonlinear Hyperbolic Problems with Entropy-Based Artificial Viscosity Stabilization

    E-Print Network [OSTI]

    Zingan, Valentin Nikolaevich

    2012-07-16T23:59:59.000Z

    Galerkin DoF degree of freedom ENO essentially non-oscillatory EoS equation of state FE finite element FEM finite element method FVM finite volume method IBVP initial boundary value problem IVP initial value problem JLG Jean-Luc Guermond KPP..., the influence of diffusion, friction, and thermal conduction is limited to a very thin region adjacent to a body surface, a so-called boundary layer, and the remainder of the flow outside this region is essentially inviscid. For flows over slender bodies...

  9. Probabilistic finite element analysis of a craniofacial finite element model Michael A. Berthaume a,b

    E-Print Network [OSTI]

    Probabilistic finite element analysis of a craniofacial finite element model Michael A. Berthaume a s t r a c t We employed a probabilistic finite element analysis (FEA) method to determine how 2011 Accepted 18 January 2012 Available online 27 January 2012 Keywords: Probabilistic analysis Finite

  10. Finite Element Analysis Skateboard Truck

    E-Print Network [OSTI]

    De, Suvranu

    Finite Element Analysis Of a Skateboard Truck #12;2 Executive Summary: Engineering is and always is an element of the `truck,' which holds the wheels. Finite Element analysis will be conducted on this piece a combination of SolidWorks (for modeling) and ABAQUS (for finite element analysis). It is evident from

  11. Process metallurgy simulation for metal drawing process optimization by using two-scale finite element method

    SciTech Connect (OSTI)

    Nakamachi, Eiji; Yoshida, Takashi; Yamaguchi, Toshihiko; Morita, Yusuke [Dept. of Biomedical Engineering, Doshisha Univ., Kyotanabe, Kyoto, 610-0394 (Japan); Kuramae, Hiroyuki [Dept. of Technol. Management, Osaka Institute of Technol., Asahiku, Osaka, 535-8585 (Japan); Morimoto, Hideo [Furukawa Electric Co. Ltd., Okano, Nishiku, Yokohama, Kanagawa, 220-0073 (Japan)

    2014-10-06T23:59:59.000Z

    We developed two-scale FE analysis procedure based on the crystallographic homogenization method by considering the hierarchical structure of poly-crystal aluminium alloy metal. It can be characterized as the combination of two-scale structure, such as the microscopic polycrystal structure and the macroscopic elastic plastic continuum. Micro polycrystal structure can be modeled as a three dimensional representative volume element (RVE). RVE is featured as by 3×3×3 eight-nodes solid finite elements, which has 216 crystal orientations. This FE analysis code can predict the deformation, strain and stress evolutions in the wire drawing processes in the macro- scales, and further the crystal texture and hardening evolutions in the micro-scale. In this study, we analyzed the texture evolution in the wire drawing processes by our two-scale FE analysis code under conditions of various drawing angles of dice. We evaluates the texture evolution in the surface and center regions of the wire cross section, and to clarify the effects of processing conditions on the texture evolution.

  12. A STABILIZED VOLUME-AVERAGING FINITE ELEMENT METHOD FOR FLOW IN POROUS MEDIA AND BINARY ALLOY SOLIDIFICATION SYSTEMS

    E-Print Network [OSTI]

    Zabaras, Nicholas J.

    A STABILIZED VOLUME-AVERAGING FINITE ELEMENT METHOD FOR FLOW IN POROUS MEDIA AND BINARY ALLOY of Mechanical and Aerospace Engineering 188 Frank H. T. Rhodes Hall Cornell University Ithaca, NY 14853 algorithm is presented for the analysis of flow in porous media and in the solidification of binary alloys

  13. Finite element methods for 3D eddy current prob-lems in bounded domains subject to realistic

    E-Print Network [OSTI]

    Rodríguez, Rodolfo

    Finite element methods for 3D eddy current prob- lems in bounded domains subject to realistic techniques to approximate eddy current prob- lems. The focus of the article is on the analysis of weak]), magnetostatics ([59, 60, 80, 81, 109, 110, 114]), eddy current ([1, 2, 3, 6, 7, 8, 9, 19, 20, 21, 23, 24, 25, 26

  14. A stabilized mixed finite element method for Darcy flow based on a multiscale decomposition of the solution

    E-Print Network [OSTI]

    Masud, Arif

    . There are many processes in nat- ure that can be modeled using the Darcy flow equations. Oil extraction to be solenoidal (i.e., a divergence free vector field). The main issue is that the equations of Darcy flow leadA stabilized mixed finite element method for Darcy flow based on a multiscale decomposition

  15. Full wave simulation of waves in ECRIS plasmas based on the finite element method

    SciTech Connect (OSTI)

    Torrisi, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania, Italy and Università Mediterranea di Reggio Calabria, Dipartimento di Ingegneria dell'Informazione, delle Infrastrutture e dell'Energia Sostenibile (DIIES), Via Graziella, I (Italy); Mascali, D.; Neri, L.; Castro, G.; Patti, G.; Celona, L.; Gammino, S.; Ciavola, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania (Italy); Di Donato, L. [Università degli Studi di Catania, Dipartimento di Ingegneria Elettrica Elettronica ed Informatica (DIEEI), Viale Andrea Doria 6, 95125 Catania (Italy); Sorbello, G. [INFN - Laboratori Nazionali del Sud, via S. Sofia 62, 95123, Catania, Italy and Università degli Studi di Catania, Dipartimento di Ingegneria Elettrica Elettronica ed Informatica (DIEEI), Viale Andrea Doria 6, 95125 Catania (Italy); Isernia, T. [Università Mediterranea di Reggio Calabria, Dipartimento di Ingegneria dell'Informazione, delle Infrastrutture e dell'Energia Sostenibile (DIIES), Via Graziella, I-89100 Reggio Calabria (Italy)

    2014-02-12T23:59:59.000Z

    This paper describes the modeling and the full wave numerical simulation of electromagnetic waves propagation and absorption in an anisotropic magnetized plasma filling the resonant cavity of an electron cyclotron resonance ion source (ECRIS). The model assumes inhomogeneous, dispersive and tensorial constitutive relations. Maxwell's equations are solved by the finite element method (FEM), using the COMSOL Multiphysics{sup ®} suite. All the relevant details have been considered in the model, including the non uniform external magnetostatic field used for plasma confinement, the local electron density profile resulting in the full-3D non uniform magnetized plasma complex dielectric tensor. The more accurate plasma simulations clearly show the importance of cavity effect on wave propagation and the effects of a resonant surface. These studies are the pillars for an improved ECRIS plasma modeling, that is mandatory to optimize the ion source output (beam intensity distribution and charge state, especially). Any new project concerning the advanced ECRIS design will take benefit by an adequate modeling of self-consistent wave absorption simulations.

  16. Finite Element Analysis of the Schroedinger Equation

    E-Print Network [OSTI]

    Avtar S. Sehra

    2007-04-17T23:59:59.000Z

    The purpose of this work is to test the application of the finite element method to quantum mechanical problems, in particular for solving the Schroedinger equation. We begin with an overview of quantum mechanics, and standard numerical techniques. We then give an introduction to finite element analysis using the diffusion equation as an example. Three numerical time evolution methods are considered: the (tried and tested) Crank-Nicolson method, the continuous space-time method, and the discontinuous space-time method.

  17. Visualization methods for high-resolution, transient, 3-D, finite element situations

    SciTech Connect (OSTI)

    Christon, M.A.

    1995-01-10T23:59:59.000Z

    Scientific visualization is the process whereby numerical data is transformed into a visual form to augment the process of discovery and understanding. Visualizing the data generated by large-scale, transient, three-dimensional finite element simulations poses many challenges due to geometric complexity, the presence of multiple materials and multiple element types, and the inherent unstructured nature of the meshes. In this paper, the direct use of finite element data structures, nodal assembly procedures, and element interpolants for volumetric adaptive surface extraction, surface rendering, vector grids and particle tracing is discussed. A brief description of a {open_quotes}direct-to-disk{close_quotes} animation system is presented, and case studies which demonstrate the use of isosurfaces, vector plots, cutting planes, reference surfaces and particle tracing are then discussed in the context of several case studies for transient incompressible viscous flow, and acoustic fluid-structure interaction simulations. An overview of the implications of massively parallel computers on visualization is presented to highlight the issues in parallel visualization methodology, algorithms. data locality and the ultimate requirements for temporary and archival data storage and network bandwidth.

  18. Use of the finite element method to predict roof collapse and subsidence resulting from the underground gasification of coal

    SciTech Connect (OSTI)

    Jegbefume, E.U.

    1981-01-01T23:59:59.000Z

    The plane strain two-dimensional finite element model employing the method of negative reaction loading was used to model roof collapse and subsidence that result from underground coal gasification. Subsidence results from a preliminary finite element model were compared with those from a theoretical model having identical assumptions. The finite element model was extended to include dry zones of rock around the cavity resulting from the underground gasification. The model was further extended to include thermal loading, thermal softening and roof collapse. Roof collapse initiation and growth were modelled by examining the portion of the roof immediately above the cavity. A creep program utilizing the method incremental strains was developed and run in conjunction with the Pafec 70 + program. The overburden rock was treated as a linear viscoelastic material having deformation parameters of the generalized Kelvin model. Finally, the finite element model was used to carry out a case history study of subsidence and roof collapse of an underground gasification field experiment in Wyoming.

  19. Scalability Analysis of Partitioning Strategies for Finite Element Graphs. \\Lambda

    E-Print Network [OSTI]

    Kumar, Vipin

    Scalability Analysis of Partitioning Strategies for Finite Element Graphs. \\Lambda Grama Y. Ananth Y. Grama ananth@cs:umn:edu Abstract Issues of partitioning Finite Element Graphs are central for implementing the Finite Element Method (FEM). Parallel formulations of finite element techniques require

  20. DOLFIN: Automated Finite Element Computing

    E-Print Network [OSTI]

    Logg, Anders; Wells, G N

    2009-01-01T23:59:59.000Z

    ´de´lec 1980]. (4) L2-conforming finite elements: (a) DGq, arbitrary degree discontinuous Lagrange elements; and (b) CR1, first degree Crouzeix–Raviart5 elements [Crouzeix and Raviart 1973]. Arbitrary combinations of the above elements may be used to define...

  1. electroseismic monitoring of co2 sequestration: a finite element ...

    E-Print Network [OSTI]

    Fabio Zyserman

    Keywords: Electroseismic Modeling, Poroelasticity, CO2 sequestration, Finite element methods. 2000 AMS ... carbon dioxide emissisons into the atmosphere.

  2. Numerical Stochastic Homogenization Method and Multiscale Stochastic Finite Element Method - A Paradigm for Multiscale Computation of Stochastic PDEs

    SciTech Connect (OSTI)

    X. Frank Xu

    2010-03-30T23:59:59.000Z

    Multiscale modeling of stochastic systems, or uncertainty quantization of multiscale modeling is becoming an emerging research frontier, with rapidly growing engineering applications in nanotechnology, biotechnology, advanced materials, and geo-systems, etc. While tremendous efforts have been devoted to either stochastic methods or multiscale methods, little combined work had been done on integration of multiscale and stochastic methods, and there was no method formally available to tackle multiscale problems involving uncertainties. By developing an innovative Multiscale Stochastic Finite Element Method (MSFEM), this research has made a ground-breaking contribution to the emerging field of Multiscale Stochastic Modeling (MSM) (Fig 1). The theory of MSFEM basically decomposes a boundary value problem of random microstructure into a slow scale deterministic problem and a fast scale stochastic one. The slow scale problem corresponds to common engineering modeling practices where fine-scale microstructure is approximated by certain effective constitutive constants, which can be solved by using standard numerical solvers. The fast scale problem evaluates fluctuations of local quantities due to random microstructure, which is important for scale-coupling systems and particularly those involving failure mechanisms. The Green-function-based fast-scale solver developed in this research overcomes the curse-of-dimensionality commonly met in conventional approaches, by proposing a random field-based orthogonal expansion approach. The MSFEM formulated in this project paves the way to deliver the first computational tool/software on uncertainty quantification of multiscale systems. The applications of MSFEM on engineering problems will directly enhance our modeling capability on materials science (composite materials, nanostructures), geophysics (porous media, earthquake), biological systems (biological tissues, bones, protein folding). Continuous development of MSFEM will further contribute to the establishment of Multiscale Stochastic Modeling strategy, and thereby potentially to bring paradigm-shifting changes to simulation and modeling of complex systems cutting across multidisciplinary fields.

  3. THE MORTAR FINITE ELEMENT METHOD FOR 3D MAXWELL EQUATIONS: FIRST RESULTS

    E-Print Network [OSTI]

    Buffa, Annalisa

    of the mortar element method to the numerical solution of Maxwell's equations (in wave form) by H; accepted for publication (in revised form) February 2, 2001; published electronically July 31, 2001. This work has been partially supported by Istituto di Analisi Numerica del Centro Nazionale delle Ricerche

  4. A multiscale preconditioner for stochastic mortar mixed finite elements

    E-Print Network [OSTI]

    Yotov, Ivan

    physical models in different subdomains. The mortar finite element method is a generalization is the stochastic collocation method [4,40,39]. It combines a finite element discretiza- tion in physical space the approximation properties of the stochastic finite element method [5,37], making it more efficient than MCS

  5. PHYSICAL REVIEW E 85, 036706 (2012) Finite-size scaling for quantum criticality using the finite-element method

    E-Print Network [OSTI]

    Kais, Sabre

    initio approximations and exact formulations. The critical parameters could be atomic nuclear charges, and an "exact" formulation using FEM. We then use the FSS approach to determine its critical nuclear chargePHYSICAL REVIEW E 85, 036706 (2012) Finite-size scaling for quantum criticality using the finite

  6. Finite Element Analysis in Functional BRIAN G. RICHMOND,1

    E-Print Network [OSTI]

    Finite Element Analysis in Functional Morphology BRIAN G. RICHMOND,1 * BARTH W. WRIGHT,1 IAN GROSSE element analysis. © 2005 Wiley-Liss, Inc. Key words: finite-element analysis; mastication; primates; biome This article reviews the fundamental principles of the finite element method and the three basic steps (model

  7. Flexible finite-element modeling of global geomagnetic depth sounding

    E-Print Network [OSTI]

    Ribaudo, Joseph Thomas

    2011-01-01T23:59:59.000Z

    OF CALIFORNIA, SAN DIEGO Flexible Finite-Element Modeling ofOF THE DISSERTATION Flexible Finite-Element Modeling ofin Chapter 2. The method is flexible enough to allow for a

  8. Numerical simulation of two-phase flow in discrete fractures using Rayleigh-Ritz finite element method

    E-Print Network [OSTI]

    Kaul, Sandeep P.

    2004-09-30T23:59:59.000Z

    Spontaneous imbibition plays a very important role in the displacement mechanism of non-wetting fluid in naturally fractured reservoirs. We developed a new 2D two-phase finite element numerical model, as available commercial simulators cannot...

  9. A Full-wave Model for Wave Propagation and Dissipation in the Inner Magnetosphere Using the Finite Element Method

    SciTech Connect (OSTI)

    Ernest Valeo, Jay R. Johnson, Eun-Hwa and Cynthia Phillips

    2012-03-13T23:59:59.000Z

    A wide variety of plasma waves play an important role in the energization and loss of particles in the inner magnetosphere. Our ability to understand and model wave-particle interactions in this region requires improved knowledge of the spatial distribution and properties of these waves as well as improved understanding of how the waves depend on changes in solar wind forcing and/or geomagnetic activity. To this end, we have developed a two-dimensional, finite element code that solves the full wave equations in global magnetospheric geometry. The code describes three-dimensional wave structure including mode conversion when ULF, EMIC, and whistler waves are launched in a two-dimensional axisymmetric background plasma with general magnetic field topology. We illustrate the capabilities of the code by examining the role of plasmaspheric plumes on magnetosonic wave propagation; mode conversion at the ion-ion and Alfven resonances resulting from external, solar wind compressions; and wave structure and mode conversion of electromagnetic ion cyclotron waves launched in the equatorial magnetosphere, which propagate along the magnetic field lines toward the ionosphere. We also discuss advantages of the finite element method for resolving resonant structures, and how the model may be adapted to include nonlocal kinetic effects.

  10. Partition-of-unity finite-element method for large scale quantum molecular dynamics on massively parallel computational platforms

    SciTech Connect (OSTI)

    Pask, J E; Sukumar, N; Guney, M; Hu, W

    2011-02-28T23:59:59.000Z

    Over the course of the past two decades, quantum mechanical calculations have emerged as a key component of modern materials research. However, the solution of the required quantum mechanical equations is a formidable task and this has severely limited the range of materials systems which can be investigated by such accurate, quantum mechanical means. The current state of the art for large-scale quantum simulations is the planewave (PW) method, as implemented in now ubiquitous VASP, ABINIT, and QBox codes, among many others. However, since the PW method uses a global Fourier basis, with strictly uniform resolution at all points in space, and in which every basis function overlaps every other at every point, it suffers from substantial inefficiencies in calculations involving atoms with localized states, such as first-row and transition-metal atoms, and requires substantial nonlocal communications in parallel implementations, placing critical limits on scalability. In recent years, real-space methods such as finite-differences (FD) and finite-elements (FE) have been developed to address these deficiencies by reformulating the required quantum mechanical equations in a strictly local representation. However, while addressing both resolution and parallel-communications problems, such local real-space approaches have been plagued by one key disadvantage relative to planewaves: excessive degrees of freedom (grid points, basis functions) needed to achieve the required accuracies. And so, despite critical limitations, the PW method remains the standard today. In this work, we show for the first time that this key remaining disadvantage of real-space methods can in fact be overcome: by building known atomic physics into the solution process using modern partition-of-unity (PU) techniques in finite element analysis. Indeed, our results show order-of-magnitude reductions in basis size relative to state-of-the-art planewave based methods. The method developed here is completely general, applicable to any crystal symmetry and to both metals and insulators alike. We have developed and implemented a full self-consistent Kohn-Sham method, including both total energies and forces for molecular dynamics, and developed a full MPI parallel implementation for large-scale calculations. We have applied the method to the gamut of physical systems, from simple insulating systems with light atoms to complex d- and f-electron systems, requiring large numbers of atomic-orbital enrichments. In every case, the new PU FE method attained the required accuracies with substantially fewer degrees of freedom, typically by an order of magnitude or more, than the current state-of-the-art PW method. Finally, our initial MPI implementation has shown excellent parallel scaling of the most time-critical parts of the code up to 1728 processors, with clear indications of what will be required to achieve comparable scaling for the rest. Having shown that the key remaining disadvantage of real-space methods can in fact be overcome, the work has attracted significant attention: with sixteen invited talks, both domestic and international, so far; two papers published and another in preparation; and three new university and/or national laboratory collaborations, securing external funding to pursue a number of related research directions. Having demonstrated the proof of principle, work now centers on the necessary extensions and optimizations required to bring the prototype method and code delivered here to production applications.

  11. A Novel High Order Time Domain Vector Finite Element Method for the Simulation of Electromagnetic Devices

    SciTech Connect (OSTI)

    Rieben, R N

    2004-07-20T23:59:59.000Z

    The goal of this dissertation is twofold. The first part concerns the development of a numerical method for solving Maxwell's equations on unstructured hexahedral grids that employs both high order spatial and high order temporal discretizations. The second part involves the use of this method as a computational tool to perform high fidelity simulations of various electromagnetic devices such as optical transmission lines and photonic crystal structures to yield a level of accuracy that has previously been computationally cost prohibitive. This work is based on the initial research of Daniel White who developed a provably stable, charge and energy conserving method for solving Maxwell's equations in the time domain that is second order accurate in both space and time. The research presented here has involved the generalization of this procedure to higher order methods. High order methods are capable of yielding far more accurate numerical results for certain problems when compared to corresponding h-refined first order methods , and often times at a significant reduction in total computational cost. The first half of this dissertation presents the method as well as the necessary mathematics required for its derivation. The second half addresses the implementation of the method in a parallel computational environment, its validation using benchmark problems, and finally its use in large scale numerical simulations of electromagnetic transmission devices.

  12. PERFORMANCE OF FULLY-COUPLED DOMAIN DECOMPOSITION PRECONDITIONERS for FINITE ELEMENT TRANSPORT / REACTION SIMULATIONS1

    E-Print Network [OSTI]

    Tuminaro, Ray S.

    element simulation of strongly coupled fluid flow, heat transfer, and mass transfer with non methods, stabilized finite element methods. 1. INTRODUCTION Modern computational fluid dynamics flows with unstructured mesh finite element methods. These flow problems are characterized by both

  13. A 3D finite element method for flexible multibody Johannes Gerstmayr, Joachim Schoberl

    E-Print Network [OSTI]

    Schoeberl, Joachim

    -blade of a helicopter or the stability of a satellite with highly flexible solar panels, see Bremer [5]. While most methods. This is applicable for the case for several engineering applications, like a gear unit

  14. A performance study of plane wave finite element methods with a Pad-type artificial boundary condition in acoustic scattering

    E-Print Network [OSTI]

    Antoine, Xavier - Institut de Mathématiques �lie Cartan, Université Henri Poincaré

    Keywords: Acoustic scattering Finite element Plane wave Pollution High-frequency Artificial boundary condition in acoustic scattering R. Kechroud a , A. Soulaimani a , X. Antoine b,c,* a �cole de Technologie are found. Indeed, it is well known that pollution effects [10­12] arise in the standard FEM, limiting its

  15. Hybrid asymptotic-finite element method for stiff two-point boundary value problems

    SciTech Connect (OSTI)

    Chin, R.C.Y.; Krasny, R.

    1983-06-01T23:59:59.000Z

    An accurate and efficient numerical method has been developed for a nonlinear stiff second order two-point boundary value problem. The scheme combines asymptotic methods with the usual solution techniques for two-point boundary value problems. A new modification of Newton's method or quasilinearization is used to reduce the nonlinear problem to a sequence of linear problems. The resultant linear problem is solved by patching local solutions at the knots or equivalently by projecting onto an affine subset constructed from asymptotic expansions. In this way, boundary layers are naturally incorporated into the approximation. An adaptive mesh is employed to achieve an error of O(1/N/sup 2/) + O(..sqrt..e). Here, N is the number of intervals and epsilon << 1 is the singular perturbation parameter. Numerical computations are presented.

  16. Submitted to: Finite Elements in Analysis and Design Finite Element Analysis over Tangled Simplicial Meshes

    E-Print Network [OSTI]

    Suresh, Krishnan

    Submitted to: Finite Elements in Analysis and Design 1 Finite Element Analysis over Tangled In modern finite element analysis (FEA), a mesh is said to be `tangled' if it contains one or more inverted are also illustrated. 1. INTRODUCTION In modern finite element analysis (FEA), the underlying mesh

  17. Finite element calculations illustrating a method of model reduction for the dynamics of structures with localized nonlinearities.

    SciTech Connect (OSTI)

    Griffith, Daniel Todd; Segalman, Daniel Joseph

    2006-10-01T23:59:59.000Z

    A technique published in SAND Report 2006-1789 ''Model Reduction of Systems with Localized Nonlinearities'' is illustrated in two problems of finite element structural dynamics. That technique, called here the Method of Locally Discontinuous Basis Vectors (LDBV), was devised to address the peculiar difficulties of model reduction of systems having spatially localized nonlinearities. It's illustration here is on two problems of different geometric and dynamic complexity, but each containing localized interface nonlinearities represented by constitutive models for bolted joint behavior. As illustrated on simple problems in the earlier SAND report, the LDBV Method not only affords reduction in size of the nonlinear systems of equations that must be solved, but it also facilitates the use of much larger time steps on problems of joint macro-slip than would be possible otherwise. These benefits are more dramatic for the larger problems illustrated here. The work of both the original SAND report and this one were funded by the LDRD program at Sandia National Laboratories.

  18. Time finite element methods for large rotational dynamics of multibody systems

    SciTech Connect (OSTI)

    Mello, F.J.; Borri, M.; Atluri, S.N. (Georgia Inst. of Tech., Atlanta, GA (USA). Center for Computational Mechanics)

    1990-01-01T23:59:59.000Z

    Weak formulations in Analytical Dynamics are developed, paralleling the variational methods in elastostatics, and including a fundamental yet novel approach for treating constraints (both holonomic and nonholonomic). A general three field approach is presented, in which the momentum balance conditions, the compatibility conditions between displacement and velocity, the constitutive relations and the displacement and momentum boundary conditions are all enforced in weak form. A primal, or kinematic formulation is developed from the general form by enforcing the compatibility conditions and displacement boundary conditions a priori. The conditional stability of the kinematic formulation is the counterpart of the locking phenomenon in elastostatics and may be avoided, either by reduced order integration, or by utilizing a mixed formulation. Toward this end, a two field mixed formulation is presented, which follows from the general form, when the constitutive relations are satisfied a priori. A general set of the constraint equations are introduced into the kinematic and mixed formulations, using a specific choice of multipliers, which results in modified variational principles. Several simple examples concerning rigid body dynamics are presented. 15 refs., 18 figs.

  19. Student Award Competition Finite Element EEG and MEG Simulations for Realistic Head Models

    E-Print Network [OSTI]

    Zhukov, Leonid

    than would be possible using linear elements alone. [1] D.S. Burnett. Finite Element Analysis: FromStudent Award Competition Finite Element EEG and MEG Simulations for Realistic Head Models be approximated by a technique such as the finite element method [1]. Most applications of the finite element

  20. Uzawa conjugate gradient method for the Stokes problem: Matlab implementation with P1-iso-P2/P1 finite element

    E-Print Network [OSTI]

    Frey, Pascal

    Uzawa conjugate gradient method for the Stokes problem: Matlab implementation with P1-iso-P2/P1 implementation of the Uzawa Conjugate gradient algorithm for the generalized Stokes problem with a P1-iso-P2/P1 for the generalized Stokes problem using the P1-iso-P2/P1 finite element. The Uzawa scheme is a decomposition

  1. Construction of Higher Order Finite Element with

    E-Print Network [OSTI]

    Kern, Michel

    ' & $ % Construction of Higher Order Finite Element with Mass Lumping Using Computer Algebra. (3D, combinatorial analysis, new third order element) 2 #12; ' & $ % Guidelines for the construction of nodes must be ~ P k unisolvent. 2. Finite element must be continuous. 3. Quadrature formula must satisfy

  2. CCM Continuity Constraint Method: A finite-element computational fluid dynamics algorithm for incompressible Navier-Stokes fluid flows

    SciTech Connect (OSTI)

    Williams, P.T.

    1993-09-01T23:59:59.000Z

    As the field of computational fluid dynamics (CFD) continues to mature, algorithms are required to exploit the most recent advances in approximation theory, numerical mathematics, computing architectures, and hardware. Meeting this requirement is particularly challenging in incompressible fluid mechanics, where primitive-variable CFD formulations that are robust, while also accurate and efficient in three dimensions, remain an elusive goal. This dissertation asserts that one key to accomplishing this goal is recognition of the dual role assumed by the pressure, i.e., a mechanism for instantaneously enforcing conservation of mass and a force in the mechanical balance law for conservation of momentum. Proving this assertion has motivated the development of a new, primitive-variable, incompressible, CFD algorithm called the Continuity Constraint Method (CCM). The theoretical basis for the CCM consists of a finite-element spatial semi-discretization of a Galerkin weak statement, equal-order interpolation for all state-variables, a 0-implicit time-integration scheme, and a quasi-Newton iterative procedure extended by a Taylor Weak Statement (TWS) formulation for dispersion error control. Original contributions to algorithmic theory include: (a) formulation of the unsteady evolution of the divergence error, (b) investigation of the role of non-smoothness in the discretized continuity-constraint function, (c) development of a uniformly H{sup 1} Galerkin weak statement for the Reynolds-averaged Navier-Stokes pressure Poisson equation, (d) derivation of physically and numerically well-posed boundary conditions, and (e) investigation of sparse data structures and iterative methods for solving the matrix algebra statements generated by the algorithm.

  3. Finite element analysis of human joints

    SciTech Connect (OSTI)

    Bossart, P.L.; Hollerbach, K.

    1996-09-01T23:59:59.000Z

    Our work focuses on the development of finite element models (FEMs) that describe the biomechanics of human joints. Finite element modeling is becoming a standard tool in industrial applications. In highly complex problems such as those found in biomechanics research, however, the full potential of FEMs is just beginning to be explored, due to the absence of precise, high resolution medical data and the difficulties encountered in converting these enormous datasets into a form that is usable in FEMs. With increasing computing speed and memory available, it is now feasible to address these challenges. We address the first by acquiring data with a high resolution C-ray CT scanner and the latter by developing semi-automated method for generating the volumetric meshes used in the FEM. Issues related to tomographic reconstruction, volume segmentation, the use of extracted surfaces to generate volumetric hexahedral meshes, and applications of the FEM are described.

  4. A PROBABILISTIC FINITE ELEMENT ANALYSIS OF RESIDUAL STRESS FORMATION

    E-Print Network [OSTI]

    Grujicic, Mica

    A PROBABILISTIC FINITE ELEMENT ANALYSIS OF RESIDUAL STRESS FORMATION IN SHRINK-FIT CERAMIC shrink fitting of the jacket over the lining is studied using a probabilistic finite element analysis structural analysis approach, known as the Advanced Mean Value (AMV) method, is used which enables

  5. Finite element approximation of coupled seismic and ...

    E-Print Network [OSTI]

    zyserman

    layer, having a thickness of about 10 nm. Finite element approximation of coupled seismic and electromagnetic waves in gas hydrate-bearing sediments – p.

  6. Finite Element Modeling of Drilling Using DEFORM

    E-Print Network [OSTI]

    Gardner, Joel D.; Dornfeld, David

    2006-01-01T23:59:59.000Z

    Vijayaraghavan, A. (2005), “Drilling of Fiber- ReinforcedFINITE ELEMENT MODELING OF DRILLING USING DEFORM J. Gardner,of Comprehensive Drilling Simulation Tool” ABSTRACT DEFORM-

  7. Smooth exterior complex-scaling, full-angular-momentum, and three-dimensional finite-element method applied to doubly excited states of helium

    SciTech Connect (OSTI)

    Elander, Nils [Department of Physics, Stockholm University, Alba Nova University Center, SE-106 91 Stockholm (Sweden); Levin, Sergey [Department of Physics, Stockholm University, Alba Nova University Center, SE-106 91 Stockholm (Sweden); Institute for Physics, St. Petersburg University, Uljanovskaya 1, St. Petersburg 198904 (Russian Federation); Institute for Theoretical Atomic and Molecular Physics at Harvard University and Smithsonian Astrophysical Observatory, 60 Garden Street, MS 14 Cambridge, Massachusetts 02138 (United States); Yarevsky, Evgeny [Institute for Physics, St. Petersburg University, Uljanovskaya 1, St. Petersburg 198904 (Russian Federation); International Solvay Institutes for Physics and Chemistry, Case Postale 231, Campus Plaine ULB Boulevard du Triomphe, Brussels 1050 (Belgium)

    2003-06-01T23:59:59.000Z

    A technique based on the total-angular-momentum representation the smooth exterior complex-scaling procedure, and the three-dimensional finite-element method, is applied to calculations of S, P, and D resonant, so-called doubly excited, states of the helium atom. The resonances are calculated with an accuracy better than 10{sup -5} a.u. The applicability of an extrapolation procedure to complex energies is analyzed.

  8. Probabilistic slope stability analysis by finite elements D.V. Griffiths*

    E-Print Network [OSTI]

    Probabilistic slope stability analysis by finite elements D.V. Griffiths* and Gordon A. Fenton. The advanced method, called the random finite element method (RFEM), uses elasto- plasticity combined nonlinear finite element methods are combined with ran- dom field generation techniques. This method, called

  9. Optimal convergence analysis for the eXtended Finite Element Serge Nicaise 1

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Optimal convergence analysis for the eXtended Finite Element Method Serge Nicaise 1 , Yves Renard 2Xtended Finite Element Method (Xfem), namely the Xfem with a cut-off function and the stan- dard Xfem finite element method, error estimates, stress intensity factors. 1 Introduction Inspired by the Pufem

  10. Benchmarking of finite element codes for radioactive material transportation packages

    SciTech Connect (OSTI)

    Ammerman, D.J.

    1996-10-01T23:59:59.000Z

    The increased power of computers and computer codes makes the use of nonlinear dynamic finite element analyses attractive for use as a tool used in the design and certification of radioactive material transportation packages. For this analysis technique to be acceptable it must be demonstrated. The technique has the ability to accurately capture the response of the packages to accident environments required by the regulations. The best method of demonstrating this ability is via a series of benchmark analyses. In this paper three benchmark problems involving significant inelastic deformations will be discussed. One of the problems has been analyzed using many different finite element codes. The other two problems involve comparison of finite element calculations to the results form physical tests. The ability of the finite element method to accurately capture the response in these three problems indicates the method should be acceptable for radioactive material transportation package design and certification.

  11. Finite element analysis of multilayer coextrusion.

    SciTech Connect (OSTI)

    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-01T23:59:59.000Z

    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.

  12. Application of Mass Lumped Higher Order Finite Elements

    SciTech Connect (OSTI)

    Chen, J.; Strauss, H. R.; Jardin, S. C.; Park, W.; Sugiyama, L. E.; Fu., G.; Breslau, J.

    2005-11-01T23:59:59.000Z

    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.

  13. adaptive finite element: Topics by E-print Network

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    properties and performance of the adaptive finite element approach to the design of nano-photonic components. Central issues are the construction of vectorial finite elements...

  14. adaptive finite elements: Topics by E-print Network

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    properties and performance of the adaptive finite element approach to the design of nano-photonic components. Central issues are the construction of vectorial finite elements...

  15. The Efficient Computation of Bounds for Functionals of Finite Element Solutions in Large Strain Elasticity

    E-Print Network [OSTI]

    Bonet, J.

    We present an implicit a-posteriori finite element procedure to compute bounds for functional outputs of finite element solutions in large strain elasticity. The method proposed relies on the existence of a potential energy ...

  16. 2.094 Finite Element Analysis of Solids and Fluids, Spring 2008

    E-Print Network [OSTI]

    Kim, Do-Nyun

    This course presents finite element theory and methods for general linear and nonlinear analyses. Reliable and effective finite element procedures are discussed with their applications to the solution of general problems ...

  17. Finite element decomposition of the human neocortex

    E-Print Network [OSTI]

    Chow, Seeling

    1998-01-01T23:59:59.000Z

    Relationships between iso-parametric curves and principal curvature directions, . . . . . . . . . . . . . . . . . . . . . . . . . . . . . , . ?. ???, . ?. . . . . . . . . . ??. ?. 25 Reparameterizing macro elements for finite element decomposition... visualization and functional analysis, building 3D surface models through serial reconstruction has received much attention in biomedical research [3], [8], [50]. Surface reconstruction techniques include volume rendering, iso-surface algorithms, parametric...

  18. Modeling Elastic Properties in Finite-Element Analysis: How Much

    E-Print Network [OSTI]

    Modeling Elastic Properties in Finite- Element Analysis: How Much Precision Is Needed to Produce analysis was investigated using a finite-element model of a Macaca fascicularis skull. Four finite-element realistically using the orthotropic elastic properties employed in analysis 4. Results suggest that finite-element

  19. Techniques for Modeling Muscle-Induced Forces in Finite Element

    E-Print Network [OSTI]

    Dumont, Elizabeth R.

    program was written to interface with a commercial finite element analysis tool to automatically apply: finite element analysis; muscle force; skull; muscle loading alogorithm, biting Finite element analysisTechniques for Modeling Muscle- Induced Forces in Finite Element Models of Skeletal Structures IAN

  20. Phase Reduction Models for Improving the Accuracy of the Finite Element Solution of

    E-Print Network [OSTI]

    Boyer, Edmond

    , acoustic scattering, short-wave problem, finite element method, pollution, accuracy, On-Surface Radiation in view of a numerical finite element solution. It is well-known that pollution error causes inaccuracies the pollution error in the finite element solution of time-harmonic scattering problems. The pollution error

  1. JOURNAL OF MATERIALS SCIENCE 38 (2003) 307 322 Crystal plasticity-based finite element analysis

    E-Print Network [OSTI]

    Grujicic, Mica

    2003-01-01T23:59:59.000Z

    JOURNAL OF MATERIALS SCIENCE 38 (2003) 307­ 322 Crystal plasticity-based finite element analysis -TiAl + 2-Ti3Al alloys has been analyzed using a finite element method. A three-dimensional rate interfaces and colony boundaries are implemented in the commercial finite element program Abaqus

  2. Nonrigid Motion Analysis Based on Dynamic Refinement of Finite Element Models

    E-Print Network [OSTI]

    Sarkar, Sudeep

    Nonrigid Motion Analysis Based on Dynamic Refinement of Finite Element Models Leonid V. Tsap finite element models. The method is based on the iterative analysis of the differences betweenÐPhysically-based vision, deformable models, nonrigid motion analysis, biomedical applications, finite element analysis. æ

  3. Finite Elements in Analysis and Design 42 (2006) 580591 www.elsevier.com/locate/finel

    E-Print Network [OSTI]

    Li, Shaofan

    2006-01-01T23:59:59.000Z

    Finite Elements in Analysis and Design 42 (2006) 580­591 www numerical Green's functions, for example bubble functions, #12;X. Liu, S. Li / Finite Elements in Analysis.elsevier.com/locate/finel A variational multiscale stabilized finite element method for the Stokes flow problem Xiaohu Liu , Shaofan Li

  4. A cracked beam finite element for rotating shaft dynamics and stability analysis

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    A cracked beam finite element for rotating shaft dynamics and stability analysis Saber El Arem Palaiseau, France Abstract In this paper, a method for the construction of a cracked beam finite element is presented. The additional flexibility due to the cracks is identified from three- dimensional finite element

  5. Finite Elements in Analysis and Design 43 (2007) 397410 www.elsevier.com/locate/finel

    E-Print Network [OSTI]

    Ghosh, Somnath

    2007-01-01T23:59:59.000Z

    Finite Elements in Analysis and Design 43 (2007) 397­410 www.elsevier.com/locate/finel Modeling introduces an extended Voronoi cell finite element model (X-VCFEM) for modeling the initiation: The extended Voronoi cell finite element method (X-VCFEM); Cohesive zone models; Interfacial debonding; Matrix

  6. Quadrilateral/hexahedral finite element mesh coarsening

    DOE Patents [OSTI]

    Staten, Matthew L; Dewey, Mark W; Scott, Michael A; Benzley, Steven E

    2012-10-16T23:59:59.000Z

    A technique for coarsening a finite element mesh ("FEM") is described. This technique includes identifying a coarsening region within the FEM to be coarsened. Perimeter chords running along perimeter boundaries of the coarsening region are identified. The perimeter chords are redirected to create an adaptive chord separating the coarsening region from a remainder of the FEM. The adaptive chord runs through mesh elements residing along the perimeter boundaries of the coarsening region. The adaptive chord is then extracted to coarsen the FEM.

  7. Optimal control of drug delivery to brain tumors for a PDE driven model using the Galerkin finite element method

    E-Print Network [OSTI]

    Hanson, Floyd B.

    Optimal control of drug delivery to brain tumors for a PDE driven model using the Galerkin finite is used to examine the optimal drug delivery to brain tumors. The PDE driven mathematical model and the drug concentration. An optimal control problem is formulated keeping in mind the primary goals

  8. A New Parallel Domain Decomposition Method for the Adaptive Finite Element Solution of Elliptic Partial Di erential Equations

    E-Print Network [OSTI]

    Bank, Randolph E.

    Partial Di#11;erential Equations Randolph E. Bank #3; and Peter K. Jimack y Abstract We present a new domain decomposition algorithm for the parallel #12;nite element solution of elliptic partial di#11 to be utilized in a well load-balanced manner. Finally, numerical evidence is presented which suggests

  9. Component mode synthesis methods applied to 3D heterogeneous core calculations, using the mixed dual finite element solver MINOS

    SciTech Connect (OSTI)

    Guerin, P.; Baudron, A. M.; Lautard, J. J. [Commissariat a l'Energie Atomique, DEN/DANS/DM2S/SERMA/LENR, CEA Saclay, 91191 Gif sur Yvette (France)

    2006-07-01T23:59:59.000Z

    This paper describes a new technique for determining the pin power in heterogeneous core calculations. It is based on a domain decomposition with overlapping sub-domains and a component mode synthesis technique for the global flux determination. Local basis functions are used to span a discrete space that allows fundamental global mode approximation through a Galerkin technique. Two approaches are given to obtain these local basis functions: in the first one (Component Mode Synthesis method), the first few spatial eigenfunctions are computed on each sub-domain, using periodic boundary conditions. In the second one (Factorized Component Mode Synthesis method), only the fundamental mode is computed, and we use a factorization principle for the flux in order to replace the higher order Eigenmodes. These different local spatial functions are extended to the global domain by defining them as zero outside the sub-domain. These methods are well-fitted for heterogeneous core calculations because the spatial interface modes are taken into account in the domain decomposition. Although these methods could be applied to higher order angular approximations - particularly easily to a SPN approximation - the numerical results we provide are obtained using a diffusion model. We show the methods' accuracy for reactor cores loaded with UOX and MOX assemblies, for which standard reconstruction techniques are known to perform poorly. Furthermore, we show that our methods are highly and easily parallelizable. (authors)

  10. Analysis of PostProcessing for Nonconforming Finite Element Solutions

    E-Print Network [OSTI]

    Schieweck, Friedhelm

    Analysis of Post­Processing for Nonconforming Finite Element Solutions F. Schieweck \\Lambda of a conforming finite element function v C h 2 V C h . One ``nice'' feature of our analysis­processing procedure where a conforming approximation is computed from a nonconforming finite element solution

  11. STATIC ANALYSIS OF SANDWICH PLATES BY FINITE ELEMENTS

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    STATIC ANALYSIS OF SANDWICH PLATES BY FINITE ELEMENTS V. MANET, W.-S. HAN AND A. VAUTRIN ´Ecole des continuity at hal-00663233,version1-26Jan2012 #12;Static analysis of sandwich plates by finite elements 3 and are presented in this paper. 2. Hybrid sandwich finite elements In this section, we present the development

  12. Modeling Masticatory Muscle Force in Finite Element Analysis

    E-Print Network [OSTI]

    Modeling Masticatory Muscle Force in Finite Element Analysis: Sensitivity Analysis Using Principal: electromyography; muscle force; mastication; primates; principal coordinates analysis; finite element analysis Our Origins, University of Arizona, Tempe, Arizona ABSTRACT Our work on a finite element model of the skull

  13. Medially Based Meshing with Finite Element Analysis of Prostate Deformation

    E-Print Network [OSTI]

    Crouch, Jessica R.

    Medially Based Meshing with Finite Element Analysis of Prostate Deformation Jessica R. Crouch1 have approached prostate imaging problems using meth- ods that incorporate finite element analysis. Assume the prostate is a linearly elastic body and compute its deformation using finite element analysis

  14. Automated Finite Element Analysis for Deformable Registration of Prostate Images

    E-Print Network [OSTI]

    1 Automated Finite Element Analysis for Deformable Registration of Prostate Images Jessica R, Marco Zaider Abstract-- Two major factors preventing the routine clinical use of finite element analysis the effort required to apply finite element analysis to image registration. Encouraging results are presented

  15. FINITE ELEMENT ANALYSIS OF A CLASS OF STRESSFREE MARTENSITIC MICROSTRUCTURES

    E-Print Network [OSTI]

    Li, Bo

    FINITE ELEMENT ANALYSIS OF A CLASS OF STRESS­FREE MARTENSITIC MICROSTRUCTURES BO LI Abstract. This work is concerned with the finite element approximation of a class of stress­free martensitic microstructures modeled by multi­well en­ ergy minimization. Finite element energy­minimizing sequences are first

  16. A finite element analysis of pneumatic-tire/sand interactions

    E-Print Network [OSTI]

    Grujicic, Mica

    A finite element analysis of pneumatic-tire/sand interactions during off-road vehicle travel M pneumatic tire and sand during off-road vehicle travel. Keywords Finite element analysis, Road vehicles and for other tire/sand combinations. Since the finite element analysis of the tire/sand interaction enables

  17. Static Analysis of Sandwich Plates by Hybrid Finite Elements

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Static Analysis of Sandwich Plates by Hybrid Finite Elements V. Manet and W.-S. Han Materials analysis tool. The subject of the present work concerns the development of hybrid sandwich finite elements into acount. The present work concerns the development of hybrid sandwich finite elements mod- elling

  18. Building Finite Element Analysis Programs in Distributed Services Environment

    E-Print Network [OSTI]

    Stanford University

    1. Building Finite Element Analysis Programs in Distributed Services Environment Jun Peng1 and Kincho H. Law2 Abstract Traditional finite element analysis (FEA) programs are typically built on a dedicated computer using the developments offered by a finite element analysis (FEA) program. Typically

  19. Pollution error in the h-version of the finite-element method and the local quality of a-posteriori error estimators

    E-Print Network [OSTI]

    Mathur, Anuj

    1994-01-01T23:59:59.000Z

    in the L-shaped domain shown in Figs. la-ld. We applied boundary condi- 1 tions consistent with the singular solution u(r, 8) = rr sin(-) in the infinite wedge (see also Chapter II) and we assumed that we are interested only in the solution in the hol...-element solution in the subdomain. In Fig. l. la we show an adaptive finite-element mesh which was gener- ated using remeshing and linear elements, over the L-shaped domain fl. For this mesh an accuracy of 1. 09% for the relative error in the energy-norm over...

  20. A finite element method with mesh adaptivity for computing vortex states in fast-rotating Bose-Einstein condensates

    E-Print Network [OSTI]

    Boyer, Edmond

    -Einstein condensates Ionut Danaila,a,b , Fr´ed´eric Hechta,b aUPMC Univ Paris 06, UMR 7598, Laboratoire Jacques Abstract Numerical computations of stationary states of fast-rotating Bose-Einstein condensates require-Einstein condensate, vortex, Sobolev gradient, descent method. 1. Introduction Recent research efforts in the field

  1. A finite element method with mesh adaptivity for computing vortex states in fastrotating BoseEinstein condensates

    E-Print Network [OSTI]

    Recanati, Catherine

    ­Einstein condensates Ionut Danaila #,a,b , Frâ??edâ??eric Hecht a,b a UPMC Univ Paris 06, UMR 7598, Laboratoire Jacques Abstract Numerical computations of stationary states of fast­rotating Bose­Einstein condensates require­Einstein condensate, vortex, Sobolev gradient, descent method. 1. Introduction Recent research e#orts in the field

  2. Finite Elements in Analysis and Design 40 (2004) 16191640 www.elsevier.com/locate/nel

    E-Print Network [OSTI]

    Ghosh, Somnath

    2004-01-01T23:59:59.000Z

    Finite Elements in Analysis and Design 40 (2004) 1619­1640 www.elsevier.com/locate/ÿnel Two scale. Raghavan et al. / Finite Elements in Analysis and Design 40 (2004) 1619­1640 performed by the method displace- ment-based ÿnite element model with a microstructural Voronoi cell FEM (VCFEM) for multi

  3. Least-squares variational principles and the finite element method: theory, formulations, and models for solid and fluid mechanics

    E-Print Network [OSTI]

    Pontaza, Juan Pablo

    2004-09-30T23:59:59.000Z

    of the weak form Galerkin procedure to problems whose solutions can be characterized as global minimizers results in global minimiza- tion of quadratic functionals, such is the case for linear elasticity problems; see [93]. In this case, given a conforming... and the success of the method for linear elasticity problems led to early attempts to extend weak form Galerkin formulations beyond problems whose solutions could be characterized as unconstrained global minimizers. For example, problems of the saddle-point type...

  4. A Quadrature Finite Element Galerkin Scheme for a Biharmonic Problem on a Rectangular Polygon

    E-Print Network [OSTI]

    Aitbayev, Rakhim

    to the OSC method (see [3]), which uses Gaussian quadrature nodes as collocation points and a finite elementA Quadrature Finite Element Galerkin Scheme for a Biharmonic Problem on a Rectangular Polygon with an underdetermined orthogonal spline collocation scheme. © 2007 Wiley Periodicals, Inc. Numer Methods Partial

  5. Impeller deflection and modal finite element analysis.

    SciTech Connect (OSTI)

    Spencer, Nathan A.

    2013-10-01T23:59:59.000Z

    Deflections of an impeller due to centripetal forces are calculated using finite element analysis. The lateral, or out of plane, deflections are an important design consideration for this particular impeller because it incorporates an air bearing with critical gap tolerances. The target gap distance is approximately 10 microns at a rotational velocity of 2500 rpm. The centripetal forces acting on the impeller cause it deflect in a concave fashion, decreasing the initial gap distance as a function of radial position. This deflection is characterized for a previous and updated impeller design for comparative purposes. The impact of design options such as material selection, geometry dimensions, and operating rotational velocity are also explored, followed by a sensitivity study with these parameters bounded by specific design values. A modal analysis is also performed to calculate the impeller's natural frequencies which are desired to be avoided during operation. The finite element modeling techniques continue to be exercised by the impeller design team to address specific questions and evaluate conceptual designs, some of which are included in the Appendix.

  6. A NONCONFORMING FINITE ELEMENT METHOD

    E-Print Network [OSTI]

    Karakashian, Ohannes

    . Karakashiany and Wadi N. Jureidiniz Abstract. Approximations to solutions of the inhomogeneous boundary value

  7. Tetrahedral hp finite elements: Algorithms and flow simulations

    SciTech Connect (OSTI)

    Sherwin, S.J.; Karniadakis, G.E. [Brown Univ., Providence RI (United States)] [Brown Univ., Providence RI (United States)

    1996-03-01T23:59:59.000Z

    We present a new discretisation for the incompressible Navier-Stokes equations that extends spectral methods to three-dimensional complex domains consisting of tetrahedral subdomains. The algorithm is based on standard concepts of hp finite elements as well as tensorial spectral elements. This new formulation employs a hierarchical/modal basis constructed from a new apex co-ordinate system which retains a generalised tensor product. These properties enable the development of computationally efficienct algorithms for use on standard finite volume unstructed meshes. A detailed analysis is presented that documents the stability and exponential convergence of the method and several flow cases are simulated and compared with analytical and experimental results. 34 refs., 28 figs., 1 tab.

  8. Summary of the Models and Methods for the FEHM Application-A Finite-Element Heat- and Mass-Transfer Code

    SciTech Connect (OSTI)

    George A. Zyvoloski; Bruce A. Robinson; Zora V. Dash; Lynn L. Trease

    1997-07-01T23:59:59.000Z

    The mathematical models and numerical methods employed by the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multi-component flow in porous media, are described. The use of this code is applicable to natural-state studies of geothermal systems and groundwater flow. A primary use of the FEHM application will be to assist in the understanding of flow fields and mass transport in the saturated and unsaturated zones below the proposed Yucca Mountain nuclear waste repository in Nevada. The component models of FEHM are discussed. The first major component, Flow- and Energy-Transport Equations, deals with heat conduction; heat and mass transfer with pressure- and temperature-dependent properties, relative permeabilities and capillary pressures; isothermal air-water transport; and heat and mass transfer with noncondensible gas. The second component, Dual-Porosity and Double-Porosity/Double-Permeability Formulation, is designed for problems dominated by fracture flow. Another component, The Solute-Transport Models, includes both a reactive-transport model that simulates transport of multiple solutes with chemical reaction and a particle-tracking model. Finally, the component, Constitutive Relationships, deals with pressure- and temperature-dependent fluid/air/gas properties, relative permeabilities and capillary pressures, stress dependencies, and reactive and sorbing solutes. Each of these components is discussed in detail, including purpose, assumptions and limitations, derivation, applications, numerical method type, derivation of numerical model, location in the FEHM code flow, numerical stability and accuracy, and alternative approaches to modeling the component.

  9. Geodesic Finite Elements of Higher Order Oliver Sander

    E-Print Network [OSTI]

    conditions. As a side product we can define geodesic finite elements for non-simplex reference elements of a simplex by the weighted average (v1, . . . , dd+1; w) := arg min qM d+1 i=1 wi dist(vi, q)2 , (1) where w are barycentric coordinates on the simplex. Based on this formula a finite element theory could be constructed

  10. Implementation of B-splines in a Conventional Finite Element Framework

    E-Print Network [OSTI]

    Owens, Brian C.

    2010-01-16T23:59:59.000Z

    The use of B-spline interpolation functions in the finite element method (FEM) is not a new subject. B-splines have been utilized in finite elements for many reasons. One reason is the higher continuity of derivatives and smoothness of B...

  11. A ThreeDimensional Finite Element Simulation for Transport of Nuclear Waste Contamination in Porous Media

    E-Print Network [OSTI]

    Ewing, Richard E.

    A Three­Dimensional Finite Element Simulation for Transport of Nuclear Waste Contamination of South Carolina, Columbia, South Carolina ABSTRACT: Model equations for transport of nuclear­waste based up on the inherent physics. A three­dimensional finite element method for nuclear waste

  12. FINITE-ELEMENT FORMULATIONS FOR PROBLEMS OF LARGE ELASTIC-PLASTIC DEFORMATION

    E-Print Network [OSTI]

    FINITE-ELEMENT FORMULATIONS FOR PROBLEMS OF LARGE ELASTIC-PLASTIC DEFORMATION R. M. MCM~EKING and J elastic-plastic flow.The method is based on Hill's variational principle for incremental deformations in a manner which allows any conventions finite element program, for "small strain" elastic-plastic analysis

  13. CONVERGENCE ANALYSIS OF A QUADRATURE FINITE ELEMENT GALERKIN SCHEME FOR A BIHARMONIC PROBLEM

    E-Print Network [OSTI]

    Aitbayev, Rakhim

    of the quadrature problem with an orthogonal spline collocation problem. Key words: biharmonic problem, finite elements, Galerkin method, Gaussian quadrature, orthogonal spline collocation AMS subject classificationCONVERGENCE ANALYSIS OF A QUADRATURE FINITE ELEMENT GALERKIN SCHEME FOR A BIHARMONIC PROBLEM RAKHIM

  14. Bending buckling of single-walled carbon nanotubes by atomic-scale finite element

    E-Print Network [OSTI]

    Jiang, Hanqing

    - vered SWNT loaded by a follower lateral force at the free end and found that it exhibited local bucklingBending buckling of single-walled carbon nanotubes by atomic-scale finite element X. Guo a , A Abstract This paper employs the atomic-scale finite element method to study bending buckling of single

  15. A finite element multigrid preconditioner for Chebyshev–collocation ...

    E-Print Network [OSTI]

    2000-03-20T23:59:59.000Z

    Keywords: Multigrid; Preconditioning; Finite element; Chebyshev–collocation. 1. .... Mi be the restriction of A on Mi; Qi,Pi :M ? Mi be the orthogonal projections.

  16. Algebraic multigrid for stabilized finite element discretizations of the Navier Stokes equation

    E-Print Network [OSTI]

    Okusanya, Tolulope Olawale, 1972 -

    2002-01-01T23:59:59.000Z

    A multilevel method for the solution of systems of equations generated by stabilized Finite Element discretizations of the Euler and Navier Stokes equations on generalized unstructured grids is described. The method is ...

  17. An Inverse Finite Element Analysis and A Parametric Study of Small Punch Tests 

    E-Print Network [OSTI]

    Xu, Zhenzhen

    2012-02-14T23:59:59.000Z

    the effects of SPT parameters including friction coefficient, punch head diameter, sample thickness, specimen scale and boundary conditions. The proposed inverse finite element (FE) method improves the accuracy of existing inverse FE methods...

  18. An Inverse Finite Element Analysis and A Parametric Study of Small Punch Tests

    E-Print Network [OSTI]

    Xu, Zhenzhen

    2012-02-14T23:59:59.000Z

    the effects of SPT parameters including friction coefficient, punch head diameter, sample thickness, specimen scale and boundary conditions. The proposed inverse finite element (FE) method improves the accuracy of existing inverse FE methods...

  19. HYDRA, A finite element computational fluid dynamics code: User manual

    SciTech Connect (OSTI)

    Christon, M.A.

    1995-06-01T23:59:59.000Z

    HYDRA is a finite element code which has been developed specifically to attack the class of transient, incompressible, viscous, computational fluid dynamics problems which are predominant in the world which surrounds us. The goal for HYDRA has been to achieve high performance across a spectrum of supercomputer architectures without sacrificing any of the aspects of the finite element method which make it so flexible and permit application to a broad class of problems. As supercomputer algorithms evolve, the continuing development of HYDRA will strive to achieve optimal mappings of the most advanced flow solution algorithms onto supercomputer architectures. HYDRA has drawn upon the many years of finite element expertise constituted by DYNA3D and NIKE3D Certain key architectural ideas from both DYNA3D and NIKE3D have been adopted and further improved to fit the advanced dynamic memory management and data structures implemented in HYDRA. The philosophy for HYDRA is to focus on mapping flow algorithms to computer architectures to try and achieve a high level of performance, rather than just performing a port.

  20. Energy conserving Anisotropic Anhysteretic Magnetic Modelling for Finite Element Analysis

    E-Print Network [OSTI]

    Jens Krause

    2012-12-20T23:59:59.000Z

    To model ferromagnetic material in finite element analysis a correct description of the constitutive relationship (BH-law) must be found from measured data. This article proposes to use the energy density function as a centrepiece. Using this function, which turns out to be a convex function of the flux density, guarantees energy conservative modelling. The magnetic field strength can be seen as a derivative with respect to the flux density. Especially for anisotropic materials (from lamination and/or grain orientation) this method has advantages. Strictly speaking this method is only valid for anhysteretic and thermodynamically stable material.

  1. Finite element analysis of heat transport in a hydrothermal zone

    SciTech Connect (OSTI)

    Bixler, N.E.; Carrigan, C.R.

    1987-01-01T23:59:59.000Z

    Two-phase heat transport in the vicinity of a heated, subsurface zone is important for evaluation of nuclear waste repository design and estimation of geothermal energy recovery, as well as prediction of magma solidification rates. Finite element analyses of steady, two-phase, heat and mass transport have been performed to determine the relative importance of conduction and convection in a permeable medium adjacent to a hot, impermeable, vertical surface. The model includes the effects of liquid flow due to capillarity and buoyancy and vapor flow due to pressure gradients. Change of phase, with its associated latent heat effects, is also modeled. The mechanism of capillarity allows for the presence of two-phase zones, where both liquid and vapor can coexist, which has not been considered in previous investigations. The numerical method employs the standard Galerkin/finite element method, using eight-node, subparametric or isoparametric quadrilateral elements. In order to handle the extreme nonlinearities inherent in two-phase, nonisothermal, porous-flow problems, steady-state results are computed by integrating transients out to a long time (a method that is highly robust).

  2. FETI-R (Part II): A regularized variational approach for Finite Element Tearing and

    E-Print Network [OSTI]

    FETI-R (Part II): A regularized variational approach for Finite Element Tearing and Interconnecting Laboratories y Albuquerque, New Mexico 87185 Abstract The coarse grid problem of the original FETI method

  3. Finite element decomposition and grid generation for brain modeling and visualization

    E-Print Network [OSTI]

    Batte, David Allan

    1997-01-01T23:59:59.000Z

    Numerical grid generation is used to provide a framework for brain and neuron visualization. Smoothing spline surfaces are fit to contour data to generate 3D solid model reconstruction of brain tissues. Finite element methods are then used...

  4. 3-D Finite Element Electromagnetic and Stress Analyses of the JET LB-SRP Divertor Element (Tungsten Lamella Design)

    E-Print Network [OSTI]

    3-D Finite Element Electromagnetic and Stress Analyses of the JET LB-SRP Divertor Element (Tungsten Lamella Design)

  5. JAC3D -- A three-dimensional finite element computer program for the nonlinear quasi-static response of solids with the conjugate gradient method; Yucca Mountain Site Characterization Project

    SciTech Connect (OSTI)

    Biffle, J.H.

    1993-02-01T23:59:59.000Z

    JAC3D is a three-dimensional finite element program designed to solve quasi-static nonlinear mechanics problems. A set of continuum equations describes the nonlinear mechanics involving large rotation and strain. A nonlinear conjugate gradient method is used to solve the equation. The method is implemented in a three-dimensional setting with various methods for accelerating convergence. Sliding interface logic is also implemented. An eight-node Lagrangian uniform strain element is used with hourglass stiffness to control the zero-energy modes. This report documents the elastic and isothermal elastic-plastic material model. Other material models, documented elsewhere, are also available. The program is vectorized for efficient performance on Cray computers. Sample problems described are the bending of a thin beam, the rotation of a unit cube, and the pressurization and thermal loading of a hollow sphere.

  6. JAC2D: A two-dimensional finite element computer program for the nonlinear quasi-static response of solids with the conjugate gradient method; Yucca Mountain Site Characterization Project

    SciTech Connect (OSTI)

    Biffle, J.H.; Blanford, M.L.

    1994-05-01T23:59:59.000Z

    JAC2D is a two-dimensional finite element program designed to solve quasi-static nonlinear mechanics problems. A set of continuum equations describes the nonlinear mechanics involving large rotation and strain. A nonlinear conjugate gradient method is used to solve the equations. The method is implemented in a two-dimensional setting with various methods for accelerating convergence. Sliding interface logic is also implemented. A four-node Lagrangian uniform strain element is used with hourglass stiffness to control the zero-energy modes. This report documents the elastic and isothermal elastic/plastic material model. Other material models, documented elsewhere, are also available. The program is vectorized for efficient performance on Cray computers. Sample problems described are the bending of a thin beam, the rotation of a unit cube, and the pressurization and thermal loading of a hollow sphere.

  7. NORTHWESTERN UNIVERSITY Finite Element Analysis of TDR Cable-Grout-Soil Mass

    E-Print Network [OSTI]

    NORTHWESTERN UNIVERSITY Finite Element Analysis of TDR Cable-Grout-Soil Mass Interaction During............................................................................................ 13 Finite Element Analysis of Cable Shear Test............................................................................................ 22 Finite Element Analysis of Grout-Cable Composite Shear Test .................................. 25

  8. 1 INTRODUCTION Nowadays, geomechanical finite-element codes

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    1 INTRODUCTION Nowadays, geomechanical finite-element codes running on modern computers allow for years. Geotechnical characterization, geophysical surveys and geomechanical modelling have been in geomechanical modelling. In particular, seismic and borehole data, obtained during previous campaigns, are used

  9. A frequency domain finite element model for tidal circulation

    E-Print Network [OSTI]

    Westerink, Joannes J.

    1985-01-01T23:59:59.000Z

    A highly efficient finite element model has been developed for the numerical prediction of depth average circulation within small scale embayments which are often characterized by irregular boundaries and bottom topography.

  10. Transactions of NAMRI/SME FINITE ELEMENT SIMULATION OF

    E-Print Network [OSTI]

    Ozel, Tugrul

    Transactions of NAMRI/SME FINITE ELEMENT SIMULATION OF Ti-6Al-4V ALLOY Turul Manufacturing #12;Transactions of NAMRI/SME 50 Volume 38, 2010 excessive heat generation. Usually above 120 m

  11. Progress in Developing Finite Element Models Replicating Flexural Graphite Testing

    SciTech Connect (OSTI)

    Robert Bratton

    2010-06-01T23:59:59.000Z

    This report documents the status of flexural strength evaluations from current ASTM procedures and of developing finite element models predicting the probability of failure. This work is covered under QLD REC-00030. Flexural testing procedures of the American Society for Testing and Materials (ASTM) assume a linear elastic material that has the same moduli for tension and compression. Contrary to this assumption, graphite is known to have different moduli for tension and compression. A finite element model was developed and demonstrated that accounts for the difference in moduli tension and compression. Brittle materials such as graphite exhibit significant scatter in tensile strength, so probabilistic design approaches must be used when designing components fabricated from brittle materials. ASTM procedures predicting probability of failure in ceramics were compared to methods from the current version of the ASME graphite core components rules predicting probability of failure. Using the ASTM procedures yields failure curves at lower applied forces than the ASME rules. A journal paper was published in the Journal of Nuclear Engineering and Design exploring the statistical models of fracture in graphite.

  12. Omega3P: A Parallel Finite-Element Eigenmode Analysis Code for Accelerator Cavities

    SciTech Connect (OSTI)

    Lee, Lie-Quan; Li, Zenghai; Ng, Cho; Ko, Kwok; /SLAC

    2009-03-04T23:59:59.000Z

    Omega3P is a parallel eigenmode calculation code for accelerator cavities in frequency domain analysis using finite-element methods. In this report, we will present detailed finite-element formulations and resulting eigenvalue problems for lossless cavities, cavities with lossy materials, cavities with imperfectly conducting surfaces, and cavities with waveguide coupling. We will discuss the parallel algorithms for solving those eigenvalue problems and demonstrate modeling of accelerator cavities through different examples.

  13. Surface photovoltage measurements and finite element modeling of SAW devices.

    SciTech Connect (OSTI)

    Donnelly, Christine

    2012-03-01T23:59:59.000Z

    Over the course of a Summer 2011 internship with the MEMS department of Sandia National Laboratories, work was completed on two major projects. The first and main project of the summer involved taking surface photovoltage measurements for silicon samples, and using these measurements to determine surface recombination velocities and minority carrier diffusion lengths of the materials. The SPV method was used to fill gaps in the knowledge of material parameters that had not been determined successfully by other characterization methods. The second project involved creating a 2D finite element model of a surface acoustic wave device. A basic form of the model with the expected impedance response curve was completed, and the model is ready to be further developed for analysis of MEMS photonic resonator devices.

  14. Introduction to Finite Element Modeling Engineering analysis of mechanical systems have been addressed by deriving differential

    E-Print Network [OSTI]

    Lin, Liwei

    Introduction to Finite Element Modeling Engineering analysis of mechanical systems have been of the FEM (don't misuse the FEM - a numerical tool) Finite Element Analysis A typical finite element quantities (e.g., strains and stresses) at selected elements Basic Theory The way finite element analysis

  15. A Novel Membrane Finite Element with Drilling Rotations

    E-Print Network [OSTI]

    Kouhia, Reijo

    A Novel Membrane Finite Element with Drilling Rotations Reijo Kouhia 1 Abstract. A new low order interpolation is used for the drill rotation #12;eld. Both triangular and quadrilateral elements are considered of freedom. 1 INTRODUCTION In-plane rotational degrees of freedom, \\drilling de- grees of freedom

  16. Phase reduction models for improving the accuracy of the finite element solution of time-harmonic scattering problems I: General

    E-Print Network [OSTI]

    Antoine, Xavier - Institut de Mathématiques �lie Cartan, Université Henri Poincaré

    Accepted 12 January 2009 Available online 21 January 2009 Keywords: Helmholtz equation Acoustic scattering Short-wave problem Finite element method Pollution Accuracy On-Surface Radiation Condition method a b- lems in view of a numerical finite element solution. It is well-known that pollution error causes

  17. CharacteristicGalerkin and Mixed Finite Element Approximation of Contamination by Compressible

    E-Print Network [OSTI]

    Chou, So-Hsiang

    for the brine, radionuclides, and heat and by a mixed finite element method for the pressure and velocity disposal buried in a salt dome. The salt dissolves to generate a brine, radioactive elements decay), the viscosity of the fluid, is dependent on â?? c, the concentration of the brine in the fluid. For simplicity we

  18. A Mixed Finite-Element Discretization of the Energy-Transport Model for Semiconductors

    E-Print Network [OSTI]

    Pietra, Paola

    A Mixed Finite-Element Discretization of the Energy-Transport Model for Semiconductors Stefan Holst #12;tting mixed #12;nite-element method is used to discretize the stationary energy. Energy-transport models describe the ow of electrons through a semi- conductor device, in uenced by di

  19. Preprint version, February 2006 to appear in Acta Numerica (2006) Finite element exterior calculus,

    E-Print Network [OSTI]

    2006-01-01T23:59:59.000Z

    problem being solved. In the finite element exterior calculus, many finite element spaces are re- vealed-dimensional subspace of the given #12;Finite element exterior calculus 3 function space where the solution isPreprint version, February 2006 to appear in Acta Numerica (2006) Finite element exterior calculus

  20. An Online Data Access System for a Finite Element Program , David Liu2

    E-Print Network [OSTI]

    Stanford University

    for a finite element analysis (FEA) program. The system incorporates a commercial off-the-shelf (COTS) database also be observed in the field of finite element analysis (FEA). Modern finite element programs1 An Online Data Access System for a Finite Element Program Jun Peng1 , David Liu2 , and Kincho H

  1. Parallel finite element algorithm for three dimensional ...

    E-Print Network [OSTI]

    methods for electromagnetic induction. INTRODUCTION. N umerical modelling of three dimensional conductivity structures in the earth has experienced a rapid ...

  2. Nonlinear Legendre Spectral Finite Elements for Wind Turbine Blade Dynamics: Preprint

    SciTech Connect (OSTI)

    Wang, Q.; Sprague, M. A.; Jonkman, J.; Johnson, N.

    2014-01-01T23:59:59.000Z

    This paper presents a numerical implementation and examination of new wind turbine blade finite element model based on Geometrically Exact Beam Theory (GEBT) and a high-order spectral finite element method. The displacement-based GEBT is presented, which includes the coupling effects that exist in composite structures and geometric nonlinearity. Legendre spectral finite elements (LSFEs) are high-order finite elements with nodes located at the Gauss-Legendre-Lobatto points. LSFEs can be an order of magnitude more efficient that low-order finite elements for a given accuracy level. Interpolation of the three-dimensional rotation, a major technical barrier in large-deformation simulation, is discussed in the context of LSFEs. It is shown, by numerical example, that the high-order LSFEs, where weak forms are evaluated with nodal quadrature, do not suffer from a drawback that exists in low-order finite elements where the tangent-stiffness matrix is calculated at the Gauss points. Finally, the new LSFE code is implemented in the new FAST Modularization Framework for dynamic simulation of highly flexible composite-material wind turbine blades. The framework allows for fully interactive simulations of turbine blades in operating conditions. Numerical examples showing validation and LSFE performance will be provided in the final paper.

  3. Study of the mixed finite volume method for Stokes and Navier-Stokes equations

    E-Print Network [OSTI]

    Droniou, Jérôme

    (without "orthogonality" assumptions as for classical finite volume methods) and in any space dimension. We methods: finite element methods, mostly used by the mathematician community (see for example [10, 11, 12, and for example to [2, 13, 8, 9] for collocated finite volume schemes). One reason for this difference of practice

  4. Materials Science and Engineering B 117 (2005) 5361 Finite element analysis-based design of a fluid-flow control nano-valve

    E-Print Network [OSTI]

    Grujicic, Mica

    2005-01-01T23:59:59.000Z

    Materials Science and Engineering B 117 (2005) 53­61 Finite element analysis-based design A finite element method-based procedure is developed for the design of molecularly functionalized nano of a fluid-flow control nano-valve. The results obtained suggest that the finite element-based procedure

  5. APPLICATIONS OF ALGEBRAIC MULTIGRID TO LARGE-SCALE FINITE ELEMENT ANALYSIS OF WHOLE BONE MICRO-MECHANICS ON THE IBM SP

    E-Print Network [OSTI]

    California at Berkeley, University of

    APPLICATIONS OF ALGEBRAIC MULTIGRID TO LARGE-SCALE FINITE ELEMENT ANALYSIS OF WHOLE BONE MICRO,5 Abstract. Accurate micro-finite element analyses of whole bones require the solution of large sets architectures. Key words. multigrid, trabecular bone, human vertebral body, finite element method, massively

  6. Abundant p-singular elements in finite classical groups

    E-Print Network [OSTI]

    Niemeyer, Alice C; Praeger, Cheryl E

    2012-01-01T23:59:59.000Z

    In 1995, Isaacs, Kantor and Spaltenstein proved that for a finite simple d-dimensional classical group G, and for any prime divisor p of |G| distinct from the characteristic, the proportion of p-singular elements (elements with order divisible by p) is at least c/d for some constant c. We define a new subfamily of p-singular elements, called p-abundant elements, which leave invariant certain "large" subspaces of the natural G-module. We find explicit upper and lower bounds for the proportion of p-abundant elements in G, and prove that it approaches a (positive) limiting value as d \\rightarrow \\infty. It turns out that the limiting proportion of p-abundant elements is at least a constant multiple of the lower bound for the proportion of all p-singular elements derived by Issacs, Kantor and Spaltenstein.

  7. Inversion of Robin coefficient by a spectral stochastic finite element approach

    SciTech Connect (OSTI)

    Jin Bangti [Department of Mathematics, Chinese University of Hong Kong, Shatin, N.T., Hong Kong (China)], E-mail: kimbtsing@yahoo.com.cn; Zou Jun [Department of Mathematics, Chinese University of Hong Kong, Shatin, N.T., Hong Kong (China)], E-mail: zou@math.cuhk.edu.hk

    2008-03-01T23:59:59.000Z

    This paper investigates a variational approach to the nonlinear stochastic inverse problem of probabilistically calibrating the Robin coefficient from boundary measurements for the steady-state heat conduction. The problem is formulated into an optimization problem, and mathematical properties relevant to its numerical computations are investigated. The spectral stochastic finite element method using polynomial chaos is utilized for the discretization of the optimization problem, and its convergence is analyzed. The nonlinear conjugate gradient method is derived for the optimization system. Numerical results for several two-dimensional problems are presented to illustrate the accuracy and efficiency of the stochastic finite element method.

  8. Diffusive mesh relaxation in ALE finite element numerical simulations

    SciTech Connect (OSTI)

    Dube, E.I.

    1996-06-01T23:59:59.000Z

    The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.

  9. Experimentally validated finite element model of electrocaloric multilayer ceramic structures

    SciTech Connect (OSTI)

    Smith, N. A. S., E-mail: nadia.smith@npl.co.uk, E-mail: maciej.rokosz@npl.co.uk, E-mail: tatiana.correia@npl.co.uk; Correia, T. M., E-mail: nadia.smith@npl.co.uk, E-mail: maciej.rokosz@npl.co.uk, E-mail: tatiana.correia@npl.co.uk [National Physical Laboratory, Hampton Road, TW11 0LW Middlesex (United Kingdom); Rokosz, M. K., E-mail: nadia.smith@npl.co.uk, E-mail: maciej.rokosz@npl.co.uk, E-mail: tatiana.correia@npl.co.uk [National Physical Laboratory, Hampton Road, TW11 0LW Middlesex (United Kingdom); Department of Materials, Imperial College London, London SW7 2AZ (United Kingdom)

    2014-07-28T23:59:59.000Z

    A novel finite element model to simulate the electrocaloric response of a multilayer ceramic capacitor (MLCC) under real environment and operational conditions has been developed. The two-dimensional transient conductive heat transfer model presented includes the electrocaloric effect as a source term, as well as accounting for radiative and convective effects. The model has been validated with experimental data obtained from the direct imaging of MLCC transient temperature variation under application of an electric field. The good agreement between simulated and experimental data, suggests that the novel experimental direct measurement methodology and the finite element model could be used to support the design of optimised electrocaloric units and operating conditions.

  10. Inverse Finite Element Modelling and Identification of Constitutive Parameters of UHS Steel Based on Gleeble Tensile Tests at High Temperature

    E-Print Network [OSTI]

    Boyer, Edmond

    identification method - associating a direct finite element model of Gleeble tests and an optimization module1 Inverse Finite Element Modelling and Identification of Constitutive Parameters of UHS Steel Based-2 ) hc heat transfer coefficient at interface between specimen and grips (W m-2 K) hth_eff effective

  11. Implementation of surface tension with wall adhesion effects in a three-dimensional finite element model for fluid flow

    E-Print Network [OSTI]

    Boyer, Edmond

    - 1 - Implementation of surface tension with wall adhesion effects in a three-dimensional finite element modelling of surface tension. The external stress vectors associated with surface tension a drop of liquid on a plane is treated. Keywords : surface tension, finite element method, average

  12. J. Sulem & M. Cerrolaza (2002): Finite element analysis of the indentation test on rocks with microstructure 1 FINITE ELEMENT ANALYSIS OF THE INDENTATION TEST

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    2002-01-01T23:59:59.000Z

    J. Sulem & M. Cerrolaza (2002): Finite element analysis of the indentation test on rocks with microstructure 1 FINITE ELEMENT ANALYSIS OF THE INDENTATION TEST ON ROCKS WITH MICROSTRUCTURE Published in-117" #12;J. Sulem & M. Cerrolaza (2002): Finite element analysis of the indentation test on rocks

  13. An AMR Capable Finite Element Diffusion Solver

    Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)

    Rd. Berkeley, CA 94720, USA E-mail: fisher47@llnl.gov Abstract. We present a novel method for the solution of the diffusion equation on a composite AMR mesh. This approach is...

  14. Method for detecting an element

    DOE Patents [OSTI]

    Blackwood, Larry G.; Reber, Edward L.; Rohde, Kenneth W.

    2007-02-06T23:59:59.000Z

    A method for detecting an element is disclosed and which includes the steps of providing a gamma-ray spectrum which depicts, at least in part, a test region having boundaries, and which has a small amount of the element to be detected; providing a calculation which detects the small amount of the element to be detected; and providing a moving window and performing the calculation within the moving window, and over a range of possible window boundaries within the test region to determine the location of the optimal test region within the gamma-ray spectrum.

  15. Experimental validation of finite element codes for welding deformations

    E-Print Network [OSTI]

    Boyer, Edmond

    Experimental validation of finite element codes for welding deformations H. M. Aarbogha,b, , M Institute for Energy Technology, N-2027 Kjeller, Norway. Abstract A single pass Metal Inert Gas welding which numerical codes quantifying welding stresses can be validated. It includes a mov- ing heat source

  16. CONTINUOUSTIME FINITE ELEMENT ANALYSIS OF MULTIPHASE FLOW IN GROUNDWATER HYDROLOGY

    E-Print Network [OSTI]

    CONTINUOUS­TIME FINITE ELEMENT ANALYSIS OF MULTIPHASE FLOW IN GROUNDWATER HYDROLOGY Zhangxin Chen­water system in groundwater hydrology is given. The system is written in a fractional flow formulation, i for an air­water system in groundwater hydrology, ff = a; w [1], [11], [26]: @(OEae ff s ff ) @t +r \\Delta

  17. Dennis, Eberhart, Dulikravich & Radons FINITE ELEMENT SIMULATION OF COOLING

    E-Print Network [OSTI]

    Dennis, Brian

    arterial supply and effectively blocks the external cooling wave from advancing to the core of the brain1 Dennis, Eberhart, Dulikravich & Radons FINITE ELEMENT SIMULATION OF COOLING OF REALISTIC 3-D Rapid cooling of the brain in the first minutes following the onset of cerebral ischemia

  18. FEWA: a Finite Element model of Water flow through Aquifers

    SciTech Connect (OSTI)

    Yeh, G.T.; Huff, D.D.

    1983-11-01T23:59:59.000Z

    This report documents the implementation and demonstration of a Finite Element model of Water flow through Aquifers (FEWA). The particular features of FEWA are its versatility and flexibility to deal with as many real-world problems as possible. Point as well as distributed sources/sinks are included to represent recharges/pumpings and rainfall infiltrations. All sources/sinks can be transient or steady state. Prescribed hydraulic head on the Dirichlet boundaries and fluxes on Neumann or Cauchy boundaries can be time-dependent or constant. Source/sink strength over each element and node, hydraulic head at each Dirichlet boundary node, and flux at each boundary segment can vary independently of each other. Either completely confined or completely unconfined aquifers, or partially confined and partially unconfined aquifers can be dealt with effectively. Discretization of a compound region with very irregular curved boundaries is made easy by including both quadrilateral and triangular elements in the formulation. Large-field problems can be solved efficiently by including a pointwise iterative solution strategy as an optional alternative to the direct elimination solution method for the matrix equation approximating the partial differential equation of groundwater flow. FEWA also includes transient flow through confining leaky aquifers lying above and/or below the aquifer of interest. The model is verified against three simple cases to which analytical solutions are available. It is then demonstrated by two examples of how the model can be applied to heterogeneous and anisotropic aquifers with transient boundary conditions, time-dependent sources/sinks, and confining aquitards for a confined aquifer of variable thickness and for a free surface problem in an unconfined aquifer, respectively. 20 references, 25 figures, 8 tables.

  19. FINITE ELEMENT ANALYSIS OF THERMAL TENSIONING TECHNIQUES MITIGATING WELD BUCKLING DISTORTION

    E-Print Network [OSTI]

    Michaleris, Panagiotis

    FINITE ELEMENT ANALYSIS OF THERMAL TENSIONING TECHNIQUES MITIGATING WELD BUCKLING DISTORTION. This paper presents a finite element analysis model of the thermal tensioning technique. A series of finite by the finite element simulations, the residual stresses of large size and high heat input welds are reduced

  20. SUBDIVISION SURFACES: A NEW PARADIGM FOR THIN-SHELL FINITE-ELEMENT ANALYSIS

    E-Print Network [OSTI]

    Desbrun, Mathieu

    SUBDIVISION SURFACES: A NEW PARADIGM FOR THIN-SHELL FINITE-ELEMENT ANALYSIS Fehmi Cirak1 , Michael for thin-shell finite-element analysis based on the use of subdivision surfaces for: i) describing by subdivision are H2 and, consequently, have a finite Kirchhoff-Love energy. The resulting finite elements

  1. A Mixed Finite Element Approximation for Compressible Flow of Contamination from Nuclear

    E-Print Network [OSTI]

    Chou, So-Hsiang

    of this system using a finite element method for the brine, radionuclides, and heat combined with a mixed finite; and ¯(â??c), the viscosity of the fluid is depen­ dent on â?? c, the concentration of the brine in the fluid + (ff L \\Gamma ff T )u i u j =juj). Brine: @ @t (OEaeâ??c) = \\Gammar \\Delta (aeâ??cu) +r \\Delta (aeE c (u

  2. Finite element analysis of axisymmetric oscillations of sessile liquid drops

    SciTech Connect (OSTI)

    Bixler, N.E.; Benner, R.E.

    1985-01-01T23:59:59.000Z

    Inviscid oscillations of sessile liquid drops are simulated by the Galerkin finite element method in conjunction with the time integrator proposed by Gresho et al. Simulations are of drops in spherical containers which are subjected to imposed oscillations of specified frequency and amplitude. Five equations govern drop response: (1) Laplace's equation for velocity potential within the drop; (2) a kinematic condition at the free surface; (3) a Bernoulli equation augmented to include gravity and capillary pressure at the free surface; (4) a kinematic condition at the solid surface; and (5) either a condition for fixed contact line or fixed contact angle. Each of these equations is modified to account for an accelerating frame of reference which moves the container. Normalized drop volume, contact angle, and gravitational Bond number are dimensionless parameters which control drop response to an imposed oscillation. Given a set of fluid properties, such as those for mercury, gravitational Bond number is uniquely defined by the container radius. Resonant frequencies and mode interaction are detected by Fourier analysis of a transient signal, such as free surface position at the pole of a spherical coordinate system. Results, especially resonant frequencies, are found to depend strongly on contact line condition. Calculation of resonant frequencies by eigenanalysis with Stewart's method is also discussed. 11 refs., 8 figs.

  3. On Smoothing Surfaces in Voxel Based Finite Element Analysis of Trabecular Bone

    E-Print Network [OSTI]

    Frey, Pascal

    On Smoothing Surfaces in Voxel Based Finite Element Analysis of Trabecular Bone Peter Arbenz on complicated domains composed of often hundreds of millions of voxel elements. The finite element analysis finite element (FE) analysis. The approach based on the FE analysis leads to linear systems of equations

  4. Time-independent hybrid enrichment for finite element solution of transient conduction–radiation in diffusive grey media

    SciTech Connect (OSTI)

    Mohamed, M. Shadi, E-mail: m.s.mohamed@durham.ac.uk [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom); Seaid, Mohammed; Trevelyan, Jon [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom)] [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom); Laghrouche, Omar [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)] [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)

    2013-10-15T23:59:59.000Z

    We investigate the effectiveness of the partition-of-unity finite element method for transient conduction–radiation problems in diffusive grey media. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary diffusion approximation to the radiation in grey media. The coupled equations are integrated in time using a semi-implicit method in the finite element framework. We show that for the considered problems, a combination of hyperbolic and exponential enrichment functions based on an approximation of the boundary layer leads to improved accuracy compared to the conventional finite element method. It is illustrated that this approach can be more efficient than using h adaptivity to increase the accuracy of the finite element method near the boundary walls. The performance of the proposed partition-of-unity method is analyzed on several test examples for transient conduction–radiation problems in two space dimensions.

  5. A finite element analysis of shrinkage stress fields in two-phase materials

    E-Print Network [OSTI]

    Hall, Frank

    1980-01-01T23:59:59.000Z

    expansion were determined as constants. A simple finite element model was then analyzed to verify the suit- ability of the code for prediction of restrained shrinkage stress fields. Cree Com lienee Test A series of constant load isothermal creep tests... was judged to be suitable for prediction of shrinkage stress f ields. 0 m ~/0 ? Q ? Beam Bending Theory ? O ? Finite Element Method U o m m a 40 60 Time, min eo 100 120 Figure 8. Vertical deflection at midspan of four point bending test...

  6. On Some Versions of the Element Agglomeration AMGe Method

    SciTech Connect (OSTI)

    Lashuk, I; Vassilevski, P

    2007-08-09T23:59:59.000Z

    The present paper deals with element-based AMG methods that target linear systems of equations coming from finite element discretizations of elliptic PDEs. The individual element information (element matrices and element topology) is the main input to construct the AMG hierarchy. We study a number of variants of the spectral agglomerate element based AMG method. The core of the algorithms relies on element agglomeration utilizing the element topology (built recursively from fine to coarse levels). The actual selection of the coarse degrees of freedom (dofs) is based on solving large number of local eigenvalue problems. Additionally, we investigate strategies for adaptive AMG as well as multigrid cycles that are more expensive than the V-cycle utilizing simple interpolation matrices and nested conjugate gradient (CG) based recursive calls between the levels. The presented algorithms are illustrated with an extensive set of experiments based on a matlab implementation of the methods.

  7. Calibration under uncertainty for finite element models of masonry monuments

    SciTech Connect (OSTI)

    Atamturktur, Sezer,; Hemez, Francois,; Unal, Cetin

    2010-02-01T23:59:59.000Z

    Historical unreinforced masonry buildings often include features such as load bearing unreinforced masonry vaults and their supporting framework of piers, fill, buttresses, and walls. The masonry vaults of such buildings are among the most vulnerable structural components and certainly among the most challenging to analyze. The versatility of finite element (FE) analyses in incorporating various constitutive laws, as well as practically all geometric configurations, has resulted in the widespread use of the FE method for the analysis of complex unreinforced masonry structures over the last three decades. However, an FE model is only as accurate as its input parameters, and there are two fundamental challenges while defining FE model input parameters: (1) material properties and (2) support conditions. The difficulties in defining these two aspects of the FE model arise from the lack of knowledge in the common engineering understanding of masonry behavior. As a result, engineers are unable to define these FE model input parameters with certainty, and, inevitably, uncertainties are introduced to the FE model.

  8. Residual stresses in weld overlay tubes: A finite element study

    SciTech Connect (OSTI)

    Taljat, B.; Zacharia, T.; Wang, X.L.; Keiser, J.R. [Oak Ridge National Lab., TN (United States). Metals and Ceramics Div.; Feng, Z. [Edison Welding Inst., Columbus, OH (United States); Jirinec, M.J. [Welding Services, Inc., Norcross, GA (United States)

    1997-01-03T23:59:59.000Z

    Residual stresses and strains in a tube with circumferential weld overlay were analyzed by the finite element (FE) method. The objective of this work was to develop and verify a FE model, to determine the magnitude and distribution of residual stresses in the weld overlay tube, and to evaluate the significance of two contributing factors to residual stress: (1) difference in material properties between tube and weld material, and (2) thermal gradients in the weld. An axisymmetric FE model was developed to simulate the circumferential two-layer welding process of alloy 625 overlay on SA210 tube. The first layer was modeled as a gas metal arc welding process with filler metal, whereas the autogenous gas tungsten arc welding process was modeled for the second layer. Neutron diffraction technique was used to experimentally determine residual elastic strains in the weld overlay tube. Comparison with the FE results shows overall good agreement. Both the experimental and FE results show high compressive stresses at the inside tube surface and high tensile stresses in the weld overlay. This suggests that weld overlay may be used to relieve tensile or produce compressive stresses at the inside tube surface, which is significant for applications where crack initiation is found at the root pass of the joining weld.

  9. Computers & Geosciences 34 (2008) 10051014 Icetools: A full Stokes finite element model for glaciers$

    E-Print Network [OSTI]

    Frey, Pascal

    2008-01-01T23:59:59.000Z

    for glaciers$ A.H. Jarosch Institute of Earth Sciences, University of Iceland, Sturlugata 7, 101 Reykjavik Detailed studies of glacier flow and deformation of ice require high-resolution numerical modeling.40.Cy; 92.40.Sn; 92.40.Vq Keywords: Glacier; Finite element method; Full Stokes; Ice 1. Introduction

  10. Finite element analysis of a darcy velocimeter in a variably saturated soil

    SciTech Connect (OSTI)

    Bixler, N.E.; Carrigan, C.R.

    1989-01-01T23:59:59.000Z

    A new device, the Darcy velocimeter, is analyzed by the finite element method. The model is of two-phase flow and heat transport in a porous medium. Both laboratory-and field-scale experiments are simulated. The results show that the device should readily detect groundwater fluxes in the 10 to 100 m/yr range.

  11. Discontinuous finite volume element discretization for coupled flowtransport problems arising in models of sedimentation

    E-Print Network [OSTI]

    Bürger, Raimund

    in models of sedimentation Raimund B¨urgera , Sarvesh Kumarb , Ricardo Ruiz-Baierc aCI2MA and Departamento The sedimentation-consolidation and flow processes of a mixture of small particles dispersed in a vis- cous fluid discontinuous finite volume element (DFVE) discretization for this model is pro- posed. The numerical method

  12. Piezoelectric theory for finite element analysis of ultrasonic motors

    SciTech Connect (OSTI)

    Emery, J.D.; Mentesana, C.P.

    1997-06-01T23:59:59.000Z

    The authors present the fundamental equations of piezoelectricity and references. They show how a second form of the equations and a second set of coefficients can be found, through inversions involving the elasticity tensor. They show how to compute the clamped permittivity matrix from the unclamped matrix. The authors list the program pzansys.ftn and present examples of its use. This program does the conversions and calculations needed by the finite element program ANSYS.

  13. 3-D Finite Element Analyses of the Egan Cavern Field

    SciTech Connect (OSTI)

    Klamerus, E.W.; Ehgartner, B.L.

    1999-02-01T23:59:59.000Z

    Three-dimensional finite element analyses were performed for the two gas-filled storage caverns at the Egan field, Jennings dome, Louisiana. The effects of cavern enlargement on surface subsidence, storage loss, and cavern stability were investigated. The finite element model simulated the leaching of caverns to 6 and 8 billion cubic feet (BCF) and examined their performance at various operating conditions. Operating pressures varied from 0.15 psi/ft to 0.9 psi/ft at the bottom of the lowest cemented casing. The analysis also examined the stability of the web or pillar of salt between the caverns under differential pressure loadings. The 50-year simulations were performed using JAC3D, a three dimensional finite element analysis code for nonlinear quasistatic solids. A damage criterion based on onset of dilatancy was used to evaluate cavern instability. Dilation results from the development of microfractures in salt and, hence, potential increases in permeability onset occurs well before large scale failure. The analyses predicted stable caverns throughout the 50-year period for the range of pressures investigated. Some localized salt damage was predicted near the bottom walls of the caverns if the caverns are operated at minimum pressure for long periods of time. Volumetric cavern closures over time due to creep were moderate to excessive depending on the salt creep properties and operating pressures. However, subsidence above the cavern field was small and should pose no problem, to surface facilities.

  14. Finite Element Simulation of Piezoelectric Wafer Active Sensors for Structural Health Monitoring with Coupled-Filed

    E-Print Network [OSTI]

    Giurgiutiu, Victor

    Finite Element Simulation of Piezoelectric Wafer Active Sensors for Structural Health Monitoring) is emerging as an effective and powerful technique in structural health monitoring (SHM). Modeling to analytical calculation and experimental data. Key words: Structural Health Monitoring, PWAS, finite element

  15. An evaluation of finite element models of stiffened plates subjected to impulsive loading

    E-Print Network [OSTI]

    Pedatzur, Omri

    2004-01-01T23:59:59.000Z

    Different finite element models are evaluated for two very common structures, a cantilever beam and a stiffened plate, subjected to impulsive loading. For the cantilever beam case, the finite element models are one, two ...

  16. Implementation of finite element analysis into the athletic shoe design process

    E-Print Network [OSTI]

    Hidalgo, Maria E. (Maria Estela), 1982-

    2004-01-01T23:59:59.000Z

    Finite element analysis is used by companies throughout the world as a substitution for manually testing prototypes. With the assistance of finite element analysis many companies and industries have decreased the time and ...

  17. An enriched finite element model with q-refinement for radiative boundary layers in glass cooling

    SciTech Connect (OSTI)

    Mohamed, M. Shadi [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)] [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom); Seaid, Mohammed; Trevelyan, Jon [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom)] [School of Engineering and Computing Sciences, University of Durham, South Road, Durham DH1 3LE (United Kingdom); Laghrouche, Omar [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)] [Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS (United Kingdom)

    2014-02-01T23:59:59.000Z

    Radiative cooling in glass manufacturing is simulated using the partition of unity finite element method. The governing equations consist of a semi-linear transient heat equation for the temperature field and a stationary simplified P{sub 1} approximation for the radiation in non-grey semitransparent media. To integrate the coupled equations in time we consider a linearly implicit scheme in the finite element framework. A class of hyperbolic enrichment functions is proposed to resolve boundary layers near the enclosure walls. Using an industrial electromagnetic spectrum, the proposed method shows an immense reduction in the number of degrees of freedom required to achieve a certain accuracy compared to the conventional h-version finite element method. Furthermore the method shows a stable behaviour in treating the boundary layers which is shown by studying the solution close to the domain boundaries. The time integration choice is essential to implement a q-refinement procedure introduced in the current study. The enrichment is refined with respect to the steepness of the solution gradient near the domain boundary in the first few time steps and is shown to lead to a further significant reduction on top of what is already achieved with the enrichment. The performance of the proposed method is analysed for glass annealing in two enclosures where the simplified P{sub 1} approximation solution with the partition of unity method, the conventional finite element method and the finite difference method are compared to each other and to the full radiative heat transfer as well as the canonical Rosseland model.

  18. FINITE WAVELET DOMAIN METHOD FOR EFFICIENT MODELING OF LAMB WAVE BASED STRUCTURAL HEALTH MONITORING

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    FINITE WAVELET DOMAIN METHOD FOR EFFICIENT MODELING OF LAMB WAVE BASED STRUCTURAL HEALTH MONITORING element exploits the advantages of wavelets for the spatial discretization of the displacement field : Wavelet-based Finite Element, Transient Response. 1 INTRODUCTION The dynamic transient response and wave

  19. Element-based concrete design with three-dimensional finite element models

    SciTech Connect (OSTI)

    O'Leary, M.; Huberty, K.; Winch, S. [Nuclear Power Technologies Div., Sargent and Lundy, 55 East Monroe, Chicago, IL 60603 (United States)

    2012-07-01T23:59:59.000Z

    A shell element based design of a typical shear wall using analytical results from a three-dimensional finite element model subjected to a combination of vertical and lateral loads is evaluated. The axial and flexural force resultants from each element for every load combination are used to calculate the required reinforcing for each element. Strength for axial loads (P) and out-of-plane flexure (M) in structural walls is determined according to the same P-M interaction procedures used for columns. After each element has been evaluated, a required reinforcing map for each face of each element in the wall is presented along with a constructible reinforcement pattern enveloping the required reinforcing. In order to determine whether the element-based approach meets the requirements of the section cut approach to design, which is typically employed in manual calculations, the total in-plane moment (M) and total vertical axial force (P) across the entire length of the wall is calculated and the P-M points are plotted on an in-plane P-M interaction diagram. It is concluded that element-based design for a structural wall ensures that reinforcement is provided where required by the three-dimensional finite element analysis while still providing sufficient reinforcing to satisfy the section cut approach to design. (authors)

  20. Detecting the Causes of Ill-Conditioning in Structural Finite Element Models

    E-Print Network [OSTI]

    Higham, Nicholas J.

    .6 of the finite element-based structural analysis package Oasys GSA was released. A new feature in this release

  1. Detecting the causes of ill-conditioning in structural finite element Ramaseshan Kannan b,a

    E-Print Network [OSTI]

    Tisseur, Francoise

    .6 of the finite element-based structural analysis package Oasys GSA was released. A new feature in this release

  2. Finite element analyses for seismic shear wall international standard problem

    SciTech Connect (OSTI)

    Park, Y.J.; Hofmayer, C.H.

    1998-04-01T23:59:59.000Z

    Two identical reinforced concrete (RC) shear walls, which consist of web, flanges and massive top and bottom slabs, were tested up to ultimate failure under earthquake motions at the Nuclear Power Engineering Corporation`s (NUPEC) Tadotsu Engineering Laboratory, Japan. NUPEC provided the dynamic test results to the OECD (Organization for Economic Cooperation and Development), Nuclear Energy Agency (NEA) for use as an International Standard Problem (ISP). The shear walls were intended to be part of a typical reactor building. One of the major objectives of the Seismic Shear Wall ISP (SSWISP) was to evaluate various seismic analysis methods for concrete structures used for design and seismic margin assessment. It also offered a unique opportunity to assess the state-of-the-art in nonlinear dynamic analysis of reinforced concrete shear wall structures under severe earthquake loadings. As a participant of the SSWISP workshops, Brookhaven National Laboratory (BNL) performed finite element analyses under the sponsorship of the U.S. Nuclear Regulatory Commission (USNRC). Three types of analysis were performed, i.e., monotonic static (push-over), cyclic static and dynamic analyses. Additional monotonic static analyses were performed by two consultants, F. Vecchio of the University of Toronto (UT) and F. Filippou of the University of California at Berkeley (UCB). The analysis results by BNL and the consultants were presented during the second workshop in Yokohama, Japan in 1996. A total of 55 analyses were presented during the workshop by 30 participants from 11 different countries. The major findings on the presented analysis methods, as well as engineering insights regarding the applicability and reliability of the FEM codes are described in detail in this report. 16 refs., 60 figs., 16 tabs.

  3. Domain decomposition for coupling finite and boundary element methods in EEG Emmanuel Olivi1, Maureen Clerc1 and Theodore Papadopoulo1

    E-Print Network [OSTI]

    Boyer, Edmond

    in the patient head. Yet, the skull anisotropy happens to be highly anisotropic, and must then be modeled.e. for the brain and the scalp). A domain decomposition (DD) framework allows to split the global system. This work presents such a coupling formulation of a 3-DD method solving iteratively a BEM for the brain

  4. Finite Element Modeling of the Deformation of a Thin Magnetoelastic Film Compared to a Membrane Model

    SciTech Connect (OSTI)

    Barham, M; White, D; Steigmann, D; Rudd, R

    2009-04-08T23:59:59.000Z

    Recently a new class of biocompatible elastic polymers loaded with small ferrous particles (magnetoelastomer) was developed at Lawrence Livermore National Laboratory. This new material was formed as a thin film using spin casting. The deformation of this material using a magnetic field has many possible applications to microfluidics. Two methods will be used to calculate the deformation of a circular magneto-elastomeric film subjected to a magnetic field. The first method is an arbitrary Lagrangian-Eulerian (ALE) finite element method (FEM) and the second is based on nonlinear continuum electromagnetism and continuum elasticity in the membrane limit. The comparison of these two methods is used to test/validate the finite element method.

  5. Wear 258 (2005) 17871793 Finite element analysis and experiments of metal/metal

    E-Print Network [OSTI]

    Sawyer, Wallace

    2005-01-01T23:59:59.000Z

    Wear 258 (2005) 1787­1793 Finite element analysis and experiments of metal/metal wear-on-disk experiments. The results from the finite element analysis were in close agreement with the block of this simulation. © 2004 Published by Elsevier B.V. Keywords: Wear modeling; Finite element analysis 1

  6. Wrinkle development analysis in thin sail-like structures using MITC shell finite elements

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Wrinkle development analysis in thin sail-like structures using MITC shell finite elements D submitted to Finite Elements in Analysis and Design September 20, 2012 hal-00733994,version1-20Sep2012 Author manuscript, published in "Finite Elements in Analysis and Design 64 (2013) 48-64" DOI : 10.1016/j

  7. Tolerance Analysis of Flexible Assemblies Using Finite Element and Spectral Analysis

    E-Print Network [OSTI]

    Tolerance Analysis of Flexible Assemblies Using Finite Element and Spectral Analysis ADCATS Report the autocorrelation function from frequency spectrum analysis to model random surface variations. Finite element BACKGROUND LITERATURE REVIEW 5 3 STATISTICAL TOLERANCE ANALYSIS USING FINITE ELEMENT ANALYSIS 3.1 Assumptions

  8. Finite element analysis of performance in the skulls of marmosets and tamarins

    E-Print Network [OSTI]

    Dumont, Elizabeth R.

    Finite element analysis of performance in the skulls of marmosets and tamarins Elizabeth R. Dumont compara- tive analyses. Key words finite element analysis; marmoset; performance; tamarin. Introduction. Using 3D finite element models of the marmoset Callithrix jacchus and the tamarin Saguinus fuscicollis

  9. Approximation Theory and Finite Element Analysis January 2014 Lecturer David Silvester

    E-Print Network [OSTI]

    Silvester, David J.

    Approximation Theory and Finite Element Analysis January 2014 Lecturer David Silvester Office Alan Analysis, CUP, 2003. · Dietrich Braess, Finite Elements: Theory, Fast Solvers and Applications in Solid Mechanics, CUP, Third Edition, 2007. · Howard Elman, David Silvester and Andy Wathen, Finite Elements

  10. Free Body Analysis, Beam Mechanics, and Finite Element Modeling of the Mandible of Alligator

    E-Print Network [OSTI]

    Free Body Analysis, Beam Mechanics, and Finite Element Modeling of the Mandible of Alligator arm mechanics, 2D and 3D beam mod- els, and three high-resolution finite element models- pared with the beam models, the Alligator finite element models exhibited less spatial variability

  11. The Finite Element Analysis of Quartz Crystal Resonators with Nonlinear Plate Equations

    E-Print Network [OSTI]

    Wang, Ji

    The Finite Element Analysis of Quartz Crystal Resonators with Nonlinear Plate Equations Ji Wang1, dujianke, wanglihong}@nbu.edu.cn, yyong@rci.rutgers.edu Abstract--The finite element analysis of quartz for typical thickness-shear vibrations of quartz crystal plates. While linear finite element analysis

  12. Finite Element Analysis of Nonlinear Thickness-shear Vibrations of AT-cut Quartz Crystal Plates

    E-Print Network [OSTI]

    Wang, Ji

    Finite Element Analysis of Nonlinear Thickness-shear Vibrations of AT-cut Quartz Crystal Plates Ji, dujianke}@nbu.edu.cn Abstract--The nonlinear finite element analysis is performed with the nonlinear a smaller size in comparison with the 3D approach. General procedure of nonlinear finite element analysis

  13. FINITE ELEMENT ANALYSIS OF A CLASS OF STRESS-FREE MARTENSITIC MICROSTRUCTURES

    E-Print Network [OSTI]

    Li, Bo

    FINITE ELEMENT ANALYSIS OF A CLASS OF STRESS-FREE MARTENSITIC MICROSTRUCTURES BO LI Abstract. This work is concerned with the finite element approximation of a class of stress-free martensitic microstructures modeled by multi-well en- ergy minimization. Finite element energy-minimizing sequences are first

  14. Comparison of Beam Theory and Finite-Element Analysis With In Vivo

    E-Print Network [OSTI]

    Comparison of Beam Theory and Finite-Element Analysis With In Vivo Bone Strain Data From structures and, more recently, finite- element (FE) analysis. In this study, we compare experimentally-Liss, Inc. Key words: bone strain; finite-element analysis; skull biome- chanics; alligator; feeding

  15. Finite Elements in Analysis and Design 41 (2004) 253265 www.elsevier.com/locate/finel

    E-Print Network [OSTI]

    Rivero, Jose Pablo Suárez

    2004-01-01T23:59:59.000Z

    Finite Elements in Analysis and Design 41 (2004) 253­265 www.elsevier.com/locate/finel The 8 favorable in finite element analysis. Furthermore, since a right prism can be subdivided into six right adaptive finite element analysis procedure has received much attention

  16. A finite element inverse analysis to assess functional improvement during the fracture healing process

    E-Print Network [OSTI]

    Miga, Michael I.

    A finite element inverse analysis to assess functional improvement during the fracture healing architecture on the FEA estimated material property metric. The finite element model inverse analysis developed i n f o Article history: Accepted 2 September 2009 Keywords: Fracture healing Finite element

  17. A Finite Element Analysis of Masticatory Stress Janine Chalk,1,2

    E-Print Network [OSTI]

    A Finite Element Analysis of Masticatory Stress Hypotheses Janine Chalk,1,2 * Brian G. Richmond,2). This study uses finite element analysis to examine the extent to which these geometric models provide loading regimes, were applied to a previously validated finite element model (FEM) of a macaque (Macaca

  18. Finite Element Analysis of Ultrasonic Processing of a Polymer-Matrix Composite

    E-Print Network [OSTI]

    Roylance, David

    Finite Element Analysis of Ultrasonic Processing of a Polymer- Matrix Composite Wenguang Zhao1 parameters, makes trial-and-error analysis of the process difficult, and this paper demonstrates how a finite-element; Ultrasonic Tape Lamination; Viscoelastic dissipation; Finite-element analysis 1. Introduction This paper

  19. NONLINEAR FINITE ELEMENT ANALYSIS OF REINFORCED CONCRETE STRUCTURES SUBJECTED TO TRANSIENT THERMAL LOADS

    E-Print Network [OSTI]

    Vecchio, Frank J.

    Abstract ii NONLINEAR FINITE ELEMENT ANALYSIS OF REINFORCED CONCRETE STRUCTURES SUBJECTED nonlinear finite element transient thermal analysis is implemented into program VecTor2©, a nonlinear's response to thermal loads, this research develops and implements a 2D nonlinear finite element transient

  20. Author's personal copy The effects of modeling simplifications on craniofacial finite element models

    E-Print Network [OSTI]

    a r t i c l e i n f o Article history: Accepted 17 March 2011 Keywords: Finite element analysis Tooth), a biomechanical system commonly studied using finite element analysis (FEA) is that part of the craniofacial in such finite element models (FEMs) depends exten- sively on the objective of the analysis and the performance

  1. Finite Element Multi-physics Modeling for Ohmic Contact of Microswitches , D. Leray1, 2

    E-Print Network [OSTI]

    Boyer, Edmond

    Finite Element Multi-physics Modeling for Ohmic Contact of Microswitches H.Liu1, 2 , D. Leray1, 2-physics modelling of electrical contact is accomplished with the finite element commercial package ANSYSTM . Two behaviour of electrical contact with rough surface included. Keywords: multi-physics, finite element

  2. Applied Finite Element Technology (3 Credits) Instructor Kamran Behdinan, MIE Dept, University of Toronto; behdinan@mie.utoronto.ca

    E-Print Network [OSTI]

    Chen, Zheng

    -F, July 7­25, 2014; Final Exam: 1-4 PM, July 26, 2014 Objective To develop an understanding of finite. References 1. D. L. Logan, The First Course in the Finite Element Method, 5th Edition, SI; Cengage Learning Project Assessment (20%) 30% Midterm Exam 10% Final Exam 40% Total 100% #12;

  3. Guaranteed Verification of Finite Element Solutions of Heat Conduction

    E-Print Network [OSTI]

    Wang, Delin

    2012-07-16T23:59:59.000Z

    Engineering iii ABSTRACT Guaranteed Verification of Finite Element Solutions of Heat Conduction. (May 2011) Delin Wang, B.E., Qingdao University of Science & Technology, China; M.S., Jilin University, China; M.E., Texas A&M University Chair of Advisory... level and ? ky kx is the characteristic thickness of the boundary layer. . . . . 37 3.5 Model problem with interface layer. The relative value of the energy norm of the error ErelU = ||eSp?h ||U / ||uEX||U? 100% versus #15; for n = 1, 2, 3...

  4. Visualizing Higher Order Finite Elements: FY05 Yearly Report.

    SciTech Connect (OSTI)

    Thompson, David; Pebay, Philippe Pierre

    2005-11-01T23:59:59.000Z

    This report contains an algorithm for decomposing higher-order finite elementsinto regions appropriate for isosurfacing and proves the conditions under which thealgorithm will terminate. Finite elements are used to create piecewise polynomialapproximants to the solution of partial differential equations for which no analyticalsolution exists. These polynomials represent fields such as pressure, stress, and mo-mentim. In the past, these polynomials have been linear in each parametric coordinate.Each polynomial coefficient must be uniquely determined by a simulation, and thesecoefficients are called degrees of freedom. When there are not enough degrees of free-dom, simulations will typically fail to produce a valid approximation to the solution.Recent work has shown that increasing the number of degrees of freedom by increas-ing the order of the polynomial approximation (instead of increasing the number offinite elements, each of which has its own set of coefficients) can allow some typesof simulations to produce a valid approximation with many fewer degrees of freedomthan increasing the number of finite elements alone. However, once the simulation hasdetermined the values of all the coefficients in a higher-order approximant, tools donot exist for visual inspection of the solution.This report focuses on a technique for the visual inspection of higher-order finiteelement simulation results based on decomposing each finite element into simplicialregions where existing visualization algorithms such as isosurfacing will work. Therequirements of the isosurfacing algorithm are enumerated and related to the placeswhere the partial derivatives of the polynomial become zero. The original isosurfacingalgorithm is then applied to each of these regions in turn.3 AcknowledgementThe authors would like to thank David Day and Louis Romero for their insight into poly-nomial system solvers and the LDRD Senior Council for the opportunity to pursue thisresearch. The authors were supported by the United States Department of Energy, Officeof Defense Programs by the Labratory Directed Research and Development Senior Coun-cil, project 90499. Sandia is a multiprogram laboratory operated by Sandia Corporation,a Lockheed-Martin Company, for the United States Department of Energy under contractDE-AC04-94-AL85000.4

  5. Current commutator anomalies in finite-element quantum electrodynamics

    E-Print Network [OSTI]

    Dean F. Miller

    1993-12-20T23:59:59.000Z

    Four-dimensional quantum electrodynamics has been formulated on a hypercubic Minkowski finite-element lattice. The equations of motion have been derived so as to preserve lattice gauge invariance and have been shown to be unitary. In addition, species doubling is avoided due to the nonlocality of the interactions. The model is used to investigate the lattice current algebra. Regularization of the current is shown to arise in a natural and nonarbitrary way. The commutators of the lattice current are calculated and shown to have the expected qualitative behavior. These lattice results are compared to various continuum calculations. (Five figures available from author.)

  6. Transport and dispersion of pollutants in surface impoundments: a finite element model

    SciTech Connect (OSTI)

    Yeh, G.T.

    1980-07-01T23:59:59.000Z

    A surface impoundment model in finite element (SIMFE) is presented to enable the simulation of flow circulations and pollutant transport and dispersion in natural or artificial lakes, reservoirs or ponds with any number of islands. This surface impoundment model consists of two sub-models: hydrodynamic and pollutant transport models. Both submodels are simulated by the finite element method. While the hydrodynamic model is solved by the standard Galerkin finite element scheme, the pollutant transport model can be solved by any of the twelve optional finite element schemes built in the program. Theoretical approximations and the numerical algorithm of SIMFE are described. Detail instruction of the application are given and listing of FORTRAN IV source program are provided. Two sample problems are given. One is for an idealized system with a known solution to show the accuracy and partial validation of the models. The other is applied to Prairie Island for a set of hypothetical input data, typifying a class of problems to which SIMFE may be applied.

  7. Generalized finite element method for Helmholtz equation

    E-Print Network [OSTI]

    Hidajat, Realino Lulie

    2009-05-15T23:59:59.000Z

    ,0) and its estimator E h;p,q u(?1,0) (d) computed using d/R = 1, with the truncation radius R =1.5. The di?erence between the curve for the Robin and the curves for the SOM, BT1, BT2 boundary conditions is the error due to the artificial boundary condition... results as the ones obtained using the Robin boundary condition. 94 5.9 Location of output points, the domain ? with truncation radius R =2.5 and the employed mesh. ............................ 94 5.10 q-convergence of e h;p,q u(A) using Sommerfeld...

  8. Control-volume mixed finite element methods

    E-Print Network [OSTI]

    1998-01-23T23:59:59.000Z

    for miscible displacement [15]. This theory extends in a ...... a scalar multiple of the tensor, so that the above double integrals can be stored and later multiplied by a ..... 3, Project Report, Reservoir Simulation Research. Corporation (1995).

  9. CONVERGENCE OF A MULTISCALE FINITE ELEMENT METHOD ...

    E-Print Network [OSTI]

    1999-04-17T23:59:59.000Z

    Mar 3, 1999 ... In the above proof, ul is used as an intermediate step towards ...... Current address: Exxon Production Research Company, P. O. Box 2189, Houston, TX 77252 ... Department of Mathematics, Purdue University, West Lafayette, ...

  10. Fracture and Fragmentation of Simplicial Finite Elements Meshes using Graphs

    SciTech Connect (OSTI)

    Mota, A; Knap, J; Ortiz, M

    2006-10-18T23:59:59.000Z

    An approach for the topological representation of simplicial finite element meshes as graphs is presented. It is shown that by using a graph, the topological changes induced by fracture reduce to a few, local kernel operations. The performance of the graph representation is demonstrated and analyzed, using as reference the 3D fracture algorithm by Pandolfi and Ortiz [22]. It is shown that the graph representation initializes in O(N{sub E}{sup 1.1}) time and fractures in O(N{sub I}{sup 1.0}) time, while the reference implementation requires O(N{sub E}{sup 2.1}) time to initialize and O(N{sub I}{sup 1.9}) time to fracture, where NE is the number of elements in the mesh and N{sub I} is the number of interfaces to fracture.

  11. Finite Element Analysis of Warpage in Laminated Aluminium Alloy Plates for Machining of Primary Aeronautic Parts

    SciTech Connect (OSTI)

    Reis, A. C.; Moreira Filho, L. A. [ITA - Aeronautic Technological Institute, IEM, Praca Mal. Eduardo Gomes, 50 - Vila das Acacias - S. J. Campos, CEP 1228-900 (Brazil); Menezes, M. A. [ITA - Aeronautic Technological Institute, IEM, Praca Mal. Eduardo Gomes, 50 - Vila das Acacias - S. J. Campos, CEP 1228-900 (Brazil); UNESP - Sao Paulo State University, Ilha Solteira Engineering Faculty, DEM, Av. Brasil, 56 - Centro - Ilha Solteira - SP, CEP 15385-000 (Brazil)

    2007-04-07T23:59:59.000Z

    The aim of this paper consists in presenting a method of simulating the warpage in 7xxx series aluminium alloy plates. To perform this simulation finite element software MSC.Patran and MSC.Marc were used. Another result of this analysis will be the influence on material residual stresses induced on the raw material during the rolling process upon the warpage of primary aeronautic parts, fabricated through machining (milling) at Embraer. The method used to determinate the aluminium plate residual stress was Layer Removal Test. The numerical algorithm Modified Flavenot Method was used to convert layer removal and beam deflection in stress level. With such information about the level and profile of residual stresses become possible, during the step that anticipate the manufacturing to incorporate these values in the finite-element approach for modelling warpage parts. Based on that warpage parameter surely the products are manufactured with low relative vulnerability propitiating competitiveness and price.

  12. Fidelity of a Finite Element Model for Longitudinal Wave Propagation in Thick Cylindrical Wave Guides

    SciTech Connect (OSTI)

    Puckett, A.D.

    2000-09-01T23:59:59.000Z

    The ability to model wave propagation in circular cylindrical bars of finite length numerically or analytically has many applications. In this thesis the capability of an explicit finite element method to model longitudinal waves in cylindrical rods with circular cross-sections is explored. Dispersion curves for the first four modes are compared to the analytical solution to determine the accuracy of various element sizes and time steps. Values for the time step and element size are determined that retain accuracy while minimizing computational time. The modeling parameters are validated by calculating a signal propagated with a broadband input force. Limitations on the applicability are considered along with modeling parameters that should be applicable to more general geometries.

  13. Large-Scale Parallel Finite Element Analysis of the Stress Singular Problems

    SciTech Connect (OSTI)

    Noriyuki Kushida; Hiroshi Okuda; Genki Yagawa [University of Tokyo (Japan)

    2002-07-01T23:59:59.000Z

    In this paper, the convergence behavior of large-scale parallel finite element method for the stress singular problems was investigated. The convergence behavior of iterative solvers depends on the efficiency of the pre-conditioners. However, efficiency of pre-conditioners may be influenced by the domain decomposition that is necessary for parallel FEM. In this study the following results were obtained: Conjugate gradient method without preconditioning and the diagonal scaling preconditioned conjugate gradient method were not influenced by the domain decomposition as expected. symmetric successive over relaxation method preconditioned conjugate gradient method converged 6% faster as maximum if the stress singular area was contained in one sub-domain. (authors)

  14. Finite element visualization in the cave virtual reality environment

    SciTech Connect (OSTI)

    Plaskacz, E.J. [Argonne National Lab., IL (United States); Kuhn, M.A. [Northern Illinois Univ., De Kalb, IL (United States)

    1996-03-01T23:59:59.000Z

    Through the use of the post-processing software, Virtual Reality visualization (VRviz), and the Cave Automatic Virtual Environment (CAVE), finite element representations can be viewed as they would be in real life. VRviz is a program written in ANSI C to translate the mathematical results generated by finite element analysis programs into a virtual representation. This virtual representation is projected into the CAVE environment and the results are animated. The animation is fully controllable. A user is able to translate the image, rotate about any axis and scale the image at any time. The user is also able to freeze the animation at any time step and control the image update rate. This allows the user to navigate around, or even inside, the image in order to effectively analyze possible failure points and redesign as necessary. Through the use of the CAVE and the real life image that is being produced by VRviz, engineers are able to save considerable time, money, and effort in the design process.

  15. A Finite Element Model for Simulation of Carbon Dioxide Sequestration

    SciTech Connect (OSTI)

    Bao, Jie; Xu, Zhijie; Fang, Yilin

    2013-11-02T23:59:59.000Z

    We present a hydro-mechanical model, followed by stress, deformation, and shear-slip failure analysis for geological sequestration of carbon dioxide (CO2). The model considers the poroelastic effects by taking into account of the two-way coupling between the geomechanical response and the fluid flow process. Analytical solutions for pressure and deformation fields were derived for a typical geological sequestration scenario in our previous work. A finite element approach is introduced here for numerically solving the hydro-mechanical model with arbitrary boundary conditions. The numerical approach was built on an open-source finite element code Elmer, and results were compared to the analytical solutions. The shear-slip failure analysis was presented based on the numerical results, where the potential failure zone is identified. Information is relevant to the prediction of the maximum sustainable injection rate or pressure. The effects of caprock permeability on the fluid pressure, deformation, stress, and the shear-slip failure zone were also quantitatively studied. It was shown that a larger permeability in caprock and base rock leads to a larger uplift but a smaller shear-slip failure zone.

  16. Finite-Element Analysis of Biting Behavior and Bone Stress in the

    E-Print Network [OSTI]

    Dumont, Elizabeth R.

    Finite-Element Analysis of Biting Behavior and Bone Stress in the Facial Skeletons of Bats biting behavior and bite force data gathered in the field with finite-element (FE) analysis. Our FE words: biting behavior; bone stress; adaptation; finite-ele- ment analysis; Chiroptera Mammal evolution

  17. The Uranium Processing Facility Finite Element Meshing Discussion

    Office of Environmental Management (EM)

    more nodes and elements (1 element 9 elements) - Changed from GTStrudl to SAP 2000. - Linear elastic computer model * DOE Project 2 - 15'-20' element size deemed...

  18. FINITE ELEMENT ANALYSIS OF BULK TRITIUM SHIPPING PACKAGE

    SciTech Connect (OSTI)

    Jordan, J.

    2010-06-02T23:59:59.000Z

    The Bulk Tritium Shipping Package was designed by Savannah River National Laboratory. This package will be used to transport tritium. As part of the requirements for certification, the package must be shown to meet the scenarios of the Hypothetical Accident Conditions (HAC) defined in Code of Federal Regulations Title 10 Part 71 (10CFR71). The conditions include a sequential 30-foot drop event, 30-foot dynamic crush event, and a 40-inch puncture event. Finite Element analyses were performed to support and expand upon prototype testing. Cases similar to the tests were evaluated. Additional temperatures and orientations were also examined to determine their impact on the results. The peak stress on the package was shown to be acceptable. In addition, the strain on the outer drum as well as the inner containment boundary was shown to be acceptable. In conjunction with the prototype tests, the package was shown to meet its confinement requirements.

  19. Finite element analyses of a linear-accelerator electron gun

    SciTech Connect (OSTI)

    Iqbal, M., E-mail: muniqbal.chep@pu.edu.pk, E-mail: muniqbal@ihep.ac.cn [Centre for High Energy Physics, University of the Punjab, Lahore 45590 (Pakistan); Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049 (China); Wasy, A. [Department of Mechanical Engineering, Changwon National University, Changwon 641773 (Korea, Republic of)] [Department of Mechanical Engineering, Changwon National University, Changwon 641773 (Korea, Republic of); Islam, G. U. [Centre for High Energy Physics, University of the Punjab, Lahore 45590 (Pakistan)] [Centre for High Energy Physics, University of the Punjab, Lahore 45590 (Pakistan); Zhou, Z. [Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049 (China)] [Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049 (China)

    2014-02-15T23:59:59.000Z

    Thermo-structural analyses of the Beijing Electron-Positron Collider (BEPCII) linear-accelerator, electron gun, were performed for the gun operating with the cathode at 1000?°C. The gun was modeled in computer aided three-dimensional interactive application for finite element analyses through ANSYS workbench. This was followed by simulations using the SLAC electron beam trajectory program EGUN for beam optics analyses. The simulations were compared with experimental results of the assembly to verify its beam parameters under the same boundary conditions. Simulation and test results were found to be in good agreement and hence confirmed the design parameters under the defined operating temperature. The gun is operating continuously since commissioning without any thermal induced failures for the BEPCII linear accelerator.

  20. TOPAZ3D. 3-D Finite Element Heat Transfer

    SciTech Connect (OSTI)

    Shapiro, A.B. [Lawrence Livermore National Lab., CA (United States)

    1992-02-24T23:59:59.000Z

    TOPAZ3D is a three-dimensional implicit finite element computer code for heat transfer analysis. TOPAZ3D can be used to solve for the steady-state or transient temperature field on three-dimensional geometries. Material properties may be temperature-dependent and either isotropic or orthotropic. A variety of time-dependent and temperature-dependent boundary conditions can be specified including temperature, flux, convection, and radiation. By implementing the user subroutine feature, users can model chemical reaction kinetics and allow for any type of functional representation of boundary conditions and internal heat generation. TOPAZ3D can solve problems of diffuse and specular band radiation in an enclosure coupled with conduction in the material surrounding the enclosure. Additional features include thermal contact resistance across an interface, bulk fluids, phase change, and energy balances.

  1. Numerical solution of 2-D nonlinear inverse heat conduction problems using finite-element techniques

    SciTech Connect (OSTI)

    Khachfe, R.A.; Jarny, Y.

    2000-02-01T23:59:59.000Z

    A general method is presented for solving different classes of nonlinear inverse heat conduction problems (IHCP) for two-dimensional, arbitrarily shaped bodies. It is based on the systematic use of a finite-element library. It is shown that, following this approach, the conjugate gradient method can be easily implemented. The method offers a very wide field of practical applications in inverse thermal analysis, while reducing very significantly the amount of work which remains specific for each particular IHPC. Two numerical experiments illustrate the influence of data errors and the iterative regularization principle.

  2. Structural Analysis of a New Design of a Glenoid Prosthesis by Finite Element Modeling

    E-Print Network [OSTI]

    Downey, John A

    1990-01-01T23:59:59.000Z

    18 Abstract Prosthesis loosening in total shoulderexcision. Present glenoid prostheses have fixation only inDesign of a Glenoid Prosthesis by Finite Element Modeling By

  3. Materials Science & Metallurgy Master of Philosophy, Materials Modelling, Course MP7, Finite Element Analysis, H. K. D. H. Bhadeshia

    E-Print Network [OSTI]

    Cambridge, University of

    Element Analysis, H. K. D. H. Bhadeshia Lecture 3: Finite Elements Steady­state heat flow through A and length L. The finite element representation consists of two nodes i and j. Heat flow in one. The temperatures T1 and T2 are maintained constant. The finite element representation consists of a single element

  4. A FINITE ELEMENT APPROACH FOR SIMULATING THE QUENCHING OF LARGE ALUMINIUM COMPRESSOR WHEELS

    E-Print Network [OSTI]

    A FINITE ELEMENT APPROACH FOR SIMULATING THE QUENCHING OF LARGE ALUMINIUM COMPRESSOR WHEELS A FINITE ELEMENT APPROACH FOR SIMULATING THE QUENCHING OF LARGE ALUMINIUM COMPRESSOR WHEELS Dirk Büche Aluminum devices such as compressor wheels are typically heat treated in order to achieve the required

  5. Equal-order finite elements with local projection stabilization for the Darcy-Brinkman equations

    E-Print Network [OSTI]

    Schieweck, Friedhelm

    Equal-order finite elements with local projection stabilization for the Darcy-Brinkman equations M. Braack and F. Schieweck Mai 2010 Abstract For the Darcy-Brinkman equations, which model porous media: Porous media flow, Darcy-Brinkman equations, Stokes, equal-order finite elements, local projection

  6. STATISTICAL AND 3D NONLINEAR FINITE ELEMENT ANALYSIS OF SCHLEGEIS DAM

    E-Print Network [OSTI]

    Balaji, Rajagopalan

    STATISTICAL AND 3D NONLINEAR FINITE ELEMENT ANALYSIS OF SCHLEGEIS DAM VICTOR SAOUMA, ERIC HANSEN is composed of two parts. First a statistical analysis of the dam crest displacement is performed, along with a prediction for the years 2000-2001. Then a 3D finite element analysis of Schlegeis dam is performed using

  7. Prediction of a Mysticete Audiogram via Finite Element Analysis of the Middle Ear

    E-Print Network [OSTI]

    Prediction of a Mysticete Audiogram via Finite Element Analysis of the Middle Ear Andrew Tubelli (Rosowski 1994). Biophysical models of the cetacean middle ear can be developed using finite element (FE for FE analysis (COMSOL, Stockholm, Sweden). _________________________ A. Tubelli (!) · A. Zosuls · D

  8. Solar Electric Generating System II finite element analysis

    SciTech Connect (OSTI)

    Dohner, J.L.; Anderson, J.R.

    1994-04-01T23:59:59.000Z

    On June 2, 1992, Landers` earthquake struck the Solar Electric Generating System II, located in Daggett, California. The 30 megawatt power station, operated by the Daggett Leasing Corporation (DLC), suffered substantial damage due to structural failures in the solar farm. These failures consisted of the separation of sliding joints supporting a distribution of parabolic glass mirrors. At separation, the mirrors fell to the ground and broke. It was the desire of the DLC and the Solar Thermal Design Assistance Center (STDAC) of Sandia National Laboratories (SNL) and to redesign these joints so that, in the event of future quakes, costly breakage will be avoided. To accomplish this task, drawings of collector components were developed by the STDAC, from which a detailed finite element computer model of a solar collector was produced. This nonlinear dynamic model, which consisted of over 8,560 degrees of freedom, underwent model reduction to form a low order nonlinear dynamic model containing only 40 degrees of freedom. This model was then used as a design tool to estimate joint dynamics. Using this design tool, joint configurations were modified, and an acceptable joint redesign determined. The results of this analysis showed that the implementation of metal stops welded to support shafts for the purpose of preventing joint separation is a suitable joint redesign. Moreover, it was found that, for quakes of Landers` magnitude, mirror breakage due to enhanced vibration in the trough assembly is unlikely.

  9. A finite element simulation of biological conversion processes in landfills

    SciTech Connect (OSTI)

    Robeck, M., E-mail: markus.robeck@uni-due.de [Department of Water and Waste Management, Building Sciences, University of Duisburg-Essen, Universitaetsstrasse 15, 45141 Essen (Germany); Ricken, T. [Institute of Mechanics/Computational Mechanics, Building Sciences, University of Duisburg-Essen, Universitaetsstrasse 15, 45141 Essen (Germany); Widmann, R. [Department of Water and Waste Management, Building Sciences, University of Duisburg-Essen, Universitaetsstrasse 15, 45141 Essen (Germany)

    2011-04-15T23:59:59.000Z

    Landfills are the most common way of waste disposal worldwide. Biological processes convert the organic material into an environmentally harmful landfill gas, which has an impact on the greenhouse effect. After the depositing of waste has been stopped, current conversion processes continue and emissions last for several decades and even up to 100 years and longer. A good prediction of these processes is of high importance for landfill operators as well as for authorities, but suitable models for a realistic description of landfill processes are rather poor. In order to take the strong coupled conversion processes into account, a constitutive three-dimensional model based on the multiphase Theory of Porous Media (TPM) has been developed at the University of Duisburg-Essen. The theoretical formulations are implemented in the finite element code FEAP. With the presented calculation concept we are able to simulate the coupled processes that occur in an actual landfill. The model's theoretical background and the results of the simulations as well as the meantime successfully performed simulation of a real landfill body will be shown in the following.

  10. FEMA: a Finite Element Model of Material Transport through Aquifers

    SciTech Connect (OSTI)

    Yeh, G.T.; Huff, D.D.

    1985-01-01T23:59:59.000Z

    This report documents the construction, verification, and demonstration of a Finite Element Model of Material Transport through Aquifers (FEMA). The particular features of FEMA are its versatility and flexibility to deal with as many real-world problems as possible. Mechanisms included in FEMA are: carrier fluid advection, hydrodynamic dispersion and molecular diffusion, radioactive decay, sorption, source/sinks, and degradation due to biological, chemical as well as physical processes. Three optional sorption models are embodied in FEMA. These are linear isotherm and Freundlich and Langmuir nonlinear isotherms. Point as well as distributed source/sinks are included to represent artificial injection/withdrawals and natural infiltration of precipitation. All source/sinks can be transient or steady state. Prescribed concentration on the Dirichlet boundary, given gradient on the Neumann boundary segment, and flux at each Cauchy boundary segment can vary independently of each other. The aquifer may consist of as many formations as desired. Either completely confined or completely unconfined or partially confined and partially unconfined aquifers can be dealt with effectively. FEMA also includes transient leakage to or from the aquifer of interest through confining beds from or to aquifers lying below and/or above.

  11. Finite element meshing approached as a global minimization process

    SciTech Connect (OSTI)

    WITKOWSKI,WALTER R.; JUNG,JOSEPH; DOHRMANN,CLARK R.; LEUNG,VITUS J.

    2000-03-01T23:59:59.000Z

    The ability to generate a suitable finite element mesh in an automatic fashion is becoming the key to being able to automate the entire engineering analysis process. However, placing an all-hexahedron mesh in a general three-dimensional body continues to be an elusive goal. The approach investigated in this research is fundamentally different from any other that is known of by the authors. A physical analogy viewpoint is used to formulate the actual meshing problem which constructs a global mathematical description of the problem. The analogy used was that of minimizing the electrical potential of a system charged particles within a charged domain. The particles in the presented analogy represent duals to mesh elements (i.e., quads or hexes). Particle movement is governed by a mathematical functional which accounts for inter-particles repulsive, attractive and alignment forces. This functional is minimized to find the optimal location and orientation of each particle. After the particles are connected a mesh can be easily resolved. The mathematical description for this problem is as easy to formulate in three-dimensions as it is in two- or one-dimensions. The meshing algorithm was developed within CoMeT. It can solve the two-dimensional meshing problem for convex and concave geometries in a purely automated fashion. Investigation of the robustness of the technique has shown a success rate of approximately 99% for the two-dimensional geometries tested. Run times to mesh a 100 element complex geometry were typically in the 10 minute range. Efficiency of the technique is still an issue that needs to be addressed. Performance is an issue that is critical for most engineers generating meshes. It was not for this project. The primary focus of this work was to investigate and evaluate a meshing algorithm/philosophy with efficiency issues being secondary. The algorithm was also extended to mesh three-dimensional geometries. Unfortunately, only simple geometries were tested before this project ended. The primary complexity in the extension was in the connectivity problem formulation. Defining all of the interparticle interactions that occur in three-dimensions and expressing them in mathematical relationships is very difficult.

  12. Masticatory biomechanics and its relevance to early hominid phylogeny: An examination of palatal thickness using finite-element analysis

    E-Print Network [OSTI]

    thickness using finite-element analysis David S. Strait a,*, Brian G. Richmond b , Mark A. Spencer c which early hominid features are functionally related to chewing. This paper uses finite-element Paranthropus synapomorphies putatively related to mastication. A finite-element model of 145,680 elements

  13. Method of lightening radiation darkened optical elements

    DOE Patents [OSTI]

    Reich, Frederich R. (Richland, WA); Schwankoff, Albert R. (W. Richland, WA)

    1980-01-01T23:59:59.000Z

    A method of lightening a radiation-darkened optical element in wich visible optical energy or electromagnetic radiation having a wavelength in the range of from about 2000 to about 20,000 angstroms is directed into the radiation-darkened optical element; the method may be used to lighten radiation-darkened optical element in-situ during the use of the optical element to transmit data by electronically separating the optical energy from the optical output by frequency filtering, data cooling, or interlacing the optic energy between data intervals.

  14. Low frequency eddy current finite element model validation and benchmark studies

    SciTech Connect (OSTI)

    Cherry, M.; Knopp, J. [Materials and Manufacturing Directorate, Air Force Research Laboratory (United States); Mooers, R.; Boehnlein, T. [University of Dayton Research Institute, Dayton, OH (United States); Aldrin, J. C. [Computational Tools, Gurnee, IL (United States); Sabbagh, H. A. [Victor Technologies, LLC, Bloomington, IN (United States)

    2011-06-23T23:59:59.000Z

    A finite element method (FEM) model was created to calculate the change in impedance of a coil due to the presence of a notch in a plate. The rectangular notches were created via electrical discharge machining (EDM) in a thick aluminum plate and were positioned at normal and oblique angles (10, 20, and 30 degrees) with respect to the vertical axis of the coil. The FEM method was chosen for this model due to its ability to solve problems in complicated geometries with the use of irregular mesh elements to discretize the solution domain. The change in impedance was calculated from the field variables in the simulation for each probe position along the parallel axis of the plate. The error between the model and the experimental data was approximately 5% for the majority of cases. The validated model was used to investigate more complex problems.

  15. A finite element approach to study cavitation instabilities in non-linear elastic solids under general loading conditions

    E-Print Network [OSTI]

    Lopez-Pamies, Oscar

    A finite element approach to study cavitation instabilities in non-linear elastic solids under a b s t r a c t This paper proposes an effective numerical method to study cavitation instabilities a single infinitesimal defect at its center. The occurrence of cavitation is identified as the event when

  16. Simulation of dynamic fracture using peridynamics, finite element modeling, and contact.

    SciTech Connect (OSTI)

    Littlewood, David John

    2010-11-01T23:59:59.000Z

    Peridynamics is a nonlocal extension of classical solid mechanics that allows for the modeling of bodies in which discontinuities occur spontaneously. Because the peridynamic expression for the balance of linear momentum does not contain spatial derivatives and is instead based on an integral equation, it is well suited for modeling phenomena involving spatial discontinuities such as crack formation and fracture. In this study, both peridynamics and classical finite element analysis are applied to simulate material response under dynamic blast loading conditions. A combined approach is utilized in which the portion of the simulation modeled with peridynamics interacts with the finite element portion of the model via a contact algorithm. The peridynamic portion of the analysis utilizes an elastic-plastic constitutive model with linear hardening. The peridynamic interface to the constitutive model is based on the calculation of an approximate deformation gradient, requiring the suppression of possible zero-energy modes. The classical finite element portion of the model utilizes a Johnson-Cook constitutive model. Simulation results are validated by direct comparison to expanding tube experiments. The coupled modeling approach successfully captures material response at the surface of the tube and the emerging fracture pattern. The coupling of peridynamics and finite element analysis via a contact algorithm has been shown to be a viable means for simulating material fracture in a high-velocity impact experiment. A combined peridynamics/finite element approach was applied to model an expanding tube experiment performed by Vogler, et al., in which loading on the tube is a result of Lexan slugs impacting inside the tube. The Lexan portion of the simulation was modeled with finite elements and a Johnson-Cook elastic-plastic material model in conjunction with an equation-of-state law. The steel tube portion of the simulation was modeled with peridynamics, an elastic-plastic material model, and a critical stretch bond damage model. The application of peridynamics to the tube portion of the model allowed the capture of the formation of cracks and eventual fragmentation of the tube. The simulation results yielded good agreement with the experimental results published by Vogler, et al., for the velocity and displacement profiles on the surface of the tube and the resulting fragment distribution. Numerical difficulties were encountered that required removal of hexahedron elements from the Lexan portion of the model over the course of the simulation. The significant number of inverted and nearly-inverted elements appearing over the course of the simulation is believed to be a result of irregularities in the contact between the Lexan and AerMet portions of the model, and was likely exacerbated by the ultra-high strength of the AerMet tube. Future simulations are planned in which the Lexan portion of the simulation is modeled with peridynamics, or with an alternative method such as smoothed particle hydrodynamics, with the goal of reducing these numerical difficulties.

  17. Spectral analysis method for detecting an element

    DOE Patents [OSTI]

    Blackwood, Larry G [Idaho Falls, ID; Edwards, Andrew J [Idaho Falls, ID; Jewell, James K [Idaho Falls, ID; Reber, Edward L [Idaho Falls, ID; Seabury, Edward H [Idaho Falls, ID

    2008-02-12T23:59:59.000Z

    A method for detecting an element is described and which includes the steps of providing a gamma-ray spectrum which has a region of interest which corresponds with a small amount of an element to be detected; providing nonparametric assumptions about a shape of the gamma-ray spectrum in the region of interest, and which would indicate the presence of the element to be detected; and applying a statistical test to the shape of the gamma-ray spectrum based upon the nonparametric assumptions to detect the small amount of the element to be detected.

  18. Version: 6/16/98 Keywords: wavy surface flow, finite element, longwave analysis, weakly-nonlinear analysis

    E-Print Network [OSTI]

    McCready, Mark J.

    Version: 6/16/98 Keywords: wavy surface flow, finite element, longwave analysis, weakly and drag are found, from finite element calculations, to increase as amplitude to approximately the third wavelength problem is solved numerically with a finite element formulation providing qualitative trends

  19. A Novel Inverse Finite Element Analysis to Assess Bone Fracture Healing in Mice Receiving Bone Marrow Mesenchymal Stem Cell Transplantation

    E-Print Network [OSTI]

    Miga, Michael I.

    A Novel Inverse Finite Element Analysis to Assess Bone Fracture Healing in Mice Receiving Bone generation, and an iterative optimization (using finite element analysis) of the fracture callus material approach includes acquisition of microCT image volumes, biomechanical testing, finite element mesh

  20. 48 2 Functional Analysis for Linear Saddle Point Problems 2.4 Examples of Pairs of Finite Element Spaces

    E-Print Network [OSTI]

    John, Volker

    48 2 Functional Analysis for Linear Saddle Point Problems 2.4 Examples of Pairs of Finite Element Spaces Violating the Discrete Inf-Sup Condition Remark 2.50. On simple pairs of finite element spaces. The simplest and most common finite element spaces are spaces of continuous functions which are piecewise

  1. Finite Element Analysis of Spherical Indentation Tests John Homans1, Sinisa Dj. Mesarovic2, Harish Radhakrishnan2

    E-Print Network [OSTI]

    Collins, Gary S.

    Finite Element Analysis of Spherical Indentation Tests John Homans1, Sinisa Dj. Mesarovic2, Harish finite element analysis program Abaqus 6.7. · Runs were repeated using a dense mesh to ensure accuracy not influence the solution. Objective · Develop a finite element model to analyze: i) Load-depth response

  2. GEO+, a finite element program on a personal computer

    E-Print Network [OSTI]

    Guillin, Clement Gerard

    1993-01-01T23:59:59.000Z

    and BW' is the difference between the largest and smallest node numbers for an element e. If it, i2, i3, and i4 are the node numbers of a 4- noded quadrilateral element e, the maximum element node number difference BWe is given by: BW'=Max( [ii...-is[, [i2-i3[, [i3-i4[, [i4-ii[, [i3-ii[, [i2-i4[) (84) One way to minimize the bandwidth is to number the nodes such that BWe is as small as possible in each element, and thereby keep the numbers of each element as close as possible. This property...

  3. A Mortar Mimetic Finite Difference Method on Non-Matching Grids

    E-Print Network [OSTI]

    Berndt, Markus

    of Petroleum and Geosystems Engineering, The University of Texas at Austin, Austin, TX 78712, U.S.A., mfw]. The methods has also been extended to locally refined meshes with hanging nodes [19], A connection between velocities in MFD methods on h2-uniform quadrilateral meshes was established. Mixed finite element (MFE

  4. Spectral element methods for elliptic problems in nonsmooth domains

    SciTech Connect (OSTI)

    Pathria, D.; Karniadakis, G.E. [Brown Univ., Providence, RI (United States)] [Brown Univ., Providence, RI (United States)

    1995-11-01T23:59:59.000Z

    Engineering applications frequently require the numerical solution of elliptic boundary value problems in irregularly shaped domains. For smooth problems spectral element methods have proved very successful, since they can accommodate fairly complicated geometries while retaining a rapid rate of convergence. Geometric singularities, however, often give rise to singular solutions. The accuracy of the spectral element methods is then degraded, and they offer no apparent advantage over low-order finite element methods. In many cases, however, the singular structure of the solution is known, and its form may be exploited by the spectral element method. Among the various ways of doing so (through supplementary basis functions, eigenfunction expansions, and graded meshes), the method of auxiliary mapping proves to be particularly effective. For certain simple cases, the problem is transformed to a coordinate system in which the solution is analytic, and exponential convergence is recovered. Even when this is not possible, the singularity is usually much weaker after mapping, so that other treatments are more effective in the new coordinate system. In this paper, the authors study different ways of treating singularities, and in particular, the method of auxiliary mapping coupled with the use of supplementary basis functions. Error estimates are presented explaining why the combined approach is more effective, and these estimates are confirmed through a number of numerical experiments for the Laplace. Poisson, and Helmholtz equations. 28 refs., 13 figs., 5 tabs.

  5. Performance of U/P finite elements subjected to nearly incompressible linear axisymmetric orthotropic conditions

    E-Print Network [OSTI]

    Bateman, Robert D

    2011-01-01T23:59:59.000Z

    Finite element analysis is a vastly expanding field which provides engineers a viable instrumentation to model and measure idealized constitutive stress strain relationships for various complex physical bodies. In the field ...

  6. TEA - a linear frequency domain finite element model for tidal embayment analysis

    E-Print Network [OSTI]

    Westerink, Joannes J.

    1984-01-01T23:59:59.000Z

    A frequency domain (harmonic) finite element model is developed for the numerical prediction of depth average circulation within small embayments. Such embayments are often characterized by irregular boundaries and bottom ...

  7. Ionic Polymer-Metal Composites: Thermodynamical Modeling and Finite Element Solution

    E-Print Network [OSTI]

    Arumugam, Jayavel

    2012-10-19T23:59:59.000Z

    cantilever beams involving pure bending which could be used in the characterization of the material parameters. A linear finite element solution along with a staggered time stepping algorithm is provided to numerically solve the governing equations...

  8. NUMERICAL ANALYSIS OF A FINITE ELEMENT SCHEME FOR THE APPROXIMATION OF HARMONIC MAPS INTO SURFACES

    E-Print Network [OSTI]

    Bartels, Soeren

    NUMERICAL ANALYSIS OF A FINITE ELEMENT SCHEME FOR THE APPROXIMATION OF HARMONIC MAPS INTO SURFACES Geometric partial differential equations and their analysis as well as numerical simulation have recently

  9. Design of a tricycle chassis using computer-aided design and finite element analysis

    E-Print Network [OSTI]

    Avila, Elliot

    2014-01-01T23:59:59.000Z

    Finite element analysis and computer-aided design are powerful tools for modeling complex systems and their responses to external stimuli. This paper explores how these techniques were employed in a highly iterative design ...

  10. On a tensor-based finite element model for the analysis of shell structures 

    E-Print Network [OSTI]

    Arciniega Aleman, Roman Augusto

    2006-04-12T23:59:59.000Z

    In the present study, we propose a computational model for the linear and nonlinear analysis of shell structures. We consider a tensor-based finite element formulation which describes the mathematical shell model in a ...

  11. A mixed-dimensional finite volume method for multiphase flow

    E-Print Network [OSTI]

    Cirpka, Olaf Arie

    ; Numerical reservoir simulation; Fractured reser- voir; Mixed-dimensional Finite Volume Method; MultigridA mixed-dimensional finite volume method for multiphase flow in fractured porous media Volker method for the fully coupled, fully implicit discretization of two-phase flow in fractured porous media

  12. A finite element viscous flow analysis in a radial turbine scroll

    E-Print Network [OSTI]

    Hill, Donald Lee

    1987-01-01T23:59:59.000Z

    A FINITE ELEMENT VISCOUS FLOW ANALYSIS IN A RADIAL TURBINE SCROLL A Thesis DONALD LEE HILL JR. Submitted to the Graduate College. of Texas ASM University in partial fulfillment of the requirement for the degree of MASTER OF SCIENCE December... 1987 Major Subject: Mechanical Engineering A FINITE ELEMENT VISCOUS FLOW ANALYSIS IN A RADIAL TURBINE SCROLL A Thesis by DONALD LEE HILL JR. Approved as to style snd content by: Dr. Erian A. Baskharone (Chairman of Conunittee) Dr. Alan B azzolo...

  13. Two dimensional finite element analysis of homogenization and failure in plain weave textile composites

    E-Print Network [OSTI]

    Kondagunta, Veeraraghava Gopal

    1993-01-01T23:59:59.000Z

    TWO DIMENSIONAL FINITE ELEMENT ANALYSIS OF HOMOGENIZATION AND FAILURE IN PLAIN WEAVE TEXTILE COMPOSITES A Thesis by VEERARAGHAVA GOPAL KONDAGUNTA Submitted to the Office of Graduate Studies of Texas ARM University In partial fulffllment... of the requirements for the degree of MASTER OF SCIENCE August 1993 Major Subject: Aerospace Engineering TWO DIMENSIONAL FINITE ELEMENT ANALYSIS OF HOMOGENIZATION AND FAILURE IN PLAIN WEAVE TEXTILE COMPOSITES A Thesis by VEERARAGHAVA GOPAL KONDAGUNTA...

  14. Evaluation of new techniques for two dimensional finite element analysis of woven composites

    E-Print Network [OSTI]

    Gundapaneni, Sitaram Chowdary

    1992-01-01T23:59:59.000Z

    EVALUATION OF NEW TECHNIQUES FOR TWO DIMENSIONAL FINITE ELEMENT ANALYSIS OF WOVEN COMPOSITES A Thesis by SITARAM CHOWDARY GUNDAPANENI Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment... of the requirement for the degree of MASTER OF SCIENCE DECEMBER 1992 Major Subject: Aerospace Engineering EVALUATION OF NEW TECHNIQUES FOR TWO DIMENSIONAL FINITE ELEMENT ANALYSIS OF WOVEN COMPOSITES A Thesis by SITARAM CHOWDARY GUNDAPANENI Approved...

  15. Multi-Dimensional Astrophysical Structural and Dynamical Analysis I. Development of a Nonlinear Finite Element Approach

    E-Print Network [OSTI]

    David L. Meier

    1998-11-10T23:59:59.000Z

    A new field of numerical astrophysics is introduced which addresses the solution of large, multidimensional structural or slowly-evolving problems (rotating stars, interacting binaries, thick advective accretion disks, four dimensional spacetimes, etc.). The technique employed is the Finite Element Method (FEM), commonly used to solve engineering structural problems. The approach developed herein has the following key features: 1. The computational mesh can extend into the time dimension, as well as space, perhaps only a few cells, or throughout spacetime. 2. Virtually all equations describing the astrophysics of continuous media, including the field equations, can be written in a compact form similar to that routinely solved by most engineering finite element codes. 3. The transformations that occur naturally in the four-dimensional FEM possess both coordinate and boost features, such that (a) although the computational mesh may have a complex, non-analytic, curvilinear structure, the physical equations still can be written in a simple coordinate system independent of the mesh geometry. (b) if the mesh has a complex flow velocity with respect to coordinate space, the transformations will form the proper arbitrary Lagrangian- Eulerian advective derivatives automatically. 4. The complex difference equations on the arbitrary curvilinear grid are generated automatically from encoded differential equations. This first paper concentrates on developing a robust and widely-applicable set of techniques using the nonlinear FEM and presents some examples.

  16. Finite element solution of multi-scale transport problems using the least squares based bubble function enrichment

    E-Print Network [OSTI]

    A. Yazdani; V. Nassehi

    2011-01-10T23:59:59.000Z

    This paper presents an optimum technique based on the least squares method for the derivation of the bubble functions to enrich the standard linear finite elements employed in the formulation of Galerkin weighted-residual statements. The element-level linear shape functions are enhanced with supplementary polynomial bubble functions with undetermined coefficients. The best least squares minimization of the residual functional obtained from the insertion of these trial functions into model equations results in an algebraic system of equations whose solution provides the unknown coefficients in terms of element-level nodal values. The normal finite element procedures for the construction of stiffness matrices may then be followed with no extra degree of freedom incurred as a result of such enrichment. The performance of the proposed method has been tested on a number of benchmark linear transport equations with the results compared against the exact and standard linear element solutions. It has been observed that low order bubble enriched elements produce more accurate approximations than the standard linear elements with no extra computational cost despite employing relatively crude mesh. However, for the solution of strongly convection or reaction dominated problems significantly higher order enrichments as well as extra mesh refinements will be required.

  17. A mixed-dimensional finite volume method for two-phase flow in fractured porous media

    E-Print Network [OSTI]

    Bastian, Peter

    . Key words: Multiphase flow; Numerical reservoir simulation; Fractured reservoir; MixedA mixed-dimensional finite volume method for two-phase flow in fractured porous media Volker-phase flow in fractured porous media. Fractures are dis- cretely modeled as lower dimensional elements

  18. Finite Element Simulation of Nonlinear Elastic Dynamics Using Cuda

    E-Print Network [OSTI]

    California at Berkeley, University of

    and e elements · A time-varying vector q(t) R3n which describes the displacement of the nodes positive definite linear system involving fext, fext, and M We focus on step 1 in this projectForce 8800 GT with 512 MB of memory. · Only uses single precision floating-point. ­ GPUs with double

  19. Finite element analysis of the Arquin-designed CMU wall under a dynamic (blast) load.

    SciTech Connect (OSTI)

    Lopez, Carlos; Petti, Jason P.

    2008-12-01T23:59:59.000Z

    The Arquin Corporation designed a CMU (concrete masonry unit) wall construction and reinforcement technique that includes steel wire and polymer spacers that is intended to facilitate a faster and stronger wall construction. Since the construction method for an Arquin-designed wall is different from current wall construction practices, finite element computer analyses were performed to estimate the ability of the wall to withstand a hypothetical dynamic load, similar to that of a blast from a nearby explosion. The response of the Arquin wall was compared to the response of an idealized standard masonry wall exposed to the same dynamic load. Results from the simulations show that the Arquin wall deformed less than the idealized standard wall under such loading conditions. As part of a different effort, Sandia National Laboratories also looked at the relative static response of the Arquin wall, results that are summarized in a separate SAND Report.

  20. Post-1906 stress recovery of the San Andreas fault system calculated from three-dimensional finite element analysis

    E-Print Network [OSTI]

    century as in the 19th century. In this paper I use finite element analysis to address the question, how element analysis Tom Parsons U.S. Geological Survey, Menlo Park, California, USA Received 23 August 2001 area seismic hazard assessment. This study presents a three-dimensional (3-D) finite element simulation

  1. COYOTE : a finite element computer program for nonlinear heat conduction problems. Part I, theoretical background.

    SciTech Connect (OSTI)

    Glass, Micheal W.; Hogan, Roy E., Jr.; Gartling, David K.

    2010-03-01T23:59:59.000Z

    The need for the engineering analysis of systems in which the transport of thermal energy occurs primarily through a conduction process is a common situation. For all but the simplest geometries and boundary conditions, analytic solutions to heat conduction problems are unavailable, thus forcing the analyst to call upon some type of approximate numerical procedure. A wide variety of numerical packages currently exist for such applications, ranging in sophistication from the large, general purpose, commercial codes, such as COMSOL, COSMOSWorks, ABAQUS and TSS to codes written by individuals for specific problem applications. The original purpose for developing the finite element code described here, COYOTE, was to bridge the gap between the complex commercial codes and the more simplistic, individual application programs. COYOTE was designed to treat most of the standard conduction problems of interest with a user-oriented input structure and format that was easily learned and remembered. Because of its architecture, the code has also proved useful for research in numerical algorithms and development of thermal analysis capabilities. This general philosophy has been retained in the current version of the program, COYOTE, Version 5.0, though the capabilities of the code have been significantly expanded. A major change in the code is its availability on parallel computer architectures and the increase in problem complexity and size that this implies. The present document describes the theoretical and numerical background for the COYOTE program. This volume is intended as a background document for the user's manual. Potential users of COYOTE are encouraged to become familiar with the present report and the simple example analyses reported in before using the program. The theoretical and numerical background for the finite element computer program, COYOTE, is presented in detail. COYOTE is designed for the multi-dimensional analysis of nonlinear heat conduction problems. A general description of the boundary value problems treated by the program is presented. The finite element formulation and the associated numerical methods used in COYOTE are also outlined. Instructions for use of the code are documented in SAND2010-0714.

  2. Finite element analysis of the distortion of a crystal monochromator from synchrotron radiation thermal loading

    SciTech Connect (OSTI)

    Edwards, W.R.; Hoyer, E.H.; Thompson, A.C.

    1985-10-01T23:59:59.000Z

    The first crystal of the Brown-Hower x-ray monochromator of the LBL-EXXON 54 pole wiggler beamline at Stanford Synchrotron Radiation Laboratory (SSRL) is subjected to intense synchrotron radiation. To provide an accurate thermal/structural analysis of the existing monochromator design, a finite element analysis (FEA) was performed. A very high and extremely localized heat flux is incident on the Si (220) crystal. The crystal, which possesses pronouncedly temperature-dependent orthotropic properties, in combination with the localized heat load, make the analysis ideally suited for finite element techniques. Characterization of the incident synchrotron radiation is discussed, followed by a review of the techniques employed in modeling the monochromator and its thermal/structural boundary conditions. The results of the finite element analysis, three-dimensional temperature distributions, surface displacements and slopes, and stresses, in the area of interest, are presented. Lastly, the effects these results have on monochromator output flux and resolution are examined.

  3. Determining matrix elements and resonance widths from finite volume: the dangerous mu-terms

    E-Print Network [OSTI]

    G. Takacs

    2011-10-10T23:59:59.000Z

    The standard numerical approach to determining matrix elements of local operators and width of resonances uses the finite volume dependence of energy levels and matrix elements. Finite size corrections that decay exponentially in the volume are usually neglected or taken into account using perturbation expansion in effective field theory. Using two-dimensional sine-Gordon field theory as "toy model" it is shown that some exponential finite size effects could be much larger than previously thought, potentially spoiling the determination of matrix elements in frameworks such as lattice QCD. The particular class of finite size corrections considered here are mu-terms arising from bound state poles in the scattering amplitudes. In sine-Gordon model, these can be explicitly evaluated and shown to explain the observed discrepancies to high precision. It is argued that the effects observed are not special to the two-dimensional setting, but rather depend on general field theoretic features that are common with models relevant for particle physics. It is important to understand these finite size corrections as they present a potentially dangerous source of systematic errors for the determination of matrix elements and resonance widths.

  4. Visualization of high resolution, three-dimensional, nonlinear finite element analyses

    SciTech Connect (OSTI)

    Christon, M.A.; Spelce, T.E.

    1992-03-24T23:59:59.000Z

    A two-pass surface extraction algorithm for adaptive finite element meshes is presented in the context of a visualization study for a particle impact and a turbine containment problem. The direct use of finite element data structures for the computation of external surfaces, surface normals, and derived physical quantities is discussed. An overview of the in-betweening algorithm which compensates for rigid body dynamics is presented with a brief discussion of a {open_quotes}direct-to-videodisk{close_quotes} animation strategy.

  5. Training course No. 2: the implementation of FEMWASTE (ORNL-5601) computer program. Final report. [Finite Element Model of WASTE transport

    SciTech Connect (OSTI)

    Yeh, G.T.

    1982-10-01T23:59:59.000Z

    This report documents a training course conducted for the US Nuclear Regulatory Commission (NRC) on the implementation of a Finite Element Model of WASTE transport through saturated-unsaturated porous media (FEMWASTE) - ORNL-5601. In addition to presenting basic program operations (Appendices A-V through A-VII), the course also covers the following topics: (1) Heuristic derivation of governing equations based on physical and chemical principles, (2) finite element derivation of FEMWASTE, (3) various numerical schemes provided by FEMWASTE, (4) FEMWASTE program structure, and (5) running of three samples problems to demonstrate various options the FEMWASTE can handle. The purpose of the training seminar is to enable NRC staff to use the model (and to be able to modify the code, if necessary) for checking information provided by a licensee, for evaluating alternative sites and designs for burial, and for comparing their results from other methods of solution.

  6. Finite element analysis and computed tomography based structural rigidity analysis of rat tibia with simulated lytic defects

    E-Print Network [OSTI]

    Vaziri, Ashkan

    Finite element analysis and computed tomography based structural rigidity analysis of rat tibia) (Damron et al., 2003; Mirels, 1989). In contrast, Computed Tomography based Structural Rigidly Ana- lysis

  7. Iterative solutions to large sparse finite element equations

    E-Print Network [OSTI]

    Wang, Hongbing

    1995-01-01T23:59:59.000Z

    direct solvers. The theory behind the convergence rate relationship and storage requirements for the preconditioned conjugate gradient methods using the diagonal scaling, incomplete Cholesky decomposition and SSOR preconditioners is explained in detail...

  8. NONCONFORMING FINITE ELEMENTS FOR REISSNER-MINDLIN PLATES

    E-Print Network [OSTI]

    Lovadina, Carlo

    to the more complex (and more inter- esting) case of shell problems. Indeed, the methods we will consider take.e. rotations and deflections) share the same nodes and degrees of freedom. The paper is organized as follows the transversal displacement, the scaled shear stresses and g a given transversal load. Finally, is the usual

  9. Finite-element modeling of subglacial cavities and related friction law

    E-Print Network [OSTI]

    Gagliardini, Olivier

    Finite-element modeling of subglacial cavities and related friction law O. Gagliardini,1 D. Cohen,2, the friction law, can be easily extended from linear to nonlinear ice rheology and is bounded even for bedrocks with locally infinite slopes. Combining our results with earlier works by others, a phenomenological friction

  10. FINITE ELEMENT MODEL-BASED STRUCTURAL HEALTH MONITORING (SHM) SYSTEMS FOR COMPOSITE MATERIAL UNDER

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    FINITE ELEMENT MODEL-BASED STRUCTURAL HEALTH MONITORING (SHM) SYSTEMS FOR COMPOSITE MATERIAL UNDER). To design a Structural Health Monitoring (SHM) system, it is important to understand phenomenologically Workshop on Structural Health Monitoring July 8-11, 2014. La Cité, Nantes, France Copyright © Inria (2014

  11. Thermo-elasto-plastic finite element analysis of quasi-state processes in Eulerian reference

    E-Print Network [OSTI]

    Michaleris, Panagiotis

    Thermo-elasto-plastic finite element analysis of quasi-state processes in Eulerian reference frames ­ Incremental scheme ­ Fine mesh along entire heat source path ­ Lengthy computer runs · Elasto-Plasticity at times ti and ti-1, respectively. Penn State University 5 #12;Elasto-Plasticity Equilibrium: r(r, t) + b

  12. Local Flaps: A Real-Time Finite Element Based Solution to the Plastic Surgery

    E-Print Network [OSTI]

    Liblit, Ben

    Local Flaps: A Real-Time Finite Element Based Solution to the Plastic Surgery Defect Puzzle fundamental challenges in plastic surgery is the alter- ation of the geometry and topology of the skin for the patient after the procedure is completed. The plastic surgeon must look at the defect created

  13. A Finite Element Algorithm of a Nonlinear Diffusive Climate Energy Balance Model

    E-Print Network [OSTI]

    Díaz, Jesús Ildefonso

    A Finite Element Algorithm of a Nonlinear Diffusive Climate Energy Balance Model R. BERMEJO,1 J. This model belongs to the category of energy balance models introduced independently by the climatologists M climate. The energy balance model we are dealing with consists of a two-dimensional nonlinear parabolic

  14. A mixed finite-element scheme of a semiconductor energy-transport model

    E-Print Network [OSTI]

    Hanke-Bourgeois, Martin

    A mixed finite-element scheme of a semiconductor energy-transport model using dual entropy which are able to deal with physical effects such as carrier heating and velocity overshoot. The energy that the Joule heating term vanishes if the dual entropy variables w1 = (µ - V )/T and w2 = -1/T are employed

  15. Finite element implementation of an eigenfunction solution for the contact pressure variation due to wear

    E-Print Network [OSTI]

    Barber, James R.

    Elsevier B.V. All rights reserved. 1. Introduction The contact pressure between two sliding bodiesFinite element implementation of an eigenfunction solution for the contact pressure variation due is influenced by the profile of the bodies, but this is altered by wear so that the pressure distribution

  16. Modeling of board-level package by Finite Element Analysis and laser interferometer measurements

    E-Print Network [OSTI]

    Cao, Wenwu

    and Vibration, School of Mechanical Engineering, Shanghai Jiao Tong University, 800 Dong Chuan Road, Shanghai be simulated by an equivalent Finite Element Analysis (FEA) model of the board-level package, in which detailed layer structure of the Print Circuit Board (PCB), signal wires and through-holes were ignored

  17. A Finite-Element Model for Simulation of Carbon Dioxide Sequestration

    SciTech Connect (OSTI)

    Bao, Jie; Xu, Zhijie; Fang, Yilin

    2014-09-01T23:59:59.000Z

    Herein, we present a coupled thermal-hydro-mechanical model for geological sequestration of carbon dioxide followed by the stress, deformation, and shear-slip failure analysis. This fully coupled model considers the geomechanical response, fluid flow, and thermal transport relevant to geological sequestration. Both analytical solutions and numerical approach via finite element model are introduced for solving the thermal-hydro-mechanical model. Analytical solutions for pressure, temperature, deformation, and stress field were obtained for a simplified typical geological sequestration scenario. The finite element model is more general and can be used for arbitrary geometry. It was built on an open-source finite element code, Elmer, and was designed to simulate the entire period of CO2 injection (up to decades) both stably and accurately—even for large time steps. The shear-slip failure analysis was implemented based on the numerical results from the finite element model. The analysis reveals the potential failure zone caused by the fluid injection and thermal effect. From the simulation results, the thermal effect is shown to enhance well injectivity, especially at the early time of the injection. However, it also causes some side effects, such as the appearance of a small failure zone in the caprock. The coupled thermal-hydro-mechanical model improves prediction of displacement, stress distribution, and potential failure zone compared to the model that neglects non-isothermal effects, especially in an area with high geothermal gradient.

  18. Finite Element Modeling of Rate-Dependent Ratcheting in Granular Materials

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Finite Element Modeling of Rate-Dependent Ratcheting in Granular Materials A. Karrech a,, A. Seibi thrust is the introduc- tion of the "Chicago" law in a continuum approach to account for the ratcheting or jointly with the ratcheting effect. Non-associated plasticity is adopted since the paper mainly targets

  19. Mortar finite element discretization of a model coupling Darcy and Stokes equations

    E-Print Network [OSTI]

    Boyer, Edmond

    Mortar finite element discretization of a model coupling Darcy and Stokes equations by C. Bernardi1 flowing on a homogeneous porous ground, we consider a system where the Darcy and Stokes equations´erons un syst`eme o`u les ´equations de Darcy et de Stokes sont coupl´ees par des conditions de raccord

  20. Mortar finite element discretization of the time dependent nonlinear Darcy's equations

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Mortar finite element discretization of the time dependent nonlinear Darcy's equations by Karima variations have been handled in [1] and [6], while time-dependent Darcy's equations have been studied in [5 equations model the flow in a porous medium, but underground porous media are most often nonhomogeneous

  1. A POSTERIORI ERROR ESTIMATE FOR THE H(div) CONFORMING MIXED FINITE ELEMENT FOR THE COUPLED

    E-Print Network [OSTI]

    Wang, Yanqiu

    for the Darcy equation, and vice versa. Special techniques usually need to be employed. In [3], a conformingA POSTERIORI ERROR ESTIMATE FOR THE H(div) CONFORMING MIXED FINITE ELEMENT FOR THE COUPLED DARCY pro- posed for the coupled Darcy-Stokes flow in [30], which imposes normal con- tinuity

  2. Hybrid spectral/finite element analysis of dynamic delamination of patterned thin films

    E-Print Network [OSTI]

    Sottos, Nancy R.

    Hybrid spectral/finite element analysis of dynamic delamination of patterned thin films Phuong Tran Accepted 10 March 2008 Available online 20 March 2008 Keywords: Thin film Adhesion Delamination Dynamic analysis is performed to investigate the dynamic edge delamination of patterned thin films from a substrate

  3. FAST AND ACCURATE FINITE ELEMENT APPROXIMATION OF WAVE MAPS INTO SPHERES

    E-Print Network [OSTI]

    . The scheme has a second order consistency property and is energy preserving and reversible. Its unconditional show in this article that a similar strategy can be applied for finite element discretizations this first-order system in time according to dtuk+1 = uk+1/2 × wk+1/2 , dtwk+1 = uk+1/2 × uk+1/2 , where dt

  4. Salt tectonics driven by differential sediment loading: Stability analysis and finite element experiments

    E-Print Network [OSTI]

    Beaumont, Christopher

    1 Salt tectonics driven by differential sediment loading: Stability analysis and finite element University of Leeds Leeds LS2 9JT United Kingdom e-mail: lykke@earth.leeds.ac.uk Short running title: Salt salt layer drives salt deformation and has a significant impact on the structural evolution

  5. Finite-element discretization of a linearized 2 -D model for lubricated oil transportation

    E-Print Network [OSTI]

    Frey, Pascal

    Finite-element discretization of a linearized 2 - D model for lubricated oil transportation V acts as a lubricant by coating the wall of the pipeline, thus preventing the oil from adhering is devoted to the numerical simulation of a linearized model for the lubricated trans- portation of heavy

  6. TIME DOMAIN ELECTROMAGNETIC SCATTERING USING FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS

    E-Print Network [OSTI]

    TIME DOMAIN ELECTROMAGNETIC SCATTERING USING FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS H. T medium by a non­harmonic microwave pulse from an antenna source in free space, and we compute of the electromagnetic interrogation problems introduced in [3]. Here we use perfectly matched layers (PMLs) as absorbing

  7. TIME DOMAIN ELECTROMAGNETIC SCATTERING USING FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS

    E-Print Network [OSTI]

    TIME DOMAIN ELECTROMAGNETIC SCATTERING USING FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS H. T medium by a non-harmonic microwave pulse from an antenna source in free space, and we compute-dimensional version of the electromagnetic interrogation problems introduced in [3]. Here we use perfectly matched

  8. A mixed finite element discretization scheme for a concrete carbonation model with concentration-dependent porosity

    E-Print Network [OSTI]

    Suciu, Nicolae

    A mixed finite element discretization scheme for a concrete carbonation model with concentration/unsaturated porous media. The special features of our problem are twofold: the reaction produces water and therefore for the case of concrete carbonation ­ one of the most important physico-chemical processes affecting

  9. Simulation of dynamics-coupling in piezoelectric tube scanners by reduced order finite element analysis

    E-Print Network [OSTI]

    Fleming, Andrew J.

    and imperfections during the manufactur- ing process of the tube on the dynamics-coupling-caused errors in openSimulation of dynamics-coupling in piezoelectric tube scanners by reduced order finite element as the compensation of dynamics-coupling effects. The present article gives a detailed description of the fully

  10. Geological formation - drill string dynamic interaction finite-element program (GEODYN). Phase 1. Theoretical description

    SciTech Connect (OSTI)

    Baird, J.A.; Apostal, M.C.; Rotelli, R.L. Jr.; Tinianow, M.A.; Wormley, D.N.

    1984-06-01T23:59:59.000Z

    The Theoretical Description for the GEODYN interactive finite-element computer program is presented. The program is capable of performing the analysis of the three-dimensional transient dynamic response of a Polycrystalline Diamond Compact Bit-Bit Sub arising from the intermittent contact of the bit with the downhole rock formations. The program accommodates nonlinear, time-dependent, loading and boundary conditions.

  11. Finite element analysis of a current density -electric eld formulation of Bean's model for superconductivity

    E-Print Network [OSTI]

    Styles, Vanessa

    Finite element analysis of a current density - electric #12;eld formulation of Bean's model a current density-electric #12;eld formulation of Bean's model for the experimental set-up of a in#12;nitely. As observed by Bossavit, [4], Bean's critical state model can be formulated as a degenerate Stefan problem

  12. A FINITE ELEMENT APPROXIMATION OF A VARIATIONAL INEQUALITY FORMULATION OF BEAN'S MODEL FOR

    E-Print Network [OSTI]

    Elliott, Charles

    A FINITE ELEMENT APPROXIMATION OF A VARIATIONAL INEQUALITY FORMULATION OF BEAN'S MODEL of a variational formulation of Bean's model for the physical configuration of an infinitely long cylindrical's axis, for x R2 . The well-known Bean critical state model can be formulated as an evolutionary

  13. Simulation of Lung Behaviour with Finite Elements : Inuence of Bio-Mechanical Parameters

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    Simulation of Lung Behaviour with Finite Elements : Inuence of Bio-Mechanical Parameters Pierre by medical needs, we propose to simu- late lung deformation and motion during respiration to track tumours. This paper presents a model of lung behaviour based on a continuous media mechan- ics model and solved

  14. A UNSTRUCTURED NODAL SPECTRAL-ELEMENT METHOD FOR ...

    E-Print Network [OSTI]

    2012-01-17T23:59:59.000Z

    [19] S. J. Sherwin and G. E. Karniadakis. A triangular spectral element method; applications to the incompressible Navier-Stokes equations. Comput. Methods ...

  15. Percutaneous Mitral Valve Dilatation: Single Balloon versus Double Balloon - A Finite Element Study

    SciTech Connect (OSTI)

    Schievano, Silvia; Kunzelman, Karyn; Nicosia, Mark; Cochran, R. P.; Einstein, Daniel R.; Khambadkone, Sachin; Bonheoffer, Philipp

    2009-01-01T23:59:59.000Z

    Background: Percutaneous mitral valve (MV) dilatation is performed with either a single balloon (SB) or double balloon (DB) technique. The aim of this study was to compare the two balloon system results using the finite element (FE) method. Methods and Results: An established FE model of the MV was modified by fusing the MV leaflet edges at commissure level to simulate a stenotic valve (orifice area=180mm2). A FE model of a 30mm SB (low-pressure, elastomeric balloon) and an 18mm DB system (high-pressure, non-elastic balloon) was created. Both SB and DB simulations resulted in splitting of the commissures and subsequent stenosis dilatation (final MV area=610mm2 and 560mm2 respectively). Stresses induced by the two balloon systems varied across the valve. At the end of inflation, SB showed higher stresses in the central part of the leaflets and at the commissures compared to DB simulation, which demonstrated a more uniform stress distribution. The higher stresses in the SB analysis were due to the mismatch of the round balloon shape with the oval mitral orifice. The commissural split was not easily accomplished with the SB due to its high compliance. The high pressure applied to the DB guaranteed the commissural split even when high forces were required to break the commissure welds. Conclusions: The FE model demonstrated that MV dilatation can be accomplished by both SB and DB techniques. However, the DB method resulted in higher probability of splitting of the fused commissures and less damage caused to the MV leaflets by overstretching.

  16. Finite Element Model of Training in the superconducting quadrupole magnet SQ02

    SciTech Connect (OSTI)

    Caspi, Shlomo; Ferracin, Paolo

    2007-11-01T23:59:59.000Z

    This paper describes the use of 3D finite element models to study training in superconducting magnets. The simulations are used to examine coil displacements when the electromagnetic forces are cycled, and compute the frictional energy released during conductor motion with the resulting temperature rise. A computed training curve is then presented and discussed. The results from the numerical computations are compared with test results of the Nb{sub 3}Sn racetrack quadrupole magnet SQ02.

  17. R.H. Saunders GS concrete growth mitigation project instrumentation and finite element analysis

    SciTech Connect (OSTI)

    Adeghe, L.; Hindy, A.; Ho, M.S. [Hydroelectric Business Unit, Ontario Hydro, Toronto, Ontario (Canada)

    1995-12-31T23:59:59.000Z

    In 1991, concrete expansion due to Alkali-Silica Reaction (ASR), was diagnosed as the cause of operational and structural problems at R.H. Saunders Generating Station. Reestablishment of contraction joints between the generating units by cutting slots was initiated in January 1993 to relieve accumulated stresses and allow further expansion without deforming the turbine-generators. An integral part of the remedial measure is the installation of an extensive instrumentation system and development of detailed finite element models of the powerhouse structures. In total, two hundred and thirty five instruments have been installed. The instruments include surface extensometers, borehole extensometers, crackmeters, stressmeters, pendulums, thermocouples, borehole convergence meters and strain gauges. All the instruments are monitored continuously by an Automatic Data Acquisition System (ADAS) which allows instrument data to be evaluated remotely. This instrumentation system is being used to collect data on the structural concrete expansion and response to slot cutting. To complement the instrumentation data, three different finite element models have been developed for use. The models range from a very detailed representation of a single powerhouse unit to a less detailed model of the sixteen-unit powerhouse. The finite element models have been calibrated to reflect measured data and subsequently used to estimate the location and frequency of future cuts.

  18. Diffusion Synthetic Acceleration for High-Order Discontinuous Finite Element SN Transport Schemes and Application to Locally Refined Unstructured Meshes

    SciTech Connect (OSTI)

    Yaqi Wang; Jean C. Ragusa

    2011-10-01T23:59:59.000Z

    Diffusion synthetic acceleration (DSA) schemes compatible with adaptive mesh refinement (AMR) grids are derived for the SN transport equations discretized using high-order discontinuous finite elements. These schemes are directly obtained from the discretized transport equations by assuming a linear dependence in angle of the angular flux along with an exact Fick's law and, therefore, are categorized as partially consistent. These schemes are akin to the symmetric interior penalty technique applied to elliptic problems and are all based on a second-order discontinuous finite element discretization of a diffusion equation (as opposed to a mixed or P1 formulation). Therefore, they only have the scalar flux as unknowns. A Fourier analysis has been carried out to determine the convergence properties of the three proposed DSA schemes for various cell optical thicknesses and aspect ratios. Out of the three DSA schemes derived, the modified interior penalty (MIP) scheme is stable and effective for realistic problems, even with distorted elements, but loses effectiveness for some highly heterogeneous configurations. The MIP scheme is also symmetric positive definite and can be solved efficiently with a preconditioned conjugate gradient method. Its implementation in an AMR SN transport code has been performed for both source iteration and GMRes-based transport solves, with polynomial orders up to 4. Numerical results are provided and show good agreement with the Fourier analysis results. Results on AMR grids demonstrate that the cost of DSA can be kept low on locally refined meshes.

  19. Modeling brittle fracture, slip weakening, and variable friction in geomaterials with an embedded strong discontinuity finite element.

    SciTech Connect (OSTI)

    Regueiro, Richard A. (University of Colorado, Boulder, CO); Borja, R. I. (Stanford University, Stanford, CA); Foster, C. D. (Stanford University, Stanford, CA)

    2006-10-01T23:59:59.000Z

    Localized shear deformation plays an important role in a number of geotechnical and geological processes. Slope failures, the formation and propagation of faults, cracking in concrete dams, and shear fractures in subsiding hydrocarbon reservoirs are examples of important effects of shear localization. Traditional engineering analyses of these phenomena, such as limit equilibrium techniques, make certain assumptions on the shape of the failure surface as well as other simplifications. While these methods may be adequate for the applications for which they were designed, it is difficult to extrapolate the results to more general scenarios. An alternative approach is to use a numerical modeling technique, such as the finite element method, to predict localization. While standard finite elements can model a wide variety of loading situations and geometries quite well, for numerical reasons they have difficulty capturing the softening and anisotropic damage that accompanies localization. By introducing an enhancement to the element in the form of a fracture surface at an arbitrary position and orientation in the element, we can regularize the solution, model the weakening response, and track the relative motion of the surfaces. To properly model the slip along these surfaces, the traction-displacement response must be properly captured. This report focuses on the development of a constitutive model appropriate to localizing geomaterials, and the embedding of this model into the enhanced finite element framework. This modeling covers two distinct phases. The first, usually brief, phase is the weakening response as the material transitions from intact continuum to a body with a cohesionless fractured surface. Once the cohesion has been eliminated, the response along the surface is completely frictional. We have focused on a rate- and state-dependent frictional model that captures stable and unstable slip along the surface. This model is embedded numerically into the element using a generalized trapezoidal formulation. While the focus is on the constitutive model of interest, the framework is also developed for a general surface response. This report summarizes the major research and development accomplishments for the LDRD project titled 'Cohesive Zone Modeling of Failure in Geomaterials: Formulation and Implementation of a Strong Discontinuity Model Incorporating the Effect of Slip Speed on Frictional Resistance'. This project supported a strategic partnership between Sandia National Laboratories and Stanford University by providing funding for the lead author, Craig Foster, during his doctoral research.

  20. PARALLEL COMPUTING OF UNSATURATED SOILS USING ELEMENT-BY-ELEMENT AND DOMAIN DEOMPOSITON METHODS

    E-Print Network [OSTI]

    Augarde, Charles

    the processors with minimum data exchanges and load balancing. The linear system of equations was solved using and primary unknowns. The mass balance of water flow leads to the following equation: ( ) 0=+ ww wr div t n as in compacted soils used as fill material. The accurate analysis of unsaturated soils by the finite element

  1. A fabrication method for integrated filter elements with inductance cancellation

    E-Print Network [OSTI]

    Perreault, David J.

    This paper outlines a fabrication method for integrated filter elements. An integrated filter element is a three- (or more) terminal device comprising a capacitor and coupled air-core magnetic windings, in which the magnetic ...

  2. Operator-adapted finite element wavelets : theory and applications to a posteriori error estimation and adaptive computational modeling

    E-Print Network [OSTI]

    Sudarshan, Raghunathan, 1978-

    2005-01-01T23:59:59.000Z

    We propose a simple and unified approach for a posteriori error estimation and adaptive mesh refinement in finite element analysis using multiresolution signal processing principles. Given a sequence of nested discretizations ...

  3. Thermo-elasto-plastic finite element analysis of quasi-state processes in Eulerian reference J. Y. Shanghvi

    E-Print Network [OSTI]

    Michaleris, Panagiotis

    Thermo-elasto-plastic finite element analysis of quasi-state processes in Eulerian reference frames-state one way coupled thermo-elasto- plastic systems has been developed. The formulation is suitable

  4. A voxel-based finite element model for the prediction of bladder deformation

    SciTech Connect (OSTI)

    Chai Xiangfei; Herk, Marcel van; Hulshof, Maarten C. C. M.; Bel, Arjan [Radiation Oncology Department, Academic Medical Center, University of Amsterdam, 1105 AZ Amsterdam (Netherlands); Radiation Oncology Department, Netherlands Cancer Institute, 1066 CX Amsterdam (Netherlands); Radiation Oncology Department, Academic Medical Center, University of Amsterdam, 1105 AZ Amsterdam (Netherlands)

    2012-01-15T23:59:59.000Z

    Purpose: A finite element (FE) bladder model was previously developed to predict bladder deformation caused by bladder filling change. However, two factors prevent a wide application of FE models: (1) the labor required to construct a FE model with high quality mesh and (2) long computation time needed to construct the FE model and solve the FE equations. In this work, we address these issues by constructing a low-resolution voxel-based FE bladder model directly from the binary segmentation images and compare the accuracy and computational efficiency of the voxel-based model used to simulate bladder deformation with those of a classical FE model with a tetrahedral mesh. Methods: For ten healthy volunteers, a series of MRI scans of the pelvic region was recorded at regular intervals of 10 min over 1 h. For this series of scans, the bladder volume gradually increased while rectal volume remained constant. All pelvic structures were defined from a reference image for each volunteer, including bladder wall, small bowel, prostate (male), uterus (female), rectum, pelvic bone, spine, and the rest of the body. Four separate FE models were constructed from these structures: one with a tetrahedral mesh (used in previous study), one with a uniform hexahedral mesh, one with a nonuniform hexahedral mesh, and one with a low-resolution nonuniform hexahedral mesh. Appropriate material properties were assigned to all structures and uniform pressure was applied to the inner bladder wall to simulate bladder deformation from urine inflow. Performance of the hexahedral meshes was evaluated against the performance of the standard tetrahedral mesh by comparing the accuracy of bladder shape prediction and computational efficiency. Results: FE model with a hexahedral mesh can be quickly and automatically constructed. No substantial differences were observed between the simulation results of the tetrahedral mesh and hexahedral meshes (<1% difference in mean dice similarity coefficient to manual contours and <0.02 cm difference in mean standard deviation of residual errors). The average equation solving time (without manual intervention) for the first two types of hexahedral meshes increased to 2.3 h and 2.6 h compared to the 1.1 h needed for the tetrahedral mesh, however, the low-resolution nonuniform hexahedral mesh dramatically decreased the equation solving time to 3 min without reducing accuracy. Conclusions: Voxel-based mesh generation allows fast, automatic, and robust creation of finite element bladder models directly from binary segmentation images without user intervention. Even the low-resolution voxel-based hexahedral mesh yields comparable accuracy in bladder shape prediction and more than 20 times faster in computational speed compared to the tetrahedral mesh. This approach makes it more feasible and accessible to apply FE method to model bladder deformation in adaptive radiotherapy.

  5. Adaptive Finite Element Modeling Techniques for the Poisson-Boltzmann Equation

    E-Print Network [OSTI]

    Michael Holst; James Andrew McCammon; Zeyun Yu; Yongcheng Zhou; Yunrong Zhu

    2010-09-30T23:59:59.000Z

    We develop an efficient and reliable adaptive finite element method (AFEM) for the nonlinear Poisson-Boltzmann equation (PBE). We first examine the regularization technique of Chen, Holst, and Xu; this technique made possible the first a priori pointwise estimates and the first complete solution and approximation theory for the Poisson-Boltzmann equation. It also made possible the first provably convergent discretization of the PBE, and allowed for the development of a provably convergent AFEM for the PBE. However, in practice the regularization turns out to be numerically ill-conditioned. In this article, we examine a second regularization, and establish a number of basic results to ensure that the new approach produces the same mathematical advantages of the original regularization, without the ill-conditioning property. We then design an AFEM scheme based on the new regularized problem, and show that the resulting AFEM scheme is accurate and reliable, by proving a contraction result for the error. This result, which is one of the first results of this type for nonlinear elliptic problems, is based on using continuous and discrete a priori pointwise estimates to establish quasi-orthogonality. To provide a high-quality geometric model as input to the AFEM algorithm, we also describe a class of feature-preserving adaptive mesh generation algorithms designed specifically for constructing meshes of biomolecular structures, based on the intrinsic local structure tensor of the molecular surface. The stability advantages of the new regularization are demonstrated using an FETK-based implementation, through comparisons with the original regularization approach for a model problem. The convergence and accuracy of the overall AFEM algorithm is also illustrated by numerical approximation of electrostatic solvation energy for an insulin protein.

  6. Static and Fatigue Analysis of Wind Turbine Blades Subject to Cold Weather Conditions Using Finite Element Analysis

    E-Print Network [OSTI]

    Victoria, University of

    Static and Fatigue Analysis of Wind Turbine Blades Subject to Cold Weather Conditions Using Finite Turbine Blades Subject to Cold Weather Conditions Using Finite Element Analysis by Patricio Andres Lillo experienced in candi- date Canadian wind turbine deployment locations. The thesis then narrows its focus

  7. Parallel 3D Finite Element Numerical Modelling of DC Electron Guns

    SciTech Connect (OSTI)

    Prudencio, E.; Candel, A.; Ge, L.; Kabel, A.; Ko, K.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; /SLAC

    2008-02-04T23:59:59.000Z

    In this paper we present Gun3P, a parallel 3D finite element application that the Advanced Computations Department at the Stanford Linear Accelerator Center is developing for the analysis of beam formation in DC guns and beam transport in klystrons. Gun3P is targeted specially to complex geometries that cannot be described by 2D models and cannot be easily handled by finite difference discretizations. Its parallel capability allows simulations with more accuracy and less processing time than packages currently available. We present simulation results for the L-band Sheet Beam Klystron DC gun, in which case Gun3P is able to reduce simulation time from days to some hours.

  8. Curved finite elements by the method of initial strains

    E-Print Network [OSTI]

    Leick, Roger Dale

    1974-01-01T23:59:59.000Z

    Bw) yz By Bz P (By Bz By Bz + By Bz/ Bx ~ Bw I/(Bu Bu Bv Bv Bw Bw) 'zx Bz+ Bx+P QBz Bx+ Bz Bx+ Bz Bx/ (3 I lj) For tais development e , c , and s are zero from the discrete zP yzl zx Kirchhoff hypothesis and hence need not be considered further... this research. Thanks is also expressed to the other two committee members who contributed to this thesis, Dr. W. E. Haisler and Dr. W. L. Perry. This research was supported by Sandia Corporation contract No. 829005. DEDICATION This thesis is dedicated...

  9. stabilized finite element methods for coupled geomechanics and ...

    E-Print Network [OSTI]

    Abstract v. Acknowledgments vii. 1 Introduction. 1. 1.1 Problem Statement . ..... recent studies have shifted focus onto coupled modeling (e.g., Settari & Price 1984, ... Water, oil and natural gas are all fluids that exist within the pore spaces and ...

  10. The Finite Volume Element Method for Diffusion ... - Purdue University

    E-Print Network [OSTI]

    Dec 6, 2003 ... SIAM Journal on Numerical Analysis is currently published by ... was introduced several years earlier in the mechanical engineering literature While ... tCourant Institute of Mathematical Sciences, New York University, 251 Mercer Street, New .... bid”) = we» — vent-2%] / mews, .... Heat Transfer, 3 (1980), pp.

  11. Coupling of finite element and boundary integral methods for ...

    E-Print Network [OSTI]

    "Peijun Li"

    2009-11-09T23:59:59.000Z

    Oct 8, 2009 ... scattering in both transverse electric and magnetic polarization cases. ..... the equivalent transmission problem (2.7) for loss medium, i.e., r > 0,.

  12. Generalized spectral decomposition method for solving stochastic finite element equations

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    or stochastic bases, have been proposed in order to reduce these computational costs. Recently, a new ap- proach times and memory requirements. The reduction of these computational costs has now become a key question and the automatic construction of reduced bases. In this paper, the concept of GSD, initially introduced for a class

  13. A DISCONTINUOUS FINITE ELEMENT METHOD FOR SOLVING A MULTIWELL PROBLEM

    E-Print Network [OSTI]

    elasticity, shape­memory alloys, materials science. AMS subject classifications. 49M07, 65K10, 65N30, 73C50 in the engineering community, along with the rise of `rational' materials as seen for instance in shape­memory alloys is the simulation of austenite­martensite transformations for shape­memory alloys, and we refer to the work of Ball

  14. Convergence of goal-oriented adaptive finite element methods

    E-Print Network [OSTI]

    Pollock, Sara

    2012-01-01T23:59:59.000Z

    v ? V T . w ? V T . The dual energy estimator is: ? T p (v, T ), v ? V T . The dual energy estimator is: ? T 2 , j (

  15. FINITE ELEMENT METHODS FOR THE SIMULATION OF INCOMPRESSIBLE POWDER FLOW

    E-Print Network [OSTI]

    Hron, Jaroslav

    configurations are provided. PHYSICAL BACKGROUD 1. Mohr-Coulomb criterion for friction The Mohr theory suggests a frictional stress model must be taken into account. This can be done using plasticity and similar theories cannot accurately describe granular flow. It is assumed that the material is incompressible, dry

  16. FLUX RECOVERY FROM PRIMAL HYBRID FINITE ELEMENT METHODS #

    E-Print Network [OSTI]

    Chou, So-Hsiang

    and Statistics, Bowling Green State University, Bowling Green, OH (chou@bgnet.bgsu.edu). The research

  17. A finite element method using singular functions for the Poisson ...

    E-Print Network [OSTI]

    2002-09-05T23:59:59.000Z

    1Department of Mathematics; Purdue University; 1395 Mathematical Science ... Contract=grant sponsor: U.S. Department of Energy by Lawrence Livermore ...... Djaoua M. Equations Int?egrales pour un Probleme Singulier dans le Plan.

  18. Experimental verification of finite element enhancement and MAC methods

    E-Print Network [OSTI]

    Ibrahim, Syed Azman

    1991-01-01T23:59:59.000Z

    for Shaft and Pipe Model 57 12. Measured Natural Frequencies and Damping Factor for Both Models . . 70 13. MAC on FE Mode Shapes - Shaft Model 14. Expanded MAC on FE Mode Shapes ? Shaft Model . . 88 88 15. MAC on FE/Test Mode Shapes ? Shaft Model Data... . . 7. Deformation of Piping System for the Cryogenic Rig at 9. 4 Hertz . . 8. Translational Vibration of the Gear Coupling Rig at 17 Hertz 9. Angular Vibration of the Gear Coupling Rig at 33 Hertz . . . . . . . . 15 41 42 43 44 45 46 10...

  19. Discontinuous Galerkin Subgrid Finite Element Method for Heterogeneous Brinkman's

    E-Print Network [OSTI]

    Lazarov, Raytcho

    made of glass or mineral wool with porosity 0.99, or open foams with porosity above 0.95, see Fig. 1 15 Fig. 1. Microstructure of industrial foams Fig. 2. Microstructure and marcostructure of mineral

  20. Parallel spectral-element direction splitting method for ...

    E-Print Network [OSTI]

    2014-07-03T23:59:59.000Z

    Jun 14, 2014 ... [10] George E. Karniadakis, Spencer J. Sherwin, Spectral/hp Element Methods for Computational Fluid Dynamics, Oxford University Press, ...

  1. NORIA - a finite element computer program for analyzing water, vapor, air, and energy transport in porous media

    SciTech Connect (OSTI)

    Bixler, N.E.

    1985-08-01T23:59:59.000Z

    NORIA is a finite element computer program that solves four nonlinear, parabolic, partial differential equations simultaneously. The four equations describe the transport of water, water vapor, air, and energy through partially saturated porous media. The numerical procedure uses the standard Galerkin finite element method to handle spatial discretization of two-dimensional domains with either planar symmetry or axisymmetry. Time integration is performed by a third-order predictor-corrector scheme that uses error estimates to automatically adjust time-step size so as to maintain uniform local time truncation error throughout the calculation. Thus, the user is not required to select time-step size except at the first time step. Nearly all material properties, such as permeability, can either be set to constant values or can be defined as functions of the dependent and independent variables by user-supplied subroutines. The gas phase is taken to be ideal. This report is intended primarily as a user`s manual but also includes discussions of the theory of two-phase transport in porous media and the numerical procedure used in NORIA. 33 refs., 10 figs., 1 tab.

  2. Determining finite volume elements for the 2D Navier-Stokes equations

    SciTech Connect (OSTI)

    Jones, D.A. (California Univ., Irvine, CA (United States). Dept. of Mathematics); Titi, E.S. (California Univ., Irvine, CA (United States). Dept. of Mathematics Cornell Univ., Ithaca, NY (United States). Mathematical Sciences Inst.)

    1991-01-01T23:59:59.000Z

    We consider the 2D Navier-Stokes equations on a square with periodic boundary conditions. Dividing the square into N equal subsquares, we show that if the asymptotic behavior of the average of solutions on these subsquares (finite volume elements) is known, then the large time behavior of the solution itself is completely determined, provided N is large enough. We also establish a rigorous upper bound for N needed to determine the solutions to the Navier-Stokes equation in terms of the physical parameters of the problem. 34 refs.

  3. Double torsion testing and finite element analysis for determining the electric fracture properties of piezoelectric ceramics

    SciTech Connect (OSTI)

    Shindo, Yasuhide; Narita, Fumio; Mikami, Masaru [Department of Materials Processing, Graduate School of Engineering, Tohoku University, Aoba-yama 6-6-02, Sendai 980-8579 (Japan)

    2005-06-01T23:59:59.000Z

    This paper presents the results of an experimental and numerical investigation in electric fracture behavior of composite [Pb(Zr,Ti)O{sub 3}] double torsion (DT) specimens. DT tests were conducted on a commercial piezoelectric ceramic bonded between two metals. Fracture loads under different electric fields were obtained from the experiment. Nonlinear three-dimensional finite element analysis was also employed to calculate the energy release rate for DT specimens based on the exact (permeable) and approximate (impermeable) crack models. The effects of applied electric field and domain switching on the energy release rate are discussed, and the model predictions are compared with the results of the experiments.

  4. Salinas - An implicit finite element structural dynamics code developed for massively parallel platforms

    SciTech Connect (OSTI)

    BHARDWAJ, MANLJ K.; REESE,GARTH M.; DRIESSEN,BRIAN; ALVIN,KENNETH F.; DAY,DAVID M.

    2000-04-06T23:59:59.000Z

    As computational needs for structural finite element analysis increase, a robust implicit structural dynamics code is needed which can handle millions of degrees of freedom in the model and produce results with quick turn around time. A parallel code is needed to avoid limitations of serial platforms. Salinas is an implicit structural dynamics code specifically designed for massively parallel platforms. It computes the structural response of very large complex structures and provides solutions faster than any existing serial machine. This paper gives a current status of Salinas and uses demonstration problems to show Salinas' performance.

  5. Finite element analysis of surface-stress effects in the Si lattice-parameter measurement

    E-Print Network [OSTI]

    Sasso, Carlo; Massa, Enrico; Mana, Giovanni; Kuetgens, Ulrich

    2013-01-01T23:59:59.000Z

    A stress exists in solids surfaces, similarly to liquids, also if the underlying bulk material is stress-free. This paper investigates the surface stress effect on the measured value of the Si lattice parameter used to determine the Avogadro constant by counting Si atoms. An elastic-film model has been used to provide a surface load in a finite element analysis of the lattice strain of the x-ray interferometer crystal used to measure the lattice parameter. Eventually, an experiment is proposed to work a lattice parameter measurement out so that there is a visible effect of the surface stress.

  6. A finite element analysis of a silicon based double quantum dot structure

    E-Print Network [OSTI]

    S. Rahman; J. Gorman; C. H. W. Barnes; D. A. Williams; H. P. Langtangen

    2006-04-06T23:59:59.000Z

    We present the results of a finite-element solution of the Laplace equation for the silicon-based trench-isolated double quantum-dot and the capacitively-coupled single-electron transistor device architecture. This system is a candidate for charge and spin-based quantum computation in the solid state, as demonstrated by recent coherent-charge oscillation experiments. Our key findings demonstrate control of the electric potential and electric field in the vicinity of the double quantum-dot by the electric potential applied to the in-plane gates. This constitutes a useful theoretical analysis of the silicon-based architecture for quantum information processing applications.

  7. Finite element analysis on the fracture of rubber toughened polymer blends

    SciTech Connect (OSTI)

    Wu, Y.; Mai, Y.W. [Univ. of Sydney, New South Wales (Australia); Wu, J. [Hong Kong Univ. of Science and Technology (Hong Kong)

    1997-12-31T23:59:59.000Z

    The effect of rubber particle volume fraction on the constitutive relation and fracture toughness of polymer blends was studied using elastic-plastic Finite Element Analysis (FEA). The effect of rubber particle cavitation on the stress-strain state at a crack tip was also investigated. Stress analysis reveals that because of the high rubber bulk modulus, the hydrostatic stress inside the rubber particle is close to that in the adjacent matrix material element. As a result, the rubber particle imposes a severe plastic constraint to the surrounding matrix and limits its plastic strain. Rubber particle cavitation can effectively release the constraint and enable large scale plastic strain to occur. Different failure criteria were used to determine the optimum rubber particle volume fraction for the polymer blends studied in this paper.

  8. Adapting a Beam-Based Rotordynamics Model to Accept a General Three-Dimensional Finite-Element Casing Model 

    E-Print Network [OSTI]

    James, Stephen M.

    2011-08-08T23:59:59.000Z

    The subject of this thesis is an extension of a two-dimensional, axisymmetric, Timoshenko-beam finite-element rotordynamic code to include a three-dimensional non-axisymmetric solid-element casing model. Axisymmetric beams are sufficient to model...

  9. Efficient Computation of Info-Gap Robustness for Finite Element Models

    SciTech Connect (OSTI)

    Stull, Christopher J. [Los Alamos National Laboratory; Hemez, Francois M. [Los Alamos National Laboratory; Williams, Brian J. [Los Alamos National Laboratory

    2012-07-05T23:59:59.000Z

    A recent research effort at LANL proposed info-gap decision theory as a framework by which to measure the predictive maturity of numerical models. Info-gap theory explores the trade-offs between accuracy, that is, the extent to which predictions reproduce the physical measurements, and robustness, that is, the extent to which predictions are insensitive to modeling assumptions. Both accuracy and robustness are necessary to demonstrate predictive maturity. However, conducting an info-gap analysis can present a formidable challenge, from the standpoint of the required computational resources. This is because a robustness function requires the resolution of multiple optimization problems. This report offers an alternative, adjoint methodology to assess the info-gap robustness of Ax = b-like numerical models solved for a solution x. Two situations that can arise in structural analysis and design are briefly described and contextualized within the info-gap decision theory framework. The treatments of the info-gap problems, using the adjoint methodology are outlined in detail, and the latter problem is solved for four separate finite element models. As compared to statistical sampling, the proposed methodology offers highly accurate approximations of info-gap robustness functions for the finite element models considered in the report, at a small fraction of the computational cost. It is noted that this report considers only linear systems; a natural follow-on study would extend the methodologies described herein to include nonlinear systems.

  10. Rolling Process Modeling Report: Finite-Element Prediction of Roll Separating Force and Rolling Defects

    SciTech Connect (OSTI)

    Soulami, Ayoub; Lavender, Curt A.; Paxton, Dean M.; Burkes, Douglas

    2014-04-23T23:59:59.000Z

    Pacific Northwest National Laboratory (PNNL) has been investigating manufacturing processes for the uranium-10% molybdenum (U-10Mo) alloy plate-type fuel for the U.S. high-performance research reactors. This work supports the Convert Program of the U.S. Department of Energy’s National Nuclear Security Administration (DOE/NNSA) Global Threat Reduction Initiative. This report documents modeling results of PNNL’s efforts to perform finite-element simulations to predict roll separating forces and rolling defects. Simulations were performed using a finite-element model developed using the commercial code LS-Dyna. Simulations of the hot rolling of U-10Mo coupons encapsulated in low-carbon steel have been conducted following two different schedules. Model predictions of the roll-separation force and roll-pack thicknesses at different stages of the rolling process were compared with experimental measurements. This report discusses various attributes of the rolled coupons revealed by the model (e.g., dog-boning and thickness non-uniformity).

  11. Automation of the matrix element reweighting method

    E-Print Network [OSTI]

    P. Artoisenet; V. Lemaître; F. Maltoni; O. Mattelaer

    2011-02-01T23:59:59.000Z

    Matrix element reweighting is a powerful experimental technique widely employed to maximize the amount of information that can be extracted from a collider data set. We present a procedure that allows to automatically evaluate the weights for any process of interest in the standard model and beyond. Given the initial, intermediate and final state particles, and the transfer functions for the final physics objects, such as leptons, jets, missing transverse energy, our algorithm creates a phase-space mapping designed to efficiently perform the integration of the squared matrix element and the transfer functions. The implementation builds up on MadGraph, it is completely automatized and publicly available. A few sample applications are presented that show the capabilities of the code and illustrate the possibilities for new studies that such an approach opens up.

  12. Method of holding optical elements without deformation during their fabrication

    DOE Patents [OSTI]

    Hed, P. Paul (80 Milaw Ct., San Ramon, CA 94583)

    1997-01-01T23:59:59.000Z

    An improved method for securing and removing an optical element to and from a blocking tool without causing deformation of the optical element. A lens tissue is placed on the top surface of the blocking tool. Dots of UV cement are applied to the lens tissue without any of the dots contacting each other. An optical element is placed on top of the blocking tool with the lens tissue sandwiched therebetween. The UV cement is then cured. After subsequent fabrication steps, the bonded blocking tool, lens tissue, and optical element are placed in a debonding solution to soften the UV cement. The optical element is then removed from the blocking tool.

  13. Finite element analysis of the effect of up-armouring on the off-road braking and sharp-turn performance

    E-Print Network [OSTI]

    Grujicic, Mica

    Finite element analysis of the effect of up-armouring on the off-road braking and sharp-mobility multi-purpose wheeled vehicle, off-road vehicle performance, finite element modelling and simulations revision for publication on 15 June 2009. DOI: 10.1243/09544070JAUTO1187 Abstract: A comprehensive finite

  14. Multigrid and Krylov Solvers for Large Scale Finite Element Groundwater Flow Simulations on Distributed Memory Parallel Platforms

    SciTech Connect (OSTI)

    Mahinthakumar, K.

    1997-01-01T23:59:59.000Z

    In this report we present parallel solvers for large linear systems arising from the finite-element discretization of the three-dimensional steady-state groundwater flow problem. Our solvers are based on multigrid and Krylov subspace methods. The parallel implementation is based on a domain decomposition strategy with explicit message passing using NX and MPI libraries. We have tested our parallel implementations on the Intel Paragon XP/S 150 supercomputer using up to 1024 parallel processors and on other parallel platforms such as SGI/Power Challenge Array, Cray/SGI Origin 2000, Convex Exemplar SPP-1200, and IBM SP using up to 64 processors. We show that multigrid can be a scalable algorithm on distributed memory machines. We demonstrate the effectiveness of parallel multigrid based solvers by solving problems requiring more than 70 million nodes in less than a minute. This is more than 25 times faster than the diagonal preconditioned conjugate gradient method which is one of the more popular methods for large sparse linear systems. Our results also show that multigrid as a stand alone solver works best for problems with smooth coefficients, but for rough coefficients it is best used as a preconditioner for a Krylov subspace method such as the conjugate gradient method. We show that even for extremely heterogeneous systems the multigrid pre-conditioned conjugate gradient method is at least 10 times faster than the diagonally preconditioned conjugate gradient method.

  15. Adaptive Finite element approximation of steady flows of incompressible fluids with implicit power-law-like rheology

    E-Print Network [OSTI]

    Christian Kreuzer; Endre Süli

    2015-03-18T23:59:59.000Z

    We develop the a posteriori error analysis of finite element approximations of implicit power-law-like models for viscous incompressible fluids. The Cauchy stress and the symmetric part of the velocity gradient in the class of models under consideration are related by a, possibly multi--valued, maximal monotone $r$-graph, with $\\frac{2d}{d+1}finite element residual, as well as the local stability of the error bound. We then consider an adaptive finite element approximation of the problem, and, under suitable assumptions, we show the weak convergence of the adaptive algorithm to a weak solution of the boundary-value problem. The argument is based on a variety of weak compactness techniques, including Chacon's biting lemma and a finite element counterpart of the Acerbi--Fusco Lipschitz truncation of Sobolev functions, introduced by L. Diening, C. Kreuzer and E. S\\"uli [Finite element approximation of steady flows of incompressible fluids with implicit power-law-like rheology. SIAM J. Numer. Anal., 51(2), 984--1015].

  16. A fast multipole boundary element method for modeling 2-D multiple crack problems with constant elements

    E-Print Network [OSTI]

    Liu, Yijun

    A fast multipole boundary element method for modeling 2-D multiple crack problems with constant 3 April 2014 Accepted 20 May 2014 Keywords: Fast multipole BEM 2-D multi-crack problems Constant elements Crack opening displacements Stress intensity factors a b s t r a c t A fast multipole boundary

  17. Finite element residual stress analysis of induction heating bended ferritic steel piping

    SciTech Connect (OSTI)

    Kima, Jong Sung [Sunchon National University, 255 Jungang-ro, Sucheon, Jeonnam (Korea, Republic of); Kim, Kyoung-Soo; Oh, Young-Jin; Chang, Hyung-Young; Park, Heung-Bae [KEPCO E and C, Co., Ltd., 188, Kumi-ro, Seongnam, Kyounggi (Korea, Republic of)

    2014-10-06T23:59:59.000Z

    Recently, there is a trend to apply the piping bended by induction heating process to nuclear power plants. Residual stress can be generated due to thermo-mechanical mechanism during the induction heating bending process. It is well-known that the residual stress has important effect on crack initiation and growth. The previous studies have focused on the thickness variation. In part, some studies were performed for residual stress evaluation of the austenitic stainless steel piping bended by induction heating. It is difficult to find the residual stresses of the ferritic steel piping bended by the induction heating. The study assessed the residual stresses of induction heating bended ferriticsteel piping via finite element analysis. As a result, it was identified that high residual stresses are generated on local outersurface region of the induction heating bended ferritic piping.

  18. A parallel high-order accurate finite element nonlinear Stokes ice sheet model and benchmark experiments

    SciTech Connect (OSTI)

    Leng, Wei [Chinese Academy of Sciences; Ju, Lili [University of South Carolina; Gunzburger, Max [Florida State University; Price, Stephen [Los Alamos National Laboratory; Ringler, Todd [Los Alamos National Laboratory,

    2012-01-01T23:59:59.000Z

    The numerical modeling of glacier and ice sheet evolution is a subject of growing interest, in part because of the potential for models to inform estimates of global sea level change. This paper focuses on the development of a numerical model that determines the velocity and pressure fields within an ice sheet. Our numerical model features a high-fidelity mathematical model involving the nonlinear Stokes system and combinations of no-sliding and sliding basal boundary conditions, high-order accurate finite element discretizations based on variable resolution grids, and highly scalable parallel solution strategies, all of which contribute to a numerical model that can achieve accurate velocity and pressure approximations in a highly efficient manner. We demonstrate the accuracy and efficiency of our model by analytical solution tests, established ice sheet benchmark experiments, and comparisons with other well-established ice sheet models.

  19. Finite Element Modeling and Validation of Residual Stresses in 304 L Girth Welds

    SciTech Connect (OSTI)

    Dike, J.J.; Ortega, A.R.; Cadden, C.H.; Rangaswamy, P. Brown, D.

    1998-06-01T23:59:59.000Z

    Three dimensional finite element simulations of thermal and mechanical response of a 304 L stainless steel pipe subjected to a circumferential autogenous gas tungsten arc weld were used to predict residual stresses in the pipe. Energy is input into the thermal model using a volumetric heat source. Temperature histories from the thermal analysis are used as loads in the mechanical analyses. In the mechanical analyses, a state variable constitutive model was used to describe the material behavior. The model accounts for strain rate, temperature, and load path histories. The predicted stresses are compared with x-ray diffraction determinations of residual stress in the hoop and circumferential directions on the outside surface of the pipe. Calculated stress profiles fell within the measured data. Reasons for observed scatter in measured stresses are discussed.

  20. Application of the discontinuous spectral Galerkin method to groundwater flow

    E-Print Network [OSTI]

    Fagherazzi, Sergio

    method uses a finite-element discretization of the groundwater flow domain with basis functions with a standard finite-element method and a mixed finite-element method. Numerical sim- ulations show that the discontinuous spectral Galerkin approximation is more efficient than the standard finite-element method (in

  1. Analysis of Refrigerant Flow and Deformation for a Flexible Short-Tube using a Finite Element Model

    E-Print Network [OSTI]

    O'Neal, D.L.; Bassiouny, R.

    reliability. Short-tubes have either a constant inner dia- meter flow channel or a tapered channel with a smallAbstract A finite element model was used to simulate single-phase flow of R-22 through flexible short-tubes. The numerical model included the fluid... in the flow area. The more flexible (5513 kPa) short-tube restricted the mass flow rate more than the most rigid (9889 kPa) short-tube used in this study. The mass flow rates estimated with the finite element model were as much as 14% higher than those from...

  2. Resolution of grain scale interactions using the Discrete Element Method

    E-Print Network [OSTI]

    Johnson, Scott M. (Scott Matthew), 1978-

    2006-01-01T23:59:59.000Z

    Granular materials are an integral part of many engineering systems. Currently, a popular tool for numerically investigating granular systems is the Discrete Element Method (DEM). Nearly all implementations of the DEM, ...

  3. Use of laboratory triaxial-creep data and finite-element analysis to predict observed creep behavior of leached salt caverns

    SciTech Connect (OSTI)

    Preece, D.S.; Stone, C.M.

    1982-08-01T23:59:59.000Z

    An increasing interest is being shown worldwide in using leached salt caverns to store oil and natural gas. A critical factor in the use of existing caverns and the design of new ones is the creep behavior of the salt surrounding the caverns. An understanding of this behavior is being gained by using laboratory triaxial creep data as material property input to finite element computer programs designed to calculate displacements and stresses due to creep. An important step in verifying these predictive methods is the comparison of field data from existing caverns with finite element analyses which incorporate the material properties and geometry of each site. This comparison has been made for caverns in the Eminence Dome (Mississippi), West Hackberry Dome (Louisiana), and Bayou Chocktaw Dome (Louisiana) with reasonably good correlation being obtained between measured and predicted volumetric response of the caverns. These comparisons are discussed in this paper.

  4. Methods and devices for fabricating and assembling printable semiconductor elements

    DOE Patents [OSTI]

    Nuzzo, Ralph G. (Champaign, IL); Rogers, John A. (Champaign, IL); Menard, Etienne (Urbana, IL); Lee, Keon Jae (Savoy, IL); Khang, Dahl-Young (Urbana, IL); Sun, Yugang (Champaign, IL); Meitl, Matthew (Champaign, IL); Zhu, Zhengtao (Urbana, IL)

    2009-11-24T23:59:59.000Z

    The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

  5. Methods and devices for fabricating and assembling printable semiconductor elements

    DOE Patents [OSTI]

    Nuzzo, Ralph G; Rogers, John A; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao

    2014-03-04T23:59:59.000Z

    The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

  6. Methods and devices for fabricating and assembling printable semiconductor elements

    DOE Patents [OSTI]

    Nuzzo, Ralph G; Rogers, John A; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao

    2013-05-14T23:59:59.000Z

    The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

  7. Methods and devices for fabricating and assembling printable semiconductor elements

    DOE Patents [OSTI]

    Nuzzo, Ralph G. (Champaign, IL); Rogers, John A. (Champaign, IL); Menard, Etienne (Durham, NC); Lee, Keon Jae (Daejeon, KR); Khang, Dahl-Young (Urbana, IL); Sun, Yugang (Champaign, IL); Meitl, Matthew (Raleigh, NC); Zhu, Zhengtao (Urbana, IL)

    2011-07-19T23:59:59.000Z

    The invention provides methods and devices for fabricating printable semiconductor elements and assembling printable semiconductor elements onto substrate surfaces. Methods, devices and device components of the present invention are capable of generating a wide range of flexible electronic and optoelectronic devices and arrays of devices on substrates comprising polymeric materials. The present invention also provides stretchable semiconductor structures and stretchable electronic devices capable of good performance in stretched configurations.

  8. A three dimensional finite element code for the analysis of damage in helically wound composite cylinders

    E-Print Network [OSTI]

    Zocher, Marvin Anthony

    1990-01-01T23:59:59.000Z

    be used for predicting quasi-static crack growth in a class ol' nonlinear viscoelastic materials under finite strain. C'lassicu/ Plate Theory: Schapery and Davidson ' have used classical plate theory (C'PT) in predicting the energy release rate (ERR... plate, shell, or beam. . 1Icthod of C'elis: Aboudi has developed a. method for determination of the effec- tive moduli of infinite elastic solids containing a, doubly periodic rectangular array of cracks. It consists of the analysis of a, unit cell...

  9. A Layerwise Shell Stiffener and Stand-Alone Curved Beam Element

    E-Print Network [OSTI]

    Kassegne, Samuel Kinde

    for the analysis of stiffened plates and shells. The finite element method has been used by a number elements. The element can also be used as a stand-alone curved or straight beam element. A finite element and Schnobrich [1] proposed a sixteen degrees of freedom isotropic beam finite element which has a displacement

  10. A mimetic finite difference method for the Stokes problem with elected edge bubbles

    SciTech Connect (OSTI)

    Lipnikov, K [Los Alamos National Laboratory; Berirao, L [DIPARTMENTO DI MATERMATICA

    2009-01-01T23:59:59.000Z

    A new mimetic finite difference method for the Stokes problem is proposed and analyzed. The unstable P{sub 1}-P{sub 0} discretization is stabilized by adding a small number of bubble functions to selected mesh edges. A simple strategy for selecting such edges is proposed and verified with numerical experiments. The discretizations schemes for Stokes and Navier-Stokes equations must satisfy the celebrated inf-sup (or the LBB) stability condition. The stability condition implies a balance between discrete spaces for velocity and pressure. In finite elements, this balance is frequently achieved by adding bubble functions to the velocity space. The goal of this article is to show that the stabilizing edge bubble functions can be added only to a small set of mesh edges. This results in a smaller algebraic system and potentially in a faster calculations. We employ the mimetic finite difference (MFD) discretization technique that works for general polyhedral meshes and can accomodate non-uniform distribution of stabilizing bubbles.

  11. Coupled Transient Finite Element Simulation of Quench in Jefferson Lab's 11 GeV Super High Momentum Spectrometer Superconducting Magnets

    SciTech Connect (OSTI)

    E. Sun, P. Brindza, S. Lassiter, M. Fowler, E. Xu

    2010-06-01T23:59:59.000Z

    This paper presents coupled transient thermal and electromagnetic finite element analysis of quench in the Q2, Q3, and dipole superconducting magnets using Vector Fields Quench code. Detailed temperature distribution within coils and aluminum force collars were computed at each time step. Both normal (quench with dump resistor) and worst-case (quench without dump resistor) scenarios were simulated to investigate the maximum temperatures. Two simulation methods were utilized, and their algorithms, implementation, advantages, and disadvantages are discussed. The first method simulated the coil using nonlinear transient thermal analysis directly linked with the transient circuit analysis. It was faster because only the coil was meshed and no eddy current was modeled. The second method simulated the whole magnet including the coil, the force collar, and the iron yoke. It coupled thermal analysis with transient electromagnetic field analysis which modeled electromagnetic fields including eddy currents within the force collar. Since eddy currents and temperature in the force collars were calculated in various configurations, segmentation of the force collars was optimized under the condition of fast discharge.

  12. Analysis of Refrigerant Flow and Deformation for a Flexible Short-Tube using a Finite Element Model 

    E-Print Network [OSTI]

    O'Neal, D.L.; Bassiouny, R.

    2004-01-01T23:59:59.000Z

    an short-tube using a finite element model International Journal of Refrig Vol 27 No 2, pp. 176-183The short-tube is used by United States heat pump manufactures as an expansion device because of its lowAnalyse de l’e´coulement d’un frigorige`ne et...

  13. Finite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil and Gas Reservoirs

    E-Print Network [OSTI]

    Patzek, Tadeusz W.

    Finite Element Solution of Nonlinear Transient Rock Damage with Application in Geomechanics of Oil the geomechanics of oil and gas reservoirs. The fragile microstructure of some rocks makes it difficult to predict. The results could be used in similar geomechanical and structural damage problems such as failure and rupture

  14. A ThreeDimensional Finite Element Simulation for Transport of Nuclear Waste Contamination in Porous Media \\Lambda

    E-Print Network [OSTI]

    Ewing, Richard E.

    A Three­Dimensional Finite Element Simulation for Transport of Nuclear Waste Contamination for transport of nuclear­waste contamination in three­dimensional porous media are presented with a description of contamination of groundwater by high­level nuclear waste and a wide variety of other sources makes a proper

  15. Finite element solution of Laplace's equation for ion-atom chambers Jacob Golde, Janine Shertzer, and Paul Oxley

    E-Print Network [OSTI]

    Oxley, Paul

    Finite element solution of Laplace's equation for ion-atom chambers Jacob Golde, Janine Shertzer for ion-atom chambers. We first consider a simplified model chamber for which an analytical solution can is carried out in an ion-atom chamber. The atomic beam first passes through a region of the cham- ber where

  16. Two-dimensional finite element simulation of fracture and fatigue behaviours of alumina microstructures for hip prosthesis

    E-Print Network [OSTI]

    Boyer, Edmond

    microstructures for hip prosthesis Kyungmok Kim¹, Bernard Forest¹, Jean Geringer¹* ¹Ecole Nationale Supérieure des and fatigue behaviours of pure alumina microstructures found at hip prosthesis. Finite element models are similar to those found at a slip zone between a femoral head and an acetabular cup of hip prosthesis

  17. An implicit finite-element model for 3D non-hydrostatic mesoscale ocean M.A. Maidana1

    E-Print Network [OSTI]

    Politècnica de Catalunya, Universitat

    An implicit finite-element model for 3D non-hydrostatic mesoscale ocean flows M.A. Maidana1 , J-dimensional, non-hydrostatic mesoscale ocean flows. The model considered here incorporates surface wind stress and the idea of using unstructured grids for modelling mesoscale ocean dynamics sounds very attractive given

  18. (a finite-element computer code for the computation of visco-plastic creeping flows with applications to tectonics)

    E-Print Network [OSTI]

    Beaumont, Christopher

    SOPALE (a finite-element computer code for the computation of visco-plastic creeping flows are bracketed like this ** 1. What is SOPALE? SOPALE is a computer code written in Fortran 77 which computes of the code history and contributors ** my comments, which may be eliminated or used at your discretion

  19. A two-dimensional finite element thermomechanical approach to a global stress-strain analysis of steel continuous casting

    E-Print Network [OSTI]

    Boyer, Edmond

    of steel continuous casting Michel Bellet, Alban Heinrich Ecole des Mines de Paris, Centre de Mise en Forme.bellet@ensmp.fr Synopsis This paper addresses the two-dimensional finite element simulation of steel continuous casting state (temperature, deformation, stresses) of steel all along the continuous casting machine. Both plane

  20. A characteristic/finite element algorithm for time-dependent 3-D advection-dominated transport using unstructured grids

    E-Print Network [OSTI]

    Mofrad, Mohammad R. K.

    number mass/heat transport in a fully developed pipe flow. Ó 2002 Elsevier Science B.V. All rights of physics and engineering. Examples are heat transfer processes in rubber extrusion, plastic castingA characteristic/finite element algorithm for time-dependent 3-D advection-dominated transport

  1. Finite element analysis of grain-matrix micro-cracking in shale within the context of a multiscale modeling

    E-Print Network [OSTI]

    Regueiro, Richard A.

    , hydrofracturing, or oil shale production. Current macro- scale and multiscale models do not account simultaneouslyFinite element analysis of grain-matrix micro-cracking in shale within the context of a multiscale-cracking in shale at grain-matrix inter- faces, assuming constituents are composed of quart silt grains

  2. Method of holding optical elements without deformation during their fabrication

    DOE Patents [OSTI]

    Hed, P.P.

    1997-04-29T23:59:59.000Z

    An improved method for securing and removing an optical element to and from a blocking tool without causing deformation of the optical element is disclosed. A lens tissue is placed on the top surface of the blocking tool. Dots of UV cement are applied to the lens tissue without any of the dots contacting each other. An optical element is placed on top of the blocking tool with the lens tissue sandwiched therebetween. The UV cement is then cured. After subsequent fabrication steps, the bonded blocking tool, lens tissue, and optical element are placed in a debonding solution to soften the UV cement. The optical element is then removed from the blocking tool. 16 figs.

  3. Method and system for high power reflective optical elements

    DOE Patents [OSTI]

    Demos, Stavros G.; Rubenchik, Alexander M.; Negres, Raluca A.

    2013-03-12T23:59:59.000Z

    A method of repairing damage in an optical element includes providing a laser system including at least one optical element having a coating layer having an incident light surface and directing a laser pulse from the laser system to impinge on the incident light surface. The method also includes sustaining damage to a portion of the incident light surface and melting the damaged portion of the incident light surface and a region adjacent to the damaged portion. The method further includes flowing material from the region adjacent the damaged portion to the damaged portion and solidifying the material in the damaged portion and the region adjacent to the damaged portion.

  4. A MULTIPOLE BOUNDARY ELEMENT METHOD FOR TWO DIMENSIONAL ELASTOSTATICS 3

    E-Print Network [OSTI]

    Yamamoto, Hirosuke

    A MULTIPOLE BOUNDARY ELEMENT METHOD FOR TWO DIMENSIONAL ELASTOSTATICS 3 YOSIHIRO YAMADAy z AND KEN. To avoid these problems, Rokhlin proposed the multipole method for the potential prob- lem. This paper requires O(N log N) work and memory. Theoretical error estimates for the multipole expansions are also

  5. Solving ODEs: Cardiac drug risk assessment The heat equation and finite elements: electrical propagation in the heart Object-oriented scientific computing: Applications

    E-Print Network [OSTI]

    Martin, Ralph R.

    Solving ODEs: Cardiac drug risk assessment The heat equation and finite elements: electrical;Solving ODEs: Cardiac drug risk assessment The heat equation and finite elements: electrical propagation in the heart Solving ODEs: Cardiac drug risk assessment #12;Solving ODEs: Cardiac drug risk assessment The heat

  6. Electrical and thermal finite element modeling of arc faults in photovoltaic bypass diodes.

    SciTech Connect (OSTI)

    Bower, Ward Isaac; Quintana, Michael A.; Johnson, Jay

    2012-01-01T23:59:59.000Z

    Arc faults in photovoltaic (PV) modules have caused multiple rooftop fires. The arc generates a high-temperature plasma that ignites surrounding materials and subsequently spreads the fire to the building structure. While there are many possible locations in PV systems and PV modules where arcs could initiate, bypass diodes have been suspected of triggering arc faults in some modules. In order to understand the electrical and thermal phenomena associated with these events, a finite element model of a busbar and diode was created. Thermoelectrical simulations found Joule and internal diode heating from normal operation would not normally cause bypass diode or solder failures. However, if corrosion increased the contact resistance in the solder connection between the busbar and the diode leads, enough voltage potentially would be established to arc across micron-scale electrode gaps. Lastly, an analytical arc radiation model based on observed data was employed to predicted polymer ignition times. The model predicted polymer materials in the adjacent area of the diode and junction box ignite in less than 0.1 seconds.

  7. Combined simulation of whirl and stick-slip phenomena using a nonlinear finite element model

    SciTech Connect (OSTI)

    Schmalhorst, B. [Baker Hughes INTEQ GmbH, Celle (Germany); Baumgart, A. [Technical Univ. Braunschweig (Germany)

    1997-07-01T23:59:59.000Z

    A time domain drillstring dynamics model allows the prediction of BHA dynamics related to torsional and two plane lateral vibrations. The model considers a nonlinear Finite Element approach that takes into account arbitrarily shaped assemblies during drilling of curved boreholes. Each nodal point represents a stabilizer, a collar connection or a wear knot moving viscously dampled inside open hole or a casing. In the case of wall contact, restitution forces from the formation act on the string and contact velocity dependent friction forces are introduced. The torsional and the lateral degrees of freedom (d.o.f.) are therefore coupled. The authors introduce a nonlinear relationship between the cutting forces at the PDC bit and the angular velocity of the bit. Bit and BHA whirl, rotational Stick-Slip phenomena and normal drilling conditions are solved numerically. Results show that the drilling process is always transient. Commonly used static models only can provide the drilling engineer with rough estimates of the loading. Evaluation of BHA failures sometimes reveal an ultra high dynamic loading of the BHA components that can be analyzed and predicted with this time domain application.

  8. Finite element analyses of continuous filament ties for masonry applications : final report for the Arquin Corporation.

    SciTech Connect (OSTI)

    Quinones, Armando, Sr. (Arquin Corporation, La Luz, NM); Bibeau, Tiffany A.; Ho, Clifford Kuofei

    2008-08-01T23:59:59.000Z

    Finite-element analyses were performed to simulate the response of a hypothetical vertical masonry wall subject to different lateral loads with and without continuous horizontal filament ties laid between rows of concrete blocks. A static loading analysis and cost comparison were also performed to evaluate optimal materials and designs for the spacers affixed to the filaments. Results showed that polypropylene, ABS, and polyethylene (high density) were suitable materials for the spacers based on performance and cost, and the short T-spacer design was optimal based on its performance and functionality. Simulations of vertical walls subject to static loads representing 100 mph winds (0.2 psi) and a seismic event (0.66 psi) showed that the simulated walls performed similarly and adequately when subject to these loads with and without the ties. Additional simulations and tests are required to assess the performance of actual walls with and without the ties under greater loads and more realistic conditions (e.g., cracks, non-linear response).

  9. A stochastic method for computing hadronic matrix elements

    DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)

    Drach, Vincent [DESY, Zeuthen (Germany); Jansen, Karl [DESY, Zeuthen (Germany); Alexandrou, Constantia [University of Cyprus, Nicosia (Cyprus); Constantinou, Martha [University of Cyprus, Nicosia (Cyprus); Dinter, Simon [DESY, Zeuthen (Germany); Hadjiyiannakou, Kyriakos [University of Cyprus, Nicosia (Cyprus); Renner, Dru B. [JLAB, Newport News, VA (United States)

    2014-01-01T23:59:59.000Z

    We present a stochastic method for the calculation of baryon 3-point functions which is an alternative to the typically used sequential method offering more versatility. We analyze the scaling of the error of the stochastically evaluated 3-point function with the lattice volume and find a favorable signal to noise ratio suggesting that the stochastic method can be extended to large volumes providing an efficient approach to compute hadronic matrix elements and form factors.

  10. Discrete Element Modeling

    SciTech Connect (OSTI)

    Morris, J; Johnson, S

    2007-12-03T23:59:59.000Z

    The Distinct Element Method (also frequently referred to as the Discrete Element Method) (DEM) is a Lagrangian numerical technique where the computational domain consists of discrete solid elements which interact via compliant contacts. This can be contrasted with Finite Element Methods where the computational domain is assumed to represent a continuum (although many modern implementations of the FEM can accommodate some Distinct Element capabilities). Often the terms Discrete Element Method and Distinct Element Method are used interchangeably in the literature, although Cundall and Hart (1992) suggested that Discrete Element Methods should be a more inclusive term covering Distinct Element Methods, Displacement Discontinuity Analysis and Modal Methods. In this work, DEM specifically refers to the Distinct Element Method, where the discrete elements interact via compliant contacts, in contrast with Displacement Discontinuity Analysis where the contacts are rigid and all compliance is taken up by the adjacent intact material.

  11. Finite element analysis and design of large diameter flexible vertical pipes subjected to incremental compacted backfill loads and creep effects

    E-Print Network [OSTI]

    Hossain, Mohammad Kabir

    1990-01-01T23:59:59.000Z

    FINITE ELEMENT ANALYSIS AND DESiGN OF LARGE DIANETER FLEXIBLE VERTICAL PIPES SUBJECTED TO INCREMENTAL COMPACTED BACKFILL LOADS AND CREEP EFFECTS A Thesis by MOHAMMAD KABIR HOSSAIN Submitted to the Office of Graduate Studies of Texas A... BACKFILL LOADS AND CREEP EFFECTS A Thesis by MOHAMMAD KABIR HOSSAIN Approved as to sty1e and content by: R be+r L. Lytton (Chair of Committee) Ozden 0. Ochoa (Member) Derek V. Morris (Member) ames T P Yao (Head of rtment) ABSTRACT Finite...

  12. Reproducing kernel element method Part III: Generalized enrichment and applications

    E-Print Network [OSTI]

    Li, Shaofan

    Reproducing kernel element method Part III: Generalized enrichment and applications Hongsheng Lu enrichment is proposed to construct the global partition polynomials or to enrich global partition polynomial. This is accomplished by either multiplying enrichment functions with the original global partition poly- nomials

  13. A benchmark study of 2D and 3D finite element calculations simulating dynamic pulse buckling tests of cylindrical shells under axial impact

    SciTech Connect (OSTI)

    Hoffman, E.L.; Ammerman, D.J.

    1993-08-01T23:59:59.000Z

    A series of tests investigating dynamic pulse buckling of a cylindrical shell under axial impact is compared to several finite element simulations of the event. The purpose of the study is to compare the performance of the various analysis codes and element types with respect to a problem which is applicable to radioactive material transport packages, and ultimately to develop a benchmark problem to qualify finite element analysis codes for the transport package design industry.

  14. FINITE ELEMENT ANALYSIS OF JNES/NUPEC SEISMIC SHEAR WALL CYCLIC AND SHAKING TABLE TEST DATA.

    SciTech Connect (OSTI)

    XU,J.; NIE, J.; HOFMAYER, C.; ALI, S.

    2007-04-12T23:59:59.000Z

    This paper describes a finite element analysis to predict the JNES/NUPEC cyclic and shaking table RC shear wall test data, as part of a collaborative agreement between the U.S. NRC and JNES to study seismic issues important to the safe operation of commercial nuclear power plant (NPP) structures, systems and components (SSC). The analyses described in this paper were performed using ANACAP reinforced concrete models. The paper describes the ANACAP analysis models and discusses the analysis comparisons with the test data. The ANACAP capability for modeling nonlinear cyclic characteristics of reinforced concrete shear wall structures was confirmed by the close comparisons between the ANACAP analysis results and the JNES/NUPEC cyclic test data. Reasonable agreement between the analysis results and the test data was demonstrated for the hysteresis loops and the shear force orbits, in terms of both the overall shape and the cycle-to-cycle comparisons. The ANACAP simulation analysis of the JNES/NUPEC shaking table test was also performed, which demonstrated that the ANACAP dynamic analysis with concrete material model is able to capture the progressive degrading behavior of the shear wall as indicated from the test data. The ANACAP analysis also predicted the incipient failure of the shear wall, reasonably close to the actual failure declared for the test specimen. In summary, the analyses of the JNES/NUPEC cyclic and shaking table RC shear wall tests presented in this paper have demonstrated the state-of-the-art analysis capability for determining the seismic capacity of RC shear wall structures.

  15. Best Practices for Finite Element Analysis of Spent Nuclear Fuel Transfer, Storage, and Transportation Systems

    SciTech Connect (OSTI)

    Bajwa, Christopher S.; Piotter, Jason; Cuta, Judith M.; Adkins, Harold E.; Klymyshyn, Nicholas A.; Fort, James A.; Suffield, Sarah R.

    2010-08-11T23:59:59.000Z

    Storage casks and transportation packages for spent nuclear fuel (SNF) are designed to confine SNF in sealed canisters or casks, provide structural integrity during accidents, and remove decay through a storage or transportation overpack. The transfer, storage, and transportation of SNF in dry storage casks and transport packages is regulated under 10 CFR Part 72 and 10 CFR Part 71, respectively. Finite Element Analysis (FEA) is used with increasing frequency in Safety Analysis Reports and other regulatory technical evaluations related to SNF casks and packages and their associated systems. Advances in computing power have made increasingly sophisticated FEA models more feasible, and as a result, the need for careful review of such models has also increased. This paper identifies best practice recommendations that stem from recent NRC review experience. The scope covers issues common to all commercially available FEA software, and the recommendations are applicable to any FEA software package. Three specific topics are addressed: general FEA practices, issues specific to thermal analyses, and issues specific to structural analyses. General FEA practices covers appropriate documentation of the model and results, which is important for an efficient review process. The thermal analysis best practices are related to cask analysis for steady state conditions and transient scenarios. The structural analysis best practices are related to the analysis of casks and associated payload during standard handling and drop scenarios. The best practices described in this paper are intended to identify FEA modeling issues and provide insights that can help minimize associated uncertainties and errors, in order to facilitate the NRC licensing review process.

  16. Two-dimensional finite element simulation of fracture and fatigue behaviours of alumina microstructures for hip prosthesis

    E-Print Network [OSTI]

    Kim, Kyungmok; Géringer, Jean; 10.1177/0954411911422843

    2012-01-01T23:59:59.000Z

    This paper describes a two-dimensional (2D) finite element simulation for fracture and fatigue behaviours of pure alumina microstructures such as those found at hip prostheses. Finite element models are developed using actual Al2O3 microstructures and a bilinear cohesive zone law. Simulation conditions are similar to those found at a slip zone in a dry contact between a femoral head and an acetabular cup of hip prosthesis. Contact stresses are imposed to generate cracks in the models. Magnitudes of imposed stresses are higher than those found at the microscopic scale. Effects of microstructures and contact stresses are investigated in terms of crack formation. In addition, fatigue behaviour of the microstructure is determined by performing simulations under cyclic loading conditions. It is shown that crack density observed in a microstructure increases with increasing magnitude of applied contact stress. Moreover, crack density increases linearly with respect to the number of fatigue cycles within a given con...

  17. Method and apparatus for diagnosing breached fuel elements

    DOE Patents [OSTI]

    Gross, K.C.; Lambert, J.D.B.; Nomura, S.

    1987-03-02T23:59:59.000Z

    The invention provides an apparatus and method for diagnosing breached fuel elements in a nuclear reactor. A detection system measures the activity of isotopes from the cover gas in the reactor. A data acquisition and processing system monitors the detection system and corrects for the effects of the cover-gas clean up system on the measured activity and further calculates the derivative curve of the corrected activity as a function of time. A plotting system graphs the derivative curve, which represents the instantaneous release rate of fission gas from a breached fuel element. 8 figs.

  18. Method and apparatus for diagnosing breached fuel elements

    DOE Patents [OSTI]

    Gross, Kenny C. (Lemont, IL); Lambert, John D. B. (Wheston, IL); Nomura, Shigeo (Mito, JP)

    1988-01-01T23:59:59.000Z

    The invention provides an apparatus and method for diagnosing breached fuel elements in a nuclear reactor. A detection system measures the activity of isotopes from the cover-gas in the reactor. A data acquisition and processing system monitors the detection system and corrects for the effects of the cover-gas clean up system on the measured activity and further calculates the derivative cure of the corrected activity as a function of time. A plotting system graphs the derivative curve, which represents the instantaneous release rate of fission gas from a breached fuel element.

  19. Finite element simulations of compositionally graded InGaN solar cells G.F. Brown a,b,n

    E-Print Network [OSTI]

    Wu, Junqiao

    , a highly conductive p-type GaN layer provides the hole contact while absorption takes place in the lowerFinite element simulations of compositionally graded InGaN solar cells G.F. Brown a,b,n , J.W. Ager Keywords: Device modeling InGaN Composition grading Heterojunction a b s t r a c t The solar power

  20. Unstructured spectral element methods of simulation of turbulent flows

    SciTech Connect (OSTI)

    Henderson, R.D. [California Inst. of Technology, Pasadena, CA (United States)] [California Inst. of Technology, Pasadena, CA (United States); Karniadakis, G.E. [Brown Univ., Providence, RI (United States)] [Brown Univ., Providence, RI (United States)

    1995-12-01T23:59:59.000Z

    In this paper we present a spectral element-Fourier algorithm for simulating incompressible turbulent flows in complex geometries using unstructured quadrilateral meshes. To this end, we compare two different interface formulations for extending the conforming spectral element method in order to allow for surgical mesh refinement and still retain spectral accuracy: the Zanolli iterative procedure and variational patching based on auxiliary {open_quotes}mortar{close_quotes} functions. We present an interpretation of the original mortar element method as a patching scheme and develop direct and iterative solution techniques that make the method efficient for simulations of turbulent flows. The properties of the new method are analyzed in detail by studying the eigenspectra of the advection and diffusion operators. We then present numerical results that illustrate the flexibility as well as the exponential convergence of the new algorithm for nonconforming discretizations. We conclude with simulation studies of the turbulent cylinder wake at Re = 1000 (external flow) and turbulent flow over riblets at Re = 3280 (internal flow). 36 refs., 29 figs., 7 tabs.

  1. A posteriori error analysis for a finite volume method

    E-Print Network [OSTI]

    Tavener, Simon

    of nuclear fuel rods, electrical impedance tomography · Pose well known challenges for accurate numerical interface method · Fictitious domain method · Ghost fluid method · (Immersed boundary method

  2. Adaptive Algebraic Multigrid for Finite Element Elliptic Equations with Random Coefficients

    SciTech Connect (OSTI)

    Kalchev, D

    2012-04-02T23:59:59.000Z

    This thesis presents a two-grid algorithm based on Smoothed Aggregation Spectral Element Agglomeration Algebraic Multigrid (SA-{rho}AMGe) combined with adaptation. The aim is to build an efficient solver for the linear systems arising from discretization of second-order elliptic partial differential equations (PDEs) with stochastic coefficients. Examples include PDEs that model subsurface flow with random permeability field. During a Markov Chain Monte Carlo (MCMC) simulation process, that draws PDE coefficient samples from a certain distribution, the PDE coefficients change, hence the resulting linear systems to be solved change. At every such step the system (discretized PDE) needs to be solved and the computed solution used to evaluate some functional(s) of interest that then determine if the coefficient sample is acceptable or not. The MCMC process is hence computationally intensive and requires the solvers used to be efficient and fast. This fact that at every step of MCMC the resulting linear system changes, makes an already existing solver built for the old problem perhaps not as efficient for the problem corresponding to the new sampled coefficient. This motivates the main goal of our study, namely, to adapt an already existing solver to handle the problem (with changed coefficient) with the objective to achieve this goal to be faster and more efficient than building a completely new solver from scratch. Our approach utilizes the local element matrices (for the problem with changed coefficients) to build local problems associated with constructed by the method agglomerated elements (a set of subdomains that cover the given computational domain). We solve a generalized eigenproblem for each set in a subspace spanned by the previous local coarse space (used for the old solver) and a vector, component of the error, that the old solver cannot handle. A portion of the spectrum of these local eigen-problems (corresponding to eigenvalues close to zero) form the coarse basis used to define the new two-level method of our interest. We illustrate the performance of this adaptive two-level procedure with a large set of numerical experiments that demonstrate its efficiency over building the solvers from scratch.

  3. Method for measuring recovery of catalytic elements from fuel cells

    DOE Patents [OSTI]

    Shore, Lawrence (Edison, NJ); Matlin, Ramail (Berkeley, NJ)

    2011-03-08T23:59:59.000Z

    A method is provided for measuring the concentration of a catalytic clement in a fuel cell powder. The method includes depositing on a porous substrate at least one layer of a powder mixture comprising the fuel cell powder and an internal standard material, ablating a sample of the powder mixture using a laser, and vaporizing the sample using an inductively coupled plasma. A normalized concentration of catalytic element in the sample is determined by quantifying the intensity of a first signal correlated to the amount of catalytic element in the sample, quantifying the intensity of a second signal correlated to the amount of internal standard material in the sample, and using a ratio of the first signal intensity to the second signal intensity to cancel out the effects of sample size.

  4. A Comparison of Least-Squares Finite Element Models with the Conventional Finite Element Models of Problems in Heat Transfer and Fluid Mechanics

    E-Print Network [OSTI]

    Nellie Rajarova,

    2010-07-14T23:59:59.000Z

    disadvantage of classical variational methods is the difficulty in constructing the approximate functions for arbitrary domains. 1.4. Review of Weighted-Residual Methods Weighted-residual methods are those in which we seek approximate solutions using a... functions. The pth order interpolation function is given by )48.2( 1 2 1 1 2 1 221 2 21 T pp npp p kp g g g f f f...

  5. Method for detection of antibodies for metallic elements

    DOE Patents [OSTI]

    Barrick, Charles W. (Golden, CO); Clarke, Sara M. (Larkspur, CO); Nordin, Carl W. (Golden, CO)

    1993-11-30T23:59:59.000Z

    An apparatus and method for detecting antibodies specific to non-protein antigens. The apparatus is an immunological plate containing a plurality of plastic projections coated with a non-protein material. Assays utilizing the plate are capable of stabilizing the non-protein antigens with detection levels for antibodies specific to the antigens on a nanogram level. A screening assay with the apparatus allows for early detection of exposure to non-protein materials. Specifically metallic elements are detected.

  6. Rotordynamic analysis with shell elements for the transfer matrix method

    E-Print Network [OSTI]

    L'Antigua, Edward Anthony

    1989-01-01T23:59:59.000Z

    Committee) Al B. elazzo1o (Member) John T. Tielki g (Member) Michael J. Rabins (Head of Department) August 1989 111 ABSTRACT Rotordynamic Analysis with Shell Elements for the Transi'er Matrix Method. (August 1989) Edward Anthony L'Antigua, B. M... Varying Bending Load with (-') = . 01, and Varying ( ? ) ratio. Table 2. Nondimensional Flexibilities (aEv ) for Cylindrical Sections (8 = 0) with Varying ( ? ) and (?-) Ratios. Table 3. Nondimensional Flexibilities (bEr = oEr ) for Cylindrical...

  7. Method for detection of antibodies for metallic elements

    DOE Patents [OSTI]

    Barrick, C.W.; Clarke, S.M.; Nordin, C.W.

    1993-11-30T23:59:59.000Z

    An apparatus and method for detecting antibodies specific to non-protein antigens. The apparatus is an immunological plate containing a plurality of plastic projections coated with a non-protein material. Assays utilizing the plate are capable of stabilizing the non-protein antigens with detection levels for antibodies specific to the antigens on a nanogram level. A screening assay with the apparatus allows for early detection of exposure to non-protein materials. Specifically metallic elements are detected. 10 figures.

  8. PetrovGalerkin Spectral Elements 1 A PetrovGalerkin Spectral Element Technique for Heterogeneous Porous

    E-Print Network [OSTI]

    in spirit to a mixed finite element method or a Petrov­Galerkin scheme. If the sequence of orthogonal commonly used for finite element approximations [14]. The technique proceeds by integrating against ``testPetrov­Galerkin Spectral Elements 1 A Petrov­Galerkin Spectral Element Technique for Heterogeneous

  9. Key Elements of and Materials Performance Targets for Highly Insulating Window Frames

    E-Print Network [OSTI]

    Gustavsen, Arild

    2012-01-01T23:59:59.000Z

    finite-element method (FEM) program, THERM (Finlayson et al. , 1998), to solve the conductive heat-transfer

  10. Towards the application of the Maximum Entropy Method to finite temperature Upsilon Spectroscopy

    E-Print Network [OSTI]

    M. Oevers; C. Davies; J. Shigemitsu

    2000-09-22T23:59:59.000Z

    According to the Narnhofer Thirring Theorem interacting systems at finite temperature cannot be described by particles with a sharp dispersion law. It is therefore mandatory to develop new methods to extract particle masses at finite temperature. The Maximum Entropy method offers a path to obtain the spectral function of a particle correlation function directly. We have implemented the method and tested it with zero temperature Upsilon correlation functions obtained from an NRQCD simulation. Results for different smearing functions are discussed.

  11. Efficient numerical methods for capacitance extraction based on boundary element method

    E-Print Network [OSTI]

    Yan, Shu

    2006-04-12T23:59:59.000Z

    algorithms based on boundary element methods (BEM) and to compute the capacitance extraction in the presence of floating dummy conductors. We propose the PHiCap algorithm, which is based on a hierarchical refinement algorithm and the wavelet transform. Unlike...

  12. Benchmark on Anisotropic Problems Numerical investigation of a mimetic finite difference method

    E-Print Network [OSTI]

    Cirpka, Olaf Arie

    Bernd Flemisch, Rainer Helmig Universität Stuttgart Institute of Hydraulic Engineering Department ourselves, we only provide a brief and rough overview in the following. For more details, we would like is given by Algorithm 1 in [BRE 05b]. By usual hybridization known from mixed finite elements, system [2

  13. A Design-Oriented Framework to Determine the Parasitic Parameters of High Frequency Magnetics in Switching Power Supplies using Finite Element Analysis Techniques

    E-Print Network [OSTI]

    Shadmand, Mohammad

    2012-07-16T23:59:59.000Z

    A DESIGN-ORIENTED FRAMEWORK TO DETERMINE THE PARASITIC PARAMETERS OF HIGH FREQUENCY MAGNETICS IN SWITCING POWER SUPPLIES USING FINITE ELEMENT ANALYSIS TECHNIQUES A Thesis by MOHAMMAD BAGHER SHADMAND Submitted to the Office... to Determine the Parasitic Parameters of High Frequency Magnetics in Switching Power Supplies using Finite Element Analysis Techniques Copyright 2012 Mohammad Bagher Shadmand A DESIGN-ORIENTED FRAMEWORK TO DETERMINE THE PARASITIC PARAMETERS OF HIGH...

  14. A non-linear constrained optimization technique for the mimetic finite difference method

    SciTech Connect (OSTI)

    Manzini, Gianmarco [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Svyatskiy, Daniil [Los Alamos National Lab. (LANL), Los Alamos, NM (United States); Bertolazzi, Enrico [Univ. of Trento (Italy); Frego, Marco [Univ. of Trento (Italy)

    2014-09-30T23:59:59.000Z

    This is a strategy for the construction of monotone schemes in the framework of the mimetic finite difference method for the approximation of diffusion problems on unstructured polygonal and polyhedral meshes.

  15. Multicomponent fluid flow by discontinuous Galerkin and mixed methods in unfractured and fractured media

    E-Print Network [OSTI]

    Firoozabadi, Abbas

    , and fractured media is presented in single phase. In the numerical model we combine the mixed finite element), finite volume (FV) or finite element (FE) methods generally do not provide satisfactory results combined the mixed finite element (MFE) method and the discontin- uous Galerkin (DG) method. The MFE method is

  16. The Distinct Element Method - Application to Structures in Jointed Rock

    SciTech Connect (OSTI)

    Morris, J.P.; Glen, L.; Blair, S.; Heuze, F.

    2001-11-30T23:59:59.000Z

    The Distinct Element Method (DEM) is a meshfree method with applications to rock mechanics, mining sciences, simulations of nuclear repositories, and the stability of underground structures. Continuum mesh-based methods have been applied successfully to many problems in geophysics. Even if the geology includes fractures and faults, when sufficiently large length scales are considered a continuum approximation may be sufficient. However, a large class of problems exist where individual rock joints must be taken into account. This includes problems where the structures of interest have sizes comparable with the block size. In addition, it is possible that while the structure may experience loads which do no measurable damage to individual blocks, some joints may fail. This may launch smaller blocks as dangerous projectiles or even cause total failure of a tunnel. Traditional grid-based continuum approaches are wholly unsuited to this class of problem. It is possible to introduce discontinuities or slide lines into existing grid-based methods, however, such limited approaches can break down when new contacts form between blocks. The distinct element method (DEM) is an alternative, meshfree approach. The DEM can directly approximate the block structure of the jointed rock using arbitrary polyhedra. Using this approach, preexisting joints are readily incorporated into the DEM model. In addition, the method detects all new contacts between blocks resulting from relative block motion. We will describe the background of the DEM and review previous application of the DEM to geophysical problems. Finally we present preliminary results from a investigation into the stability of underground structures subjected to dynamic loading.

  17. Multi-element probabilistic collocation method in high dimensions

    SciTech Connect (OSTI)

    Foo, Jasmine [Division of Applied Mathematics, Brown University, 182 George St., Box F, Providence, RI 02912 (United States); Karniadakis, George Em [Division of Applied Mathematics, Brown University, 182 George St., Box F, Providence, RI 02912 (United States)], E-mail: gk@dam.brown.edu

    2010-03-01T23:59:59.000Z

    We combine multi-element polynomial chaos with analysis of variance (ANOVA) functional decomposition to enhance the convergence rate of polynomial chaos in high dimensions and in problems with low stochastic regularity. Specifically, we employ the multi-element probabilistic collocation method MEPCM and so we refer to the new method as MEPCM-A. We investigate the dependence of the convergence of MEPCM-A on two decomposition parameters, the polynomial order {mu} and the effective dimension {nu}, with {nu}<=}{mu} for monotonic convergence of the method. We also employ MEPCM-A to obtain error bars for the piezometric head at the Hanford nuclear waste site under stochastic hydraulic conductivity conditions. Finally, we compare the cost of MEPCM-A against Monte Carlo in several hundred dimensions, and we find MEPCM-A to be more efficient for up to 600 dimensions for a specific multi-dimensional integration problem involving a discontinuous function.

  18. Brake rotor design and comparison using finite element analysis : an investigation in topology optimization

    E-Print Network [OSTI]

    Domond, Kenneth

    2010-01-01T23:59:59.000Z

    An Introduction to the Conjugate Gradient Method Without thedescent method and the conjugate gradient method. The onlymethod and the conjugate gradient method can be found in

  19. Effect of substrate compliance on the global unilateral post-buckling of coatings: AFM observations and finite element calculations

    SciTech Connect (OSTI)

    Parry, G. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex, (France)]. E-mail: guillaume.parry@etu.univ-poitiers.fr; Colin, J. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex (France); Coupeau, C. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex (France); Foucher, F. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex (France); Cimetiere, A. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex (France); Grilhe, J. [Laboratoire de Metallurgie Physique, UMR 6630, du CNRS, Universite de Poitiers, BP 30179, 86962 Futuroscope Cedex (France)

    2005-01-10T23:59:59.000Z

    The post-critical regime of straight-sided wrinkles on compliant substrates of polycarbonate has been observed by atomic force microscope and investigated by means of finite element simulations. The effect of coupling between the film and its substrate has revealed a global buckling phenomenon, characterized by critical loads lower than those found in the case of a rigid substrate. Characteristic shapes of the buckled structure have been also found to spread over a region wider than the delaminated zone itself. A law relating the film deflexion to the stress has finally been established for any film/substrate system.

  20. LOWER TRANSIENT STRESSES IN AN AORTIC VALVE LEAFLET WITH OBLIQUE REINFORCEMENT FIBERS: A FINITE ELEMENT STUDY

    E-Print Network [OSTI]

    LOWER TRANSIENT STRESSES IN AN AORTIC VALVE LEAFLET WITH OBLIQUE REINFORCEMENT FIBERS: A FINITE Cambridge, MA, USA INTRODUCTION Native aortic valve leaflets contain a network of strong collagen fibers the stresses of cyclic loading over billions of heart beats. Replacement valves that consist of animal tissues

  1. A fast multipole boundary element method for solving the thin plate bending problem

    E-Print Network [OSTI]

    Liu, Yijun

    A fast multipole boundary element method for solving the thin plate bending problem S. Huang, Y: Fast multipole method Boundary element method Thin plate bending problem a b s t r a c t A fast multipole boundary element method (BEM) for solving large-scale thin plate bending problems is presented

  2. Implementation of a Finite Di erence Method on a Custom Computing Platform

    E-Print Network [OSTI]

    erence method a two-dimensional heat transfer modeling system using a Splash-2 Con gurable ComputingImplementation of a Finite Di erence Method on a Custom Computing Platform Kevin J. Paar and Peter involvingcomplex boundary conditions. The nite di erence method can be used to perform heat and shock analysis

  3. REMARKS ON THE MAXIMUM ENTROPY METHOD APPLIED TO FINITE TEMPERATURE LATTICE QCD.

    SciTech Connect (OSTI)

    UMEDA, T.; MATSUFURU, H.

    2005-07-25T23:59:59.000Z

    We make remarks on the Maximum Entropy Method (MEM) for studies of the spectral function of hadronic correlators in finite temperature lattice QCD. We discuss the virtues and subtlety of MEM in the cases that one does not have enough number of data points such as at finite temperature. Taking these points into account, we suggest several tests which one should examine to keep the reliability for the results, and also apply them using mock and lattice QCD data.

  4. Remarks on the Maximum Entropy Method applied to finite temperature lattice QCD

    E-Print Network [OSTI]

    Takashi Umeda; Hideo Matsufuru

    2005-10-05T23:59:59.000Z

    We make remarks on the Maximum Entropy Method (MEM) for studies of the spectral function of hadronic correlators in finite temperature lattice QCD. We discuss the virtues and subtlety of MEM in the cases that one does not have enough number of data points such as at finite temperature. Taking these points into account, we suggest several tests which one should examine to keep the reliability for the results, and also apply them using mock and lattice QCD data.

  5. A new Cement to Glue non-conforming Grids with Robin interface conditions: the finite element case

    E-Print Network [OSTI]

    Japhet, Caroline; Nataf, Frédéric

    2007-01-01T23:59:59.000Z

    We design and analyze a new non-conforming domain decomposition method based on Schwarz type approaches that allows for the use of Robin interface conditions on non-conforming grids. The method is proven to be well posed, and the iterative solver to converge. The error analysis is performed in 2D piecewise polynomials of low and high order and extended in 3D for $P_1$ elements. Numerical results in 2D illustrate the new method.

  6. Two dimensional finite element analysis of homogenization and failure in plain weave textile composites 

    E-Print Network [OSTI]

    Kondagunta, Veeraraghava Gopal

    1993-01-01T23:59:59.000Z

    contains one repeating pattern of fabric design. The elastic analysis of one such structure yields all the 3D moduli of the composite. There have also been a few papers in the area of damage mechanics of textile composites. Li et. al. [21] used 2D finite... are considered. The understanding gained from this two dimen- sional study can be used to extend these techniques for a complete 3D analysis. The high degree of inhomogeneity in textile composites due to the coarse microstructure makes the O ~ eel 0 Cl CD...

  7. Mixed finite element discretization of a model for organic pollution in waters

    E-Print Network [OSTI]

    Paris-Sud XI, Université de

    of the DO and its flux at the boundary. Currently, no oxygen flux is generated by the environmental medium the Biochemical Oxygen Demand and Dissolved Oxygen to which dispersion terms are added. In this work, we propose. The central element of such problems is the oxygen. The main tracers currently used are the density b

  8. Higher-Order Spectral/HP Finite Element Technology for Structures and Fluid Flows

    E-Print Network [OSTI]

    Vallala, Venkat Pradeep

    2013-06-20T23:59:59.000Z

    This study deals with the use of high-order spectral/hp approximation functions in the ?nite element models of various nonlinear boundary-value and initial-value problems arising in the ?elds of structural mechanics and ?ows of viscous...

  9. New finite elements with embedded strong discontinuities to model failure of three-dimensional continua

    E-Print Network [OSTI]

    Kim, Jongheon

    2013-01-01T23:59:59.000Z

    An- alytical Methods in Geomechanics, 31(2):239–259. DvorkinAnalytical Methods in Geomechanics, 28:609–632. Oliver J,and Analytical Methods in Geomechanics, 2:99–128. Ventura G,

  10. Wave propagation across acoustic / Biot's media: a finite-difference method

    E-Print Network [OSTI]

    Boyer, Edmond

    Wave propagation across acoustic / Biot's media: a finite-difference method Guillaume Chiavassa1 Marseille, France. Abstract. Numerical methods are developed to simulate the wave propagation in het- erogeneous 2D fluid / poroelastic media. Wave propagation is described by the usual acoustics equations (in

  11. A Mortar Mimetic Finite Difference Method on NonMatching Grids

    E-Print Network [OSTI]

    Berndt, Markus

    , and Department of Petroleum and Geosystems Engineering, The University of Texas at Austin, Austin, TX 78712, U intersecting at each vertex [15]. The methods has also been extended to locally refined meshes with hanging for the normal velocities in MFD methods on h 2 ­uniform quadrilateral meshes was established. Mixed finite

  12. Generalized cell-centered finite volume methods: application to two-phase flow in porous media

    E-Print Network [OSTI]

    techniques have been used in nuclear engineering for the diffusion approximation and are called nodal methods hand, cell-centered finite volume methods are widely used by engineers and scien- tists who have., (John Wiley, Chichester, England,1997), pp. 231-241. #12;resulting schemes may be different from what

  13. A finite element model for transient thermal/structural analysis of large composite space structures

    E-Print Network [OSTI]

    Lutz, James Delmar

    1986-01-01T23:59:59.000Z

    Versus Time at Solar and Antisolar Points For Cross-Section Number 2 84 28 Temperature Versus Time at Solar and Antisolar Points f' or Cross-Section Number 3 . . . . . . . . . . . . . 85 29 Thermal Noments Versus Time 87 30 Thermal Axial Forces... ar bitrar y domain xi v Superscripts conduction element iteration number flux radiation thermal Subscripts initial reference solar unstrained state INTRODVCTION In the next few decades large space structures will be placed into earth...

  14. Enhanced finite element analysis using the PATRAN-G pre- and post processor program

    E-Print Network [OSTI]

    Murry, Margaret Lois

    1988-01-01T23:59:59.000Z

    better and it was concluded that color graphics saved time in the verification of the model. Clarke and Muller conducted a study in which they broke down the cost of a typical analytical solution into twenty percent computing costs and eighty percent... plots, orthographic plots, view control, window control, a color table, plot labeling, and solid contour plots. Other PATRAN-G features are element connectivity optimization, interactive editing of the model, neutral file interfacing, and a large help...

  15. Numerical study on transient heat transfer under soil with plastic mulch in agriculture applications using a nonlinear finite element model

    E-Print Network [OSTI]

    De Castro, Carlos Armando

    2011-01-01T23:59:59.000Z

    In this paper is developed a simple mathematical model of transient heat transfer under soil with plastic mulch in order to determine with numerical studies the influence of different plastic mulches on the soil temperature and the evolutions of temperatures at different depths with time. The governing differential equations are solved by a Galerkin Finite Element Model, taking into account the nonlinearities due to radiative heat exchange between the soil surface, the plastic mulch and the atmosphere. The model was validated experimentally giving good approximation of the model to the measured data. Simulations were run with the validated model in order to determine the optimal combination of mulch optical properties to maximize the soil temperature with a Taguchi's analysis, proving that the material most used nowadays in Colombia is not the optimal and giving quantitative results of the properties the optimal mulch must possess.

  16. Convergence analysis of sectional methods for solving aggregation population balance equations: The fixed pivot technique

    E-Print Network [OSTI]

    Magdeburg, Universität

    element methods [11, 30, 32], finite volume methods [3, 12, 13], stochastic methods [8, 9, 26], moment methods [34] and sectional methods [17, 18, 22, 23]. A large variety of finite element methods, finite volume methods, weighted residuals, the method of orthogonal collocation and Galerkin's method

  17. Anisotropic mesh adaptation for solution of finite element problems using hierarchical edge-based error estimates

    SciTech Connect (OSTI)

    Lipnikov, Konstantin [Los Alamos National Laboratory; Agouzal, Abdellatif [UNIV DE LYON; Vassilevski, Yuri [Los Alamos National Laboratory

    2009-01-01T23:59:59.000Z

    We present a new technology for generating meshes minimizing the interpolation and discretization errors or their gradients. The key element of this methodology is construction of a space metric from edge-based error estimates. For a mesh with N{sub h} triangles, the error is proportional to N{sub h}{sup -1} and the gradient of error is proportional to N{sub h}{sup -1/2} which are optimal asymptotics. The methodology is verified with numerical experiments.

  18. Development of geometrically-nonlinear finite element analysis for marine risers

    E-Print Network [OSTI]

    Haas, Mark Edward

    1987-01-01T23:59:59.000Z

    an element that can work with high-aspect ratio structures. ACKNOWLEDGMENTS I wish to thank Dr. Y. K. "Jack" Lou for his patient persistence, for his strong encouragement, and for deep insight into this reses, rch. Dr. Lou stepped in as the committee... members, Dr. Mogens Henriksen, Dr. Harry A. Hogan, Dr. James R. Morgan, Dr. James K. Nelson, and Dr. John M. Niedziewicki, for their time, support, and criticisms (although not always taken in the proper context). Special thanks go to Dr. Hogan...

  19. Simplified modeling methods for mechanically fastened connections in flight structures

    E-Print Network [OSTI]

    Brewer, Brett Andrew

    2012-01-01T23:59:59.000Z

    load sharing, failure modes, and basic design are outlined such that later developments in the finite element method

  20. Damage Identification of a Composite Beam Using Finite Element Model Updating

    E-Print Network [OSTI]

    Moaveni, Babak; He, Xianfei; Conte, Joel P; De Callafon, Raymond A.

    2008-01-01T23:59:59.000Z

    Shevitz, D. W. (1996), Damage identification in structuresreview of vibration-based damage identification methods, TheM. , & Samman, M. M. (1991), Damage detection from changes

  1. NORIA-SP: A finite element computer program for analyzing liquid water transport in porous media; Yucca Mountain Site Characterization Project

    SciTech Connect (OSTI)

    Hopkins, P.L.; Eaton, R.R.; Bixler, N.E.

    1991-12-01T23:59:59.000Z

    A family of finite element computer programs has been developed at Sandia National Laboratories (SNL) most recently, NORIA-SP. The original NORIA code solves a total of four transport equations simultaneously: liquid water, water vapor, air, and energy. Consequently, use of NORIA is computer-intensive. Since many of the applications for which NORIA is used are isothermal, we decided to ``strip`` the original four-equation version, leaving only the liquid water equation. This single-phase version is NORIA-SP. The primary intent of this document is to provide the user of NORIA-SP an accurate user`s manual. Consequently, the reader should refer to the NORIA manual if additional detail is required regarding the equation development and finite element methods used. The single-equation version of the NORIA code (NORIA-SP) has been used most frequently for analyzing various hydrological scenarios for the potential underground nuclear waste repository at Yucca Mountain in western Nevada. These analyses are generally performed assuming a composite model to represent the fractured geologic media. In this model the material characteristics of the matrix and the fractures are area weighted to obtain equivalent material properties. Pressure equilibrium between the matrix and fractures is assumed so a single conservation equation can be solved. NORIA-SP is structured to accommodate the composite model. The equations for water velocities in both the rock matrix and the fractures are presented. To use the code for problems involving a single, nonfractured porous material, the user can simply set the area of the fractures to zero.

  2. Higher-Order, Space-Time Adaptive Finite Volume Methods: Algorithms, Analysis and Applications

    SciTech Connect (OSTI)

    Minion, Michael

    2014-04-29T23:59:59.000Z

    The four main goals outlined in the proposal for this project were: 1. Investigate the use of higher-order (in space and time) finite-volume methods for fluid flow problems. 2. Explore the embedding of iterative temporal methods within traditional block-structured AMR algorithms. 3. Develop parallel in time methods for ODEs and PDEs. 4. Work collaboratively with the Center for Computational Sciences and Engineering (CCSE) at Lawrence Berkeley National Lab towards incorporating new algorithms within existing DOE application codes.

  3. High-order Time-stepping for Galerkin and Collocation Methods Based on Component-wise

    E-Print Network [OSTI]

    Lambers, James

    to other spatial discretiza- tions, such as finite elements and meshless methods, extends the applicability--spectral methods, finite elements, meshless methods, Gaussian quadrature, block Lanczos I. INTRODUCTION The rapid with existing methods that are applicable to such problems, including finite element methods and meshless

  4. Superconvergence of finite element approximations by a computer-based approach

    E-Print Network [OSTI]

    Gangaraj, Srihari Kumar

    1994-01-01T23:59:59.000Z

    . Superconvcrgence poi?ts I'or I, hc & lass of "harmonic" solutious of the cquatio?s of pla?c elasl i& i ty: Cubi& i& n?dipitv square elements The superconvergence points arc le& atcd at the intersection of thc contours C a?' (P; f, T), i = 1, . . . , 4... of &he mat('ri;&I orl, l&ot ropy. Sup( r(o?v& rgcuce points fol the r& curl]po?eut of thc Ih&s, &I&('u) I'o&. tl&& ?lais of "harmonic" solutions of th&. cquatio? ol' or&i&olropi& h( &I, -&. ond?ction: Cubic serendipity ele&?cnts. /cro-co?to?rs a?d s...

  5. On simple and accurate finite element models for nonlinear bending analysis of beams and plates

    E-Print Network [OSTI]

    Urthaler Lapeira, Yetzirah Yksya

    2007-09-17T23:59:59.000Z

    ¡B ¶ ¡C3x¡C4 (2.11) where the superscript \\U00 denotes the quantities belonging to either the Timoshenko beam theory (fi = 0) or the Reddy-Bickford beam theory (fi 6= 0). The slope U has a difierent meaning for each theory, as deflned below U = 8... unifled beam element, the axial displacement ^u0(^x) and transverse de ection ^w0(^x) are interpolated as ^u0(^x) = 2X j=1 ^ujˆj(^x) ^w0(^x) = 4X j=1 ¢j`j(^x); ^u1 = ^u0( ^x1); ^u2 = ^u0( ^x2) (2.22a) ¢1 = ^w0( ^x1); ¢2 = ( ^x1); ¢3 = ^w0( ^x2); ¢4...

  6. Dynamic pulse buckling of cylindrical shells under axial impact: A benchmark study of 2D and 3D finite element calculations

    SciTech Connect (OSTI)

    Hoffman, E.L.; Ammerman, D.J.

    1995-04-01T23:59:59.000Z

    A series of tests investigating dynamic pulse buckling of a cylindrical shell under axial impact is compared to several 2D and 3D finite element simulations of the event. The purpose of the work is to investigate the performance of various analysis codes and element types on a problem which is applicable to radioactive material transport packages, and ultimately to develop a benchmark problem to qualify finite element analysis codes for the transport package design industry. During the pulse buckling tests, a buckle formed at each end of the cylinder, and one of the two buckles became unstable and collapsed. Numerical simulations of the test were performed using PRONTO, a Sandia developed transient dynamics analysis code, and ABAQUS/Explicit with both shell and continuum elements. The calculations are compared to the tests with respect to deformed shape and impact load history.

  7. Least squares based finite element formulations and their applications in fluid mechanics

    E-Print Network [OSTI]

    Prabhakar, Vivek

    2009-05-15T23:59:59.000Z

    are proposed, implemented in a computational scheme, and evaluated in this study. The continuity equation is treated as a constraint on the velocity field and the constraint is enforced using the penalty method. These penalty based formulations produce accurate...

  8. Stability Analysis for a Saline Boundary Layer Formed by Uniform Up ow Using Finite Elements

    E-Print Network [OSTI]

    Eindhoven, Technische Universiteit

    density in the deeper underground and high salt water density at the boundary layer), gravitation plays and Darcy's law. In this report we #12;rst give an overview of semi-analytical methods to analyse

  9. Power and efficiency limits for internal combustion engines via methods of finite-time thermodynamics

    E-Print Network [OSTI]

    Berry, R. Stephen

    Power and efficiency limits for internal combustion engines via methods of finite publication 17 June 1993) Analytical expressionsfor the upper bounds of power and efficiency of an internal and expensiveto compute and analyze.2If we are interestedin maximum power output or in maximum effi- ciency

  10. A stable finite difference method for the elastic wave equation on complex geometries with free surfaces

    E-Print Network [OSTI]

    Appelö, Daniel

    and periodic boundary conditions. The fully discrete version of the method conserves a discrete energy; curvilinear grids ; finite differences; stability; energy estimate; seismic wave propagation 1 Introduction, such as gas pipes, wave guides, railroad rails and disc brakes. In the vast majority of wave propagation

  11. FINITE VOLUME METHODS APPLIED TO THE COMPUTATIONAL MODELLING OF WELDING PHENOMENA

    E-Print Network [OSTI]

    Taylor, Gary

    1 FINITE VOLUME METHODS APPLIED TO THE COMPUTATIONAL MODELLING OF WELDING PHENOMENA Gareth A.Taylor@brunel.ac.uk ABSTRACT This paper presents the computational modelling of welding phenomena within a versatile numerical) and Computational Solid Mechanics (CSM). With regard to the CFD modelling of the weld pool fluid dynamics, heat

  12. Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water

    E-Print Network [OSTI]

    Kurganov, Alexander

    Finite-Volume-Particle Methods for Models of Transport of Pollutant in Shallow Water Alina Chertock of shallow water equations and the pollutant propagation is described by a transport equation. The idea and the pollution computations: the shallow water equations are numerically integrated using a #12;nite- volume

  13. A MT System from Turkmen to Turkish Employing Finite State and Statistical Methods

    E-Print Network [OSTI]

    Yanikoglu, Berrin

    between close language pairs can be relatively easier and can still benefit from simple(r) paradigms in MT with a disambiguation post-processing stage based on statistical language models. The very productive inflectionalA MT System from Turkmen to Turkish Employing Finite State and Statistical Methods A. Cüneyd TANTU

  14. A fourth-order symplectic finite-difference time-domain (FDTD) method for light scattering and a 3D Monte Carlo code for radiative transfer in scattering systems

    E-Print Network [OSTI]

    Zhai, Pengwang

    2009-06-02T23:59:59.000Z

    meter. 60 20 Geometry of a scattering event. . . . . . . . . . . . . . . . . . . . . . 63 21 An example of the atmosphere model used in the 3D Monte Carlo code for the vector radiative transfer systems. Inhomogeneous layers are divided into voxels... cases can be solved analytically. Several popular numerical methods include the T-matrix method [15, 16, 17, 18, 19], finite-element method [20, 21], finite-difference time-domain(FDTD)method[22,23,24,25,26,27,28,29,30,31,32], point-matching method [33...

  15. Finite Elements for a Beam System With Nonlinear Contact Under Periodic Excitation

    E-Print Network [OSTI]

    Boyer, Edmond

    . Rousselet Abstract Solar arrays are structures which are connected to satellites; during launch, they are in a folded position and submitted to high vibrations. In order to save mass, the flexibility of the panels investigation of nonlinear normal mode (NNM) can be a new method to describe nonlinear behaviors, this work

  16. Finite element study of multi-modal vibration damping for thermal barrier coating applications

    E-Print Network [OSTI]

    barriers for gas turbine blades, there is a growing interest in devel- oping multifunctional coatings methods is critical in many tech- nologies and especially so for gas turbine engines since many the spinning up of the turbine to avoid resonances. One of the advances made in gas turbine technologies, has

  17. Madabhushi & Haigh: The pull-out behaviour of piles. 1 FINITE ELEMENT ANALYSIS OF PILE FOUNDATIONS SUBJECT TO PULL-OUT

    E-Print Network [OSTI]

    Haigh, Stuart

    Madabhushi & Haigh: The pull-out behaviour of piles. 1 FINITE ELEMENT ANALYSIS OF PILE FOUNDATIONS.K.HAIGH, Cambridge University Engineering Department, UK ABSTRACT: The behaviour of piles in sand under tensile-out resistance with soil density, pile size and speed of pull was investigated in an attempt to predict pile

  18. Finite element simulation of laser spot welding A. De, S. K. Maiti, C. A. Walsh and H. K. D. H. Bhadeshia

    E-Print Network [OSTI]

    Cambridge, University of

    Finite element simulation of laser spot welding A. De, S. K. Maiti, C. A. Walsh and H. K. D. H laser spot welding, taking into account the temperature dependence of the physical properties and latent formation during laser spot welding, although the `double ellipsoidal' representation requires an a priori

  19. Finite element analysis of floatation of rectangular tunnels following earthquake induced liquefaction

    E-Print Network [OSTI]

    Madabhushi, S. S. C.; Madabhushi, S. P. G.

    2014-10-29T23:59:59.000Z

    of the 15th World Conference on Earthquake Engineering, Lisbon, Paper ID 3705. [18] Chan, A.H.C., (1988), A generalised fully coupled effective stress based computer procedure for problems in Geomechanics, SWANDYNE User Manual, Swansea, UK. [19... . 7, 667–691. [27] Pastor, M., Zienkiewicz, O.C. and Leung, K.H., (1985). Simple model for transient soil loading in earthquake analysis. II: Non-associative models for sands, Int. Journal of Numerical and Analytical Methods in Geomechanics, Vol.9...

  20. System and method for reproducibly mounting an optical element

    DOE Patents [OSTI]

    Eisenbies, Stephen; Haney, Steven

    2005-05-31T23:59:59.000Z

    The present invention provides a two-piece apparatus for holding and aligning the MEMS deformable mirror. The two-piece apparatus comprises a holding plate for fixedly holding an adaptive optics element in an overall optical system and a base spatially fixed with respect to the optical system and adapted for mounting and containing the holding plate. The invention further relates to a means for configuring the holding plate through adjustments to each of a number of off-set pads touching each of three orthogonal plane surfaces on the base, wherein through the adjustments the orientation of the holding plate, and the adaptive optics element attached thereto, can be aligned with respect to the optical system with six degrees of freedom when aligning the plane surface of the optical element. The mounting system thus described also enables an operator to repeatedly remove and restore the adaptive element in the optical system without the need to realign the system once that element has been aligned.