Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuous material properties”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for inverse problems: Determining material properties for heterogeneous systems

Here, we determine spatially varying discontinuous material properties using a domain-decomposition based physics-informed neural networks (PINNs) framework named the Adaptive Interface-PINNs or AdaI-PINNs (Roy et al., 2024). We propose the use of distinct neural networks for the field variables and material properties within each material, utilizing adaptive activation functions. While the neural networks across different materials share the same weights and biases, their activation functions are uniquely tailored using a hyperparameter that influences the slope of the activation function. The proposed framework is tested on several one-dimensional and two-dimensional benchmark examples, and its performance is compared with conventional PINNs and existing domain-decomposition PINNs frameworks, namely, the Multi-domain physics-informed neural network (M-PINN), and the eXtended physics-informed neural networks (XPINNs). The results demonstrate that the proposed approach can determine randomly distributed discontinuous material properties with an L 2 error of $\mathscr{O}$ (10 -3 ) for the material property and the root-mean-square error of $\mathscr{O}$ (10 -3 ) for the primary variable while the other approaches yield errors that are approximately two orders of magnitude larger (that is, $\mathscr{O}$ (10 -1 )). Moreover, the spatial distribution of material properties obtained using the proposed framework is in close agreement with the true distribution, whereas the other approaches fare much worse. Additionally, the proposed approach is approximately 40% faster than its competitors, indicating its potential as a robust alternative for solving inverse problems in heterogeneous materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Physics-informed neural networks for heterogeneous poroelastic media

This study presents a novel physics-informed neural network (PINN) framework for modeling poroelasticity in heterogeneous media with material interfaces. The approach introduces a composite neural network (CoNN) where separate neural networks predict displacement and pressure variables for each material. While sharing identical activation functions, these networks are independently trained for all other parameters. To address challenges posed by heterogeneous material interfaces, the CoNN is integrated with the Interface-PINNs (I-PINNs) framework (Sarma et al., Comput. Methods Appl. Mech. Eng. 429: 117135, 2024), allowing different activation functions across material interfaces. Further, this ensures accurate approximation of discontinuous solution fields and gradients. Performance and accuracy of this combined architecture were evaluated against the conventional PINNs approach, a single neural network (SNN) architecture, and the eXtended PINNs (XPINNs) framework through two one-dimensional benchmark examples with discontinuous material properties. The results show that the proposed CoNN with I-PINNs architecture achieves an RMSE that is two orders of magnitude better than the conventional PINNs approach and is at least 40 times faster than the SNN framework. Compared to XPINNs, the proposed method achieves an RMSE at least one order of magnitude better and is 40% faster.

42 ENGINEERING↗
DOE OSTI2025

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Modeling diffusion and types I-V sorption of water vapor in heterogeneous systems

Here, a sorption and diffusion model is developed using a mobile-immobile decomposition of the vapor concentration. Henry’s mode represents the mobile species which diffuses through the solid, while Langmuir and pooling modes represent possible immobilization processes that cause nonlinearities in sorption capacity. The model can simulate the five classical sorption types and can incorporate additional dynamics compared to models based on thermodynamic equilibrium. Additionally, the framework allows for discontinuities in material properties encountered in multi-material systems. The resulting non-linear coupled equations are solved by employing a finite element method to discretize in space, and a backward Euler method to discretize in time. The discretized system of equations is solved via a Newton-type iteration scheme at each time-step. Four different materials are parameterized for the model, which are then used in examples to demonstrate the model’s ability to capture complex sorption processes in multi-material systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Elastic Wave Propagation in Curvilinear Coordinates with Mesh Refinement Interfaces by a Fourth Order Finite Difference Method

In this work, we develop a fourth order accurate finite difference method for the three dimensional elastic wave equation in isotropic media with the piecewise smooth material property. In our model, the material property can be discontinuous at curved interfaces. The governing equations are discretized in second order form on curvilinear meshes by using a fourth order finite difference operator satisfying a summation-by-parts property. The method is energy stable and high order accurate. The highlight is that mesh sizes can be chosen according to the velocity structure of the material so that computational efficiency is improved. At the mesh refinement interfaces with hanging nodes, physical interface conditions are imposed by using ghost points and interpolation. With a fourth order predictor-corrector time integrator, the fully discrete scheme is energy conserving. Numerical experiments are presented to verify the fourth order convergence rate and the energy conserving property.

58 GEOSCIENCES↗
DOE OSTI2021

Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times

Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.

42 ENGINEERING↗
DOE OSTI2025

An innovative radial gradient material design using hot isostatic pressing for applications in extreme environments

Functionally graded materials (FGMs) are highly advanced continuous or discontinuous structures whose structural and material properties vary along a singular geometric dimension either in the axial or radial direction. Here, the radial gradient FGM design makes for an optimal structural design to incorporate a bi-metallic structure with a copper-based high entropy alloy (Cu-HEA) with good mechanical properties and high irradiation resistance, and Chromium (Cr) with great corrosion resistance. This study focuses on the experimental design of a metal powder loading mechanism to fabricate a bi-metallic radial gradient structure using Cu-HEA and 99.9 % pure Cr metal powders. The powder loading strategy uses custom-designed concentric cylindrical dividers to separate the individual compositions. Two benchtop trial runs were performed for design optimization. The optimized design was then implemented to eventually load the HEA and Cr powders for consolidation via powder metallurgy hot isostatic pressing (PM-HIP). The electron microscopy analysis reveals the successful fabrication of the radial gradient structure with the chemical mapping analysis, demonstrating the gradual composition shift from the HEA at the center to the pure-Cr at the periphery via a three-step gradient.

High Entropy Alloys (HEAs)↗
DOE OSTI2024

Reinforcement induced microcracking during the conversion of polymer-derived ceramics

This study investigates the role of discontinuous reinforcement on the microcracking of a preceramic polymer matrix during polymer to ceramic conversion. The material system comprised a carbosilane-based, preceramic polymer reinforced with mullite particles. The preceramic resin slurry was formulated for photopolymerization on a digital light projection 3D printer. Micro X-ray computed tomography, using a synchrotron light source, monitored a printed composite part during pyrolysis. The conversion of the carbosilane matrix into Si(O)C generated microcracks in the matrix that radially extended from particles. The likelihood of microcrack formation positively correlated with particle size. The largest particles (>10 4 µm 3 volume or >60 µm side length) always abutted microcracks, while the matrix surrounding smaller particles (<5 µm) was free of microcracks. Numerical calculations showed the particles resist the shrinking matrix during its conversion. This created tensile stresses within the matrix. The driving force for microcrack growth was also analyzed with regard to the influence of matrix material parameters, particle form-factor and crack length. These results reveal the need to consider the form factor of discontinuous reinforcements, with regard to the material properties of the polymer-derived ceramic, in order to reduce or eliminate defects in the final part.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Shock Stabilization of the HLLC Riemann Solver for the Carbuncle Instability

The HLLC approximate Riemann solver improves upon the HLL Riemann solver by resolving contact discontinuities. This is a particularly desirable property for multi-material codes in which problems usually contain material interfaces. However, the HLLC solver is known to suffer from the carbuncle phenomenon, a numerical instability most apparent at grid-aligned shocks in multi-dimensional simulations. Many problems of interest, including high energy-density physics applications, require the accurate resolution of both material interfaces and hydrodynamic shocks. A variety of methods have been developed to cure this instability, with varying degrees of complexity. The objective of this work is to describe a simple approach to modify the HLLC Riemann solver and prevent the carbuncle instability. Here, the method is then demonstrated for assorted two-dimensional test problems known to exhibit the shock instability. The performance of the new solver is compared with that of the standard HLL and HLLC Riemann solvers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Interface-aware Methods for Diffusion Equations [Slides]

In this talk, we will present interface-aware methods for diffusion equation. The focus is on the necessity of interface-aware, mixed cells, interface reconstruction, discontinuity of material diffusion coefficients, order of accuracy in time, property for large time steps, and simpler interface-aware methods. Numerical examples are provided to demonstrate the features on these aspects of the methods

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Segmentation method comparison for residual fiber length measurement across tiled microscopy images

Fiber length distribution (FLD), in part, governs mechanical properties in discontinuous fiber composites, yet manual measurement methods limit the high-throughput characterization needed for materials design optimization. This study compares deep learning segmentation approaches for automated FLD measurement in large-field microscopy, evaluating how method choice affects the microstructural descriptors used in structure-property-processing relationships. A critical challenge is that high-resolution microscopy images (10,000×10,000 pixels) must be tiled for deep learning analysis, fragmenting fibers at boundaries. We demonstrate that segmentation method proves crucial for measurement accuracy. For example, instance segmentation with Slicing Aided Hyper Inference (SAHI) preserves individual fiber integrity across tiles while semantic segmentation prioritizes speed. Comparing against manual measurement of extracted carbon fibers, YOLOv11-SAHI matched manual ground truth (238 μm weighted mean) with 40x speedup (4.5 vs 167 minutes per image). U-Net provides rapid quantification although it is at the cost of reduced accuracy due only reliably measuring stand-alone fibers. Our comparative analysis reveals that instance segmentation with SAHI better preserves length measurements while semantic segmentation prioritizes speed, providing empirical guidance for method selection. The characterization provides essential inputs for mechanical property prediction models and inverse design workflows, accelerating composite materials development cycles.

Additive manufacturing↗
DOE OSTI2022

Discontinuous Aligned Carbon Fiber Intermediates for Automotive and Related Applications

This work focused on preferentially aligning discontinuous carbon fibers in wet-laid or air-laid processes. It is well known that aligned fibers provides higher directional strength and stiffness. Discontinuous fibers further allow higher degree of draw and formability as the gaps in the fibers allow for higher material movement. The current processes are limited in their ability to align carbon fibers during processing. The aligned fibers have several benefits - (a) in applications where chopped fibers can replace continuous fibers for targeted strength and stiffness metrics, but at a substantially reduced cost; (b) they can tolerate deeper draws than continuous fiber composites in thermo-stamping and compression molding processes; (c) they can be tailored for pultrusion and unidirectional applications. Although pultrusion is primarily a process that adopts continuous fibers, stitch bonded entangled discontinuous fibers can provide unique intermediates. This is analogous to natural coir fibers which get aligned and entangled to produce ropes/rods for example, (d) they can be processed in cross-ply and multi-directional formats, like composite laminates. In this work Neenah Paper partnered with IACMI, UT and ORNL to evaluate structure-process-property relationships with Zoltek carbon fiber. A few process parameters such as machine speed, weight basis, fiber length, effect of fiber sizing, direction of mat lay-up etc. were investigated. The produced mats were converted to thermoplastic composite laminates using polyamide 6 (PA6, nylon) resin. The specific objective of this project is to produce a wet-laid nonwoven carbon fiber mat with a high degree of unidirectional fiber alignment, using discontinuous carbon fibers. The report provides details about the processing, characterization, and lower-upper bound properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗
DOE OSTI2022

Epitaxial Sc x Al 1- x N on GaN exhibits attractive high-K dielectric properties

Epitaxial Sc x Al 1- x N thin films of ~100 nm thickness grown on metal polar GaN substrates are found to exhibit significantly enhanced relative dielectric permittivity (ε r ) values relative to AlN. ε r values of ~17–21 for Sc mole fractions of 17%–25% ( x = 0.17–0.25) measured electrically by capacitance–voltage measurements indicate that Sc x Al 1- x N has the largest relative dielectric permittivity of any existing nitride material. Since epitaxial Sc x Al 1- x N layers deposited on GaN also exhibit large polarization discontinuity, the heterojunction can exploit the in situ high-K dielectric property to extend transistor operation for power electronics and high-speed microwave applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Terahertz near-field imaging of sidewall-induced losses in superconducting qubits

Correlating superconducting qubit performance with advanced materials analysis is a key strategy for improving coherence. Existing diagnostics for key properties, such as dielectric loss, structural discontinuity, and interface heterogeneity, often rely on destructive electron microscopy or low-throughput millikelvin measurements. Here, in this study, we demonstrate noninvasive terahertz (THz) nano-imaging/spectroscopy of encapsulated niobium transmon qubits as a high-throughput proxy for performance evaluation. We identify large variations in sidewall near-field signals, implicating sidewall loss and discontinuity as major coherence limiters, and also use THz hyperspectral line scans to probe dielectric responses and field participation at Al junction interfaces.

Kim, Richard H. J. [Ames Laboratory (AMES), Ames, ↗
DOE OSTI2022

Performance of discontinuity-free components produced by additive turning computer aided manufacturing strategy

Computer aided manufacturing (CAM) techniques for directed energy deposition (DED) affect the material properties of the manufactured component based on the scan strategy used. In this study, we investigate the material characteristics of turning-style toolpath strategies to generate axisymmetric components with additive manufacturing (AM), referred to in this research as additive turning. This novel approach leverages existing CAM technology for turning, where the component rotates around a stationary cutting tool, to generate toolpath trajectories for DED with varying wall-thicknesses and controlled deposition angles. This strategy allows for entire components to be deposited in one continuous deposition, resulting in reduced cycle-time and improved material usage efficiency compared to conventional AM strategies where the beam is switched off at the end of every layer. Results from this study show that the use of additive turning can produce over 99 % dense components with less variation and anisotropy in texture and hardness, as well as a lower variation in elongation to failure when compared to conventional strategies. This research highlights that various CAM strategies could be deployed for AM to improve process efficiency or enable localized control over part performance.

36 MATERIALS SCIENCE↗
DOE OSTI2020

HDG-1 Graphite Preirradiation Data Package Report

This report documents all pre-irradiation examination material-property measurement data for graphite specimens that are going to be used within the first high dose graphite (HDG) -1 irradiation capsule. The two new HDG capsules signify a major change to the AGC Experiment. HDG-1 and HDG-2 will replace the last two Advanced Graphite Creep (AGC) capsules (AGC-5 and AGC-6) which were designed to irradiate graphite at the extreme upper operational temperatures for a very-high-temperature reactor (VHTR) design, 1100°C. These very high temperature AGC-5 and AGC-6 capsules have been repurposed to re-irradiated specimens (from AGC-2, AGC-3, and AGC-4) at the lower temperatures of 600°C and 800°C. HDG-1 will be irradiated at 600°C and HDG-2 will be irradiated at 800°C. By re-irradiating the previous AGC specimens a total maximum neutron dose of around 15 dpa (displacements per atom) can be achieved for all major graphite grades at irradiation temperatures of 600°C and 800°C. Specimens in the HDG-1 capsule are made up of previously irradiated specimens from the AGC-2 capsule and unirradiated specimens prepared for the now discontinued AGC-5 capsule. Utilizing the irradiated specimens, a maximum neutron dose of around 15 dpa is anticipated. These new maximum dose levels will provide irradiated material property data over a total neutron dose range of 1-15 dpa at a temperature of 600°C when combined with the previous AGC-1 and AGC-2 irradiation data. This will provide quantitative data necessary for predicting the irradiation behavior and operating performance of new nuclear graphite grades for use within high temperature reactor designs. Similar to previous AGC test trains, HDG-1 includes the major graphite grades (IG-110, NBG-17, NBG-18, PCEA, and 2114) as well as adding the very fine-grain grade IG-430 which is of interest to the Molten Salt Reactor (MSR) designs. Also new to the HDG-1 capsule are 90 smaller geometry specimens designated as pencil specimens. These specimens take up only one third the space of a standard creep size specimen. This increased number of specimens will enhance property measurement statistics because they will provide 3 times the control specimen data at a position that would otherwise only have a single measurement.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

HDG-1 Graphite Preirradiation Data Package Report

This report documents all pre-irradiation examination material-property measurement data for graphite specimens that are going to be used within the first high dose graphite (HDG) -1 irradiation capsule. The two new HDG capsules signify a major change to the AGC Experiment. HDG-1 and HDG-2 will replace the last two Advanced Graphite Creep (AGC) capsules (AGC-5 and AGC-6) which were designed to irradiate graphite at the extreme upper operational temperatures for a very-high-temperature reactor (VHTR) design, 1100°C. These very high temperature AGC-5 and AGC-6 capsules have been repurposed to re-irradiated specimens (from AGC-2, AGC-3, and AGC-4) at the lower temperatures of 600°C and 800°C. HDG-1 will be irradiated at 600°C and HDG-2 will be irradiated at 800°C. By re-irradiating the previous AGC specimens a total maximum neutron dose of around 15 dpa (displacements per atom) can be achieved for all major graphite grades at irradiation temperatures of 600°C and 800°C. Specimens in the HDG-1 capsule are made up of previously irradiated specimens from the AGC-2 capsule and unirradiated specimens prepared for the now discontinued AGC-5 capsule. Utilizing the irradiated specimens, a maximum neutron dose of around 15 dpa is anticipated. These new maximum dose levels will provide irradiated material property data over a total neutron dose range of 1-15 dpa at a temperature of 600°C when combined with the previous AGC-1 and AGC-2 irradiation data. This will provide quantitative data necessary for predicting the irradiation behavior and operating performance of new nuclear graphite grades for use within high temperature reactor designs. Similar to previous AGC test trains, HDG-1 includes the major graphite grades (IG-110, NBG-17, NBG-18, PCEA, and 2114) as well as adding the very fine-grain grade IG-430 which is of interest to the Molten Salt Reactor (MSR) designs. Also new to the HDG-1 capsule are 90 smaller geometry specimens designated as pencil specimens. These specimens take up only one third the space of a standard creep size specimen. This increased number of specimens will enhance property measurement statistics because they will provide 3 times the control specimen data at a position that would otherwise only have a single measurement.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗