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.

SourceResearch connectionFrom yearThrough year

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

At least 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

NASA NTRS1999

A Hysteresis Model for Piezoceramic Materials

This paper addresses the modeling of nonlinear constitutive relations and hysteresis inherent to piezoceramic materials at moderate to high drive levels. Such models are, necessary to realize the, full potential of the materials in high performance control applications, and a necessary prerequisite is the development of techniques which permit control implementation. The approach employed here is based on the qualification of reversible and irreversible domain wall motion in response to applied electric fields. A comparison with experimental data illustrates that because the resulting ODE model is physics-based, it can be employed for both characterization and prediction of polarization levels throughout the range of actuator operation. Finally, the ODE formulation is amenable to inversion which facilitates the development of an inverse compensator for linear control design.

Smith, Ralph C.↗

NASA NTRS1985

Impact Crater Scaling Laws

The effects of the impactor size and velocity, of the material properties of the impactor and impacted body, and of the gravitational field strength were examined. The dependence of the coupling parameter on the Gruniesen parameter and the us/up slope values that characterize a Tillotson material model were tested. Whether a code calculation can recover the scaling exponents that can be derived theoretically for an idealized perfectly-porous material that crushes at zero strength and subsequently is incompressible was investigated. A model that was close to the idealized case was generated in a form suitable for code calculations and utilized for a series of one dimensional calculations. The code calculation demonstrated the theoretical coupling parameter. Calculations of impacts using a physically real porous material model seem to indicate a coupling parameter of the same form applicable to a nonporous material, distinctly different from that expected from other theoretical and experimental results for porous materials. The one dimensional calculations are being extended to two dimensional impacts.

Holsapple, K. A.↗

NASA NTRSDate not supplied

Guided Design of Alloys

Using multiscale materials modeling, we are designing stronger alloys for the airspace applications. We focus on improving mechanical properties of materials at operating temperature. We employ basic science and fundamental theory for guided design of the next-generation alloys. Among the examples, we consider precipitated superalloys with a local phase transformation strengthening. We acknowledge funding by NASA’s Aeronautics Research Mission Directorate (ARMD) via Transformational Tools and Technologies (TTT) Project.

multiscale↗

NASA NTRSDate not supplied

Guided Design of Alloys: Part 2

Multiscale↗

DOE OSTI2020

Tensor Basis Gaussian Process Models of Hyperelastic Materials

In this study, we develop Gaussian process regression (GPR) models of isotropic hyperelastic material behavior. First, we consider the direct approach of modeling the components of the Cauchy stress tensor as a function of the components of the Finger stretch tensor in a Gaussian process. We then consider an improvement on this approach that embeds rotational invariance of the stress-stretch constitutive relation in the GPR representation. This approach requires fewer training examples and achieves higher accuracy while maintaining invariance to rotations exactly. Finally, we consider an approach that recovers the strain-energy density function and derives the stress tensor from this potential. Although the error of this model for predicting the stress tensor is higher, the strain-energy density is recovered with high accuracy from limited training data. The approaches presented here are examples of physics-informed machine learning. They go beyond purely data-driven approaches by embedding the physical system constraints directly into the Gaussian process representation of materials models.

36 MATERIALS SCIENCE↗

NASA NTRS2019

A Multiscale Nonlocal Progressive Damage Model for Composite Materials

In this paper, the advantages of a nonlocal progressive damage formulation are described and demonstrated. An approximation of the nonlocal formulation was implemented coupled with the MAT162 composite damage model as a User defined material model in the LS DYNA environment. A comparison of the local model and the nonlocal model is simulated for an 8-ply laminate under tension is carried for increasing mesh densities. The results show the regularization achieved by nonlocal models by providing mesh independent results.

Kodagali, Karan↗

NASA NTRS2018

Development of a Mesoscale Finite Element Constitutive Model for Fiber Kinking

A mesoscale finite element material model is proposed to analyze structures that fail by the fiber kinking damage mode. To evaluate the assumptions of the mesoscale model, the results were compared with those of a high-fidelity micromechanical model. A direct comparison between the two models shows remarkable correlation, indicating that the key features of the fiber kinking phenomenon are appropriately accounted for in the mesoscale model. The mesoscale model is applied to structural analysis cases to demonstrate the capabilities of the model. A verification study is conducted with an unnotched compression specimen and preliminary validation is demonstrated with a notched compression specimen. The results show that the model is successful at representing the kinematics of fiber kinking while at the same time highlighting the need for further verification and validation.

Bergan, Andrew C.↗

NASA NTRSDate not supplied

Application of the MAT 213 Composite Impact Model to NASA Problems of Interest

As composite materials are gaining increased use in aircraft components where impact resistance under high-energy impact conditions is important (such as the turbine engine fan case), the need for accurate material models to simulate the deformation, damage, and failure response of polymer matrix composites under impact conditions is becoming more critical.

Polymer Matrix Composites↗

NASA NTRSDate not supplied

Progress on Quantitative Infrared Thermography at the NASA Langley Aerothermodynamic Laboratory

A quantitative infrared thermography technique is being developed for accurate and high-resolution measurements of surface temperature and heat flux in the hypersonic wind tunnels at NASA Langley Research Center. A hemisphere test campaign was carried out in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel to assess the current performance of the technique and to identify remaining challenges. Measurements were obtained for three different model materials over a range of freestream Reynolds numbers. The raw infrared images were converted to temperature via a radiometric calibration of the camera sensor. Corrections were made for losses due to test article emissivity and window transmissivity. A preliminary uncertainty analysis yielded an estimated accuracy within 4 K at moderate viewing angles, similar to that of a standard thermocouple. Surface heat flux was obtained from the measured surface temperatures via a finite-difference solution of the one-dimensional heat equation that modeled thermal variation of the material properties. Comparisons of the measured heating to both computational and theoretical predictions demonstrated agreement within the expected uncertainty of the technique. The heat transfer measurement near the stagnation point was found to be highly sensitive to the test article emissivity, which establishes the need for improved characterization of the wind tunnel model materials for high temperature applications.

infrared thermography↗

NASA NTRS1994

Comparison Of Models Of Metal-Matrix Composites

Report presents comparative review of four mathematical models of micromechanical behaviors of fiber/metal-matrix composite materials. Models differ in various details, all based on properties of fiber and matrix constituent materials, all involve square arrays of fibers continuous and parallel and all assume complete bonding between constituents. Computer programs implementing models used to predict properties and stress-vs.-strain behaviors of unidirectional- and cross-ply laminated composites made of boron fibers in aluminum matrices and silicon carbide fibers in titanium matrices. Stresses in fiber and matrix constituent materials also predicted.

Bigelow, C. A.↗

NASA NTRSDate not supplied

Computational Investigation of Oxidative Etch Pitting in FiberForm and Its Impact on Material Properties

Oxidation-driven carbon erosion does not occur uniformly but rather through the development of localized etch pits at active surface sites. These active sites form due to atomic defects on the carbon surface, making them significantly more reactive than the surrounding, non-defective areas. As a result, these sites are the first to react during ablation, leading to their removal. This process creates new defects in neighboring atoms, increasing their reactivity and causing localized carbon removal around these active sites. In this way, the highly reactive defective areas serve as nucleation points for the formation and growth of etch pits, which can have adverse effects on structural integrity of FiberForm. To better understand how these etch pits impact the material properties of carbon fiber microstructures, we have developed a new capability within the direct simulation Monte Carlo (DSMC) framework to capture the etch pit formation process. This capability, integrated into the DSMC code SPARTA (Stochastic Parallel Rarefied-gas Time-accurate Analyzer), models material removal in the presence of active sites, leading to the formation of etch pits. The current work focuses on studying the effects of these etch pits on the material properties of FiberForm, a widely used base material in thermal protection systems (TPS). The microstructure of virgin FiberForm, obtained via X-ray microtomography, is imported into SPARTA to generate the ablated geometries with etch pits. These modified microstructures are then analyzed using the Porous Microstructure Analysis (PuMA) software to compute various material properties, including elasticity, thermal conductivity, and permeability. We investigate the variation of these properties due to the complex surface topology changes caused by etch pit formation. Additionally, we compare the effects of pitting with the conventional model of shrinking fibers, traditionally used to simulate the ablation of carbon structures. Significant differences emerge between the two approaches. Consequently, this physically realistic model of material removal through etch pit formation offers improved accuracy in predicting the degradation of carbon-based TPS during oxidation. It also provides insights into other mechanisms, such as spallation, where chunks of material are removed into the flow due to etch pit growth. Ultimately, this model enhances our understanding of failure modes in these materials during ablation.

PuMA↗

DOE OSTI2018

American Society of Composites, 32nd Technical Conference

This paper will present a general methodology by which weave draping manufacturing simulation results can be utilized to include the effects of weave draping and scissor angle in a structural multiscale simulation. While the methodology developed is general in nature, this paper will specifically demonstrate the methodology applied to a truncated pyramid, utilizing manufacturing simulation weave draping results from ESI PAM-FORM, and multiscale simulation using Altair Multiscale Designer (MDS) and OptiStruct. From a multiscale simulation perspective, the weave draping manufacturing simulation results will be used to develop a series of woven unit cells which cover the range of weave scissor angles existing within the part. For each unit cell, a multiscale material model will be developed, and applied to the corresponding spatial locations within the structural simulation mesh. In addition, the principal material orientation will be mapped from the wave draping manufacturing simulation mesh to the structural simulation mesh using Altair HyperMesh mapping technology. Results of the coupled simulation will be compared and verified against experimental data of the same available via General Motors (GM) Department of Energy (DOE) project.

42 ENGINEERING↗

DOE OSTI2021

Efficient Failure Probability Calculations and Modeling Interface Debonding in TRISO Particles with BISON

This document constitutes completion of the NEAMS milestone, which is titled Demonstrate efficiency improvements in TRISO failure probability calculations and the effect of debonding on fission product transport in TRISO fuel particles in Bison . In this report we present the development of the following: (1) the anisotropic elasticity model of pyrolytic carbon and extended material models using the local coordinate system for aspherical particle geometry; (2) a capability of modeling interface debonding in TRISO particles; and (3) an efficient high-fidelity approach to calculate the failure probability of TRISO particles. Benchmark problems and AGR-2 tests are used to verify the models and demonstrate these new capabilities.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Shrinkage-induced deformations and creep of structural concrete: 1-year measurements and numerical prediction

Highlights: • Extensive experimental study on drying shrinkage, creep and microcracking of concrete • All specimens prepared from a single batch of ordinary-strength structural concrete • 1st year of measurements of (not only) non-uniformly drying beams with span up to 3-m • The database is downloadable from free-to-use research data repository • Modified MPS model for concrete creep used in blind prediction of all experiments The material models for creep and shrinkage operating on the material point level in FEM are usually intended for challenging complex applications and structures, where the average cross-sectional approach does not suffice. The identification of the growing number of material parameters induced by increasing model capabilities relies on very specific and narrow-oriented yet interconnected experiments which are scarce. The presented comprehensive experiments aim to provide a clearer image of the complicated interaction among the basic phenomena: drying, shrinkage, creep, and microcracking. The cornerstone of this ongoing research is a unique set of 30 partially-sealed unreinforced concrete beams with span 1.75–3.0 m subjected to drying. To minimize material variation, all specimens in this study were cast from a single batch of ordinary strength structural concrete with slag-blended binder. The resulting experimental database will be suitable both for validation and development of the constitutive models.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Digital Twin Model for Advanced Manufacture of a Rotating Detonation Engine Injector

A digital twin material model (DTMM) of an additive manufacturing (AM) process was created to advance the state of the art in rotating detonation engine (RDE) injector design. Current RDE injectors are designed with large pressure drops, enabling a stable and repeatable combustion process. However, this comes at the cost of system efficiency. For the technology to transition to commercial fossil-based power generation, it is important to develop injectors with reduced flow losses. Low-loss injectors are difficult to design and manufacture with conventional manufacturing techniques. AM enables new design options, but the AM manufacturing process must be thoroughly understood to result in a robust design. A DTMM provides the necessary insight by defining the cause-effect relationships between process parameters, microstructure features, and properties. Therefore, a DTMM to support the design and manufacturing process was developed and applied to the design of a new additively manufactured low-loss injector. The injector combustion behavior was characterized through hot-fire tests, and mechanical performance was compared to the DTMM predictions. The two project goals were the successful development of the DTMM and the demonstration of an improved RDE injector design. The RDE injector design and DTMM developments occurred on parallel but dependent paths. The injector was designed to reduce pressure drop by increasing the cross-sectional flow area ratio between the injector air passages and the combustor annulus. This resulted in less structural material, raising the concern that thin members would be susceptible to high-cycle fatigue (HCF) under the periodic loading inherent to an RDE. It was most important for the DTMM to predict behavior in these features; therefore, the injector design concept guided the material thicknesses used in fatigue tests. The DTMM development started by manufacturing a series of coupons over the range of possible AM process variations. A design-of-experiment approach was used to select which process variable combinations gave the most efficient coverage relevant to the injector design space. The microstructure in each of these coupons was characterized, and then computational methods were used to create a numerical model of the correlation between process variables and microstructure. Next, a set of HCF samples were tested to calibrate existing models that map microstructure to HCF performance. Together, these two links formed the DTMM that calculates HCF behavior from AM process variables. Two injector prototypes were additively manufactured. The first injector design strategy aggressively pursued low-loss performance by substantially increasing the oxidizer flow area. The combination of manufacturing lead times and the fatigue testing schedule meant that the DTMM was not available when building this first prototype. Therefore, its process parameters were chosen based on a manual review of the available coupon data. This prototype was built successfully and evaluated in 58 combustion tests. Sustained detonation was achieved with remarkably reduced pressure loss, and some tests even displayed pressure loss characteristics similar to conventional gas turbine combustors. This achieved the project goal of improving RDE injector design. The second injector was manufactured according to the optimized parameters predicted by the DTMM. The flow area modifications of this injector were less aggressive than the first injector since demonstrating low pressure loss was not an objective of the second hot-fire test series. Rather, the test objective was to cause high cycle fatigue failure in the part due to periodic loading from the rotating detonation wave. The observed number of cycles to failure was to be compared to the number predicted by the DTMM and thereby assess the utility of the DTMM in component design. However, the required level of vibration was not obtained during combustion. Therefore, high cycle fatigue was not experienced in the hot-fire tests of the second injector. Fatigue data was obtained by further testing the second injector in a conventional HCF test apparatus. The injector demonstrated HCF strength above the DTMM prediction. In fact, it did not fail and testing was only discontinued due to reaching the end of the period of performance. This points to some success in the project’s primary goal of successfully developing and applying the DTMM to a component design. Implementing the DTMM recommendations for optimal processing parameters led to a part with acceptable properties. The DTMM was also shown to be an efficient correlator of data and to provide insight into the relationship between process settings, microstructure, and property performance. However, the failure of the DTMM prediction to match the experimental result of the injector fatigue test also points to the need to include significantly more data in the model development. In this project, coupons made with identical processing parameters exhibited drastically different properties from each other and from the injector part, which clearly influences the accuracy of a model that predicts performance based on parameters. Uncertainties in the build process must be quantified to develop more robust models. A denser and broader matrix of coupon process and geometry variations, several repeated builds of every point, more in-situ build process measurements, and direct observation of tensile and HCF sample microstructure (as opposed to separate microstructure specimens) are recommendations to improve future AM modeling efforts.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2019

Application of Electronics Packaging Fundamentals to Photovoltaic Interconnects and Packaging

The energy capacity of photovoltaic (PV) installations worldwide has increased significantly in the last decade, increasing the demand for rigorous investigation of the physical phenomena causing degradation and failures in PV modules. As PV reliability science develops, established methods and approaches from longer-standing industries can inform and expedite PV reliability advances. This work demonstrates how thermomechanical solder bond fatigue models derived for electronics packaging applications can be applied to both standard and emerging PV interconnect designs. This expertise cannot be directly translated to PV however, as target PV module lifetimes are significantly longer than most electronics and the modules must withstand the natural climate wherever they are deployed. Verification of analytical reliability models is then an additional challenge, due to model size, test time, and climate variability. Furthermore, unconventional materials, such as low-temperature solders, are now being integrated into PV designs, for which appropriate material models must be selected and material parameters derived. In this paper, the current state of PV interconnect research is explored, with an emphasis on the experimental and simulation approaches and models being used for this work. Recent results for standard module architectures as well as emerging interconnect schemes are discussed.

14 SOLAR ENERGY↗

DOE OSTI2020

Inelastic equation of state for solids

In this study, a complete inelastic equation of state (IEOS) for solids is developed based on a superposition of thermodynamic energy potentials. The IEOS allows for a tensorial stress state by including an isochoric hyperelastic Helmholtz potential in addition to the zero-kelvin isotherm and lattice vibration energy contributions. Inelasticity is introduced through the nonlinear equations of finite strain plasticity which utilize the temperature dependent Johnson–Cook yield model. Material failure is incorporated into the model by a coupling of the damage history variable to the energy potentials. The numerical evaluation of the IEOS requires a nonlinear solution of stress, temperature and history variables associated with elastic trial states for stress and temperature. The model is implemented into the ALEGRA shock and multi-physics code and the applications presented include single element deformation paths, the Taylor anvil problem and an energetically driven thermo-mechanical problem.

42 ENGINEERING↗