Search NASA⌕ Search

SEARCH · Search NASA

Results for “material parameters”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Summary of Geometric Parameters and Their Effects on Performance of U-10Mo Fuel Plates

A monolithic plate-type fuel system has been under development to convert high-performance test reactors from highly enriched uranium to low-enrichment uranium fuels and is now moving into the qualification phase, a predecessor to the timely conversion of the target reactors. To qualify this fuel system, the plates must meet safety standards and perform well in a reactor. The plates must maintain mechanical integrity, exhibit geometric stability, and have stable and predictable in-reactor behavior. The requirement to maintain mechanical integrity under normal operating conditions is primarily demonstrated by successful testing. However, each high-performance reactor employs a distinct design, resulting in distinct plate geometries, with unique features, attributes, irregularities, and tolerances. Due to the abundance of such distinct geometric varieties, a single “generic” plate geometry capturing all extremes is not achievable. It is also impractical to test each of these proposed designs in a reactor. This limitation necessitates cautious evaluations since the thermomechanical response of a plate with a certain geometry may not be representative for a plate with a significantly different geometry. To address concerns related to in-reactor plate performance, large set of sensitivity studies were performed. These parametric studies aimed to better understand irradiation performance, while evaluating the sensitivity of results to various modeling inputs, including geometric, operational, and material parameters. This work studied selected geometric parameters based on provided fuel specifications and performed a series of parametric simulations. The resulting temperature, displacement, and stress strains were comparatively evaluated to determine the effects of various geometric parameters. This draft provides a “high-level summary” of our parametric sensitivity studies and summarizes the key findings.

42 ENGINEERING↗
DOE OSTI2022

Evaluation of accuracy and convergence of numerical coupling approaches for poroelasticity benchmark problems

Accurate modeling of subsurface flow and transport processes is vital as the prevalence of subsurface activities such as carbon sequestration, geothermal recovery, and nuclear waste disposal increases. Computational modeling of these problems leverages poroelasticity theory, which describes coupled fluid flow and mechanical deformation. Although fully coupled monolithic schemes are accurate for coupled problems, they can demand significant computational resources for large problems. In this work, a fixed stress scheme is implemented into the Sandia Sierra Multiphysics toolkit. Two implementation methods, along with the fully coupled method, are verified with one-dimensional (1D) Terzaghi, 2D Mandel, and 3D Cryer sphere benchmark problems. The impact of a range of material parameters and convergence tolerances on numerical accuracy and efficiency was evaluated. Overall the fixed stress schemes achieved acceptable numerical accuracy and efficiency compared to the fully coupled scheme. However, the accuracy of the fixed stress scheme tends to decrease with low permeable cases, requiring the finer tolerance to achieve a desired numerical accuracy. For the fully coupled scheme, high numerical accuracy was observed in most of cases except a low permeability case where an order of magnitude finer tolerance was required for accurate results. Finally, a two-layer Terzaghi problem and an injection–production well system were used to demonstrate the applicability of findings from the benchmark problems for more realistic conditions over a range of permeability. Simulation results suggest that the fixed stress scheme provides accurate solutions for all cases considered with the proper adjustment of the tolerance. Furthermore, this work clearly demonstrates the robustness of the fixed stress scheme for coupled poroelastic problems, while a cautious selection of numerical tolerance may be required under certain conditions with low permeable materials.

58 GEOSCIENCES↗
DOE OSTI2021

Deep Learning Analysis of Polaritonic Wave Images

Deep learning (DL) is an emerging analysis tool across the sciences and engineering. Encouraged by the successes of DL in revealing quantitative trends in massive imaging data, we applied this approach to nanoscale deeply subdiffractional images of propagating polaritonic waves in complex materials. Utilizing the convolutional neural network (CNN), we developed a practical protocol for the rapid regression of images that quantifies the wavelength and the quality factor of polaritonic waves. Using simulated near-field images as training data, the CNN can be made to simultaneously extract polaritonic characteristics and material parameters in a time scale that is at least 3 orders of magnitude faster than common fitting/processing procedures. The CNN-based analysis was validated by examining the experimental near-field images of charge-transfer plasmon polaritons at graphene/α-RuCl3 interfaces. Our work provides a general framework for extracting quantitative information from images generated with a variety of scanning probe methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quasiparticle conductance in spin valve Josephson structures

We study the quasiparticle current in clean ferromagnetic Josephson structures of the form $S_1/F_1/N/F_2/S_2$, where $\textit{S, F,}$ and $\textit{N}$ denote superconducting, ferromagnetic, or normal layers, respectively. Our focus is on the structure of the conductance $\textit{G}$ as a function of bias $\textit{V$, emphasizing the subgap region. We use a fully self-consistent numerical method, coupled to a transfer matrix procedure to extract $\textit{G(V)}$. We choose material parameters appropriate to experimentally realized Co Cu Nb structures. In this work, we find a resonance peak structure as a function of the intermediate layer thickness and of the misalignment angle $\phi$ between $\textit{F}$ layers. To understand this resonance structure, we develop an approximate analytic method. For experimentally relevant thicknesses, the conductance has multiple subgap peaks, which oscillate in position between low and critical bias positions. These oscillations occur in both $\phi$ and the layer thicknesses. We compare our results with those obtained for the spin valve structures ($F_1/N/F_2/S_2$) and discuss the implications of our results for the fabrication of spin Josephson devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

On plasticity-enhanced interfacial toughness in bonded joints

Here, the performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from -60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γ o depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γ o and σ*/σ yb with the value Γ/Γ o increasing rapidly above a threshold value of σ*/σ yb . The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γ o during micrometer-scale crack-growth when σ*/σ yb = 2 (a reasonable choice for σ*/σ yb ). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Diverse simulations of time-resolved photoluminescence in thin-film solar cells: A SnO 2 /CdSe y Te 1-y case study

Time-resolved photoluminescence (TRPL) is widely used to measure carrier lifetime in thin-film solar cell absorbers. However, the injection dependence of data and frequent non-exponential decay shapes complicate the interpretation. Here, we develop a numerical model to simulate injection-dependent TRPL measurements in a SnO 2 /CdSe y Te 1-y case study solar cell structure, considering parameters of interest to researchers in industry and academia. Previous simulations have shown that, in low injection, excess electrons and holes injected by the laser pulse are rapidly separated in the electric field formed by the pn junction. As a result, at early times, the PL signal can decay faster than the Shockley-Read-Hall lifetime in the absorber bulk (t_bulk). Prior simulations have shown that the charge stored in the junction can slowly leak out to affect decays at late times. However, it has not been clear if and to what degree charge storage can affect the slopes extracted from TRPL decays - t_2 - commonly cited as the TRPL-measured lifetime. Here, we show that charge storage can, in some cases, result in t_2 values that substantially overestimate t_bulk. Previous simulations indicate that high injection conditions can screen the junction field and minimize charge separation. Here, we show that continued injection increases can drive down t_2 below t_bulk as radiative recombination becomes dominant. We catalog charge storage and radiative recombination impacts for a diverse set of material parameters and compare results to double-heterostructure models.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,S↗
DOE OSTI2020

PBX9502 PAGOSA/SURF Calibration for 100 µm Grid Size

There are two methods to calibrate the SURF reactive burn explosives model.The first is to calibrate using Pop plot data. Because there are only two main parameters in SURF, called A and B, there is just enough flexibility in SURF to determine the Pop plot line, i.e., equivalent information to a slope and a single point on the line. However, due to the linearity of the Pop plot being based on logarithmic scales, a slight deviation off the Pop plot (though the fit to the Pop plot might be very good overall) can make considerable difference when comparing with the velocity profile data. The second calibration method thusly focuses strongly on matching entire velocity profile data. This is to say that individual velocity profiles extracted from embedded gages (in the case of PBX 9502, we refer to R. Gustavsen’s experimental data) may fit poorly, while the overall fit to the Pop plot may seem reasonable. In addition, there are other factors that could produce slightly different calibrations for SURF, while leading to an overall acceptable fit to data. One of those may be the physics code used and here we use Pagosa. Another is mesh size dependence. A third one stems from fitting a reactive burn model being done in conjunction with an equation of state (EOS) of the reactants and one for the products – different EOS can result in slightly different parameterizations. Finally, there are several secondary material parameters in SURF, and the overall fit will depend on the entire parameterization. The goal of this report is to clarify some of our previously made statements about SURF calibration, as well as to redo a systematic calibration using 100 µm grid resolution.

42 ENGINEERING↗
DOE OSTI2021

Dislocation transmission across Σ3{112} incoherent twin boundary: a combined atomistic and phase-field study

Grain boundaries (GBs) in polycrystalline materials act as impediments to dislocation motion and result in strengthening. Understanding slip transmission through GBs, specifically twin boundaries, is essential to understand the plastic deformation behavior of polycrystalline fcc materials. Here the interaction between a glide dislocation and Σ3{112} incoherent twin boundary (ITB) in copper is investigated using a combined atomistic and mesoscale approach. The material parameters and structure of the GB in the mesoscale phase field dislocation dynamics (PFDD) model are informed from Molecular Statics (MS) simulations. The structural unit of the ITB consists of an array of three partial dislocations. The interaction between a glide dislocation impinging on each of the GB partial dislocations is investigated using both PFDD and Molecular Dynamics (MD) with two boundary conditions. Transmission planes predicted by both PFDD and MD (NVT) are in agreement, and show that not all transmission events are direct. Critical transmission stresses predicted by PFDD are in the range of 276 MPa to 1380 MPa, while MD predictions are in the range from 100 MPa to 700 MPa. The PFDD and MD predictions of slip transmission are explained using dislocation theory based on isotropic linear elasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Neutron scattering off Weyl semimetals

We present how to detect type-I Weyl nodes in a material by inelastic neutron scattering. Such an experiment first of all allows one to determine the dispersion of the Weyl fermions. We extend the reasoning to produce a quantitative test of the Weyl equation, taking into account realistic anisotropic properties. These anisotropies are mostly contained in the form of the emergent magnetic moment of the excitations, which determines how they couple to the neutrons. Although there are many material parameters, we find several quantitative predictions that are universal and demonstrate that the excitations are described by solutions to the Weyl equation. The anisotropic coupling between electrons and neutrons implies that even fully unpolarized neutrons can reveal the spin-momentum locking of the Weyl fermions because the neutrons will couple to some components of the Weyl fermion pseudospin more strongly. On the other hand, in an experiment with polarized neutrons, the scattered neutron beam remains fully polarized in a direction that varies as a function of momentum transfer (within the range of validity of the Weyl equation). This allows measurement of the chirality of Weyl fermions for inversion-symmetric nodes. Furthermore, in this work we estimate that the scattering rate may be large enough for such experiments to be practical; in particular, the magnetic moment may be larger than the ordinary Bohr magneton, compensating for a small density of states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Smart material based multilayered microbeam structures for spatial self-deployment and reconfiguration: A residual stress approach

Alleviation of the potentially damaging effects induced by residual stresses was comprehensively investigated in previous research. Here, this paper, however, presents a spatially self-deployable and reconfigurable multilayered microbeam which takes advantage of residual stresses and shape memory effects. Reconfigurable mechanism of a typical four-layered microbeam composed of Pt\Ni 50 Ti 50 \Ni 50 Ti 50 \Pt is introduced, followed by analytical modeling of the maximum distance of the self-deployed gap as functions of variable structural and material parameters, including compressive residual stress in Pt layers and tensile residual stress in Ni 50 Ti 50 layers. Analytical solutions given by the static model agree well with the results obtained via finite element models (FEMs). Fabrication, characterization, and in-situ experiments were carried out to validate the feasibility of deployment of the as-released four-layered microbeam. The maximum distance of the gap was measured to be 41.39 μm at 20 °C, which could be increased to 51.73 μm thanks to controllable reconfiguration driven by shape memory effects. Theoretical analysis of such self-deployment and reconfiguration suggested a tensile residual stress increase by 52 MPa in Ni 50 Ti 50 layers. The multilayered microbeam structure with capabilities of self-deployment and reconfiguration offers great potential for various emerging applications, such as micro robotics, medical drug delivery devices, and intelligent chip scale spacecraft.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quantum efficiency of InGaN–GaN multi-quantum well solar cells: Experimental characterization and modeling

InGaN-based multi-quantum well (MQW) solar cells are promising devices for photovoltaics (e.g., for tandem solar cells and concentrator systems), space applications, and wireless power transfer. In order to improve the efficiency of these devices, the factors limiting their efficiency and stability must be investigated in detail. Due to the complexity of a MQW structure, compared with a simple pn junction, modeling the spectral response of these solar cells is not straightforward, and ad hoc methodologies must be implemented. In this paper, we propose a model, based on material parameters and closed-formula equations, that describes the shape of the quantum efficiency of InGaN/GaN MQW solar cells, by taking into account the layer thickness, the temperature dependence of the absorption coefficient, and quantum confinement effects. We demonstrate (i) that the proposed model can effectively reproduce the spectral response of the cells; in addition, (ii) we prove that the bulk p-GaN layer absorbs radiation, but the carriers photogenerated in this region do not significantly contribute to device current. Finally, we show that (iii) by increasing the temperature, there is a redshift of the absorption edge due to bandgap narrowing, which can be described by Varshni law and is taken into account by the model, and a lowering in the extraction efficiency due to the increase in recombination (mostly Shockley–Read–Hall) inside the quantum wells, which is also visible by decreasing light intensity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

A molecular dynamics survey of grain boundary energy in uranium dioxide and cerium dioxide

Uranium dioxide (UO 2 ) is the primary fuel material that is used in current nuclear reactors. As one of the most fundamental material parameters, grain boundary (GB) energy strongly influences many fuel properties, and the influences depend on the characters and properties of individual GBs. Using molecular dynamics simulations, a high throughput survey of GB energy in UO 2 was carried out for the purpose of elucidating the roles of GB geometry such as misorientation and inclination, as well as the bonding nature of UO 2 , in affecting GB energy. GB energies in CeO 2 were calculated as well for comparison with UO2toinvestigate the generality of GB energy anisotropy in fluorite phase oxides. The results show significant GB energy anisotropy in both UO 2 and CeO 2 that is associated with the cubic symmetry of the fluorite structure. More interestingly, the GB anisotropy is found to be dependent not only on the crystal structure but also the ionic bonding. As such, the GB energy anisotropy in fluorite oxides has significant differences compared with that in fcc metals. We report that the data obtained and the increased knowledge on GB anisotropy will facilitate GB engineering for nuclear fuels with improved properties

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Durable and Low-Cost Fractal Structured Multifunctional Coatings for Next-Generation CSP

Fractal coatings on various low-cost and high-temperature metallic substrates such as SS316, SS347, In800H, In740H, In625, Haynes 230, etc. were fabricated employing our patented novel electrodeposition process at various deposition potentials and deposition times, as well as chemical etching. The best combination of the materials and processing conditions were co-optimized for the reliable and repeatable fabrication of the multifunctional coatings. Surface morphologies were examined using SEM and AFM, and fractal characterization was done using a 3D optical surface profiler. Fractal dimension was used to identify successful fabrication and material parameters. The mechanical durability of the coatings was evaluated using several standardized tests, such as sand abrasion (ASTM D968) and adhesion (ASTM D3359), to fully characterize and demonstrate the mechanical durability of the fabricated surfaces. HTF immersion tests on the fractal multifunctional coatings were conducted based on a modification of the ASTM D870, which is established for static water immersion testing. To calculate the corrosion rates and descaled weight loss of the samples subjected to thermal exposure for different time durations (100 h, 300 h, 500 h, and optionally 750 h) to estimate the corrosion rate. While protecting the substrates from corrosion is an important goal, an allied goals is also to examine the post-corrosion temperature stability and properties of the molten salts themselves due to the dissolution of the constituents of the alloys and the corrosion products from the coatings that cause “fouling” of the molten salts. Fouling characteristics of the different molten salt HTFs exposed to coated and uncoated substrates were investigated. Finally, considering a prototype molten salt/sCO2 heat exchanger, cost/performance model is developed to demonstrate that the coatings on low-cost alloys such as stainless steel provide lower levelized cost of heat exchanger compared to high nickel content Ha230.

14 SOLAR ENERGY↗
DOE OSTI2023

MHz free electron laser x-ray diffraction and modeling of pulsed laser heated diamond anvil cell

A new diamond anvil cell experimental approach has been implemented at the European x-ray Free Electron Laser, combining pulsed laser heating with MHz x-ray diffraction. In this report we use this setup to determine liquidus temperatures under extreme conditions, based on the determination of time-resolved crystallization. The focus is on a Fe-Si-O ternary system, relevant for planetary cores. This time-resolved diagnostic is complemented by a finite-element model, reproducing temporal temperature profiles measured experimentally using streaked optical pyrometry. This model calculates the temperature and strain fields by including (i) pressure and temperature dependencies of material properties, and (ii) the heat-induced thermal stress, including feedback effect on material parameter variations. Making our model more realistic, these improvements are critical as they give 7000 K temperature differences compared to previous models. Laser intensities are determined by seeking minimal deviation between measured and modeled temperatures. Combining models and streak optical pyrometry data extends temperature determination below detection limit. The presented approach can be used to infer the liquidus temperature by the appearance of SiO 2 diffraction spots. In addition, temperatures obtained by the model agree with crystallization temperatures reported for Fe–Si alloys. Our model reproduces the planetary relevant experimental conditions, providing temperature, pressure, and volume conditions. Those predictions are then used to determine liquidus temperatures at experimental timescales where chemical migration is limited. This synergy of novel time-resolved experiments and finite-element modeling pushes further the interpretation capabilities in diamond anvil cell experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Revisiting the Role of Mobile Ion Size for the Activation Barrier of Charge Transport in Single Ion Conducting Polymers

The development of polymer electrolytes for energy storage and conversion technologies requires a fundamental understanding of the material parameters controlling the energy barriers for ion transport. In glassy polymers, these activation barriers are usually extracted by using Arrhenius procedures. However, our recent studies of single ion conducting polymers reveal that this traditional Arrhenius description provides anomalously small prefactors, an issue that is widely disregarded in the literature. Here, to address it, an alternative approach was introduced for extracting the effective activation barrier E by imposing a physically valid limit (∼10 –13 s) for the ion hopping attempt time. Two consequences of this recent approach are that in polymer electrolytes (i) E is significantly (∼30–40%) lower than the barrier estimated using traditional Arrhenius fits and (ii) E displays significant temperature dependence even in the glassy state. Under these circumstances, we are revisiting the role of ion size, glass transition temperature T g , and dielectric and elastic constants of the polymer matrix on the effective activation barriers for ion transport as extracted using the recently proposed approach and highlight the differences from the picture that previously emerged based on simplistic apparent Arrhenius analysis. To this end, we investigate cation transport in three families of single ion conducting trifluoromethane sulfonimide-based polyanions with varied T g . Our results indicate that E decreases with mobile cation size and is highly sensitive to changes in the dielectric permittivity of the matrix, even for large cations. These insights call for revisions of many earlier results based on apparent Arrhenius fits as the proposed approach can provide more accurate guidance for the design of polymer electrolytes with enhanced ionic conductivity.

energy barrier↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗