Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Fracture Capabilities in Grizzly with the extended Finite Element Method (X-FEM)

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

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Finite element modeling of biomass hopper flow

Hoppers are widely used biomass handling devices that channel bulk biomass from storage to subsequent handling equipment. Jenike’s longstanding approach, based on the Mohr-Coulomb model, has been successfully used to design hoppers handling cohesionless granular materials such as grains and other agricultural produces. However, designing a hopper to ensure reliable biomass flow is found to be challenging due to cohesion, irregular particle shape, and bulk material elastoplasticity. This study aims to address the biomass handling engineering challenge with alternative constitutive material models concerning the flow behavior of bulk solids. Finite element modeling is an approach that allows for implementation of different material models, whose underlying constitutive theories assist in investigating the origin and manifestation of bulk mechanical behavior of granular materials. This study focuses on the incipient gravity hopper flow of two types of biomass feedstocks, i.e., ground corn stover and Douglas fir wood. Three widely used constitutive material models, i.e., Mohr-Coulomb model, modified Cam-Clay model, and Drucker-Prager/Cap model, are implemented. Using the flow pattern represented by the volume of biomass exhibiting more than 7% of axial strain (Kamath and Puri, 1999), the finite element model predicts that the bulk corn stover particulate material forms an arch, which represents a hampered transition from the static state to the dynamic flow-state out of the hopper, whereas bulk Douglas fir wood particulate material develops a reliable mass flow pattern. A laboratory scale hopper was used to experimentally determine the biomass flow conditions, which were subsequently compared with the predicted onset of flow by a finite element model (FEM). The developed FEM was found to correctly predict the initiation of mass flow for the milled Douglas fir wood, whereas corn stover was predicted to establish a strong core flow suggesting an unreliable handling characteristic. This observation aligns with the reported poor handling of milled corn stover.

09 BIOMASS FUELS↗
DOE OSTI2024

A Comprehensive Review on Finite Element Analysis of Laser Shock Peening

Laser shock peening (LSP) is a formidable cold working surface treatment that provides high-energy precision to enhance the mechanical properties of materials. This paper delves into the intricacies of the LSP process, offering insights into its methodology and the simulation thereof through the finite element method. This review critically examines various points, such as laser energy, overlapping of shots, effect of LSP on residual stress, effect of LSP on grain refinement, and algorithms for simulation extrapolated from finite element analyses conducted by researchers, shedding light on the nuanced considerations integral to this technique. As the significance of LSP continues to grow, the collective findings underscore its potential as a transformative technology for fortifying materials against mechanical stress and improving their overall performance and longevity. The discourse encapsulates the evolving landscape of the LSP, emphasizing the pivotal role played by finite element analysis in advancing our understanding and application of this innovative surface treatment.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Variational Asymptotic Homogenization of Finitely Deformed Viscoelastic-Viscoplastic Composites

The objective of this paper is to develop a constitutive model for finitely deformed viscoelastic-viscoplastic materials and a micromechanics approach to homogenizing composites consisting of such materials. The development of the constitutive model involves establishing a thermodynamic framework based on finite strain theory, developing a viscoelasticity and a viscoplasticity model based on the thermodynamic framework, developing a radial return algorithm based on a classic framework, and deriving a closed-from incremental constitutive relation in the spatial configuration. The development of the micromechanics approach involves pulling-back the above constitutive relation to the material configuration, formulating a variational statement with the resulting constitutive relation, discretizing variational statement in a finite-dimensional space, and solving the discretized variational statement using an Euler–Newton predictor–corrector method. The constitutive model is calibrated via monotonic uniaxial compression tests on a polymer, and the calibrated model is validated by comparing its predictions with the cyclic test data. It is shown capable of characterizing viscoelasticity, viscoplasticity, and complex loading paths. The micromechanics approach’s capabilities are demonstrated through homogenizing a unidirectional fiber-reinforced composite, subjected to uniaxial, biaxial, and shear loading, at different strain rates. It is demonstrated to be capable of handling rate dependence and complex loading paths. The present framework can be further improved by implementing more sophisticated viscoelasticity and viscoplasticity models in future work.

Finite element analysis↗
DOE OSTI2019

Electromagnetics Simulations with Vector-valued Finite Elements in MOOSE

Vector-valued finite elements have become an essential part of modern electromagnetic (EM) numerical simulation. This talk will describe the development and initial validation of a MOOSE-based EM solver application, ELK (Electromagnetics Library for Kinetics and fluids), that utilizes the recently added vector-valued finite element system. The primary motivation for the development of ELK was to study EM wave propagation in radio-frequency (rf) driven systems, such as those encountered in the study of antenna structures, waveguides, and plasma physics. In simulating these scenarios with traditional scalar-valued finite elements, boundary conditions on tangential field components can sometimes be difficult to define, and field discontinuities encountered with certain geometric or structural features (e.g., abrupt corners and inhomogeneous materials) can often lead to non-physical results. Vector-valued elements can largely mitigate and even eliminate these difficulties. ELK verification and validation efforts thus far have focused on standard wave propagation benchmarks for antennas and waveguides, and the impact of vector-valued elements on the implementation and results of these simulations will be described. Progress using ELK for plasma physics simulation, where inhomogeneous simulation domains are commonplace, will also be discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Pellet cladding mechanical interaction modeling using the extended finite element method

Fracturing of light water reactor (LWR) UO{sub 2} pellet has multiple important effects on fuel performance. One particularly important concern is that cracks in the fuel cause elevated stresses in the cladding when pellet cladding mechanical interaction (PCMI) occurs. The extended finite element method (XFEM) is a powerful method to represent arbitrary propagating discrete cracks in finite element models. The use of XFEM has been previously demonstrated for modeling propagating discrete cracks in the BISON fuel performance code. This paper demonstrates an initial application of XFEM to model stress concentrations induced by fuel fractures at the fuel/cladding interface during PCMI. This is demonstrated on a study of a pre-defined stationary crack in a 2D cross-section model of a LWR fuel rod. The results from a model with a discrete crack defined with XFEM compare favorably with the results from a model with the same crack geometry defined in the finite element mesh. This study focuses on benchmarking the use of XFEM for PCMI modeling with a stationary crack, but this technique will be readily extended in the future to consider PCMI in conjunction with arbitrary, physics-driven crack propagation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

An implementation of a high-order generalized finite difference method for solving the time-harmonic cold plasma wave equation in toroidal geometry

A high-order physics-informed meshless finite difference numerical technique is introduced for solving the time-harmonic cold plasma wave equation in toroidal geometries, presenting a novel application of the generalized finite difference (GFD) method to plasma wave simulations. The algorithm employs an irregular distribution of computational points, with local point density informed by the shortest wavelength derived from the cold plasma dispersion relation. Numerical stability and robustness are addressed using regularization techniques. The algorithm, implemented for two spatial dimensions, solves for the wave electric field and is demonstrated to achieve convergence rates of $\mathcal{O}$($\mathcal{h}$ $\mathcal{P}$ )⁠. Verification tests reproduce plane wave solutions, and example simulations of ion cyclotron resonance heating and electron cyclotron resonance heating demonstrate its capability, approaching realistic tokamak plasma scenarios. This work contributes to laying a foundation for the GFD method to be used in more sophisticated, optimized, and physically realistic full-wave simulations in time-harmonic plasma wave research.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Modeling curved interfaces without element‐partitioning in the extended finite element method

Summary In this paper, we model holes and material interfaces (weak discontinuities) in two‐dimensional linear elastic continua using the extended finite element method on higher‐order (spectral) finite element meshes. Arbitrary parametric curves such as rational Bézier curves and cubic Hermite curves are adopted in conjunction with the level set method to represent curved interfaces. Efficient computation of weak form integrals with polynomial integrands is realized via the homogeneous numerical integration scheme—a method that uses Euler's homogeneous function theorem and Stokes' theorem to reduce integration to the boundary of the domain. Numerical integration on cut elements requires the evaluation of a one‐dimensional integral over a parametric curve, and hence, the need to partition curved elements is eliminated. To improve stiffness matrix conditioning, ghost penalty stabilization and the Jacobi preconditioner are used. For material interface problems, we develop an enrichment function that captures weak discontinuities on spectral meshes. Taken together, we show through numerical experiments that these advances deliver optimal algebraic rates of convergence with h ‐refinement ( p =1,2,…,5) and exponential rates of convergence with p ‐refinement ( p =1,2,…,7) for elastostatic problems with holes and material inclusions on Cartesian p th‐order spectral finite element meshes.

Chin, Eric B.↗
DOE OSTI2014

Iterative structure of finite loop integrals

In this paper we develop further and refine the method of differential equations for computing Feynman integrals. In particular, we show that an additional iterative structure emerges for finite loop integrals. As a concrete non-trivial example we study planar master integrals of light-by-light scattering to three loops, and derive analytic results for all values of the Mandelstam variables s and t and the mass m. We start with a recent proposal for defining a basis of loop integrals having uniform transcendental weight properties and use this approach to compute all planar two-loop master integrals in dimensional regularization. We then show how this approach can be further simplified when computing finite loop integrals. This allows us to discuss precisely the subset of integrals that are relevant to the problem. We find that this leads to a block triangular structure of the differential equations, where the blocks correspond to integrals of different weight. We explain how this block triangular form is found in an algorithmic way. Another advantage of working in four dimensions is that integrals of different loop orders are interconnected and can be seamlessly discussed within the same formalism. We use this method to compute all finite master integrals needed up to three loops. Finally, we remark that all integrals have simple Mandelstam representations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Thermodynamics and CP-odd transport in holographic QCD with finite magnetic field

We consider a bottom-up holographic model of QCD at finite temperature T and magnetic field B, and study dependence of thermodynamics and CP-odd transport on these variables. As the magnetic field couples to the flavor sector only, one should take the Veneziano limit where the number of flavors and colors are large while their ratio is kept fixed. We investigate the corresponding holographic background in the approximation where the ratio of flavors to colors is finite but small. We demonstrate that B-dependence of the entropy of QCD is in qualitative agreement with the recent lattice studies. Finally we study the CP-odd transport properties of this system. In particular, we determine the Chern-Simons decay rate at finite B and T, that is an important ingredient in the Chiral Magnetic Effect.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Finite Convergence of a Subgradient Projections Method with Expanding Controls

We study finite convergence of the modified cyclic subgradient projections (MCSP) algorithm for the convex feasibility problem (CFP) in the Euclidean space. Expanding control sequences allow the indices of the sets of the CFP to re-appear and be used again by the algorithm within windows of iteration indices whose lengths are not constant but may increase without bound. Motivated by another development in finitely convergent sequential algorithms that has a significant real-world application in the field of radiation therapy treatment planning, we show that the MCSP algorithm retains its finite convergence when used with an expanding control that is repetitive and fulfills an additional condition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Mean-Variance Problems for Finite Horizon Semi-Markov Decision Processes

This paper deals with a mean-variance problem for finite horizon semi-Markov decision processes. The state and action spaces are Borel spaces, while the reward function may be unbounded. The goal is to seek an optimal policy with minimal finite horizon reward variance over the set of policies with a given mean. Using the theory of N-step contraction, we give a characterization of policies with a given mean and convert the second order moment of the finite horizon reward to a mean of an infinite horizon reward/cost generated by a discrete-time Markov decision processes (MDP) with a two dimension state space and a new one-step reward/cost under suitable conditions. We then establish the optimality equation and the existence of mean-variance optimal policies by employing the existing results of discrete-time MDPs. We also provide a value iteration and a policy improvement algorithms for computing the value function and mean-variance optimal policies, respectively. In addition, a linear program and the dual program are developed for solving the mean-variance problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Few-Body Bound States and Resonances in Finite Volume

Since the pioneering work of Lüscher in the 1980s it is well known that considering quantum systems in finite volume, specifically, finite periodic boxes, can be used as a powerful computational tool to extract physical observables. While this formalism has been worked out in great detail in the two-body sector, much effort is currently being invested into deriving analogous relations for systems with more constituents. This work is relevant not only for nuclear physics, where lattice methods are now able to calculate few- and many-nucleon states, but also for other fields such as simulations of cold atoms. This article discusses recent progress regarding the extraction of few-body bound-state and resonance properties from finite-volume calculations of systems with an arbitrary number of constituents.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A Discrete Nonlinear Schrödinger-Type Hierarchy, Its Finite-Dimensional Reduction Analysis, and the Numerical Integration Scheme

We investigate the procedures of discretization of the integrable nonlinear Schrödinger dynamical system, well known as the Ablowitz–Ladik equation, the corresponding symplectic structures, and the finite-dimensional invariant reductions. We develop an efficient scheme of invariant reduction of the corresponding infinite system of ordinary differential equations to an equivalent finite system of ordinary differential equations with respect to the evolution parameter. We construct a finite set of recurrence algebraic regular relations that allows one to generate solutions of the discrete nonlinear Schrödinger dynamical system and discuss the related functional spaces of solutions. Finally, we analyze the Fourier-transform approach to the study of the set of solutions of the discrete nonlinear Schrödinger dynamical system and its functional-analytic aspects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Transformation of double-half inverse Gaussian hollow beams into superposition of finite Airy beams using an optical Airy transform

In the present work, we investigate theoretically the transformation of a double-half inverse Gaussian hollow (DHIGH) beam into a superposition of finite Airy beams by using an optical Airy transform system. The analytical expression of the generated finite Airy beam is derived by means of the generalized Huygens–Fresnel integral diffraction. Numerical calculations show that the optical Airy transform system converts the central dark spot of the hollow incident beam into peak intensity. It is demonstrated that the characteristics of the transverse intensity distribution of the output finite Airy beam are determined by the phase pattern’s constants of the Airy transform system and the radius size of the incident DHIGH beam. Moreover, it is shown that the profile of the central peak intensity of the generated beam can be controlled by means of the radius size of the incident hollow DHIGH beam.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Compact finite difference relaxation method for chaotic and hyperchaotic initial value systems

In this paper, we present a new application of higher order compact finite differences to solve nonlinear initial value problems exhibiting chaotic behaviour. The method involves dividing the domain of the problem into multiple sub-domains, with each sub-domain integrated using higher order compact finite difference schemes. The nonlinearity is dealt using a Gauss–Seidel-like relaxation. The method is, therefore, referred to as the multi-domain compact finite difference relaxation method (MD-CFDRM). In this new application, the MD-CFDRM is used to solve famous chaotic systems and hyperchaotic systems. The main advantage of the new approach is that it offers better accuracy on coarser grids which significantly improves the computational speed of the method. The results are compared with spectral-based multi-domain method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Towards space from Hilbert space: finding lattice structure in finite-dimensional quantum systems

We report that field theories place one or more degrees of freedom at every point in space. Hilbert spaces describing quantum field theories, or their finite-dimensional discretizations on lattices, therefore have large amounts of structure: they are isomorphic to the tensor product of a smaller Hilbert space for each lattice site or point in space. Local field theories respecting this structure have interactions which preferentially couple nearby points. The emergence of classicality through decoherence relies on this framework of tensor-product decomposition and local interactions. We explore the emergence of such lattice structure from Hilbert-space considerations alone. We point out that the vast majority of finite-dimensional Hilbert spaces cannot be isomorphic to the tensor product of Hilbert-space subfactors that describes a lattice theory. A generic Hilbert space can only be split into a direct sum corresponding to a basis of state vectors spanning the Hilbert space, we consider setups in which the direct sum is naturally decomposed into two pieces. We define a notion of direct-sum locality which characterizes states and decompositions compatible with Hamiltonian time evolution. Finally, we illustrate these notions for a toy model that is the finite-dimensional discretization of the quantum-mechanical double-well potential. We discuss their relevance in cosmology and field theory, especially for theories which describe a landscape of vacua with different spacetime geometries.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗