Search NASA⌕ Search

SEARCH · Search NASA

Results for “microstructure scale 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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Neptune inbound bow shock

Voyager 2 crossed the inbound or upstream Neptunian bow shock at 1430 spacecraft event time on August 24 in 1989 (Belcher et al., 1989). The plasma and magnetic field measurements allow us to study the solar wind interaction with the outermost gas giant. To fully utilize all of the spacecraft observations, an improved nonlinear least squares, 'Rankine-Hugoniot' magnetohydrodynamic shock-fitting technique has been developed (Szabo, 1994). This technique is applied to the Neptunian data set. We find that the upstream bow shock normal points nearly exactly toward the Sun consistent with any reasonable large-scale model of the bow shock for a near subsolar crossing. The shock was moving outward with a speed of 14 +/- 12 km/s. The shock can be characterized as a low beta, high Mach number, strong quasi-perpendicular shock. Finally, the shock microstructure features are resolved and found to scale well with theoretical expectations.

Szabo, Adam↗
DOE OSTI2023

Effect of sample size on the maximum value distribution of fatigue driving forces in metals and alloys

An analytical framework is presented to predict the effects of sample size on the maximum value distribution (MVD) of the driving forces for fatigue crack formation in metals and alloys. The distribution of the maximum driving force for fatigue crack formation over the domain follows the generalized extreme value theory in the limit as the domain size increases to infinity. Here, a simulation-based analysis of microstructure influences on fatigue resistance for polycrystalline metals and alloys is very costly, and reaching those limits is intractable. This work models the MVD of Fatigue Indicator Parameters (FIPs), which serve as surrogate measures for the driving force for fatigue crack formation, at finite sample sizes prior to their convergence to a limiting extreme value distribution. Large-scale crystal plasticity finite element (CPFE) simulations of FCC Al 7075-T6 with microstructure realizations of various sizes are incorporated to calibrate and evaluate the developed framework, and a total of ∼6.5 million grains of Al 7075-T6 are examined. The calibrated analytical solution agrees well with the brute force Monte Carlo simulation results extracted from the CPFE simulations. Furthermore, the developed formulation can predict the MVD of FIPs for different sample sizes using a size-dependent parameter, and it is capable of accurately extrapolating the MVD of FIPs for much larger microstructure sample sizes than the size used for its calibration.

Crystal plasticity↗
DOE OSTI2025

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗
DOE OSTI2021

Multiscale thermal properties prediction in the Multiphysics Object Oriented Simulation Environment (MOOSE) via a general Boltzmann solver [Poster]

The project objectives were: (1) Engineering-scale fuel performance modeling relies on accurate thermal properties; (2) Thermal properties are inherently multi-scale, arising from atomistic processes and interactions with a material’s microstructure; (3) Heat transport in solids via conduction occurs through transport and scattering of electrons and phonons; (4) Use the Boltzmann transport equation (BTE) to predict the macroscopic behavior of a materials system in terms of the microscopic dynamics of its heat carriers; and (5) This project establishes a new MOOSE (Multiphysics Object Oriented Simulation Environment) module, Boltzmann, dedicated to phonon and thermal electron transport.

36 MATERIALS SCIENCE↗
NASA NTRS2015

Multiscale Modeling of Ceramic Matrix Composites

Results of multiscale modeling simulations of the nonlinear response of SiC/SiC ceramic matrix composites are reported, wherein the microstructure of the ceramic matrix is captured. This micro scale architecture, which contains free Si material as well as the SiC ceramic, is responsible for residual stresses that play an important role in the subsequent thermo-mechanical behavior of the SiC/SiC composite. Using the novel Multiscale Generalized Method of Cells recursive micromechanics theory, the microstructure of the matrix, as well as the microstructure of the composite (fiber and matrix) can be captured.

ceramic matrix composites↗
NASA NTRS2019

Pot-ability Assessment of Litz Wires for High Power Density Electric Motor

An alternative process technique, namely vacuum-assisted axial injection potting (VaAIP), has been developed to pot the Litz wires in the stator winding of high power density electric motors for the future electrified aircrafts. Initial trials of the process showed significant improvement in potting quality with less voids, thus potential improvement in thermal management of the motors. As an initial effort of pot-ability assessment, microstructures, 2-D and 3-D, of the Litz wires including dimensions and distribution of conductor filament, coating, and open spaces; packing patterns; shape/configuration changes of each bundles or the overall cross-sections per degree of twist were determined and quantified successfully. The microstructure analyses were performed not only for effective potting process development but also for more realistic electro-thermal modeling solutions. This paper will present results of the microstructure analyses, potentials of the VaAIP process from the trials, and future plans for scale-up and implementation of the process into a full-scale prototype stator winding.

Electrified aircraft↗
DOE OSTI2022

Time-dependent THMC properties and microstructural evolution of damaged rocks in excavation damage zone

Modeling coupled thermo-hydro-mechanical-chemical (THMC) processes in host rocks near high-level nuclear waste (HLW) repositories at various time scales is an extremely challenging task. The current study integrates experimental, theoretical, and numerical methods in assessing the evolution of excavation damage zone (EDZ) over time and its implication on the long-term migration of hazardous species. Argillite and rock salt and are the focus of this study. The first part of the report presents a novel time-dependent directional microcrack damage theory for generic brittle rocks. It features detailed statistical description of the microcracks within a damaged solid, permitting a direct upscaling of microscale processes such as crack growth kinetics, crack closure/opening, sliding friction to explain the macroscopic creep, nonlinear elasticity, shear dilation, and loading-unloading hysteresis. This provides a basic platform for describing the anisotropic mechanical and transport properties of damaged rocks during excavation and subsequent THMC loadings. The model is validated and numerically implemented to Finite Element (FE) package ABAQUS through the user-defined material (UMAT) interface and have demonstrated great potential in resolving the time-dependent and anisotropic evolution of damage in the EDZ. The second part of the report presents a multi-scale experimental effort in characterizing the thermal, hydraulic, and mechanical properties of Mancos shale and Avery Island salt. For the Mancos shale, triaxial compression tests are performed at different confining pressures and temperatures to probe its thermomechanical properties relevant to HLW repositories. The obtained stress-strain data are interpreted using the proposed directional damage theory. Post-test specimens are subjected to gas permeability tests to reveal the correlation between permeability and the degree of microcracking. At microscale, temperature-controlled nanoindentation tests are performed and found a linear correlation between fracture toughness and elastic modulus from 25°C to 300°C. For Avery Island salt, we have designed and manufactured a novel relative-humidity controlled uniaxial creep device. Long-term creep tests at low stresses (< 5 MPa) are performed at different levels of relative humidity (RH). Besides confirming the much higher creep rates as one would expect through extrapolating the high-stress creep data, the results reveal that the steady-state creep rate of rock salt is strongly dependent on the ambient RH, an aspect that is often neglected in the literature. Both behaviors can be attributed to the pressure-solution creep mechanism which dominates at low stress and high RH levels. The third part of the report explores a set of numerical strategies in modeling the THMC behavior of porous geomaterials. A fully implicit, monolithic FE solution that can flexibly interface with different material models and coupling mechanisms for THMC problems is developed and verified through the ABAQUS user-defined element (UEL) interface. The scheme is used to study the THM response of a hypothetical HLW storage site with reference to an existing in-situ heater test. Strategies for integrating the UEL and the microcrack UMAT are suggested. Another numerical endeavor of this study is to implement a higher-order asymptotic homogenization method to account for the heterogeneous porous structures. The same method is then extended to perform microstructure-informed thermo-mechanical modeling of generalized continua. The above outcomes of this project provide a strong thrust towards enhancing the fundamental understanding and modeling capability of the long-term evolution of host rocks in EDZ, thus helping achieve the design goal of 1-million-year isolation of high-level nuclear wastes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Thixotropic spectra and Ashby-style charts for thixotropy

There is no universal model for thixotropy, and comparing thixotropic effects between different fluids is a subtle yet challenging problem. We introduce a generalized (model-insensitive) framework for comparing thixotropic properties based on thixotropic spectra. A superposition of exponential stress modes distributed over thixotropic time scales is used to quantify buildup and breakdown times and mode strengths in response to step-change input. This mathematical framework is tested with several experimental step-shear rate data on colloidal suspensions. Low-dimensional metrics based on moments of the distribution reveal characteristic average thixotropic properties, which are visualized on Ashby-style diagrams. This method outlines a framework for describing thixotropy across a diverse range of microstructures, supporting scientific studies as well as material selection for engineering design applications.

Mechanics↗
NASA NTRS1995

3D thin film microstructures for space microrobots

Micromechanisms of locomotion and a manipulator with an external skeleton like the structure of an insect are proposed. Several micro-sized models were built on silicon wafers by using polysilicon for rigid plates and polyimide for elastic joints. Due to scale effects, friction in micromechanical components is dominant as compared to the inertial forces because friction is proportional to L(exp 2) while mass is proportional to L(exp 3). Therefore, to ensure efficient motion, rotational joint that exhibits rubbing should be avoided. In this paper, paper models of a robot leg and a micro-manipulator are presented to show structures with external skeletons and elastic joints. Then the large scale implementation using plastic plates, springs, and solenoids is demonstrated. Since the assembly technique is based on paper folding, it is compatible with thin film micro-fabrication and integrated circuit (IC) planar processes. Finally, several micromechanisms were fabricated on silicon wafers to demonstrate the feasibility of building a 3D microstructure from a single planar structure that can be used for space microrobots.

Shimoyama, Isao↗
DOE OSTI2024

A non-isothermal breakage-damage model for plastic-bonded granular materials incorporating temperature, pressure, and rate dependencies

Plastic-bonded granular materials (PBM) are widely used in industrial sectors, including building construction, abrasive applications, and defense applications such as plastic-bonded explosives. The mechanical behavior of PBM is highly nonlinear, irreversible, rate dependent, and temperature sensitive governed by various micromechanical attributions such as grain crushing and binder damage. This paper presents a thermodynamically consistent, microstructure-informed constitutive model to capture these characteristic behaviors of PBM. Key features of the model include a breakage internal variable to upscale the grain-scale information to the continuum level and to predict grain size evolution under mechanical loading. In addition, a damage internal state variable is introduced to account for the damage, deterioration, and debonding of the binder matrix upon loading. Temperature is taken as a fundamental external state variable to handle non-isothermal loading paths. The proposed model is able to capture with good accuracy several important aspects of the mechanical properties of PBM, such as pressure-dependent elasticity, pressure-dependent yield strength, brittle-to-ductile transition, temperature dependency, and rate dependency in the post-yielding regime. Furthermore, the model is validated against multiple published datasets obtained from confined and unconfined compression tests, covering various PBM compositions, confining pressures, temperatures, and strain rates.

Breakage↗
DOE OSTI2024

Stress localization investigation of additively manufactured GRCop-42 thin-wall structure

A full-field crystal plasticity (CP) framework is presented for the GRCop-42 alloy to study microscopic mechanical behavior and local stress heterogeneities. The microstructures of additively manufactured (AM) materials are often unique relative to conventionally processed materials, and the local thermal histories drive these differences during the build process. These thermal histories depend on the process parameters (laser power, scan speed, and scan strategy) and the part geometry. Prior research has shown that the mechanical properties of thin-walled structures can vary significantly with wall thickness due to changes in the thermal boundary conditions during manufacturing. It is, therefore, desirable to perform CP simulations based on the phenomenological constitutive model to predict the local mechanical responses induced by microstructural heterogeneities. This work generates representative microstructures based on experimentally collected grain information (i.e., texture) for grain scale stress analysis, and the material constitutive parameters are calibrated using the experimental mechanical testing data. Here, we specifically investigated the effect of crystallographic texture and grain morphologies on the size-dependent mechanical properties of AM GRCop-42. The selection of appropriate material properties for implementing an effective free surface boundary condition and the influence of adjacent buffer layers are also discussed. Analysis of local field results reveals a strong correlation between stress localization and the initial grain orientation. However, no significant relationship between the misorientation of the individual adjacent grains and the average misorientation is observed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Characterization of Radiation Induced Defects Across Scales

The Whole Powder Pattern Modelling (WPPM) approach of diffraction line profile analysis (LPA) was applied for the first time on PuO 2 in order to understand the evolution of the microstructural properties as a function of process conditions. Crystallite size distributions and defect concentrations were determined by WPPM on three different PuO 2 samples - including one with a bimodal distribution and compared to the result of more common LPA analyses used in PuO 2 literature studies (e.g. Scherrer equation, double Voigt method). The conclusion was that a lognormal distribution of spherical crystallites with dislocations most closely represent the microstructure of PuO 2 produced by heating Pu 2 (C 2 O 4 ) 3 (H 2 O) 6 •3H 2 O between 650°C and 850°C. The volume-weighted mean column length (Lvol) parameter obtained from the double Voigt method for the monomodal samples were the same as the Lvol obtained from the WPPM within uncertainty. The two LPA methods did not agree on the Lvol parameters of the bimodal distribution sample with only the WPPM getting values that were consistent with the individual components. The WPPM can potentially provide more insight into future studies of the influences on the microstructure properties of PuO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Image-Based Digital Twin for Assessing the Coupled Electro-Chemo-Mechanical Behavior of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is 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 discrete nature of the images' pixel points is employed throughout the model and approximation construction. The mechanical model is verified using an image-based microstructure under tensile loading. A transient electro-chemo-mechanical coupled simulation is performed to evaluate the potential micro-cracking induced degradation of the battery cathode material subjected to charge/discharge cycling.

battery degradation↗
DOE OSTI2023

High fidelity model of directed energy deposition: Laser-powder-melt pool interaction and effect of laser beam profile on solidification microstructure

Metal additive manufacturing technologies keep receiving a great deal of interest as well as strong requests to develop methods to link the process science to printed parts performance and understand how to overcome inherent limitations. Here, a high-fidelity model based on the multiphysics ALE3D code was developed to reproduce the directed energy deposition process down to the powder scale. This includes resolving the laser-powder-melt pool interactions (powder impingement and incorporation into melt pool, hydrodynamics flow condition and laser absorption inefficiencies) as well as the resulting solidification microstructure. This micrometer scale digital twin captured the effect of powder incorporation process and powder flow rate on porosity. Furthermore, it was used to explore how a ring laser beam profile instead of the standard Gaussian laser profile could decrease the thermal gradient along the solidification front in the melt pool, which in turn can increase propensity for more desirable equiaxed grains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phase Stability of Yttrium Titanium Oxide

This proposed study aims at elevating Idaho National Laboratory (INL)’s portfolio in nuclear material research on two major fronts: 1) expanding INL’s material characterization capabilities to provide holistic microstructural and microchemical characterization ranging from atomic scale to millimeter scale, and 2) filling the knowledge gap on the effect of microstructural evolution on material phase stability at extreme conditions resembling those inside nuclear reactors, using Yttrium-Titanium-Oxide (YTO) as a model material. Previous research indicated the kinetics of phase transition of YTO materials varies between bulk samples and powder samples, suggesting that microstructural characteristics such as defect, micro-strain, and crystallite size could play a significant role in phase transition. The effect of crystallite size on the thermodynamics of phase transition has been well documented for other ceramic materials at nanometer crystallite size, with potential impacts from surface and strain energy and defects. In this proposed study, we will perform microstructural (micro-strain, crystallite size, defects, phase) and microchemical (elemental mapping) analyses through in-situ high temperature measurements using X-ray diffraction (XRD; at millimeter scale for powder and bulk samples) and transmission electron microscopy (TEM, at nanometer scale and atomic scale for powder and thin foil samples), to provide a comprehensive view of the effect of microstructural and microchemical evolution on phase transition and phase stability of YTO materials. These findings will also establish a roadmap and provide new insights on the phase stability of other materials that are critical to nuclear energy and other renewable energy infrastructure. Besides advancing our understanding of this fundamental material science question, this project will also enable the development of state-of-the-art material characterization capabilities at INL/MFC for microstructural and microchemical analysis at high temperature conditions. The new collaboration between the High-Resolution Materials Characterization group at INL and the Electron Microscopy group at the University of Florida (UF) in this project will open the possibilities for more fruitful and synergetic cutting-edge nuclear material research in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Diffusion, atomic transport, and ordering in Al-Zr alloys: FCC and liquid phases

Additive manufacturing of materials with controlled microstructure demands knowledge of atomic scale properties near the solid-liquid transition state. Many of these properties are not affordable by experimental techniques and computer modeling is the possible solution to the problem. In this paper, we present the results of an extended atomistic study of intrinsic atomic transport due to vacancy diffusion in FCC and L12 solid phases and diffusion in the liquid phase of Al-Zr alloys. A deceleration of the overall self-diffusion was observed when Zr was added to Al. The effect was stronger in the solid and weaker in the liquid. Additionally, the effect was strongly temperature dependent in the solid phases, but not in the liquid. Atomic transport was chemically biased: transport of Zr atoms was significantly slower than that of Al atoms, and this bias effect was stronger in the solid phases. The overall diffusion and chemical ordering processes in the liquid state were five to six orders in magnitude faster than in the solid. Chemical short-range order parameters in the liquid saturated at values close to those in the ordered L12 structure of Al3Zr. Chemical and structural ordering in the solid phases was negligible over the modeled microsecond time scale. Here, the results are discussed in view of optimizing additive manufacturing parameters for the controlled formation of metastable L1 2 precipitates.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A comparative study of solute trapping in Fe-(33–45 at%) Cu alloys manufactured by laser directed energy deposition

The high cooling rates in laser directed energy deposition (DED-LB) of alloys lead to substantial amounts of solute trapping as solute atoms cannot diffuse away from the solid/liquid interface before it advances. In some concentrated alloys, this results in supersaturated solid phases that form nanoscale hierarchical microstructures when the solute atoms precipitate out during reheating from subsequent laser passes. We choose the Iron-Copper (Fe-Cu) binary alloy as model system as it is chemically homogeneous in the liquid phase and has negligible solid solubility at room temperature. Two alloys with nominal compositions in atomic (at.) %, Fe67Cu33 and Fe55Cu45, were manufactured using DED-LB. Scanning transmission electron microscopy (STEM), energy dispersive spectroscopy (EDS) and wavelength dispersive spectroscopy (WDS) were used to characterize the nanostructures and heterogeneous chemical compositions. A non-equilibrium solute partitioning model was used to compute the supersaturated chemistries of the constituent phases and validated with experimentally measured compositions. The measured phase compositions of the two alloys were very similar, at roughly 12 at.% Cu and 4 at.% Fe in the α(bcc)-Fe and ε(fcc)-Cu phases respectively, despite having different processing parameters and mechanical behavior. This indicates that the total thermal history, that depends on both the processing parameters and laser scan pattern, plays a stronger role on the final microstructure evolution than just the initial quantity of trapped solute. In addition, we find that current non-equilibrium solute partitioning models applied on the continuum scale fall short of predicting accurate quantitative phase compositions in concentrated alloys, although the qualitative trends are captured correctly.

36 MATERIALS SCIENCE↗
DOE OSTI2022

304L Can Crush Validation Studies

Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.

36 MATERIALS SCIENCE↗