Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Normal Distribution”

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 721 records · Page 40

Wasserstein normalized autoencoder for anomaly detection

A novel anomaly detection algorithm is presented. The Wasserstein normalized autoencoder (WNAE) is a normalized probabilistic model that minimizes the Wasserstein distance between the learned probability distribution—a Boltzmann distribution where the energy is the reconstruction error of the autoencoder (AE)—and the distribution of the training data. This algorithm has been developed and applied to the identification of semivisible jets—conical sprays of visible standard model (SM) particles and invisible dark matter states—with the CMS experiment at the CERN LHC. Trained on jets of particles from simulated SM processes, the WNAE is shown to learn the probability distribution of the input data in a fully unsupervised fashion, such that it effectively identifies new physics jets as anomalies. The model exhibits stable, convergent training and recovers strong classification performance for a wide range of signals against the selected background process, for which a standard AE fails because of outlier reconstruction. In addition, the model improves upon standard normalized autoencoders while remaining fully agnostic to the signal. The WNAE directly tackles the problem of outlier reconstruction, a common failure mode of autoencoders in anomaly detection tasks.

Hayrapetyan, Aram [Yerevan Phys. Inst.]↗

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗

Scaling laws for the additive manufacturing

The evaluation of simple thermal model of selective laser melting (SLM) process shows that the temperature distribution in the sample is characterized by two dimensionless parameters: normalized enthalpy and the ratio of dwell time to the thermal diffusion time. Here we demonstrated that the melt depth data taken for different machines for different materials collapsed in one curve, making possible to rescale the optimal processing parameters between the different materials and machines. The melt pool depth, width and the length are the universal function of these two parameters. Within the operational range of parameters for SLM these functions can be interpolated by the simple algebraic expressions given the possibility to reduce the calculation of the melt pool characteristic to spreadsheet model.

36 MATERIALS SCIENCE↗

Predicting properties of gas and solid streams by intrinsic kinetics of fast pyrolysis of wood

Pyrolysis has the potential to create a biocrude oil from biomass sources that can be used as fuel or as feedstock for subsequent upgrading to hydrocarbon fuels or other chemicals. The product distribution/composition, however, is linked to the biomass source. This work investigates the products formed from pyrolysis of woody biomass with a previously developed chemical kinetics model. Different woody feedstocks reported in prior literature are placed on a common basis (moisture, ash, fixed carbon free) and normalized by initial elemental composition through ultimate analysis. Observed product distributions over the full devolatilization range are explored, reconstructed by the model, and verified with independent experimental data collected with a microwave-assisted pyrolysis system. These trends include production of permanent gas (CO, CO 2 ), char, and condensable (oil, water) species. Elementary compositions of these streams are also investigated. As a result, close agreement between literature data, model predictions, and independent experimental data indicate that the proposed model/method is able to predict the ideal distribution from fast pyrolysis given reaction temperature, residence time, and feedstock composition.

09 BIOMASS FUELS↗

Large Scale Study of Ligand–Protein Relative Binding Free Energy Calculations: Actionable Predictions from Statistically Robust Protocols

The accurate and reliable prediction of protein–ligand binding affinities can play a central role in the drug discovery process as well as in personalized medicine. Of considerable importance during lead optimization are the alchemical free energy methods that furnish an estimation of relative binding free energies (RBFE) of similar molecules. Recent advances in these methods have increased their speed, accuracy, and precision. This is evident from the increasing number of retrospective as well as prospective studies employing them. However, such methods still have limited applicability in real-world scenarios due to a number of important yet unresolved issues. Here, we report the findings from a large data set comprising over 500 ligand transformations spanning over 300 ligands binding to a diverse set of 14 different protein targets which furnish statistically robust results on the accuracy, precision, and reproducibility of RBFE calculations. We use ensemble-based methods which are the only way to provide reliable uncertainty quantification given that the underlying molecular dynamics is chaotic. These are implemented using TIES (Thermodynamic Integration with Enhanced Sampling). Results achieve chemical accuracy in all cases. Ensemble simulations also furnish information on the statistical distributions of the free energy calculations which exhibit non-normal behavior. We find that the “enhanced sampling” method known as replica exchange with solute tempering degrades RBFE predictions. We also report definitively on numerous associated alchemical factors including the choice of ligand charge method, flexibility in ligand structure, and the size of the alchemical region including the number of atoms involved in transforming one ligand into another. Our findings provide a key set of recommendations that should be adopted for the reliable application of RBFE methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crowd cluster data in the USA for analysis of human response to COVID-19 events and policies

We provide data on daily social contact intensity of clusters of people at different types of Points of Interest (POI) by zip code in Florida and California. This data is obtained by aggregating fine-scaled details of interactions of people at the spatial resolution of 10 m, which is then normalized as a social contact index. We also provide the distribution of cluster sizes and average time spent in a cluster by POI type. This data will help researchers perform fine-scaled, privacy-preserving analysis of human interaction patterns to understand the drivers of the COVID-19 epidemic spread and mitigation. Current mobility datasets either provide coarse-level metrics of social distancing, such as radius of gyration at the county or province level, or traffic at a finer scale, neither of which is a direct measure of contacts between people. We use anonymized, de-identified, and privacy-enhanced location-based services (LBS) data from opted-in cell phone apps, suitably reweighted to correct for geographic heterogeneities, and identify clusters of people at non-sensitive public areas to estimate fine-scaled contacts.

60 APPLIED LIFE SCIENCES↗

Theoretical analysis of the background intensity distribution in X-ray Birefringence Imaging using synchrotron bending-magnet radiation

In the recently developed technique of X-ray Birefringence Imaging, molecular orientational order in anisotropic materials is studied by exploiting the birefringence of linearly polarized X-rays with energy close to an absorption edge of an element in the material. In the experimental setup, a vertically deflecting high-resolution double-crystal monochromator is used upstream from the sample to select the appropriate photon energy, and a horizontally deflecting X-ray polarization analyzer, consisting of a perfect single crystal with a Bragg reflection at Bragg angle of approximately 45°, is placed downstream from the sample to measure the resulting rotation of the X-ray polarization. However, if the experiment is performed on a synchrotron bending-magnet beamline, then the elliptical polarization of the X-rays out of the electron orbit plane affects the shape of the output beam. Also, because the monochromator introduces a correlation between vertical position and photon energy to the X-ray beam, the polarization analyzer does not select the entire beam, but instead selects a diagonal stripe, the slope of which depends on the Bragg angles of the monochromator and the polarization analyzer. In the present work, the final background intensity distribution is calculated analytically because the phase space sampling methods normally used in ray traces are too inefficient for this setup. X-ray Birefringence Imaging data measured at the Diamond Light Source beamline B16 agree well with the theory developed here.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A forward modeling approach to analyzing galaxy clustering with S IM BIG

We present cosmological constraints from a simulation-based inference (SBI) analysis of galaxy clustering from the SimBIG forward modeling framework. SimBIG leverages the predictive power of high-fidelity simulations and provides an inference framework that can extract cosmological information on small nonlinear scales. In this work, we apply SimBIG to the Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy sample and analyze the power spectrum, P ℓ (k), to k max = 0.5 h/Mpc. We construct 20,000 simulated galaxy samples using our forward model, which is based on 2,000 high-resolution Quijote N -body simulations and includes detailed survey realism for a more complete treatment of observational systematics. We then conduct SBI by training normalizing flows using the simulated samples and infer the posterior distribution of ΛCDM cosmological parameters: Ω m , Ω b , h, n s , σ 8 . We derive significant constraints on Ω m and σ 8 , which are consistent with previous works. Our constraint on σ 8 is 27% more precise than standard P ℓ analyses because we exploit additional cosmological information on nonlinear scales beyond the limit of current analytic models, k > 0.25 h/Mpc. This improvement is equivalent to the statistical gain expected from a standard P ℓ analysis of galaxy sample ~ 60% larger than CMASS. While we focus on P ℓ in this work for validation and comparison to the literature, SimBIG provides a framework for analyzing galaxy clustering using any summary statistic. We expect further improvements on cosmological constraints from subsequent SimBIG analyses of summary statistics beyond P ℓ .

79 ASTRONOMY AND ASTROPHYSICS↗

Characteristics of He ii proximity profiles

The proximity profile in the spectra of z≈3 quasars, where fluxes extend blueward of the He ii Lyα wavelength 304 (1+z) Å, is one of the most important spectral features in the study of the intergalactic medium (IGM). Based on the Hubble Space Telescope spectra of 24 He ii quasars, we find that the majority of them display a proximity profile, corresponding to an ionization radius as large as 20 Mpc in the source's rest frame. In comparison with those in the H i spectra of the quasars at z ≈ 6, the He ii proximity effect is more prominent and is observed over a considerably longer period of reionization. The He ii proximity zone sizes decrease at higher redshifts, particularly at z>3.3. This trend is similar to that for H i, signaling an onset of He ii reionization at z≳4. For quasar SDSS1253+6817 (z = 3.48), the He ii absorption trough displays a gradual decline and serves as a good case for modeling the He ii reionization. To model such a broad profile requires a quasar radiation field whose energy distribution between 4 and 1 Rydberg is considerably harder than normally assumed. The UV continuum of this quasar is indeed exceptionally steep, and the He ii ionization level in the quasar vicinity is higher than the average level in the IGM. These results are evidence that a very hard EUV continuum from this quasar produces a large ionized zone around it. Distinct exceptions are the two brightest He ii quasars at z ≈ 2.8, for which no significant proximity profile is present, probably implying that they are very young.

79 ASTRONOMY AND ASTROPHYSICS↗

Plasma-arc lamp high heat flux cycling exposure of neutron irradiated tungsten materials

Thick plate, unalloyed W was neutron irradiated in the High Flux Isotope Reactor (HFIR) at 550 °C to a fast fluence of 1.24 × 10 25 n m -2 E > 0.1 MeV (~0.24 dpa). Unirradiated and irradiated specimens of the material were high heat flux (HHF) tested in the Plasma Arc Lamp (PAL) facility. The PAL uses a high-power photon source to provide a broad and even heat distribution on the sample surface. To simulate on/off cycling of normal operating plasma, the samples were exposed to approximately 800 cycles at 4.73 MW m -2 absorbed heat flux (incident heat fluxes of 10.95 MW m -2 ). Additionally, after PAL exposure, slight changes were observed on the surfaces of the samples with SEM. The samples showed some annealing in the near surface polished region, but they were all below the damage threshold for cracking or other destructive features. The PAL has a large parameter space for future testing. The use of the HFIR and PAL to sequentially expose neutron irradiated samples to HHF will be a powerful tool for understanding materials behavior in a fusion-like environment.

36 MATERIALS SCIENCE↗

γ decay from the quasicontinuum of 197,198 Au

Here, the average electromagnetic dipole response of levels in the quasicontinuum of 197,198 Au has been measured using ( 3 He, 3 He') and (d,p) reactions. The extracted γ-ray strength functions have been normalized according to three model assumptions for the nuclear spin distribution. An enhancement in the energy region E γ = 3.0–6.5 MeV is observed for both isotopes. The E1 component of such excess of strength is studied in detail for 198 Au and is interpreted as the pygmy dipole resonance with an energy centroid of 5.9(1) MeV and exhausts about 1% of the total integrated strength. The pygmy dipole resonance is shown to have a significant impact on the calculated 197 Au(n,γ) 198 Au cross section.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Constraints on dark photon dark matter from Lyman- α forest simulations and an ultrahigh signal-to-noise quasar spectrum

The ultralight dark photon is a well-motivated, hypothetical dark matter candidate. In a dilute plasma, they can resonantly convert into photons, and heat up the intergalactic medium between galaxies. In this work, we explore the dark photon dark matter parameter space by comparing synthetic Lyman- α forest data from cosmological hydrodynamical simulations to observational data from VLT/UVES of the quasar HE0940-1050 ( z em = 3.09 ). We use a novel flux normalization technique that targets underdense gas, reshaping the flux probability distribution. Not only do we place robust constraints on the kinetic mixing parameter of dark photon dark matter, but notably our findings suggest that this model can still reconcile simulated and observed Doppler parameter distributions of z ∼ 0 Lyman- α lines, as seen by HST/COS. This work opens new pathways for the use of the Lyman- α forest to explore new physics, and can be extended to other scenarios such as primordial black hole evaporation, dark matter decay, and annihilation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Wake loss and energy spread factor of the LEReC Booster cavity caused by short range wake field

LEReC project uses a DC photoemission gun with multi-alkali (CsK 2 Sb or NaK 2 Sb) cathode [1]. To get 24 mm “flat-top” distribution, 32 Gaussian laser bunches with 0.6 mm rms length are stacked together with 0.75 mm distance [2]. In this case one cannot simply use a 1 cm rms length Gaussian/step/delta bunch for short range wake field simulation since a 0.6 mm bunch contains frequency much higher than the 1 cm bunch. A short range wake field simulation was done using CST Particle Studio™ with 0.6 mm rms Gaussian bunch at the speed of light, and this result was compared with the result for 1 cm rms Gaussian bunch in Figure 1, from where one notice that the wake potential for the 0.6 mm bunch is ~10 times higher than that of the 1 cm bunch. The wake potential of the 0.6 mm bunch, as well as the charge distribution, was then “shift and stack” every 0.75 mm, the normalized results are shown in Figure 2. The wake loss factor (WLF) is the integration of the product of wake potential and normalized bunch charge, and the energy spread factor (ESF) is the rms deviation from the average energy loss. It is calculated by summing the weighted squares of the differences and taking the square root of the sum. These two factors were then divided by β 2 for 1.6 MV beam energy. The wake loss factor is at 0.86 V/pC and energy spread factor is at 0.54 V/pC rms. With 100 pC electron bunch, the energy spread inter-bunch is 54 V rms.

43 PARTICLE ACCELERATORS↗

Subchannel codes: CTF and VIPRE-01

In this chapter, two state-of-the-art and well-used subchannel codes are presented: CTF and VIPRE-01. As part of the Consortium for Advanced Simulation of Light Water Reactors, CTF is being used for single- and two-phase flow analysis under light-water reactor normal and anticipated transient operating conditions. Accurate determination of flow distribution, pressure drop, and void content is crucial for predicting margins to thermal crisis and ensuring more efficient nuclear power plant performance. In preparation for the intended applications, CTF is extensively validated against data from experimental facilities. Meanwhile, the licensed commercial code VIPRE-01 is used to generate a baseline set of simulations for the targeted tests and the solution parameters are compared to CTF results. The closure models in CTF that are in the greatest need of improvement are highlighted and discussed.

Zhao, Xingang↗

Bayesian Monte Carlo Evaluation Framework for Imperfect Data and Models [Abstract]

Nuclear data evaluation methods conventionally make the following assumptions: prior and posterior probability distribution functions (PDFs) of all model parameters and data are normal (Gaussian); the linear approximation is sufficiently accurate for minimization of a cost function (even for non-linear models); and that both the model (of, e.g., neutron cross section) and experimental data (including their covariance data) are perfect. These assumptions are inherent to the well-known generalized linear least squares (GLLS) minimization method commonly used for evaluations of resolved resonance region (RRR) neutron cross sections. However, these assumptions are often not justified due to the presence of non-normal PDFs, non-linear models (e.g. R -matrix formalism), and inherent imperfections in data and models (e.g. discrepant data sets, discrepancies between the previous evaluation and newly measured data, or imperfect covariance data). We remove the said assumptions in a mathematical framework of Bayes’ theorem, and implement it using the Metropolis-Hastings Monte Carlo method. Parameters of a new kind are introduced to parameterize inherent imperfections, e.g. , any discrepancies between the theoretical model and measured data. These new parameters enable evaluators to quantify their expert judgement about any discrepancies or imperfections in a reproducible manner. We demonstrate the framework with an ongoing evaluation of 233 U in the eV region using the ENDF-B/VIII library and transmission data measured by Guber, et al. , and compare the posterior parameters to those obtained by conventional evaluation methods.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fe Coated Optical Fiber for Distributed Corrosion Monitoring in Soil and Aqueous Environment

Natural gas pipeline corrosion represents a substantial cost and safety concern during normal operations. Effective monitoring of corrosion is important to detect and mitigate pipeline risks before corrosion-related catastrophic events happen. Here, we describe the use of Fe-coated optical fiber sensors (OFS) for distributed corrosion monitoring where Fe acts as a corrosion proxy. By using an optical backscatterring reflectometer (OBR), corrosion was monitored based on the increase in backscattered intensity amplitude of the light being passed as Fe undergoes corrosion. The Fe-coated OFSs were prepared with a film thickness between 25-225 nm by an electroless plating approach and corroded by a CO2-saturated acidic electrolyte. The corrosion rate was approximately 2.5 mm/yr for Fe with a film thickness of 30 nm and increased with film thickness. Backscattering-based corrosion rates were supported by visible light transmission measurements, which have been previously demonstrated. Additionally, a correlation between the intermittent transmission and residual Fe film thickness during corrosion of Fe was established. Corrosion was measurable out to > 100 m which is a significant improvement over our previous work showing corrosion sensing at < 10 m. The corrosion sensor was also tested in soil at >1 ft depth. Fe-coated OFS was loaded into a Draka cable which provides mechanical support during fiber deployment and installation into either acidified soil, top-soil, or sandy soil.

aqueous environment↗

HfZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-zirconium (Hf-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions.

36 MATERIALS SCIENCE↗

HfTa_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys hafnium-tantalum (Hf-Ta). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Hf and Ta. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect BCC lattice sites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,040 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗