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 55 records · Page 3

Direct Measurement of Dust Attenuation in z approx. 1.5 Star-Forming Galaxies from 3D-HST: Implications for Dust Geometry and Star Formation Rates

The nature of dust in distant galaxies is not well understood, and until recently few direct dust measurements have been possible. We investigate dust in distant star-forming galaxies using near-infrared grism spectra of the 3D-HST survey combined with archival multi-wavelength photometry. These data allow us to make a direct comparison between dust towards star-forming regions (measured using Balmer decrements) and the integrated dust properties (derived by comparing spectral energy distributions [SEDs] with stellar population and dust models) for a statistically significant sample of distant galaxies. We select a sample of 163 galaxies between 1.36< or = z< or = 1.5 with H(alpha) SNR > or = 5 and measure Balmer decrements from stacked spectra. First, we stack spectra in bins of integrated stellar dust attenuation, and find that there is extra dust extinction towards star-forming regions (AV,HII is 1.81 times the integrated AV, star), though slightly lower than found for low-redshift starburst galaxies. Next, we stack spectra in bins of specific star formation rate (log sSFR), star formation rate (log SFR), and stellar mass (logM∗). We find that on average AV,HII increases with SFR and mass, but decreases with increasing sSFR. The amount of extra extinction also decreases with increasing sSFR and decreasing stellar mass. Our results are consistent with the two-phase dust model - in which galaxies contain both a diffuse and a stellar birth cloud dust component - as the extra extinction will increase once older stars outside the star-forming regions become more dominant. Finally, using our Balmer decrements we derive dust-corrected H(alpha) SFRs, and find evidence that SED fitting produces incorrect SFRs if very rapidly declining SFHs are included in the explored parameter space. Subject headings: dust, extinction- galaxies: evolution- galaxies: high-redshift

Direct↗

ULYSSES: Universal LeptogeneSiS Equation Solver

ULYSSES is a python package that calculates the baryon asymmetry produced from leptogenesis in the context of a type-I seesaw mechanism. The code solves the semi-classical Boltzmann equations for points in the model parameter space as specified by the user. We provide a selection of predefined Boltzmann equations as well as a plugin mechanism for externally provided models of leptogenesis. Furthermore, the ULYSSES code provides tools for multi-dimensional parameter space exploration. The emphasis of the code is on user flexibility and rapid evaluation. It is publicly available at https://github.com/earlyuniverse/ulysses

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Exploring DFT+U parameter space with a Bayesian calibration assisted by Markov chain Monte Carlo sampling

Abstract The density-functional theory is widely used to predict the physical properties of materials. However, it usually fails for strongly correlated materials. A popular solution is to use the Hubbard correction to treat strongly correlated electronic states. Unfortunately, the values of the Hubbard U and J parameters are initially unknown, and they can vary from one material to another. In this semi-empirical study, we explore the U and J parameter space of a group of iron-based compounds to simultaneously improve the prediction of physical properties (volume, magnetic moment, and bandgap). We used a Bayesian calibration assisted by Markov chain Monte Carlo sampling for three different exchange-correlation functionals (LDA, PBE, and PBEsol). We found that LDA requires the largest U correction. PBE has the smallest standard deviation and its U and J parameters are the most transferable to other iron-based compounds. Lastly, PBE predicts lattice parameters reasonably well without the Hubbard correction.

36 MATERIALS SCIENCE↗

Numerical simulations of thermal instabilities in stratified gases. II - Exploration of the parameter space

The temporal evolution of density perturbations in an initially hydrostatic isothermal atmosphere consisting of an optically thin radiating compressible plasma is studied. Numerical techniques are used to describe the nonlinear evolution of the perturbations, and the relative equilibrium between dynamic and thermal instabilities as governed by three independent control parameters are examined, namely, the initial density contrast of the perturbation, the ratio of the local buoyancy oscillation period to the local radiative cooling time, and the ratio of the perturbation radius to the local scaleheight. Four orders of magnitude of initial density contrasts and ratios of buoyancy and cooling times, and one order of magnitude of the bubble dimensions are explored. Well-defined oscillations were found to occur in a limited parameter range, and thermal instability to occur even within secondary condensations deriving from the bubble fragmentation.

Reale, F.↗

Exploration of the parameter space of quasisymmetric stellarator vacuum fields through adjoint optimisation

Optimising stellarators for quasisymmetry leads to strongly reduced collisional transport and energetic particle losses compared with unoptimised configurations. Although stellarators with precise quasisymmetry have been obtained in the past, it remains unclear how broad the parameter space is where good quasisymmetry may be achieved. We study the range of aspect ratios and rotational transform values for which stellarators with excellent quasisymmetry on the boundary can be obtained. A large number of Fourier harmonics is included in the boundary representation, which is made computationally tractable by the use of adjoint methods to enable fast gradient-based optimisation and by the direct optimisation of vacuum magnetic fields, which converge more robustly compared with solutions from magnetohydrostatics. Several novel configurations are presented, including stellarators with record levels of quasisymmetry on a surface, three field period quasiaxisymmetric stellarators with substantial magnetic shear, and compact quasisymmetric stellarators at low aspect ratios similar to tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ignition and Burn in a Hybrid Nuclear Fuel for a Pulsed Rocket Engine

This work explores the parameter space of a cylindrical hybrid nuclear fuel in support of a z-pinch driven pulsed fission-fusion (PuFF) engine. 0D power balance and 1D burn wave calculations have been performed to explore the parameter space of hybrid cylindrical nuclear fuels. The boundary of minimal initial conditions needed to reach breakeven are found within the context of the model. The effects of initial conditions upon the yield at the end of burn wave expansion are also determined. The model is used to examine the minimum initial energies predicted to result in yields of a few MJ. The goal of this work is to guide fuel design, inform future models and experiments as well as look for the most efficient parameter space for ignition in a z-pinch driven hybrid nuclear reaction. The impact of initial parameters upon ignition, burn, and gain are discussed. It is found that a hybrid cylindrical target may breakeven with initial energy near the axis of approximately 4 to 6 MJ in lithium deuteride/uranium 235 fuel. It is also found that magnetic field can lower the threshold further of which the magnitude changes across the parameter space. The dual fusion/fission reactions are found to boost each other leading to lower initial driving energies needed to reach breakeven in the hybrid cylindrical nuclear fuel.

fission↗

The Dynamics of Ions on Phased Radio-frequency Carpets in High Pressure Gases and Application for Barium Tagging in Xenon Gas Time Projection Chambers

Radio-frequency (RF) carpets with ultra-fine pitches are examined for ion transport in gases at atmospheric pressures and above. We develop new analytic and computational methods for modeling RF ion transport at densities where dynamics are strongly influenced by buffer gas collisions. An analytic description of levitating and sweeping forces from phased arrays is obtained, then thermodynamic and kinetic principles are used to calculate ion loss rates in the presence of collisions. This methodology is validated against detailed microscopic SIMION simulations. We then explore a parameter space of special interest for neutrinoless double beta decay experiments: transport of barium ions in xenon at pressures from 1 to 10 bar. Our computations account for molecular ion formation and pressure dependent mobility as well as finite temperature effects. We discuss the challenges associated with achieving suitable operating conditions, which lie beyond the capabilities of existing devices, using presently available or near-future manufacturing techniques.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Non-decoupling new particles

We initiate the study of a new class of beyond the Standard Model states that we call “Loryons.” They have the defining characteristic of being non-decoupling, in the sense that their physical mass is dominated by a contribution from the vacuum expectation value of the Higgs boson. The stakes are high: the discovery of a Loryon would tell us that electroweak symmetry must be non-linearly realized in the effective field theory of the Standard Model. Loryons have their masses bounded from above by perturbative unitarity considerations and thus define a finite parameter space for exploration. After providing a complete catalog of Loryon representations under mild assumptions, we turn to examining the constraints on the parameter space from Higgs couplings measurements, precision electroweak tests, and direct collider searches. We show that most fermionic candidates are already ruled out (with some notable exceptions), while much of the scalar Loryon parameter space is still wide open for discovery.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Geo-Statistical Approach to Estimating Asteroid Exploration Parameters

NASA's vision for space exploration calls for a human visit to a near earth asteroid (NEA). Potential human operations at an asteroid include exploring a number of sites and analyzing and collecting multiple surface samples at each site. In this paper two approaches to formulation and scheduling of human exploration activities are compared given uncertain information regarding the asteroid prior to visit. In the first approach a probability model was applied to determine best estimates of mission duration and exploration activities consistent with exploration goals and existing prior data about the expected aggregate terrain information. These estimates were compared to a second approach or baseline plan where activities were constrained to fit within an assumed mission duration. The results compare the number of sites visited, number of samples analyzed per site, and the probability of achieving mission goals related to surface characterization for both cases.

probability↗

Optimization of Plug Nozzles Using CFD Informed Design Space Querying

This paper presents an axisymmetric CFD approach based on Reynolds-averaged Navier-Stokes (RANS) calculations for the shape optimization of plug nozzle contours for both supersonic cruise and landing/takeoff (LTO) conditions. Three key external plug design parameters were defined as independent variables in the optimization scheme, while nozzle performance was evaluated using thrust coefficient (Cfg) and discharge coefficient (Cd) as the primary dependent variables. By varying these parameters, a wide range of axisymmetric conic section plug profiles can be achieved, simplifying the problem. The paper also discusses the supersonic plug nozzle optimization approach, design space seeding, meshing, and CFD simulations for the optimization process. The results demonstrate that the approach is effective in exploring the parameter space and reducing computational expenses, making it a viable option for optimizing nozzle plug contours. Using this optimization toolchain, the study yielded preliminary designs for LTO, cruise, and LTO-cruise balanced plug nozzles.

nozzles↗

Optimization of Plug Nozzles Using CFD Informed Design Space Querying

This paper presents an axisymmetric CFD approach based on Reynolds-averaged Navier-Stokes (RANS) calculations for the shape optimization of plug nozzle contours for both supersonic cruise and landing/takeoff (LTO) conditions. Three key external plug design parameters were defined as independent variables in the optimization scheme, while nozzle performance was evaluated using thrust coefficient (Cfg) and discharge coefficient (Cd) as the primary dependent variables. By varying these parameters, a wide range of axisymmetric conic section plug profiles can be achieved, simplifying the problem. The paper also discusses the supersonic plug nozzle optimization approach, design space seeding, meshing, and CFD simulations for the optimization process. The results demonstrate that the approach is effective in exploring the parameter space and reducing computational expenses, making it a viable option for optimizing nozzle plug contours. Using this optimization toolchain, the study yielded preliminary designs for LTO, cruise, and LTO-cruise balanced plug nozzles.

CFD↗

Guided search for desired functional responses via Bayesian optimization of generative model: Hysteresis loop shape engineering in ferroelectrics

Advances in theoretical modeling across multiple disciplines have yielded generative models capable of high veracity in predicting macroscopic functional responses of materials emerging as a result of complex non-local interactions. Correspondingly, of interest is the inverse problem of finding the model parameter that will yield desired macroscopic responses, such as stress–strain curves, ferroelectric hysteresis loops, etc. Here, we suggest and implement Gaussian process based methods that allow to effectively sample the degenerate parameter space of a complex non-local model to output regions of parameter space which yield desired functionalities. In this work, we discuss the specific adaptation of the acquisition function and sampling function to make the process efficient and balance the efficient exploration of parameter space for multiple possible minima and exploitation to densely sample the regions of interest where target behaviors are optimized. Additionally, this approach is illustrated via the hysteresis loop engineering in ferroelectric materials but can be adapted to other functionalities and generative models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nested active learning for efficient model contextualization and parameterization: pathway to generating simulated populations using multi-scale computational models

There is increasing interest in the use of mechanism-based multi-scale computational models (such as agent-based models (ABMs)) to generate simulated clinical populations in order to discover and evaluate potential diagnostic and therapeutic modalities. The description of the environment in which a biomedical simulation operates (model context) and parameterization of internal model rules (model content) requires the optimization of a large number of free parameters. In this work, we utilize a nested active learning (AL) workflow to efficiently parameterize and contextualize an ABM of systemic inflammation used to examine sepsis. Contextual parameter space was examined using four parameters external to the model’s rule set. The model’s internal parameterization, which represents gene expression and associated cellular behaviors, was explored through the augmentation or inhibition of signaling pathways for 12 signaling mediators associated with inflammation and wound healing. We have implemented a nested AL approach in which the clinically relevant (CR) model environment space for a given internal model parameterization is mapped using a small Artificial Neural Network (ANN). The outer AL level workflow is a larger ANN that uses AL to efficiently regress the volume and centroid location of the CR space given by a single internal parameterization. We have reduced the number of simulations required to efficiently map the CR parameter space of this model by approximately 99%. In addition, we have shown that more complex models with a larger number of variables may expect further improvements in efficiency.

97 MATHEMATICS AND COMPUTING↗

Anatomy of the electroweak phase transition for dark sector induced baryogenesis

We investigate the electroweak phase transition patterns for a recently proposed baryogenesis model with CP violation originated in the dark sector. The model includes a complex scalar singlet-Higgs boson portal, a U(1) l gauge lepton symmetry with a Z ' gauge boson portal and a fermionic dark matter particle. We find a novel thermal history of the scalar sector, featuring a Z 2 breaking singlet vacuum in the early Universe driven by a dark Yukawa coupling, that induces a one-step strongly first order electroweak phase transition. We explore the parameter space that generates the observed matter-antimatter asymmetry and dark matter relic abundance, while being consistent with constraints from electric dipole moment, collider searches, and dark matter direct detection bounds. The complex singlet can be produced via the Higgs portal and decays into Standard Model particles after traveling a certain distance. We explore the reach for long-lived singlet scalars at the 13 TeV Large Hadron Collider with \( \mathcal{L} \) = 139 fb - 1 and show its impact on the parameter space of the model. Setting aside currently unresolved theoretical uncertainties, we estimate the gravitational wave signatures detectable at future observatories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Hypothesis Learning in Automated Experiment: Application to Combinatorial Materials Libraries

Machine learning is rapidly becoming an integral part of experimental physical discovery via automated and high-throughput synthesis, and active experiments in scattering and electron/probe microscopy. This, in turn, necessitates the development of active learning methods capable of exploring relevant parameter spaces with the smallest number of steps. In this work, an active learning approach based on conavigation of the hypothesis and experimental spaces is introduced. This is realized by combining the structured Gaussian processes containing probabilistic models of the possible system's behaviors (hypotheses) with reinforcement learning policy refinement (discovery). This approach closely resembles classical human-driven physical discovery, when several alternative hypotheses realized via models with adjustable parameters are tested during an experiment. This approach is demonstrated for exploring concentration-induced phase transitions in combinatorial libraries of Sm-doped BiFeO 3 using piezoresponse force microscopy, but it is straightforward to extend it to higher-dimensional parameter spaces and more complex physical problems once the experimental workflow and hypothesis generation are available.

36 MATERIALS SCIENCE↗

Recursive Branching Simulated Annealing Algorithm

This innovation is a variation of a simulated-annealing optimization algorithm that uses a recursive-branching structure to parallelize the search of a parameter space for the globally optimal solution to an objective. The algorithm has been demonstrated to be more effective at searching a parameter space than traditional simulated-annealing methods for a particular problem of interest, and it can readily be applied to a wide variety of optimization problems, including those with a parameter space having both discrete-value parameters (combinatorial) and continuous-variable parameters. It can take the place of a conventional simulated- annealing, Monte-Carlo, or random- walk algorithm. In a conventional simulated-annealing (SA) algorithm, a starting configuration is randomly selected within the parameter space. The algorithm randomly selects another configuration from the parameter space and evaluates the objective function for that configuration. If the objective function value is better than the previous value, the new configuration is adopted as the new point of interest in the parameter space. If the objective function value is worse than the previous value, the new configuration may be adopted, with a probability determined by a temperature parameter, used in analogy to annealing in metals. As the optimization continues, the region of the parameter space from which new configurations can be selected shrinks, and in conjunction with lowering the annealing temperature (and thus lowering the probability for adopting configurations in parameter space with worse objective functions), the algorithm can converge on the globally optimal configuration. The Recursive Branching Simulated Annealing (RBSA) algorithm shares some features with the SA algorithm, notably including the basic principles that a starting configuration is randomly selected from within the parameter space, the algorithm tests other configurations with the goal of finding the globally optimal solution, and the region from which new configurations can be selected shrinks as the search continues. The key difference between these algorithms is that in the SA algorithm, a single path, or trajectory, is taken in parameter space, from the starting point to the globally optimal solution, while in the RBSA algorithm, many trajectories are taken; by exploring multiple regions of the parameter space simultaneously, the algorithm has been shown to converge on the globally optimal solution about an order of magnitude faster than when using conventional algorithms. Novel features of the RBSA algorithm include: 1. More efficient searching of the parameter space due to the branching structure, in which multiple random configurations are generated and multiple promising regions of the parameter space are explored; 2. The implementation of a trust region for each parameter in the parameter space, which provides a natural way of enforcing upper- and lower-bound constraints on the parameters; and 3. The optional use of a constrained gradient- search optimization, performed on the continuous variables around each branch s configuration in parameter space to improve search efficiency by allowing for fast fine-tuning of the continuous variables within the trust region at that configuration point.

Bolcar, Matthew↗

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↗