Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element”

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 973 records · Page 54

Electromagnetic Extended Finite Elements for High-Fidelity Multimaterial Problems LDRD Final Report

Surface effects are critical to the accurate simulation of electromagnetics (EM) as current tends to concentrate near material surfaces. Sandia EM applications, which include exploding 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 impractical 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 empirical 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 satisfy 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.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

A blended continuous–discontinuous finite element method for solving the multi-fluid plasma model

The multi-fluid plasma model represents electrons, multiple ion species, and multiple neutral species as separate fluids that interact through short-range collisions and long-range electromagnetic fields. The model spans a large range of temporal and spatial scales, which renders the model stiff and presents numerical challenges. To address the large range of timescales, a blended continuous and discontinuous Galerkin method is proposed, where the massive ion and neutral species are modeled using an explicit discontinuous Galerkin method while the electrons and electromagnetic fields are modeled using an implicit continuous Galerkin method. This approach is able to capture large-gradient ion and neutral physics like shock formation, while resolving high-frequency electron dynamics in a computationally efficient manner. The details of the Blended Finite Element Method (BFEM) are presented. The numerical method is benchmarked for accuracy and tested using two-fluid one-dimensional soliton problem and electromagnetic shock problem. The results are compared to conventional finite volume and finite element methods, and demonstrate that the BFEM is particularly effective in resolving physics in stiff problems involving realistic physical parameters, including realistic electron mass and speed of light. The benefit is illustrated by computing a three-fluid plasma application that demonstrates species separation in multi-component plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Discontinuous finite element space-angle treatment of the first order linear Boltzmann transport equation with magnetic fields: Application to MRI-guided radiotherapy

Purpose: The advent of magnetic resonance imaging (MRI) guided radiotherapy systems demands the incorporation of the magnetic field into dose calculation algorithms of treatment planning systems. This is due to the fact that the Lorentz force of the magnetic field perturbs the path of the relativistic electrons, hence altering the dose deposited by them. Building on the previous work, the authors have developed a discontinuous finite element space-angle treatment of the linear Boltzmann transport equation to accurately account for the effects of magnetic fields on radiotherapy doses. Methods: The authors present a detailed description of their new formalism and compare its accuracy to GEANT4 Monte Carlo calculations for magnetic fields parallel and perpendicular to the radiation beam at field strengths of 0.5 and 3 T for an inhomogeneous 3D slab geometry phantom comprising water, bone, and air or lung. The accuracy of the authors’ new formalism was determined using a gamma analysis with a 2%/2 mm criterion. Results: Greater than 98.9% of all points analyzed passed the 2%/2 mm gamma criterion for the field strengths and orientations tested. The authors have benchmarked their new formalism against Monte Carlo in a challenging radiation transport problem with a high density material (bone) directly adjacent to a very low density material (dry air at STP) where the effects of the magnetic field dominate collisions. Conclusions: A discontinuous finite element space-angle approach has been proven to be an accurate method for solving the linear Boltzmann transport equation with magnetic fields for cases relevant to MRI guided radiotherapy. The authors have validated the accuracy of this novel technique against GEANT4, even in cases of strong magnetic field strengths and low density air.

61 RADIATION PROTECTION AND DOSIMETRY↗

Simulation of thin slot spirals and dual circular patch antennas using the finite element method with mixed elements

This semi-annual report describes progress up to mid-January 1995. The report contains five sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. Of significance is the development of decomposition schemes which separate the different regions of the antenna volume. Substantial effort was devoted to improving the feed model in the context of the finite element method (FEM). Finally, an innovative scheme for truncating finite element meshes is presented.

Gong, Jian↗

NASTRAN thermal analyzer: A general purpose finite element heat transfer computer program

The program not only can render temperature distributions in solids subjected to various thermal boundary conditions, including effects of diffuse-gray thermal radiation, but is fully compatible in capacity and in the finite-element model representation with that of its structural counterpart in the NASTRAN system. The development history of the finite-element approach for determining temperatures is summarized. The scope of analysis capability, program structure, features, and limitations are given with the objective of providing NASTRAN users with an overall veiw of the NASTRAN thermal analyzer.

Lee, H.↗

Application of a finite element method for computing grazing incidence wave structure in an impedance tube - Comparison with experiment

The acoustic performance of a liner specimen, in a grazing incidence impedance tube, is analyzed using a finite element method. The liner specimen was designed to be a locally reacting, two-degree-of-freedom type with the resistance and reactance provided by perforated facesheets and compartmented cavities. Measured and calculated wave structures are compared for both normal and grazing incidence from 0.3 to 1.2 kHz. A finite element algorithm was incorporated into an optimization loop in order to predict liner grazing incidence impedance from measured SWR and null position data. Results suggest that extended reaction effects may have been responsible for differences between normal and grazing incidence impedance estimates.

Lester, H. C.↗

Transient conduction-radiation analysis of an absolute active cavity radiometer using finite elements

A NASA-developed finite element-based model of a generic active cavity radiometer (ACR) has been developed in order to study the dependence on operating temperature of the closed-loop and open-loop transient response of the instrument. Transient conduction within the sensing element is explored, and the transient temperature distribution resulting from the application of a time-varying radiative boundary condition is calculated. The results verify the prediction that operation of an ACR at cryogenic temperatures results in large gains in frequency response.

Mahan, J. R.↗

Finite-element analysis and multibody dynamics issues in rotorcraft dynamic analysis

There is general agreement that the development of effective rotorcraft analysis software will require the use of modern computational mechanics methodologies, especially finite element analysis and multibody dynamics. This paper examines the analysis of rotorcraft dynamics from the perspective of these methodologies. First, a general discussion of rotorcraft analysis and modeling is presented. Then, a hierarchy of rotorcraft analyses is presented, ranging from simple to complex kinematics, where it is shown that in comprehensive rotorcraft software, finite element analysis must be augmented by multibody dynamics in order to properly analyze large motions of rotorcraft components. Finally, a review of multibody dynamics is presented to further familiarize the rotorcraft community with this technology.

Ruzicka, Gene C.↗

Scattering and radiation analysis of three-dimensional cavity arrays via a hybrid finite element method

A hybrid numerical technique is presented for a characterization of the scattering and radiation properties of three-dimensional cavity arrays recessed in a ground plane. The technique combines the finite element and boundary integral methods and invokes Floquet's representation to formulate a system of equations for the fields at the apertures and those inside the cavities. The system is solved via the conjugate gradient method in conjunction with the Fast Fourier Transform (FFT) thus achieving an O(N) storage requirement. By virtue of the finite element method, the proposed technique is applicable to periodic arrays comprised of cavities having arbitrary shape and filled with inhomogeneous dielectrics. Several numerical results are presented, along with new measured data, which demonstrate the validity, efficiency, and capability of the technique.

Jin, Jian-Ming↗

Finite element analyses of single particle crushing tests incorporating computed tomography imaging and damage mechanics

Fracture of individual grains within a granular assembly affect the strength and deformation behavior of granular materials under large stress states. Interparticle forces, leading to contact stresses and ultimately fracture initiation, are influenced by particle morphology. A numerical method that incorporates both particle fracture and morphology can provide a more accurate model of a granular material’s macro-scale response. This research investigates an approach to modeling fracture initiation and propagation within individual grains, utilizing high resolution X-ray computed tomography (CT) imaging to consider grain morphology and an explicit finite element code optimized for high performance computing which incorporates damage mechanics. Using precisely measured force-displacement curves from single particle crushing tests of manufactured quartz spheres, the method is validated, and the effects of fracture properties on grain fragmentation are determined. This approach is then used to simulate single particle crushing of Ottawa sand, with grain morphology incorporated through discretized CT images. Here the effect of grain orientation on the force and applied displacement at catastrophic splitting is also explored. The proposed discrete particle based finite element framework can be applied to modeling granular assemblies, considering the effect of individual particle shape and fracture on the assembly’s deformation response through well calibrated numerical simulations.

42 ENGINEERING↗

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↗

CoCrMo cellular structures made by Electron Beam Melting studied by local tomography and finite element modelling

The work focuses on the structural and mechanical characterization of Co-Cr-Mo cellular samples with cubic pore structure made by Electron Beam Melting (EBM). X-ray tomography was used to characterize the architecture of the sample. High resolution images were also obtained thanks to local tomography in which the specimen is placed close to the X-ray source. These images enabled to observe some defects due to the fabrication process: small pores in the solid phase, partially melted particles attached to the surface. Then, in situ compression tests were performed in the tomograph. The images of the deformed sample show a progressive buckling of the vertical struts leading to final fracture. The deformation initiated where the defects were present in the strut i.e. in regions with reduced local thickness. The finite element modelling confirmed the high stress concentrations of these weak points leading to the fracture of the sample. - Highlights: • CoCrMo samples fabricated by Electron Beam Melting (EBM) process are considered. • X-ray Computed Tomography is used to observe the structure of the sample. • The mechanical properties are tested thanks to an in situ test in the tomograph. • A finite element model is developed to model the mechanical behaviour.

36 MATERIALS SCIENCE↗

An arbitrary Lagrangian–Eulerian finite element formulation for a poroelasticity problem stemming from mixture theory

In this paper, a finite element formulation is developed for a poroelastic medium consisting of an incompressible hyperelastic skeleton saturated by an incompressible fluid. The governing equations stem from mixture theory and the application is motivated by the study of interstitial fluid flow in brain tissue. The formulation is based on the adoption of an arbitrary Lagrangian–Eulerian (ALE) perspective. We focus on a flow regime in which inertia forces are negligible. Finally, the stability and convergence of the formulation is discussed, and numerical results demonstrate agreement with the theory.

97 MATHEMATICS AND COMPUTING↗

A comparison of classical mechanics models and finite element simulation of elastically tailored wing boxes

Structural tailoring concepts were developed to create wings with elastically produced camber for the purpose of increasing lift during takeoff conditions. Simple models based upon enhancements to the thin walled composite beam theory of Rehfield were developed to investigate prospects for elastic tailoring of the chordwise deformation of wing structures. The purpose here is to provide a comparison of the theoretical results with a finite element model for the bending method of producing camber. Finite element correlation studies were completed for two cases: a bonded unstiffened structural box, and a bolted unstiffened structural box. Results from these studies show an error of less than one percent for the bonded case and less than six percent for the bolted case in predicting camber curvature for the structural box. Examination of the results shows that the theory is very accurate for the cases studied and will provide an excellent basis for conducting further tailoring studies.

Rehfield, Lawrence W.↗

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

The DEAL.II finite element library: Design, features, and insights

DEAL.II is a state-of-the-art finite element library focused on generality, dimension-independent programming, parallelism, and extensibility. In this paper, we outline its primary design considerations and its sophisticated features such as distributed meshes, h p -adaptivity, support for complex geometries, and matrix-free algorithms. But DEAL.II is more than just a software library: It is also a diverse and worldwide community of developers and users, as well as an educational platform. We therefore also discuss some of the technical and social challenges and lessons learned in running a large community software project over the course of two decades.

97 MATHEMATICS AND COMPUTING↗

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗