SEARCH · Engineering Papers
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.
In this paper, an analytical method for comparing the performance of an ideal forecasting algorithm to the best statistical methods in terms of average data return is derived.
We describe the near real-time transient-source discovery engine for the intermediate Palomar Transient Factory (iPTF), currently in operations at the Infrared Processing and Analysis Center (IPAC), Caltech. We coin this system the IPAC/iPTF Discovery Engine (or IDE). We review the algorithms used for PSF-matching, image subtraction, detection, photometry, and machine-learned (ML) vetting of extracted transient candidates. We also review the performance of our ML classifier. For a limiting signal-to-noise ratio of 4 in relatively unconfused regions, bogus candidates from processing artifacts and imperfect image subtractions outnumber real transients by approximately equal to 10:1. This can be considerably higher for image data with inaccurate astrometric and/or PSF-matching solutions. Despite this occasionally high contamination rate, the ML classifier is able to identify real transients with an efficiency (or completeness) of approximately equal to 97% for a maximum tolerable false-positive rate of 1% when classifying raw candidates. All subtraction-image metrics, source features, ML probability-based real-bogus scores, contextual metadata from other surveys, and possible associations with known Solar System objects are stored in a relational database for retrieval by the various science working groups. We review our efforts in mitigating false-positives and our experience in optimizing the overall system in response to the multitude of science projects underway with iPTF.
The acoustic analogy introduced by Lighthill to study jet noise is now over 50 years old. In the present paper, Lighthill s Acoustic Analogy is revisited together with a brief evaluation of the state-of-the-art of the subject and an exploration of the possibility of further improvements in jet noise prediction from analytical methods, computational fluid dynamics (CFD) predictions, and measurement techniques. Experimental Particle Image Velocimetry (PIV) data is used both to evaluate turbulent statistics from Reynolds-averaged Navier-Stokes (RANS) CFD and to propose correlation models for the Lighthill stress tensor. The NASA Langley Jet3D code is used to study the effect of these models on jet noise prediction. From the analytical investigation, a retarded time correction is shown that improves, by approximately 8 dB, the over-prediction of aft-arc jet noise by Jet3D. In experimental investigation, the PIV data agree well with the CFD mean flow predictions, with room for improvement in Reynolds stress predictions. Initial modifications, suggested by the PIV data, to the form of the Jet3D correlation model showed no noticeable improvements in jet noise prediction.
In striking contrast to equilibrium systems, inertia can profoundly alter the structure of active systems. Here, we demonstrate that driven systems can exhibit effective equilibrium-like states with increasing particle inertia, despite rigorously violating the fluctuation–dissipation theorem. Increasing inertia progressively eliminates motility-induced phase separation and restores equilibrium crystallization for active Brownian spheres. This effect appears to be general for a wide class of active systems, including those driven by deterministic time-dependent external fields, whose nonequilibrium patterns ultimately disappear with increasing inertia. The path to this effective equilibrium limit can be complex, with finite inertia sometimes acting to accentuate nonequilibrium transitions. The restoration of near equilibrium statistics can be understood through the conversion of active momentum sources to passive-like stresses. Unlike truly equilibrium systems, the effective temperature is now density dependent, the only remnant of the nonequilibrium dynamics. Further, this density-dependent temperature can in principle introduce departures from equilibrium expectations, particularly in response to strong gradients. Our results provide additional insight into the effective temperature ansatz while revealing a mechanism to tune nonequilibrium phase transitions.
Determining wildland fire impacts on streamflow can be problematic as the hydrology in burned watersheds is influenced by post-fire weather conditions. Here, this study presents a quantile-based analytical framework for assessing fire impacts on low and peak daily flow magnitudes, while accounting for post-fire weather influences. This framework entails (a) the bootstrap method to compute the relative change in the post-fire annual flow and weather statistics, (b) double mass analysis to detect if post-fire baseflow and quick-flow yield ratios are significantly altered, and (c) a quantile regression method to parse fire effects on flow at a specific quantile. We illustrate the applicability of this analytical framework using 44 western US streams with at least 5% of their watershed area burned. Results indicate that large, high-severity burns in upland watersheds can raise the streamflow magnitude at the 0. 05 th and 0. 95 th quantiles for at least the five post-fire years. Quantile regression results show that the median fire-related increase in flow for the five post-fire years can be up to 5,000% (Standard Error; S.E. < 2%) at the 0. 05 th E quantile and 161% (S.E. < 10%) at the 0. 95 th quantile. The fire-related increase in flow was often pronounced at the 0. 05 th quantile for streams in the Pacific Northwest and California regions. The difference in fire effects on flow (at both quantiles) across streams was related to post-fire weather, pyrology, physiography, and land cover. The proposed analytical framework can be useful for detecting and quantifying fire effects on the low and peak stream flows in burned watersheds without overlapping disturbances.
The modeling of a spatially incoherent laser beam remains a central problem of the parametric instabilities in the context of inertial confinement fusion. This letter gives a simplified and comprehensive overview of the recent analytical developments regarding the modeling of these laser beams and a comparison with a dedicated experiment. Our model accounts for the first time for the statistical standard deviation of the gain and accurately captures the entanglement between wave mixing processes and the speckle correlations thus resolving the longstanding contradictions between the random phase approximation and the model of independent speckles. It is successfully compared to a recent laser beam spray experiment and the associated paraxial simulations, demonstrating that backscattering predictions require accounting for the beam spray. Furthermore, our framework thus provides a way to evaluate and guide the analysis of parametric instabilities in high laser energy experiments.
The purpose of this document is to summarize environmental monitoring data collected at the Remote-Handled Low-Level Waste (RHLLW) Disposal Facility during the first four years of facility operations (FY 2019 through FY 2022). This summary provides a “baseline” condition for the facility, which can be used to distinguish contaminant releases from the RHLLW Disposal Facility from pre existing contamination as well as potential future releases from other sources (i.e., upgradient aquifer sources), as specified in the facility monitoring plan. Sufficient data has been collected over the first four years of operation of the RHLLW Disposal Facility to establish baseline conditions for future compliance monitoring of aquifer wells and performance monitoring of vadose zone lysimeters. Except for seven elevated tritium measurements from lysimeter HFEF-South believed to have been impacted by waste disposals, all data was deemed appropriate for establishing baseline concentrations. Monitoring of the aquifer detected indicator analytes gross alpha and gross beta and target analytes tritium and C-14. There were a few C-14 detections prior to an increase in the required detection level (RDL) after the first round of sampling. I-129 and Tc-99 were not detected above RDLs in aquifer samples. Lysimeter samples detected indicator analytes gross alpha and gross beta and target analyte tritium. Target analytes C-14, I-129, and Tc-99 were not detected above RDLs in lysimeter samples. Except for tritium, baseline concentrations of indicator and target analytes measured during the first four years of RHLLW facility operations are expected to represent conditions through the expected 20 year operating period of the facility. Elevated tritium in the aquifer, a result of past discharges of tritium at upgradient facilities, showed a decline over the first four years of operations, and levels are expected to continue to decline with time as a result of dilution and decay. Current tritium levels in lysimeter samples may be elevated due to tritium in the water applied during construction and infiltration testing of the facility. As a result, tritium concentrations in lysimeters may decline with time as a result of dilution and decay. Statistical measures (i.e., mean, standard deviation, and 99% upper confidence level [UCL]) were calculated for all detected analytes using the four years of concentration data. These measures will be used to evaluate future monitoring results and to demonstrate the facility is performing, or not performing, as established in the facility performance assessment (PA). Future monitoring results will also be combined with the baseline data in this report to further establish temporal trends and natural variability in the measurements. Based on the data collected during the baseline period, it is recommended the gross alpha action level for performance monitoring increase from 10 to 20 pCi/L due to several baseline measurements exceeding the preliminary action level of 10. An action level of 20 pCi/L is slightly greater than the 99% UCL, is protective of the aquifer, and would reduce unnecessary sampling. It is also recommended that tritium be added to the lysimeter analyte list for post-baseline period monitoring. Tritium, while not a dose concern, is a good tracer that can provide valuable information on water flow in and around the RHLLW Disposal Facility. Recommendations for a tritium action level are provided if considered necessary.
Abstract The local structure of a transition metal (TM) ion is a function of cation elements and valence states. More than that, in this work, by employing a trove of first‐principles data of TM oxides, the local structures of TM cations are statistically analyzed to extract detailed information about cation site preference, bond length, site structural distortion, and cation magnetization. It is found that cation radius alone poorly describes the local structure of a transition metal oxide, while the statistics of coordination number as well as the TMO bond length distribution, especially that of the 3d TMs, can provide comprehensive knowledge for understanding the behavior of TM elements. Based on these statistics, the interplay of site distortion due to the Jahn–Teller effect, cation site similarity, and a new set of ionic radii are all obtained to chart the “persona” of transition metal ions in solids.
This work presents a generalized approach for analytical method optimization that branches the gap between techniques historically employed and accurate modern optimization techniques suitable for various applications. The novelty of the described strategy is the utilization of multivariate, multiobjective optimization with Karush-Kuhn-Tucker conditions to bound the optimization space to solutions within the physical limitations of instrumentation. Briefly, the basic steps outlined in this paper are to (1) determine the objective(s) that should be maximized or minimized based on the goals of the analytical application, (2) conduct a screening experiment, (3) perform ANOVA to determine the parameters which have a statistically significant effect on the objective, (4) conduct an experiment (e.g., Box-Behnken design) to collect data for fitting the objective equation, and (5) determine the physical constraints of the parameters and solve the Lagrangian to determine the optimal method parameters. A broad approach to optimization target selection allows for robust method tuning to develop improved data sets amenable for chemometrics and machine learning algorithm development. Gas chromatography-mass spectrometry was selected as a use case due to its broad use across scientific fields and time-consuming method development involving numerous parameters. In conclusion, this strategy can reduce the cost of research, improve data quality, and enable the rapid development of new analytical technique.
Protective vaccines and reliable diagnostics are essential tools for controlling viral diseases. However, the efficacy of these tools can be diminished by mutations in viral genomes. The delay between the emergence of new viral strains and the redesign of vaccines and diagnostics allows for continued viral transmission. Is it possible to address this challenge by computationally predicting viral genome sequence evolution? Can we “future-proof” vaccines and diagnostics by targeting both current and anticipated future sequence variants? While predicting viral evolution is still an unsolved, “grand challenge” problem in biology, the large, and rapidly growing, number of SARS-CoV-2 genome sequences provide an opportunity to quantify the ability of machine learning to predict viral genome sequence evolution. Towards this end, we have developed a simple computational model for predicting viral evolution at the level of individual nucleotides. The key metric for quantifying the per-base, prediction accuracy for viral evolution is the Mann-Whitney U statistic (or, equivalently, the area under the receiver operator curve). Since the Mann-Whitney U statistic is not a differentiable function, existing deep leaning packages (like Pytorch and Keras/TensorFlow) are not useful, as they require that the accuracy metric/objective function be analytically differentiable with respect to the model parameters. To overcome this challenge, we have implemented custom software, “FutureTense”, that can train a machine learning model by maximizing the non-differentiable Mann-Whitney U statistic. This software trains a machine learning model by exploring along the direction of the discrete gradient of the Mann-Whitney U statistic in the model parameter space. Parallel computing and genome sequence-specific optimizations are used to accelerate model training. The resulting machine learning model learns the observed high C->U mutation rates in the SARS-CoV-2 genome (which are potentially induced by host defenses) and provides prediction accuracies that are significantly better than one would expect from random chance. While predicting viral evolution is still quite far from a solved problem, the surprising performance of this simple model gives hope that the accuracy of predicting viral genome evolution can be further increased by more sophisticated approaches.
We study the fluctuations of subsystem entropies in closed quantum many-body systems after thermalization. Using a combination of analytics and numerics for both random quantum circuits and Hamiltonian dynamics, we find that the statistics of such entropy fluctuations is drastically different than in the classical setting. For instance, shortly after a system thermalizes, the probability of entropy fluctuations for a subregion is suppressed in the dimension of the Hilbert space of the complementary subregion. This suppression becomes increasingly stringent as a function of time, ultimately depending on the exponential of the Hilbert space dimension, until extremely late times when the amount of suppression saturates. We also use our results to estimate the total number of rare fluctuations at large timescales. We find that the “Boltzmann brain” paradox is largely ameliorated in quantum many-body systems, in contrast with the classical setting.
Studies of many-body non-Hermitian parity-time (PT)-symmetric quantum systems are attracting a lot of interest due to their relevance in research areas ranging from quantum optics and continuously monitored dynamics to Euclidean wormholes in quantum gravity and dissipative quantum chaos. While a symmetry classification of non-Hermitian systems leads to 38 universality classes, we show that, under certain conditions, PT-symmetric systems are grouped into 24 universality classes. We identify 14 of them in a coupled two-site Sachdev-Ye-Kitaev (SYK) model and confirm the classification by spectral analysis using exact diagonalization techniques. Intriguingly, in 4 of these 14 universality classes, AIII ν , BDI ν † , BDI + + ν , and CI − − ν , we identify a basis in which the SYK Hamiltonian has a block structure in which some blocks are rectangular, with ν ∈ N the difference between the number of rows and columns. We show analytically that this feature leads to the existence of ν robust purely eigenvalues, whose level statistics follow the predictions of Hermitian random matrix theory for classes A, AI, BDI, and CI, respectively. We have recently found that this ν is a topological invariant, so these classes are topological. By contrast, nontopological real eigenvalues display a crossover between Hermitian and non-Hermitian level statistics. Similarly to the case of Lindbladian dynamics, the reduction of universality classes leads to unexpected results, such as the absence of Kramers degeneracy in a given sector of the theory. Another novel feature of the classification scheme is that different sectors of the PT-symmetric Hamiltonian may have different symmetries. Published by the American Physical Society 2024
It is shown from flux transfer event (FTE) occurrence statistics, observed as a function of MLT by the ISEE satellites, that recent 2-dimensional analytic theories of the effects of pulsed Petschek reconnection predict FTEs to contribute between 50 and 200 kV to the total reconnection voltage when the magnetosheath field points southward. The upper limit (200 kV) allows the possibility that FTEs provide all the antisunward transport of open field lines into the tail lobe. This range is compared with the voltages associated with series of FTEs signatures, as inferred from ground-based observations, which are in the range 10-60 kV. We conclude that the contribution could sometimes be made by a series of single, large events; however, the voltage is often likely to be contributed by several FTEs at different MLT.
The Sachdev–Ye–Kitaev (SYK) model is a rare example of a strongly-interacting system that is analytically tractable. Tractability arises because the model is largely structureless by design and therefore artificial: while the interaction is restricted to two-body terms, interaction matrix elements are “randomized” and therefore the corresponding interaction operator does not commute with the local density. Unlike conventional density–density-type interactions, the SYK-interaction is, in this sense, not integrable. We here investigate a variant of the (complex) SYK model, which restores this integrability. It features a randomized single-body term and a density–density-type interaction. We present numerical investigations suggesting that the model exhibits two integrable phases separated by several intermediate phases including a chaotic one. The chaotic phase carries several characteristic SYK-signatures including in the spectral statistics and the frequency scaling of the Green’s function and therefore should be adiabatically connected to the non-Fermi liquid phase of the original SYK model. Thus, our model Hamiltonian provides a bridge from the SYK-model towards microscopic realism.
The increasing precision of cosmology data in the modern era is calling for methods to allow the extraction of non-Gaussian information using tools beyond two-point statistics. The marked power spectrum has the potential to extract beyond two-point information in a computationally efficient way while using much of the infrastructure already available for the power spectrum. In this work we explore the marked power spectrum from an analytical perspective. In particular, we explore a low-order polynomial for the mark that allows us to better control the theoretical uncertainties and we show that with minimal new degrees of freedom the analytical results match measurements from N-body simulations for both the matter field and biased tracers in redshift space. Finally, we show that even within the limited forms of mark that we consider, there are degeneracies that can be broken by inclusion of the marked auto-spectrum or the cross-spectrum with the unmarked field. I n conclusion, we discuss future theoretical developments that would enable us to apply this approach to survey data.
Owing to its excellent mechanical properties and stability under high temperature and neutron irradiation conditions, SiC/SiC composites have emerged as a promising material for light water reactors (LWRs) in the development of accident-tolerant fuel (ATF) systems. Structural integrity and retention of hermeticity are two crucial requirements for SiC/SiC claddings during normal operations, and both of them are closely related to the proportional limit stress (PLS) of the material. Understanding the behavior of SiC/SiC composites under multiaxial stress states and developing a probabilistic approach for evaluating the structural vulnerability are of paramount importance for reliability-based analysis and design of SiC/SiC composite claddings. So far, there has been very limited effort towards experimental and analytical investigations of probabilistic failure of SiC/SiC claddings. This critical knowledge gap motivates this research. A probabilistic failure criterion for SiC/SiC composites under multi-axial loading is developed, and this criterion is incorporated into reliability analysis of the structural integrity of SiC/SiC fuel cladding. The research consists of two parts: 1) experimental investigation of multiaxial failure behavior of SiC/SiC composites, and 2) theoretical modeling of time-dependent probabilistic failure of SiC/SiC cladding. In the experimental investigation, the PLS is determined through the examination of stress-strain response, the acoustic emission measurement, as well as the X-ray computed tomography. The theoretical framework is derived by combin- ing the finite weakest-link statistical model and the subcritical damage growth model. This theoretical model captures the time-dependent failure mechanism of the material, which has a major consequence for predicting the lifetime distribution of the cladding. Meanwhile, the model also predicts that the failure statistics of the cladding depends strongly on the cladding length. The results of the multiaxial experiments reveal the level of statistical variation of the PLS of SiC/SiC materials under different stress states. The theoretical model provides a robust analytical tool for extrapolation of small-scale laboratory test results to the behavior of full-scale claddings. These findings establish a scientific foundation for the development of reliability-based design of SiC/SiC fuel claddings, which will play an essential role in improving the structural safety and integrity of LWRs.
Aims.In this work, we perform the first cosmological parameter analysis of the fourth release of Kilo Degree Survey (KiDS-1000) data with second- and third-order shear statistics. This paper builds on a series of studies aimed at describing the roadmap to third-order shear statistics. Methods.We derived and tested a combined model of the second-order shear statistic, namely, the COSEBIs and the third-order aperture mass statistics 〈ℳ ap 3 〉 in a tomographic set-up. We validated our pipeline withN-body mock simulations of the KiDS-1000 data release. To model the second- and third-order statistics, we used the latest version of HMCODE2020 for the power spectrum and BIHALOFITfor the bispectrum. Furthermore, we used an analytic description to model intrinsic alignments and hydro-dynamical simulations to model the effect of baryonic feedback processes. Lastly, we decreased the dimension of the data vector significantly by considering only equal smoothing radii for the 〈ℳ ap 3 〉 part of the data vector. This makes it possible to carry out a data analysis of the KiDS-1000 data release using a combined analysis of COSEBIs and third-order shear statistics. Results.We first validated the accuracy of our modelling by analysing a noise-free mock data vector, assuming the KiDS-1000 error budget, finding a shift in the maximum of the posterior distribution of the matter density parameter, ΔΩ m < 0.02 σ Ω m , and of the structure growth parameter, ΔS 8 < 0.05 σ S 8 . Lastly, we performed the first KiDS-1000 cosmological analysis using a combined analysis of second- and third-order shear statistics, where we constrained Ω m = 0.248 −0.055 +0.062 andS 8 = σ 8 √(Ω m /0.3 )= 0.772 ± 0.022. The geometric average on the errors of Ω m andS 8 of the combined statistics decreases, compared to the second-order statistic, by a factor of 2.2.
This study reports a simple modification to a commercial resin 3D printer that significantly reduces the amount of prepolymer material needed for the production of extraction sorbents. Here, the modified printing platform is demonstrated in the printing of two imidazolium-based ionic liquid (IL) monomers. Two geometries resembling a blade-type polymeric ionic liquid (PIL) sorbent used in thin-film microextraction and a fiber-type sorbent used in solid-phase microextraction (SPME) were printed. The SPME PIL sorbents were used to extract 10 organic contaminants, including plasticizers, antimicrobial agents, UV filters, and pesticides, from water followed by high-performance liquid chromatographic (HPLC) analysis. To compare the extraction performance of the SPME sorbents, seven fibers printed with the same prepolymer composition from the same printing batch as well as different batches were evaluated. The results revealed highly reproducible extraction efficiencies for all tested sorbents with no statistical difference in their extraction performance. Method validation showed acceptable linearity ($R$ 2 > 0.92) for all analytes with limits of detection and limits of quantification ranging from 0.13 to 45 μg L –1 and 0.43 to 150 μg L –1 , respectively.