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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Dynamically morphing microchannels in liquid crystal elastomer coatings containing disclinations

Liquid crystal elastomers (LCEs) hold a major promise as a versatile material platform for smart soft coatings, since their orientational order can be predesigned to program a desired dynamic profile. In this work, we introduce temperature-responsive dynamic coatings based on LCEs with arrays of singular defects-disclinations that run parallel to the surface. The disclinations form in response to antagonistic patterns of the molecular orientation at the top and bottom surfaces, imposed by the plasmonic mask photoalignment. Upon heating, an initially flat LCE coating develops linear microchannels located above each disclination. The stimulus that causes a non-flat profile of LCE coatings upon heating is the activation force induced by the gradients of molecular orientation around disclinations. To describe the formation of microchannels and their thermal response, we adopt a Frank-Oseen model of disclinations in patterned director field and propose a linear elasticity theory to connect the complex spatially-varying molecular orientation to the displacements of the LCE. The thermo-responsive surface profiles predicted by the theory and by the finite element modeling are in good agreement with the experimental data; in particular, higher gradients of molecular orientation produce a stronger modulation of the coating profile. The elastic theory and the finite element simulations allow us to estimate the material parameter that characterizes the elastomer coating’s response to the thermal activation. The disclinations containing LCEs show potential as soft dynamic coatings with a predesigned responsive surface profile.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Effect of lattice orientation on compressive properties of selective laser sintered nylon lattice coupons

As AM lattices become more popular in medical devices, it is important to consider how lattice design parameters affect the mechanical integrity and performance of a device. This research investigated the effect of lattice orientation on the compressive mechanical response of five common lattice geometries: Hexagonal Honeycomb (Hex), Diamond, Voronoi Tessellation Method (VTM), sheet-based Gyroid, and Face-Centered-Cubic (FCC). Samples were tested at two relative densities and eight orientations, printed in Nylon 12 (PA2200) on an EOS P396. The mechanical response was compared to simulated finite element analysis (FEA) response for Hex and FCC lattices. Hex lattices displayed significant orientation dependance. Diamond lattices demonstrated some orientation dependance when at a lower relative density. Gyroid, FCC, and VTM all displayed minimal orientation dependance. Increasing lattice relative density appeared to reduce anisotropy among tested orientations of all lattice geometries except Hex. FEA simulations were able to produce a trend response that matched the mechanical response trends based on the investigated orientations. Furthermore, these results support the need for multi axial lattice mechanical test considerations. It also bolsters the need for lattice FEA models to have a sufficient comparator when identifying “worst-case” lattice test scenarios.

36 MATERIALS SCIENCE↗

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↗

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 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 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↗

DOE OSTI2023

Residual Stresses Affect Cell Fragment Movement

Predictive modeling tools such as the finite element method can be of tremendous help in assessing the reliability and long-term performance of photovoltaic modules. In order to obtain accurate results, the proper modeling of materials and manufacturing processes are of utmost importance. Module fabrication introduces thermo-mechanical stresses inside the module laminate, which need to be accounted for as residual stresses in finite element simulations. We found that cell fragment movement and crack opening displacements of fractured silicon cells within modules are affected by those residual stresses. Cell cracking remains a challenging topic in assessing the reliability and durability of damaged modules. Hence, accurately quantifying the separation and movement between cell fragments creates the foundation for establishing reliable lifetime and performance assessments of fractured silicon modules. Here, we present a modeling approach that uses upper and lower bounds to accurately account for the residual stresses introduced by the module lamination process. We designed a four-point flexure coupon test of a laminated, fractured silicon strip to validate our numerical results and found good agreement between our modeling methodology and the experimental data. Finally, we discuss the implications of the residual stresses on the normal crack opening and metallization wear-out of fractured silicon cells.

14 SOLAR ENERGY↗

DOE OSTI2020

Intrinsic thermal interfacial resistance measurement in bonded metal–polymer foils

Heat conduction through bonded metal–polymer interfaces often limits the overall heat transfer in electronic packaging, batteries, and heat recovery systems. To design the thermal circuit in such systems, it is essential to measure the thermal interfacial resistance (TIR) across ~1 µm to 100 µm junctions. Previously reported TIR of metal–polymer junctions utilize ASTM E1530-based two-block systems that measure the TIR by applying pressure across the interface through external heating and cooling blocks. Here, we report a novel modification of the ASTM-E1530 technique that employs integrated heaters and sensors to provide an intrinsic TIR measurement of an adhesively bonded metal–polymer junction. We design the measurement technique using finite element simulations to either passively suppress or actively compensate the lateral heat diffusion through the polymer, which can minimize the systematic error to ≲5%. Through proof-of-concept experiments, we report the TIR of metal–polymer interfaces made from DuPont’s Pyralux double-side copper-clad laminates, commonly used in flexible printed circuit boards. Our TIR measurement errors are <10%. We highlight additional sources of errors due to non-idealities in the experiment and discuss possible ways to overcome them. Our measurement technique is also applicable to interfaces that are electrically insulating such as adhesively joined metal–metal junctions and sputter-coated or welded metal–polymer junctions. Altogether, the technique is capable of measuring TIR ≳10 –5 m 2 KW –1 in bonded metal–polymer foils and can be tailored for in situ measurements in flexible electronics, circuit packaging, and other hybrid metal–polymer systems.

47 OTHER INSTRUMENTATION↗