Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

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 73 records · Page 4

Lighting up the LHC with Dark Matter

We show that simultaneously explaining dark matter and the observed value of the muon’s magnetic dipole moment may lead to yet unexplored photon signals at the LHC. We consider the Minimal Supersymmetric Standard Model with electroweakino masses in the few-to-several hundred GeV range, and opposite sign of the Bino mass parameter with respect to both the Higgsino and Wino mass parameters. In such region of parameter space, the spin-independent elastic scattering cross section of a Bino-like dark matter candidate in direct detection experiment is suppressed by cancellations between different amplitudes, and the observed dark matter relic density can be realized via Bino-Wino co-annihilation. Moreover, the observed value of the muon’s magnetic dipole moment can be explained by Bino and Wino loop contributions. Interestingly, “radiative” decays of Wino-like neutralinos into the lightest neutralino and a photon are enhanced, whereas decays into leptons are suppressed. While these decay patterns weaken the reach of multi-lepton searches at the LHC, the radiative decay opens a new window for probing dark matter at the LHC through the exploration of parameter space regions beyond those currently accessible. To complement the current electroweakino searches, we propose searching for a single (soft) photon plus missing transverse energy, accompanied by a hard initial state radiation jet.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Film Thickness and Composition Effects in Symmetric Ternary Block Copolymer/Homopolymer Blend Films: Domain Spacing and Orientation

We report thickness and homopolymer molar mass dependencies in thin film blends of a compositionally symmetric 75 kg/mol lamellar polystyrene-block-poly(methyl methacrylate) diblock copolymer (PS-b-PMMA) with near-equal volume fractions of PS and PMMA homopolymers on “neutral” surfaces that promote a vertical domain orientation. PS and PMMA homopolymers with equal molar masses of 1, 3, and 22 kg/mol are blended in a 1:1 mass ratio. To efficiently explore this parameter space, we prepared combinatorial samples with gradient homopolymer concentrations using electrospray deposition (ESD) and gradient film thicknesses using flow coating; these samples were subsequently characterized by grazing-incidence small-angle X-ray scattering and scanning electron microscopy, respectively. Our results demonstrate that the added homopolymers increase or reduce domain spacings with respect to the neat (unblended) block copolymer depending on their molar masses, in line with previous reports; a simple heuristic to estimate the domain scaling in ternary blends based on the homopolymer molar mass agrees well with data for the thinnest blend films studied here. More surprisingly, we observe that the measured domain spacing in blends also depends on film thickness, with thicker films exhibiting larger domain spacings than thinner films at the same homopolymer concentration. Furthermore, thickness gradients at fixed homopolymer concentrations reveal a change in lamellae orientation from vertical to horizontal as film thickness is increased beyond the nominal lamellar domain spacing. Coarse-grained molecular dynamics simulations indicate that this change in orientation is induced by homopolymer segregation to film surfaces. Finally, tying these effects together, we hypothesize that the observed dependence of domain spacing on film thickness is a consequence of changes in the vertical homopolymer composition profile in thicker films during lamellae reorientation and the conformational asymmetry between PMMA and PS segments.

36 MATERIALS SCIENCE↗

Influence of transient baroclinic eddies on planetary-scale waves

Low-order models with a mean zonal flow forced by zonally symmetric heating and one planetary-scale wave forced by topography have multiple flow equilibria. Two stable equilibria are characterized by strong zonal flow with low wave amplitude (high index circulation) and weak zonal flow with a high wave amplitude (low index circulation) fixed in phase with the topography. These two states presumably represent normal zonal circulation and a blocking configuration, respectively. When a shorter, baroclinically unstable wave is dded to this low-order model, the planetary-scale wave no longer stays in stable equilibrium states. Instead, the long wave remains in weather regimes or preferred regions in phase space. These regimes are, in general, different from the equilibria of the model with the planetary-scale wave only. Thus, the short unstable wave adds some randomness to planetary-scale circulation and changes its position with respect to the planetary-scale topography. The Reinhold and Pierrehumbert model is further explored here by adding a long-wave in the thermal forcing and exploring wider parameter space. When the symmetric thermal forcing is weak so that the short wave is stable, the planetary scale wave is fixed in amplitude and phase by the asymmetric forcing. However, when the symmetric forcing is increased to an unstable level, the planetary-scale wave becomes less organized. With asymmetric thermal forcing only (i.e., no topography) the long-wave is randomly distributed through phase space. Thermal forcing seems to be less effective in organizing weather regimes than topographic forcing.

Branscome, L. E.↗

Reinforcement learning based automated history matching for improved hydrocarbon production forecast

History matching aims to find a numerical reservoir model that can be used to predict the reservoir performance. An engineer and model calibration (data inversion) method are required to adjust various parameters/properties of the numerical model in order to match the reservoir production history. In this study, we develop deep neural networks within the reinforcement learning framework to achieve automated history matching that will reduce engineers’ efforts, human bias, automatically and intelligently explore the parameter space, and remove the need of large set of labeled training data. To that end, a fast-marching-based reservoir simulator is encapsulated as an environment for the proposed reinforcement learning. The deep neural-network-based learning agent interacts with the reservoir simulator within reinforcement learning framework to achieve the automated history matching. Reinforcement learning techniques, such as discrete Deep Q Network and continuous Deep Deterministic Policy Gradients, are used toth, used to train the learning agents. The continuous actions enable the Deep Deterministic Policy Gradients to explore more states at each iteration in a a learning episode; consequently, a better history matching is achieved using this algorithm as compared to Deep Q Network. For simplified dual-target composite reservoir models, the best history-matching performances of the discrete and continuous learning methods in terms of normalized root mean square errors are 0.0447 and 0.0038, respectively. Furthermore, our study shows that continuous action space achieved by the deep deterministic policy gradient drastically outperforms deep Q network.

42 ENGINEERING↗

Testing Dissipative Magnetosphere Model Light Curves and Spectra with Fermi Pulsars

We explore the emission properties of a dissipative pulsar magnetosphere model introduced by Kalapotharakos et al. comparing its high-energy light curves and spectra, due to curvature radiation, with data collected by the Fermi LAT. The magnetosphere structure is assumed to be near the force-free solution. The accelerating electric field, inside the light cylinder (LC), is assumed to be negligible, while outside the LC it rescales with a finite conductivity (sigma). In our approach we calculate the corresponding high-energy emission by integrating the trajectories of test particles that originate from the stellar surface, taking into account both the accelerating electric field components and the radiation reaction forces. First, we explore the parameter space assuming different value sets for the stellar magnetic field, stellar period, and conductivity. We show that the general properties of the model are in a good agreement with observed emission characteristics of young gamma-ray pulsars, including features of the phase-resolved spectra. Second, we find model parameters that fit each pulsar belonging to a group of eight bright pulsars that have a published phase-resolved spectrum. The sigma values that best describe each of the pulsars in this group show an increase with the spin-down rate (E˙ ) and a decrease with the pulsar age, expected if pair cascades are providing the magnetospheric conductivity. Finally, we explore the limits of our analysis and suggest future directions for improving such models.

acceleration of particles↗

Multi-objective Bayesian optimization of ferroelectric materials with interfacial control for memory and energy storage applications

Optimization of materials’ performance for specific applications often requires balancing multiple aspects of materials’ functionality. Even for the cases where a generative physical model of material behavior is known and reliable, this often requires search over multidimensional function space to identify low-dimensional manifold corresponding to the required Pareto front. In this work, we introduce the multi-objective Bayesian optimization (MOBO) workflow for the ferroelectric/antiferroelectric performance optimization for memory and energy storage applications based on the numerical solution of the Ginzburg–Landau equation with electrochemical or semiconducting boundary conditions. MOBO is a low computational cost optimization tool for expensive multi-objective functions, where we update posterior surrogate Gaussian process models from prior evaluations and then select future evaluations from maximizing an acquisition function. Using the parameters for a prototype bulk antiferroelectric (PbZrO 3 ), we first develop a physics-driven decision tree of target functions from the loop structures. We further develop a physics-driven MOBO architecture to explore multidimensional parameter space and build Pareto-frontiers by maximizing two target functions jointly—energy storage and loss. This approach allows for rapid initial materials and device parameter selection for a given application and can be further expanded toward the active experiment setting. The associated notebooks provide both the tutorial on MOBO and allow us to reproduce the reported analyses and apply them to other systems (https://github.com/arpanbiswas52/MOBO_AFI_Supplements).

36 MATERIALS SCIENCE↗

Sensitivity of two-dimensional model predictions of ozone response to stratospheric aircraft: An update

The Goddard Space Flight Center (GSFC) two-dimensional model of stratospheric photochemistry and dynamics has been used to calculate the O3 response to stratospheric aircraft (high-speed civil transport (HSCT)) emissions. The sensitivity of the model O3 response was examined for systematic variations of five parameters and two reaction rates over a wide range, expanding on calculations by various modeling groups for the NASA High Speed Research Program and the World Meteorological Organization. In all, 448 model runs were required to test the effects of variations in the latitude, altitude, and magnetitude of the aircraft emissions perturbation, the background chlorine levels, the background sulfate aerosol surface area densities, and the rates of two key reactions. No deviation from previous conclusions concerning the response of O3 to HSCTs was found in this more exhaustive exploration of parameter space. Maximum O3 depletions occur for high-altitude, low altitude HSCT perturbations. Small increases in global total O3 can occur for low-altitude, high-altitude injections. Decreasing aerosol surface area densities and background chlorine levels increases the sensitivity of model O3 to the HSCT perturbations. The location of the aircraft emissions is the most important determinant of the model response. Response to the location of the HSCT emissions is not changed qualitatively by changes in background chlorine and aerosol loading. The response is also not very sensitive to changes in the rates of the reactions NO + HO2 yields NO2 + OH and HO2 + O3 yields OH + 2O2 over the limits of their respective uncertainties. Finally, levels of lower stratospheric HO(sub x) generally decrease when the HSCT perturbation is included, even though there are large increases in H2O due to the perturbation.

Considine, David B.↗

Airfoil Computational Fluid Dynamics - 2k shapes, 25 AoA's, 3 Re numbers

This dataset contains aerodynamic quantities - including flow field values (momentum, energy, and vorticity) and summary values (coefficients of lift, drag, and momentum) - for 1,830 airfoil shapes computed using the HAM2D CFD (computational fluid dynamics) model. The airfoil shapes were designed using the separable shape tensor parameterization that encodes two-dimensional shapes as elements of the Grassmann manifold. This data-driven approach learns two independent spaces of parameter from a collection of sample airfoils. The first captures large-scale, linear perturbations, and the second defines small-scale, higher-order perturbations. For this dataset, we used the G2Aero database of over 19,000 airfoil shapes to learn a parameter space that captured a wide array of shape characteristics. We sampled airfoil designs over both parameter spaces to explore the full range of possible shape variations. The aerodynamic quantities for the generated airfoil were obtained using the HAM2D code, which is a finite-volume Reynolds-averaged Navier-Stokes (RANS) flow solver. We employ a fifth-order WENO scheme for spatial reconstruction with Roe's flux difference scheme for inviscid flux and second-order central differencing for viscous flux. A preconditioned GMRES method is applied for implicit integration. The Spalart-Allmaras 1-eq turbulence model is used for the turbulence closure, and the Medida-Baeder 2-eq transition model is applied to account for the effects of laminar turbulent transition. The airfoil grid is generated with a total of 400 points on the airfoil surface, the initial wall-normal spacing of y+ = 1, and an outer boundary located at 300 chord lengths away from the wall. The CFD simulations are performed at a freestream Mach number of 0.1, for or three different Reynolds' numbers (3M, 6M, and 9M), and for 25 angles of attack from -4 deg. to 20 deg. with 1 degree increments. Across all these various parameters, this dataset includes the results from over 250,000 CFD simulations. The simulations were performed using the Bridges-2 system at the Pittsburgh Supercomputing Center in February 2023 as part of the INTEGRATE project funded by the Advanced Research Projects Agency - Energy, in the U.S. Department of Energy. The data was collected, reformatted, and preprocessed for this OEDI submission in July 2023 under the Foundational AI for Wind Energy project funded by the U.S. Department of Energy Wind Energy Technologies Office. This dataset is intended to serve as a benchmark against which new artificial intelligence (AI) or machine learning (ML) tools may be tested. Baseline AI/ML methods for analyzing this dataset have been implemented, and a link to their repository containing those models has been provided. The .h5 data file structure can be found in the GitHub Repository resource under explore_airfoil_2k_data.ipynb.

2k↗

Autotuning PolyBench benchmarks with LLVM Clang/Polly loop optimization pragmas using Bayesian optimization

Here, we develop a ytopt autotuning framework that leverages Bayesian optimization to explore the parameter space search and compare four different supervised learning methods within Bayesian optimization and evaluate their effectiveness. We select six of the most complex PolyBench benchmarks and apply the newly developed LLVM Clang/Polly loop optimization pragmas to the benchmarks to optimize them. We then use the autotuning framework to optimize the pragma parameters to improve their performance. The experimental results show that our autotuning approach outperforms the other compiling methods to provide the smallest execution time for the benchmarks syr2k, 3mm, heat-3d, lu, and covariance with two large datasets in 200 code evaluations for effectively searching the parameter spaces with up to 170,368 different configurations. We find that the Floyd-Warshall benchmark did not benefit from autotuning. To cope with this issue, we provide some compiler option solutions to improve the performance. Then we present loop autotuning without a user's knowledge using a simple mctree autotuning framework to further improve the performance of the Floyd-Warshall benchmark. We also extend the ytopt autotuning framework to tune a deep learning application.

79 ASTRONOMY AND ASTROPHYSICS↗

Search for a new heavy scalar particle decaying into a Higgs boson and a new scalar singlet in final states with one or two light leptons and a pair of τ -leptons with the ATLAS detector

A search for a new heavy scalar particle X decaying into a Standard Model (SM) Higgs boson and a new singlet scalar particle S is presented. The search uses a proton-proton (pp) collision data sample with an integrated luminosity of 140 fb -1 recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV with the ATLAS detector at the Large Hadron Collider. The most sensitive mass parameter space is explored in X mass ranging from 500 to 1500 GeV, with the corresponding S mass in the range 200–500 GeV. The search selects events with two hadronically decaying τ-lepton candidates from H → τ + τ - decays and one or two light leptons (ℓ = e, μ) from S → VV (V = W, Z) decays while the remaining V boson decays hadronically or to neutrinos. A multivariate discriminant based on event kinematics is used to separate the signal from the background. No excess is observed beyond the expected SM background and 95% confidence level upper limits between 72 fb and 542 fb are derived on the cross-section σ(pp → X → SH) assuming the same SM-Higgs boson-like decay branching ratios for the S → VV decay. Upper limits on the visible cross-sections σ(pp → X → SH → WWττ) and σ(pp → X → SH → ZZττ) are also set in the ranges 3–26 fb and 6–33 fb, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring color-octet scalar parameter space in minimal R-symmetric models

In this work we study the collider phenomenology of color-octet scalars (sgluons) in minimal supersymmetric models endowed with a global continuous R symmetry. We systematically catalog the significant decay channels of scalar and pseudoscalar sgluons and identify novel features that are natural in these models. These include decays in nonstandard diboson channels, such as to a gluon and a photon; three-body decays with considerable branching fractions; and long-lived particles with displaced vertex signatures. We also discuss the single and pair production of these particles and show that they can evade existing constraints from the Large Hadron Collider, to varying extents, in large regions of reasonable parameter space. We find, for instance, that a 725 GeV scalar and a 350 GeV or lighter pseudoscalar can still be accommodated in realistic scenarios.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Formation of transient high-β plasmas in a magnetized, weakly collisional regime

We present experimental data providing evidence for the formation of transient (~20 μs) plasmas that are simultaneously weakly magnetized (i.e. Hall magnetization parameter ωτ > 1) and dominated by thermal pressure (i.e. ratio of thermal-to-magnetic pressure β > 1). Particle collisional mean free paths are an appreciable fraction of the overall system size. These plasmas are formed via the head-on merging of two plasmas launched by magnetized coaxial guns. The ratio λ gun = μ 0 I gun /ψ gun of gun current I gun to applied magnetic flux ψ gun is an experimental knob for exploring the parameter space of β and ωτ. These experiments were conducted on the Big Red Ball at the Wisconsin Plasma Physics Laboratory. The transient formation of such plasmas can potentially open up new regimes for the laboratory study of weakly collisional, magnetized, high-β plasma physics; processes relevant to astrophysical objects and phenomena; and novel magnetized plasma targets for magneto-inertial fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grad–Shafranov equilibria via data-free physics informed neural networks

A large number of magnetohydrodynamic (MHD) equilibrium calculations are often required for uncertainty quantification, optimization, and real-time diagnostic information, making MHD equilibrium codes vital to the field of plasma physics. In this paper, we explore a method for solving the Grad–Shafranov equation by using physics-informed neural networks (PINNs). For PINNs, we optimize neural networks by directly minimizing the residual of the partial differential equation as a loss function. We show that PINNs can accurately and effectively solve the Grad–Shafranov equation with several different boundary conditions, making it more flexible than traditional solvers. This method is flexible as it does not require any mesh and basis choice, thereby streamlining the computational process. We also explore the parameter space by varying the size of the model, the learning rate, and boundary conditions to map various tradeoffs such as between reconstruction error and computational speed. Additionally, we introduce a parameterized PINN framework, expanding the input space to include variables such as pressure, aspect ratio, elongation, and triangularity in order to handle a broader range of plasma scenarios within a single network. Parameterized PINNs could be used in future work to solve inverse problems such as shape optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Foundations of machine learning for low-temperature plasmas: methods and case studies

Abstract Machine learning (ML) and artificial intelligence have proven to be an invaluable tool in tackling a vast array of scientific, engineering, and societal problems. The main drivers behind the recent proliferation of ML in practically all aspects of science and technology can be attributed to: (a) improved data acquisition and inexpensive data storage; (b) exponential growth in computing power; and (c) availability of open-source software and resources that have made the use of state-of-the-art ML algorithms widely accessible. The impact of ML on the field of low-temperature plasmas (LTPs) could be particularly significant in the emerging applications that involve plasma treatment of complex interfaces in areas ranging from the manufacture of microelectronics and processing of quantum materials, to the LTP-driven electrification of the chemical industry, and to medicine and biotechnology. This is primarily due to the complex and poorly-understood nature of the plasma-surface interactions in these applications that pose unique challenges to the modeling, diagnostics, and predictive control of LTPs. As the use of ML is becoming more prevalent, it is increasingly paramount for the LTP community to be able to critically analyze and assess the concepts and techniques behind data-driven approaches. To this end, the goal of this paper is to provide a tutorial overview of some of the widely-used ML methods that can be useful, amongst others, for discovering and correlating patterns in the data that may be otherwise impractical to decipher by human intuition alone, for learning multivariable nonlinear data-driven prediction models that are capable of describing the complex behavior of plasma interacting with interfaces, and for guiding the design of experiments to explore the parameter space of plasma-assisted processes in a systematic and resource-efficient manner. We illustrate the utility of various supervised, unsupervised and active learning methods using LTP datasets consisting of commonly-available, information-rich measurements (e.g. optical emission spectra, current–voltage characteristics, scanning electron microscope images, infrared surface temperature measurements, Fourier transform infrared spectra). All the ML demonstrations presented in this paper are carried out using open-source software; the datasets and codes are made publicly available. The FAIR guiding principles for scientific data management and stewardship can accelerate the adoption and development of ML in the LTP community.

Physics↗

Interstellar gas heating by primordial black holes

Interstellar gas heating is a powerful cosmology-independent observable for exploring the parameter space of primordial black holes (PBHs) formed in the early Universe that could constitute part of the dark matter (DM). We provide a detailed analysis of the various aspects for this observable, such as PBH emission mechanisms. Using observational data from the Leo T dwarf galaxy, we constrain the PBH abundance over a broad mass-range, M PBH ~ O(1)M ⊙ – 10 7 M ⊙ , relevant for the recently detected gravitational wave signals from intermediate-mass BHs. Here, we also consider PBH gas heating of systems with bulk relative velocity with respect to the DM, such as Galactic clouds.

79 ASTRONOMY AND ASTROPHYSICS↗

MeV dark matter with MeV dark photons in Abelian kinetic mixing theories

Here we consider the cosmology and phenomenology of a dark photon portal to a simple dark sector consisting of a single, light, fermionic dark matter particle species with mass in the MeV range. We entertain three possible kinetic mixing structures of a new Abelian gauge group U(1) dark with the visible sector through U(1) e.m. , U(1) Y and T[SU(2)_L]. We assume the dark photon to be massive and around the MeV scale, thus close to the mass scale of the dark matter candidate. We compute the dark matter relic density via freeze-out and freeze-in, entertaining the additional possibility of a late inflationary period that could dilute the dark matter yield of heavy candidates, and (ii) additional production modes, for models with under-abundant thermal production. We explore the parameter space compatible with a variety of experimental and astrophysical bounds, and discuss prospects for discovery with new CMB probes and MeV gamma-ray telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating resource-efficient neutron/gamma classification ML models targeting eFPGAs

There has been considerable interest and resulting progress in implementing machine learning (ML) models in hardware over the last several years from the particle and nuclear physics communities. A big driver has been the release of the Python package, hls4ml, which has enabled porting models specified and trained using Python ML libraries to register transfer level (RTL) code. So far, the primary end targets have been commercial field-programmable gate arrays (FPGAs) or synthesized custom blocks on application specific integrated circuits (ASICs). However, recent developments in open-source embedded FPGA (eFPGA) frameworks now provide an alternate, more flexible pathway for implementing ML models in hardware. These customized eFPGA fabrics can be integrated as part of an overall chip design. In general, the decision between a fully custom, eFPGA, or commercial FPGA ML implementation will depend on the details of the end-use application. In this work, we explored the parameter space for eFPGA implementations of fully-connected neural network (fcNN) and boosted decision tree (BDT) models using the task of neutron/gamma classification with a specific focus on resource efficiency. We used data collected using an AmBe sealed source incident on Stilbene, which was optically coupled to an OnSemi J-series silicon photomultiplier (SiPM) to generate training and test data for this study. We investigated relevant input features and the effects of bit-resolution and sampling rate as well as trade-offs in hyperparameters for both ML architectures while tracking total resource usage. The performance metric used to track model performance was the calculated neutron efficiency at a gamma leakage of 10 -3 . The results of the study will be used to aid the specification of an eFPGA fabric, which will be integrated as part of a test chip.

47 OTHER INSTRUMENTATION↗

The landscape of disc outflows from black hole–neutron star mergers

ABSTRACT We investigate mass ejection from accretion discs formed in mergers of black holes (BHs) and neutron stars (NSs). The third observing run of the LIGO/Virgo interferometers provided BH–NS candidate events that yielded no electromagnetic (EM) counterparts. The broad range of disc configurations expected from BH–NS mergers motivates a thorough exploration of parameter space to improve EM signal predictions. Here we conduct 27 high-resolution, axisymmetric, long-term hydrodynamic simulations of the viscous evolution of BH accretion discs that include neutrino emission/absorption effects and post-processing with a nuclear reaction network. In the absence of magnetic fields, these simulations provide a lower limit to the fraction of the initial disc mass ejected. We find a nearly linear inverse dependence of this fraction on disc compactness (BH mass over initial disc radius). The dependence is related to the fraction of the disc mass accreted before the ouflow is launched, which depends on the disc position relative to the innermost stable circular orbit. We also characterize a trend of decreasing ejected fraction and decreasing lanthanide/actinide content with increasing disc mass at fixed BH mass. This trend results from a longer time to reach weak freezout and an increasingly dominant role of neutrino absorption at higher disc masses. We estimate the radioactive luminosity from the disc outflow alone available to power kilonovae over the range of configurations studied, finding a spread of two orders of magnitude. For most of the BH–NS parameter space, the disc outflow contribution is well below the kilonova mass upper limits for GW190814.

79 ASTRONOMY AND ASTROPHYSICS↗