Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element Method”

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 OSTI2025

Effects of Material Selection on the Remanufacturability of a Swashplate in an Axial Piston Pump

A swashplate is a critical component in an axial piston pump for motion transformation that undergoes cyclic loading. Traditionally, swashplates are not designed to have multiple lives of operation. According to industry experts, a remanufacturability analysis of swashplates suggests that the scrap rate is high. Among the primary reasons for this high scrap rate is insufficient material on the running surface, especially for subtractive recovery. Therefore, the selection of the proper material and its amount plays a significant role in the fatigue behavior, the remanufacturability of the swashplate, as well as its cost. Analyses addressing fatigue life, remanufacturability, and cost of the swashplate were carried out in this paper through a comparative study utilizing various materials. Recommendations on the material selection that balance the failure consequences, remanufacturability, and cost are provided to guide designers in making informed decisions for alternative design options.

42 ENGINEERING↗
DOE OSTI2023

Simulation studies of a full-ring, CZT SPECT system for whole-body imaging of 99m Tc and 177 Lu

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. With the goal of filling this technological gap, we have designed a cadmium zinc telluride (CZT) full-ring SPECT scanner instrumented with a broad-energy tungsten collimator. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT does not rely on the dual- and the triple-head rotational SPECT standard paradigm, enabling a larger effective area in each scan to increase the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized phantoms filled with 99m Tc and 177 Lu. Our SPECT design can resolve 7.9 mm rods for 99m Tc (140 keV) and 9.5 mm for 177 Lu (208 keV) in a hot-rod Derenzo phantom with a 3-min exposure and reach an image contrast of 78% for 99m Tc and 57% for 177 Lu using the NEMA IQ phantom with a 6-min exposure. Our modified scatter correction shows an improved contrast-recovery ratio compared to a standard correction. In this paper, we demonstrate the good performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99m Tc imaging of brain perfusion and 123 I imaging of dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design provides similar IQ and contrast to the commercial full-ring SPECT VERITON for 99m Tc. Regarding 177 Lu imaging of the 208 keV emissions, our design provides similar contrast to that of other state-of-the-art SPECTs with a significant reduction in exposure. In conclusion, the high sensitivity and extended energy range up to 250 keV makes our SPECT design a promising alternative for clinical imaging and theranostics of emerging radionuclides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

A Generalized Backward Euler algorithm for the numerical integration of a viscous breakage model

This paper discusses the formulation and the numerical performance of a fully implicit algorithm used to integrate a rate-dependent model defined within a breakage mechanics framework. For this purpose, a Generalized Backward Euler (GBE) algorithm has been implemented according to two different linearization strategies: The former is derived by a direct linearization of the constitutive equations, while the latter introduces rate effects through a consistency parameter. Additionally, the accuracy and efficiency of the GBE algorithm have been investigated by (1) performing material point analyses and (2) solving initial boundary value problems. In both cases, the overall performance of the underlying algorithm is inspected for a range of loading rates, thus simulating comminution from slow to fast dynamic problems. As the viscous response of the breakage model can be recast through a viscous nucleus function, the presented algorithm can be considered as a general framework to integrate constitutive equations relying on the overstress approach typical of Perzyna-like viscoplastic models.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Generalized radial-return mapping algorithm for anisotropic von Mises plasticity framed in material eigenspace

A computationally efficient integration algorithm for anisotropic plasticity is proposed, which is identified as a generalization of the radial-return mapping algorithm to anisotropy. The algorithm is based upon formulation within the eigenspace of a material anisotropy tensor associated with anisotropic quadratic von Mises (J 2 ) plasticity (also called Hill plasticity), for which it is shown to ensure that the flow rule remains associative, ie, the normality condition is satisfied. Extension of the algorithm to include anisotropic elasticity (anisotropic elastoplasticity) is further provided, made possible by the identification of a certain fourth-order material tensor dependent on both the elastic and plastic anisotropy. The derivation of the fully elastoplastically anisotropic algorithm involves further complexity, but the resulting algorithm is shown to closely resemble the purely plastically anisotropic one once the appropriate eigenspace is identified. The proposed generalized radial-return algorithm is compared to a classical closest-point projection algorithm, for which it is shown to provide considerable advantage in computational cost. In conclusion, the efficiency, accuracy, and robustness of the algorithm are demonstrated through various illustrative test cases and in the finite element simulation of Taylor impact tests on tantalum.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2024

A large deformation multiphase continuum mechanics model for shock loading of soft porous materials

A large deformation, coupled finite-element (FE) model is developed to simulate the multiphase response of soft porous materials subjected to high strain-rate loading. The approach is based on the theory of porous media (TPM) at large deformations. Simplifications to the one-dimensional regime studied in the numerical simulations follow. An overview of several different time integration schemes is presented for the purpose of solving the nonlinear dynamic coupled balance of momenta (mixture and fluid) and balance of mass of the mixture equations. Numerical examples are presented for (i) verification against closed-form analytical solutions assuming small loads, (ii) demonstrating large deformation effects at high strain-rate, and (iii) showing differences in deformations between a single-phase elastodynamics model with occluded compressible pore fluid and a multiphase poroelastodynamics model at high strain-rate. The multiphase model shows that the relative motion of the pore fluid significantly dampens the deformation response of the solid skeleton as compared to the single-phase model, and makes it possible to extract quantitative values for the stresses of the different constituents, thereby allowing one to form preliminary conclusions about the onset of damage in the solid skeleton. The novelty of the current work is developing a multiphase, large deformation, mixture theory numerical model for high strain-rate loading of soft porous materials. It was discovered that explicit, adaptive time-stepping Runge–Kutta schemes offer high accuracy at relatively low cost when compared to traditional implicit or explicit central difference time-stepping schemes for shock-like loadings. Here, shock viscosity is added to the mixture momentum balance equation to regularize the shock front, and a stabilization term is added to the mixture mass balance equation to stabilize equal order interpolation finite elements for the coupled finite element solution of multiphase materials.

Engineering↗
DOE OSTI2024

A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures

Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.

42 ENGINEERING↗
DOE OSTI2020

Linear viscoelastic characterization of electrically conductive adhesives used as interconnect in photovoltaic modules

Electrically conductive adhesives (ECAs) are incorporated into recent designs of photovoltaic (PV) modules and replace the traditional metallic solders as interconnects. This transition depicts a significant material change, and a proper understanding of the interconnects' mechanical response has not yet been established. However, such an understanding is necessary to (a) identify the driving forces for module degradation and failure, (b) allow for module design optimization, and (c) enable accurate lifetime predictions. This study summarizes the framework for the mechanical materials characterization and modeling of ECAs for PV applications. Only high-fidelity material models are able to capture the rate and temperature dependency of the ECA interconnect and allow for accurate modeling of the materials response. Furthermore, a linear viscoelastic representation is found to describe the mechanical response of the ECAs sufficiently well. The effects of curing conditions and environmental exposure are investigated, and material models for a variety of ECAs are reported and prepared for the use in numerical simulations. Overall, a finite element simulation of a generic submodel of a shingled cell module is used to highlight the need for high-fidelity material models and demonstrates the error made in the predicted stress states by using less sophisticated models.

14 SOLAR ENERGY↗
DOE OSTI2023

Multiscale numerical investigation of ratchet growth damage effects in PBX 9502

This paper presents results of numerical experiments conducted on the high explosive PBX 9502 to investigate how recently observed grain-scale damage mechanisms of ratchet growth affect uniaxial compression measurements. Simulations are multiscale in the sense of directly resolving grains, pores, cracks, and grain-interfaces based upon scanning electron microscope (SEM) images of damaged and undamaged samples. The combined finite-discrete element method (FDEM) is utilized to resolve both grain-scale microfracture and elastoplastic deformation of solid grains. Pristine (undamaged) and damaged microstructures are compared in simulation of unconfined compression tests of the same material from the literature. Here, the simulation results show the observed microscale mechanisms of damage, specifically microfracture predominantly around and sometimes through grains and crack-associated pore growth, can well-explain the effective degradation of strength and stiffness observed in the laboratory measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Thermomechanical Modeling of High-Temperature Bonded Interface Materials

In an automotive power electronics package, bonded interface materials play a critical role in providing a smooth heat dissipation path from the device to the coolant. However, these materials are susceptible to failure when the package gets exposed to repeated changes in its thermal environment and internal temperature variations. Under these thermal conditions, the coefficient of thermal expansion mismatch between the different component layers results in thermally induced stresses and consequently the delamination of the material from its adjoining surfaces or crack initiation and propagation within the material. Thermomechanical modeling offers a path to study the complex nature of these bonded materials, guides their design and selection, and, in conjunction with experimental results, plays an important role in predicting their lifetime. In this chapter, two modeling strategies - similar in approach and implementation but different in the underlying theory - are discussed in the context of high-temperature bonded materials. In addition, a short review of sintered silver is provided, as it is widely seen as a potentially promising and reliable bonded material for high-temperature applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design

Elastic meta-materials are those whose unique properties come from their micro-architecture, rather than, e.g., from their chemistry. The introduction of such architecture, which is increasingly able to be fabricated due to advances in additive manufacturing, expands the design domain and enables improved design, from the most complex multi-physics design problems to the simple compliance design problem that is our focus. Unfortunately, concurrent design of both the micro-scale and the macroscale is computationally very expensive when the former can vary spatially, particularly in three dimensions. Instead, we provide simple, accurate surrogate models of the homogenized linear elastic response of the isotruss, the octet truss, and the ORC truss based on high-fidelity continuum finite element analyses. These surrogate models are relatively accurate over the full range of relative densities, in contrast to analytical models in the literature, which we show lose accuracy as relative density increases. The surrogate models are also simple to implement, which we demonstrate by modifying Sigmund’s 99-line code to solve a three-dimensional, multiscale compliance design problem with spatially varying relative density. Here, we use this code to generate examples in both two and three dimensions that illustrate the advantage of elastic meta-materials over structures with a single length scale, i.e., those without micro-architectures.

42 ENGINEERING↗
DOE OSTI2018

Optimal Control of a Degenerate PDE for Surface Shape

Controlling the shapes of surfaces provides a novel way to direct self-assembly of colloidal particles on those surfaces and may be useful for material design. This motivates the investigation of an optimal control problem for surface shape in this paper. Specifically, we consider an objective (tracking) functional for surface shape with the prescribed mean curvature equation in graph form as a state constraint. The control variable is the prescribed curvature. We prove existence of an optimal control, and using improved regularity estimates, we show sufficient differentiability to make sense of the first order optimality conditions. This allows us to rigorously compute the gradient of the objective functional for both the continuous and discrete (finite element) formulations of the problem. Numerical results are shown to illustrate the minimizers and optimal controls on different domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Examination of Mean Stress Calculation Approaches in Rock Mechanics

The mean stress, as a fundamental statistical property of a group of stress data, is essential for stress variability characterisation. However, currently in rock mechanics, the mean stress is customarily and erroneously calculated by separately averaging the principal stress magnitudes and orientations. In order to draw the attention of our community to the appropriate approach for stress variability characterisation, here we compare the customary scalar/vector mean with that obtained by the mathematically rigorous tensorial approach—the Euclidean mean. Calculation of mean stress using both a small group of actual in situ stress measurement results and a large group of simulated stress data (obtained using the combined finite– discrete element method, FEMDEM) demonstrates that the two approaches yield different results. Further investigation of these results shows that the scalar/vector approach may yield non-unique and non-orthogonal mean principal stresses, and these may deviate significantly from those of the Euclidean mean. Our calculations and comparisons reveal that the scalar/ vector approach is deficient because it processes the principal stress magnitudes and orientations separately as independent quantities and ignores the connection between them. Conversely, the tensorial approach appropriately averages the tensors that simultaneously carry not only the information of stress magnitudes and orientations, but also the inherent relations between them. Therefore, arbitrarily using scalar/vector mean stress of in situ stress measurements as input in further rock engineering analyses may yield significantly erroneous results. Finally, we advise that stress data should be statistically processed in a tensorial manner using tensors referred to a common Cartesian coordinate system.

58 GEOSCIENCES↗
DOE OSTI2018

Impacts of hydrogeological characteristics on groundwater-level changes induced by earthquakes

Changes in groundwater level during earthquakes have been reported worldwide. In this study, field observations of co-seismic groundwater-level changes in wells under different aquifer conditions and sampling intervals due to near-field earthquake events in Taiwan are presented. Sustained changes, usually observed immediately after earthquakes, are found in the confined aquifer. Oscillatory changes due to the dynamic strain triggered by passing earthquake waves can only be recorded by a high-frequency data logger. While co-seismic changes recover rapidly in an unconfined aquifer, they can sustain for months or longer in a confined aquifer. Three monitoring wells with long-term groundwater-level data were examined to understand the association of co-seismic changes with local hydrogeological conditions. The finite element software ABAQUS is used to simulate the pore-pressure changes induced by the displacements due to fault rupture. The calculated co-seismic change in pore pressure is related to the compressibility of the formation. The recovery rate of the change is rapid in the unconfined aquifer due to the hydrostatic condition at the water table, but slow in the confined aquifer due to the less permeable confining layer. Fracturing of the confining layer during earthquakes may enhance the dissipation of pore pressure and induce the discharge of the confined aquifer. The study results indicated that aquifer characteristics play an important role in determining groundwater-level changes during and after earthquakes.

58 GEOSCIENCES↗
DOE OSTI2018

Hydro-mechanically coupled finite-element analysis of the stability of a fractured-rock slope using the equivalent continuum approach: a case study of planned reservoir banks in Blaubeuren, Germany

Construction works associated with the building of reservoirs in mountain areas can damage the stability of adjacent valley slopes. Seepage processes caused by the filling and drawdown operations of reservoirs also affect the stability of the reservoir banks over time. The presented study investigates the stability of a fractured-rock slope subjected to seepage forces in the lower basin of a planned pumped-storage hydropower (PSH) plant in Blaubeuren, Germany. The investigation uses a hydro-mechanically coupled finite-element analyses. For this purpose, an equivalent continuum model is developed by using a representative elementary volume (REV) approach. To determine the minimum required REV size, a large number of discrete fracture networks are generated using Monte Carlo simulations. These analyses give a REV size of 28 × 28 m, which is sufficient to represent the equivalent hydraulic and mechanical properties of the investigated fractured-rock mass. The hydro-mechanically coupled analyses performed using this REV size show that the reservoir operations in the examined PSH plant have negligible effect on the adjacent valley slope.

58 GEOSCIENCES↗
DOE OSTI2018

Experimental Validation of the Earthquake Resistance of a Full-Scale Model of RBMK-1000 Graphite Stack with Post-Repair Characteristics

In providing justification for the operational safety of RBMK-1000 power units with graphite stack distortion, earthquake resistance tests were performed on a full-scale model of RBMK graphite stack under seismic platform conditions. The tests were performed for three mutually perpendicular axes of the model successively for seismic loads corresponding to magnitudes 6, 7 and 8 on the MSK-64 scale. The tests confirmed the earthquake resistance of the graphite stack for prescribed loads. The results were used to verify the FEMGR code (Kurchatov Institute) and the LEGAK-DF finite-element software system (RFYaTs – VNIIEF) for purposes of computational validation of the earthquake resistance of the post-repair restored RMBK-1000 core.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Analysis of the Resonance Phenomena in the Operation of the Lower Bureya HPP During an Earthquake

To justify that hydroelectric power plants operating in seismic areas are safe, it is necessary to analyze the resonance phenomena in the generating unit. The analysis was performed using a finite-element model of the powerhouse-foundation-equipment system of the Lower Bureya HPP. The resonance phenomena occurring when the disturbance frequency at the elevation of the stator supports is close to the first natural frequency and speed of the water turbine were studied.

42 ENGINEERING↗
DOE OSTI2018

Study of Dissipative Losses in AC-Biased Mo/Au Bilayer Transition-Edge Sensors

We are developing kilo-pixel arrays of transition-edge sensors (TESs) for the X-ray Integral Field Unit on ESA’s Athena observatory. Previous measurements of AC-biased Mo/Au TESs have highlighted a frequency-dependent loss mechanism that results in broader transitions and worse spectral performance compared to the same devices measured under DC bias. In order to better understand the nature of this loss, we are now studying TES pixels in different geometric configurations. We present measurements on devices of different sizes and with different metal features used for noise mitigation and X-ray absorption. Our results show how the loss mechanism is strongly dependent upon the amount of metal in close proximity to the sensor and can be attributed to induced eddy current coupling to these features. We present a finite element model that successfully reproduces the magnitude and geometry dependence of the losses. Using this model, we present mitigation strategies that should reduce the losses to an acceptable level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗