Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Best practices for photovoltaic performance loss rate calculations

Abstract The performance loss rate (PLR) is a vital parameter for the time-dependent assessment of photovoltaic (PV) system performance and health state. Although this metric can be calculated in a relatively straightforward manner, it is challenging to achieve accurate and reproducible results with low uncertainty. Furthermore, the temporal evolution of PV system performance is usually nonlinear, but in many cases a linear evaluation is preferred as it simplifies the assessment and it is easier to evaluate. As such, the search for a robust and reproducible calculation methodology providing reliable linear PLR values across different types of systems and conditions has been the focus of many research activities in recent years. In this paper, the determination of PV system PLR using different pipelines and approaches is critically evaluated and recommendations for best practices are given. As nonlinear PLR assessments are fairly new, there is no consent on how to calculate reliable values. Several promising nonlinear approaches have been developed recently and are presented as tools to evaluate the PV system performance in great detail. Furthermore, challenges are discussed with respect to the PLR calculation but also opportunities for differentiating individual performance losses from a generic PLR value having the potential of enabling actionable insights for maintenance.

14 SOLAR ENERGY↗
DOE OSTI2023

Recipes for when physics fails: recovering robust learning of physics informed neural networks

Abstract Physics-informed neural networks (PINNs) have been shown to be effective in solving partial differential equations by capturing the physics induced constraints as a part of the training loss function. This paper shows that a PINN can be sensitive to errors in training data and overfit itself in dynamically propagating these errors over the domain of the solution of the PDE. It also shows how physical regularizations based on continuity criteria and conservation laws fail to address this issue and rather introduce problems of their own causing the deep network to converge to a physics-obeying local minimum instead of the global minimum. We introduce Gaussian process (GP) based smoothing that recovers the performance of a PINN and promises a robust architecture against noise/errors in measurements. Additionally, we illustrate an inexpensive method of quantifying the evolution of uncertainty based on the variance estimation of GPs on boundary data. Robust PINN performance is also shown to be achievable by choice of sparse sets of inducing points based on sparsely induced GPs. We demonstrate the performance of our proposed methods and compare the results from existing benchmark models in literature for time-dependent Schrödinger and Burgers’ equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Differences in substrate use linked to divergent carbon flow during litter decomposition

ABSTRACT Discovering widespread microbial processes that create variation in soil carbon (C) cycling within ecosystems may improve soil C modeling. Toward this end, we screened 206 soil communities decomposing plant litter in a common garden microcosm environment and examined features linked to divergent patterns of C flow. C flow was measured as carbon dioxide (CO2) and dissolved organic carbon (DOC) from 44-days of litter decomposition. Two large groups of microbial communities representing ‘high’ and ‘low’ DOC phenotypes from original soil and 44-day microcosm samples were down-selected for fungal and bacterial profiling. Metatranscriptomes were also sequenced from a smaller subset of communities in each group. The two groups exhibited differences in average rate of CO2 production, demonstrating that the divergent patterns of C flow arose from innate functional constraints on C metabolism, not a time-dependent artefact. To infer functional constraints, we identified features – traits at the organism, pathway or gene level – linked to the high and low DOC phenotypes using RNA-Seq approaches and machine learning approaches. Substrate use differed across the high and low DOC phenotypes. Additional features suggested that divergent patterns of C flow may be driven in part by differences in organism interactions that affect DOC abundance directly or indirectly by controlling community structure.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Characterizing the uppermost 100 m structure of the San Jacinto fault zone southeast of Anza, California, through joint analysis of geological, topographic, seismic and resistivity data

We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (<100 m) fault zone at this structurally complex site. Geological mapping at the surface shows the CF has three main subparallel strands within a <100 m zone with varying degrees of rock damage. These strands intersect units of banded gneiss and tonalite, and various sedimentary units. Near the surface, the weathered but more intact tonalite and gneiss to the southwest have relatively high VP. The low-lying flat sedimentary basins around the two southwestern-most CF strands and elevated damaged gneiss to the northeast have lowest VP <500 m s–1. The high relief of the northeast gneiss unit may in part be explained by its extensive damage and inferred increased relative rock uplift. Resistivity imaging shows the unconsolidated dry basin sediments (maximum >1300 Ohm.m) contrasted against the compacted fine-grained (potentially wet) materials within the CF core and the Bautista Formation (minimum <40 Ohm.m), which is slightly elevated above the flat basins. The inverse relationship between VP (increases) and resistivity (decreases) in the uppermost ~15 m can be characterized as log–log linear with slopes of –2.6 to –4. At depths >30 m, the velocity heterogeneity near the surface merges into larger-scale structures that are generally slower on the northeast side of the CF core compared to the southwest side (as much as ~40 per cent reduction in average VP). A previous study revealed a 20–37 per cent variability in peak ground velocities across the SGB site from local earthquakes. The upper end of that range is associated with the near-surface unconsolidated sedimentary basins and northeast damaged gneiss unit. Preliminary analysis of time-dependent topography mostly shows effects of changing vegetation and anthropogenic activity.

Geochemistry & Geophysics↗
DOE OSTI2023

En echelon faults reactivated by wastewater disposal near Musreau Lake, Alberta

We use machine-learning and cross-correlation techniques to enhance earthquake detectability by two magnitude units for the earthquake sequence near Musreau Lake, Alberta, which is induced by wastewater disposal. This deep catalogue reveals a series of en echelon ~N–S oriented strike-slip faults that are favourably oriented for reactivation. These faults require only ~0.6 MPa overpressure for triggering to occur. Earthquake activity occurs in bursts, or episodes; episodes restricted to the largest fault tend to have earthquakes starting near the southern end (distant from injectors) and progressing northwards (towards the injectors). While most events are concentrated along these ~N–S oriented faults, we also delineate smaller faults. Together, these findings suggest pore pressure as the triggering mechanism, where a time-dependent increase in pore pressure likely caused these faults to progressively reawaken. Analysis of the ‘next record-breaking event’, a statistical model that forecasts the sequencing of earthquake magnitudes, suggests that the next largest event would be M L ~4.3. The seismically illuminated length of the largest fault indicates potential magnitudes as large as M w 5.3.

58 GEOSCIENCES↗
DOE OSTI2022

A Process-based Model with Temperature, Water, and Lab-derived Data Improves Predictions of Daily Culex pipiens/restuans Mosquito Density

While the number of human cases of mosquito-borne diseases has increased in North America in the last decade, accurate modeling of mosquito population density has remained a challenge. Longitudinal mosquito trap data over the many years needed for model calibration, and validation is relatively rare. In particular, capturing the relative changes in mosquito abundance across seasons is necessary for predicting the risk of disease spread as it varies from year to year. We developed a discrete, semi-stochastic, mechanistic process-based mosquito population model that captures life-cycle egg, larva, pupa, adult stages, and diapause for Culex pipiens (Diptera, Culicidae) and Culex restuans (Diptera, Culicidae) mosquito populations. This model combines known models for development and survival into a fully connected age-structured model that can reproduce mosquito population dynamics. Mosquito development through these stages is a function of time, temperature, daylight hours, and aquatic habitat availability. The time-dependent parameters are informed by both laboratory studies and mosquito trap data from the Greater Toronto Area. The model incorporates city-wide water-body gauge and precipitation data as a proxy for aquatic habitat. This approach accounts for the nonlinear interaction of temperature and aquatic habitat variability on the mosquito life stages. We demonstrate that the full model predicts the yearly variations in mosquito populations better than a statistical model using the same data sources. This improvement in modeling mosquito abundance can help guide interventions for reducing mosquito abundance in mitigating mosquito-borne diseases like West Nile virus.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The effect of saturated thermal conduction on clouds in a hot plasma

We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High-energy neutrino emission from magnetized jets of rapidly rotating protomagnetars

ABSTRACT Relativistic jets originating from protomagnetar central engines can lead to long duration gamma-ray bursts (GRBs) and are considered potential sources of ultra-high-energy cosmic rays and secondary neutrinos. We explore the propagation of such jets through a broad range of progenitors, from stars which have shed their envelopes to supergiants which have not. We use a semi-analytical spin-down model for the strongly magnetized and rapidly rotating protoneutron star (PNS) to investigate the role of central engine properties such as the surface dipole field strength, initial rotation period, and jet opening angle on the interactions and dynamical evolution of the jet-cocoon system. With this model, we determine the properties of the relativistic jet, the mildly relativistic cocoon, and the collimation shock in terms of system parameters such as the time-dependent jet luminosity, injection angle, and density profile of the stellar medium. We also analyse the criteria for a successful jet breakout, the maximum energy that can be deposited into the cocoon by the relativistic jet, and structural stability of the magnetized outflow relative to local instabilities. Lastly, we compute the high-energy neutrino emission as these magnetized outflows burrow through their progenitors. Precursor neutrinos from successful GRB jets are unlikely to be detected by IceCube, which is consistent with the results of previous works. On the other hand, we find that high-energy neutrinos may be produced for extended progenitors like blue and red supergiants, and we estimate the detectability of neutrinos with next generation detectors such as IceCube-Gen2.

Astronomy & Astrophysics↗
DOE OSTI2025

Quantifying the impact of the Si/O interface in CCSN explosions using the Force Explosion Condition

ABSTRACT The explosion mechanism of a core-collapse supernova (SN) is a complex interplay between neutrino heating and cooling (including the effects of neutrino-driven convection), the gravitational potential, and the ram pressure of the infalling material. To analyse the post-bounce phase of an SN, one can use the generalized force explosion condition (FEC+), which succinctly formalizes the interplay among these four phenomena in an analytical condition, consistent with realistic simulations. In this paper, we use the FEC+ to study the post-bounce phase of 341 spherically symmetric simulations, where convection is included through a time-dependent mixing-length approach. We find that the accretion of the Si/O interface through the expanding shock can significantly change the outcome of the SN by driving the FEC+ above the explosion threshold. We systematically explore this by (i) artificially smoothing the pre-SN density profile, and (ii) artificially varying the mixing length. In both cases, we find that large-enough density contrasts at the Si/O interface lead to successful shock revival only if the FEC+ is already close to the explosion threshold. Furthermore, we find that the accretion of the Si/O interface has a substantial effect on the critical condition for SN explosions, contributing between 5 and 15 per cent, depending on how pronounced the density contrast at the interface is. Earlier studies showed that convection affects the critical condition by 25–30 per cent, which demonstrates that the accretion of the Si/O interface through the shock can play a nearly comparable role in influencing shock dynamics.

Boccioli, Luca (ORCID:000000024819310X)↗
DOE OSTI2026

Dynamics and observational signatures of core-collapse supernovae with central engines: hydrodynamics simulations with Monte Carlo post-processing

A long-lived central engine embedded in expanding supernova ejecta can alter the dynamics and observational signatures of the event, producing an unusually luminous, energetic, and/or rapidly evolving transient. We use 2D hydrodynamics simulations to study the effect of a central energy source, varying the amount, rate, and isotropy of the energy deposition. We post-process the results with a time-dependent Monte Carlo radiation transport code to extract observational signatures. The engine excavates a bubble at the centre of the ejecta, which becomes Rayleigh–Taylor unstable. Sufficiently powerful engines are able to break through the edge of the bubble and accelerate, shred, and compositionally mix the entire ejecta. The breakout of the engine-driven wind occurs at distinct rupture points, and the outflowing high-velocity gas may eventually give rise to radio emission. The dynamical impact of the engine leads to faster rising optical light curves, with photon escape facilitated by the faster expansion of the ejecta and the opening of low-density channels. For models with strong engines, the spectra are initially hot and featureless, but later evolve to resemble those of broad-line Ic supernovae. Under certain conditions, line emission from ionized, low-velocity material near the centre of the ejecta may be able to escape and produce narrow emission similar to that seen in interacting supernovae. We discuss how variability in the engine energy reservoir and injection rate could give rise to a heterogeneous set of events spanning multiple observational classes, including the fast blue optical transients, broad-line Ic supernovae, and superluminous supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
DOE OSTI2023

Dynamic and transient processes in warm dense matter

In this paper, we discuss some of the key challenges in the study of time-dependent processes and non-equilibrium behaviour in warm dense matter. We outline some of the basic physics concepts that have underpinned the definition of warm dense matter as a subject area in its own right and then cover, in a selective, non-comprehensive manner, some of the current challenges, pointing along the way to topics covered by the papers presented in this volume.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Direct structural retrieval from gas-phase ultrafast diffraction data using a genetic algorithm

Ultrafast scattering techniques such as ultrafast electron diffraction and ultrafast x-ray diffraction have been utilized to elucidate the structural dynamics, reaction intermediates, and final products in molecular reactions following photoexcitation. The time-dependent structures are typically not directly retrieved from the experimental data, but they rely on comparison with calculations. The genetic algorithm (GA), a global optimization strategy, can be used to retrieve the molecular structures directly from diffraction patterns without any theoretical input. However, the robustness of the GA with respect to real experimental conditions such as a limited momentum transfer range, noise, and artifacts has not been studied in detail. In this work, we characterize the performance of the GA with simulated data that mimic realistic experimental conditions. We have developed and implemented a variant of the GA specific to diffraction measurements which performs better in the presence of imperfect data compared to the standard implementation of the GA. We demonstrate this method with both synthetic data and experimental ultrafast electron diffraction data on the UV-induced photodissociation of trifluoroiodomethane (C⁢F 3⁡ I) molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Accessing the gluon momentum fraction of nucleons through the gradient flow

We calculate the gluon momentum fraction of the nucleon using lattice QCD, with a nonperturbative renormalization technique based on the gradient flow. The gluon momentum fraction is determined on a single Wilson-clover ensemble using 𝑁 𝑓 =2 +1 flavors with pion mass 358 MeV and lattice spacing 0.094 fm. We employ the variational method to reduce excited-state contamination and apply the distillation framework to ensure a large operator basis. To reduce systematic uncertainties, we apply Bayesian model averaging to all fit procedures. We apply matching coefficients to the flow-time dependent lattice results to recover the gluon momentum fraction in the $\overline{MS}$-scheme at 2 GeV. Our final result is ⟨𝑥⟩ 𝑔 ⁢(𝜇 =2 GeV) =0.482⁢(35), where we quote only statistical uncertainties.

Lattice QCD↗
DOE OSTI2026

Thresholded quantum sensing with a frustrated Kitaev trimer

We investigate the response of a Ramsey interferometric quantum sensor based on a frustrated, three-spin system (a Kitaev trimer) to a classical time-dependent field (signal). The system eigenspectrum is symmetric about a critical point, |𝒃| = 0, with four of the spectral components varying approximately linearly with the magnetic field and four exhibiting a nonlinear dependence. Under the adiabatic approximation and for appropriate initial states, we show that the sensor's response to a zero-mean signal is such that below a threshold, |𝒃| <⁢ 𝑏 th , the sensor does not respond to the signal, whereas above the threshold, the sensor acts as a detector that the signal has occurred. This thresholded response is approximately omnidirectional. Moreover, when deployed in an entangled multisensor configuration, the sensor achieves sensitivity at the Heisenberg limit. Such detectors could be useful both as stand-alone units for signal detection above a noise threshold and in two- or three-dimensional arrays, analogous to a quantum bubble chamber, for applications such as particle track detection and long-baseline telescopy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Finite-volume formalism for physical processes with an electroweak loop integral

This study investigates finite-volume effects in physical processes that involve the combination of long-range hadronic matrix elements with electroweak loop integrals. We adopt the approach of implementing the electroweak part as the infinite-volume version, which is denoted as the EW ∞ method in this work. A general approach is established for correcting finite-volume effects in cases where the hadronic intermediate states are dominated by either a single particle or two particles. For the single-particle case, this work derives the infinite volume reconstruction method from a new perspective. For the two-particle case, we provide the correction formulas for power-law finite-volume effects and unphysical terms with exponentially divergent time dependence. The finite-volume formalism developed in this study has broad applications, including the QED corrections in various processes and the two-photon exchange contribution in 𝐾 𝐿 → 𝜇 + ⁢𝜇 − or 𝜂 → 𝜇 + ⁢𝜇 − decays.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗