Search NASA⌕ Search

SEARCH · Search NASA

Results for “Level set method”

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 199 records · Page 11

Level-set Modeling Simulations of Chemical Vapor Infiltration for Ceramic Matrix Composites Manufacturing

Silicon-carbide (SiC) reinforced ceramic matrix composites (CMCs) are a key enabling technology to reduce fuel consumption and emissions of gas turbine engines. In one manufacturing approach, chemical vapor infiltration (CVI) is limited to only coating SiC fibers. The preform is then fabricated using a lay-up of basic plys or 2-D woven sheets composed of the precoated fibers. At the other extreme, CVI is used to completely densify a 3-D woven preform shaped almost like the gas turbine component itself. The latter approach is more suitable for highly engineered components which sit directly in the gas path of the engine, for example, a high pressure turbine blade. In this case, the geometry is necessarily complex for aerodynamic, stress, and lifing (multi-physics) requirements. Presently, optimizing the CVI-dominated manufacturing approach is largely by trial-and-error. In this work, a first-principles modeling of CVI is performed to realize optimization of SiC/SiC CMC manufacturing. The modeling is based on a level-set framework to describe the interface between the vapor and solid phases. A finite-difference numerical scheme using an immersed boundary method is developed for fixed, structured meshes. Massively parallel direct numerical simulations (DNS) of CVI through fiber-woven geometries are performed using one-step chemistry, and over a range of Thiele moduli. Illustrative applications of the resulting large DNS data sets are given, including the development of fiber-weave specific infiltration models and structure functions for mean-field (porous media) Computational Fluid Dynamics (CFD) simulations of CVI.

36 MATERIALS SCIENCE↗

Propagating sample variance uncertainties in redshift calibration: simulations, theory, and application to the COSMOS2015 data

ABSTRACT Cosmological analyses of galaxy surveys rely on knowledge of the redshift distribution of their galaxy sample. This is usually derived from a spectroscopic and/or many-band photometric calibrator survey of a small patch of sky. The uncertainties in the redshift distribution of the calibrator sample include a contribution from shot noise, or Poisson sampling errors, but, given the small volume they probe, they are dominated by sample variance introduced by large-scale structures. Redshift uncertainties have been shown to constitute one of the leading contributions to systematic uncertainties in cosmological inferences from weak lensing and galaxy clustering, and hence they must be propagated through the analyses. In this work, we study the effects of sample variance on small-area redshift surveys, from theory to simulations to the COSMOS2015 data set. We present a three-step Dirichlet method of resampling a given survey-based redshift calibration distribution to enable the propagation of both shot noise and sample variance uncertainties. The method can accommodate different levels of prior confidence on different redshift sources. This method can be applied to any calibration sample with known redshifts and phenotypes (i.e. cells in a self-organizing map, or some other way of discretizing photometric space), and provides a simple way of propagating prior redshift uncertainties into cosmological analyses. As a worked example, we apply the full scheme to the COSMOS2015 data set, for which we also present a new, principled SOM algorithm designed to handle noisy photometric data. We make available a catalogue of the resulting resamplings of the COSMOS2015 galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Benchmarking the extraction of statistical neutron capture cross sections on short-lived nuclei for applications using the β -Oslo method

Numerous scientific fields including astrophysics, nuclear power, and nuclear forensics require a knowledge of basic nuclear properties for large numbers of short-lived, radioactive isotopes far removed from stable nuclei. Neutron-capture cross sections are one such piece of nuclear data where direct measurements are not possible and theoretical predictions can vary by orders of magnitude. A recently developed indirect technique for inferring neutron capture rates, the β-Oslo method, has been introduced but not compared against a known, directly measured neutron capture cross section. To provide this benchmark, two indirect methods based on β decay and charged-particle reactions were used to extract the nuclear level density and γ-ray strength function of 51 Ti. The nuclear level density and γ-ray strength function from the two data sets were found to be equivalent and were used to extract the neutron capture cross section of 50 Ti which agrees with previous direct measurements at high neutron energies. Here, the results demonstrate the validity of the β-Oslo method for extracting neutron capture cross sections of short-lived nuclei and provide a sufficiently small uncertainty to be used in various applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multimodal parameter spaces of a complex multi-channel neuron model

One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.

97 MATHEMATICS AND COMPUTING↗

Blinding scheme for the scale-dependence bias signature of local primordial non-Gaussianity for DESI 2024

The next generation of spectroscopic surveys is expected to achieve an unprecedented level of accuracy in the measurement of cosmological parameters. To avoid confirmation bias and thereby improve the reliability of these results, blinding procedures become a standard practice in the cosmological analyses of such surveys. Blinding is especially crucial when the impact of observational systematics is important relative to the cosmological signal, and a detection of that signal would have significant implications. This is the case for local primordial non-gaussianity, as probed by the scale-dependent bias of the galaxy power spectrum at large scales that are heavily sensitive to the dependence of the target selection on the imaging quality, known as imaging systematics. We propose a blinding method for the scale-dependent bias signature of local primordial non-gaussianity at the density field level which consists in generating a set of weights for the data that replicate the scale-dependent bias. The applied blinding is predictable, and can be straightforwardly combined with other catalog-level blinding procedures that have been designed for the baryon acoustic oscillation and redshift space distortion signals. The procedure is validated through simulations that replicate data from the first year of observation of the Dark Energy Spectroscopic Instrument, but may find applications to other upcoming spectroscopic surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Anomaly detection in collider physics via factorized observables

To maximize the discovery potential of high-energy colliders, experimental searches should be sensitive to unforeseen new physics scenarios. This goal has motivated the use of machine learning for unsupervised anomaly detection. In this paper, we introduce a new anomaly detection strategy called : factorized observables for regressing conditional expectations. Our approach is based on the inductive bias of factorization, which is the idea that the physics governing different energy scales can be treated as approximately independent. Assuming factorization holds separately for signal and background processes, the appearance of nontrivial correlations between low- and high-energy observables is a robust indicator of new physics. Under the most restrictive form of factorization, a machine-learned model trained to identify such correlations will in fact converge to the optimal new physics classifier. We test on a benchmark anomaly detection task for the Large Hadron Collider involving collimated sprays of particles called jets. By teasing out correlations between the kinematics and substructure of jets, our method can reliably extract percent-level signal fractions. This strategy for uncovering new physics adds to the growing toolbox of anomaly detection methods for collider physics with a complementary set of assumptions. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Application of Point Precipitation Frequency Estimates to Watersheds

This report documents work sponsored by the U.S. Nuclear Regulatory Commission (NRC) at the Oak Ridge National Laboratory (ORNL) as part of the RES project, “Application of Point Precipitation Frequency Estimates to Watersheds.” This project was implemented as part of the Probabilistic Flood Hazard Assessment (PFHA) Research Program. The objective of the PFHA Research Program is to develop tools and guidance on the use of PFHA methods to risk-inform NRC’s licensing of new facilities as well as licensing and oversight of currently operating facilities as they relate to flooding hazards. Many nuclear power plants (NPPs) are located on or near rivers so riverine flooding hazards need to be considered in their design and operation. Probabilistic riverine flood models are important tools for realistic assessment of flooding risks. However, these models require areal estimates of the depth, duration, and frequency of rainfall distributed over the watershed, which are not often available. Point precipitation frequency estimates are more widely available. For example, the National Oceanic and Atmospheric Administration (NOAA) has published NOAA Atlas 14, which provides point precipitation frequency estimates for 5-minute through 60-day durations at average recurrence intervals of 1-year through 1,000-year. The research documented in this report addresses areal reduction factors (ARFs), which can be used to convert the widely available point precipitation frequency estimates, to estimates of areal precipitation frequency over a watershed. The most widely used ARF source is Technical Paper 29 (TP-29) published by the then U.S. Weather Bureau in 1958. However, both the methods and the underlying precipitation data used to produce TP-29 are seriously out of date. For example, due to the small gauge network available at the time of TP-29’s compilation, ARF estimates developed are only for watersheds smaller than about 400 square miles. Due to the relatively short record lengths of precipitation data available, frequency considerations could not be accurately determined. Other factors such as regional climate and seasonality were not addressed. Several newer methods have been published since TP-29 was developed and both the type and quantity of precipitation data have increased significantly, along with computational resources and analytical tools such as geographic information systems. This report reviewed and assessed the available precipitation products and methods for conducting ARF analysis. The work applied up-to-date precipitation data products and analysis methods with a novel watershed-based approach to investigate how ARF estimates vary across different methods, data sources, geographical locations, return periods, and seasons. The overall findings reported here regarding basic ARF trends are in line with other recent studies showing that ARFs decrease with increasing area, increase with increasing duration, and decrease with increasing return period. This study found significant differences among the available ARF methods. This work also found a strong geographical variability across different US hydrologic regions, suggesting that the ARF are specific to regional climate patterns and geographical characteristics and should not be applied arbitrarily to other locations. The results also reveal the importance of data record length, especially for high return level ARFs. The work reported in NUREG/CR-7271 will assist NRC staff in assessing different classes of ARF methods in conjunction with available rainfall data sets. It will also support the development of guidance for application of point precipitation data in PFHAs. It should be noted that the ARF values presented in this report for any location or region were developed for the purposes of comparing methods and investigating the factors that influence ARFs. They should not be considered official and should not be used in leu of a site-specific analysis.

54 ENVIRONMENTAL SCIENCES↗

MOSAiC Radar b1 Processing: Corrections, Calibrations, and Processing Report

The U.S. Department of Energy’s (DOE) Atmospheric Radiation Measurement (ARM) user facility deployed many instruments on board a German ice breaker, the Research Vessel (RV) Polarstern, for one year from October 2019 to October 2020. The purpose of the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign was to study the decline in the sea-ice pack around the North Pole, and what factors may be at play. After the campaign ended, efforts were undertaken to provide a calibrated radar data set for future studies. MOSAiC presented new challenges to this process, as existing methodologies often were not applicable for the frozen environment with little to no ground clutter, and concurrent engineering updates or calibrations could not be accomplished once the ship set off. Like the previous Cloud, Aerosol, and Complex Terrain Interactions (CACTI) and Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE) documentation (Hardin et al. 2020, Matthews et al. 2023), the data correction and calibration process is known in ARM as creating a “b1” datastream. This means that the radar datastreams have been well characterized to the best possible quality. This report will detail the status of the raw “a1” level data sets during the MOSAiC campaign, the corrections that were applied to create the “b1” data files, and the details of the applied methods.

54 ENVIRONMENTAL SCIENCES↗

Federal Metering Guidance

The Federal Metering Guidance documents the metering provisions in EO 13834, Efficient Federal Operations and other metering related statutory requirements. This document provides the meter-related reporting requirements for Federal agencies, specifies a method to verify the existence of currently installed building-level meters, provides guidelines on prioritizing the implementation of advanced energy and water meters, sets the schedule for agencies to update their five-year metering plans, provides the minimum requirements that agencies must meet when updating their five-year metering plans along with recommend elements for an effective plan.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Open-Access Nuclear Knowledge and Innovation Toolkit to Make Use of Nuclear World and Increase Collaboration

Nuclear engineering began as an urgent weapons development program in the early 1940s, involving work in enrichment, reactor design, shielding design, and isotope production. This work led to investigations of peaceful uses of nuclear technology, including nuclear power production and applications of radiation. Such peaceful technologies, often derived from experience gained through the weapons program, often maintained some level or attitude of secrecy. As the field continued to develop, that secrecy led to different specializations developing with different sets of terminology and methods of computation. For example, radiation safety used buildup and shielding factors while reactor designers use diffusivity and cross sections. Impressive amounts of work went into theory and application development in each sub-field over the next few decades, with much of the work left without adequate documentation or with important documentation made only within the organization. This situation has led to a loss of knowledge and experience, a lack of collaboration across nuclear sub-fields, and immense amounts of time and effort spent on redundant work. This paper proposes an open-access system that combines preserving and passing on the knowledge in all nuclear sub-fields with tools to make that knowledge useful, referred to as the Nuclear Knowledge and Innovation Toolkit (NuKIT). Knowledge would be captured in user-entered and expert-verified databases of reactors, applications of radiation, isotope production methods, and other pertinent nuclear information complete with technical information and tags that can be used to organize and investigate the data. Tools include community tools like forums and the ability to create collaborative projects; visualization and presentation tools; and an online programming system that allows code sharing, collaborative editing, database manipulation, and integration with all other tools. Finally, NuKIT is envisioned to include commercialization tools to allow freelance work and to compensate data entry/moderation work to help encourage participation. Eventually, integration with industry-standard tools like MCNP, SCALE, etc. is expected. This paper will present some proposed work in the context of its use for isotopes and radiation application, as well as possible work to be done in the context of the Combined Materials Experimental Toolkit (CoMET) being developed by the Nuclear Science User Facilities (NSUF) program. It should be noted that any work done towards this system, by any group, will advance the full system. It should also be noted that similar work previously done can be integrated into NuKIT.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Structure of Offshore Low-Level Jet Turbulence and Implications to Mesoscale-to-Microscale Coupling

This paper explores realistic nonstationary atmospheric boundary layer (ABL) turbulence arising from nonstationarity at the mesoscale, particularly within offshore low-level jets with implications to offshore wind farms, using high-fidelity multiscale large-eddy simulations (LES). To this end, we analyzed the single-point turbulence statistical structure of a North-Atlantic offshore LLJ event simulated using high-resolution LES (AMR-Wind). The nonstationary LLJ is simulated using a mesoscale-to-microscale coupled (MMC) simulation procedure involving data assimilation of mesoscale velocity and temperature data from the Weather Research and Forecasting (WRF) model. Unlike the assimilation of mesoscale velocity data into the LES, the direct assimilation of temperature profiles had a strong impact on turbulence stratification, thereby causing erroneous predictions of turbulence both above and within the jet layer. Various approaches to mitigate this effect have resulted in multiple (four) variants of this MMC strategy. Outcomes from this work clearly show that the turbulence within the low-level jet is a strong function of the MMC approach as the turbulence structure within the low-level jet is dependent on the flux of residual turbulence from outside the jet, which in turn depends on the temperature forcing history. Additionally, the turbulence predicted by all these different methods (as well as the observation data) show similar deviations from equilibrium as evidenced by comparisons with idealized atmospheric turbulence structure obtained using the same numerical method. In general, we observe that the predicted LLJ turbulence tends to differ from canonical ABL turbulence with comparable shear. Particularly, the combination of shear and turbulence observed in such nonstationary low-level turbulence cannot be matched using equilibrium settings and therefore, represents a critical use-case for both testing and leveraging meso–micro coupling strategies.

17 WIND ENERGY↗

Best practices for building energy codes compliance

Building energy codes are a proven policy mechanism to achieve economy-wide energy savings in buildings. In 2021, the International Energy Agency Energy in Buildings and Communities Programme (IEA EBC) Building Energy Codes Working Group (BECWG) set out to survey the diverse codes compliance practices across BECWG member countries. While codes vary in format and approach, attaining the energy reduction potential of building energy codes requires effective implementation and compliance. This is a large undertaking for countries, as effective compliance checking requires adequate resources, technical knowledge, capacity, and strong institutions. Many nations face the same compliance issues, such as requiring faster and easier methods to verify codes and coordinating among numerous stakeholders and levels of government. This paper explores the need for stronger institutional approaches to enforce building energy codes that will lead to code compliance. Specifically, this report sets out to address the question of what practices result in effective building energy code compliance in selected BECWG member countries.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Investigation of dose reduction in cardiac perfusion SPECT via optimization and choice of the image reconstruction strategy

Background: We investigated the extent to which the administered dose (activity) level can be reduced without sacrificing diagnostic accuracy for three reconstruction strategies for SPECT-myocardial perfusion imaging (MPI). Methods: We optimized the parameters of the three reconstruction strategies for perfusion-defect detection over a range of simulated administered dose levels using a set of hybrid studies (derived from 190 subjects) consisting of clinical SPECT-MPI data modified to contain realistic simulated lesions. The optimized strategies we considered are filtered backprojection (FBP) with no correction for degradations, ordered-subsets expectation-maximization (OS-EM) with attenuation correction (AC), scatter correction (SC), and resolution correction (RC), and OS-EM with scatter and resolution correction only. Each study was evaluated using a total perfusion deficit (TPD) score computed by the Quantitative Perfusion SPECT (QPS) software package. We conducted a receiver operating characteristics (ROC) study based on the TPD scores for each dose level and reconstruction strategy. Results: For FBP, the achieved optimum values of the area under the ROC curve (AUC) at 100%, 50%, 25%, and 12.5% of standard dose were 0.75, 0.74, 0.72, and 0.70, respectively, compared to 0.81, 0.79, 0.76, and 0.74 for OS-EM with AC–SC–RC and 0.78, 0.77, 0.74, 0.72 for OS-EM with SC–RC. Conclusions: Our results suggest that studies reconstructed by OS-EM with AC–SC–RC could possibly be reduced, on average, to 25% of the originally administered dose without causing diagnostic accuracy (AUC) to decrease below that of FBP.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Automated experimental design of safe rampdowns via probabilistic machine learning

Abstract Typically the rampdown phase of a shot consists of a decrease in current and injected power and optionally a change in shape, but there is considerable flexibility in the rate, sequencing, and duration of these changes. On the next generation of tokamaks it is essential that this is done safely as the device could be damaged by the stored thermal and electromagnetic energy present in the plasma. This works presents a procedure for automatically choosing experimental rampdown designs to rapidly converge to an effective rampdown trajectory. This procedure uses probabilistic machine learning methods paired with acquisition functions taken from Bayesian optimization. In a set of 2022 experiments at DIII-D, the rampdown designs produced by our method maintained plasma control down to substantially lower current and energy levels than are typically observed. The actions predicted by the model significantly improved as the model was able to explore over the course of the experimental campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Conditional deep generative models for simultaneous simulation and reconstruction of entire events

We extend the particle-flow neural assisted simulations (arnassus) framework of fast simulation and reconstruction to entire collider events. In particular, we use two generative artificial intelligence tools, continuous normalizing flows and diffusion models, to create a set of reconstructed particle-flow objects conditioned on truth-level particles from CMS Open Simulations. While previous work focused on jets, our updated methods now can accommodate all particle-flow objects in an event along with particle-level attributes like particle type and production vertex coordinates. This approach is fully automated, entirely written in Python, and GPU-compatible. Using a variety of physics processes at the LHC, we show that the extended arnassus is able to generalize beyond the training dataset and outperforms the standard, public tool elphes.

Dreyer, Etienne [Weizmann Institute of Science, Re↗

Search for a Sub-eV Sterile Neutrino Using Daya Bay’s Full Dataset

This Letter presents results of a search for the mixing of a sub-eV sterile neutrino with three active neutrinos based on the full data sample of the Daya Bay Reactor Neutrino Experiment, collected during 3158 days of detector operation, which contains 5.55 × 106 reactor $\overline{v}$ e candidates identified as inverse beta-decay interactions followed by neutron capture on gadolinium. The analysis benefits from a doubling of the statistics of our previous result and from improvements of several important systematic uncertainties. No significant oscillation due to mixing of a sub-eV sterile neutrino with active neutrinos was found. Exclusion limits are set by both Feldman-Cousins and CLs methods. Light sterile neutrino mixing with sin 2⁡ 2⁢θ 14 ≳ 0.01 can be excluded at 95% confidence level in the region of 0.01 eV 2 ≲ |Δ⁢$m$$^{2}_{41}$| ≲ 0.1 eV 2 . This result represents the world-leading constraints in the region of 2 × 10 –4 eV 2 ≲ |Δ⁢$m$$^{2}_{41}$| ≲ 0.2 eV 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hemiptera phylogenomic resources: Tree‐based orthology prediction and conserved exon identification

Abstract High‐throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1‐to‐1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1‐to‐1 orthogroups based on tree‐based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree–species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci , a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

Owen, Christopher L.↗

A Wave Energy Converter Design Load Case Study

This article presents an example by which design loads for a wave energy converter (WEC) might be estimated through the various stages of the WEC design process. Unlike previous studies, this study considers structural loads, for which, an accurate assessment is crucial to the optimization and survival of a WEC. Three levels of computational fidelity are considered. The first set of design load approximations are made using a potential flow frequency-domain boundary-element method with generalized body modes. The second set of design load approximations are made using a modified version of the linear-based time-domain code WEC-Sim. The final set of design load simulations are realized using computational fluid dynamics coupled with finite element analysis to evaluate the WEC’s loads in response to both regular and focused waves. This study demonstrates an efficient framework for evaluating loads through each of the design stages. In comparison with experimental and high-fidelity simulation results, the linear-based methods can roughly approximate the design loads and the sea states at which they occur. The high-fidelity simulations for regular wave responses correspond well with experimental data and appear to provide reliable design load data. The high-fidelity simulations of focused waves, however, result in highly nonlinear interactions that are not predicted by the linear-based most-likely extreme response design load method.

16 TIDAL AND WAVE POWER↗