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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Machine learning simulation of finite element analysis in augmented reality

Media, method and system for approximating a finite element analysis texture map for an object. To accomplish this, the object is converted to a computer generated model and finite element analysis is performed for a plurality of different simulated inputs to generate a plurality of simulated mappings. Each simulated mapping is converted into a simulated texture map. A machine learning model is trained on the simulated inputs and simulated texture maps to generate a texture map which approximates a finite element analysis. The machine learning model receives a user input and generates the texture map therefrom. The texture map is then wrapped to the object and displayed.

Scherer, Derek Carl↗
DOE OSTI2016

Finite Volume Method Based Neutronic Solvers for Steady and Transient Analysis of Molten Salt Reactors

In recent years, the coupled neutronics and thermal hydraulics analysis becomes more and more important, because it helps to identify the most relevant safety issues without conservative assumptions. In addition, such coupled simulations are essential for the generation IV advanced nuclear power systems such as molten salt reactor (MSR), because the presence of strong coupling between different phenomena, like the transport of delayed neutron precursors by the fluid flow?. Due to the different discretization methods used in reactor physics analysis and thermal-hydraulics analysis, it is hard to couple different kinds of codes. This work presents an attempt to solve the multigroup neutron diffusion equation with the finite volume method which is quite widely used for the numerical simulation of a variety of applications involving fluid flow and heat and mass transfer. The unstructured finite volume method (FVM) based on open source C++ library OpenFOAM was chosen as the development platform?. The objective of this paper is to concisely present the main features of the developed codes in terms of the theoretical background and numerical approach and as can be seen in the following sections, the calculation results obtained by the NDSFoam which solves the steady-state neutron diffusion equation and NDTFoam which solves the transient-state neutron diffusion equation are discussed and compared to the benchmark problems. In this paper, two new neutronic solvers using finite volume method have been developed to perform the neutron physics analysis. The solvers are based on OpenFOAM, and have been tested with several steady state and transient-state benchmark problems. By verifications of the codes, it is proved that the development of NDSFoam and NDTFoam are correct. In our future work, both of these two solvers will be coupled with the precompiled CFD solvers in OpenFOAM to analyze the molten salt fast reactors. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2018

Finite-Difference Methods for Fractional Differential Equations of Order 1/2

In this work, we study approximations of solutions of fractional differential equations of order 1/2. We present a new method of approximation and obtain the order of convergence. The presentation is given within the abstract framework of a semidiscrete approximation scheme, which includes finite-element methods, finite-difference schemes, and projection methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

High-order finite element method for atomic structure calculations

We introduce featom, an open source code that implements a high-order finite element solver for the radial Schrödinger, Dirac, and Kohn-Sham equations. The formulation accommodates various mesh types, such as uniform or exponential, and the convergence can be systematically controlled by increasing the number and/or polynomial order of the finite element basis functions. The Dirac equation is solved using a squared Hamiltonian approach to eliminate spurious states. Here, to address the slow convergence of the $κ=±1$ states due to divergent derivatives at the origin, we incorporate known asymptotic forms into the solutions. We achieve a high level of accuracy (10 -8 Hartree) for total energies and eigenvalues of heavy atoms such as uranium in both Schrödinger and Dirac Kohn-Sham solutions. We provide detailed convergence studies and computational parameters required to attain commonly required accuracies. Finally, we compare our results with known analytic results as well as the results of other methods. In particular, we calculate benchmark results for atomic numbers (Z) from 1 to 92, verifying current benchmarks. We demonstrate significant speedup compared to the state-of-the-art shooting solver dftatom. An efficient, modular Fortran 2008 implementation, is provided under an open source, permissive license, including examples and tests, wherein particular emphasis is placed on the independence (no global variables), reusability, and generality of the individual routines.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Assessing a mini-application as a performance proxy for a finite element method engineering application

The performance of a large-scale, production-quality science and engineering application (‘app’) is often dominated by a small subset of the code. Even within that subset, computational and data access patterns are often repeated, so that an even smaller portion can represent the performance-impacting features. If application developers, parallel computing experts, and computer architects can together identify this representative subset and then develop a small mini-application (‘miniapp’) that can capture these primary performance characteristics, then this miniapp can be used to both improve the performance of the app as well as provide a tool for co-design for the high-performance computing community. However, a critical question is whether a miniapp can effectively capture key performance behavior of an app. This study provides a comparison of an implicit finite element semiconductor device modeling app on unstructured meshes with an implicit finite element miniapp on unstructured meshes. The goal is to assess whether the miniapp is predictive of the performance of the app. Finally, single compute node performance will be compared, as well as scaling up to 16,000 cores. Results indicate that the miniapp can be reasonably predictive of the performance characteristics of the app for a single iteration of the solver on a single compute node.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Automatically adaptive stabilized finite elements and continuation analysis for compaction banding in geomaterials

Under compressive creep, viscoplastic solids experiencing internal mass transfer processes can accommodate singular cnoidal wave solutions as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the material's capacity to diffusive internal perturbations, leading to localized failure features (e.g., cracks and compaction bands). These cnoidal waves, generally found in fluids, have strong nonlinearities that produce periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, practical simulation tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This automatic refinement method provides an error estimate that drives mesh adaptivity, a crucial feature for the problem at hand. We compare the performance of this adaptive strategy against analytical and standard finite element solutions. We then investigate the sensitivity of the diffusivity ratio, the parameter controlling the process, and identify multiple possible solutions with several stress peaks. We also show the evolution of the spacing between peaks for all solutions as a function of that parameter.

42 ENGINEERING↗
DOE OSTI2017

From cells to laminate: probing and modeling residual stress evolution in thin silicon photovoltaic modules using synchrotron X-ray micro-diffraction experiments and finite element simulations

Fracture of silicon crystalline solar cells has recently been observed in increasing percentages especially in solar photovoltaic (PV) modules involving thinner silicon solar cells (<200 μm). Many failures due to fracture have been reported from the field because of environmental loading (snow, wind, etc.) as well as mishandling of the solar PV modules (during installation, maintenance, etc.). However, a significantly higher number of failures have also been reported during module encapsulation (lamination) indicating high residual stress in the modules and thus more prone to cell cracking. Here in this paper we report through the use of synchrotron X-ray submicron diffraction coupled with physics-based finite element modeling, the complete residual stress evolution in mono-crystalline silicon solar cells during PV module integration process. For the first time, we unravel the reason for the high stress and cracking of silicon cells near soldered inter-connects. Our experiments revealed a significant increase of residual stress in the silicon cell near the solder joint after lamination. Moreover, our finite element simulations show that this increase of stress during lamination is a result of highly localized bending of the cell near the soldered inter-connects. Further, the synchrotron X-ray submicron diffraction has proven to be a very effective way to quantitatively probe mechanical stress in encapsulated silicon solar cells. Thus, this technique has ultimately enabled these findings leading to the enlightening of the role of soldering and encapsulation processes on the cell residual stress. This model can be further used to suggest methodologies that could lead to lower stress in encapsulated silicon solar cells, which are the subjects of our continued investigations.

14 SOLAR ENERGY↗
DOE OSTI2019

deal.II Implementation of a Weak Galerkin Finite Element Solver for Darcy Flow

This paper presents a weak Galerkin (WG) finite element solver for Darcy flow and its implementation on the deal.II platform. The solver works for quadrilateral and hexahedral meshes in a unified way. It approximates pressure by Q-type degree $k(≥ 0)$ polynomials separately defined in element interiors and on edges/faces. In this work, numerical velocity is obtained in the unmapped Raviart-Thomas space $RT_{[k]}$ via postprocessing based on the novel concepts of discrete weak gradients. The solver is locally mass-conservative and produces continuous normal fluxes. The implementation in deal.II allows polynomial degrees up to 5. Numerical experiments show that our new WG solver performs better than the classical mixed finite element methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Scalability Tests of a Finite Element Code on Hundreds of Thousands Cores and Heterogeneous Architecture

A multi-scale finite element method code, msFEM, is tested on Jaguar and Nebulae, two petaflops computers that were listed as #1 and #2 on the Top500 list of June 2010 at the time of the tests. The flat MPI version of msFEM is scaled from 20K up to 200K CPU cores on Jaguar, delivering 70% parallel efficiency at the 200K cores with a finite element model of eight millions of degrees of freedom. GPU versions, in both double precision and mixed precision coded through MPI+OpenMP+CUDA hybrid programming, 900 GPU nodes on Jaguar and 1500 GPU nodes on Nebulae, achieving remarkable 90 + % parallel efficiency on the systems. The mixed-precision GPU version delivers further 1.5 times of speedup over the fully double precision version with no significant implementational cost. The large-scale tests support that the msFEM runs efficiently on petaflops computers and is highly potential for domain applications at extreme-scale.

finite element method↗
DOE OSTI2017

Finite element formulation with embedded weak discontinuities for strain localization under dynamic conditions

Here, we present an explicit finite element formulation designed for the treatment of strain localization under highly dynamic conditions. We also used a material stability analysis to detect the onset of localization behavior. Finite elements with embedded weak discontinuities are employed with the aim of representing subsequent localized deformation accurately. The formulation and its algorithmic implementation are described in detail. Numerical results are presented to illustrate the usefulness of this computational framework in the treatment of strain localization under highly dynamic conditions, and to examine its performance characteristics in the context of two-dimensional plane-strain problems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2019

Fatigue Crack Propagation of Nickel-Based Superalloy: Experiments and Simulations with Extended Finite Element Method

Numerical simulation based on extended finite element method was employed to investigate the fatigue crack propagation of nickel-based superalloy at room temperature. Experimental tests on compact tension specimens have performed to obtain fatigue crack propagation parameters in Paris region. The extended finite element method has presented a new approach to solve the stress intensity factors and can effectively predict crack propagation without re-meshing at crack tip. The simulation results are in good accordance with experimental data in real 3D cases.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Mixed finite-element method for multi-term time-fractional diffusion and diffusion-wave equations

In this paper, numerical theory based on the mixed finite-element method and finite difference analog of the Caputo fractional derivative for multi-term time-fractional diffusion equations and diffusion-wave equations is analyzed. The unconditional stability and convergence results are proved for the two resulting fully discrete schemes. Finally, the obtained results are supported by numerical experiments carried out for some test problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

A front-fixing finite element method for pricing American options under regime-switching jump-diffusion models

A front-fixing finite element method is applied to solve the partial integro-differential equations (PIDEs) arising in pricing American options under Markov-modulated jump-diffusion models with the free boundaries feature. For this purpose, we used a front-fixing method to transform the pricing problem into a nonlinear parabolic integro-differential equation on a fixed domain. Then the variational form of the resulting problem is solved by a finite element method. Under some appropriate assumptions, we establish the stability of the method and illustrate some numerical results to examine the rate of convergence of the proposed method for the pricing problem and compare its accuracy to some recent works on pricing American options under regime-switching jump-diffusion models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Crank–Nicolson WSGI difference scheme with finite element method for multi-dimensional time-fractional wave problem

In this article, a second-order Crank–Nicolson weighted and shifted Grünwald integral (WSGI) time-discrete scheme combined with finite element method is studied for finding the numerical solution of the multi-dimensional time-fractional wave equation. The time-fractional wave equation with Caputo-fractional derivative is transformed into the time-fractional integral equation by integral technique, then the second-order Crank–Nicolson finite element scheme with a time second-order WSGI discrete approximation for fractional integral is formulated. The analysis of stability is made, then a priori error estimate is given by making use of the conversion technique between the WSGI formula and fractional integral. At the end of the article, some numerical examples covering two two-dimensional cases with rectangular element and triangular element and a three-dimensional case with tetrahedral element are shown to test and verify our theoretical results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

HOSS: an implementation of the combined finite-discrete element method

Nearly thirty years since its inception, the combined finite-discrete element method (FDEM) has made remarkable strides in becoming a mainstream analysis tool within the field of Computational Mechanics. FDEM was developed to effectively “bridge the gap” between two disparate Computational Mechanics approaches known as the finite and discrete element meth-ods. At Los Alamos National Laboratory (LANL) researchers developed the Hybrid Optimization Software Suite (HOSS) as a hybrid multi-physics platform, based on FDEM, for the simulation of solid material behavior complemented with the latest technological enhancements for full fluid–solid interaction. Furthermore, in HOSS, several newly developed FDEM algorithms have been implemented that yield more accurate material deformation formulations, inter-particle interaction solvers, and fracture and fragmentation solutions. Additionally, an explicit computational fluid dynamics solver and a novel fluid–solid interaction algorithms have been fully integrated (as opposed to coupled) into the HOSS’ solid mechanical solver, allowing for the study of an even wider range of problems. Advancements such as this are leading HOSS to become a tool of choice for multi-physics problems. Finally, HOSS has been successfully applied by a myriad of researchers for analysis in rock mechanics, oil and gas industries, engineering application (structural, mechanical and biomedical engineering), mining, blast loading, high velocity impact, as well as seismic and acoustic analysis. This paper intends to summarize the latest development and application efforts for HOSS.

58 GEOSCIENCES↗
DOE OSTI2024

Efficient sensitivity analysis of the thermal profile in powder bed fusion of metals using hypercomplex automatic differentiation finite element method

Rapid cyclic temperature fluctuation occurring in powder bed fusion of metals using a laser beam (PBF-LB/M) influences the formation of flaws in printed parts. Consequently, there is a pressing need to enhance the quality of printed parts by developing innovative methodologies that can predict thermal histories and help uncover the intricate relationships between process parameters and thermal profiles. Sensitivity Analysis (SA) emerges as an essential tool for this, offering the potential for process optimization and enhanced quality control. Nonetheless, conventional SA methodologies often incur in excessive computational costs and potential numerical approximation errors. Here, to address this technical challenge, we present a novel method for SA that integrates the HYPercomplex-based Automatic Differentiation (HYPAD) technique with transient thermal simulations conducted via the finite element method (FEM). Leveraging this methodology, we efficiently and accurately perform SA for PBF-LB/M processes in a post-processing step. Compared to traditional methods like Finite Differences (FD), HYPAD-FEM required 96 % less computational time for obtaining sensitivities for 22 process parameters, under a comparative study conducted within the context of the 2018–02 AM benchmark of the National Institute of Standards and Technology. In summary, HYPAD-FEM offers superior efficiency and accuracy in SA over conventional methods, delivering the best sensitivity of a model without the need for step-size selection and problem or parameter-based implementations.

36 MATERIALS SCIENCE↗