Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite element simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Thermal Energy Storage for Dispatchable Geothermal Power

Dispatchable generation refers to sources of electricity that can be turned on or off, or can adjust their power output accordingly to an order. Geothermal plants are usually used for base-load power rather than dispatchable power. Flexibility is a challenge for a geothermal plant because the heat source does not naturally respond quickly to changes in demand. The current state-of-the-art approach to address this challenge is to maintain the geothermal fluid flows at relatively steady rates while providing a bypass around the generation equipment as needed. Under part load, some of the geothermal fluid is pumped to the surface, bypasses the generation equipment, and is re-injected to the ground without extracting any useful enthalpy. This approach is robust but necessarily incurs high parasitic power draws at part load due to the constant, full-flow pumping power requirements. There are even more challenges and unknowns for steam and flash geothermal plants, which were the focus for this project. As a process improvement to this existing configuration of a dispatchable geothermal plant, we examined the use of a phase-change thermal-energy storage (TES) unit that will provide some “inertia” against temperature fluctuations. The TES phase change material (PCM) could be in the form of eutectic molten nitrate salts, organic sugar alcohols, or the like based on the required temperature. This configuration allowed the power plant to meet the ramp rate requirements while allowing the geothermal resource to drawn down at the duty-cycle-averaged extraction rate rather than the full-rated-output extraction rate. The project team performed an experimental and computational evaluation of a peak shaving technical approach that used thermal storage to switch between two nominal operating points: full rated power and 50% power. In the initial requirements definition phase, the team decided to examine the PCM/TES-enhanced geothermal plant in the context of a 40-MW two-stage flash plant operating in mid-Southern California. Steady-state process simulations were used to understand the equipment sizing and the energetics. Dynamic process simulations were used to understand the response rate of the system to changes in requested power output from the utility. The technical approach was also guided by a techno-economic analysis (TEA) and a high-level analysis of the reservoir. Recognizing that performance of the thermal storage was a key risk element, high-fidelity finite-element simulations of the PCM-TES heat exchanger were performed to better understand its performance and temporal response. These simulations were validated with a subscale PCM-TES experimental apparatus that mimicked the geometry of the full-scale system. The key outcomes of this Phase I project were that: 1) the PCM-TES heat exchanger met the target for gravimetric energy storage density, 2) the dispatchable geothermal plant met its target for response rate, 3) the peak shaving approach reduces the number of production wells and does not disturb the wellfield, and 4) the incentives for properly remunerating dispatchable operation were identified.

15 GEOTHERMAL ENERGY↗
NASA NTRS1997

Ultrasonic Nondestructive Characterization of Adhesive Bonds

Qualitative measurements of adhesion or binding forces can be accomplished, for example, by using the reflection coefficient of an ultrasound or by using thermal waves (Light and Kwun, 1989, Achenbach and Parikh, 1991, and Bostrom and wickham, 1991). However, a quantitative determination of binding forces is rather difficult. It has been observed that higher harmonics of the fundamental frequency are generated when an ultrasound passes through a nonlinear material. It seems that such non-linearity can be effectively used to characterize the bond strength. Several theories have been developed to model this nonlinear effect (Adler and Nagy, 1991; Achenbach and Parikh, 1991; Parikh and Achenbach, 1992; and Hirose and Kitahara, 1992; Anastasi and Roberts, 1992). Based on a microscopic description of the nonlinear interface binding force, a quantitative method was presented by Pangraz and Arnold (1994). Recently, Tang, Cheng and Achenbach (1997) made a comparison between the experimental and simulated results based on this theoretical model. A water immersion mode-converted shear wave through-transmission setup was used by Berndt and Green (1997) to analyze the nonlinear acoustic behavior of the adhesive bond. In this project, the nonlinear responses of an adhesive joint was investigated through transmission tests of ultrasonic wave and analyzed by the finite element simulations. The higher order harmonics were obtained in the tests. It is found that the amplitude of higher harmonics increases as the aging increases, especially the 3dorder harmonics. Results from the numerical simulation show that the material nonlinearity does indeed generate higher order harmonics. In particular, the elastic-perfect plastic behavior generates significant 3rd and 5th order harmonics.

Qu, Jianmin↗
DOE OSTI2025

Simulations of the fault current limiting operation of a long-length REBCO CORC ® superconducting cable cooled by helium gas

Conductor-on-round-core (CORC®) cables composed of rare-earth-barium-copper-oxide high-temperature superconducting (HTS) tapes are of great interest for power transmission applications due to their many advantages such as high power density, light weight, and low loss. Closed circulation loops of cryogenic helium gas can be used to cool HTS cables down to low temperatures to significantly improve their current-carrying capacity. Coupled circuit-electromagnetic-thermal finite element simulations implemented in the COMSOL Multiphysics package were developed, validated, and then used for simulating the fault current limiting (FCL) performance and the cooling processes of an 8-layer CORC® cable cooled with a flow of cryogenic helium gas. In the simulations, the temperature dependence of the electrical and thermal properties of all component materials is implemented for improved accuracy. To overcome computational challenges caused by the considerable difference in geometrical scales (i.e. few-µm-thick HTS layers versus 10 m-long HTS cable), the model is divided into two separate simulations. The first simulation is performed on the transverse cross-section of the cable to calculate the electric field, heating power and temperature rise in each component of a CORC® cable during FCL operation. The heating power calculated in the first simulation is transferred to the second model to simulate the cooling of a 10 m-long cable after the fault is cleared. The effect of the helium gas flow rate on the cooling process is also investigated to develop strategic approaches for optimizing cooling systems for HTS cables with FCL capability. The simulations indicated that a 40 ms fault with a voltage drop of 20 V m −1 along the cable can result in a temperature increase from 60 K to about 165 K inside the cable, and it takes about 500 s to cool the cable back to nearly 60 K with a flow of cold helium gas at a rate of 5 g s −1 .

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Dopant Compensation within the Intrinsic Poly-Si Isolation Region in Poly-Si/SiOx Passivated IBC Si Solar Cells

We report on the effect of dopant compensation within intrinsic poly-Si regions between p- and n-type fingers of poly-Si/SiOx passivated interdigitated back contact (IBC) solar cells using intrinsic poly-Si as the isolation region between the doped poly-Si fingers. First, we show that dopants from the doped fingers contaminate the intrinsic gap, resulting in doping of the entire intrinsic gap and overlap of the dopant tails from each finger. Next, we show that despite this doping across the gap, shunting between the doped fingers does not occur. We show that this is a result of trap-assisted compensation creating a highly resistive intrinsic region, preventing shunt. We simulate shunt resistance across the gap based on local carrier concentration and deep trap density. We show that trap defects within the poly-Si enhance compensation between the dopant tails. We experimentally confirm these predictions by scanning spreading resistance microscopy of the gap showing ~20 um domain with resistivity ~10^7 ..omega.. cm. Additionally, Kelvin probe force microscopy are compared to finite element simulations which result in the same approximate shape for potential profile, indicating diode behavior across the isolation region. These results demonstrate the powerful effect that trap defects have within the poly-Si isolation region and suggest that precision patterning is not as essential as once thought.

dopant compensation↗
DOE OSTI2019

Numerical Investigation of the Effect of Rolling on the Localized Stress and Strain Induction for Wire + Arc Additive Manufactured Structures

Cold rolling can be used in-process or post-process to improve microstructure, mechanical properties and residual stress in directed-energy-deposition techniques, such as the high deposition rate wire + arc additive manufacturing (WAAM) process. Finite element simulations of the rolling process are employed to investigate the effect of rolling parameters, in particular rolling load and roller profile radius on the residual stress field as well as plastic strain distribution for the profiled roller. The results show the response to rolling of commonly used structural metals in WAAM, i.e., AA2319, S335JR steel and Ti-6Al-4V, taking into account the presence of residual stresses. The rolling load leads to changes in the location and the maximum value of the compressive residual stresses, as well as the depth of the compressive residual stresses. However, the roller profile radius only changes the maximum value of these compressive residual stresses. Changing the rolling load influences the equivalent plastic strain close to the top surface of the wall as well as in deeper areas, whereas the influence of the roller profile radius is negligible. The plastic strain distribution is virtually unaffected by the initial residual stresses prior to rolling. Finally, design curves were generated from the simulations for different materials, suggesting ideal rolling load and roller profile combinations for microstructural improvement requiring certain plastic strains at a specific depth of the additive structure.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Application of the Quadrupole Method for Simulation of Passive Thermography

Passive thermography has been shown to be an effective method for in-situ and real time nondestructive evaluation (NDE) to measure damage growth in a composite structure during cyclic loading. The heat generation by subsurface flaw results in a measurable thermal profile at the surface. This paper models the heat generation as a planar subsurface source and calculates the resultant temperature profile at the surface using a three dimensional quadrupole. The results of the model are compared to finite element simulations of the same planar sources and experimental data acquired during cyclic loading of composite specimens.

Winfree, William P.↗
DOE OSTI2021

Learning constitutive relations using symmetric positive definite neural networks

In this work, we present a new neural-network architecture, called the Cholesky-factored symmetric positive definite neural network (SPD-NN), for modeling constitutive relations in computational mechanics. Instead of directly predicting the stress of the material, the SPD-NN trains a neural network to predict the Cholesky factor of the tangent stiffness matrix, based on which the stress is calculated in incremental form. As a result of this special structure, SPD-NN weakly imposes convexity on the strain energy function, satisfies the second order work criterion (Hill's criterion) and time consistency for path-dependent materials, and therefore improves numerical stability, especially when the SPD-NN is used in finite element simulations. Depending on the types of available data, we propose two training methods, namely direct training for strain and stress pairs and indirect training for loads and displacement pairs. We demonstrate the effectiveness of SPD-NN on hyperelastic, elasto-plastic, and multiscale fiber-reinforced plate problems from solid mechanics. The generality and robustness of SPD-NN make it a promising tool for a wide range of constitutive modeling applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

A new model for the effective thermal conductivity of polycrystalline solids

We introduce a novel model for the effective thermal conductivity of polycrystalline solids based on the thin-interface description of grain boundaries (GBs). In contrast to existing models, our new model treats a GB as an autonomous “phase” with its own thermal conductivity. The Kapitza resistance/conductance of a thin interface is then derived in terms of the interface thermal conductivity and width. In turn, the effective thermal conductivity of polycrystals is derived in terms of grain size, grain and GB conductivities, and GB width. This treatment allows the model to simulate the change of the Kapitza resistance/conductance with segregation/doping, GB structure/phase transition, or GB decohesion. Moreover, since the model assumes a finite width for GBs, it is expected to give better predictions than its sharp-interface-based counterparts for nanoscale grains. The predictions of the new model deviate from the corresponding ones from existing models by 1%–100% as the grain size approaches the GB width. High-fidelity finite-element simulations were conducted to validate the predictions of the new model. These simulations proved the higher accuracy of the new model. We also discuss how to generalize this treatment to other types of interfaces in heterogeneous materials. The advantages and limitations of the new model are summarized, and some future directions are highlighted.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Multistage Axial Hydro-Forging Sequence: A New Forming Approach for Manufacturing of Double-Stepped Tubes

To manufacture the large-expansion-ratio double-stepped tube without thinning, a new approach known as multistage axial hydro-forging sequence has been proposed in this study. The whole forming process is divided into five stages. Compared with traditional tube hydroforming, the thickness can be increased by axial compression deformation to reduce the thinning ratio and enhance the forming capacity of the tube in hydro-forging stage and double-stepped hydro-forging stage, respectively. The proposed multistage axial hydro-forging sequence process was implemented by experiments along with the finite element simulation to validate its feasibility. The suitable loading path to prevent defects such as wrinkle and thinning was discussed. The results showed that double-stepped tube with large expansion ratio can be produced successfully with movable dies’ design. There is no thinning in the final tube, particularly in all stepped area, and a better thickness distribution could be achieved. As a result, the proposed multistage axial hydro-forging sequence can be a feasible forming approach for the fabrication of large-expansion-ratio double-stepped tubes with high reliability.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Modular machine learning-based elastoplasticity: Generalization in the context of limited data

The development of highly accurate constitutive models for materials that undergo path-dependent processes continues to be a complex challenge in computational solid mechanics. Challenges arise both in considering the appropriate model assumptions and from the viewpoint of data availability, verification, and validation. Recently, data-driven modeling approaches have been proposed that aim to establish stress-evolution laws that avoid user-chosen functional forms by relying on machine learning representations and algorithms. However, these approaches not only require a significant amount of data but also need data that probes the full stress space with a variety of complex loading paths. Furthermore, they rarely enforce all necessary thermodynamic principles as hard constraints. Hence, they are in particular not suitable for low-data or limited-data regimes, where the first arises from the cost of obtaining the data and the latter from the experimental limitations of obtaining labeled data, which is commonly the case in engineering applications. In this work, we discuss a hybrid framework that can work on a variable amount of data by relying on the modularity of the elastoplasticity formulation where each component of the model can be chosen to be either a classical phenomenological or a data-driven model depending on the amount of available information and the complexity of the response. The method is tested on synthetic uniaxial data coming from simulations as well as cyclic experimental data for structural materials. The discovered material models are found to not only interpolate well but also allow for accurate extrapolation in a thermodynamically consistent manner far outside the domain of the training data. This ability to extrapolate from limited data was the main reason for the early and continued success of phenomenological models and the main shortcoming in machine learning-enabled constitutive modeling approaches. Training aspects and details of the implementation of these models into Finite Element simulations are discussed and analyzed.

42 ENGINEERING↗
DOE OSTI2021

Assessment Approach to Advanced Fuel Models

This document constitutes completion of the NEAMS milestone M3MS-21IN0201016, which is titled: Assessment of advanced fuel improvements. In this report we present: (1) simulations of separate effects creep tests for uranium alloy and HT-9; (2) examples of standardizing and updating metallic fuel assessments to include the latest models and key output figures of merit, which is coordinated with LANL to better evaluate advanced models; (3) examples of how the FIPD database at ANL is integrated into metallic fuel assessment cases; and finally (4) documentation of fuel-specific figures of merit and the BISON review process for assessments with best practices for preparing and reviewing finite element simulations in the appendix.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRSDate not supplied

Preliminary 3D Simulations of Anomalous Thermocouple Coupling to Surface Electric Fields

Anomalous measurements from subsurface thermocouples have been observed during re-entry of the Orion capsule (EFT-1 plugs). These measurements are considered anomalous because the amplitude and frequency of the inferred temperature modulations are too large to be due to plausible heating processes. Consequently, it is believed the anomalies are caused by electrical interference generated by the re-entry plasma. Quantifying these effects requires several steps: modeling the external plasma flows, modeling the plasma-surface interaction, and finally modeling the surface-thermocouple coupling. This paper addresses this last step using 3D finite element simulations of a wire embedded in a weakly conductive medium.

thermocouple anomaly↗
DOE OSTI2017

The in-depth residual strain heterogeneities due to an indentation and a laser shock peening for Ti-6Al-4V titanium alloy

In this paper, heterogeneity of the through-thickness residual strain due to the laser shock peening (LSP) in comparison with that due to the indentation was studied in Ti-6Al-4V alloy samples. The latter is almost a quasi-static process while the former features extremely high strain-rate deformation. The synchrotron based high-energy X-ray diffraction was employed to investigate the through-thickness residual strain distribution. The studied two samples after the two processing procedures share the following features: (i) the pressure affected depths are both ~2 mm and (ii) the largest magnitudes of the compressive residual strains parallel to the surface are ~4,000 με. However, the pit depth for the indentation sample is ~9 times larger than that for the LSP. The position featuring the largest magnitude of the compressive residual strain is in the sub-surface for the indentation while it is in the surface for the LSP. Results of the elastic-visco-plastic finite element simulation for the indentation indicate that the position featuring the maximum accumulative plastic shear as defined in this paper corresponds to the location with the largest magnitude of compressive residual strain. To validate this finding, full width at half maximum (FWHM) of the X-ray diffraction profile, which is proportional to the level of the plastic deformation, is also studied. Finally, it is found that positions with the largest FWHM indeed correspond to the largest magnitude of compressive residual strain for both tests.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Quantifying leakage fields at ionic grain boundaries using off-axis electron holography

The electrical properties of interfaces in semiconductors and ionic conductors are immensely important in a wide range of applications. Electron holography is ideally suited for the direct measurement of the electrostatic potential of such interfaces. A key challenge with this approach is the contribution of the leakage field from the sample to the observed electron phase shift. This leakage field cannot be a priori independently determined and can cause an overestimation of the phase shift. In this work, we use finite element simulations to compute the three-dimensional electrostatic potential in the vicinity of an interface associated with a given interfacial charge density distribution. We then evaluate the predicted phase shift and demonstrate that the leakage field strongly affects the recovery of the projected interface potential. From the difference between the true potential and uncorrected, recovered potential, we propose a method to correct for this effect. We then demonstrate the application of this methodology to the analysis of experimental off-axis electron holography data acquired from the grain boundaries in lightly doped ceria.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamic Tensile Characterization of Thin-Sheet Brittle Metallic Materials

Refractory metals are favorable materials in applications where high strength and ductility are needed at elevated temperatures. In some cases, operating temperatures may be near the melting point of the material. However, as temperature drops, refractory metals typically undergo a significant mechanical response change - ductile-to-brittle transition. These materials may be subjected to high strain rate loading at an ambient temperature state, such as an impact or crash. Knowledge of the high rate material properties are essential for design as well as simulation of impact events. The high rate stress-strain behavior of brittle metallic materials at ambient temperature is rarely studied because of experimental challenges, particularly when failure is involved. Failure typically occurs within the non-gage section of the material, which invalidates any collected stress-strain information. Here, a method to determine a specimen geometry which will produce failures in the gage section is presented. Pure tungsten in thin-sheet form was used as a trial material to select a specimen geometry for high rate Kolsky tension bar experiments. A finite element simulation was conducted to derive a strain correction for more accurate results. The room temperature stress-strain behavior of pure tungsten at a strain rate of 24 s –1 is presented. The outcome of this experimental technique can be applied to other brittle materials for dynamic tensile characterization.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Microstructure-Based Estimation of Strength and Ductility Distributions for \(\alpha +\beta \) Titanium Alloys

Titanium alloys are processed to develop a wide range of microstructure configurations and therefore material properties. While these properties are typically measured experimentally, a framework for property prediction could greatly enhance alloy design and manufacturing. In this work, a microstructure-sensitive framework is presented for the prediction of strength and ductility as well as estimates of the bounds in variability for these properties. The framework explicitly considers distributions of microstructure via new approaches for instantiation of structure in synthetic samples. The parametric evaluation strategy, including the finite element simulation package FEpX, is used to create and test virtual polycrystalline samples to evaluate the variability bounds of mechanical properties in Ti-6Al-4V. Critical parameters for the property evaluation framework are provided by measurements of single crystal properties and advanced characterization of microstructure and slip system strengths in 2D and 3D. Property distributions for yield strength and ductility are presented, along with the validation and verification steps undertaken. Comparisons between strain localization and slip activity in virtual samples and in experimental grain-scale strain measurements are also discussed.

36 MATERIALS SCIENCE↗
NASA NTRS2018

Probabilistic Sensitivity Analyses Incorporating Surrogate Models for the Structural Response of a Representative Inflatable Space Structure

This study concentrates on results from probabilistic analysis and numerical simulation tools to identify parameter sensitivities and assess surrogate model suitability for a novel inflatable airlock concept, specifically the Non-Axisymmetric Inflatable Pressure Structure (NAIPS) under NASA's Minimalistic Advanced Softgoods Hatch (MASH) Program. The current studies extended the demonstration of probabilistic analysis tools to identify parameter sensitivities by: incorporating nonlinear material load-strain property curves and orthotropic representation of the fabric; and exploring the use of surrogate models to facilitate rapid evaluations as new information becomes available, as well as enable optimizations. The airlock model has been simplified to focus on the end dome section. Brief overviews of the finite element model and associated surrogate models are provided along with the probabilistic sensitivity analysis approaches. Responses of interest include the sensitivities of various loads in the structural elements of the NAIPS to material properties, cord lengths, inflation pressure and friction between softgoods components. The use of surrogate models was explored to facilitate rapid sensitivity evaluations as parameter ranges change. The completion of the sensitivity studies improved understanding of the dependence of load responses to several uncertain parameters and confidence in the ability to use surrogate models to represent the finite element simulations.

Lyle, Karen H.↗
DOE OSTI2024

Effect of Part Size, Displacement Rate, and Aging on Compressive Properties of Elastomeric Parts of Different Unit Cell Topologies Formed by Vat Photopolymerization Additive Manufacturing

Due to its ability to achieve geometric complexity at high resolution and low length scales, additive manufacturing (AM) has increasingly been used for fabricating cellular structures (e.g., foams and lattices) for a variety of applications. Specifically, elastomeric cellular structures offer tunability of compliance as well as energy absorption and dissipation characteristics. However, there are limited data available on compression properties for printed elastomeric cellular structures of different designs and testing parameters. In this work, the authors evaluate how unit cell topology, part size, the rate of compression, and aging affect the compressive response of polyurethane-based simple cubic, body-centered, and gyroid structures formed by vat photopolymerization AM. Finite element simulations incorporating hyperelastic and viscoelastic models were used to describe the data, and the simulated results compared well with the experimental data. Of the designs tested, only the parts with the body-centered unit cell exhibited differences in stress–strain responses at different part sizes. Of the compression rates tested, the highest displacement rate (1000 mm/min) often caused stiffer compressive behavior, indicating deviation from the quasi-static assumption and approaching the intermediate rate response. The cellular structures did not change in compression properties across five weeks of aging time, which is desirable for cushioning applications. This work advances knowledge on the structure–property relationships of printed elastomeric cellular materials, which will enable more predictable compressive properties that can be traced to specific unit cell designs.

36 MATERIALS SCIENCE↗