Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “methods: statistical”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Characterizing the effect of hypersonic boundary layer turbulence on antenna performance: A computational approach

The degradation of antenna performance during hypersonic re-entry is a well known phenomenon that can lead to complete radio blackout. Recent additions to the Empire code establish it as a tool for the study and analysis of the problem. Coupling to the Sandia Parallel Aerodynamics and Reentry Code (SPARC) enables the electromagnetic analysis of realistic re-entry plasma profiles. The geometric flexibility afforded by both Empire and SPARC allow the consideration of arbitrary vehicle and antenna configurations. We have used this tool to study antenna performance during re-entry when the boundary layer becomes turbulent. A concise description of line-of-sight transmissions, which employs advanced statistical methods, was developed. New insights into the low altitude reflectometer readings of RAM-C2 are offered. Techniques for the reconstruction of the re-entry plasma profile from reflectometer data were explored.

42 ENGINEERING↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Feasibility of an Accelerometer-Based Structural Health Monitoring System for the LANL Blast Tube

A modeling- and simulation-based study was conducted on the feasibility of implementing an accelerometer-based SHM system on the Los Alamos National Laboratory blast tube. A blast tube experiment was modeled using the Abaqus explicit finite element solver. A custom user subroutine was written to apply test-like pressure loading to the inside surface of the blast tube. The subroutine applies analytically defined pressure loads derived from tracer output taken from a Compressible Flow Computational Fluid Dynamics Solver model of the blast tube. Five unique versions of the model were created: an undamaged reference model at 65°F was used as the baseline and compared to equivalent models at 10°F and 100°F. These three models were compared to models with small damage at the reference temperature. The two types of damage considered were a radial (circumferential) crack in the main tube body and a longitudinal crack in the supports. Acceleration outputs were extracted from accelerometer bodies included in the model and were post processed using a variety of standard SHM techniques. Different potential features signaling failure were extracted and compared using statistical methods in the time and frequency domains. A method was identified that clearly shows that differences in structural response resulting from the modeled damage can be differentiated from the structural response resulting from changing environmental conditions. However, the amount of damage applied to create observable differences in the accelerometer data was so large that simpler methods of damage detection would be more cost effective in locating damage.

42 ENGINEERING↗
DOE OSTI2025

Advancing Neutrino Simulation Modeling with MARLEY: Insights from the NNSA-MSIIP Internship

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. However, there is no experimental data to validate the current event generator MARLEY (Model of Argon Reaction Low Energy Yields) that can model supernova neutrinos. To validate the model, I developed a new muon capture feature within the MARLEY simulation framework by coding key functions in C++, Python, and ROOT. I generated and analyzed one million simulated muon capture events and comparing the results to experimental data. To improve the model, I worked on optimizing the model’s parameters to improve its precision using statistical methods.

Wong, Baker [Fermilab]↗
DOE OSTI2019

Pore Formation in the Laminar Structure of Wear-Resistant Modified Powder Coatings

The open porosity of wear-resistant modified coatings is investigated, and their macrostructure is analyzed. A statistical method is proposed for describing the formation of porosity in the laminar structure of a powder coating. The matrix of transition probabilities is determined for a uniform Markov chain. The porosity of a coating on a smooth substrate is calculated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessing the Contribution of ECa and NDVI in the Delineation of Management Zones in a Vineyard

Precision fertilization implies the need to identify the variability of soil fertility, which is costly and time-consuming. Remotely measured data can be a solution. Using this strategy, a study was conducted, in a vineyard, to delineate different management zones using two indicators: apparent soil electrical conductivity (ECa) and normalized difference vegetation index (NDVI). To understand the contribution of each indicator, three scenarios were used for zone definition: (1) using only NDVI, (2) only ECa, or (3) using a combination of the two. Then the differences in soil fertility between these zones were assessed using simple statistical methods. The results indicate that the most beneficial strategy is the combined use of the two indicators, as it allowed the definition of three distinct zones regarding important soil variables and crop nutrients, such as soil total nitrogen, Mg2+ cation, exchange acidity, and effective cation exchange capacity, and some relevant cation ratios. This strategy also allowed the identification of an ionic unbalance in the soil chemistry, due to an excess of Mg2+, that was harming crop health, as reported by NDVI. This also impacted ECa and NDVI relationship, which was negative in this study. Overall, the results demonstrate the advantages of using remotely sensed data, mainly more than one type of sensing data, and suggest a high potential for differential crop fertilization and soil management in the study area.

Esteves, Catarina (ORCID:0000000327205906)↗
DOE OSTI2021

Spatial Interpolation of Air Pollutant and Meteorological Variables in Central Amazonia

The Amazon Rainforest is highlighted by the global community both for its extensive vegetation cover that constantly suffers the effects of anthropic action and for its substantial biodiversity. This dataset presents data of meteorological variables from the Amazon Rainforest region with a spatial resolution of 0.001° in latitude and longitude, resulting from an interpolation process. The original data were obtained from the GoAmazon 2014/5 project, in the Atmospheric Radiation Measurement (ARM) repository, and then processed through mathematical and statistical methods. The dataset presented here can be used in experiments in the field of Data Science, such as training models for predicting climate variables or modeling the distribution of species.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Information Length Concept Applied to Plasma Turbulence

A methodology to study statistical properties of anomalous transport in fusion plasma is investigated. Three time traces generated by the full-f gyrokinetic code GKNET are analyzed for this purpose. The time traces consist of heat flux as a function of the radial position, which is studied in a novel manner using statistical methods. The simulation data exhibit transport processes with both medium and long correlation length along the radius. A typical example of a phenomenon with long correlation length is avalanches. In order to investigate the evolution of the turbulent state, two basic configurations are studied, one flux-driven and one gradient-driven with decaying turbulence. The information length concept in tandem with Boltzmann–Gibbs and Tsallis entropy is used in the investigation. It is found that the dynamical states in both flux-driven and gradient-driven cases are surprisingly similar, but the Tsallis entropy reveals differences between them. This indicates that the types of probability distribution function are nevertheless quite different since the higher moments are significantly different.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

THE THIRD FERMI GBM GAMMA-RAY BURST CATALOG: THE FIRST SIX YEARS

Since its launch in 2008, the Fermi Gamma-ray Burst Monitor (GBM) has triggered and located on average approximately two γ -ray bursts (GRBs) every three days. Here, we present the third of a series of catalogs of GRBs detected by GBM, extending the second catalog by two more years through the middle of 2014 July. The resulting list includes 1405 triggers identified as GRBs. The intention of the GBM GRB catalog is to provide information to the community on the most important observables of the GBM-detected GRBs. For each GRB, the location and main characteristics of the prompt emission, the duration, peak flux, and fluence are derived. The latter two quantities are calculated for the 50–300 keV energy band where the maximum energy release of GRBs in the instrument reference system is observed, and also for a broader energy band from 10 to 1000 keV, exploiting the full energy range of GBM's low-energy [Nai[Tl)] detectors. Using statistical methods to assess clustering, we find that the hardness and duration of GRBs are better fit by a two-component model with short-hard and long-soft bursts than by a model with three components. Furthermore, information is provided on the settings and modifications of the triggering criteria and exceptional operational conditions during years five and six in the mission. This third catalog is an official product of the Fermi GBM science team, and the data files containing the complete results are available from the High-Energy Astrophysics Science Archive Research Center.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A NEW REFERENCE CHEMICAL COMPOSITION FOR TMC-1

Recent detections of complex organic molecules in dark clouds have rekindled interest in the astrochemical modeling of these environments. Because of its relative closeness and rich molecular complexity, TMC-1 has been extensively observed to study the chemical processes taking place in dark clouds. We use local thermodynamical equilibrium radiative transfer modeling coupled with a Bayesian statistical method which takes into account outliers to analyze the data from the Nobeyama spectral survey of TMC-1 between 8 and 50 GHz. We compute the abundance relative to molecular hydrogen of 57 molecules, including 19 isotopologues in TMC-1 along with their associated uncertainty. The new results are in general agreement with previous abundance determination from Ohishi and Kaifu and the values reported in the review from Agúndez and Wakelam. However, in some cases, large opacity and low signal to noise effects allow only upper or lower limits to be derived, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Sub-band Image Reconstruction Using Differential Chromatic Refraction

Refraction by the atmosphere causes the positions of sources to depend on the airmass through which an observation was taken. This shift is dependent on the underlying spectral energy of the source and the filter or bandpass through which it is observed. Wavelength-dependent refraction within a single passband is often referred to as differential chromatic refraction (DCR). With a new generation of astronomical surveys undertaking repeated observations of the same part of the sky over a range of different airmasses and parallactic angles, DCR should be a detectable and measurable astrometric signal. In this paper we introduce a novel procedure that takes this astrometric signal and uses it to infer the underlying spectral energy distribution of a source; we solve for multiple latent images at specific wavelengths via a generalized deconvolution procedure built on robust statistics. We demonstrate the utility of such an approach for estimating a partially deconvolved image, at higher spectral resolution than the input images, for surveys such as the Large Synoptic Survey Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Photometric and Kinematic Study of the Open Clusters SAI 44 and SAI 45

We carry out a detailed photometric and kinematic study of the poorly studied sparse open clusters SAI 44 and SAI 45 using ground-based BVR {sub c} I {sub c} data supplemented by archival data from Gaia eDR3 and Pan-STARRS. The stellar memberships are determined using a statistical method based on Gaia eDR3 kinematic data, and we found 204 members in SAI 44 while only 74 members are identified in SAI 45. The average distances to SAI 44 and SAI 45 are calculated to be 3670 ± 184 and 1668 ± 47 pc. The logarithmic age of the clusters are determined to be 8.82 ± 0.10 and 9.07 ± 0.10 yr for SAI 44 and SAI 45, respectively. The color–magnitude diagram of SAI 45 hosts an extended main-sequence turnoff (eMSTO). The apparent age spread is found to be similar to the apparent age spread predicted on the basis of the age spread and cluster age relation predicted by rotation models. This indicates that eMSTO is a stellar evolution rather than star formation phenomenon in SAI 45. We conclude that eMSTO in SAI 45 is mainly caused by the different rotation rates of stars as the SYCLIST synthetic population with different rotation rates was able to reproduce the observed eMSTO, and stars in the red part of the eMSTO were preferentially concentrated in the inner region, which again hints at different rotations being the reason for the extension in the upper MS. This finding supports the theory attributing the origin of eMSTO to the different rotations of eMSTO stars. The mass function slopes are obtained as -2.24 ± 0.66 and -2.58 ± 3.20 in the mass rages 2.426–0.990 M{sub ⊙} and 2.167–1.202 M{sub ⊙} for SAI 44 and SAI 45, respectively. SAI 44 exhibits the signature of mass segregation while we found weak evidence of mass segregation in SAI 45 possibly due to tidal stripping. The dynamical relaxation times of these clusters indicate that both clusters are in a dynamically relaxed state. Using the AD-diagram method, the apex coordinates are found to be ([[Formula]], -[[Formula]]) for SAI 44 and (-[[Formula]], -[[Formula]]) for SAI 45. The average space velocity components of the clusters SAI 44 and SAI 45 are calculated in units of km s{sup -1} as (-15.14 ± 3.90, -19.43 ± 4.41, -20.85 ± 4.57) and (28.13 ± 5.30, -9.78 ± 3.13, -19.59 ± 4.43), respectively.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

From Neutron Star Observables to the Equation of State. II. Bayesian Inference of Equation of State Pressures

One of the key goals of observing neutron stars is to infer the equation of state (EoS) of the cold, ultradense matter in their interiors. Here, we present a Bayesian statistical method of inferring the pressures at five fixed densities, from a sample of mock neutron star masses and radii. We show that while five polytropic segments are needed for maximum flexibility in the absence of any prior knowledge of the EoS, regularizers are also necessary to ensure that simple underlying EoS are not over-parameterized. For ideal data with small measurement uncertainties, we show that the pressure at roughly twice the nuclear saturation density, ρ{sub sat}, can be inferred to within 0.3 dex for many realizations of potential sources of uncertainties. The pressures of more complicated EoS with significant phase transitions can also be inferred to within ∼30%. We also find that marginalizing the multi-dimensional parameter space of pressure to infer a mass–radius relation can lead to biases of nearly 1 km in radius, toward larger radii. Using the full, five-dimensional posterior likelihoods avoids this bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Determination of Dark Matter Halo Mass from Dynamics of Satellite Galaxies

We propose that the mass of a dark matter halo can be inferred from the dynamical status of its satellite galaxies. Using nine dark matter simulations of halos like the Milky Way (MW), we find that the present-day substructures in each halo follow a characteristic distribution in the phase space of orbital binding energy and angular momentum, and that this distribution is similar from halo to halo, but has an intrinsic dependence on the halo formation history. We build this distribution directly from the simulations for a specific halo and extend the result to halos of similar formation history but different masses by scaling. The mass of an observed halo can then be estimated by maximizing the likelihood in comparing the measured kinematic parameters of its satellite galaxies with these distributions. We test the validity and accuracy of this method with mock samples taken from the simulations. Using the positions, radial velocities, and proper motions of nine tracers and assuming observational uncertainties comparable to those of MW satellite galaxies, we find that the halo mass can be recovered to within $\sim 40 \% $. The accuracy can be improved to within ~25% if 30 tracers are used. Yet, the dependence of the phase-space distribution on the halo formation history sets a minimum uncertainty of $\sim 20 \% $ that cannot be reduced by using more tracers. We believe that this minimum uncertainty also applies to any mass determination for a halo when the phase-space information of other kinematic tracers is used.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Quantitative Constraints on the Reionization History from the IGM Damping Wing Signature in Two Quasars at z > 7

During reionization, neutral hydrogen in the intergalactic medium (IGM) imprints a damping wing absorption feature on the spectrum of high-redshift quasars. A detection of this signature provides compelling evidence for a significantly neutral universe, and enables measurements of the hydrogen neutral fraction x HI (z) at that epoch. Obtaining reliable quantitative constraints from this technique, however, is challenging due to stochasticity induced by the patchy inside-out topology of reionization, degeneracies with quasar lifetime, and the unknown unabsorbed quasar spectrum close to rest-frame Lyα. Here, we combine a large-volume semi-numerical simulation of reionization topology with 1D radiative transfer through high-resolution hydrodynamical simulations of the high-redshift universe to construct models of quasar transmission spectra during reionization. Our state-of-the-art approach captures the distribution of damping wing strengths in biased quasar halos that should have reionized earlier, as well as the erosion of neutral gas in the quasar environment caused by its own ionizing radiation. Combining this detailed model with our new technique for predicting the quasar continuum and its associated uncertainty, we introduce a Bayesian statistical method to jointly constrain the neutral fraction of the universe and the quasar lifetime from individual quasar spectra. We apply this methodology to the spectra of the two quasars with the highest redshifts known, ULAS J1120+0641 and ULAS J1342+0928, and measure volume-averaged neutral fractions {x HI }(Z = 7.09) = 0.48 $+0.26\atop{-0.26}$and {x HI }(Z = 7.54) = 0.60 $+0.20\atop{-0.23}$ (posterior medians and 68% credible intervals) when marginalized over quasar lifetimes of 10 3 ≤ t q ≤ 10 8 yr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Transport of Protostellar Cosmic Rays in Turbulent Dense Cores

We report recent studies have suggested that low-energy cosmic rays (CRs) may be accelerated inside molecular clouds by the shocks associated with star formation. We use a Monte Carlo transport code to model the propagation of CRs accelerated by protostellar accretion shocks through protostellar cores. We calculate the CR attenuation and energy losses and compute the resulting flux and ionization rate as a function of both radial distance from the protostar and angular position. We show that protostellar cores have nonuniform CR fluxes that produce a broad range of CR ionization rates, with the maximum value being up to two orders of magnitude higher than the radial average at a given distance. In particular, the CR flux is focused in the direction of the outflow cavity, creating a "flashlight" effect and allowing CRs to leak out of the core. The radially averaged ionization rates are less than the measured value for the Milky Way of ζ ≈ 10 -16 s -1 ; however, within r ≈ 0.03 pc from the protostar, the maximum ionization rates exceed this value. We show that variation in the protostellar parameters, particularly in the accretion rate, may produce ionization rates that are a couple of orders of magnitude higher or lower than our fiducial values. Finally, we use a statistical method to model unresolved subgrid magnetic turbulence in the core. We show that turbulence modifies the CR spectrum and increases the uniformity of the CR distribution but does not significantly affect the resulting ionization rates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

MicroRNAs as diagnostic and prognostic biomarkers of age-related macular degeneration: advances and limitations

A main cause of vision loss in the elderly is age-related macular degeneration (AMD). Among the cellular, biochemical, and molecular changes linked to this disease, inflammation and angiogenesis appear as being crucial in AMD pathogenesis and progression. There are two forms of the disease: dry AMD, accounting for 80–90% of cases, and wet AMD. The disease usually begins as dry AMD associated with retinal pigment epithelium and photoreceptor degeneration, whereas wet AMD is associated with choroidal neovascularization resulting in severe vision impairment. The new vessels are largely malformed, leading to blood and fluid leakage within the disrupted tissue, which provokes inflammation and scar formation and results in retinal damage and detachment. MicroRNAs are dysregulated in AMD and may facilitate the early detection of the disease and monitoring disease progression. Two recent reviews of microRNAs in AMD had indicated weaknesses or limitations in four earlier investigations. Studies in the last three years have shown considerable progress in overcoming some of these concerns and identifying specific microRNAs as biomarkers for AMD. Further large-scale studies are warranted using appropriate statistical methods to take into account gender and age disparity in the study populations and confounding factors such as smoking status.

60 APPLIED LIFE SCIENCES↗