Peridynamic Multiscale Finite Element Methods
Costa, Timothy; Bond, Stephen D.; Littlewood, David John; Moore, Stan Gerald
2015-12-01
The problem of computing quantum-accurate design-scale solutions to mechanics problems is rich with applications and serves as the background to modern multiscale science research. The prob- lem can be broken into component problems comprised of communicating across adjacent scales, which when strung together create a pipeline for information to travel from quantum scales to design scales. Traditionally, this involves connections between a) quantum electronic structure calculations and molecular dynamics and between b) molecular dynamics and local partial differ- ential equation models at the design scale. The second step, b), is particularly challenging since the appropriate scales of molecular dynamic and local partial differential equation models do not overlap. The peridynamic model for continuum mechanics provides an advantage in this endeavor, as the basic equations of peridynamics are valid at a wide range of scales limiting from the classical partial differential equation models valid at the design scale to the scale of molecular dynamics. In this work we focus on the development of multiscale finite element methods for the peridynamic model, in an effort to create a mathematically consistent channel for microscale information to travel from the upper limits of the molecular dynamics scale to the design scale. In particular, we first develop a Nonlocal Multiscale Finite Element Method which solves the peridynamic model at multiple scales to include microscale information at the coarse-scale. We then consider a method that solves a fine-scale peridynamic model to build element-support basis functions for a coarse- scale local partial differential equation model, called the Mixed Locality Multiscale Finite Element Method. Given decades of research and development into finite element codes for the local partial differential equation models of continuum mechanics there is a strong desire to couple local and nonlocal models to leverage the speed and state of the
The generalized finite element method applied to the dynamic...
Office of Scientific and Technical Information (OSTI)
Title: The generalized finite element method applied to the dynamic response of heterogeneous media. Authors: Robbins, Joshua ; Voth, Thomas E. Publication Date: 2013-02-01 OSTI ...
Microscopic description of fission dynamics: finite element method...
Office of Scientific and Technical Information (OSTI)
Title: Microscopic description of fission dynamics: finite element method resolution of ... This site is a product of DOE's Office of Scientific and Technical Information (OSTI) and ...
Modular Finite Element Methods Library Version: 1.0
Energy Science and Technology Software Center (OSTI)
2010-06-22
MFEM is a general, modular library for finite element methods. It provides a variety of finite element spaces and bilinear/linear forms in 2D and 3D. MFEM also includes classes for dealing with various types of meshes and their refinement.
Kim, S.
1994-12-31
Parallel iterative procedures based on domain decomposition techniques are defined and analyzed for the numerical solution of wave propagation by finite element and finite difference methods. For finite element methods, in a Lagrangian framework, an efficient way for choosing the algorithm parameter as well as the algorithm convergence are indicated. Some heuristic arguments for finding the algorithm parameter for finite difference schemes are addressed. Numerical results are presented to indicate the effectiveness of the methods.
Parallel, adaptive finite element methods for conservation laws
Biswas, R.; Devine, K.D.; Flaherty, J.E. Rensselaer Polytechnic Institute, Troy, NY )
1994-01-01
We construct parallel finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. A posteriori estimates of spatial errors are obtained by a p-refinement technique using superconvergence at Radau points. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We compare results using different limiting schemes and demonstrate parallel efficiency through computations on an NCUBE/2 hypercube. We also present results using adaptive h- and p-refinement to reduce the computational cost of the method.
Bochev, Pavel Blagoveston
2011-06-01
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.
Code verification for the eXtended Finite Element Method (XFEM...
Office of Scientific and Technical Information (OSTI)
Title: Code verification for the eXtended Finite Element Method (XFEM): the compound cohesionless impact problem. Abstract not provided. Authors: Niederhaus, John Henry ; Voth, ...
An implementation analysis of the linear discontinuous finite element method
Becker, T. L.
2013-07-01
This paper provides an implementation analysis of the linear discontinuous finite element method (LD-FEM) that spans the space of (l, x, y, z). A practical implementation of LD includes 1) selecting a computationally efficient algorithm to solve the 4 x 4 matrix system Ax = b that describes the angular flux in a mesh element, and 2) choosing how to store the data used to construct the matrix A and the vector b to either reduce memory consumption or increase computational speed. To analyze the first of these, three algorithms were selected to solve the 4 x 4 matrix equation: Cramer's rule, a streamlined implementation of Gaussian elimination, and LAPACK's Gaussian elimination subroutine dgesv. The results indicate that Cramer's rule and the streamlined Gaussian elimination algorithm perform nearly equivalently and outperform LAPACK's implementation of Gaussian elimination by a factor of 2. To analyze the second implementation detail, three formulations of the discretized LD-FEM equations were provided for implementation in a transport solver: 1) a low-memory formulation, which relies heavily on 'on-the-fly' calculations and less on the storage of pre-computed data, 2) a high-memory formulation, which pre-computes much of the data used to construct A and b, and 3) a reduced-memory formulation, which lies between the low - and high-memory formulations. These three formulations were assessed in the Jaguar transport solver based on relative memory footprint and computational speed for increasing mesh size and quadrature order. The results indicated that the memory savings of the low-memory formulation were not sufficient to warrant its implementation. The high-memory formulation resulted in a significant speed advantage over the reduced-memory option (10-50%), but also resulted in a proportional increase in memory consumption (5-45%) for increasing quadrature order and mesh count; therefore, the practitioner should weigh the system memory constraints against any
Unconstrained paving and plastering method for generating finite element meshes
Staten, Matthew L.; Owen, Steven J.; Blacker, Teddy D.; Kerr, Robert
2010-03-02
Computer software for and a method of generating a conformal all quadrilateral or hexahedral mesh comprising selecting an object with unmeshed boundaries and performing the following while unmeshed voids are larger than twice a desired element size and unrecognizable as either a midpoint subdividable or pave-and-sweepable polyhedra: selecting a front to advance; based on sizes of fronts and angles with adjacent fronts, determining which adjacent fronts should be advanced with the selected front; advancing the fronts; detecting proximities with other nearby fronts; resolving any found proximities; forming quadrilaterals or unconstrained columns of hexahedra where two layers cross; and establishing hexahedral elements where three layers cross.
A Method for Connecting Dissimilar Finite Element Meshes in Three Dimensions
Dohrmann, C.R.; Heinstein, M.W.; Key, S.W.
1998-11-12
A method is presented for connecting dissimilar finite element meshes in three dimensions. The method combines the concept of master and slave surfaces with the uniform strain approach for surface, corrections finite elements- By modifyhg the are made to element formulations boundaries of elements on the slave such that first-order patch tests are passed. The method can be used to connect meshes which use different element types. In addition, master and slave surfaces can be designated independently of relative mesh resolutions. Example problems in three-dimensional linear elasticity are presented.
Domain decomposition based iterative methods for nonlinear elliptic finite element problems
Cai, X.C.
1994-12-31
The class of overlapping Schwarz algorithms has been extensively studied for linear elliptic finite element problems. In this presentation, the author considers the solution of systems of nonlinear algebraic equations arising from the finite element discretization of some nonlinear elliptic equations. Several overlapping Schwarz algorithms, including the additive and multiplicative versions, with inexact Newton acceleration will be discussed. The author shows that the convergence rate of the Newton`s method is independent of the mesh size used in the finite element discretization, and also independent of the number of subdomains into which the original domain in decomposed. Numerical examples will be presented.
Taylor, G.; Dong, C.; Sun, S.
2010-03-18
A mathematical model for contaminant species passing through fractured porous media is presented. In the numerical model, we combine two locally conservative methods, i.e. mixed finite element (MFE) and the finite volume methods. Adaptive triangle mesh is used for effective treatment of the fractures. A hybrid MFE method is employed to provide an accurate approximation of velocities field for both the fractures and matrix which are crucial to the convection part of the transport equation. The finite volume method and the standard MFE method are used to approximate the convection and dispersion terms respectively. The model is used to investigate the interaction of adsorption with transport and to extract information on effective adsorption distribution coefficients. Numerical examples in different fractured media illustrate the robustness and efficiency of the proposed numerical model.
A massively parallel adaptive finite element method with dynamic load balancing
Devine, K.D.; Flaherty, J.E.; Wheat, S.R.; Maccabe, A.B.
1993-05-01
We construct massively parallel, adaptive finite element methods for the solution of hyperbolic conservation laws in one and two dimensions. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. The resulting method is of high order and may be parallelized efficiently on MIMD computers. We demonstrate parallel efficiency through computations on a 1024-processor nCUBE/2 hypercube. We also present results using adaptive p-refinement to reduce the computational cost of the method. We describe tiling, a dynamic, element-based data migration system. Tiling dynamically maintains global load balance in the adaptive method by overlapping neighborhoods of processors, where each neighborhood performs local load balancing. We demonstrate the effectiveness of the dynamic load balancing with adaptive p-refinement examples.
A massively parallel adaptive finite element method with dynamic load balancing
Devine, K.D.; Flaherty, J.E.; Wheat, S.R.; Maccabe, A.B.
1993-12-31
The authors construct massively parallel adaptive finite element methods for the solution of hyperbolic conservation laws. Spatial discretization is performed by a discontinuous Galerkin finite element method using a basis of piecewise Legendre polynomials. Temporal discretization utilizes a Runge-Kutta method. Dissipative fluxes and projection limiting prevent oscillations near solution discontinuities. The resulting method is of high order and may be parallelized efficiently on MIMD computers. They demonstrate parallel efficiency through computations on a 1024-processor nCUBE/2 hypercube. They present results using adaptive p-refinement to reduce the computational cost of the method, and tiling, a dynamic, element-based data migration system that maintains global load balance of the adaptive method by overlapping neighborhoods of processors that each perform local balancing.
Analysis of Mode III Elastodynamic Cracked Plane using the Fractal Two-Level Finite Element Method
Fan, J.; Lee, Y. Y.; Leung, A. Y. T.
2010-05-21
In this study, the fractal two-level finite element method, which has mainly been used for static cracked plane problems, is applied to the cracked plane problem. Using the transformation process in the proposed method, the infinite dimension of the finite element matrices that are assembled for a singular region is made finite in terms of the dynamics stress intensity factors directly, and thus the computational time can be reduced significantly. The Newmark time integration scheme is then used to obtain the dynamic stress intensity factors. The results from the proposed method are in reasonable agreement with those of classical methods. The main drawback of the time integration scheme is that numerical oscillations are induced in some cases.
Parameters of lossy cavity resonators calculated by the finite element method
Groiss, S.; Bardi, I.; Biro, O.; Preis, K.; Richter, K.R.
1996-05-01
The calculation of the resonance frequency and quality factor of closed or aperture coupled cavity resonators with volume and wall losses by an edge finite element method is discussed. An efficient solver is developed to solve the complex nonlinear eigenvalue problem. The effect of the roughness of the walls on the quality factor is taken approximately into account.
SQA of finite element method (FEM) codes used for analyses of pit storage/transport packages
Russel, E.
1997-11-01
This report contains viewgraphs on the software quality assurance of finite element method codes used for analyses of pit storage and transport projects. This methodology utilizes the ISO 9000-3: Guideline for application of 9001 to the development, supply, and maintenance of software, for establishing well-defined software engineering processes to consistently maintain high quality management approaches.
Mixed-RKDG Finite Element Methods for the 2-D Hydrodynamic Model for Semiconductor Device Simulation
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Chen, Zhangxin; Cockburn, Bernardo; Jerome, Joseph W.; Shu, Chi-Wang
1995-01-01
In this paper we introduce a new method for numerically solving the equations of the hydrodynamic model for semiconductor devices in two space dimensions. The method combines a standard mixed finite element method, used to obtain directly an approximation to the electric field, with the so-called Runge-Kutta Discontinuous Galerkin (RKDG) method, originally devised for numerically solving multi-dimensional hyperbolic systems of conservation laws, which is applied here to the convective part of the equations. Numerical simulations showing the performance of the new method are displayed, and the results compared with those obtained by using Essentially Nonoscillatory (ENO) finite difference schemes. FrommoreÂ Â» the perspective of device modeling, these methods are robust, since they are capable of encompassing broad parameter ranges, including those for which shock formation is possible. The simulations presented here are for Gallium Arsenide at room temperature, but we have tested them much more generally with considerable success.Â«Â less
A 3D finite element ALE method using an approximate Riemann solution
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Chiravalle, V. P.; Morgan, N. R.
2016-08-09
Arbitrary Lagrangianâ€“Eulerian finite volume methods that solve a multidimensional Riemann-like problem at the cell center in a staggered grid hydrodynamic (SGH) arrangement have been proposed. This research proposes a new 3D finite element arbitrary Lagrangianâ€“Eulerian SGH method that incorporates a multidimensional Riemann-like problem. Here, two different Riemann jump relations are investigated. A new limiting method that greatly improves the accuracy of the SGH method on isentropic flows is investigated. A remap method that improves upon a well-known mesh relaxation and remapping technique in order to ensure total energy conservation during the remap is also presented. Numerical details and test problemmoreÂ Â» results are presented.Â«Â less
Partition-of-unity finite-element method for large scale quantum...
Office of Scientific and Technical Information (OSTI)
process using modern partition-of-unity (PU) techniques in finite element analysis. ... Country of Publication: United States Language: English Subject: 75 CONDENSED MATTER ...
Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)
Dolbow, John; Zhang, Ziyu; Spencer, Benjamin; Jiang, Wen
2015-09-01
Efforts are underway to develop fracture mechanics capabilities in the Grizzly code to enable it to be used to perform deterministic fracture assessments of degraded reactor pressure vessels (RPVs). A capability was previously developed to calculate three-dimensional interaction- integrals to extract mixed-mode stress-intensity factors. This capability requires the use of a finite element mesh that conforms to the crack geometry. The eXtended Finite Element Method (X-FEM) provides a means to represent a crack geometry without explicitly fitting the finite element mesh to it. This is effected by enhancing the element kinematics to represent jump discontinuities at arbitrary locations inside of the element, as well as the incorporation of asymptotic near-tip fields to better capture crack singularities. In this work, use of only the discontinuous enrichment functions was examined to see how accurate stress intensity factors could still be calculated. This report documents the following work to enhance Grizzlyâ€™s engineering fracture capa- bilities by introducing arbitrary jump discontinuities for prescribed crack geometries; X-FEM Mesh Cutting in 3D: to enhance the kinematics of elements that are intersected by arbitrary crack geometries, a mesh cutting algorithm was implemented in Grizzly. The algorithm introduces new virtual nodes and creates partial elements, and then creates a new mesh connectivity; Interaction Integral Modifications: the existing code for evaluating the interaction integral in Grizzly was based on the assumption of a mesh that was fitted to the crack geometry. Modifications were made to allow for the possibility of a crack front that passes arbitrarily through the mesh; and Benchmarking for 3D Fracture: the new capabilities were benchmarked against mixed- mode three-dimensional fracture problems with known analytical solutions.
Permeability computation on a REV with an immersed finite element method
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-04
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.
An h-adaptive finite element method for turbulent heat transfer
Carriington, David B [Los Alamos National Laboratory
2009-01-01
A two-equation turbulence closure model (k-{omega}) using an h-adaptive grid technique and finite element method (FEM) has been developed to simulate low Mach flow and heat transfer. These flows are applicable to many flows in engineering and environmental sciences. Of particular interest in the engineering modeling areas are: combustion, solidification, and heat exchanger design. Flows for indoor air quality modeling and atmospheric pollution transport are typical types of environmental flows modeled with this method. The numerical method is based on a hybrid finite element model using an equal-order projection process. The model includes thermal and species transport, localized mesh refinement (h-adaptive) and Petrov-Galerkin weighting for the stabilizing the advection. This work develops the continuum model of a two-equation turbulence closure method. The fractional step solution method is stated along with the h-adaptive grid method (Carrington and Pepper, 2002). Solutions are presented for 2d flow over a backward-facing step.
Koteras, J.R.
1993-07-01
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.
Finite element modeling of lightning
Hwang, C.C.; Huang, S.R.; Bor, S.S.
1995-12-31
In this paper the transmission line equation which describes the transient voltage and current distributions of a lightning stroke is employed. Finite element method is used to derive the element equations and one-dimensional linear elements are used to discretize the field region. The implicit Newmark time integration technique is used to convert the resulting second-order ordinary differential equations into a set of recurrence equations which are then solved at each time step. Numerical example is included and discussed.
Numerical simulation and design of a fluxset sensor by finite element method
Preis, K.; Bardi, I.; Biro, O.; Richter, K.R.; Pavo, J.; Gasparics, A.; Ticar, I.
1998-09-01
A 3D model of a fluxset sensor serving to measure magnetic fields arising in Eddy Current Nondestructive Testing applications is analyzed by the finite element method. The voltage induced in the pick-up coil is obtained by computing the flux of the core of the sensor for several values of the exciting current at various external fields. It is shown that the time shift of the ensuing voltage impulse depends linearly on the external field in a wide range. The behavior of the sensor is furthermore simulated in a real nondestructive testing arrangement consisting of an exciting coil located above a conducting plate with a crack.
Method and apparatus for connecting finite element meshes and performing simulations therewith
Dohrmann, Clark R.; Key, Samuel W.; Heinstein, Martin W.
2003-05-06
The present invention provides a method of connecting dissimilar finite element meshes. A first mesh, designated the master mesh, and a second mesh, designated the slave mesh, each have interface surfaces proximal the other. Each interface surface has a corresponding interface mesh comprising a plurality of interface nodes. Each slave interface node is assigned new coordinates locating the interface node on the interface surface of the master mesh. The slave interface surface is further redefined to be the projection of the slave interface mesh onto the master interface surface.
Energy Science and Technology Software Center (OSTI)
2006-03-08
MAPVAR-KD is designed to transfer solution results from one finite element mesh to another. MAPVAR-KD draws heavily from the structure and coding of MERLIN II, but it employs a new finite element data base, EXODUS II, and offers enhanced speed and new capabilities not available in MERLIN II. In keeping with the MERLIN II documentation, the computational algorithms used in MAPVAR-KD are described. User instructions are presented. Example problems are included to demonstrate the operationmoreÂ Â» of the code and the effects of various input options. MAPVAR-KD is a modification of MAPVAR in which the search algorithm was replaced by a kd-tree-based search for better performance on large problems.Â«Â less
Energy Science and Technology Software Center (OSTI)
2005-06-26
Exotxt is an analysis code that reads finite element results data stored in an exodusII file and generates a file in a structured text format. The text file can be edited or modified via a number of text formatting tools. Exotxt is used by analysis to translate data from the binary exodusII format into a structured text format which can then be edited or modified and then either translated back to exodusII format or tomoreÂ Â» another format.Â«Â less
Energy Science and Technology Software Center (OSTI)
2005-05-07
CONEX is a code for joining sequentially in time multiple exodusll database files which all represent the same base mesh topology and geometry. It is used to create a single results or restart file from multiple results or restart files which typically arise as the result of multiple restarted analyses. CONEX is used to postprocess the results from a series of finite element analyses. It can join sequentially the data from multiple results databases intomoreÂ Â» a single database which makes it easier to postprocess the results data.Â«Â less
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-04-14
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale mediummoreÂ Â» property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.Â«Â less
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-08-15
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
Gao, Kai; Fu, Shubin; Gibson, Richard L.; Chung, Eric T.; Efendiev, Yalchin
2015-04-14
It is important to develop fast yet accurate numerical methods for seismic wave propagation to characterize complex geological structures and oil and gas reservoirs. However, the computational cost of conventional numerical modeling methods, such as finite-difference method and finite-element method, becomes prohibitively expensive when applied to very large models. We propose a Generalized Multiscale Finite-Element Method (GMsFEM) for elastic wave propagation in heterogeneous, anisotropic media, where we construct basis functions from multiple local problems for both the boundaries and interior of a coarse node support or coarse element. The application of multiscale basis functions can capture the fine scale medium property variations, and allows us to greatly reduce the degrees of freedom that are required to implement the modeling compared with conventional finite-element method for wave equation, while restricting the error to low values. We formulate the continuous Galerkin and discontinuous Galerkin formulation of the multiscale method, both of which have pros and cons. Applications of the multiscale method to three heterogeneous models show that our multiscale method can effectively model the elastic wave propagation in anisotropic media with a significant reduction in the degrees of freedom in the modeling system.
Lee, D. W.; Joo, H. G.
2013-07-01
The methods and performance of a three-dimensional S{sub n} transport code employing the Discontinuous Finite Element Method (DFEM) and the Coarse Mesh Finite Difference (CMFD) formulation are presented. The mesh generator GMSH and a post processing visualization tool Visit are combined with the code for flexible geometry processing and versatile visualization. The CMFD method for DFEM Sn applications is formulated and the performance of the CMFD acceleration of eigenvalue calculations is demonstrated for a simple set of neutron transport problems. (authors)
Full wave simulation of waves in ECRIS plasmas based on the finite element method
Torrisi, G.; Mascali, D.; Neri, L.; Castro, G.; Patti, G.; Celona, L.; Gammino, S.; Ciavola, G.; Di Donato, L.; Sorbello, G.; Isernia, T.
2014-02-12
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.
Evaluation of Test Methods for Triaxially Braided Composites using a Meso-Scale Finite Element Model
Zhang, Chao
2015-10-01
The characterization of triaxially braided composite is complicate due to the nonuniformity of deformation within the unit cell as well as the possibility of the freeedge effect related to the large size of the unit cell. Extensive experimental investigation has been conducted to develop more accurate test approaches in characterizing the actual mechanical properties of the material we are studying. In this work, a meso-scale finite element model is utilized to simulate two complex specimens: notched tensile specimen and tube tensile specimen, which are designed to avoid the free-edge effect and free-edge effect induced premature edge damage. The full field strain data is predicted numerically and compared with experimental data obtained by Digit Image Correlation. The numerically predicted tensile strength values are compared with experimentally measured results. The discrepancy between numerically predicted and experimentally measured data, the capability of different test approaches are analyzed and discussed. The presented numerical model could serve as assistance to the evaluation of different test methods, and is especially useful in identifying potential local damage events.
Solution of elastic-plastic stress analysis problems by the P-version of the finite element method
Szabo, B.A.; Holzer, S.M.; Actis, R.L.
1995-12-31
The solution of small-strain elastic-plastic stress analysis problems by the p-version of the finite element method is discussed. The formulation is based on the deformation theory of plasticity and the displacement method. Practical realization of controlling discretization errors for elastic-plastic problems is the main focus of the paper. Numerical examples, which include comparisons between the deformation and incremental theories of plasticity under tight control of discretization errors, are presented.
X. Frank Xu
2010-03-30
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
Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Prudencio, E.; Schussman, G.; Uplenchwar, R.; Ko, K.; /SLAC
2009-06-19
Over the past years, SLAC's Advanced Computations Department (ACD), under SciDAC sponsorship, has developed a suite of 3D (2D) parallel higher-order finite element (FE) codes, T3P (T2P) and Pic3P (Pic2P), aimed at accurate, large-scale simulation of wakefields and particle-field interactions in radio-frequency (RF) cavities of complex shape. The codes are built on the FE infrastructure that supports SLAC's frequency domain codes, Omega3P and S3P, to utilize conformal tetrahedral (triangular)meshes, higher-order basis functions and quadratic geometry approximation. For time integration, they adopt an unconditionally stable implicit scheme. Pic3P (Pic2P) extends T3P (T2P) to treat charged-particle dynamics self-consistently using the PIC (particle-in-cell) approach, the first such implementation on a conformal, unstructured grid using Whitney basis functions. Examples from applications to the International Linear Collider (ILC), Positron Electron Project-II (PEP-II), Linac Coherent Light Source (LCLS) and other accelerators will be presented to compare the accuracy and computational efficiency of these codes versus their counterparts using structured grids.
Regnier, D.; Verriere, M.; Dubray, N.; Schunck, N.
2015-11-30
In this study, we describe the software package FELIX that solves the equations of the time-dependent generator coordinate method (TDGCM) in NN-dimensions (N â‰¥ 1) under the Gaussian overlap approximation. The numerical resolution is based on the Galerkin finite element discretization of the collective space and the Crankâ€“Nicolson scheme for time integration. The TDGCM solver is implemented entirely in C++. Several additional tools written in C++, Python or bash scripting language are also included for convenience. In this paper, the solver is tested with a series of benchmarks calculations. We also demonstrate the ability of our code to handle a realistic calculation of fission dynamics.
Blacker, Teddy D.
1994-01-01
An automatic quadrilateral surface discretization method and apparatus is provided for automatically discretizing a geometric region without decomposing the region. The automated quadrilateral surface discretization method and apparatus automatically generates a mesh of all quadrilateral elements which is particularly useful in finite element analysis. The generated mesh of all quadrilateral elements is boundary sensitive, orientation insensitive and has few irregular nodes on the boundary. A permanent boundary of the geometric region is input and rows are iteratively layered toward the interior of the geometric region. Also, an exterior permanent boundary and an interior permanent boundary for a geometric region may be input and the rows are iteratively layered inward from the exterior boundary in a first counter clockwise direction while the rows are iteratively layered from the interior permanent boundary toward the exterior of the region in a second clockwise direction. As a result, a high quality mesh for an arbitrary geometry may be generated with a technique that is robust and fast for complex geometric regions and extreme mesh gradations.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Regnier, D.; Verriere, M.; Dubray, N.; Schunck, N.
2015-11-30
In this study, we describe the software package FELIX that solves the equations of the time-dependent generator coordinate method (TDGCM) in NN-dimensions (N â‰¥ 1) under the Gaussian overlap approximation. The numerical resolution is based on the Galerkin finite element discretization of the collective space and the Crankâ€“Nicolson scheme for time integration. The TDGCM solver is implemented entirely in C++. Several additional tools written in C++, Python or bash scripting language are also included for convenience. In this paper, the solver is tested with a series of benchmarks calculations. We also demonstrate the ability of our code to handle amoreÂ Â» realistic calculation of fission dynamics.Â«Â less
Tao, Liang; McCurdy, C.W.; Rescigno, T.N.
2008-11-25
We show how to combine finite elements and the discrete variable representation in prolate spheroidal coordinates to develop a grid-based approach for quantum mechanical studies involving diatomic molecular targets. Prolate spheroidal coordinates are a natural choice for diatomic systems and have been used previously in a variety of bound-state applications. The use of exterior complex scaling in the present implementation allows for a transparently simple way of enforcing Coulomb boundary conditions and therefore straightforward application to electronic continuum problems. Illustrative examples involving the bound and continuum states of H2+, as well as the calculation of photoionization cross sections, show that the speed and accuracy of the present approach offer distinct advantages over methods based on single-center expansions.
Partition-of-unity finite-element method for large scale quantum...
Office of Scientific and Technical Information (OSTI)
... OF FREEDOM; MOLECULAR DYNAMICS METHOD; OPTIMIZATION; PARALLEL PROCESSING; QUANTUM MECHANICS Word Cloud More Like This Full Text preview image File size NAView Full Text ...
Nguyen, Dang Van; Li, Jing-Rebecca; Grebenkov, Denis; Le Bihan, Denis
2014-04-15
The complex transverse water proton magnetization subject to diffusion-encoding magnetic field gradient pulses in a heterogeneous medium can be modeled by the multiple compartment Blochâ€“Torrey partial differential equation (PDE). In addition, steady-state Laplace PDEs can be formulated to produce the homogenized diffusion tensor that describes the diffusion characteristics of the medium in the long time limit. In spatial domains that model biological tissues at the cellular level, these two types of PDEs have to be completed with permeability conditions on the cellular interfaces. To solve these PDEs, we implemented a finite elements method that allows jumps in the solution at the cell interfaces by using double nodes. Using a transformation of the Blochâ€“Torrey PDE we reduced oscillations in the searched-for solution and simplified the implementation of the boundary conditions. The spatial discretization was then coupled to the adaptive explicit Rungeâ€“Kuttaâ€“Chebyshev time-stepping method. Our proposed method is second order accurate in space and second order accurate in time. We implemented this method on the FEniCS C++ platform and show time and spatial convergence results. Finally, this method is applied to study some relevant questions in diffusion MRI.
2-d Finite Element Code Postprocessor
Energy Science and Technology Software Center (OSTI)
1996-07-15
ORION is an interactive program that serves as a postprocessor for the analysis programs NIKE2D, DYNA2D, TOPAZ2D, and CHEMICAL TOPAZ2D. ORION reads binary plot files generated by the two-dimensional finite element codes currently used by the Methods Development Group at LLNL. Contour and color fringe plots of a large number of quantities may be displayed on meshes consisting of triangular and quadrilateral elements. ORION can compute strain measures, interface pressures along slide lines, reaction forcesmoreÂ Â» along constrained boundaries, and momentum. ORION has been applied to study the response of two-dimensional solids and structures undergoing finite deformations under a wide variety of large deformation transient dynamic and static problems and heat transfer analyses.Â«Â less
3-D Finite Element Code Postprocessor
Energy Science and Technology Software Center (OSTI)
1996-07-15
TAURUS is an interactive post-processing application supporting visualization of finite element analysis results on unstructured grids. TAURUS provides the ability to display deformed geometries and contours or fringes of a large number of derived results on meshes consisting of beam, plate, shell, and solid type finite elements. Time history plotting is also available.
Williams, P.T.
1993-09-01
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.
Finite element model for heat conduction in jointed rock masses
Gartling, D.K.; Thomas, R.K.
1981-01-01
A computatonal procedure for simulating heat conduction in a fractured rock mass is proposed and illustrated in the present paper. The method makes use of a simple local model for conduction in the vicinity of a single open fracture. The distributions of fractures and fracture properties within the finite element model are based on a statistical representation of geologic field data. Fracture behavior is included in the finite element computation by locating local, discrete fractures at the element integration points.
Finite element analysis of multilayer coextrusion.
Hopkins, Matthew Morgan; Schunk, Peter Randall; Baer, Thomas A.; Mrozek, Randy A.; Lenhart, Joseph Ludlow; Rao, Rekha Ranjana; Collins, Robert; Mondy, Lisa Ann
2011-09-01
Multilayer coextrusion has become a popular commercial process for producing complex polymeric products from soda bottles to reflective coatings. A numerical model of a multilayer coextrusion process is developed based on a finite element discretization and two different free-surface methods, an arbitrary-Lagrangian-Eulerian (ALE) moving mesh implementation and an Eulerian level set method, to understand the moving boundary problem associated with the polymer-polymer interface. The goal of this work is to have a numerical capability suitable for optimizing and troubleshooting the coextrusion process, circumventing flow instabilities such as ribbing and barring, and reducing variability in layer thickness. Though these instabilities can be both viscous and elastic in nature, for this work a generalized Newtonian description of the fluid is used. Models of varying degrees of complexity are investigated including stability analysis and direct three-dimensional finite element free surface approaches. The results of this work show how critical modeling can be to reduce build test cycles, improve material choices, and guide mold design.
Hugger, J.
1995-12-31
When the finite element solution of a variational problem possesses certain super convergence properties, it is possible very inexpensively to obtain a correction term providing an additional order of approximation of the solution. The correction can be used for error estimation locally or globally in whatever norm is preferred, or if no error estimation is wanted it can be used for postprocessing of the solution to improve the quality. In this paper such a correction term is described for the general case of n dimensional, linear or nonlinear problems. Computational evidence of the performance in one space dimension is given with special attention to the effects of the appearance of singularities and zeros of derivatives in the exact solution.
Quadrilateral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Benzley, Steven E
2012-10-16
Techniques for coarsening a quadrilateral mesh are described. These techniques include identifying a coarsening region within the quadrilateral mesh to be coarsened. Quadrilateral elements along a path through the coarsening region are removed. Node pairs along opposite sides of the path are identified. The node pairs along the path are then merged to collapse the path.
Application of Mass Lumped Higher Order Finite Elements
Chen, J.; Strauss, H. R.; Jardin, S. C.; Park, W.; Sugiyama, L. E.; G. Fu; Breslau, J.
2005-11-01
There are many interesting phenomena in extended-MHD such as anisotropic transport, mhd, 2-fluid effects stellarator and hot particles. Any one of them challenges numerical analysts, and researchers are seeking for higher order methods, such as higher order finite difference, higher order finite elements and hp/spectral elements. It is true that these methods give more accurate solution than their linear counterparts. However, numerically they are prohibitively expensive. Here we give a successful solution of this conflict by applying mass lumped higher order finite elements. This type of elements not only keep second/third order accuracy but also scale closely to linear elements by doing mass lumping. This is especially true for second order lump elements. Full M3D and anisotropic transport models are studied.
Cheng, C.Z.
1988-12-01
A nonvariational ideal MHD stability code (NOVA) has been developed. In a general flux coordinate (/psi/, theta, /zeta/) system with an arbitrary Jacobian, the NOVA code employs Fourier expansions in the generalized poloidal angle theta and generalized toroidal angle /zeta/ directions, and cubic-B spline finite elements in the radial /psi/ direction. Extensive comparisons with these variational ideal MHD codes show that the NOVA code converges faster and gives more accurate results. An extended version of NOVA is developed to integrate non-Hermitian eigenmode equations due to energetic particles. The set of non-Hermitian integro-differential eigenmode equations is numerically solved by the NOVA-K code. We have studied the problems of the stabilization of ideal MHD internal kink modes by hot particle pressure and the excitation of ''fishbone'' internal kink modes by resonating with the energetic particle magnetic drift frequency. Comparisons with analytical solutions show that the values of the critical ..beta../sub h/ from the analytical theory can be an order of magnitude different from those computed by the NOVA-K code. 24 refs., 11 figs., 1 tab.
Quadrilateral/hexahedral finite element mesh coarsening
Staten, Matthew L; Dewey, Mark W; Scott, Michael A; Benzley, Steven E
2012-10-16
A technique for coarsening a finite element mesh ("FEM") is described. This technique includes identifying a coarsening region within the FEM to be coarsened. Perimeter chords running along perimeter boundaries of the coarsening region are identified. The perimeter chords are redirected to create an adaptive chord separating the coarsening region from a remainder of the FEM. The adaptive chord runs through mesh elements residing along the perimeter boundaries of the coarsening region. The adaptive chord is then extracted to coarsen the FEM.
Finite Element Interface to Linear Solvers
Energy Science and Technology Software Center (OSTI)
2005-03-18
Sparse systems of linear equations arise in many engineering applications, including finite elements, finite volumes, and others. The solution of linear systems is often the most computationally intensive portion of the application. Depending on the complexity of problems addressed by the application, there may be no single solver capable of solving all of the linear systems that arise. This motivates the desire to switch an application from one solver librwy to another, depending on themoreÂ Â» problem being solved. The interfaces provided by solver libraries differ greatly, making it difficult to switch an application code from one library to another. The amount of library-specific code in an application Can be greatly reduced by having an abstraction layer between solver libraries and the application, putting a common "face" on various solver libraries. One such abstraction layer is the Finite Element Interface to Linear Solvers (EEl), which has seen significant use by finite element applications at Sandia National Laboratories and Lawrence Livermore National Laboratory.Â«Â less
Finite Element Heat & Mass Transfer Code
Energy Science and Technology Software Center (OSTI)
1996-10-10
FEHM is a numerical simulation code for subsurface transport processes. It models 3-D, time-dependent, multiphase, multicomponent, non-isothermal, reactive flow through porous and fractured media. It can accurately represent complex 3-D geologic media and structures and their effects on subsurface flow and transport. Its capabilities include flow of gas, water, and heat; flow of air, water, and heat; multiple chemically reactive and sorbing tracers; finite element/finite volume formulation; coupled stress module; saturated and unsaturated media; andmore »double porosity and double porosity/double permeability capabilities.« less
Finite Element Heat & Mass Transfer Code
Energy Science and Technology Software Center (OSTI)
1996-10-10
FEHM is a numerical simulation code for subsurface transport processes. It models 3-D, time-dependent, multiphase, multicomponent, non-isothermal, reactive flow through porous and fractured media. It can accurately represent complex 3-D geologic media and structures and their effects on subsurface flow and transport. Its capabilities include flow of gas, water, and heat; flow of air, water, and heat; multiple chemically reactive and sorbing tracers; finite element/finite volume formulation; coupled stress module; saturated and unsaturated media; andmoreÂ Â» double porosity and double porosity/double permeability capabilities.Â«Â less
DYNA3D: a finite element program for supercomputers
Benson, D.J.; Hallquist, J.O.
1986-01-01
DYNA3D is an explicit three-dimensional finite element code for analyzing the large deformation dynamic response of inelastic solids and structures. A contact-impact algorithm permits gaps and sliding along material interfaces with friction. By a specialization of this algorithm, such interfaces can be rigidly tied to admit variable zoning without the need of transition regions. Spatial discretization is achieved by the use of 8-node solid elements, 2-node beam elements, 4-node shell elements, 8-node solid shell elements, and rigid bodies. The equations-of-motion are integrated in time by the central difference method. 8 refs., 15 figs.
Finite Element Results Visualization for Unstructured Grids
Energy Science and Technology Software Center (OSTI)
1996-07-15
GRIZ is a general-purpose post-processing application supporting interactive visualization of finite element analysis results on unstructured grids. In addition to basic pseudocolor renderings of state variables over the mesh surface, GRIZ provides modern visualization techniques such as isocontours and isosurfaces, cutting planes, vector field display, and particle traces. GRIZ accepts both command-line and mouse-driven input, and is portable to virtually any UNIX platform which provides Motif and OpenGl libraries.
ExodusII Finite Element Data Model
Energy Science and Technology Software Center (OSTI)
2005-05-14
EXODUS II is a model developed to store and retrieve data for finite element analyses. It is used for preprocessing (problem definition), postprocessing (results visualization), as well as code to code data transfer. An EXODUS II data file is a random access, machine independent, binary file that is written and read via C, C++, or Fortran library routines which comprise the Application Programming Interface. (exodus II is based on netcdf)
TAURUS. 3-d Finite Element Code Postprocessor
Whirley, R.G.
1991-05-01
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D (ESTSC 139), DYNA3D (ESTSC 138), TACO3D (ESTSC 287), TOPAZ3D (ESTSC 231), and GEMINI (ESTSC 455) and plots contours, time histories,and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
TAURUS. 3-d Finite Element Code Postprocessor
Whirley, R.G.
1992-03-03
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D (ESTSC 139), DYNA3D (ESTSC 138), TACO3D (ESTSC 287), TOPAZ3D (ESTSC 231), and GEMINI (ESTSC 455) and plots contours, time histories,and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
TAURUS. 3-D Finite Element Code Postprocessor
Whirley, R.G.
1984-05-01
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D, DYNA3D, TACO3D, TOPAZ3D, and GEMINI and plots contours, time histories,and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
TAURUS. 3-D Finite Element Code Postprocessor
Whirley, R.G.
1993-11-30
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D, DYNA3D, TACO3D, TOPAZ3D, and GEMINI and plots contours, time histories,and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
TAURUS. 3-D Finite Element Code Postprocessor
Whirley, R.G.
1992-03-03
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D, DYNA3D, TACO3D, TOPAZ3D, and GEMINI and plots contours, time histories,and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
TAURUS. 3-D Finite Element Code Postprocessor
Kennedy, T.
1992-03-03
TAURUS reads the binary plot files generated by the LLNL three-dimensional finite element analysis codes, NIKE3D, DYNA3D, TACO3D, TOPAZ3D, and GEMINI and plots contours, time histories, and deformed shapes. Contours of a large number of quantities may be plotted on meshes consisting of plate, shell, and solid type elements. TAURUS can compute a variety of strain measures, reaction forces along constrained boundaries, and momentum. TAURUS has three phases: initialization, geometry display with contouring, and time history processing.
Modal Substructuring of Geometrically Nonlinear Finite-Element Models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Kuether, Robert J.; Allen, Matthew S.; Hollkamp, Joseph J.
2016-02-01
The efficiency of a modal substructuring method depends on the component modes used to reduce each subcomponent model. Methods such as Craigâ€“Bampton have been used extensively to reduce linear finite-element models with thousands or even millions of degrees of freedom down orders of magnitude while maintaining acceptable accuracy. A novel reduction method is proposed here for geometrically nonlinear finite-element models using the fixed-interface and constraint modes of the linearized system to reduce each subcomponent model. The geometric nonlinearity requires an additional cubic and quadratic polynomial function in the modal equations, and the nonlinear stiffness coefficients are determined by applying amoreÂ Â» series of static loads and using the finite-element code to compute the response. The geometrically nonlinear, reduced modal equations for each subcomponent are then coupled by satisfying compatibility and force equilibrium. This modal substructuring approach is an extension of the Craigâ€“Bampton method and is readily applied to geometrically nonlinear models built directly within commercial finite-element packages. The efficiency of this new approach is demonstrated on two example problems: one that couples two geometrically nonlinear beams at a shared rotational degree of freedom, and another that couples an axial spring element to the axial degree of freedom of a geometrically nonlinear beam. The nonlinear normal modes of the assembled models are compared with those of a truth model to assess the accuracy of the novel modal substructuring approach.Â«Â less
The Uranium Processing Facility (UPF) Finite Element Meshing Discussion
Office of Energy Efficiency and Renewable Energy (EERE)
The Uranium Processing Facility (UPF) Finite Element Meshing Discussion Loring Wyllie Arne Halterman Degenkolb Engineers, San Francisco
2-D Finite Element Heat Conduction
Energy Science and Technology Software Center (OSTI)
1989-10-30
AYER is a finite element program which implicitly solves the general two-dimensional equation of thermal conduction for plane or axisymmetric bodies. AYER takes into account the effects of time (transient problems), in-plane anisotropic thermal conductivity, a three-dimensional velocity distribution, and interface thermal contact resistance. Geometry and material distributions are arbitrary, and input is via subroutines provided by the user. As a result, boundary conditions, material properties, velocity distributions, and internal power generation may be mademoreÂ Â» functions of, e.g., time, temperature, location, and heat flux.Â«Â less
Guerin, P.; Baudron, A. M.; Lautard, J. J.
2006-07-01
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)
3-D Finite Element Heat Transfer
Energy Science and Technology Software Center (OSTI)
1992-02-01
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 functionalmoreÂ Â» 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.Â«Â less
HYDRA, A finite element computational fluid dynamics code: User manual
Christon, M.A.
1995-06-01
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.
Rivard, MJ; Ghadyani, HR; Bastien, AD; Lutz, NN; Hepel, JT
2015-06-15
Purpose: Noninvasive image-guided breast brachytherapy delivers conformal HDR Ir-192 brachytherapy treatments with the breast compressed, and treated in the cranial-caudal and medial-lateral directions. This technique subjects breast tissue to extreme deformations not observed for other disease sites. Given that, commercially-available software for deformable image registration cannot accurately co-register image sets obtained in these two states, a finite element analysis based on a biomechanical model was developed to deform dose distributions for each compression circumstance for dose summation. Methods: The model assumed the breast was under planar stress with values of 30 kPa for Youngâ€™s modulus and 0.3 for Poissonâ€™s ratio. Dose distributions from round and skin-dose optimized applicators in cranial-caudal and medial-lateral compressions were deformed using 0.1 cm planar resolution. Dose distributions, skin doses, and dose-volume histograms were generated. Results were examined as a function of breast thickness, applicator size, target size, and offset distance from the center. Results: Over the range of examined thicknesses, target size increased several millimeters as compression thickness decreased. This trend increased with increasing offset distances. Applicator size minimally affected target coverage, until applicator size was less than the compressed target size. In all cases, with an applicator larger or equal to the compressed target size, > 90% of the target covered by > 90% of the prescription dose. In all cases, dose coverage became less uniform as offset distance increased and average dose increased. This effect was more pronounced for smaller target-applicator combinations. Conclusions: The model exhibited skin dose trends that matched MC-generated benchmarking results and clinical measurements within 2% over a similar range of breast thicknesses and target sizes. The model provided quantitative insight on dosimetric treatment variables over
Time-dependent finite-element models of phase-change problems with moving heat sources
Westerberg, K.W. ); Wiklof, C. ); Finlayson, B.A. . Dept. of Chemical Engineering)
1994-03-01
A mathematical model is developed for melting of a multilayered medium while a heat source traverses one boundary. The finite-element method uses moving meshes, front-tracking using spines, an automatic time-step algorithm, and an efficient solution of the linearized equations. A novel solution method allows the fixed-mesh code to work unchanged but allows a moving mesh in other problems. The finite-element method is applied when the heater mesh moves with respect to the multilayered medium mesh. The same technique allows parallel processing for finite-element codes. The model is applied to several test problems and then to the title problem.
Progress in Developing Finite Element Models Replicating Flexural Graphite Testing
Robert Bratton
2010-06-01
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.
Surface photovoltage measurements and finite element modeling of SAW devices.
Donnelly, Christine
2012-03-01
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.
CUERVO: A finite element computer program for nonlinear scalar transport problems
Sirman, M.B.; Gartling, D.K.
1995-11-01
CUERVO is a finite element code that is designed for the solution of multi-dimensional field problems described by a general nonlinear, advection-diffusion equation. The code is also applicable to field problems described by diffusion, Poisson or Laplace equations. The finite element formulation and the associated numerical methods used in CUERVO are outlined here; detailed instructions for use of the code are also presented. Example problems are provided to illustrate the use of the code.
Biffle, J.H.
1993-02-01
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.
Diffusive mesh relaxation in ALE finite element numerical simulations
Dube, E.I.
1996-06-01
The theory for a diffusive mesh relaxation algorithm is developed for use in three-dimensional Arbitary Lagrange/Eulerian (ALE) finite element simulation techniques. This mesh relaxer is derived by a variational principle for an unstructured 3D grid using finite elements, and incorporates hourglass controls in the numerical implementation. The diffusive coefficients are based on the geometric properties of the existing mesh, and are chosen so as to allow for a smooth grid that retains the general shape of the original mesh. The diffusive mesh relaxation algorithm is then applied to an ALE code system, and results from several test cases are discussed.
Experimentally validated finite element model of electrocaloric multilayer ceramic structures
Smith, N. A. S. E-mail: maciej.rokosz@npl.co.uk Correia, T. M. E-mail: maciej.rokosz@npl.co.uk; Rokosz, M. K. E-mail: maciej.rokosz@npl.co.uk
2014-07-28
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.
Assessment of a hybrid finite element and finite volume code for turbulent incompressible flows
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Xia, Yidong; Wang, Chuanjin; Luo, Hong; Christon, Mark; Bakosi, Jozsef
2015-12-15
Hydra-TH is a hybrid finite-element/finite-volume incompressible/low-Mach flow simulation code based on the Hydra multiphysics toolkit being developed and used for thermal-hydraulics applications. In the present work, a suite of verification and validation (V&V) test problems for Hydra-TH was defined to meet the design requirements of the Consortium for Advanced Simulation of Light Water Reactors (CASL). The intent for this test problem suite is to provide baseline comparison data that demonstrates the performance of the Hydra-TH solution methods. The simulation problems vary in complexity from laminar to turbulent flows. A set of RANS and LES turbulence models were used in themoreÂ Â» simulation of four classical test problems. Numerical results obtained by Hydra-TH agreed well with either the available analytical solution or experimental data, indicating the verified and validated implementation of these turbulence models in Hydra-TH. Where possible, we have attempted some form of solution verification to identify sensitivities in the solution methods, and to suggest best practices when using the Hydra-TH code.Â«Â less
Assessment of a hybrid finite element and finite volume code for turbulent incompressible flows
Xia, Yidong; Wang, Chuanjin; Luo, Hong; Christon, Mark; Bakosi, Jozsef
2015-12-15
Hydra-TH is a hybrid finite-element/finite-volume incompressible/low-Mach flow simulation code based on the Hydra multiphysics toolkit being developed and used for thermal-hydraulics applications. In the present work, a suite of verification and validation (V&V) test problems for Hydra-TH was defined to meet the design requirements of the Consortium for Advanced Simulation of Light Water Reactors (CASL). The intent for this test problem suite is to provide baseline comparison data that demonstrates the performance of the Hydra-TH solution methods. The simulation problems vary in complexity from laminar to turbulent flows. A set of RANS and LES turbulence models were used in the simulation of four classical test problems. Numerical results obtained by Hydra-TH agreed well with either the available analytical solution or experimental data, indicating the verified and validated implementation of these turbulence models in Hydra-TH. Where possible, we have attempted some form of solution verification to identify sensitivities in the solution methods, and to suggest best practices when using the Hydra-TH code.
Finite Element Solver for Fission Dynamics
Energy Science and Technology Software Center (OSTI)
2015-01-30
FELIX is a physics computer code used to model fission fragment mass distributions in a fully quantum-mechanical, misroscopic framework that only relies on our current knowledge of nuclear forces. It is an implementation of the time-dependent generator coordinate method (TDGCM), which simulates the dynamics of a collective quantum wave-packet assuming the motion is adiabatic. In typical applications of the TDGCM, the nuclear collective wavepacket is obtained as a superposition of wavefunctions obtained by solving themoreÂ Â» Hartree-Fock-Bogoliubov equations of nuclear density functional theory (DFT). The program calculates at each time step the coefficients of that superposition.Â«Â less
FUEL ELEMENT FABRICATION METHOD
Hix, J.N.; Cooley, G.E.; Cunningham, J.E.
1960-05-31
A method is given for assembling and fabricating a fuel element comprising a plurality of spaced parallel fuel plates of a bowed configuration supported by and between a pair of transperse aluminum side plates. In this method, a brasing alloy is preplated on one surface of the aluminum side plates in the form of a cladding or layer-of uniform thickness. Grooves are then cut into the side plates through the alloy layer and into the base aluminum which results in the utilization of thinner aluminum side plates since a portion of the necessary groove depth is supplied by the brazing alloy.
Finite-element analysis of an epoxy-curing process
Gartling, D K; Hickox, C E; Nunziato, J W
1983-01-01
A finite element numerical procedure is used to study the curing of an epoxy compound. The problem involves the gelation of an incompressible liquid due to an exothermic chemical reaction. Nonuniform temperature fields produce buoyancy-driven fluid motions that interact with the solidifying material. The numerical simulations provide temperature histories and the progression of the gel front that are compared with experimental data.
Nonlinear Legendre Spectral Finite Elements for Wind Turbine Blade Dynamics: Preprint
Wang, Q.; Sprague, M. A.; Jonkman, J.; Johnson, N.
2014-01-01
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.
An Annotated Reference Guide to the Finite-Element Interface Specification Version 1.0
Alan B. Williams; Ivan J. Otero; Kyran D. Mish; Lee M. Tayor; Robert L. Clay
1999-04-01
The Finite-Element Interface (FEI) specification provides a layered abstraction that permits finite-element analysis codes to utilize various linear-algebra solution packages with minimal concern for the internal details of the solver modules. Alternatively, this interface can be viewed as a way for solver developers to provide solution services to finite-element clients without having to embed finite-element abstractions within their solver libraries. The purpose of this document is to provide some level of documentation between the bare interface specification itself, which consists only of C/C++ header files, and the full documentation suite that supports the interface definition by providing considerable detail as to its design and implementation. This document primarily provides the ''how'' of calling the interface member functions, so that programmers can readily learn how to utilize the interface implementation without having to consider all the details contained in the interface's definition, design, and motivation. The interface specification is presented three times in this document, each time with an increasing level of detail. The first presentation provides a general overview of the calling sequence, in order to acquaint the programmer with a basic introduction to how the interface is used to ''train'' the underlying solver software on the particular finite-element problem that is to be solved. The second pass through the interface definition provides considerable detail on each method, including specific considerations as to the structure of the underlying data, and an exposition of potential pitfalls that may occur as a byproduct of either the finite-element modeling process, or of the use of the associated interface implementation. Finally, a third description of the interface is given implicitly via the discussion of sample problems that provide concrete examples of the use of the finite-element interface.
A Moving Window Technique in Parallel Finite Element Time Domain Electromagnetic Simulation
Lee, Lie-Quan; Candel, Arno; Ng, Cho; Ko, Kwok; ,
2010-06-07
A moving window technique for the finite element time domain (FETD) method is developed to simulate the propagation of electromagnetic waves induced by the transit of a charged particle beam inside large and long structures. The window moving along with the beam in the computational domain adopts high-order finite-element basis functions through p refinement and/or a high-resolution mesh through h refinement so that a sufficient accuracy is attained with substantially reduced computational costs. Algorithms to transfer discretized fields from one mesh to another, which are the key to implementing a moving window in a finite-element unstructured mesh, are presented. Numerical experiments are carried out using the moving window technique to compute short-range wakefields in long accelerator structures. The results are compared with those obtained from the normal FETD method and the advantages of using the moving window technique are discussed.
Residual stresses in weld overlay tubes: A finite element study
Taljat, B.; Zacharia, T.; Wang, X.L.; Keiser, J.R.; Feng, Z.; Jirinec, M.J.
1997-01-03
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.
Beam and Truss Finite Element Verification for DYNA3D
Rathbun, H J
2007-07-16
The explicit finite element (FE) software program DYNA3D has been developed at Lawrence Livermore National Laboratory (LLNL) to simulate the dynamic behavior of structures, systems, and components. This report focuses on verification of beam and truss element formulations in DYNA3D. An efficient protocol has been developed to verify the accuracy of these structural elements by generating a set of representative problems for which closed-form quasi-static steady-state analytical reference solutions exist. To provide as complete coverage as practically achievable, problem sets are developed for each beam and truss element formulation (and their variants) in all modes of loading and physical orientation. Analyses with loading in the elastic and elastic-plastic regimes are performed. For elastic loading, the FE results are within 1% of the reference solutions for all cases. For beam element bending and torsion loading in the plastic regime, the response is heavily dependent on the numerical integration rule chosen, with higher refinement yielding greater accuracy (agreement to within 1%). Axial loading in the plastic regime produces accurate results (agreement to within 0.01%) for all integration rules and element formulations. Truss elements are also verified to provide accurate results (within 0.01%) for elastic and elastic-plastic loading. A sample problem to verify beam element response in ParaDyn, the parallel version DYNA3D, is also presented.
Goudreau, G.L.
1983-05-01
The confrontation of the finite element technology with the awesome number crunching required for the nonlinear problem has forced a new assessment. The finite element community discovered that for impact, wave propagation and even slower impulsively driven dynamic problems, the finite difference hydrocodes were much more efficient, albeit requiring contorted mesh topologies to model practical geometries and fifty thousand time steps to follow a quasi-static process or reach static equilibrium. The experience of engineering analysis at the Lawrence Livermore National Laboratory has straddled both worlds, and the intent of this lecture is to review choices in the context of the two and three dimensional implicit and explicit Lagrangian codes developed in our Methods Development Group. The explicit DYNA2D and DYNA3D and implicit NIKE2D and NIKE3D of Hallquist form the focus of this discussion.
Finite element analyses for seismic shear wall international standard problem
Park, Y.J.; Hofmayer, C.H.
1998-04-01
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.
Finite-element analysis of the deformation of thin Mylar films...
Office of Scientific and Technical Information (OSTI)
Finite-element analysis of the deformation of thin Mylar films due to measurement forces. Citation Details In-Document Search Title: Finite-element analysis of the deformation of ...
An enriched finite element model with q-refinement for radiative boundary layers in glass cooling
Mohamed, M. Shadi; Seaid, Mohammed; Trevelyan, Jon; Laghrouche, Omar
2014-02-01
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.
On the validation of seismic imaging methods: Finite frequency...
Office of Scientific and Technical Information (OSTI)
On the validation of seismic imaging methods: Finite ... Our study suggests that finite-frequency approaches to ... Country of Publication: United States Language: English ...
A Finite Element Model for Simulation of Carbon Dioxide Sequestration
Bao, Jie; Xu, Zhijie; Fang, Yilin
2013-11-02
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.
Barham, M; White, D; Steigmann, D; Rudd, R
2009-04-08
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.
Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report.
Siefert, Christopher; Bochev, Pavel Blagoveston; Kramer, Richard Michael Jack; Voth, Thomas Eugene; Cox, James
2014-09-01
Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include explod- ing bridge wires for detonator design, electromagnetic launch of flyer plates for material testing and gun design, lightning blast-through for weapon safety, electromagnetic armor, and magnetic flux compression generators, all require accurate resolution of surface effects. These applications operate in a large deformation regime, where body-fitted meshes are im- practical and multimaterial elements are the only feasible option. State-of-the-art methods use various mixture models to approximate the multi-physics of these elements. The em- pirical nature of these models can significantly compromise the accuracy of the simulation in this very important surface region. We propose to substantially improve the predictive capability of electromagnetic simulations by removing the need for empirical mixture models at material surfaces. We do this by developing an eXtended Finite Element Method (XFEM) and an associated Conformal Decomposition Finite Element Method (CDFEM) which sat- isfy the physically required compatibility conditions at material interfaces. We demonstrate the effectiveness of these methods for diffusion and diffusion-like problems on node, edge and face elements in 2D and 3D. We also present preliminary work on h -hierarchical elements and remap algorithms.
Smith, Jovanca J.; Bishop, Joseph E.
2013-11-01
This report summarizes the work performed by the graduate student Jovanca Smith during a summer internship in the summer of 2012 with the aid of mentor Joe Bishop. The projects were a two-part endeavor that focused on the use of the numerical model called the Lattice Discrete Particle Model (LDPM). The LDPM is a discrete meso-scale model currently used at Northwestern University and the ERDC to model the heterogeneous quasi-brittle material, concrete. In the first part of the project, LDPM was compared to the Karagozian and Case Concrete Model (K&C) used in Presto, an explicit dynamics finite-element code, developed at Sandia National Laboratories. In order to make this comparison, a series of quasi-static numerical experiments were performed, namely unconfined uniaxial compression tests on four varied cube specimen sizes, three-point bending notched experiments on three proportional specimen sizes, and six triaxial compression tests on a cylindrical specimen. The second part of this project focused on the application of LDPM to simulate projectile perforation on an ultra high performance concrete called CORTUF. This application illustrates the strengths of LDPM over traditional continuum models.
Shadid, J.N.; Moffat, H.K.; Hutchinson, S.A.; Hennigan, G.L.; Devine, K.D.; Salinger, A.G.
1996-05-01
The theoretical background for the finite element computer program, MPSalsa, is presented in detail. MPSalsa is designed to solve laminar, low Mach number, two- or three-dimensional incompressible and variable density reacting fluid flows on massively parallel computers, using a Petrov-Galerkin finite element formulation. The code has the capability to solve coupled fluid flow, heat transport, multicomponent species transport, and finite-rate chemical reactions, and to solver coupled multiple Poisson or advection-diffusion- reaction equations. The program employs the CHEMKIN library to provide a rigorous treatment of multicomponent ideal gas kinetics and transport. Chemical reactions occurring in the gas phase and on surfaces are treated by calls to CHEMKIN and SURFACE CHEMKIN, respectively. The code employs unstructured meshes, using the EXODUS II finite element data base suite of programs for its input and output files. MPSalsa solves both transient and steady flows by using fully implicit time integration, an inexact Newton method and iterative solvers based on preconditioned Krylov methods as implemented in the Aztec solver library.
Wakefield Simulation of CLIC PETS Structure Using Parallel 3D Finite Element Time-Domain Solver T3P
Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; Ko, K.; Syratchev, I.; /CERN
2009-06-19
In recent years, SLAC's Advanced Computations Department (ACD) has developed the parallel 3D Finite Element electromagnetic time-domain code T3P. Higher-order Finite Element methods on conformal unstructured meshes and massively parallel processing allow unprecedented simulation accuracy for wakefield computations and simulations of transient effects in realistic accelerator structures. Applications include simulation of wakefield damping in the Compact Linear Collider (CLIC) power extraction and transfer structure (PETS).
FINITE ELEMENT ANALYSIS OF BULK TRITIUM SHIPPING PACKAGE
Jordan, J.
2010-06-02
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.
Finite element analyses of a linear-accelerator electron gun
Iqbal, M. E-mail: muniqbal@ihep.ac.cn; Wasy, A.; Islam, G. U.; Zhou, Z.
2014-02-15
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.
Three-dimensional finite element analysis of a bolted joint
Lin, H.; Reddy, H.; McKee, R.B.
1995-11-01
A flat steel flange clamped by 14 bolts was modeled by three dimensional finite element analysis to evaluate the apparent stiffness in compression of the flange. The bolt load was applied via an extremely stiff washer to eliminate the effect of bolt head deflection on the results. The dimensionless stiffness was found to vary with bolt aspect ratio d/L in a nearly linear fashion, with slope equal to that predicted by a simple cylindrical model. The conical frustum model, with a half angle of 32.8 degrees, predicted stiffnesses within a few percent of those measured. The interflange pressure along the bolt circle was calculated for a common bolt spacing and a range of aspect ratios.
Multi-element microelectropolishing method
Lee, Peter J.
1994-01-01
A method is provided for microelectropolishing a transmission electron microscopy nonhomogeneous multi-element compound foil. The foil is electrolyzed at different polishing rates for different elements by rapidly cycling between different current densities. During a first portion of each cycle at a first voltage a first element electrolyzes at a higher current density than a second element such that the material of the first element leaves the anode foil at a faster rate than the second element and creates a solid surface film, and such that the solid surface film is removed at a faster rate than the first element leaves the anode foil. During a second portion of each cycle at a second voltage the second element electrolyzes at a higher current density than the first element, and the material of the second element leaves the anode foil at a faster rate than the first element and creates a solid surface film, and the solid surface film is removed at a slower rate than the second element leaves the foil. The solid surface film is built up during the second portion of the cycle, and removed during the first portion of the cycle.
Multi-element microelectropolishing method
Lee, P.J.
1994-10-11
A method is provided for microelectropolishing a transmission electron microscopy nonhomogeneous multi-element compound foil. The foil is electrolyzed at different polishing rates for different elements by rapidly cycling between different current densities. During a first portion of each cycle at a first voltage a first element electrolyzes at a higher current density than a second element such that the material of the first element leaves the anode foil at a faster rate than the second element and creates a solid surface film, and such that the solid surface film is removed at a faster rate than the first element leaves the anode foil. During a second portion of each cycle at a second voltage the second element electrolyzes at a higher current density than the first element, and the material of the second element leaves the anode foil at a faster rate than the first element and creates a solid surface film, and the solid surface film is removed at a slower rate than the second element leaves the foil. The solid surface film is built up during the second portion of the cycle, and removed during the first portion of the cycle. 10 figs.
A phenomenological finite element model of part building in the stereolithography process
Chambers, R.S.; Guess, T.R.; Hinnerichs, T.D.
1995-03-01
The finite element method has been used to develop the framework for a tool that can be used to model the structural deformation arising from the stereolithography build process. Such a tool when fully developed can facilitate numerical studies aimed at evaluating build parameters and build styles. Although the current software makes no attempt to capture all the physics of the process, provisions for three important build features have been made: (1) laser path history including scanning rate and depth of cure, (2) structural linkage, and (3) time varying material behavior. For demonstration purposes, a three dimensional finite element code was modified to include a phenomenological material model of solidification. The model was based on cure shrinkage and stress relaxation data collected from in-situ tests on individual strands drawn using 3D Systems` stereolithography apparatus (SLA-250). To depict the directed path of solidification within layers, a finite element birthing scheme was conceived to activate elements along the predetermined coordinate path of the laser. Structural linkage was enforced by joining element strands of layers when laser paths connect or overlap, respectively. A limited number of analyses have been performed to contrast simple build styles.
Bailey, T S; Chang, J H; Warsa, J S; Adams, M L
2010-12-22
We present a new spatial discretization of the discrete-ordinates transport equation in two-dimensional Cartesian (X-Y) geometry for arbitrary polygonal meshes. The discretization is a discontinuous finite element method (DFEM) that utilizes piecewise bi-linear (PWBL) basis functions, which are formally introduced in this paper. We also present a series of numerical results on quadrilateral and polygonal grids and compare these results to a variety of other spatial discretizations that have been shown to be successful on these grid types. Finally, we note that the properties of the PWBL basis functions are such that the leading-order piecewise bi-linear discontinuous finite element (PWBLD) solution will satisfy a reasonably accurate diffusion discretization in the thick diffusion limit, making the PWBLD method a viable candidate for many different classes of transport problems.
Electrogalvanic finite element analysis of partially protected marine structures
Kasper, R.G.; April, M.G.
1983-05-01
A general finite element modeling (FEM) procedure for calculating electrogalvanic field responses due to multiple anodic/cathodic interactions has been developed. Within any defined electrolytic continuum the exact geometry and location of anodes, cathodes, and paint surfaces can now be realistically incorporated in the mathematical model formulation. The anodic/cathodic interactions in the conductive electrolyte are predicted by the application of classical d.c. electric field theory for conductive continuums in conjunction with widely accepted laboratory oxidation/reduction responses for the electrodes. The electrogalvanic fields in the electrolyte are calculated using the scalar Poisson equation whereby traditional boundary conditions are prescribed in the far field of the electrolyte. In the near field of the anodes, cathodes, and the painted metallic substrate, complex boundary conditions are enforced based on empirical polarization curves and paint impedance values. The ionic current in the electrolyte leaving the anode and arriving at the cathodes are mathematically constrained to sum to zero over the metallic surface (spatial Kirchoff's law). Three example problems (three dimensional) of an iron bar in salt water will demonstrate the electrochemical field effects due to relative anodic/cathodic spacing, anodic/cathodic geometries, and infinite paint impedance. A fourth boundary value problem (two dimensional) will illustrate the effect of a paint discontinuity (unprotected area) in the center of the otherwise painted (finite impedance) metallic substrate. The proposed electrochemical modeling procedure has the advantage of concurrently treating numerous essential electrochemical parameters for any geometric anodic/cathodic configurations consistent with the physical laws of electrochemistry.
Parallel 3D Finite Element Particle-in-Cell Simulations with Pic3P
Candel, A.; Kabel, A.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; Ko, K.; Ben-Zvi, I.; Kewisch, J.; /Brookhaven
2009-06-19
SLAC's Advanced Computations Department (ACD) has developed the parallel 3D Finite Element electromagnetic Particle-In-Cell code Pic3P. Designed for simulations of beam-cavity interactions dominated by space charge effects, Pic3P solves the complete set of Maxwell-Lorentz equations self-consistently and includes space-charge, retardation and boundary effects from first principles. Higher-order Finite Element methods with adaptive refinement on conformal unstructured meshes lead to highly efficient use of computational resources. Massively parallel processing with dynamic load balancing enables large-scale modeling of photoinjectors with unprecedented accuracy, aiding the design and operation of next-generation accelerator facilities. Applications include the LCLS RF gun and the BNL polarized SRF gun.
Finite Element Simulations to Explore Assumptions in Kolsky Bar Experiments.
Crum, Justin
2015-08-05
The chief purpose of this project has been to develop a set of finite element models that attempt to explore some of the assumptions in the experimental set-up and data reduction of the Kolsky bar experiment. In brief, the Kolsky bar, sometimes referred to as the split Hopkinson pressure bar, is an experimental apparatus used to study the mechanical properties of materials at high strain rates. Kolsky bars can be constructed to conduct experiments in tension or compression, both of which are studied in this paper. The basic operation of the tension Kolsky bar is as follows: compressed air is inserted into the barrel that contains the striker; the striker accelerates towards the left and strikes the left end of the barrel producing a tensile stress wave that propogates first through the barrel and then down the incident bar, into the specimen, and finally the transmission bar. In the compression case, the striker instead travels to the right and impacts the incident bar directly. As the stress wave travels through an interface (e.g., the incident bar to specimen connection), a portion of the pulse is transmitted and the rest reflected. The incident pulse, as well as the transmitted and reflected pulses are picked up by two strain gauges installed on the incident and transmitted bars as shown. By interpreting the data acquired by these strain gauges, the stress/strain behavior of the specimen can be determined.
ELLIPT2D: A Flexible Finite Element Code Written Python
Pletzer, A.; Mollis, J.C.
2001-03-22
The use of the Python scripting language for scientific applications and in particular to solve partial differential equations is explored. It is shown that Python's rich data structure and object-oriented features can be exploited to write programs that are not only significantly more concise than their counter parts written in Fortran, C or C++, but are also numerically efficient. To illustrate this, a two-dimensional finite element code (ELLIPT2D) has been written. ELLIPT2D provides a flexible and easy-to-use framework for solving a large class of second-order elliptic problems. The program allows for structured or unstructured meshes. All functions defining the elliptic operator are user supplied and so are the boundary conditions, which can be of Dirichlet, Neumann or Robbins type. ELLIPT2D makes extensive use of dictionaries (hash tables) as a way to represent sparse matrices.Other key features of the Python language that have been widely used include: operator over loading, error handling, array slicing, and the Tkinter module for building graphical use interfaces. As an example of the utility of ELLIPT2D, a nonlinear solution of the Grad-Shafranov equation is computed using a Newton iterative scheme. A second application focuses on a solution of the toroidal Laplace equation coupled to a magnetohydrodynamic stability code, a problem arising in the context of magnetic fusion research.
Method for detecting an element
Blackwood, Larry G.; Reber, Edward L.; Rohde, Kenneth W.
2007-02-06
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.
SEACAS Theory Manuals: Part III. Finite Element Analysis in Nonlinear Solid Mechanics
Laursen, T.A.; Attaway, S.W.; Zadoks, R.I.
1999-03-01
This report outlines the application of finite element methodology to large deformation solid mechanics problems, detailing also some of the key technological issues that effective finite element formulations must address. The presentation is organized into three major portions: first, a discussion of finite element discretization from the global point of view, emphasizing the relationship between a virtual work principle and the associated fully discrete system, second, a discussion of finite element technology, emphasizing the important theoretical and practical features associated with an individual finite element; and third, detailed description of specific elements that enjoy widespread use, providing some examples of the theoretical ideas already described. Descriptions of problem formulation in nonlinear solid mechanics, nonlinear continuum mechanics, and constitutive modeling are given in three companion reports.
Finite element meshing approached as a global minimization process
WITKOWSKI,WALTER R.; JUNG,JOSEPH; DOHRMANN,CLARK R.; LEUNG,VITUS J.
2000-03-01
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
FEHMN 1.0: Finite element heat and mass transfer code; Revision 1
Zyvoloski, G.; Dash, Z.; Kelkar, S.
1992-05-01
A computer code is described which can simulate non-isothermal multi-phase multicomponent flow in porous media. It is applicable to natural-state studies of geothermal systems and groundwater flow. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved sing the finite element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat and mass transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. A summary of the equations in the model and the numerical solution procedure are provided in this report. A user`s guide and sample problems are also included. The FEHMN (Finite Element Heat and Mass Nuclear) code, described in this report, is a version of FEHM (Finite Element Heat and Mass, Zyvoloski et al., 1988) developed for the Yucca Mountain Site Characterization Project (YMP). The main use of FEHMN will be to assist in the understanding of flow fields in the saturated zone below the potential Yucca Mountain repository.
Grid Generator for Two, Three-dimensional Finite Element Subsurface Flow Models
Energy Science and Technology Software Center (OSTI)
1993-04-28
GRIDMAKER serves as a preprocessor for finite element models in solving two- and three-dimensional subsurface flow and pollutant transport problems. It is designed to generate three-point triangular or four-point quadrilateral elements for two-dimensional domains and eight-point hexahedron elements for three-dimensional domains. A two-dimensional domain of an aquifer with a variable depth layer is treated as a special case for depth-integrated two-dimensional, finite element subsurface flow models. The program accommodates the need for aquifers with heterogeneousmoreÂ Â» systems by identifying the type of material in each element.Â«Â less
Hexahedral finite element mesh coarsening using pillowing technique
Staten, Matthew L.; Woodbury, Adam C.; Benzley, Steven E.; Shepherd, Jason F.
2012-06-05
A techniques for coarsening a hexahedral mesh is described. The technique includes identifying a coarsening region within a hexahedral mesh to be coarsened. A boundary sheet of hexahedral elements is inserted into the hexahedral mesh around the coarsening region. A column of hexahedral elements is identified within the boundary sheet. The column of hexahedral elements is collapsed to create an extraction sheet of hexahedral elements contained within the coarsening region. Then, the extraction sheet of hexahedral elements is extracted to coarsen the hexahedral mesh.
Coupled discrete element and finite volume solution of two classical soil mechanics problems
Chen, Feng; Drumm, Eric; Guiochon, Georges A
2011-01-01
One dimensional solutions for the classic critical upward seepage gradient/quick condition and the time rate of consolidation problems are obtained using coupled routines for the finite volume method (FVM) and discrete element method (DEM), and the results compared with the analytical solutions. The two phase flow in a system composed of fluid and solid is simulated with the fluid phase modeled by solving the averaged Navier-Stokes equation using the FVM and the solid phase is modeled using the DEM. A framework is described for the coupling of two open source computer codes: YADE-OpenDEM for the discrete element method and OpenFOAM for the computational fluid dynamics. The particle-fluid interaction is quantified using a semi-empirical relationship proposed by Ergun [12]. The two classical verification problems are used to explore issues encountered when using coupled flow DEM codes, namely, the appropriate time step size for both the fluid and mechanical solution processes, the choice of the viscous damping coefficient, and the number of solid particles per finite fluid volume.
A high-order staggered finite-element vertical discretization for non-hydrostatic atmospheric models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Guerra, Jorge E.; Ullrich, Paul A.
2016-06-01
Atmospheric modeling systems require economical methods to solve the non-hydrostatic Euler equations. Two major differences between hydrostatic models and a full non-hydrostatic description lies in the vertical velocity tendency and numerical stiffness associated with sound waves. In this work we introduce a new arbitrary-order vertical discretization entitled the staggered nodal finite-element method (SNFEM). Our method uses a generalized discrete derivative that consistently combines the discontinuous Galerkin and spectral element methods on a staggered grid. Our combined method leverages the accurate wave propagation and conservation properties of spectral elements with staggered methods that eliminate stationary (2Î”x) modes. Furthermore, high-order accuracy alsomoreÂ Â» eliminates the need for a reference state to maintain hydrostatic balance. In this work we demonstrate the use of high vertical order as a means of improving simulation quality at relatively coarse resolution. We choose a test case suite that spans the range of atmospheric flows from predominantly hydrostatic to nonlinear in the large-eddy regime. Lastly, our results show that there is a distinct benefit in using the high-order vertical coordinate at low resolutions with the same robust properties as the low-order alternative.Â«Â less
A high-order staggered finite-element vertical discretization for non-hydrostatic atmospheric models
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Guerra, Jorge E.; Ullrich, Paul A.
2016-06-01
Atmospheric modeling systems require economical methods to solve the non-hydrostatic Euler equations. Two major differences between hydrostatic models and a full non-hydrostatic description lies in the vertical velocity tendency and numerical stiffness associated with sound waves. In this work we introduce a new arbitrary-order vertical discretization entitled the staggered nodal finite-element method (SNFEM). Our method uses a generalized discrete derivative that consistently combines the discontinuous Galerkin and spectral element methods on a staggered grid. Our combined method leverages the accurate wave propagation and conservation properties of spectral elements with staggered methods that eliminate stationary (2Î”x) modes. Furthermore, high-order accuracy alsomoreÂ Â» eliminates the need for a reference state to maintain hydrostatic balance. In this work we demonstrate the use of high vertical order as a means of improving simulation quality at relatively coarse resolution. We choose a test case suite that spans the range of atmospheric flows from predominantly hydrostatic to nonlinear in the large-eddy regime. Our results show that there is a distinct benefit in using the high-order vertical coordinate at low resolutions with the same robust properties as the low-order alternative.Â«Â less
The Uranium Processing Facility Finite Element Meshing Discussion
Office of Environmental Management (EM)
2 - 15'-20' element size deemed adequate as the building geometry was rather straightforward. - Linear dynamic computer model October 25, 2011 8 Department of Energy -...
XU, J.; COSTANTINO, C.; HOFMAYER, C.
2006-06-26
PAPER DISCUSSES COMPUTATIONS OF SEISMIC INDUCED SOIL PRESSURES USING FINITE ELEMENT MODELS FOR DEEPLY EMBEDDED AND OR BURIED STIFF STRUCTURES SUCH AS THOSE APPEARING IN THE CONCEPTUAL DESIGNS OF STRUCTURES FOR ADVANCED REACTORS.
METHOD OF MAKING FUEL ELEMENTS
Bean, C.H.; Macherey, R.E.
1959-12-01
A method is described for fabricating fuel elements, particularly for enclosing a plate of metal with a second metal by inserting the plate into an aperture of a frame of a second plate, placing a sheet of the second metal on each of opposite faces of the assembled plate and frame, purging with an inert gas the air from the space within the frame and the sheets while sealing the seams between the frame and the sheets, exhausting the space, purging the space with air, re-exhausting the spaces, sealing the second aperture, and applying heat and pressure to bond the sheets, the plate, and the frame to one another.
A phase field - Finite element approach to model the interaction between
Office of Scientific and Technical Information (OSTI)
phase transformations and plasticity in shape memory alloys (Journal Article) | SciTech Connect phase field - Finite element approach to model the interaction between phase transformations and plasticity in shape memory alloys Citation Details In-Document Search This content will become publicly available on March 22, 2018 Title: A phase field - Finite element approach to model the interaction between phase transformations and plasticity in shape memory alloys Authors: Paranjape, Harshad M.
Mimetic finite difference method for the stokes problem on polygonal meshes
Lipnikov, K; Beirao Da Veiga, L; Gyrya, V; Manzini, G
2009-01-01
Various approaches to extend the finite element methods to non-traditional elements (pyramids, polyhedra, etc.) have been developed over the last decade. Building of basis functions for such elements is a challenging task and may require extensive geometry analysis. The mimetic finite difference (MFD) method has many similarities with low-order finite element methods. Both methods try to preserve fundamental properties of physical and mathematical models. The essential difference is that the MFD method uses only the surface representation of discrete unknowns to build stiffness and mass matrices. Since no extension inside the mesh element is required, practical implementation of the MFD method is simple for polygonal meshes that may include degenerate and non-convex elements. In this article, we develop a MFD method for the Stokes problem on arbitrary polygonal meshes. The method is constructed for tensor coefficients, which will allow to apply it to the linear elasticity problem. The numerical experiments show the second-order convergence for the velocity variable and the first-order for the pressure.
Finite-element analysis of earing using non-quadratic yield surfaces
Logan, R.W.
1995-06-18
During deep draw cupping, the phenomenon known as earing may occur as the cup wall is formed, resulting in a periodic variation of cup wall height around the perimeter of the finished cup. This is generally due to planar anisotropy of flow in rolled sheet product. It is generally observed that the anisotropy parameter R will vary in the plane of the sheet when ears are observed in cupping, with a parameter {Delta}R describing the variation of R in the plane of the sheet. For many common textures in face-centered and body-centered materials, the ears form relative to the sheet rolling direction at 0{degrees} and 90{degrees} around the perimeter if {Delta}R>0, and at -45{degrees} and +45{degrees} if {Delta}R<0. There is extensive experimental evidence that ear height shows a linear correlation with {Delta}R/R, but attempts to duplicate this using the finite-element method are highly dependent on both the methodology and yield surface used. It was shown previously that using a coarse mesh and the quadratic Hill yield surface tends to greatly under predict earing. In this study, we have used two different finite-element codes developed at LLNL to examine the predicted earing using both quadratic Hill and alternative non-quadratic yield surfaces. These results are compared to experimental data and conclusions drawn about the most desirable closed-form yield surfaces to duplicate the observed earing phenomena.
Kolotilina, L.; Nikishin, A.; Yeremin, A.
1994-12-31
The solution of large systems of linear equations is a crucial bottleneck when performing 3D finite element analysis of structures. Also, in many cases the reliability and robustness of iterative solution strategies, and their efficiency when exploiting hardware resources, fully determine the scope of industrial applications which can be solved on a particular computer platform. This is especially true for modern vector/parallel supercomputers with large vector length and for modern massively parallel supercomputers. Preconditioned iterative methods have been successfully applied to industrial class finite element analysis of structures. The construction and application of high quality preconditioners constitutes a high percentage of the total solution time. Parallel implementation of high quality preconditioners on such architectures is a formidable challenge. Two common types of existing preconditioners are the implicit preconditioners and the explicit preconditioners. The implicit preconditioners (e.g. incomplete factorizations of several types) are generally high quality but require solution of lower and upper triangular systems of equations per iteration which are difficult to parallelize without deteriorating the convergence rate. The explicit type of preconditionings (e.g. polynomial preconditioners or Jacobi-like preconditioners) require sparse matrix-vector multiplications and can be parallelized but their preconditioning qualities are less than desirable. The authors present results of numerical experiments with Factorized Sparse Approximate Inverses (FSAI) for symmetric positive definite linear systems. These are high quality preconditioners that possess a large resource of parallelism by construction without increasing the serial complexity.
ROMERO,VICENTE J.; SWILER,LAURA PAINTON; GIUNTA,ANTHONY A.
2000-04-25
This paper examines the modeling accuracy of finite element interpolation, kriging, and polynomial regression used in conjunction with the Progressive Lattice Sampling (PLS) incremental design-of-experiments approach. PLS is a paradigm for sampling a deterministic hypercubic parameter space by placing and incrementally adding samples in a manner intended to maximally reduce lack of knowledge in the parameter space. When combined with suitable interpolation methods, PLS is a formulation for progressive construction of response surface approximations (RSA) in which the RSA are efficiently upgradable, and upon upgrading, offer convergence information essential in estimating error introduced by the use of RSA in the problem. The three interpolation methods tried here are examined for performance in replicating an analytic test function as measured by several different indicators. The process described here provides a framework for future studies using other interpolation schemes, test functions, and measures of approximation quality.
Cofer, W.F.
1992-03-01
The microplane concrete material model is based upon assumptions regarding the behavior of the material components. At any point, the response to the strain tensor on arbitrarily oriented surfaces is considered. Simple, softening stress-strain relationships are assumed in directions perpendicular and parallel to the surfaces. The macroscopic material behavior is then composed of the sum of the effects. The implementation of this model into the explicit, nonlinear, dynamic finite element program, DYNA3D, is described. To avoid the spurious mesh sensitivity that accompanies material failure, a weighted integral strain averaging approach is used to ensure that softening is nonlocal. This method is shown to be effective for limiting the failure zone in a concrete rod subjected to an impulse loading. 36 refs., 7 figs.
Krabicka, J.; Yan, Y.
2009-08-15
Electrostatic sensors are used in certain industries for the flow measurement of pneumatically conveyed solids. However, despite various advances that have been made in recent years, relatively little information is known about the exact nature of the electrostatic charge induced onto the sensor electrode due to moving particles, which is dependent on electrode geometry, particle distribution, and particle velocity. This paper presents a novel approach to the study of the charge induced onto electrostatic sensors based on fitting a Lorentzian curve to the results of a finite-element model of the electrostatic sensor and pipeline. The modeling method is validated by comparing the modeling results of a nonintrusive circular electrode with an established analytical solution. The modeling results are used for in-depth analysis and informed design of a particular sensor configuration.
Parallel performance of a preconditioned CG solver for unstructured finite element applications
Shadid, J.N.; Hutchinson, S.A.; Moffat, H.K.
1994-06-01
A parallel unstructured finite element (FE) implementation designed for message passing machines is described. This implementation employs automated problem partitioning algorithms for load balancing unstructured grids, a distributed sparse matrix representation of the global finite element equations and a parallel conjugate gradient (CG) solver. In this paper a number of issues related to the efficient implementation of parallel unstructured mesh applications are presented. These include the differences between structured and unstructured mesh parallel applications, major communication kernels for unstructured CG solvers, automatic mesh partitioning algorithms, and the influence of mesh. partitioning metrics on parallel performance. Initial results are presented for example finite element (FE) heat transfer analysis applications on a 1024 processor nCUBE 2 hypercube. Results indicate over 95% scaled efficiencies are obtained for some large problems despite the required unstructured data communication.
Parallel performance of a preconditioned CG solver for unstructured finite element applications
Shadid, J.N.; Hutchinson, S.A.; Moffat, H.K.
1994-12-31
A parallel unstructured finite element (FE) implementation designed for message passing MIMD machines is described. This implementation employs automated problem partitioning algorithms for load balancing unstructured grids, a distributed sparse matrix representation of the global finite element equations and a parallel conjugate gradient (CG) solver. In this paper a number of issues related to the efficient implementation of parallel unstructured mesh applications are presented. These include the differences between structured and unstructured mesh parallel applications, major communication kernels for unstructured CG solvers, automatic mesh partitioning algorithms, and the influence of mesh partitioning metrics on parallel performance. Initial results are presented for example finite element (FE) heat transfer analysis applications on a 1024 processor nCUBE 2 hypercube. Results indicate over 95% scaled efficiencies are obtained for some large problems despite the required unstructured data communication.
DYNA3D Finite Element Analysis of Steam Explosion Loads on a Pedestal Wall Design
Noble, C R
2007-01-18
The objective of this brief report is to document the ESBWR pedestal wall finite element analyses that were performed as a quick turnaround effort in July 2005 at Lawrence Livermore National Laboratory and describe the assumptions and failure criteria used for these analyses [Ref 4]. The analyses described within are for the pedestal wall design that included an internal steel liner. The goal of the finite element analyses was to assist in determining the load carrying capacity of the ESBWR pedestal wall subjected to an impulsive pressure generated by a steam explosion.
On the validation of seismic imaging methods: Finite frequency or ray theory?
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Maceira, Monica; Larmat, Carene; Porritt, Robert W.; Higdon, David M.; Rowe, Charlotte A.; Allen, Richard M.
2015-01-23
We investigate the merits of the more recently developed finite-frequency approach to tomography against the more traditional and approximate ray theoretical approach for state of the art seismic models developed for western North America. To this end, we employ the spectral element method to assess the agreement between observations on real data and measurements made on synthetic seismograms predicted by the models under consideration. We check for phase delay agreement as well as waveform cross-correlation values. Based on statistical analyses on S wave phase delay measurements, finite frequency shows an improvement over ray theory. Random sampling using cross-correlation values identifiesmoreÂ Â» regions where synthetic seismograms computed with ray theory and finite-frequency models differ the most. Our study suggests that finite-frequency approaches to seismic imaging exhibit measurable improvement for pronounced low-velocity anomalies such as mantle plumes.Â«Â less
Conventional {delta}f-particle simulations of electromagnetic perturbations with finite elements
Mishchenko, Alexey; Hatzky, Roman; Koenies, Axel
2004-12-01
The possibility of electromagnetic particle-in-cell simulations with a conventional {delta}f approach is shown in slab geometry using finite elements. Both the ion-temperature-gradient driven mode and the shear Alfven wave are reproduced and benchmarked with the analytical linear dispersion relation. Particularly, the Alfven wave is simulated successfully at the limit k{sub perpendicular}{yields}0.
Finite Element Code For 3D-Hydraulic Fracture Propagation Equations (3-layer).
Energy Science and Technology Software Center (OSTI)
1992-03-24
HYFRACP3D is a finite element program for simulation of a pseudo three-dimensional fracture geometries with a two-dimensional planar solution. The model predicts the height, width and winglength over time for a hydraulic fracture propagating in a three-layered system of rocks with variable rock mechanics properties.
A Finite-Element Model for Simulation of Carbon Dioxide Sequestration
Bao, Jie; Xu, Zhijie; Fang, Yilin
2014-09-01
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.
Software Library for Storing and Retrieving Mesh and Results of Finite Element
Energy Science and Technology Software Center (OSTI)
1997-07-07
EXOII is a model developed to store and retrieve data for finite element analyses. It is used for preprocessing (problem definition), postprocessing (results visualization), as well as code to code transfer. An EXOII data file is a random access, machine independent, binary file that is written and read via C, C++, or Fortran library routines which comprise the Application Programming Interface (API).
Maginot, P. G.; Ragusa, J. C.; Morel, J. E.
2013-07-01
We examine several possible methods of mass matrix lumping for discontinuous finite element discrete ordinates transport using a Lagrange interpolatory polynomial trial space. Though positive outflow angular flux is guaranteed with traditional mass matrix lumping in a purely absorbing 1-D slab cell for the linear discontinuous approximation, we show that when used with higher degree interpolatory polynomial trial spaces, traditional lumping does yield strictly positive outflows and does not increase in accuracy with an increase in trial space polynomial degree. As an alternative, we examine methods which are 'self-lumping'. Self-lumping methods yield diagonal mass matrices by using numerical quadrature restricted to the Lagrange interpolatory points. Using equally-spaced interpolatory points, self-lumping is achieved through the use of closed Newton-Cotes formulas, resulting in strictly positive outflows in pure absorbers for odd power polynomials in 1-D slab geometry. By changing interpolatory points from the traditional equally-spaced points to the quadrature points of the Gauss-Legendre or Lobatto-Gauss-Legendre quadratures, it is possible to generate solution representations with a diagonal mass matrix and a strictly positive outflow for any degree polynomial solution representation in a pure absorber medium in 1-D slab geometry. Further, there is no inherent limit to local truncation error order of accuracy when using interpolatory points that correspond to the quadrature points of high order accuracy numerical quadrature schemes. (authors)
Efficient large-scale finite-element computations in a CRAY environment
Goudreau, G.L.; Bailey, R.A.; Hallquist, J.O.; Murray, R.C.; Sackett, S.J.
1983-06-01
The Lawrence Livermore National Laboratory engineering computational experience on the CRAY-1 is highlighted in the context of our large general purpose solid and structural mechanics codes. DYNA2D and DYNA3D are explicit large deformation inelastic Lagrangian codes with one point elements and hourglass control. NIKE2D and NIKE3D are implicit codes of comparable continuum formulation but use two point constant pressure elements and an optimized linear equation solver. NIKE3D has a finite rotation plastic resultant shell element. The new general purpose linear elastic structures code GEMINI is also illustrated for large static and eigenvalue analysis. 19 references.
Lin, T.; Wang, H.
1995-12-31
The swift improvement of computational capabilities enables us to apply finite element methods to simulate more and more problems arising from various applications. A fundamental question associated with finite element simulations is their accuracy. In other words, before we can make any decisions based on the numerical solutions, we must be sure that they are acceptable in the sense that their errors are within the given tolerances. Various estimators have been developed to assess the accuracy of finite element solutions, and they can be classified basically into two types: a priori error estimates and a posteriori error estimates. While a priori error estimates can give us asymptotic convergence rates of numerical solutions in terms of the grid size before the computations, they depend on certain Sobolev norms of the true solutions which are not known, in general. Therefore, it is difficult, if not impossible, to use a priori estimates directly to decide whether a numerical solution is acceptable or a finer partition (and so a new numerical solution) is needed. In contrast, a posteriori error estimates depends only on the numerical solutions, and they usually give computable quantities about the accuracy of the numerical solutions.
A knowledge-based approach to the adaptive finite element analysis
Haghighi, K.; Kang, E.
1995-12-31
An automatic and knowledge-based finite element mesh generator (INTELMESH), which makes extensive use of interactive computer graphics techniques, has been developed. INTELMESH is designed for planar domains and axisymmetric 3-D structures of elasticity and heat transfer subjected to mechanical and thermal loading. It intelligently identifies the critical regions/points in the problem domain and utilizes the new concepts of substructuring and wave propagation to choose the proper mesh size for them. INTELMESH generates well-shaped triangular elements by applying triangulation and Laplacian smoothing procedures. The adaptive analysis involves the initial finite element analysis and an efficient a-posteriori error analysis and estimation. Once a problem is defined, the system automatically builds a finite element model and analyzes the problem through an automatic iterative process until the error reaches a desired level. It has been shown that the proposed approach which initiates the process with an a-priori, and near optimum mesh of the object, converges to the desired accuracy in less time and at less cost.
Method of lightening radiation darkened optical elements
Reich, Frederich R.; Schwankoff, Albert R.
1980-01-01
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.
METHOD FOR MAKING FUEL ELEMENTS
Kates, L.W.; Campbell, R.W.; Heartel, R.H.W.
1960-08-01
A method is given for making zirconium-clad uranium wire. A tube of zirconium is closed with a zirconium plug, after which a chilled uranium core is inserted in the tube to rest against the plug. Additional plugs and cores are inserted alternately as desired. The assembly is then sheathed with iron, hot worked to the desired size, and the iron sheath removed.
Finite element strategies for the efficient analysis and evaluation of solar collector structures
Koteras, J. R.
1980-01-01
Concentrating or reflecting structures for solar energy systems must be evaluated as to their structural integrity and optical performance. Computer studies can be used as an integral part of these evaluations. The computer studies make use of finite element structural codes coupled with post-processors that calculate optical data. If the analysis of a solar structure is to be carried out in an efficient manner, these computer codes must have certain capabilities. A number of solar energy projects at Sandia National Laboratories have made extensive use of finite element analyses. The analyses have been useful in evaluating design concepts which hold promise for large scale use in solar energy projects. Analysis procedures have been developed for some structures so that evaluations can be carried out in a straightforward manner.
Gen Purpose 1-D Finite Element Network Fluid Flow Heat Transfer System Simulator
Energy Science and Technology Software Center (OSTI)
1993-08-02
SAFSIM (System Analysis Flow Simulator) is a FORTRAN computer program to simulate the integrated performance of systems involving fluid mechanics, heat transfer, and reactor dynamics. SAFSIM provides sufficient versatility to allow the engineering simulation of almost any system, from a backyard sprinkler system to a clustered nuclear reactor propulsion system. In addition to versatility, speed and robustness are primary SAFSIM development goals. SAFSIM contains three basic physics modules: (1) a one-dimensional finite element fluid mechanicsmoreÂ Â» module with multiple flow network capability; (2) a one-dimensional finite element structure heat transfer module with multiple convection and radiation exchange capability; and (3) a point reactor dynamics module with reactivity feedback and decay heat capability. SAFSIM can be used for compressible and incompressible, single-phase, multicomponent flow systems.Â«Â less
FEHMN 1.0: Finite element heat and mass transfer code
Zyvoloski, G.; Dash, Z.; Kelkar, S.
1991-04-01
A computer code is described which can simulate non-isothermal multiphase multicomponent flow in porous media. It is applicable to natural-state studies of geothermal systems and ground-water flow. The equations of heat and mass transfer for multiphase flow in porous and permeable media are solved using the finite element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat and mass transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. A summary of the equations in the model and the numerical solution procedure are provided in this report. A user`s guide and sample problems are also included. The main use of FEHMN will be to assist in the understanding of flow fields in the saturated zone below the proposed Yucca Mountain Repository. 33 refs., 27 figs., 12 tabs.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Barabash, R. I.; Agarwal, V.; Koric, S.; Jasiuk, I.; Tischler, J. Z.
2016-01-01
TmoreÂ Â» he depth-dependent strain partitioning across the interfaces in the growth direction of the NiAl/Cr(Mo) nanocomposite between the Cr and NiAl lamellae was directly measured experimentally and simulated using a finite element method (FEM). Depth-resolved X-ray microdiffraction demonstrated that in the as-grown state both Cr and NiAl lamellae grow along the 111 direction with the formation of as-grown distinct residual ~0.16% compressive strains for Cr lamellae and ~0.05% tensile strains for NiAl lamellae.hree-dimensional simulations were carried out using an implicit FEM. First simulation was designed to study residual strains in the composite due to cooling resulting in formation of crystals. Strains in the growth direction were computed and compared to those obtained from the microdiffraction experiments. Second simulation was conducted to understand the combined strains resulting from cooling and mechanical indentation of the composite. Numerical results in the growth direction of crystal were compared to experimental results confirming the experimentally observed trends.Â«Â less
Spectral analysis method for detecting an element
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-12
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.
Wakefield Computations for the CLIC PETS using the Parallel Finite Element Time-Domain Code T3P
Candel, A; Kabel, A.; Lee, L.; Li, Z.; Ng, C.; Schussman, G.; Ko, K.; Syratchev, I.; /CERN
2009-06-19
In recent years, SLAC's Advanced Computations Department (ACD) has developed the high-performance parallel 3D electromagnetic time-domain code, T3P, for simulations of wakefields and transients in complex accelerator structures. T3P is based on advanced higher-order Finite Element methods on unstructured grids with quadratic surface approximation. Optimized for large-scale parallel processing on leadership supercomputing facilities, T3P allows simulations of realistic 3D structures with unprecedented accuracy, aiding the design of the next generation of accelerator facilities. Applications to the Compact Linear Collider (CLIC) Power Extraction and Transfer Structure (PETS) are presented.
Finite element simulation of magnetic detection of creep damage at seam welds
Sablik, M.J.; Rubin, S.W.; Jiles, D.C.; Kaminski, D.; Bi, Y.
1996-09-01
Using appropriately modified magnetization curves for each element of creep-damaged material, a finite element calculation has been carried out to simulate magnetic detection of non-uniform creep damage around a seam weld in a 2.25 Cr 1 Mo steam pipe. The magnetization curves for the creep-damaged elements were obtained from an earlier model for the magnetic effect of a uniformly creep-damaged material as given by Chen, et al. In the finite element calculation, a magnetic C-core with primary and secondary coils was placed with its pole pieces flush against the specimen in the vicinity of the weld. The secondary emf was shown to be reduced when creep damage was present inside the pipe wall at the cusp of weld and in the vicinity of the cusp. The C-core detected the creep damage best if it completely spanned the weld seam width. Also, the current in the primary needed to be such that the C-core was not magnetically saturated.
Laser optical element mounting arrangement and method
Whitehouse, D.R.
1988-10-11
This patent describes a method of mounting and sealing an optical element of a laser to a mounting support member of the laser comprising the steps of providing an optical element mounting surface the the support member with a smooth surface which is sufficiently smooth to form a vacuum seal with a surface of the optical element, the vacuum seal preventing at least to a large and acceptable degree, movement of gases through the seal between the optical element and the mounting support member as a result of a gas pressure differential across the seal during operation of the laser and contacting the surface of the optical element with the mounting surface of the support member, at ambient temperature to form the vacuum seal between the surfaces without the use of a sealant between the surfaces, and therein the vacuum seal between the surfaces can be broken and remade without use of a sealant.
Method and apparatus for staking optical elements
Woods, Robert O.
1988-01-01
A method and apparatus for staking two optical elements together in order to retain their alignment is disclosed. The apparatus includes a removable adaptor made up of first and second adaptor bodies each having a lateral slot in their front and side faces. The adaptor also includes a system for releasably attaching each adaptor body to a respective optical element such that when the two optical elements are positioned relative to one another the adaptor bodies are adjacent and the lateral slots therein are aligned to form key slots. The adaptor includes keys which are adapted to fit into the key slots. A curable filler material is employed to retain the keys in the key slots and thereby join the first and second adaptor bodies to form the adaptor. Also disclosed is a method for staking together two optical elements employing the adaptor of the present invention.
Method and apparatus for staking optical elements
Woods, Robert O.
1988-10-04
A method and apparatus for staking two optical elements together in order to retain their alignment is disclosed. The apparatus includes a removable adaptor made up of first and second adaptor bodies each having a lateral slot in their front and side faces. The adaptor also includes a system for releasably attaching each adaptor body to a respective optical element such that when the two optical elements are positioned relative to one another the adaptor bodies are adjacent and the lateral slots therein are aligned to form key slots. The adaptor includes keys which are adapted to fit into the key slots. A curable filler material is employed to retain the keys in the key slots and thereby join the first and second adaptor bodies to form the adaptor. Also disclosed is a method for staking together two optical elements employing the adaptor of the present invention.
Gartling, D.K.
1996-05-01
User instructions are given for the finite element, electromagnetics program, TORO II. The theoretical background and numerical methods used in the program are documented in SAND95-2472. The present document also describes a number of example problems that have been analyzed with the code and provides sample input files for typical simulations. 20 refs., 34 figs., 3 tabs.
Regueiro, Richard A.; Borja, R. I.; Foster, C. D.
2006-10-01
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
Improving the selection of bolted joint parameters by using results from finite-element analysis
Ziada, H.H.
1992-10-01
Finite-element techniques were used to determine the combination of parameters that would produce the optimum design for the bolted joint. A finite-element model was constructed from a combination of two-dimensional isoparametric axisymmetric solid elements and gap elements. The gap elements represented the physical contact between the plate and bolt head and the plate and nut. The model was used to determine the stress patterns, displacements, and contact areas (separation diameters) across the bolted joint. Computed results for the contact area showed the following: the contact area is determined by joint design, not bolt Toad magnitude; joint thickness, plate thickness ratio, and bolt-head dimensions have a pronounced effect on the contact area; and the maximum contact area occurs with plates of equal thickness. The change in separation diameters (D) can be presented in a simple straight-line expression D/d = 1 + 1.1 (L/d), where L = joint thickness and d = bolt diameter. The stress results indicated that neither load nor stress distribution under the bolt head or nut is constant or uniform. Bolt-head dimensions, joint thickness, and plate thickness ratio markedly affect the developed stresses, with the maximum stress value occurring at the bolt-plate interface. Beyond the bolt head, maximum stress develops at the interface plane of the bolted plates, with that stress decreasing as the plate thickness ratio increases reaching its lowest value at the equal-thickness ratio. These results indicate that design stiffness formulas used heretofore are not adequate. Results provided by this analysis should prove more effective in achieving optimum bolt design for bolted joints.
Shephard, M.S.; Dey, S.; Georges, M.K.
1995-12-31
Specific issues associated with the automatic generation of finite element meshes for curved geometric domains axe considered. A review of the definition of when a triangulation is a valid mesh, a geometric triangulation, for curved geometric domains is given. Consideration is then given to the additional operations necessary to maintain the validity of a mesh when curved finite elements are employed. A procedure to control the mesh gradations based on the curvature of the geometric model faces is also given.
Three-dimensional finite element impact analysis of a nuclear waste truck cask
Miller, J.D.
1985-05-01
A three-dimensional finite element impact analysis of a hypothetical accident event for the preliminary design of a shipping cask to be used to transport radioactive waste by standard tractor-semitrailer truck is presented. The dynamic structural analysis code DYNA3D, run on Sandia's Cray-1 computer, was used to calculate the effects of the closure-end of the cask impacting a rigid, frictionless surface on an edge of its external impact limiter after a 30-foot fall. The center of gravity of the 304 stainless steel and depleted uranium cask was assumed to be directly above the impact point. An elastic-plastic material constitutive model was used to calculate the nonlinear response of the cask components to the transient loading. Results from the calculations show the cask sustained large localized deformations. However, these were almost entirely confined to the impact limiters built into the cask. The closure sections were determined to remain intact and leakage would not be expected after the event. Interactive color computer graphics were used throughout the analysis, proving to be extremely helpful for generation and verification of the geometry and boundary conditions of the finite element model and for interpretation of the analysis results. 12 refs., 29 figs., 4 tabs.
Finite element thermomechanical simulation of droplets impacting on a rigid substrate
Bertagnolli, M.; Marchese, M.; Jacucci, G.; St. Doltsinis, I.; Noelting, S.
1994-12-31
The plasma spray process is a convenient way to coat a piece of material with a layer of another material, to protect the first from thermal shock or environmental degradation. Finite Element simulation techniques (FEM) for the spreading process of a ceramic liquid droplet impacting on a flat cold surface have been developed. The goal of the present investigation is (1) to predict the geometrical form of the splat as a function of process parameters, such as initial temperature and velocity, and (2) to follow the thermal field developing in the droplet up to solidification. A non-linear finite element procedure has been extended in order to model the complex physical phenomena involved in the impact process. The dynamic motion of the viscous melt in the drops as constrained by elastic surface tensions in interaction with the developing contact with the target, ultimately has been coupled to transient thermal phenomena accounting also for the solidification of the material. In a first model description, spherical particles of liquid ceramic of given temperature and velocity impact on a flat, cool rigid surface. The deformation of the splat geometry as well as the evolution of the thermal field within the splat are followed up to the final state and require adaptive discretization techniques. The authors discuss an utilization of the proposed model in correlating flattening degrees with the initial process parameters.
Efficient Computation of Info-Gap Robustness for Finite Element Models
Stull, Christopher J.; Hemez, Francois M.; Williams, Brian J.
2012-07-05
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.
Soulami, Ayoub; Lavender, Curt A.; Paxton, Dean M.; Burkes, Douglas
2014-04-23
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).
Microscopic description of fission dynamics: finite element method...
Office of Scientific and Technical Information (OSTI)
Report Number(s): LLNL-PROC-678472 DOE Contract Number: AC52-07NA27344 Resource Type: Conference Resource Relation: Conference: Presented at: Wonder 2015: Fourth International ...
An Implementation of the Generalized Finite Element Method for...
Office of Scientific and Technical Information (OSTI)
Abstract not provided. Authors: Robbins, Joshua ; Voth, Thomas E. Publication Date: 2009-08-01 OSTI Identifier: 1142075 Report Number(s): SAND2009-5068C 507910 DOE Contract Number: ...
Stochastic finite element methods for partial differential equations...
Office of Scientific and Technical Information (OSTI)
DOI: 10.1017S0962492914000075 Select the DOI to obtain a copy of this journal article from the publisher. Find in Google Scholar Find in Google Scholar Search WorldCat Search ...
METHOD OF MAKING WIRE FUEL ELEMENTS
Zambrow, J.L.
1960-08-01
A method is given for making a nuclear reactor fuel element in the form of a uranium-bearing wire clad with zirconium. A uranium bar is enclosed in a zirconium sheath which is coated with an oxide of magnesium, beryllium, or zirconium. The sheathed bar is then placed in a steel tube and reduced to the desired diameter by swaging at 800 to 900 deg C, after which the steel and oxide are removed.
Some aspects of sheet forming simulation using explicit finite element techniques
Whirley, R.G.; Engelmann, B.E.; Logan, R.W.
1992-07-28
This paper discusses some aspects of the application of the explicit nonlinear finite element code DYNA3D to the numerical simulation of sheet forming processes. Although DYNA3D is intrinisically a code for transient dynamic response, it has been successfully applied to the quasistatic sheet forming problem by slow application of the loads. Although somewhat nontraditional, this approach has proven very robust in application to large-scale problems. This paper discusses the components of a successful sheet forming analysis, including constitutive modeling, contact definition, guidelines for load application, and some suggestions for solution evaluation. Further, this paper illustrates the use of DYNA3D to compute the Limiting Draw Ratio as a function of material anisotropy and compares the computed results for 3-D forming problems with experimental data.
GRIZ: Finite element analysis results visualization for unstructured grids. User manual
Dovey, D.J.; Spelce, T.E.
1993-10-01
GRIZ supports interactive visualization of finite element analysis results on unstructured grids. GRIZ is a general-purpose post-processing application which is designed to work with a variety of an analysis codes. Currently, GRIZ is capable of calculating and displaying derived variables for the DYNA3D, NIKE3D and TOPAZ3D analysis codes. GRIZ reads in data files in the ``MDG plotfile`` format. GRIZ provides support for modern 3D visualization techniques such as isosurface display, cutting planes and display of vector data. GRIZ also incorporates the ability to animate data over time and to store animation frames to a video disk. GRIZ is designed to utilize the capabilities of modern graphics workstations which provide hardware support for 3D graphics, thereby giving the user as much interactive performance as possible. This should make it easier for analysts to explore and interrogate their analysis results.
dgtoexo2: A Distorted Grid Output File to Exodus II Finite Element Database Conversion Utility
Moffat, H.K.
1998-12-01
This report describes how to obtain publication-quality graphics from distorted grid electronic structure codes using the combination of the conversion utility, dgtoexo2, and mustafa, an AVS Express application. dgtoexo2 converts scalar function results from a format applicable to distorted grid codes into the Exodus II unstructured finite element data representation. nmstafa can read Exodus II files and use the AVS Express engine to visualize data on unix and Windows NT platforms. Though not designed for the purpose, the dgtoexo2/EXOdUS II/mustafa combination is sufficiently versatile to provide for the specialized graphics needs of electronic structure codes. The combination also scales well, producing robust performance for problems involving millions of grid points.
3-D Finite Element Analysis of Induction Logging in a Dipping Formation
EVERETT,MARK E.; BADEA,EUGENE A.; SHEN,LIANG C.; MERCHANT,GULAMABBAS A.; WEISS,CHESTER J.
2000-07-20
Electromagnetic induction by a magnetic dipole located above a dipping interface is of relevance to the petroleum well-logging industry. The problem is fully three-dimensional (3-D) when formulated as above, but reduces to an analytically tractable one-dimensional (1-D) problem when cast as a small tilted coil above a horizontal interface. The two problems are related by a simple coordinate rotation. An examination of the induced eddy currents and the electric charge accumulation at the interface help to explain the inductive and polarization effects commonly observed in induction logs from dipping geological formations. The equivalence between the 1-D and 3-D formulations of the problem enables the validation of a previously published finite element solver for 3-D controlled-source electromagnetic induction.
Finite element residual stress analysis of induction heating bended ferritic steel piping
Kima, Jong Sung; Kim, Kyoung-Soo; Oh, Young-Jin; Chang, Hyung-Young; Park, Heung-Bae
2014-10-06
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.
Finite element analysis of three TVA dams with alkali-aggregate reaction
Grenoble, B.A.; Meisenheimer, J.K.; Wagner, C.D.; Newell, V.A.
1995-12-31
Three large Tennessee Valley Authority (TVA) dams are currently experiencing problems caused by alkali-aggregate reaction (AAR). Since the fall of 1990, engineers in Stone & Webster`s Denver, Colorado office have been working with TVA to evaluate how AAR is affecting the dams and to identify measures for controlling the adverse effects of the concrete growth. This paper provides an overview of how finite element analysis is being used to understand the affects of AAR on these structures and to evaluate alternatives for minimizing the adverse effects of the concrete growth. Work on Hiwassee Dam is essentially complete, while that on the Chickamauga and Fontana Projects is still in progress. Consequently, this paper will focus primarily on Hiwassee Dam. The ongoing work on the other two projects will only be discussed briefly.
White, D; Fasenfest, B; Rieben, R; Stowell, M
2006-09-08
We are concerned with the solution of time-dependent electromagnetic eddy current problems using a finite element formulation on three-dimensional unstructured meshes. We allow for multiple conducting regions, and our goal is to develop an efficient computational method that does not require a computational mesh of the air/vacuum regions. This requires a sophisticated global boundary condition specifying the total fields on the conductor boundaries. We propose a Biot-Savart law based volume-to-surface boundary condition to meet this requirement. This Biot-Savart approach is demonstrated to be very accurate. In addition, this approach can be accelerated via a low-rank QR approximation of the discretized Biot-Savart law.
Mostafa, Salwa; Lavrik, Nickolay V; Bannuru, Thirumalesh; Rajic, Slobodan; Islam, Syed K; Datskos, Panos G; Hunter, Scott Robert
2011-01-01
A self resonating bimorph cantilever structure for fast temperature cycling in a pyroelectric energy harvester has been modeled using a finite element method. The effect of constituting material properties and system parameters on the frequency and magnitude of temperature cycling and the efficiency of energy recycling using the proposed structure has been investigated. Results show that thermal contact conductance and heat source temperature play a key role in dominating the cycling frequency and efficiency of energy recycling. An optimal solution for the most efficient energy scavenging process has been sought by studying the performance trend with different variable parameters such as thermal contact conductance, heat source temperature, device aspect ratio and constituent materials of varying thermal conductivity and expansion coefficients.
Finite element analysis of an extended end-plate connection using the T-stub approach
Muresan, Ioana Cristina; Balc, Roxana
2015-03-10
Beam-to-column end-plate bolted connections are usually used as moment-resistant connections in steel framed structures. For this joint type, the deformability is governed by the deformation capacity of the column flange and end-plate under tension and elongation of the bolts. All these elements around the beam tension flange form the tension region of the joint, which can be modeled by means of equivalent T-stubs. In this paper a beam-to-column end-plate bolted connection is substituted with a T-stub of appropriate effective length and it is analyzed using the commercially available finite element software ABAQUS. The performance of the model is validated by comparing the behavior of the T-stub from the numerical simulation with the behavior of the connection as a whole. The moment-rotation curve of the T-stub obtained from the numerical simulation is compared with the behavior of the whole extended end-plate connection, obtained by numerical simulation, experimental tests and analytical approach.
Fluorescence photon migration by the boundary element method
Fedele, Francesco; Eppstein, Margaret J. . E-mail: maggie.eppstein@uvm.edu; Laible, Jeffrey P.; Godavarty, Anuradha; Sevick-Muraca, Eva M.
2005-11-20
The use of the boundary element method (BEM) is explored as an alternative to the finite element method (FEM) solution methodology for the elliptic equations used to model the generation and transport of fluorescent light in highly scattering media, without the need for an internal volume mesh. The method is appropriate for domains where it is reasonable to assume the fluorescent properties are regionally homogeneous, such as when using highly specific molecularly targeted fluorescent contrast agents in biological tissues. In comparison to analytical results on a homogeneous sphere, BEM predictions of complex emission fluence are shown to be more accurate and stable than those of the FEM. Emission fluence predictions made with the BEM using a 708-node mesh, with roughly double the inter-node spacing of boundary nodes as in a 6956-node FEM mesh, match experimental frequency-domain fluorescence emission measurements acquired on a 1087 cm{sup 3} breast-mimicking phantom at least as well as those of the FEM, but require only 1/8 to 1/2 the computation time.
METHOD OF PREPARING A CERAMIC FUEL ELEMENT
Ross, W.T.; Bloomster, C.H.; Bardsley, R.E.
1963-09-01
A method is described for preparing a fuel element from -325 mesh PuO/ sub 2/ and -20 mesh UO/sub 2/, and the steps of screening --325 mesh UO/sub 2/ from the -20 mesh UO/sub 2/, mixing PuO/sub 2/ with the --325 mesh UO/sub 2/, blending this mixture with sufficient --20 mesh UO/sub 2/ to obtain the desired composition, introducing the blend into a metal tube, repeating the procedure until the tube is full, and vibrating the tube to compact the powder are included. (AEC)
Method for locating defective nuclear fuel elements
Lawrie, W.E.; White, N.W.; Womack, R.E.
1982-02-02
Defects in nuclear fuel elements are ascertained and located within an assembled fuel assembly by ultrasonic means. In a typical embodiment of the invention, an ultrasonic search unit is positioned within the fuel assembly opposite the lower plenum of the fuel element to be tested. An ultrasonic pulse is radially projected into the element. Defective fuel elements are ascertained by ultrasonic reflection measurements.
Method of monitoring stored nuclear fuel elements
Borloo, E.; Buergers, W.; Crutzen, S.; Vinche, C.
1983-05-24
To monitor a nuclear fuel element or fuel elements located in a store, e.g. a pond in a swimming pool reactor, the store is illuminated ultrasonically using one or more transducers transmitting ultrasonic signals in one or more predetermined directions to obtain an output which, because it depends on the number and relative location of the fuel elements in the store, and the structure of the store itself is distinctive to the fuel elements or elements stored therein. From this distinctive output is derived an identity unique to the stored fuel element or elements and a reference signal indicative of the whole structure when intact, the reference signal and identity being recorded. Subsequent ultrasonic testing of the store and its contents under identical operating conditions produces a signal which is compared to the recorded reference signal and if different therefrom reveals the occurrence of tampering with the store and/or the fuel element or elements.
Soulami, Ayoub; Choi, Kyoo Sil; Liu, Wenning N.; Sun, Xin; Khaleel, Mohammad A.
2009-04-01
Recently, several studies conducted by automotive industry revealed the tremendous advantages of Advanced High Strength Steels (AHSS). TRansformation Induced Plasticity (TRIP) steel is one of the typical representative of AHSS. This kind of materials exhibits high strength as well as high formability. Analyzing the crack behaviour in TRIP steels is a challenging task due to the microstructure level inhomogeneities between the different phases (Ferrite, Bainite, Austenite, Martensite) that constitute these materials. This paper aims at investigating the fracture resistance of TRIP steels. For this purpose, a micromechanical finite element model is developed based on the actual microstructure of a TRIP 800 steel. Uniaxial tensile tests on TRIP 800 sheet notched specimens were also conducted and tensile properties and R-curves (Resistance curves) were determined. The comparison between simulation and experimental results leads us to the conclusion that the method using microstructure-based representative volume element (RVE) captures well enough the complex behavior of TRIP steels. The effect of phase transformation, which occurs during the deformation process, on the toughness is observed and discussed.
Electrical and thermal finite element modeling of arc faults in photovoltaic bypass diodes.
Bower, Ward Isaac; Quintana, Michael A.; Johnson, Jay
2012-01-01
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.
Christian, Joshua M.; Ho, Clifford Kuofei
2010-04-01
Understanding the effects of gravity and wind loads on concentrating solar power (CSP) collectors is critical for performance calculations and developing more accurate alignment procedures and techniques. This paper presents a rigorous finite-element model of a parabolic trough collector that is used to determine the impact of gravity loads on bending and displacements of the mirror facets and support structure. The geometry of the LUZ LS-2 parabolic trough collector was modeled using SolidWorks, and gravity-induced loading and displacements were simulated in SolidWorks Simulation. The model of the trough collector was evaluated in two positions: the 90{sup o} position (mirrors facing upward) and the 0{sup o} position (mirrors facing horizontally). The slope errors of the mirror facet reflective surfaces were found by evaluating simulated angular displacements of node-connected segments along the mirror surface. The ideal (undeformed) shape of the mirror was compared to the shape of the deformed mirror after gravity loading. Also, slope errors were obtained by comparing the deformed shapes between the 90{sup o} and 0{sup o} positions. The slope errors resulting from comparison between the deformed vs. undeformed shape were as high as {approx}2 mrad, depending on the location of the mirror facet on the collector. The slope errors resulting from a change in orientation of the trough from the 90{sup o} position to the 0{sup o} position with gravity loading were as high as {approx}3 mrad, depending on the location of the facet.
MPSalsa a finite element computer program for reacting flow problems. Part 2 - user`s guide
Salinger, A.; Devine, K.; Hennigan, G.; Moffat, H.
1996-09-01
This manual describes the use of MPSalsa, an unstructured finite element (FE) code for solving chemically reacting flow problems on massively parallel computers. MPSalsa has been written to enable the rigorous modeling of the complex geometry and physics found in engineering systems that exhibit coupled fluid flow, heat transfer, mass transfer, and detailed reactions. In addition, considerable effort has been made to ensure that the code makes efficient use of the computational resources of massively parallel (MP), distributed memory architectures in a way that is nearly transparent to the user. The result is the ability to simultaneously model both three-dimensional geometries and flow as well as detailed reaction chemistry in a timely manner on MT computers, an ability we believe to be unique. MPSalsa has been designed to allow the experienced researcher considerable flexibility in modeling a system. Any combination of the momentum equations, energy balance, and an arbitrary number of species mass balances can be solved. The physical and transport properties can be specified as constants, as functions, or taken from the Chemkin library and associated database. Any of the standard set of boundary conditions and source terms can be adapted by writing user functions, for which templates and examples exist.
Efficient multilevel finite-element approach to three-dimensional phase-change problems
Lee, R.T.; Chiou, W.Y.
1997-01-01
A finite-element (FE) formulation suitable for a multigrid algorithm in solving three-dimensional phase-change problems is described. This formulation is based on the averaged specific heat model. The algorithm has been proved to be very useful for large problems where the computational complexity can be reduced from O(n{sup 3}) to O(n ln n) with high storage efficiency in a personal computer. To evaluate the accuracy of the present algorithm, the numerical results for larger slender ratio are compared with previous analytical solutions. Results show that the numerical solutions at the symmetric surface of the long axis are in very good agreement with the two-dimensional exact solutions for slender ratio = 5. The magnitudes of time steps and freezing-temperature intervals are insensitive to the maximal and average absolute errors when the time step is less than 0.01 s. Consequently, a larger time step can be used to save computing time and retain the same order of accuracy. This algorithm is also available for pure metals and alloys that exhibit a very large or small (or zero) freezing-temperature interval.
Sondak, D.; Shadid, J. N.; Oberai, A. A.; Pawlowski, R. P.; Cyr, E. C.; Smith, T. M.
2015-04-29
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and high Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Sondak, D.; Shadid, J. N.; Oberai, A. A.; Pawlowski, R. P.; Cyr, E. C.; Smith, T. M.
2015-04-29
New large eddy simulation (LES) turbulence models for incompressible magnetohydrodynamics (MHD) derived from the variational multiscale (VMS) formulation for finite element simulations are introduced. The new models include the variational multiscale formulation, a residual-based eddy viscosity model, and a mixed model that combines both of these component models. Each model contains terms that are proportional to the residual of the incompressible MHD equations and is therefore numerically consistent. Moreover, each model is also dynamic, in that its effect vanishes when this residual is small. The new models are tested on the decaying MHD Taylor Green vortex at low and highmoreÂ Â» Reynolds numbers. The evaluation of the models is based on comparisons with available data from direct numerical simulations (DNS) of the time evolution of energies as well as energy spectra at various discrete times. Thus a numerical study, on a sequence of meshes, is presented that demonstrates that the large eddy simulation approaches the DNS solution for these quantities with spatial mesh refinement.Â«Â less
A finite integration method for conformal, structured-grid, electromagnetic simulation
Cooke, S.J. . E-mail: simon.cooke@nrl.navy.mil; Shtokhamer, R.; Mondelli, A.A.; Levush, B.
2006-06-10
We describe a numerical scheme for solving Maxwell's equations in the frequency domain on a conformal, structured, non-orthogonal, multi-block mesh. By considering Maxwell's equations in a volume parameterized by dimensionless curvilinear coordinates, we obtain a set of tensor equations that are a continuum analogue of common circuit equations, and that separate the metrical and metric-free parts of Maxwell's equations and the material constitutive relations. We discretize these equations using a new formulation that treats the electric field and magnetic induction using simple basis-function representations to obtain a discrete form of Faraday's law of induction, but that uses finite integral representations for the displacement current and magnetic field to obtain a discrete form of Ampere's law, as in the finite integration technique [T. Weiland, A discretization method for the solution of Maxwell's equations for six-component fields, Electron. Commun. (AE U) 31 (1977) 116; T. Weiland, Time domain electromagnetic field computation with finite difference methods, Int. J. Numer. Model: Electron. Netw. Dev. Field 9 (1996) 295-319]. We thereby derive new projection operators for the discrete tensor material equations and obtain a compact numerical scheme for the discrete differential operators. This scheme is shown to exhibit significantly reduced numerical dispersion when compared to the standard linear finite element method. We take advantage of the mesh structure on a block-by-block basis to implement these numerical operators efficiently, and achieve computational speed with modest memory requirements when compared to explicit sparse matrix storage. Using the Jacobi-Davidson [G.L.G. Sleijpen, H.A. van der Vorst, A Jacobi-Davidson iteration method for linear eigenvalue problems, SIAM J. Matrix Anal. Appl. 17 (2) (1996) 401-425; S.J. Cooke, B. Levush, Eigenmode solution of 2-D and 3-D electromagnetic cavities containing absorbing materials using the Jacobi
Iskander, S.K.
1981-02-01
Two finite element (FE) approaches were used to calculate opening mode I stress intensity factors (K/sub I/) in two- or three-dimensional (2-D and 3-D) problems for the Heavy-Section Steel Technology (HSST) program. For problems that can be modeled in two dimensions, two techniques were used. One of these may be termed an ''energy release rate'' technique, and the other is based on the classical near-tip displacement and stress field equations. For three-dimensional problems, only the latter technique was used. In the energy release technique, K/sub I/ is calculated as the change in potential energy of the structure due to a small change in crack length. The potential energy is calculated by the FE method but without completely solving the system of linear equations for the displacements. Furthermore, the system of linear equations is only slightly perturbed by the change in crack length and, therefore, many computations need not be repeated for the second structure with the slight change in crack length. Implementation of these last two items has resulted in considerable savings in the calculation of K/sub I/ as compared to two complete FE analyses. These ideas are incorporated in the FMECH code. The accuracy of the methods has been checked by comparing the results of the two approaches with each other and with closed form solutions. It is estimated that the accuracy of the results is about +-5%.
Lehtikangas, O.; Tarvainen, T.; Kim, A.D.; Arridge, S.R.
2015-02-01
The radiative transport equation can be used as a light transport model in a medium with scattering particles, such as biological tissues. In the radiative transport equation, the refractive index is assumed to be constant within the medium. However, in biomedical media, changes in the refractive index can occur between different tissue types. In this work, light propagation in a medium with piece-wise constant refractive index is considered. Light propagation in each sub-domain with a constant refractive index is modeled using the radiative transport equation and the equations are coupled using boundary conditions describing Fresnel reflection and refraction phenomena on the interfaces between the sub-domains. The resulting coupled system of radiative transport equations is numerically solved using a finite element method. The approach is tested with simulations. The results show that this coupled system describes light propagation accurately through comparison with the Monte Carlo method. It is also shown that neglecting the internal changes of the refractive index can lead to erroneous boundary measurements of scattered light.
MODELING OF HIGH SPEED FRICTION STIR SPOT WELDING USING A LAGRANGIAN FINITE ELEMENT APPROACH
Miles, Michael; Karki, U.; Woodward, C.; Hovanski, Yuri
2013-09-03
Friction stir spot welding (FSSW) has been shown to be capable of joining steels of very high strength, while also being very flexible in terms of controlling the heat of welding and the resulting microstructure of the joint. This makes FSSW a potential alternative to resistance spot welding (RSW) if tool life is sufficiently high, and if machine spindle loads are sufficiently low so that the process can be implemented on an industrial robot. Robots for spot welding can typically sustain vertical loads of about 8kN, but FSSW at tool speeds of less than 3000 rpm cause loads that are too high, in the range of 11-14 kN. Therefore, in the current work tool speeds of 3000 rpm and higher were employed, in order to generate heat more quickly and to reduce welding loads to acceptable levels. The FSSW process was modeled using a finite element approach with the ForgeÂ® software package. An updated Lagrangian scheme with explicit time integration was employed to model the flow of the sheet material, subjected to boundary conditions of a rotating tool and a fixed backing plate [3]. The modeling approach can be described as two-dimensional, axisymmetric, but with an aspect of three dimensions in terms of thermal boundary conditions. Material flow was calculated from a velocity field which was two dimensional, but heat generated by friction was computed using a virtual rotational velocity component from the tool surface. An isotropic, viscoplastic Norton-Hoff law was used to model the evolution of material flow stress as a function of strain, strain rate, and temperature. The model predicted welding temperatures and the movement of the joint interface with reasonable accuracy for the welding of a dual phase 980 steel.
Method and system for processing optical elements using magnetorheological finishing
Menapace, Joseph Arthur; Schaffers, Kathleen Irene; Bayramian, Andrew James; Molander, William A
2012-09-18
A method of finishing an optical element includes mounting the optical element in an optical mount having a plurality of fiducials overlapping with the optical element and obtaining a first metrology map for the optical element and the plurality of fiducials. The method also includes obtaining a second metrology map for the optical element without the plurality of fiducials, forming a difference map between the first metrology map and the second metrology map, and aligning the first metrology map and the second metrology map. The method further includes placing mathematical fiducials onto the second metrology map using the difference map to form a third metrology map and associating the third metrology map to the optical element. Moreover, the method includes mounting the optical element in the fixture in an MRF tool, positioning the optical element in the fixture; removing the plurality of fiducials, and finishing the optical element.
User's Manual for the FEHM Application-A Finite-Element Heat- and Mass-Transfer Code
George A. Zyvoloski; Bruce A. Robinson; Zora V. Dash; Lynn L. Trease
1997-07-07
This document is a manual for the use of the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multicomponent flow in porous media. 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 equations of heat and mass transfer for multiphase flow in porous and permeable media are solved in the FEHM application by using the finite-element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat- and mass-transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. In fact, FEHM is capable of describing flow that is dominated in many areas by fracture and fault flow, including the inherently three-dimensional flow that results from permeation to and from faults and fractures. The code can handle coupled heat and mass-transfer effects, such as boiling, dryout, and condensation that can occur in the near-field region surrounding the potential repository and the natural convection that occurs through Yucca Mountain due to seasonal temperature changes. The code is also capable of incorporating the various adsorption mechanisms, ranging from simple linear relations to nonlinear isotherms, needed to describe the very complex transport processes at Yucca Mountain. This report outlines the uses and capabilities of the FEHM application, initialization of code variables, restart procedures, and error processing. The report describes all the data files, the input data
Validation of Finite-Element Models of Persistent-Current Effects in Nb3Sn Accelerator Magnets
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Wang, X.; Ambrosio, G.; Chlachidze, G.; Collings, E. W.; Dietderich, D. R.; DiMarco, J.; Felice, H.; Ghosh, A. K.; Godeke, A.; Gourlay, S. A.; et al
2015-01-06
Persistent magnetization currents are induced in superconducting filaments during the current ramping in magnets. The resulting perturbation to the design magnetic field leads to field quality degradation, in particular at low field where the effect is stronger relative to the main field. The effects observed in NbTi accelerator magnets were reproduced well with the critical-state model. However, this approach becomes less accurate for the calculation of the persistent-current effects observed in Nb3Sn accelerator magnets. Here a finite-element method based on the measured strand magnetization is validated against three state-of-art Nb3Sn accelerator magnets featuring different subelement diameters, critical currents, magnet designsmoreÂ Â» and measurement temperatures. The temperature dependence of the persistent-current effects is reproduced. Based on the validated model, the impact of conductor design on the persistent current effects is discussed. The performance, limitations and possible improvements of the approach are also discussed.Â«Â less
FUEL ELEMENT AND METHOD OF PREPARATION
Kingston, W.E.
1961-04-25
A nuclear fuel element in the form of a wire is reported. A bar of uranium is enclosed in a thin layer of aluminum and the composite is sheathed in beryllium, zirconium, or stainnless steel. The sheathed article is then drawn to wire form, heated to alloy the aluminum with both uranium and sheath, and finally cold worked.
Optical Method for Automated Real Time Control of Elemental Compositio...
Broader source: All U.S. Department of Energy (DOE) Office Webpages (Extended Search)
Optical Method for Automated Real Time Control of Elemental Composition, Distribution, and Film Thickness in CIGS Solar Cell Production National Renewable Energy Laboratory Contact ...
A non-linear constrained optimization technique for the mimetic finite difference method
Manzini, Gianmarco; Svyatskiy, Daniil; Bertolazzi, Enrico; Frego, Marco
2014-09-30
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.
Woodward, Carol S.; Gardner, David J.; Evans, Katherine J.
2015-01-01
Efficient solutions of global climate models require effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a step size governed by accuracy of the processes of interest rather than by stability of the fastest time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton's method is applied to solve these systems. Each iteration of the Newton's method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but this Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix-vector products. The matrix-vector multiply can also be approximated by a finite difference approximation which may introduce inaccuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite difference approximations of these matrix-vector products for climate dynamics within the spectral element shallow water dynamical core of the Community Atmosphere Model.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Woodward, Carol S.; Gardner, David J.; Evans, Katherine J.
2015-01-01
Efficient solutions of global climate models require effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a step size governed by accuracy of the processes of interest rather than by stability of the fastest time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton's method is applied to solve these systems. Each iteration of the Newton's method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but thismoreÂ Â» Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix-vector products. The matrix-vector multiply can also be approximated by a finite difference approximation which may introduce inaccuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite difference approximations of these matrix-vector products for climate dynamics within the spectral element shallow water dynamical core of the Community Atmosphere Model.Â«Â less
Gardner, David; Woodward, Carol S.; Evans, Katherine J
2015-01-01
Efficient solution of global climate models requires effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a time step dictated by accuracy of the processes of interest rather than by stability governed by the fastest of the time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton s method is applied for these systems. Each iteration of the Newton s method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but this Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix-vector products. The matrix-vector multiply can also be approximated by a finite-difference which may show a loss of accuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite-difference approximations of these matrix-vector products for climate dynamics within the spectral-element based shallow-water dynamical-core of the Community Atmosphere Model (CAM).
User`s manual for the FEHM application -- A finite-element heat- and mass-transfer code
Zyvoloski, G.A.; Robinson, B.A.; Dash, Z.V.; Trease, L.L.
1997-07-01
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 equations of heat and mass transfer for multiphase flow in porous and permeable media are solved in the FEHM application by using the finite-element method. The permeability and porosity of the medium are allowed to depend on pressure and temperature. The code also has provisions for movable air and water phases and noncoupled tracers; that is, tracer solutions that do not affect the heat- and mass-transfer solutions. The tracers can be passive or reactive. The code can simulate two-dimensional, two-dimensional radial, or three-dimensional geometries. In fact, FEHM is capable of describing flow that is dominated in many areas by fracture and fault flow, including the inherently three-dimensional flow that results from permeation to and from faults and fractures. The code can handle coupled heat and mass-transfer effects, such as boiling, dryout, and condensation that can occur in the near-field region surrounding the potential repository and the natural convection that occurs through Yucca Mountain due to seasonal temperature changes. This report outlines the uses and capabilities of the FEHM application, initialization of code variables, restart procedures, and error processing. The report describes all the data files, the input data, including individual input records or parameters, and the various output files. The system interface is described, including the software environment and installation instructions.
Penjweini, R; Zhu, T
2015-06-15
Purpose: The pleural volumes will deform during surgery portion of the pleural photodynamic therapy (PDT) of lung cancer when the pleural cavity is opened. This impact the delivered dose when using highly conformal treatment techniques. In this study, a finite element-based (FEM) deformable image registration is used to quantify the anatomical variation between the contours for the pleural cavities obtained in the operating room and those determined from pre-surgery computed tomography (CT) scans. Methods: An infrared camera-based navigation system (NDI) is used during PDT to track the anatomical changes and contour the lung and chest cavity. A series of CTs of the lungs, in the same patient, are also acquired before the surgery. The structure contour of lung and the CTs are processed and contoured in Matlab and MeshLab. Then, the contours are imported into COMSOL Multiphysics 5.0, where the FEM-based deformable image registration is obtained using the deformed mesh - moving mesh (ALE) model. The NDI acquired lung contour is considered as the reference contour, and the CT contour is used as the target one, which will be deformed. Results: The reconstructed three-dimensional contours from both NDI and CT can be converted to COMSOL so that a three-dimensional ALE model can be developed. The contours can be registered using COMSOL ALE moving mesh model, which takes into account the deformation along x, y and z-axes. The deformed contour has good matches to the reference contour after the dynamic matching process. The resulting 3D deformation map can be used to obtain the locations of other critical anatomic structures, e.g., heart, during surgery. Conclusion: Deformable image registration can fuse images acquired by different modalities. It provides insights into the development of phenomenon and variation in normal anatomical structures over time. The initial assessments of three-dimensional registration show good agreement.
Method of holding optical elements without deformation during their fabrication
Hed, P. Paul
1997-01-01
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.
METHOD OF PREPARING A FUEL ELEMENT
Noland, R.A.; Stone, C.C.
1959-09-01
A method is presented for joining sections of uranium rod which are clad in a corrosion-resistant material of higher melting point than uranium. The method includes the steps of preferentially- etching the uranium in the ends of the sections to be joined to depress the level of the uranium slighily below that of the cladding, bringing the ends to be joined together, applying axial pressure to the joint, and heli-arc welding the joint while rotating the joint, to fuse the cladding to a uniform depth of about 50% of the thickness of the cladding.
Linker, K.L.; Rawlinson, K.S.; Smith, G.
1991-10-01
The Department of Energy's Solar Thermal Program has as one of its program elements the development and evaluation of conversion device technologies applicable to dish-electric systems. The primary research and development combines a conversion device (heat engine), solar receiver, and generator mounted at the focus of a parabolic dish concentrator. The Stirling-cycle heat engine was identified as the conversion device for dish-electric with the most potential for meeting the program's goals for efficiency, reliability, and installed cost. To advance the technology toward commercialization, Sandia National Laboratories has acquired a Stirling Thermal Motors, Inc., kinematic Stirling engine, STM4-120, for evaluation. The engine is being bench-tested at Sandia's Engine Test Facility and will be combined later with a solar receiver for on-sun evaluation. This report presents the engine characteristics, finite element analyses of critical engine components, test system layout, instrumentation, and preliminary performance results from the bench test.
Soulami, Ayoub; Lavender, Curt A.; Paxton, Dean M.; Burkes, Douglas
2015-06-15
Pacific Northwest National Laboratory (PNNL) has been investigating manufacturing processes for the uranium-10% molybdenum alloy plate-type fuel for high-performance research reactors in the United States. This work supports the U.S. Department of Energy National Nuclear Security Administrationâ€™s Office of Material Management and Minimization Reactor Conversion Program. This report documents modeling results of PNNLâ€™s efforts to perform finite-element simulations to predict roll-separating forces for various rolling mill geometries for PNNL, Babcock & Wilcox Co., Y-12 National Security Complex, Los Alamos National Laboratory, and Idaho National Laboratory. The model developed and presented in a previous report has been subjected to further validation study using new sets of experimental data generated from a rolling mill at PNNL. Simulation results of both hot rolling and cold rolling of uranium-10% molybdenum coupons have been compared with experimental results. The model was used to predict roll-separating forces at different temperatures and reductions for five rolling mills within the National Nuclear Security Administration Fuel Fabrication Capability project. This report also presents initial results of a finite-element model microstructure-based approach to study the surface roughness at the interface between zirconium and uranium-10% molybdenum.
METHOD OF PREPARING A FUEL ELEMENT FOR A NUCLEAR REACTOR
Roake, W.E.; Evans, E.A.; Brite, D.W.
1960-06-21
A method of preparing a fuel element for a nuclear reactor is given in which an internally and externally cooled fuel element consisting of two coaxial tubes having a plurality of integral radial ribs extending between the tubes and containing a powdered fuel material is isostatically pressed to form external coolant channels and compact the powder simultaneously.
Methods and devices for fabricating and assembling printable semiconductor elements
Nuzzo, Ralph G.; Rogers, John A.; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao
2009-11-24
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.
Methods and devices for fabricating and assembling printable semiconductor elements
Nuzzo, Ralph G.; Rogers, John A.; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao
2011-07-19
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.
Methods and devices for fabricating and assembling printable semiconductor elements
Nuzzo, Ralph G; Rogers, John A; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao
2014-03-04
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.
Methods and devices for fabricating and assembling printable semiconductor elements
Nuzzo, Ralph G; Rogers, John A; Menard, Etienne; Lee, Keon Jae; Khang, Dahl-Young; Sun, Yugang; Meitl, Matthew; Zhu, Zhengtao
2013-05-14
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.
Hoover, C G; DeGroot, A J; Sherwood, R J
2000-06-01
ParaDyn is a parallel version of the DYNA3D computer program, a three-dimensional explicit finite-element program for analyzing the dynamic response of solids and structures. The ParaDyn program has been used as a production tool for over three years for analyzing problems which range in size from a few tens of thousands of elements to between one-million and ten-million elements. ParaDyn runs on parallel computers provided by the Department of Energy Accelerated Strategic Computing Initiative (ASCI) and the Department of Defense High Performance Computing and Modernization Program. Preprocessing and post-processing software utilities and tools are designed to facilitate the generation of partitioned domains for processors on a massively parallel computer and the visualization of both resultant data and boundary data generated in a parallel simulation. This manual provides a brief overview of the parallel implementation; describes techniques for running the ParaDyn program, tools and utilities; and provides examples of parallel simulations.
Goudreau, G.L.; Hallquist, J.O.
1981-10-01
The state of Lagrangian hydrocodes for computing the large deformation dynamic response of inelastic continuua is reviewed in the context of engineering computation at the Lawrence Livermore National Laboratory, USA, and the DYNA2D/DYNA3D finite elements codes. The emphasis is on efficiency and computational cost. The simplest elements with explicit time integration. The two-dimensional four node quadrilateral and the three-dimensional hexahedron with one point quadrature are advocated as superior to other more expensive choices. Important auxiliary capabilities are a cheap but effective hourglass control, slidelines/planes with void opening/closure, and rezoning. Both strain measures and material formulation are seen as a homogeneous stress point problem and a flexible material subroutine interface admits both incremental and total strain formulation, dependent on internal energy or an arbitrary set of other internal variables. Vectorization on Class VI computers such as the CRAY-1 is a simple exercise for optimally organized primitive element formulations. Some examples of large scale computation are illustrated, including continuous tone graphic representation.
Method of holding optical elements without deformation during their fabrication
Hed, P.P.
1997-04-29
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.
Higher-Order, Space-Time Adaptive Finite Volume Methods: Algorithms, Analysis and Applications
Minion, Michael
2014-04-29
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.
Validation of Finite-Element Models of Persistent-Current Effects in Nb_{3}Sn Accelerator Magnets
Wang, X.; Ambrosio, G.; Chlachidze, G.; Collings, E. W.; Dietderich, D. R.; DiMarco, J.; Felice, H.; Ghosh, A. K.; Godeke, A.; Gourlay, S. A.; Marchevsky, M.; Prestemon, S. O.; Sabbi, G.; Sumption, M. D.; Velev, G. V.; Xu, X.; Zlobin, A. V.
2015-01-06
Persistent magnetization currents are induced in superconducting filaments during the current ramping in magnets. The resulting perturbation to the design magnetic field leads to field quality degradation, in particular at low field where the effect is stronger relative to the main field. The effects observed in NbTi accelerator magnets were reproduced well with the critical-state model. However, this approach becomes less accurate for the calculation of the persistent-current effects observed in Nb_{3}Sn accelerator magnets. Here a finite-element method based on the measured strand magnetization is validated against three state-of-art Nb3Sn accelerator magnets featuring different subelement diameters, critical currents, magnet designs and measurement temperatures. The temperature dependence of the persistent-current effects is reproduced. Based on the validated model, the impact of conductor design on the persistent current effects is discussed. The performance, limitations and possible improvements of the approach are also discussed.
Method and system for high power reflective optical elements
Demos, Stavros G.; Rubenchik, Alexander M.; Negres, Raluca A.
2013-03-12
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.
METHOD OF FORMING A FUEL ELEMENT FOR A NUCLEAR REACTOR
Layer, E.H. Jr.; Peet, C.S.
1962-01-23
A method is given for preparing a fuel element for a nuclear reactor. The method includes the steps of sandblasting a body of uranium dioxide to roughen the surface thereof, depositing a thin layer of carbon thereon by thermal decomposition of methane, and cladding the uranium dioxide body with zirconium by gas pressure bonding. (AEC)
Whirley, R.G.
1991-05-01
This report is the User Manual for the 1991 version of DYNA3D, and also serves as an interim User Guide. DYNA3D is a nonlinear, explicit, finite element code for analyzing the transient dynamic response of three-dimensional solids and structures. The code is fully vectorized and is available on several computer platforms. DYNA3D includes solid, shell, beam, and truss elements to allow maximum flexibility in modeling physical problems. Many material models are available to represent a wide range of material behavior, including elasticity, plasticity, composites, thermal effects, and rate dependence. In addition, DYNA3D has a sophisticated contact interface capability, including frictional sliding and single surface contact. Rigid materials provide added modeling flexibility. A material model driver with interactive graphics display is incorporated into DYNA3D to permit accurate modeling of complex material response based on experimental data. Along with the DYNA3D Example Problem Manual, this document provides the information necessary to apply DYNA3D to solve a wide range of engineering analysis problems. 73 refs., 49 figs.
Linker, K.L.; Rawlinson, K.S.; Smith, G.
1991-10-01
The Department of Energy`s Solar Thermal Program has as one of its program elements the development and evaluation of conversion device technologies applicable to dish-electric systems. The primary research and development combines a conversion device (heat engine), solar receiver, and generator mounted at the focus of a parabolic dish concentrator. The Stirling-cycle heat engine was identified as the conversion device for dish-electric with the most potential for meeting the program`s goals for efficiency, reliability, and installed cost. To advance the technology toward commercialization, Sandia National Laboratories has acquired a Stirling Thermal Motors, Inc., kinematic Stirling engine, STM4-120, for evaluation. The engine is being bench-tested at Sandia`s Engine Test Facility and will be combined later with a solar receiver for on-sun evaluation. This report presents the engine characteristics, finite element analyses of critical engine components, test system layout, instrumentation, and preliminary performance results from the bench test.
Nuclear fuel elements and method of making same
Schweitzer, Donald G.
1992-01-01
A nuclear fuel element for a high temperature gas nuclear reactor that has an average operating temperature in excess of 2000.degree. C., and a method of making such a fuel element. The fuel element is characterized by having fissionable fuel material localized and stabilized within pores of a carbon or graphite member by melting the fissionable material to cause it to chemically react with the carbon walls of the pores. The fissionable fuel material is further stabilized and localized within the pores of the graphite member by providing one or more coatings of pyrolytic carbon or diamond surrounding the porous graphite member so that each layer defines a successive barrier against migration of the fissionable fuel from the pores, and so that the outermost layer of pyrolytic carbon or diamond forms a barrier between the fissionable material and the moderating gases used in an associated high temperature gas reactor. The method of the invention provides for making such new elements either as generally spherically elements, or as flexible filaments, or as other relatively small-sized fuel elements that are particularly suited for use in high temperature gas reactors.
A stochastic method for computing hadronic matrix elements
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Alexandrou, Constantia; Constantinou, Martha; Dinter, Simon; Drach, Vincent; Jansen, Karl; Hadjiyiannakou, Kyriakos; Renner, Dru B.
2014-01-24
In this study, 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.
(Environmental and geophysical modeling, fracture mechanics, and boundary element methods)
Gray, L.J.
1990-11-09
Technical discussions at the various sites visited centered on application of boundary integral methods for environmental modeling, seismic analysis, and computational fracture mechanics in composite and smart'' materials. The traveler also attended the International Association for Boundary Element Methods Conference at Rome, Italy. While many aspects of boundary element theory and applications were discussed in the papers, the dominant topic was the analysis and application of hypersingular equations. This has been the focus of recent work by the author, and thus the conference was highly relevant to research at ORNL.
A Massively Parallel Sparse Eigensolver for Structural Dynamics Finite Element Analysis
Day, David M.; Reese, G.M.
1999-05-01
Eigenanalysis is a critical component of structural dynamics which is essential for determinating the vibrational response of systems. This effort addresses the development of numerical algorithms associated with scalable eigensolver techniques suitable for use on massively parallel, distributed memory computers that are capable of solving large scale structural dynamics problems. An iterative Lanczos method was determined to be the best choice for the application. Scalability of the eigenproblem depends on scalability of the underlying linear solver. A multi-level solver (FETI) was selected as most promising for this component. Issues relating to heterogeneous materials, mechanisms and multipoint constraints have been examined, and the linear solver algorithm has been developed to incorporate features that result in a scalable, robust algorithm for practical structural dynamics applications. The resulting tools have been demonstrated on large problems representative of a weapon's system.
Morris, J; Johnson, S
2007-12-03
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.
METHOD OF PREPARING A FUEL ELEMENT FOR A NUCLEAR REACTOR
Hauth, J.J.; Anicetti, R.J.
1962-12-01
A method is described for preparing a fuel element for a nuclear reactor. According to the patent uranium dioxide is compacted in a metal tabe by directlng intense sound waves at the tabe prior to tamp packing or vibration compaction of the powder. (AEC)
Finite element analysis of the I-40 bridge over the Rio Grande
Farrar, C.R.; Duffey, T.A.; Goldman, P.A.; Jauregui, D.V.; Vigil, J.S.
1996-01-01
In the 1960s and 1970s numerous bridges were built in the US with a design similar to those on Interstate 40 (I-40) over the Rio Grande in Albuquerque, New Mexico. The Federal Highway Administration (FHWA) and the National Science Foundation (NSF) have provided funds to New Mexico State University (NMSU) through the New Mexico State Highway and Transportation Department (NMSH and TD) and The Alliance For Transportation Research (ATR) for evaluation and testing of the existing fracture-critical bridges over the Rio Grande. Because the I-40 Bridges over the Rio Grande were to be razed during the summer of 1993, the investigators were able to introduce damage into the structure in order to test various damage identification methods and to observe the changes in load paths through the structure caused by the cracking. To support this research effort, NMSU contracted Los Alamos National Laboratory (LANL) to perform experimental modal analyses, and to develop experimentally verified numerical models of the bridge. A previous report (LA-12767-MS) summarizes the results of the experimental modal analyses. This report summarizes the numerical analyses of the bridges and compares the results of these analyses to the experimental results.
Method and apparatus for diagnosing breached fuel elements
Gross, K.C.; Lambert, J.D.B.; Nomura, S.
1987-03-02
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.
Method and apparatus for diagnosing breached fuel elements
Gross, Kenny C.; Lambert, John D. B.; Nomura, Shigeo
1988-01-01
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.
Finite-element analysis of transverse compressive and thermal loads on Nb3Sn wires with voids
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Zhai, Y.; D'Hauthuille, L.; Barth, C.; Senatore, C.
2016-02-29
High-field superconducting magnets play a very important role in many large-scale physics experiments, particularly particle colliders and fusion confinement devices such as Large Hadron Collider (LHC) and International Thermonuclear Experimental Reactor (ITER). The two most common superconductors used in these applications are NbTi and Nb3Sn. Nb3Sn wires are favored because of their significantly higher Jc (critical current density) for higher field applications. The main disadvantage of Nb3Sn is that the superconducting performance of the wire is highly strain sensitive and it is very brittle. This strain sensitivity is strongly influenced by two factors: plasticity and cracked filaments. Cracks are inducedmoreÂ Â» by large stress concentrators that can be traced to the presence of voids in the wire. We develop detailed 2-D and 3-D finite-element models containing wire filaments and different possible distributions of voids in a bronze-route Nb3Sn wire. We apply compressive transverse loads for various cases of void distributions to simulate the stress and strain response of a Nb3Sn wire under the Lorentz force. Furthermore, this paper improves our understanding of the effect voids have on the Nb3Sn wire's mechanical properties, and in so, the connection between the distribution of voids and performance degradation such as the correlation between irreversible strain limit and the void-induced local stress concentrations.Â«Â less
Method for measuring recovery of catalytic elements from fuel cells
Shore, Lawrence; Matlin, Ramail
2011-03-08
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.
Wang, Wei; Shu, Chi-Wang; Yee, H.C.; Sjögreen, Björn
2012-01-01
A new high order finite-difference method utilizing the idea of Harten ENO subcell resolution method is proposed for chemical reactive flows and combustion. In reaction problems, when the reaction time scale is very small, e.g., orders of magnitude smaller than the fluid dynamics time scales, the governing equations will become very stiff. Wrong propagation speed of discontinuity may occur due to the underresolved numerical solution in both space and time. The present proposed method is a modified fractional step method which solves the convection step and reaction step separately. In the convection step, any high order shock-capturing method can be used. In the reaction step, an ODE solver is applied but with the computed flow variables in the shock region modified by the Harten subcell resolution idea. For numerical experiments, a fifth-order finite-difference WENO scheme and its anti-diffusion WENO variant are considered. A wide range of 1D and 2D scalar and Euler system test cases are investigated. Studies indicate that for the considered test cases, the new method maintains high order accuracy in space for smooth flows, and for stiff source terms with discontinuities, it can capture the correct propagation speed of discontinuities in very coarse meshes with reasonable CFL numbers.
Method for detection of antibodies for metallic elements
Barrick, Charles W.; Clarke, Sara M.; Nordin, Carl W.
1993-11-30
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.
Method for detection of antibodies for metallic elements
Barrick, C.W.; Clarke, S.M.; Nordin, C.W.
1993-11-30
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.
Tezaur, I. K.; Perego, M.; Salinger, A. G.; Tuminaro, R. S.; Price, S. F.
2015-04-27
This paper describes a new parallel, scalable and robust finite element based solver for the first-order Stokes momentum balance equations for ice flow. The solver, known as Albany/FELIX, is constructed using the component-based approach to building application codes, in which mature, modular libraries developed as a part of the Trilinos project are combined using abstract interfaces and template-based generic programming, resulting in a final code with access to dozens of algorithmic and advanced analysis capabilities. Following an overview of the relevant partial differential equations and boundary conditions, the numerical methods chosen to discretize the ice flow equations are described, along with their implementation. The results of several verification studies of the model accuracy are presented using (1) new test cases for simplified two-dimensional (2-D) versions of the governing equations derived using the method of manufactured solutions, and (2) canonical ice sheet modeling benchmarks. Model accuracy and convergence with respect to mesh resolution are then studied on problems involving a realistic Greenland ice sheet geometry discretized using hexahedral and tetrahedral meshes. Also explored as a part of this study is the effect of vertical mesh resolution on the solution accuracy and solver performance. The robustness and scalability of our solver on these problems is demonstrated. Lastly, we show that good scalability can be achieved by preconditioning the iterative linear solver using a new algebraic multilevel preconditioner, constructed based on the idea of semi-coarsening.
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
Tezaur, I. K.; Perego, M.; Salinger, A. G.; Tuminaro, R. S.; Price, S. F.
2015-04-27
This paper describes a new parallel, scalable and robust finite element based solver for the first-order Stokes momentum balance equations for ice flow. The solver, known as Albany/FELIX, is constructed using the component-based approach to building application codes, in which mature, modular libraries developed as a part of the Trilinos project are combined using abstract interfaces and template-based generic programming, resulting in a final code with access to dozens of algorithmic and advanced analysis capabilities. Following an overview of the relevant partial differential equations and boundary conditions, the numerical methods chosen to discretize the ice flow equations are described, alongmoreÂ Â» with their implementation. The results of several verification studies of the model accuracy are presented using (1) new test cases for simplified two-dimensional (2-D) versions of the governing equations derived using the method of manufactured solutions, and (2) canonical ice sheet modeling benchmarks. Model accuracy and convergence with respect to mesh resolution are then studied on problems involving a realistic Greenland ice sheet geometry discretized using hexahedral and tetrahedral meshes. Also explored as a part of this study is the effect of vertical mesh resolution on the solution accuracy and solver performance. The robustness and scalability of our solver on these problems is demonstrated. Lastly, we show that good scalability can be achieved by preconditioning the iterative linear solver using a new algebraic multilevel preconditioner, constructed based on the idea of semi-coarsening.Â«Â less
Adaptive finite volume methods for time-dependent P.D.E.S.
Ware, J.; Berzins, M.
1995-12-31
The aim of adaptive methods for time-dependent p.d.e.s is to control the numerical error so that it is less than a user-specified tolerance. This error depends on the spatial discretization method, the spatial mesh, the method of time integration and the timestep. The spatial discretization method and positioning of the spatial mesh points should attempt to ensure that the spatial error is controlled to meet the user`s requirements. It is then desirable to integrate the o.d.e. system in time with sufficient accuracy so that the temporal error does not corrupt the spatial accuracy or the reliability of the spatial error estimates. This paper is concerned with the development of a prototype algorithm of this type, based on a cell-centered triangular finite volume scheme, for two space dimensional convection-dominated problems.
High-order finite-volume methods for hyperbolic conservation laws on mapped multiblock grids
DOE Public Access Gateway for Energy & Science Beta (PAGES Beta)
McCorquodale, P. W.; Colella, P.; Dorr, M. R.; Hittinger, J. A. F.
2015-01-13
We present an approach to solving hyperbolic conservation laws by finite-volume methods on mapped multiblock grids, extending the approach of Colella, Dorr, Hittinger, and Martin (2011) [10] for grids with a single mapping. We consider mapped multiblock domains for mappings that are conforming at inter-block boundaries. By using a smooth continuation of the mapping into ghost cells surrounding a block, we reduce the inter-block communication problem to finding an accurate, robust interpolation into these ghost cells from neighboring blocks. Lastly, we demonstrate fourth-order accuracy for the advection equation for multiblock coordinate systems in two and three dimensions.
Dash, Z.V.; Robinson, B.A.; Zyvoloski, G.A.
1997-07-01
The requirements, design, and verification and validation of the software used in the FEHM application, a finite-element heat- and mass-transfer computer code that can simulate nonisothermal multiphase multicomponent flow in porous media, are described. The test of the DOE Code Comparison Project, Problem Five, Case A, which verifies that FEHM has correctly implemented heat and mass transfer and phase partitioning, is also covered.
Methods for removing transuranic elements from waste solutions
Slater, S.A.; Chamberlain, D.B.; Connor, C.; Sedlet, J.; Srinivasan, B.; Vandegrift, G.F.
1994-11-01
This report outlines a treatment scheme for separating and concentrating the transuranic (TRU) elements present in aqueous waste solutions stored at Argonne National Laboratory (ANL). The treatment method selected is carrier precipitation. Potential carriers will be evaluated in future laboratory work, beginning with ferric hydroxide and magnetite. The process will result in a supernatant with alpha activity low enough that it can be treated in the existing evaporator/concentrator at ANL. The separated TRU waste will be packaged for shipment to the Waste Isolation Pilot Plant.
High-Fidelity RF Gun Simulations with the Parallel 3D Finite Element Particle-In-Cell Code Pic3P
Candel, A; Kabel, A.; Lee, L.; Li, Z.; Limborg, C.; Ng, C.; Schussman, G.; Ko, K.; /SLAC
2009-06-19
SLAC's Advanced Computations Department (ACD) has developed the first parallel Finite Element 3D Particle-In-Cell (PIC) code, Pic3P, for simulations of RF guns and other space-charge dominated beam-cavity interactions. Pic3P solves the complete set of Maxwell-Lorentz equations and thus includes space charge, retardation and wakefield effects from first principles. Pic3P uses higher-order Finite Elementmethods on unstructured conformal meshes. A novel scheme for causal adaptive refinement and dynamic load balancing enable unprecedented simulation accuracy, aiding the design and operation of the next generation of accelerator facilities. Application to the Linac Coherent Light Source (LCLS) RF gun is presented.
Modeling Slag Penetration and Refractory Degradation Using the Finite Element Method
Johnson, Kenneth I.; Williford, Ralph E.; Matyas, Josef; Pilli, Siva Prasad; Sundaram, S. K.; Korolev, Vladimir N.
2008-09-01
Refractory degradation due to slag penetration can significantly reduce the service life of gasifier refractory linings. This paper describes a modeling approach that was developed to predict refractory spalling as a function of operating temperature, coal feedstock and refractory type. The model simulates the coupled thermal, diffusion, and mechanical interactions of coal slag with refractory ceramics. The heat transfer and slag diffusion solutions are directly coupled through a temperature-dependent effective diffusivity for slag penetration. The effective diffusivity is defined from slag penetration tests conducted in our laboratories on specific coal slag and refractory combinations. Chemically-induced swelling of the refractory and the build-up of mechanical stresses are functions of the slag penetration. The model results are compared with analytical spalling models and validated by experimental data in order to develop an efficient refractory degradation model for implementation in a systems level gasifier model. The ultimate goal of our research is to provide a tool that will help optimize gasifier performance by balancing conversion efficiency with refractory life.
Zhou, Xiangyu; Ghione, Giovanni; Bertazzi, Francesco Goano, Michele; Bellotti, Enrico
2014-07-21
We present a multiband envelope-function model for wurtzite nanostructures based on a rigorous numerical procedure to determine operator ordering and band parameters from nonlocal empirical pseudopotential calculations. The proposed approach, implemented within a finite-element scheme, leads to well-posed, numerically stable envelope equations that accurately reproduce full-Brillouin-zone subband dispersions of quantum systems. Although demonstrated here for III-nitride nonlocal empirical pseudopotentials, the model provides a general theoretical framework applicable to ab initio electronic structures of wurtzite semiconductors.
Richgels, M A; Biffle, J H
1980-09-01
ALGEBRA is a program that allows the user to process output data from finite-element analysis codes before they are sent to plotting routines. These data take the form of variable values (stress, strain, and velocity components, etc.) on a tape that is both the output tape from the analyses code and the input tape to ALGEBRA. The ALGEBRA code evaluates functions of these data and writes the function values on an output tape that can be used as input to plotting routines. Convenient input format and error detection capabilities aid the user in providing ALGEBRA with the functions to be evaluated. 1 figure.
Anand, S.C.; Pandit, A.
1983-06-01
In the investigation, a Galerkin finite element model in two dimensions is developed to study the phenomena of mass transfer in porous media. In particular, the problems of the saltwater encroachment in coastal aquifers and the transport of hazardous wastes in groundwater environment are studied for a wide range of aquifer parameters. The coupled governing partial differential equations are nondimensionalized and solved for a two-dimensional, saturated aquifer in the vertical plane for both steady state and transient conditions using an iterative solution procedure. The flow transport is represented either in terms of the stream function or the freshwater hydraulic head.
Method of determining lanthanidies in a transition element host
De Kalb, Edward L.; Fassel, Velmer A.
1976-02-03
A phosphor composition contains a lanthanide activator element within a host matrix having a transition element as a major component. The host matrix is composed of certain rare earth phosphates or vanadates such as YPO.sub.4 with a portion of the rare earth replaced with one or more of the transition elements. On X-ray or other electromagnetic excitation, trace lanthanide impurities or additives within the phosphor are spectrometrically determined from their characteristic luminescence.
System and method for reproducibly mounting an optical element
Eisenbies, Stephen; Haney, Steven
2005-05-31
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.
Method of fabricating reflection-mode EUV diffraction elements
Naulleau, Patrick P.
2002-01-01
Techniques for fabricating a well-controlled, quantized-level, engineered surface that serves as substrates for EUV reflection multilayer overcomes problems associated with the fabrication of reflective EUV diffraction elements. The technique when employed to fabricate an EUV diffraction element that includes the steps of: (a) forming an etch stack comprising alternating layers of first and second materials on a substrate surface where the two material can provide relative etch selectivity; (b) creating a relief profile in the etch stack wherein the relief profile has a defined contour; and (c) depositing a multilayer reflection film over the relief profile wherein the film has an outer contour that substantially matches that of the relief profile. For a typical EUV multilayer, if the features on the substrate are larger than 50 nm, the multilayer will be conformal to the substrate. Thus, the phase imparted to the reflected wavefront will closely match that geometrically set by the surface height profile.