Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material Modeling”

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 1,747 records · Page 97

Review of the synergies between computational modeling and experimental characterization of materials across length scales

With the increasing interplay between experimental and computational approaches at multiple length scales, new research directions are emerging in materials science and computational mechanics. Such cooperative interactions find many applications in the development, characterization and design of complex material systems. This manuscript provides a broad and comprehensive overview of recent trends in which predictive modeling capabilities are developed in conjunction with experiments and advanced characterization to gain a greater insight into structure–property relationships and study various physical phenomena and mechanisms. The focus of this review is on the intersections of multiscale materials experiments and modeling relevant to the materials mechanics community. After a general discussion on the perspective from various communities, the article focuses on the latest experimental and theoretical opportunities. Emphasis is given to the role of experiments in multiscale models, including insights into how computations can be used as discovery tools for materials engineering, rather than to “simply” support experimental work. This is illustrated by examples from several application areas on structural materials. In conclusion this manuscript ends with a discussion on some problems and open scientific questions that are being explored in order to advance this relatively new field of research.

3-D microstructure characterization↗

Multiscale Modeling of Carbon Fiber Reinforced Polymer (CFRP) for Integrated Computational Materials Engineering Process

In this work, a multiscale modeling framework for CFRP is introduced to study hierarchical structure of CFRP. Four distinct scales are defined: nanoscale, microscale, mesoscale, and macroscale. Information at lower scales can be passed to higher scale, which is beneficial for studying effect of constituents on macroscale part’s mechanical property. This bottom-up modeling approach enables better understanding of CFRP from finest details. Current study focuses on microscale and mesoscale. Representative volume element is used at microscale and mesoscale to model material’s properties. At microscale, unidirection CFRP (UD) RVE is used to study properties of UD. The UD RVE can be modeled with different volumetric fraction to encounter non-uniform fiber distribution in CFRP part. Such consideration is important in modeling uncertainties at microscale level. Currently, we identified volumetric fraction as the only uncertainty parameters in UD RVE. To measure effective material properties of UD RVE, periodic boundary conditions (PBC) are applied to UD RVE to ensure convergence of obtained properties. Properties of UD is directly used at mesoscale woven RVE modeling, where each yarn is assumed to have same properties as UD. Within woven RVE, there can be many potential uncertainties parameters to consider for a physical modeling of CFRP. Currently, we will consider fiber misalignment within yarn and angle between wrap and weft yarns. PBC is applied to woven RVE to calculate its effective material properties. The effect of uncertainties are investigated quantitatively by Gaussian process. Preliminary results of UD and Woven study are analyzed for efficacy of the RVE modeling. This work is considered as the foundation for future multiscale modeling framework development for ICME project.

Composite, ICME, Multiscale Modeling↗

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗

Kernel Enriched Meshfree Multiphysics Degradation Modeling of Energy Storage Materials

Energy storage materials exhibit strong electro-chemo-mechanical coupling and highly anisotropic material properties, contributing to the formation and propagation of micro-cracking during charge/discharge cycling and ultimately diminishing performance and service life. With microstructural images supplied by the National Laboratory of the Rockies (NLR), pixel-based meshfree model construction by the reproducing kernel particle method (RKPM) is used to represent the complex material microstructures that dictate the coupled physics of these systems. Traditional electro-chemo-mechanical models rely on mesh-based finite element methods, which can lead to difficulties in meshing such complex geometries and capturing crack propagation due to mesh dependency. The first kernel enrichment discussed will be the interface modified reproducing kernel (IM-RK) [1, 2], constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The IM-RK is especially useful for areas in which a known discontinuity-type is expected a priori. The second kernel enrichment to be discussed is a neural network-enhanced reproducing kernel (NN-RK) [3, 4], which is introduced to effectively model non-obvious damage and crack propagation in the material microstructures; the location, orientation, and solution transition near a localization are automatically captured by superimposed block-level NN optimizations. This NN enrichment approach allows for effective modeling of localizations via a fixed background discretization, relieving tedious efforts for adaptive refinement in traditional mesh-based methods. Applications to the heterogeneous microstructures of Li-ion battery cathodes will be presented to demonstrate the effectiveness of the proposed methods. NN-RK is additionally used to inform how crack opening and closure in turn affect the electro-chemo-mechanical responses in the material microstructure. References: [1] Wang, Y., Baek, J., Tang, Y. et al. "Support vector machine guided reproducing kernel particle method for image-based modeling of microstructures," Comput Mech 73, 907-942 (2024). https://doi.org/10.1007/s00466-023-02394-9. [2] Susuki, K., Allen, J. & Chen, J. S.. "Image-based modeling of coupled electro-chemo-mechanical behavior of Li-ion battery cathode using an interface-modified reproducing kernel particle method," Engineering with Computers (2024). https://doi.org/10.1007/s00366-024-02016-9. [3] Baek, J., Chen, J. S., Susuki, K., "Neural Network enhanced Reproducing Kernel Particle Method for Modeling Localizations," International Journal for Numerical Methods in Engineering, Vol. 123, 4422-4454 (2022). https://doi.org/10.1002/nme.7040.

97 MATHEMATICS AND COMPUTING↗

Assessing Alignment Between a Top-Down LCA Model and a Bottom-Up Supply Chain Model for Data Gap Identification and Prioritization

This work presents a method for evaluating the alignment between the results of bottom-up and top-down models and using the evaluation to enhance the bottom-up model and underlying data set. Disaggregation is used to enhance the embodied energy results of top-down models, which are then compared to bottom-up model results. Using the disaggregated top-down results makes it possible to identify specific data gaps and deficiencies within the bottom-up model and its data set. The top-down results will also provide information around which data gaps in the bottom-up model are contributing the most error to the results; in this way, data gaps can be identified as either critical or non-critical, and the effort needed to acquire and add data can be prioritized. Although the proposed method could be applied to any physical quantity or environmental impact that can be quantified with both bottom-up models (typically in physical units) and top-down models (typically in monetary units), embodied energy is the focus of this work. The disaggregated energy baseline is calculated from the throughput vector x, which gives the total economic activity required from each sector to provide a final demand. Using the direct requirements matrix, each element in x is disaggregated into a vector of inputs to the corresponding sector. Energy data from the interventions matrix B is then used to convert the monetary flows of the disaggregated x into physical energy flows. For each energy flow, the origin sector, the destination sector, and the energy carrier type will all be known. These energetic flows are analogous but not identical to direct energy inputs that would appear in a unit process used in life cycle analysis and other bottom-up methods, and thus can be used to evaluate the alignment between the top-down results and the results of a bottom-up model. To demonstrate and further develop the proposed method, embodied energy for chemicals assessed in the 2015 U.S. DOE bandwidth study on energy use within the chemical sector will be calculated using a top-down and a bottom-up model. The Materials Flows through Industry (MFI) modeling tool developed at the National Renewable Energy Laboratory is used as the bottom-up model, and the U.S. EEIO model, recently published by the U.S. Environmental Protection Agency, is used as the top-down model. U.S. EEIO and MFI results will be evaluated to determine how well the model results align. Where the two models are not in good alignment, the disaggregation method described above will be applied to the top-down model to generate a more detailed set of results that will aid in identifying data gaps in the MFI database. Using only MFI, data gaps can be identified in terms of which sectors have few or no unit processes in the database, but no information is available regarding the impact of the missing data on MFI results. The disaggregated top-down results will provide this information on which gaps in the MFI data are critical and have a significant impact on the results, and which are less important.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 Mark-I Station Blackout: Part II – Uncertainty Analysis

Single case comparisons between severe accident simulations can provide detailed insights into severe accident model behavior, however, they cannot offer insights into model uncertainty, sensitivity to uncertain parameters, or underlying model biases.Here in this analysis, the single case benchmark comparison of the MELCOR material interaction models for a station blackout (SBO) scenario of a boiling water reactor (BWR) using representative Fukushima Daiichi Unit 1 boundary conditions is expanded to include an uncertainty analysis. As part of this uncertainty analysis, 1200 simulations are performed for each material interaction model (2400 total), with random sampling of 14 uncertain MELCOR input parameters. Input parameters are selected for their impact on models representing core degradation processes. These include candling, fuel rod failure, debris quenching and dryout. The analysis performed here is not a traditional “best-estimate” uncertainty analysis that uses best-estimate parameters or identifies best-estimate figure of merit distributions. Instead, it is an exploratory uncertainty analysis that identifies and interrogates underlying model form biases of the two material interaction models (eutectics and interactive materials models). Uniform distributions are applied to all uncertain parameters to ensure coverage of the model parameter uncertainty space. Key findings from this study include underlying model form biases exhibited by material interaction models, and notable differences in accident progression outcomes between the material interaction models. This uncertainty study extends and confirms the conclusions from the first part of this study, which compared the impact of material interaction modeling on simulation of a short-term station blackout scenario with representative Fukushima Daiichi Unit I boundary conditions. In particular, this study confirms that the eutectics model generally exhibits accelerated degradation and failure of fuel components, the core plate, and the lower head. The eutectics model also has a tendency to exhibit a greater degree of core degradation, greater debris mass formation, and larger debris mass ejection. Finally, the eutectics model exhibits higher maximum temperatures for fuel, cladding, particulate debris, oxidic molten pool, and metallic molten pool components than the interactive materials model; interactive materials model simulations exhibit a soft “limitation” on maximum temperatures that is related to the temperature at which material relocation occurs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Integrated modeling of boron powder injection for real-time plasma-facing component conditioning

An integrated modeling framework for investigating the application of solid boron (B) powder injection for real-time surface conditioning of plasma-facing components (PFCs) in tokamak environments is presented. Utilizing the DIII-D impurity powder dropper (IPD) setup, this study simulates B powder injection scenarios ranging from milligrams to tens of milligrams per second, corresponding to boron flux rates of 10 20 –10 21 B/s in standard L-mode conditions. The comprehensive modeling approach combines EMC3-EIRENE for simulating the deuterium plasma background and the Dust Injection Simulator (DIS) for the ablation and transport of the boron powder particles. EMC3 trace impurity fluid modeling results show substantial boron transport to the inboard lower divertor, predominantly influenced by the main ion plasma flow. The dependency on powder particle size (5-250 µm) was found to be insignificant for the scenario considered. The effects of erosion and redeposition were considered to reconcile the discrepancies with experimental observations, which saw substantial deposition on the outer divertor plasma-facing components. For this purpose, the WallDYN3D code was updated to include boron sources within the plasma domain and integrated into the modeling framework. The mixed-material migration modeling shows evolving boron deposition patterns, suggesting the formation of mixed B-C layers or predominantly B coverage depending on the powder mass flow rate. While the modeling outcomes at lower B injection rates tend to align with DIII-D experimental observations, the prediction of near-pure boron layers at higher rates has yet to be experimentally verified in the carbon environment of the DIII-D tokamak. The extensive reach of boron layers found in the modeling suggests the need for modeling that encompasses the entire wall geometry for more accurate experimental correlations. This integrated approach sets a precedent for analyzing and applying real-time in-situ boron coating techniques in advanced tokamak scenarios, potentially extendable to ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Design and Analysis of a Mutual Inductance Level Sensor for Liquid Metals

Here, this article describes the design and analysis of an electromagnetic level sensor for use in high-temperature liquid metal systems. The mutual inductance level sensor (MILS) described in this work was fabricated using two single-conductor mineral insulated cables wrapped in a bifilar fashion around a stainless steel tube core and was housed in an isolating thimble that preserved the pressure boundary of the test vessel. Two sensor variations were fabricated that differ only in active length, 1016 and 1778 mm. Experimental data were collected using the 1016-mm sensor (MILS-MKII-040) that demonstrated a sensitivity of 9.2 μ V/mm in a room temperature testing stand that used solid aluminum as a surrogate for liquid metal. Experimental data were collected using the 1778-mm sensor (MILS-MKII-070) that demonstrated a sensitivity of 6.9 μ V/mm in the high-temperature (300 ° C) liquid sodium environment at the mechanisms engineering test loop (METL) of Argonne National Laboratory. The sensor performance was found to be repeatable over the course of several months, with roughly ±1% deviation from nominal output. Finite element models were developed in COMSOL Multiphysics that fully describe each test setup, and the models were validated using experimental data. The validated COMSOL models were used to perform an array of analyses that examined the performance of the sensor in differing environments. Maximizing the coil diameter inside the isolating thimble was found to maximize the signal and sensitivity of the sensor. An optimal operating frequency was found near 1000 Hz using both experimental data and COMSOL. The influence of a metallic thimble surrounding the sensor and a metallic sensor core was quantified and found to be negligible at the optimal operating frequency. The sensitivity of the sensor was quantified when monitoring the level of additional liquid metals. These include lead, lead-bismuth eutectic (LBE), sodium-potassium alloy (NaK), and lithium (in addition to sodium). The sensitivities were quantified using liquid metal properties at 350 ° C and 650 ° C. The geometry of the test stand model, all material properties used in the model, and the results are presented in a manner that allows the reader can replicate the model and perform additional analyses.

COMSOL↗

Low-Loss, High-Flux Saturation Magnetic Cores for Power Converters

As new and alternative energy sources are added to the energy grid, efficient and low-loss power converters are essential to the storage of energy from these intermittent power sources. Higher operating frequency can reduce the size and weight of transformers. However, current magnetic cores limit the operating frequency because core losses increase with frequency. Therefore, new magnetic core materials are needed with higher resistivity (lower eddy-current losses), lower coercivity (lower hysteresis losses), and higher permeability and flux density (higher power density). These new materials must be integrated into isolated power converters designed to optimize efficiency based on the transformer performance. To address this need, Mainstream pursued a two-pronged approach in Phase I to develop novel magnetic materials and model a planar transformer design. We irradiated silicon steel with a high-energy electron beam in an attempt to induce nanocrystalline grain formation. The nanocrystalline grains were expected to increase resistivity to reduce eddy-current losses and decrease coercivity to reduce hysteresis losses. We also modeled the transformer performance through Simulation Program with Integrated Circuit Emphasis (SPICE) and analytical models to establish the effect of magnetic properties of the core on the performance of the transformer. There still are hurdles to be overcome in Phase II. The electron beam processing was successful in increasing the resistivity and the flux saturation density. However, the resistivity is still too low for high-frequency magnetic cores, because the eddy-current losses would be prohibitively high. Additionally, the coercivity increased for the processed samples leading to higher hysteresis losses. In comparing our SPICE and analytical models to the planar transformer measured values and impedance plots, good agreement was seen for primary and secondary short-circuit impedance measurements and for primary-to-secondary capacitance measurements. The SPICE model deviated from measured values, though, when measuring the transformer open-circuit impedance, indicating that we might need to add more model complexity to the transformer model. Phase II research will continue to build on the transformer model and integrating novel magnetic materials into high-frequency transformers.

15 GEOTHERMAL ENERGY↗

Advanced Physics Models for Particle-Particle Interactions (UCD final report)

The LDRD "Advanced Physics Models for Particle-Particle Interactions” collaborative research with Professor Jennifer Sinclair Curtis and PhD student Chanh Nguyen at the University of California, Davis involved meso-scale, discrete element method (DEM) modeling of micron-sized particles impacting a substrate. Using DEM, physical insight into the details of particle response upon impact can be obtained. Such details provide the type of information necessary to enhance existing continuum or macro-scale simulation tools so that modeling for material strength and particle breakage models is based on particle physics versus empirical relationships. While molecular dynamics simulations of atoms can also provide physical insight into the particle response to impact, the scale of particles that are treated is on the order of angstroms, not microns. This study gives evidence to the ability of the DEM simulation technique to provide such fundamental, detailed insight for micron-sized particles as a function of particle properties and impact velocity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Additive Manufacturing Evaporative Casting

Traditional lost foam casting has been around for decades. The process uses foam forms blown by an injection molding like process or CNC milled into the desired shape, then pouring molten metal over the foam pattern to create a metallic object. Additive Manufacturing Evaporative Casting (AMEC) is a new process that eliminates foam forms by using additive manufacturing to 3D print the desired cast geometry. This not only saves time and money but allows for more advanced and complex designs that can’t be achieved by carving foam. Additionally, AMEC doesn’t require molds and tooling like other casting and foundry options. Because the AMEC process is new, extensive testing is needed to develop a better understanding of the process to minimize defects, quantify material properties, and start computer modeling for the process. This CRADA (collaborative research and development agreement) between ORNL and Skuld seeks to improve the process, develop a computational model, characterize material properties, and explore new applications.

36 MATERIALS SCIENCE↗

International SAMPE Symposium and Exhibition, 31st, Los Angeles, CA, April 7-10, 1986, Proceedings

The present conference on the development status of advanced structural materials considers topics arising in such areas as automated structural manufacturing, advanced material and structure design techniques, environmental effects on materials, composite matrix processing, computer modeling for materials and processes, materials development trends in Europe and in Japan, fiber and whisker reinforcement development status, and novel thermoplastic materials and their applications. Also discussed are pressure-sensitive adhesive systems, materials suitable for space applications, polyimide resin systems, electronic materials, novel resin chemistries, ceramic and metallic systems, and the impact performance of state-of-the-art materials.

Bauer, J. L.↗

Status of NASA Research on Projectile Shape Effects-Impact Simulations

Numerical evaluation of a Whipple shield has yielded a ballistic model for cylindrical projectiles. A numerical model of a Whipple shield covered with a thermal-blanket that is representative of operational shields has been developed using CTH: The shield model has been evaluated against existing test data; The projectile model has been developed from preexisting models for graphite-epoxy materials. The model has been used to identify projectile characteristics at the ballistic limit for spherical and cylindrical projectiles; Considered normal impacts of the shield at 7 kilometers per second; Varied the angle between the axis-of-symmetry and velocity vector; Varied the length to the diameter ratio over a broad range. The results have been consolidated into a generalized model that can be adapted to existing spherical models

Miller, J.↗

Rapid Tools for an AFP Manufacturing Defects Assessment Framework

This work formulates an automated fiber placement (AFP) defects assessment framework. Such framework assumes AFP manufacturing processes are able to identify manufacturing defects via automated inspection systems and intends to provide rapid analysis tools to create a defect assessment loop. The defect assessment loop defines an automated analysis process during AFP manufacturing to minimize the number of manual repairs on the part, thus accelerating AFP manufacturing of composite structures. This defects assessment loop advantages automated inspection data to evaluate the influence of encountered AFP defects on a ply-by-ply basis. Structural evaluation is based on strength criteria using local-global finite element models. Local models affect the global part model via material property reductions. Schemes to reduce material properties using the local models are the main enablers of this framework’s automated assessment capability. Finally, we discuss the main technical challenges to realize the feasibility of this framework.

Bahamonde Jacome, Luis G.↗

Microstructure Characterization and Modeling for Improved Electrode Design

This presentation describes research work led by NREL with team members from Argonne National Laboratory and Texas A&M University in microstructure analysis, modeling and validation under DOE's Computer-Aided Engineering of Batteries (CAEBAT) program. The goal of the project is to close the gaps between CAEBAT models and materials research by creating predictive models that can be used for electrode design.

33 ADVANCED PROPULSION SYSTEMS↗

A microstructure-guided constitutive modeling approach for random heterogeneous materials: Application to structural binders

This paper presents a microstructure-guided modeling approach to predict the effective elastic response of heterogeneous materials, and demonstrates its application toward two highly heterogeneous, uncon- ventional structural binders, i.e., iron carbonate and fly ash geopolymer. Microstructural information from synchrotron X-ray tomography (XRT) and intrinsic elastic properties of component solid phases from statistical nanoindentation are used as the primary inputs. The virtual periodic 3D microstructure reconstructed using XRT, along with periodic boundary conditions is used as a basis for strain- controlled numerical simulation scheme in the linear elastic range to predict the elastic modulus as well as the stresses in the microstructural phases. The elastic modulus of the composite material predicted from the microstructure-based constitutive modeling approach correlates very well with experimental measurements for both the materials considered. This technique efficiently links the microstructure to mechanical properties of interest and helps develop material design guidelines for novel heterogeneous composites

36 MATERIALS SCIENCE↗

Verification of a Finite Element Model for Pyrolyzing Ablative Materials

Ablating thermal protection system (TPS) materials have been used in many reentering spacecraft and in other applications such as rocket nozzle linings, fire protection materials, and as countermeasures for directed energy weapons. The introduction of the finite element model to the analysis of ablation has arguably resulted in improved computational capabilities due the flexibility and extended applicability of the method, especially to complex geometries. Commercial finite element codes often provide enhanced capability compared to custom, specially written programs based on versatility, usability, pre- and post-processing, grid generation, total life-cycle costs, and speed.

ablating materials↗