Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability 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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

The role of memory in non-genetic inheritance and its impact on cancer treatment resistance

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Empirical estimation of densities in NaCl-KCl-UCl 3 and NaCl-KCl-YCl 3 molten salts using Redlich-Kister expansion

Densities of molten KCl-NaCl-UCl 3 and KCl-NaCl-YCl 3 ternary systems have been estimated using a multidimensional Redlich-Kister model. Temperature and composition dependent Redlich-Kister functions have been used to generate binary interaction parameters in the outlined ternary salt systems. These binary interactions have been used in the extrapolation to ternary system densities. The results of the density extrapolations by Muggianu interpolation scheme provide agreement within 2–3% for the NaCl-KCl-YCl 3 liquids and 11% in NaCl-KCl-UCl 3 liquids compared to the available experimental data. Modeling NaCl-KCl-UCl 3 molten phase density with a ternary interaction parameter improved the agreement within 4%. Thermophysical modeling used in this study has shown promising results for use in other material properties, such as viscosity, thermal conductivity, and heat capacity of the molten salts. Lastly, the outlined modeling method applied in these specific molten salt ternaries can be used for quaternary or higher multicomponent molten salt systems.

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

Thermo4PFM: Facilitating Phase-field simulations of alloys with thermodynamic driving forces

Phase-field modeling is a popular front-tracking approach used to model solidification. Its time-evolution equations are often coupled to alloy composition and/or thermal diffusion in high-resolution multiphysics approaches. Materials thermodynamic properties tabulated in CALPHAD databases can be used for phase-field modeling to parameterize bulk energies of alloys. In addition, they can be naturally integrated into models such as the Kim-Kim-Suzuki (KKS) model where driving forces depend on the differences between chemical potentials of co-existing phases. In that case, a small system of coupled nonlinear equations needs to be solved at every point in space where the phase-field order parameter is to be updated and evolved in time. Here we present Thermo4PFM, a solver for the KKS equations for binary and ternary alloys, with two or three phases, and parameterized with CALPHAD models. Thermo4PFM is open source, written in C++, and can take advantage of Graphics Processing Units (GPU) accelerators. Using OpenMP offload capabilities for C++ classes, an excellent performance is demonstrated on GPU using the LLVM compiler. CALPHAD data is read from simple JSON files using an open source parser from the boost library.

36 MATERIALS SCIENCE↗
DOE OSTI2026

The Role of Normal Stress and Shear Stress Heterogeneity in the Inferred Depth‐Independence of Stress Drop

Earthquake stress drops are inferred to be independent of source depth, contradicting linear scaling predictions for earthquakes as frictional stick‐slip instabilities that assume increasing fault normal stress due to overburden. Here, we examine the scaling between averaged stress drops and increasing normal stress for simulated earthquake sequences in continuum rate‐and‐state fault models. The models produce a weaker dependence of stress drop on normal stress than the linearity of simple friction, which can be well‐fit by a sublinear power‐law. This result is more prominent when the fault dimension is much larger than nucleation scales. In such cases, the averaged behavior of ruptures is dominated by rupture propagation conditions, reflecting more heterogeneous shear stress conditions. As natural faults can be considerably larger than the smallest earthquakes they host, such a weaker scaling between averaged rupture conditions and normal stress may partially explain the lack of an inferred depth‐dependence of earthquake stress drops.

earthquake source↗
DOE OSTI2024

Relativistic corrections to the correlated basis function effective nuclear Hamiltonian

We discuss the inclusion of relativistic boost corrections into the correlated basis function effective nuclear Hamiltonian, derived from a realistic model of two- and three-nucleon interactions using the formalism of correlated basis functions and the cluster expansion technique. Different procedures to take into account the effects of boost interactions are compared on the basis of the ability to reproduce the nuclear matter equation of state obtained from accurate quantum many-body calculations. Furthermore, the results of our study show that the repulsive contribution of the boost interaction significantly depends on the underlying model of the nonrelativistic potential. On the other hand, the dominant relativistic correction turns out to be the corresponding reduction of the strength of repulsive three-nucleon interactions, leading to a significant softening of the equation of state of nuclear matter at supranuclear densities.

Neutron stars & pulsars↗
DOE OSTI2025

Anomalous softening of 3D printed elastomeric foam irradiated under compressive strain

Elastomeric foam is an essential component in many industrial and technological settings, primarily as thermal insulators and as positional/mechanical support cushions. In particular, silicone foam is utilized in harsh environments due to exceptional thermal and chemical stability. Under service conditions within certain applications such material gets exposed to a high dosage of gamma radiation, which can permanently alter the material’s structural and mechanical response properties. Most studies on gamma-exposure under inert or oxidative atmosphere indicate hardening of silicone foam, which is attributed to radiation-induced enhancement in chemical cross-linking. Here we report two contrasting effects depending on whether (non-oxidative) radiation exposure is carried out with the foam under zero or finite compressive strain. While in the former case we observe radiation-hardening consistent with previous studies, in the latter case (50% porous foam under 30% uniaxial compression) we see a monotonic decrease in Young’s modulus with increasing dosage, although solvent swelling experiments on the constituent rubber indicate a net increase in cross-link density independent of the state of strain. We quantitatively model all dose-dependent data using the Ogden Hyperfoam strain-energy function within the framework of Tobolsky two-network scheme and attribute the above anomaly to a combined effect of radiation-induced thickness change (compression set) and inherent nonlinearity in the foam’s stress-strain response.

Coarse-grained models↗
DOE OSTI2020

Effect of Cr(III) Adsorption on the Dissolution of Boehmite Nanoparticles in Caustic Solution

Incorporation of relatively minor impurity metals onto metal (oxy)hydroxides can strongly impact solubility. In complex highly alkaline multicomponent tank waste such as that at the Hanford Nuclear Reservation, tests indicate that the surface area-normalized dissolution rate of boehmite ($γ$-AlOOH) nanomaterials is at least an order of magnitude lower than that predicted for the pure phase. In this work we examine the dissolution kinetics of boehmite coated by adsorbed Cr(III), which adheres at saturation coverages as sparse chemisorbed monolayer clusters. Using 40 nm boehmite nanoplates as a model system, temperature-dependent dissolution rates of pure versus Cr(III)-adsorbed boehmite showed that the initial rate for the latter is consistently several times lower, with an apparent activation energy 16 kJ·mol -1 higher. Although the surface coverage is only around 50%, solution analytics coupled to multimethod solids characterization reveals a phyicochemical armoring effect by adsorbed Cr(III) that substantially reduces the number of dissolution-active sites on particle surfaces. Such findings could help improve kinetics models of boehmite and/or metal ion adsorbed boehmite nanomaterials, ultimately providing a stronger foundation for the development of more robust complex radioactive liquid waste processing strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A redefinition of the halo boundary leads to a simple yet accurate halo model of large-scale structure

ABSTRACT We present a model for the halo–mass correlation function that explicitly incorporates halo exclusion and allows for a redefinition of the halo boundary in a flexible way. We assume that haloes trace mass in a way that can be described using a single scale-independent bias parameter. However, our model exhibits scale-dependent biasing due to the impact of halo-exclusion, the use of a ‘soft’ (i.e. not infinitely sharp) halo boundary, and differences in the one halo term contributions to ξhm and ξmm. These features naturally lead us to a redefinition of the halo boundary that lies at the ‘by eye’ transition radius from the one-halo to the two-halo term in the halo–mass correlation function. When adopting our proposed definition, our model succeeds in describing the halo–mass correlation function with $\approx 2{{\ \rm per\ cent}}$ residuals over the radial range 0.1 h−1 Mpc < r < 80 h−1 Mpc, and for halo masses in the range 1013 h−1 M⊙ < M < 1015 h−1 M⊙. Our proposed halo boundary is related to the splashback radius by a roughly constant multiplicative factor. Taking the 87 percentile as reference we find rt/Rsp ≈ 1.3. Surprisingly, our proposed definition results in halo abundances that are well described by the Press–Schechter mass function with δsc = 1.449 ± 0.004. The clustering bias parameter is offset from the standard background-split prediction by $\approx 10{{\ \rm per\ cent}}\!-\!15{{\ \rm per\ cent}}$. This level of agreement is comparable to that achieved with more standard halo definitions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Formation of SMBH seeds in Population III star clusters through collisions: the importance of mass loss

ABSTRACT Runaway collisions in dense clusters may lead to the formation of supermassive black hole (SMBH) seeds, and this process can be further enhanced by accretion, as recent models of SMBH seed formation in Population III star clusters have shown. This may explain the presence of SMBHs already at high redshift, z > 6. However, in this context, mass loss during collisions was not considered and could play an important role for the formation of the SMBH seed. Here, we study the effect of mass loss, due to collisions of protostars, in the formation and evolution of a massive object in a dense primordial cluster. We consider both constant mass-loss fractions as well as analytic models based on the stellar structure of the collision components. Our calculations indicate that mass loss can significantly affect the final mass of the possible SMBH seed. Considering a constant mass loss of 5 per cent for every collision, we can lose between 60–80 per cent of the total mass that is obtained if mass loss were not considered. Using instead analytical prescriptions for mass loss, the mass of the final object is reduced by 15–40 per cent, depending on the accretion model for the cluster we study. Altogether, we obtain masses of the order of $10^4\, \mathrm{M}_{\odot }$, which are still massive enough to be SMBH seeds.

Alister Seguel, P. J.↗
DOE OSTI2024

Rapid assessment of the creep rupture life of metals: A model enabling experimental design

Prediction of the creep rupture life of engineering metals is critical for qualification and design of new materials. The use of long-term creep tests and the need to quantify the performance variability in a priori similar systems hinder the rapid creep assessment of a given material. Therefore, it is essential to develop methods that can extrapolate the long-term performance of alloys and the associated variability from short-term experiments. To this end, this study introduces a new model which enables the estimation of the rupture life of a material for a given stress and temperature. This model relies on two components. First, a new relation for the minimum creep rate (MCR) of materials is introduced. It includes a stress dependent stress exponent allowing the model to capture the variation of MCR across a wide range of temperatures and stresses. Second, employing the Monkman-Grant (MG) law, we establish a relation between stress, temperature and creep rupture life. Together, these two elements yield a new closed-form mathematical expression for the Larson Miller parameter as a function of stress and temperature. This expression captures the creep rupture time for many metals (Gr91, Copper, Gr122 and 347H) and compares favorably with alternate empirical approaches. The model is then used to assess the minimum duration of creep rates necessary to qualify the material up to 100000h. Furthermore, it is found that depending on the material system, creep tests as few as five limited to 5000 h for steels (Gr91, Gr122, 347H) and 100 h for copper are sufficient to model creep lifetimes. Finally, using a Bayesian inference-based approach to calibrate the model, we demonstrate that variability in rupture life can be captured via the quantification of the uncertainty in the model parameters and extrapolated from a limited number of short to moderately short creep tests; thereby paving the way for accelerated creep testing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Aggregation of Leslie matrix models with application to ten diverse animal species

Some grouping is necessary when constructing a Leslie matrix model because it involves discretizing a continuous process of births and deaths. Here, the level of grouping is determined by the number of age classes and frequency of sampling. It is largely unknown what is lost or gained by using fewer age classes, and I address this question using aggregation theory. I derive an aggregator for a Leslie matrix model using weighted least squares, determine what properties an aggregated matrix inherits from the original matrix, evaluate aggregation error, and measure the influence of aggregation on asymptotic and transient behaviors. To gauge transient dynamics, I employ reactivity of the standardized Leslie matrix. I apply the aggregator to 10 Leslie models developed for animal populations drawn from a diverse set of species. Several properties are inherited by the aggregated matrix: (a) it is a Leslie matrix; (b) it is irreducible whenever the original matrix is irreducible; (c) it is primitive whenever the original matrix is primitive; and (d) its stable population growth rate and stable age distribution are consistent with those of the original matrix if the least squares weights are equal to the original stable age distribution. In the application, depending on the population modeled, when the least squares weights do not follow the stable age distribution, the stable population growth rate of the aggregated matrix may or may not be approximately consistent with that of the original matrix. Transient behavior is lost with high aggregation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗
DOE OSTI2021

Self-induced spin-orbit torques in metallic ferromagnets

Herein we present a phenomenological theory of spin-orbit torques in a metallic ferromagnet with spin-relaxing boundaries. The model is rooted in the coupled diffusion of charge and spin in the bulk of the ferromagnet, where we account for the anomalous Hall effects as well as the anisotropic magnetoresistance in the corresponding constitutive relations for both charge and spin sectors. The diffusion equations are supplemented with suitable boundary conditions reflecting the spin-sink capacity of the environment. In inversion-asymmetric heterostructures, the uncompensated spin accumulation exerts a dissipative torque on the order parameter, giving rise to a current-dependent linewidth in the ferromagnetic resonance with a characteristic angular dependence. We compare our model to recent spin-torque ferromagnetic resonance measurements, illustrating how rich self-induced spin-torque phenomenology can arise even in simple magnetic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Microscopic study of the coupled-wire construction and plausible realization in spin-dependent optical lattices

Coupled-wire constructions offer particularly simple and powerful models to capture the essence of strongly correlated topological phases of matter. They often rely on effective theories valid in the low-energy and strong-coupling limits, which impose severe constraints on the physical systems where they could be realized. In this work, we investigate the microscopic relevance of a class of coupled-wire models and their possible experimental realization in cold-atom experiments. We connect with earlier results and prove the emergence of fractional quantum Hall states in the limit of strong interwire tunneling. Contrary to previous studies relying on renormalization group arguments, our microscopic approach exposes the connection between coupled-wire constructions and model wave functions in continuum Landau levels. Then, we use exact-diagonalization methods to investigate the appearance of these fractional quantum Hall states in more realistic settings. We examine the parameter regimes where these strongly correlated phases arise, and provide a way to detect their appearance in cold-atom experiments through standard time-of-flight measurements. Motivated by this experimental probe, we finally propose a realization of our model with cold atoms in spin-dependent optical lattices. Our estimates show that the previous fractional quantum Hall phases lie within experimentally accessible parameter regimes, giving a viable route toward their experimental study.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

On small black holes in string theory

Abstract We discuss the worldsheet sigma-model whose target space is thed+1 dimensional Euclidean Schwarzschild black hole. We argue that in the limit where the Hawking temperature of the black hole,T, approaches the Hagedorn temperature,T H , it can be described in terms of a generalized version of the Horowitz-Polchinski effective theory. Ford≥ 6, where the Horowitz-Polchinski EFT [1, 2] does not have suitable solutions, the modified effective Lagrangian allows one to study the black hole CFT in an expansion in powers ofd− 6 andT H −T. AtT=T H , the sigma model is non-trivial for alld> 6. It exhibits an enhanced SU(2) symmetry, and is described by a non-abelian Thirring model with a radially dependent coupling. The resulting picture connects naturally to the results of [3–5], that relate Schwarzschild black holes in flat spacetime at largedto the two dimensional black hole. We also discuss an analogous open string system, in which the black hole is replaced by a system of two separated D-branes connected by a throat. In this system, the asymptotic separation of the branes plays the role of the inverse temperature. At the critical separation, the system is described by a Kondo-type model, which again exhibits an enhanced SU(2) symmetry. At larged, the brane system gives rise to the hairpin brane [6].

Physics↗
DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗
DOE OSTI2024

Detachment scalings derived from 1D scrape-off-layer simulations

Fusion power plants will require detachment to mitigate sputtering and keep divertor heat fluxes at tolerable levels. Controlling detachment on these devices may require the use of real-time scrape-off-layer modeling to complement the limited set of available diagnostics. In this work, we use the configurable Hermes-3 edge modeling framework to perform time-dependent, fixed-fraction-impurity 1D detachment simulations. Although currently far from real-time, these simulations are used to investigate time-dependent effects and the minimum physics set required for control-relevant modeling. We show that these simulations reproduce the expected rollover of the target ion flux — a typical characteristic of detachment onset. We also perform scans of the input heat flux and impurity concentration and show that the steady-state results closely match the scalings predicted by the 0D time-independent Lengyel–Goedheer model. This allows us to indirectly compare to SOLPS simulations, which find a similar scaling but a lower value for the impurity concentration required for detachment for given upstream conditions. We use this result to suggest a series of improvements for the Hermes simulations, and finally show simulations demonstrating the impact of time-dependence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Emergence of the isotropic Kitaev honeycomb lattice $α$–RuCl 3 and its magnetic properties

Here, we present a comprehensive investigation of the crystal and magnetic structures of the van der Waals antiferromagnet $α$–RuCl 3 using single crystal x-ray and neutron diffraction. The crystal structure at room temperature is a monoclinic (C2/m). However, with decreasing temperature, a remarkable first-order structural phase transition is observed, leading to the emergence of a rhombohedral ($R\bar3$) structure characterized by three-fold rotational symmetry forming an isotropic honeycomb lattice. On further cooling, a zigzag-type antiferromagnetic order develops below $T_{\textrm{N}} = 6\sim6.6$ K. The critical exponent of the magnetic order parameter was determined to be $\beta = 0.11(1)$, which is close to the two-dimensional Ising model. Additionally, the angular dependence of the magnetic critical field of the zigzag antiferromagnetic order for the polarized ferromagnetic phase reveals a six-fold rotational symmetry within the ab–plane. These findings reflect the symmetry associated with the Ising-like bond-dependent Kitaev spin interactions and underscore the universality of the Kitaev interaction-dominated antiferromagnetic system.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗