Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data logging”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental and modeling studies of PR and ND oxalate solubility to high ionic strengths: Insight into actinide(III) oxalates

Actinide oxalates are chemical compounds important to nuclear industry, ranging from actinide separation in waste reprocessing, to production of specialty actinides, and to disposal of high level nuclear waste (HLW) and spent nuclear fuel (SNF). In this study, the solubility constants for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O by performing solubility experiments in HNO 3 and mixtures of HNO 3 and H 2 C 2 O 4 at 23.0 ± 0.2 °C have been determined. The targeted starting materials, Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O, were successfully synthesized at room temperature using PrCl 3 , NdCl 3 and oxalic acid as the source metrials. Then, we utilized the targeted solubility-controlling phases to conduct solubility measurements. There was no phase change over the entire periods of experiments, demonstrating that Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O were the solubility-controlling phases in our respective experiments. Based on our experimental data, we have developed a thermodynamic model for Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O in the mixtures of HNO 3 and H 2 C 2 O 4 to high ionic strengths. The model for Pr 2 (C 2 O 4 ) 3 •10H 2 O reproduces well the reported experimental data for Pu 2 (C 2 O 4 ) 3 •10H 2 O, which are not utilized for the model development, demonstrating that Pr(III) is an excellent analog for Pu(III). Similarly, the model for Nd 2 (C 2 O 4 ) 3 •10H 2 O reproduces the solubility of Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The Pitzer model was used for the calculation of activity coefficients. Based on the published, well established model for dissociation constants for oxalic acid and stability constants for actinide-oxalate complexes [i.e., AmC 2 O 4 + , and Am(C 2 O 4 ) 2 - ] to high ionic strengths, we have obtained the solubility constants (log 10 K 0 ) for the following reactions at 25 °C, Pr 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Pr 3+ + 3C 2 O 4 2- + 10H 2 O(l). Nd 2 (C 2 O 4 ) 3 •10H 2 O ⇌ 2Nd 3+ + 3C 2 O 4 2- + 10H 2 O(l). to be -30.82 ± 0.30 (2σ), and - 31.14 ± 0.35 (2σ), respectively. These values for can be directly applied to Pu 2 (C 2 O 4 ) 3 •10H 2 O, Am 2 (C 2 O 4 ) 3 •10H 2 O and Cm 2 (C 2 O 4 ) 3 •10H 2 O. The model established for actinide oxalates by this study provides the needed knowledge with regard to solubilities of actinide/REE oxalates at various ionic strengths, and is expected to find applications in many fields, including the geological disposal of nuclear waste and the mobility of REE under the surface conditions, as Pr 2 (C 2 O 4 ) 3 •10H 2 O and Nd 2 (C 2 O 4 ) 3 •10H 2 O can be regarded as the pure Pr and Nd end-members of deveroite, a recently discovered natural REE oxalate with the following stoichiometry, (Ce 1.01 Nd 0.33 La 0.32 Pr 0.11 Y 0.11 Sm 0.01 Pb 0.04 U 0.03 Th 0.01 Ca 0.04 ) 2.01 (C 2 O 4 ) 2.99 •9.99H 2 O. Regarding its importance in the geological disposal of nuclear waste, Am 2 (C 2 O 4 ) 3 •10H 2 O/Pu 2 (C 2 O 4 ) 3 •10H 2 O/Cm 2 (C 2 O 4 ) 3 •10H 2 O can be the source-term phase for actinides, as demonstrated by the instance in the disposal in clay/shale formations. This is exemplified by the stability of Am 2 (C 2 O 4 ) 3 •10H 2 O in comparison with Am(OH) 3 (am), Am(OH) 3 (s) and AmCO 3 (OH)(s) under the relevant geological repository conditions.

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

A broad-host-range event detector: expanding and quantifying performance between Escherichia coli and Pseudomonas species

Modern microbial biodesign relies on the principle that well-characterized genetic parts can be reused and reconfigured for different functions. However, this paradigm has only been successful in a limited set of hosts, mostly comprised from common lab strains of Escherichia coli . It is clear that new applications such as chemical sensing and event logging in complex environments will benefit from new host chassis. This study quantitatively compared how the same chemical event logger performed across four strains and three different microbial species. An integrase-based sensor and memory device was operated by two representative soil Pseudomonads— Pseudomonas fluorescens SBW25 and Pseudomonas putida DSM 291. Quantitative comparisons were made between these two non-traditional hosts and two benchmark E. coli chassis including the probiotic Nissle 1917 and common cloning strain DH5α. The performance of sensor and memory components changed according to each host, such that a clear chassis effect was observed and quantified. These results were obtained via fluorescence from reporter proteins that were transcriptionally fused to the integrase and downstream recombinant region and via data-driven kinetic models. The Pseudomonads proved to be acceptable chassis for the operation of this event logger, which outperformed the common E. coli DH5α in many ways. This study advances an emerging frontier in synthetic biology that aims to build broad-host-range devices and understand the context by which different species can execute programmable genetic operations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Wellbore Stability and Mud Loss Management in Geothermal Drilling: Optimizing Mud Weight to Mitigate Tensile Wellbore Fracturing at The Geysers, California

As part of a U.S. Department of Energy (DOE) Geothermal Technologies Office-funded initiative, Geysers Power Company, LLC, a subsidiary of Calpine Corporation, has been working to enhance drilling performance at the world’s largest geothermal field, The Geysers, in northern California. In a recent drilling operation of the GDC-36 well, excessive mud losses were encountered, initially addressed through repeated but largely ineffective cement plugging. Ultimately, the most effective strategy was to drill blind through the loss zones, made feasible by the high rate of penetration (ROP) achieved with PDC bits, allowing significant progress before the mud tanks were depleted and water-sensitive argillic formation layers could collapse. In response to these challenges, the project team explored alternative methods to minimize downtime and risks associated with cement plugging and continuous mud loss and to contemplate the driving mechanisms for the losses. Wellbore imaging using Formation MicroImager (FMI) and Ultrasonic Borehole Imager (UBI) tools revealed longitudinal tensile fractures, which were attributed to mud weights exceeding the minimum circumferential stress resulting from the native stress field and formation pressure. This study examines the mud losses encountered and leverages wellbore imaging data to understand the mechanisms behind mud induced tensile fracturing in specific rock facies. Understanding fracture behavior across different lithologies is crucial, as fractures within the reservoir can enhance steam migration throughout the system. The reservoir at The Geysers lies within the Mesozoic Franciscan Assemblage, a tectonic mélange formed by subduction. It consists of metamorphosed turbidite sandstone (greywacke) and mudstone (argillite), oceanic upper crust (including greenstone and chert), and serpentinized ultramafic rocks - each exhibiting distinct geomechanical fracturing properties. The structural fabric of the Franciscan Assemblage was shaped by low-angle Mesozoic thrust faulting and later overprinted by sub-vertical strike-slip structures related to the Pacific-North American plate boundary. A wellbore stability model was developed using core measurements and logs to simulate fracturing scenarios during drilling under varying stress conditions. These simulations guided the development of an optimized mud weight management strategy that should enable adaptive adjustments during drilling, reducing the likelihood of tensile fracturing and mud losses, ultimately improving operational efficiency.

15 GEOTHERMAL ENERGY↗
DOE OSTI2021

encore : an O ( N g2) estimator for galaxy N -point correlation functions

ABSTRACT We present a new algorithm for efficiently computing the N-point correlation functions (NPCFs) of a 3D density field for arbitrary N. This can be applied both to a discrete spectroscopic galaxy survey and a continuous field. By expanding the statistics in a separable basis of isotropic functions built from spherical harmonics, the NPCFs can be estimated by counting pairs of particles in space, leading to an algorithm with complexity $\mathcal {O}(N_\mathrm{g}^2)$ for Ng particles, or $\mathcal {O}(N_\mathrm{FFT}\log N_\mathrm{FFT})$ when using a Fast Fourier Transform with NFFT grid-points. In practice, the rate-limiting step for N > 3 will often be the summation of the histogrammed spherical harmonic coefficients, particularly if the number of radial and angular bins is large. In this case, the algorithm scales linearly with Ng. The approach is implemented in the encore code, which can compute the 3PCF, 4PCF, 5PCF, and 6PCF of a BOSS-like galaxy survey in ${\sim}100$ CPU-hours, including the corrections necessary for non-uniform survey geometries. We discuss the implementation in depth, along with its GPU acceleration, and provide practical demonstration on realistic galaxy catalogues. Our approach can be straightforwardly applied to current and future data sets to unlock the potential of constraining cosmology from the higher point functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

B/PS bulges in DESI Legacy edge-on galaxies – I. Sample building

ABSTRACT We present the biggest up-to-date sample of edge-on galaxies with boxy/peanut-shaped (B/PS) bulges and X-structures. The sample was prepared using images from the Dark Energy Spectroscopic Instrument (DESI) Legacy catalogue and contains about 2000 galaxies. To find suitable candidates in catalogue, we made the assumption that the residues (original images minus model) of galaxies with B/PS bulges should exhibit a characteristic X-shape. Galaxies with such features were selected by eye and then used as input data for a neural network training, which was applied to a bigger sample of edge-on galaxies. Using the available data and the photometric models from the literature, we investigated the observational and statistical properties of the sample created. Comparing the B/D ratios for galaxies with and without B/PS bulges, we found that the B/D ratio for galaxies from our sample is statistically higher, with typical values in the range ≈0.2–0.5 depending on the decomposition procedure. We studied how the opening angles φ of the X-structure and the length of its rays are distributed in the formed sample and found them to be consistent with previous measurements and predictions from N-body models, e.g. $\varphi \gtrsim 25{^\circ }$, but measured here for a much larger number of galaxies. We found a sharp increase in the B/PS bulge fraction for stellar masses log M⋆ ≳ 10.4, but for edge-on galaxies, which complements the results of previous works. The sample can be used in future work to test various bar models and their relationship with B/PS bulges, as well as to study their stability and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Stellar population models based on the SDSS-IV MaStar library of stellar spectra – I. Intermediate-age/old models

ABSTRACT We use the first release of the SDSS/MaStar stellar library comprising ∼9000, high S/N spectra, to calculate integrated spectra of stellar population models. The models extend over the wavelength range 0.36–1.03 µm and share the same spectral resolution ($R\sim 1800$) and flux calibration as the SDSS-IV/MaNGA galaxy data. The parameter space covered by the stellar spectra collected thus far allows the calculation of models with ages and chemical composition in the range $\rm {\mathit{ t}\gt 200 \,Myr, -2 \lt = [Z/H] \lt = + 0.35}$, which will be extended as MaStar proceeds. Notably, the models include spectra for dwarf main-sequence stars close to the core H-burning limit, as well as spectra for cold, metal-rich giants. Both stellar types are crucial for modelling λ > 0.7 µm absorption spectra. Moreover, a better parameter coverage at low metallicity allows the calculation of models as young as 500 Myr and the full account of the blue horizontal branch phase of old populations. We present models adopting two independent sets of stellar parameters (Teff, log g, [Z/H]). In a novel approach, their reliability is tested ‘on the fly’ using the stellar population models themselves. We perform tests with Milky Way and Magellanic Clouds globular clusters, finding that the new models recover their ages and metallicities remarkably well, with systematics as low as a few per cent for homogeneous calibration sets. We also fit a MaNGA galaxy spectrum, finding residuals of the order of a few per cent comparable to the state-of-art models, but now over a wider wavelength range.

Maraston, C.↗
DOE OSTI2021

Atomic Iron and Nickel in the Coma of C/1996 B2 (Hyakutake): Production Rates, Emission Mechanisms, and Possible Parents

Abstract Two papers recently reported the detection of gaseous nickel and iron in the comae of over 20 comets from observations collected over two decades, including interstellar comet 2I/Borisov. To evaluate the state of the laboratory data in support of these identifications, we reanalyzed archived spectra of comet C/1996 B2 (Hyakutake), one of the nearest and brightest comets of the past century, using a combined experimental and computational approach. We developed a new, many-level fluorescence model that indicates that the fluorescence emissions of Fe I and Ni I vary greatly with heliocentric velocity. Combining this model with laboratory spectra of an Fe-Ni plasma, we identified 22 lines of Fe I and 14 lines of Ni I in the spectrum of Hyakutake. Using Haser models, we estimate the nickel and iron production rates as Q Ni = (2.6–4.1) × 10 22 s −1 and Q Fe = (0.4–2.8) × 10 23 s −1 . From derived column densities, the Ni/Fe abundance ratio log 10 [Ni/Fe] = −0.15 ± 0.07 deviates significantly from solar abundance ratios, and it is consistent with the ratios observed in solar system comets. Possible production and emission mechanisms are analyzed in the context of existing laboratory measurements. Based on the observed spatial distributions, excellent fluorescence model agreement, and Ni/Fe ratio, our findings support an origin consisting of a short-lived unknown parent followed by fluorescence emission. Our models suggest that the strong heliocentric velocity dependence of the fluorescence efficiencies can provide a meaningful test of the physical process responsible for the Fe I and Ni I emission.

Bromley, S. J. (ORCID:0000000321108152)↗
DOE OSTI2022

Risk matrix for legacy wells within the Area of Review (AoR) of Carbon Capture & Storage (CCS) projects

The success of CCS depends on the capacity, injectivity, and confinement by the storage medium. Thousands of wells drilled over the past century with the intention to find the trapped oil and gas may penetrate the containment seals. These wells may provide leakage pathways for the CO 2 to escape and contaminate the underground sources of drinking water (USDW) or reach the surface in the worst-case scenario. Identifying the risky wells penetrating the containment seals and predicting their current as well as future well integrity is the most challenging task when limited data is available. Here, this paper proposes a unique methodology for risk assessment of the wells penetrating the containment seals based on the proximity of these wells from the proposed injection location, mechanical integrity, and accessibility of these wells over the lifecycle of the CCS project. This helps in identifying the wells which need immediate attention from the wells that need little to no attention. It also highlights the corrective actions necessary for the success of the CCS project as well as help estimate the approximate cost required to perform the corrective actions. This methodology focuses on all wells (producers, injectors, orphan, abandoned, water, stratigraphic, etc) while the majority of the studies found in literature focused on wells with sustained casing pressure (SCP) reports and cement bond logs (CBL). The proposed risk matrix, if applied to future CCS projects across the globe, will uniformly categorize the wells within the Area of Review (AoR).

58 GEOSCIENCES↗
DOE OSTI2026

Ride Pingo to Transit

In this project, we developed an on-demand microtransit first- and last-mile service. To integrate the service with fixed-route transit, we developed a feature called Transit Connect that prioritized riders’ on-time arrival at the transit station over other service requirements. We first prototyped service and related algorithms in a simulated environment, and then piloted the service in the city of Kent, Washington. Our algorithm incorporates request-specific hard drop-off deadlines to ensure timely arrivals for transit transfers. In the pilot, these constraints were obtained from GTFS Realtime data to accurately determine the schedule of the transit and the location of the stations. This approach introduced the ability to accept or decline new requests based on the timing of transit connections for these new requests and connection status of onboarding customers. The pilot (called “Ride Pingo to Transit”) deployed a fleet of three 14-person vans, ran from September 2021 to March 2023, and served a total of 21,329 trips. Transit Connect was offered for drop-offs at both Kent Station and the Kent Valley hub. In total, 2,844 such trips were completed. This dataset was collected from our pilot, which includes the following: - Requests: List of all trip requests, including those that were actually served and those not materialized. - Fleet: Daily vehicle service logs. - Service details: Daily vehicle stop logs (boarding and alighting). - Trip types: First mile, last mile, or point-to-point. ![pingo to transit](pingo-to-transit.jpg)

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

TESS Delivers Five New Hot Giant Planets Orbiting Bright Stars from the Full-frame Images

We present the discovery and characterization of five hot and warm Jupiters—TOI-628 b (TIC 281408474; HD 288842), TOI-640 b (TIC 147977348), TOI-1333 b (TIC 395171208, BD+47 3521A), TOI-1478 b (TIC 409794137), and TOI-1601 b (TIC 139375960)—based on data from NASA’s Transiting Exoplanet Survey Satellite (TESS). The five planets were identified from the full-frame images and were confirmed through a series of photometric and spectroscopic follow-up observations by the TESS Follow-up Observing Program Working Group. The planets are all Jovian size (R {sub P} = 1.01–1.77 R {sub J}) and have masses that range from 0.85 to 6.33 M {sub J}. The host stars of these systems have F and G spectral types (5595 ≤ T {sub eff} ≤ 6460 K) and are all relatively bright (9.5 < V < 10.8, 8.2 < K < 9.3), making them well suited for future detailed characterization efforts. Three of the systems in our sample (TOI-640 b, TOI-1333 b, and TOI-1601 b) orbit subgiant host stars (log g < 4.1). TOI-640 b is one of only three known hot Jupiters to have a highly inflated radius (R {sub P} > 1.7 R {sub J}, possibly a result of its host star’s evolution) and resides on an orbit with a period longer than 5 days. TOI-628 b is the most massive, hot Jupiter discovered to date by TESS with a measured mass of 6.31{sub −0.30}{sup +0.28} M {sub J} and a statistically significant, nonzero orbital eccentricity of e = 0.074{sub −0.022}{sup +0.021}. This planet would not have had enough time to circularize through tidal forces from our analysis, suggesting that it might be remnant eccentricity from its migration. The longest-period planet in this sample, TOI-1478 b (P = 10.18 days), is a warm Jupiter in a circular orbit around a near-solar analog. NASA’s TESS mission is continuing to increase the sample of well-characterized hot and warm Jupiters, complementing its primary mission goals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reassessing long-standing meteorological records: an example using the national hottest day in Ireland

This analysis highlights the potential value in reanalysing early national meteorological records from around the world. These were oftentimes measured via techniques that preceded standardisation of instrumentation and methods of observation and thus could be subject to considerable biases and uncertainties. This analysis uses the techniques pioneered by WMO record assessment teams. The highest currently recognised air temperature (33.3 °C) ever recorded in the Republic of Ireland was logged at Kilkenny Castle in 1887. The original observational record however no longer exists. Given that Ireland is now the only country in Europe to have a national heat record that was set in the 19th century, a reassessment of the verity of this record is both timely and valuable. The present analysis undertakes a fundamental reassessment of the plausibility of the 1887 temperature record using methods similar to those used to assess various weather extremes under WMO auspices over recent years. Specifically, we undertake an inter-station reassessment using sparse available records and make recourse to the new and improved 20CRv3 sparse-input reanalysis product. Neither surrounding available stations nor the reanalysis offer substantive support for the Kilkenny record of 33.3 °C being correct. Moreover, recent data rescue efforts have uncovered several earlier extreme values, one of which exceeds the Kilkenny value (33.5 °C on 16 July 1876 recorded at the Phoenix Park). However, the sparsity of early observational networks, a distinct lack of synoptic support from 20CRv3 for many of the extreme heat values, and the fact that these measurements were obtained using non-standard exposures lead us to conclude that there is grossly insufficient evidence to support any of these 19th century extremes as robust national heat record candidates. Data from the early 20th century onwards benefit from a denser network of stations undertaking measurements in a more standardised manner, many under the direct auspices of Met Éireann and its predecessors, adhering to WMO guidance and protocols. This enables more robust cross-checking of records. We argue that the Met Éireann-recognised 20th century heat record from Boora in 1976 is verified as the most plausible robust national temperature record based upon the synoptic situation and comparisons with nearby neighbouring stations. This measurement of 32.5°C thus likely constitutes the highest reliably recorded temperature measurement in the Republic of Ireland. Ultimately, the formal decision on any reassessment and reassignment of the national record rests with the national meteorological service, Met Éireann.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

The importance of round-robin validation when assessing machine-learning-based vertical extrapolation of wind speeds

The extrapolation of wind speeds measured at a meteorological mast to wind turbine rotor heights is a key component in a bankable wind farm energy assessment and a significant source of uncertainty. Industry-standard methods for extrapolation include the power-law and logarithmic profiles. The emergence of machine-learning applications in wind energy has led to several studies demonstrating substantial improvements in vertical extrapolation accuracy in machine-learning methods over these conventional power-law and logarithmic profile methods. In all cases, these studies assess relative model performance at a measurement site where, critically, the machine-learning algorithm requires knowledge of the rotor-height wind speeds in order to train the model. This prior knowledge provides fundamental advantages to the site-specific machine-learning model over the power-law and log profiles, which, by contrast, are not highly tuned to rotor-height measurements but rather can generalize to any site. Furthermore, there is no practical benefit in applying a machine-learning model at a site where winds at the heights relevant for wind energy production are known; rather, its performance at nearby locations (i.e., across a wind farm site) without rotor-height measurements is of most practical interest. To more fairly and practically compare machine-learning-based extrapolation to standard approaches, we implemented a round-robin extrapolation model comparison, in which a random-forest machine-learning model is trained and evaluated at different sites and then compared against the power-law and logarithmic profiles. We consider 20 months of lidar and sonic anemometer data collected at four sites between 50 and 100 km apart in the central United States. We find that the random forest outperforms the standard extrapolation approaches, especially when incorporating surface measurements as inputs to include the influence of atmospheric stability. When compared at a single site (the traditional comparison approach), the machine-learning improvement in mean absolute error was 28 % and 23 % over the power-law and logarithmic profiles, respectively. Using the round-robin approach proposed here, this improvement drops to 20 % and 14 %, respectively. These latter values better represent practical model performance, and we conclude that round-robin validation should be the standard for machine-learning-based wind speed extrapolation methods.

17 WIND ENERGY↗
DOE OSTI2023

Hierarchical off-diagonal low-rank approximation of Hessians in inverse problems, with application to ice sheet model initialization

Obtaining lightweight and accurate approximations of discretized objective functional Hessians in inverse problems governed by partial differential equations (PDEs) is essential to make both deterministic and Bayesian statistical large-scale inverse problems computationally tractable. The cubic computational complexity of dense linear algebraic tasks, such as Cholesky factorization, that provide a means to sample Gaussian distributions and determine solutions of Newton linear systems is a computational bottleneck at large-scale. These tasks can be reduced to log-linear complexity by utilizing hierarchical off-diagonal low-rank (HODLR) matrix approximations. In this work, we show that a class of Hessians that arise from inverse problems governed by PDEs are well approximated by the HODLR matrix format. In particular, we study inverse problems governed by PDEs that model the instantaneous viscous flow of ice sheets. In these problems, we seek a spatially distributed basal sliding parameter field such that the flow predicted by the ice sheet model is consistent with ice sheet surface velocity observations. Here, we demonstrate the use of HODLR Hessian approximation to efficiently sample the Laplace approximation of the posterior distribution with covariance further approximated by HODLR matrix compression. Computational studies are performed which illustrate ice sheet problem regimes for which the Gauss–Newton data-misfit Hessian is more efficiently approximated by the HODLR matrix format than the low-rank (LR) format. We then demonstrate that HODLR approximations can be favorable, when compared to global LR approximations, for large-scale problems by studying the data-misfit Hessian associated with inverse problems governed by the first-order Stokes flow model on the Humboldt glacier and Greenland ice sheet.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Asteroseismological analysis of the polluted ZZ Ceti star G 29 – 38 with TESS

G29 – 38 (TIC 422526868) is one of the brightest (V = 13.1) and closest (d = 17.51 pc) pulsating white dwarfs with a hydrogen-rich atmosphere (DAV/ZZ Ceti class). It was observed by the TESS spacecraft in sectors 42 and 56. The atmosphere of G 29 – 38 is polluted by heavy elements that are expected to sink out of visible layers on short time-scales. Here, the photometric TESS data set spans ~51 d in total, and from this, we identified 56 significant pulsation frequencies, that include rotational frequency multiplets. In addition, we identified 30 combination frequencies in each sector. The oscillation frequencies that we found are associated with g-mode pulsations, with periods spanning from ~ 260 to ~ 1400 s. We identified rotational frequency triplets with a mean separation δν ℓ = 1 of 4.67 μHz and a quintuplet with a mean separation δν ℓ = 2 of 6.67 μHz, from which we estimated a rotation period of about 1.35 ± 0.1 d. We determined a constant period spacing of 41.20 s for ℓ = 1 modes and 22.58 s for ℓ = 2 modes. We performed period-to-period fit analyses and found an asteroseismological model with M * /M ⊙ = 0.632 ± 0.03, K, and log g = 8.048 ± 0.005 (with a hydrogen envelope mass of M H ~ 5.6 × 10 –5 M * ), in good agreement with the values derived from spectroscopy. We obtained an asteroseismic distance of 17.54 pc, which is in excellent agreement with that provided by Gaia (17.51 pc).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robust Parameter Design on Dual Stochastic Response Models With Constrained Bayesian Optimization

In engineering system design, minimizing the variations of the quality measurements while guaranteeing their overall quality up to certain levels, namely the robust parameter design (RPD), is crucial. Recent works have dealt with the design of a system whose response-control variables relationship is a deterministic function with a complex shape and function evaluation is expensive. In this work, we propose a Bayesian optimization method for the RPD of stochastic functions. Dual stochastic response models are carefully designed for stochastic functions. The heterogeneous variance of the sample mean is addressed by the predictive mean of the log variance surrogate model in a two-step approach. We establish an acquisition function that favors exploration across the feasible and optimality-improvable regions to effectively and efficiently solve the stochastic constrained optimization problem. Further, the performance of our proposed method is demonstrated by the extensive numerical and case studies. Note to Practitioners-Many manufacturing processes involve undesirable variations, which create variations in the final products. For example, many emerging manufacturing processes, such as nanomanufacturing, involve complex physical and chemical dynamics and transformation, creating variations in the manufacturing output. In such processes, it is crucial to design the manufacturing processes or products so that they have minimum variations in their quality. Meanwhile, it is also important to maintain the overall quality of the designed processes or products. Furthermore, acquiring data from many advanced manufacturing processes is often very costly, especially in the designing stage. In this work, we propose a data-driven method that automatically finds the best setting of manufacturing processes or products with the minimum variations of quality and a given constraint on the average quality satisfied. Our proposed method is used before conducting every experiment; It analyzes the historical data from previous experiments and provides a setting to be used in the next experiment. Our proposed method efficiently utilizes the historical data, and thus finds the best robust setting by conducting only a small number of experiments.

42 ENGINEERING↗
DOE OSTI2026

Block encoding of the three-dimensional heterogeneous Poisson equation with application to fracture flow

Quantum linear system (QLS) algorithms offer the potential to solve large-scale linear systems exponentially faster than classical methods. However, applying QLS algorithms to real-world problems remains challenging due to issues such as state preparation, data loading, and efficient information extraction. In this work, we study the feasibility of applying QLS algorithms to solve discretized three-dimensional (3D) heterogeneous Poisson equations, with specific examples relating to groundwater flow through geologic fracture networks. We explicitly construct a block encoding for the 3D heterogeneous Poisson matrix by leveraging the sparse local structure of the discretized operator. While classical solvers benefit from preconditioning, we show that block encoding the system matrix and preconditioner separately does not improve the effective condition number that dominates the QLS run-time. This differs from classical approaches where the preconditioner and the system matrix can often be implemented independently. Nevertheless, due to the structure of the problem in three dimensions, the quantum algorithm achieves a run-time of 𝑂⁡(𝑁 2/3 polylog 𝑁 ⋅log (1/𝜖)), outperforming the best classical methods (with run times of 𝑂⁡(𝑁⁢log 𝑁 ⋅log (1/𝜖))) and offering exponential memory savings. These results highlight both the promise and limitations of QLS algorithms for practical scientific computing, and point to effective condition-number reduction as a key barrier in achieving quantum advantages.

58 GEOSCIENCES↗
DOE OSTI2020

Characterization and identification of HPC applications at leadership computing facility

High Performance Computing (HPC) is an important method for scientific discovery via large-scale simulation, data analysis, or artificial intelligence. Leadership-class supercomputers are expensive, but essential to run large HPC applications. The Petascale era of supercomputers began in 2008, with the first machines achieving performance in excess of one petaflops, and with the advent of new supercomputers in 2021 (e.g., Aurora, Frontier), the Exascale era will soon begin. However, the high theoretical computing capability (i.e., peak FLOPS) of a machine is not the only meaningful target when designing a supercomputer, as the resources demand of applications varies. A deep understanding of the characterization of applications that run on a leadership supercomputer is one of the most important ways for planning its design, development and operation. In order to improve our understanding of HPC applications, user demands and resource usage characteristics, we perform correlative analysis of various logs for different subsystems of a leadership supercomputer. This analysis reveals surprising, sometimes counter-intuitive patterns, which, in some cases, conflicts with existing assumptions, and have important implications for future system designs as well as supercomputer operations. For example, our analysis shows that while the applications spend significant time on MPI, most applications spend very little time on file I/O. Combined analysis of hardware event logs and task failure logs show that the probability of a hardware FATAL event causing task failure is low. Combined analysis of control system logs and file I/O logs reveals that pure POSIX I/O is used more widely than higher level parallel I/O. Based on holistic insights of the application gained through combined and co-analysis of multiple logs from different perspectives and general intuition, we engineer features to "fingerprint" HPC applications. We use t-SNE (a machine learning technique for dimensionality reduction) to validate the explainability of our features and finally train machine learning models to identify HPC applications or group those with similar characteristic. To the best of our knowledge, this is the first work that combines logs on file I/O, computing, and inter-node communication for insightful analysis of HPC applications in production.

Liu, Zhengchun↗
DOE OSTI2022

AmeriFlux US-Cst Crossett Experimental Forest

This is the AmeriFlux version of the carbon flux data for the site US-Cst Crossett Experimental Forest. Site Description - The study takes place in the Crossett Experimental Forest (CEF, 32°2′ N, 91°57′ W), which was established in 1934 with the objective of developing effective management protocols for loblolly and shortleaf pine forests in the Gulf Coastal Plain. Mean annual temperature and precipitation are 17.6 °C and 1,410 mm, respectively. Most of the soils on and near the CEF are silt loams (primarily Glossaquic Fragiudalfs). The tower is situated in the the old-growth pine management sector of the CEF, which emphasizes management for old-growth-like conditions, with the goal of emulating the structural and compositional attributes of historic forests from this region using a combination of timber harvesting, prescribed fire, and other tools. This includes periodic commercial logging of intermediate-sized loblolly pines to support further restoration treatments, such as prescribed fire, hardwood and exotic species competition control, and the increase of plant species more compatible with fire (e.g., shortleaf pine and understory grasses)

Novick, Kim↗