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 739 records · Page 41

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗

Sensitivities of time-dependent temperature profile predictions for NSTX with the multi-mode model

The Multi-Mode Model (MMM) for turbulent transport was applied to a large set of well-analyzed discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate its sensitivities to a wide range of plasma conditions. MMM calculations were performed for hundreds of milliseconds in each discharge by performing time-dependent predictive simulations with the 1.5D tokamak integrated modeling code TRANSP. A closely related study (Lestz et al 2025 Plasma Phys. Control. Fusion 67 105029) concluded that MMM predicted electron and ion temperature profiles that were in reasonable agreement with NSTX observations, generally outperforming a different reduced transport model, TGLF. This finding motivates the more thorough investigation of the characteristics of the MMM predictions conducted in this work. The simulations with MMM have electron energy transport dominated by electron temperature gradient modes in the examined discharges with relatively low plasma β (ratio of kinetic plasma pressure to magnetic field pressure) and high collisionality, transitioning to a mixture of different modes for higher β and lower collisionality. The thermal ion diffusivity predicted by MMM is much smaller than the neoclassical contribution, in line with previous experimental analysis of NSTX. Nonetheless, the electron and ion temperature profiles are coupled via collisional energy exchange and thus sensitive to which transport channels are predicted. The time-dependent simulations with MMM are robust to the simulation start time, converging to remarkably similar temperature profiles later during the discharge. MMM typically overpredicts confinement relative to NSTX observations, leading to the prediction of overly steep temperature profiles. Plasmas with spatially broader temperature profiles, higher plasma β, and longer energy confinement times tend to be predicted by MMM with better agreement with the experiment. As a result, these findings provide useful context for understanding the regime-dependent tendencies of MMM in anticipation of self-consistent, time-dependent predictive simulations of NSTX-U discharges with these same modeling tools.

MMM↗

Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment. Dataset.

To simulate the effects of mycorrhizal fungi on soil organic matter cycling, we incorporated mycorrhizal processes into the Carbon, Organisms, Rhizosphere, and Protection in the Soil Environment (CORPSE) model to develop a new soil model Myco-CORPSE. The new model was calibrated and evaluated against soil measurements taken at temperate forests in New Hampshire (NH) and Georgia (GA). A series of scenario analysis were also conducted to explore the conditions under which ectomycorrhizal (ECM) N acquisition processes can induce different soil C accumulation in ECM systems compared to arbuscular (AM) systems.In this data package, we included:-The Python codes of the standard Myco-CORPSE model we developed: "Standard Myco_CORPSE python codes.zip". The main program is the "gradient_sim.py" which calculates the bulk soil microbes and CN content along a user defined gradient of clay, soil temperature, soil moisture and mycorrhizal dominance, and relies on two subprograms "CORPSE_deriv.py" and "CORPSE_integrate.py". "CORPSE_deriv.py" calculated the changes in all simulated soil stock within every time step and "CORPSE_integrate.py" integrate the changes in all simulated soil stock within simulated time period. The program "Plot.py" is used to plot the major outputs produced by the main program "gradient_sim.py".-The modified Python codes of Myco-CORPSE models with site-level environmental inputs (in NH and GA) used to conduct simulations in NH and GA sites: "NH_GA model simulations.zip". -The Python codes used to evaluate the Myco-CORPSE simulation outputs in NH and GA sites against site-level measurements: "Plot NH_GA simulation against measurements.zip". It includes both the evaluation Python code, the model outputs on NH and GA sites, and the measured soil properties in both sites.-The modified Python codes of Myco-CORPSE models "Scenario analysis_model simulations.zip" that is used to conduct scenario analysis of how different litter properties, mycorrhizal fungal traits, climate, and seasonal variation in temperature and vegetation phenology impact the mycorrhizal effects on soil CN properties. The sub file folder "Scenario analysis_litter traits" contains the codes for scenario analysis of different litter properties; The sub file folder "Scenario analysis_ECM types" contains the codes for scenario analysis of different ECM fungal traits; The sub file folder "Scenario analysis_climate&seasonality" contains the codes for scenario analysis of different climate and seasonalities;-"Scenario analysis_model results and plotting codes.zip" contains all the output files from the the scenario analysis of Myco-CORPSE model as described above and the plotting codes used to the generate the heatmaps and scatterplots shown in the manuscript "Ectomycorrhizal effects on decomposition are highly dependent on fungal traits, climate, and litter properties: A model-based assessment"The majority of the model outputs did not have specific geographic information or temporal coverage because the analysis we conducted are mainly hypothetical model simulations. We only provided geographic description, coordinates and temporal coverage for those soil measurements which we used for model evaluations (included in the "Plot NH_GA simulation against measurements.zip").

54 ENVIRONMENTAL SCIENCES↗

Correlation between the gas-phase metallicity and ionization parameter in extragalactic H II regions

The variations of the metallicity and ionization parameter in H II regions are usually thought to be the dominant factors that produce the variations we see in the observed emission line spectra. There is an increasing amount of evidence that these two quantities are physically correlated, although the exact form of this correlation is debatable in the literature. Simulated emission line spectra from photoionized clouds provide important clues about the physical conditions of H II regions and are frequently used for deriving metallicities and ionization parameters. Through a systematic investigation on the assumptions and methodology used in applying photoionization models, we find that the derived correlation has a strong dependence on the choice of model parameters. On the one hand, models that give consistent predictions over multiple emission-line ratios yield a positive correlation between the metallicity and ionization parameter for the general population of H II regions or star-forming galaxies. On the other hand, models that are inconsistent with the data locus in high-dimensional line ratio space yield discrepant correlations when different subsets of line ratios are used in the derivation. The correlation between the metallicity and ionization parameter has a secondary dependence on the surface density of the star formation rate (SFR), with the higher SFR regions showing a higher ionization parameter but weaker correlations. The existence of the positive correlation contradicts the analytical wind-driven bubble model for H II regions. We explore assumptions in both dynamical models and photoionization models, and conclude that there is a potential bias associated with the geometry. However, this is still insufficient to explain the correlation. Mechanisms that suppress the dynamical influence of stellar winds in realistic H II regions might be the key to solving this puzzle, though more sophisticated combinations of dynamical models and photoionization models to test are required.

79 ASTRONOMY AND ASTROPHYSICS↗

Simplified spin dependence in dark matter direct detection

The interactions of dark matter with Standard Model particles can be systematically studied in the language of effective field theories. We investigate dark matter interactions with Standard Model particles, including spin-dependent interactions, for direct detection experiments and demonstrate that, although the scattering rate generally depends on multiple types of material response functions, certain linear combinations of these material response functions vanish if the initial and final electronic states share the same Hamiltonian. We also find that several other response functions vanish in parity-symmetric materials, making these systems as simple as isotropic detectors in some respects. Finally, we present the scattering rate for an anisotropic, possibly chiral detector, for generic dark matter-electron spin interactions. These relations reduce the number of independent response functions needed, thereby simplifying the computational complexity for a broad class of dark matter models. Our results provide a complete and efficient toolkit for analyzing electron recoil signals in diverse detector materials.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Impact of wildfires on regional ozone and PM 2.5 : Considering the light absorption of Brown carbon

The influence of wildfire brown carbon (BrC) in moderating the formation of O 3 and fine particulate matter (PM 2.5 ) due to its absorption of ultraviolet (UV) radiation is, for the first time, investigated using a modified community multiscale air quality (CMAQ) model. The wavelength-dependent imaginary refractive index of the organic carbon from wildfires, which are needed by the CMAQ model for its inline photolysis rate calculation, are generated experimentally from wood burning aerosols from a combustion chamber. On high emission days of the Bastrop County Complex fire in Texas in early September 2011, BrC UV absorption reduces the daytime average NO 2 photolysis rate by up to 38% in the core region of the wildfire plume compared to the base scenario that does not consider BrC absorption. Consequently, O 3 production in the smoke plume is suppressed compared to the base scenario due to reduced HO x radical concentrations. Further, in the core region, the predicted O 3 increase due to wildfire reaches as high as 47–123 ppb without considering BrC absorption, the predicted increase of O 3 is 5–15% lower when BrC absorption is considered. Similarly, considering the BrC UV absorption leads to approximately 1% (or ~2–3 µgm -3 ) lower estimation of the wildfire emissions’ impact on total PM 2.5 . This change is small because secondary aerosols, which are the components affected by BrC absorption, only account for a small fraction of the total PM 2.5 in wildfire impacted regions in this study. In addition, our study shows that assumptions about aerosol mixing state (core-shell vs. homogeneous) in the inline photolysis rate calculation would not significantly affect out assessment of the impact of wildfire BrC light absorption on O 3 and PM 2.5 .

54 ENVIRONMENTAL SCIENCES↗

Spectrally-dependent positive and negative photo-plastic effects in CdZnTe and CdZnTeSe

We studied the Vickers microhardness HV0.025 of CdZnTe and CdZnTeSe samples upon illumination with quasi-monochromatic light with wavelengths ranging between 750 and 1540 nm. We observed that light with a wavelength close to or above the bandgap of the samples results in a positive photo-plastic effect (hardening), and light in the spectral range well below the band gap can lead both to positive or negative effects depending on the parameters of deep levels in the bandgap. We performed measurements of spectrally dependent photoconductivity and compared them with HV0.025. Based on the analysis, we propose a complex qualitative model explaining the spectral dependence of HV0.025 by the charging and neutralization of deep levels by optical transitions between the valence and conduction bands and between the bands and deep levels in the band gap. Tuning of parameters of deep levels could be used to control the sign and magnitude of the photo-plastic effect.

36 MATERIALS SCIENCE↗

Sparsity dependence of Krylov state complexity in the SYK model

We study the Krylov state complexity of the Sachdev-Ye-Kitaev (SYK) model for 𝑁 ≤28 Majorana fermions with 𝑞-body fermion interaction with 𝑞 =4, 6, 8 for a range of sparse parameter 𝑘 that controls the number of remaining terms in the original SYK model after sparsification. The critical value of 𝑘 below which the model ceases to be holographic, denoted 𝑘 𝑐 , has been subject of several recent investigations. Using Krylov complexity as a probe, we find that the peak value of complexity does not change as we increase 𝑘 beyond 𝑘 ≥ 𝑘 min at large temperatures. We argue that this behavior is related to the change in the holographic nature of the Hamiltonian in the sparse SYK-type models such that the model is holographic for all 𝑘 ≥ 𝑘 min ≈ 𝑘 𝑐 . Our results provide a novel way to determine 𝑘 𝑐 in SYK-type models.

gauge-gravity dualitites↗

Evolutionary NAS in Light of Model Stability for Accurate Continual Learning

Continual learning, the capability to learn new knowledge from streaming data without forgetting the previous knowledge, is a critical requirement for dynamic learning systems, especially for emerging edge devices such as self-driving cars and drones. However, continual learning is still facing the catastrophic forgetting problem. Previous work illustrate that model performance on continual learning is not only related to the learning algorithms but also strongly dependent on the inherited model, i.e., the model where continual learning starts. The better stability of the inherited model, the less catastrophic forgetting and thus, the inherited model should be elaborately selected. Inspired by this finding, we develop an evolutionary neural architecture search (ENAS) algorithm that emphasizes the Stability of the inherited model, namely ENAS-S. ENAS-S aims to find optimal architectures for accurate continual learning on edge devices. On CIFAR-10 and CIFAR-100, we present that ENAS-S achieves competitive architectures with lower catastrophic forgetting and smaller model size when learning from a data stream, as compared with handcrafted DNNs.

Du, Xiaocong↗

A workflow for automatic generation and efficient refinement of individual pressure-dependent networks

Manual chemical kinetic model construction often requires modelists to at least implicitly guess all possible decay pathways and their relative fluxes for all important species. We present a workflow and associated tool for automatic generation and re finement of pressure-dependent networks that should enable modelists to efficiently and comprehensively identify decay pathways and estimate respective parameters for chemical species. This tool combines the capabilities of the Reaction Mechanism Generator (RMG) software to generate possible reaction paths, determine thermochemistry, approximate rate coefficients and estimate frequencies with the capabilities of the Arkane software to use quantum chemical parameters to compute thermochemistry and pressure dependent rates. A flux-based algorithm is used to decide which isomers to add to the network. Isomers added to the network are reacted to form other channels. When enough of the flux is accounted for by the isomers and bimolecular product channels, the generation process terminates to yield a comprehensive network. Network sensitivity analysis is applied to the network to identify the most important wells and barriers that can be refined using quantum chemistry calculations. Iterative refinement can then be used to achieve accuracy. The comprehensive network can be easily reduced using energy and flux-based algorithms to the most important channels and isomers in the network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spread Spectrum Time Domain Reflectometry (SSTDR) Digital Twin Simulation of Photovoltaic Systems for Fault Detection and Location

Utilizing spread spectrum time domain reflectometry (SSTDR) to detect, locate, and characterize faults in photovoltaic (PV) systems is examined in this paper. We present a method to obtain the model parameters that are needed to produce digital twin SSTDR responses for PV systems. The digital twin SSTDR responses could be used to predict faults within the PV systems. Here, the model parameters are the reflection and transmission coefficients at each impedance discontinuity in the PV system along with the propagation coefficients across each PV cable segment. We obtain model parameter by applying inverse modeling techniques to experimental SSTDR data associated with PV systems. Our model parameters can be used in any digital twin simulation method for modeling reflectometry in frequency-dependent and complex loads. For validation, we used the model parameters in a graph network simulation engine and adapted it to be used for SSTDR digital twin simulations in PV systems. We produced simulations for 0 to 10 PV modules connected in series. We also simulated SSTDR responses for open circuit disconnections in a PV setup containing 10 PV modules in series. Results show that all but one simulated disconnect locations match experimental disconnection locations of the same setup with an error of less than 5%.

14 SOLAR ENERGY↗

Modeling Needs for Power Semiconductor Devices and Power Electronics Systems: Preprint

As energy systems move towards wide-spread electrification, penetration of power semiconductor devices and power electronics continues to grow at a rapid pace. This paper attempts to review the state of the art in power semiconductor modeling and the existing state of modeling in general when in it comes to design of devices and systems. The modeling needs for semiconductor devices vary depending on the end goal, and the level of abstraction needed towards model formulation also changes with the size and complexity of system to be designed. Finally, key aspects of thermal design and modeling for power electronics are also discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Enhanced accuracy through ensembling of randomly initialized auto-regressive models for dynamical systems

Computational mechanics simulations using traditional finite element methods (FEM) require prohibitively expensive computational resources for real-time engineering applications, design optimization, and digital twin implementations. While machine learning (ML) surrogate models offer significant computational speedups, autoregressive ML models for time-dependent mechanical systems suffer from error accumulation that compromises long-term prediction reliability - a critical concern for engineering applications where accuracy over extended time horizons is essential for safety and performance assessments. Here, we propose a deep ensemble framework specifically designed to address this challenge in computational mechanics applications, where multiple ML surrogate models with random weight initializations are trained in parallel and their predictions aggregated during inference. This approach leverages statistical diversity to maximize information gain from a fixed set of training data and to mitigate error propagation, while maintaining the computational efficiency that makes ML surrogates attractive for engineering practice. We validate the framework on three representative problems spanning critical areas of computational mechanics: stress field evolution in heterogeneous microstructures under complex loading (relevant to advanced materials design and composite analysis), planetary-scale shallow water dynamics (applicable to environmental and geotechnical engineering), and Gray-Scott reaction-diffusion systems (relevant to mass transport and chemical process engineering). Across all test cases, the ensemble approach demonstrates consistent error reduction of 15-33% compared to individual models. The codes for this work are available on GitHub (https://github.com/Graham-Brady-Research-Group/AutoregressiveEnsemble_SpatioTemporal_Evolution).

autoregressive prediction↗

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗

A multi-model assessment of the early last deglaciation (PMIP4 LDv1): a meltwater perspective

Abstract. The last deglaciation (∼20–11 ka BP) is a period of a major, long-term climate transition from a glacial to interglacial state that features multiple centennial- to decadal-scale abrupt climate variations whose root cause is still not fully understood. To better understand this time period, the Paleoclimate Modelling Intercomparison Project (PMIP) has provided a framework for an internationally coordinated endeavour in simulating the last deglaciation whilst encompassing a broad range of models. Here, we present a multi-model intercomparison of 17 transient simulations of the early part of the last deglaciation (∼20–15 ka BP) from nine different climate models spanning a range of model complexities and uncertain boundary conditions and forcings. The numerous simulations available provide the opportunity to better understand the chain of events and mechanisms of climate changes between 20 and 15 ka BP and our collective ability to simulate them. We conclude that the amount of freshwater forcing and whether it follows the ice sheet reconstruction or induces an inferred Atlantic meridional overturning circulation (AMOC) history, heavily impacts the deglacial climate evolution for each simulation rather than differences in the model physics. The course of the deglaciation is consistent between simulations except when the freshwater forcing is above 0.1 Sv – at least 70 % of the simulations agree that there is warming by 15 ka BP in most places excluding the location of meltwater input. For simulations with freshwater forcings that exceed 0.1 Sv from 18 ka BP, warming is delayed in the North Atlantic and surface air temperature correlations with AMOC strength are much higher. However, we find that the state of the AMOC coming out of the Last Glacial Maximum (LGM) also plays a key role in the AMOC sensitivity to model forcings. In addition, we show that the response of each model to the chosen meltwater scenario depends largely on the sensitivity of the model to the freshwater forcing and other aspects of the experimental design (e.g. CO2 forcing or ice sheet reconstruction). The results provide insight into the ability of our models to simulate the first part of the deglaciation and how choices between uncertain boundary conditions and forcings, with a focus on freshwater fluxes, can impact model outputs. We can use these findings as helpful insight in the design of future simulations of this time period.

Snoll, Brooke (ORCID:0000000230914394)↗

Status Update: Deposition Modeling For SNF Canister CISCC

This report fulfills the M3 milestones M3SF-21PN010207025 & M3SF-20PN0102070412. During fiscal year (FY) 2020, Pacific Northwest National Laboratory (PNNL) worked to further develop the FY 2019 deposition and particle tracking models. This status report outlines these efforts and presents the progress made so far. Model development work is ongoing and is planned to continue in FY 2021. The FY 2020 model development included work on: Wind Effects. Model development and sensitivity studies, investigated how wind direction and speed affect deposition; Brownian Motion. Implementing Brownian Motion into existing models; Particle Size Variability. Depending on the particle composition, the diameter of the particle may vary with changes in relative humidity. Models were developed to analyze this; Multiphase and Fluid Film Modeling. Investigating canister surface wetting and drying, and how this effects overall deposition. Models were developed to analyze this; Difusophoresis. Performing initial work to implement diffusiophoresis into the existing models; Turbophoresis. Performing initial work to implement turbophoresis into the existing models. Much of this work will continue into FY21. The authors present initial results and discus current and future work.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

One-to-one aeroservoelastic validation of operational loads and performance of a 2.8 MW wind turbine model in OpenFAST

Abstract. This article presents a validation study of the popular aeroservoelastic code suite OpenFAST leveraging weeks of measurements obtained during normal operation of a 2.8 MW land-based wind turbine. Measured wind conditions were used to generate one-to-one turbulent flow fields (i.e., comparing simulation to measurement in 10 min increments, or bins) through unconstrained and constrained assimilation methods using the kinematic turbulence generators TurbSim and PyConTurb. A total of 253 bins of 10 min of normal turbine operation were selected for analysis, and a statistical comparison in terms of performance and loads is presented. We show that successful validation of the model was not strongly dependent on the type of inflow assimilation method used for mean quantities of interest, which had median modeling errors per wind-speed interval generally within 5 %–10 % of the measurement. The type of inflow assimilation method did have a larger effect on the fatigue predictions for blade-root flapwise and tower-base fore–aft quantities, which surprisingly saw larger errors from the assumed higher-fidelity assimilation methods. Avenues for further work are discussed and include possible improvements to the aerodynamic, structural, and controller modeling that may offer insight on the origin of the up to ∼ 40 % median overprediction of fatigue for these quantities.

17 WIND ENERGY↗