SEARCH · Engineering Papers

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

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 OSTI2014

Understanding Atom Probe Tomography of Oxide-Supported Metal Nanoparticles by Correlation with Atomic Resolution Electron Microscopy and Field Evaporation Simulation

Metal-dielectric composite materials, specifically metal nanoparticles supported on or embedded in metal oxides, are widely used in catalysis. The accurate optimization of such nanostructures warrants the need for detailed three-dimensional characterization. Atom probe tomography is uniquely capable of generating sub-nanometer structural and compositional data with part-per-million mass sensitivity, but there are reconstruction artifacts for composites containing materials with strongly differing fields of evaporation, as for oxide-supported metal nanoparticles. By correlating atom probe tomography with scanning transmission electron microscopy for Au nanoparticles embedded in an MgO support, deviations from an ideal topography during evaporation are demonstrated directly, and correlated with compositional errors in the reconstructed data. Finite element simulations of the field evaporation process confirm that protruding Au nanoparticles will evolve on the tip surface, and that evaporation field variations lead to an inaccurate assessment of the local composition, effectively lowering the spatial resolution of the final reconstructed dataset. Cross-correlating the experimental data with simulations results in a more detailed understanding of local evaporation aberrations during APT analysis of metal-oxide composites, paving the way towards a more accurate three-dimensional characterization of this technologically important class of materials.

Environmental Molecular Sciences Laboratory↗

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 OSTI2015

Hyperelastic tension of graphene

In this paper, we investigate the hyperelastic tensile behaviour of single layer graphene sheets (SLGSs). A one-term incompressible Ogden-type hyperelastic model is chosen to describe the mechanical response of C-C bonds. By establishing equality between the Ogden strain-energy and the variation of the Tersoff-Brenner interatomic potential, three different geometries of SLGSs are studied under tensile loading. We compute the Young's modulus, the finite-deformation Poisson's ratio, ultimate strains, total reactions, and the variation of the potential energy per carbon atom for large strains. Numerical simulations are compared with results obtained by molecular mechanics and molecular dynamics simulations, finite elements, continuum mechanics theory, and experiments. Our predictions are validated, revealing the potential predictive capabilities of the present hyperelastic framework for the analysis of graphene in the context of infinitesimal and large deformations. The good agreement found between our calculations and the published data suggests that graphene may be described as a hyperelastic material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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 OSTI2014

Ignition probability of polymer-bonded explosives accounting for multiple sources of material stochasticity

Accounting for the combined effect of multiple sources of stochasticity in material attributes, we develop an approach that computationally predicts the probability of ignition of polymer-bonded explosives (PBXs) under impact loading. The probabilistic nature of the specific ignition processes is assumed to arise from two sources of stochasticity. The first source involves random variations in material microstructural morphology; the second source involves random fluctuations in grain-binder interfacial bonding strength. The effect of the first source of stochasticity is analyzed with multiple sets of statistically similar microstructures and constant interfacial bonding strength. Subsequently, each of the microstructures in the multiple sets is assigned multiple instantiations of randomly varying grain-binder interfacial strengths to analyze the effect of the second source of stochasticity. Critical hotspot size-temperature states reaching the threshold for ignition are calculated through finite element simulations that explicitly account for microstructure and bulk and interfacial dissipation to quantify the time to criticality (t{sub c}) of individual samples, allowing the probability distribution of the time to criticality that results from each source of stochastic variation for a material to be analyzed. Two probability superposition models are considered to combine the effects of the multiple sources of stochasticity. The first is a parallel and series combination model, and the second is a nested probability function model. Results show that the nested Weibull distribution provides an accurate description of the combined ignition probability. The approach developed here represents a general framework for analyzing the stochasticity in the material behavior that arises out of multiple types of uncertainty associated with the structure, design, synthesis and processing of materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2019

Model for Characterization and Optimization of Spectrally Selective Structures to Reduce the Operating Temperature and Improve the Energy Yield of Photovoltaic Modules

Many existing commercially manufactured photovoltaic modules include a cover layer of glass, commonly coated with a single layer antireflection coating (ARC) to reduce reflection losses. As many common photovoltaic cells, including c-Si, CdTe, and CIGS, decrease in efficiency with increasing temperature, a more effective coating would increase reflection of sub-bandgap light while still acting as an antireflection coating for higher energy photons. The sub-bandgap reflection would reduce parasitic sub-bandgap absorption and therefore reduce operating temperature. This reduction under realistic outdoor conditions would lead to an increase in annual energy yield of a photovoltaic module beyond what is achieved by a single layer ARC. However, calculating the actual increase in energy yield provided by this approach is difficult without using time-consuming simulation. Here, we present a time-independent matrix model which can quickly determine the percentage change in annual energy yield of a module with a spectrally selective mirror by comparison to a baseline module with no mirror. The energy benefit is decomposed into a thermal component from temperature reduction and an optical component from increased transmission of light above the bandgap and therefore increased current generation. Time-independent matrix model calculations are based on real irradiance conditions that vary with geographic location and module tilt angle. The absolute predicted values of energy yield improvement from the model are within 0.1% of those obtained from combined ray-tracing and time-dependent finite-element simulations and compute 1000x faster. Uncertainty in the model result is primarily due to effects of wind speed on module temperature. Optimization of the model result produces a 13-layer and a 20-layer mirror, which increase annual module energy yield by up to 4.0% compared to a module without the mirror, varying depending on the module location and tilt angle. Finally, we analyze how spectrally selective mirrors affect the loss pathways of the photovoltaic module.

14 SOLAR ENERGY↗

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↗

DOE OSTI2026

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE↗