Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constraining”

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 847 records · Page 47

Clarity through the neutrino fog: constraining new forces in dark matter detectors

The PandaX-4T and XENONnT experiments present indications of Coherent Elastic Neutrino Nucleus Scattering (CEνNS) from 8 B solar neutrinos at 2.6σ and 2.7σ, respectively. This constitutes the first observation of the neutrino “floor” or “fog”, an irreducible background that future dark matter searches in terrestrial detectors will have to contend with. Here, we first discuss the contributions from neutrino–electron scattering and from the Migdal effect in the region of interest of these experiments, and we argue that they are non-negligible. Second, we make use of the recent PandaX-4T and XENONnT data to derive novel constraints on light scalar and vector mediators coupling to neutrinos and quarks. We demonstrate that these experiments already provide world-leading laboratory constraints on new light mediators in some regions of parameter space.

neutrino interactions↗

Constraining the Chiral Magnetic Effect with charge-dependent azimuthal correlations in Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 2.76 and 5.02 TeV

Systematic studies of charge-dependent two- and three-particle correlations in Pb-Pb collisions at $ \sqrt{s_{\mathrm{NN}}} $ = 2.76 and 5.02 TeV used to probe the Chiral Magnetic Effect (CME) are presented. These measurements are performed for charged particles in the pseudorapidity ( η ) and transverse momentum ( p T ) ranges | η | < 0 . 8 and 0 . 2 < p T < 5 GeV/ c . A significant charge-dependent signal that becomes more pronounced for peripheral collisions is reported for the CME-sensitive correlators γ 1, 1 =$\langle4cos( φ α + φ β - 2Ψ 2 )$\rangle$ and γ 1, - 3 = $\langle$cos( φ α - 3 φ β + 2Ψ 2 )$\rangle$. The results are used to estimate the contribution of background effects, associated with local charge conservation coupled to anisotropic flow modulations, to measurements of the CME. A blast-wave parametrisation that incorporates local charge conservation tuned to reproduce the centrality dependent background effects is not able to fully describe the measured γ 1 , 1 . Finally, the charge and centrality dependence of mixed-harmonics three-particle correlations, of the form γ 1, 2 = $\langle$cos( φ α + 2 φ β - 3Ψ 3 )$\rangle$, which are insensitive to the CME signal, verify again that background contributions dominate the measurement of γ 1 , 1 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Adaptive mesh refinement in stress-constrained topology optimization

We present a topology structural optimization framework with adaptive mesh refinement and stress-constraints. Finite element approximation and geometry representation benefit from such refinement by enabling more accurate stress field predictions and greater resolution of the optimal structural boundaries. We combine a volume fraction filter to impose a minimum design feature size, the RAMP penalization to generate “black-and-white designs” and a RAMP-like stress definition to resolve the “stress singularity problem.” Regions with stress concentrations dominate the optimized design. As such, rigorous simulations are required to accurately approximate the stress field. To achieve this goal, we invoke a threshold operation and mesh refinement during the optimization. We do so in an optimal fashion, by applying adaptive mesh refinement techniques that use error indicators to refine and coarsen the mesh as needed. In this way, we obtain more accurate simulations and greater resolution of the design domain. We present results in two dimensions to demonstrate the efficiency of our method.

97 MATHEMATICS AND COMPUTING↗

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗

A geochemical approach to constraining the formation of glassy fallout debris from nuclear tests

Reprocessed earth material is a glassy nuclear fallout debris from near-surface nuclear tests. A geochemical approach to analysis of glassy fallout is uniquely suited to determine the means of reprocessing and shed light on the mechanisms of fallout formation. An improved understanding of fallout formation is of interest both for its potential to guide post-detonation nuclear forensic investigations and in the context of possible affinities between glassy debris and other glasses generated by high-energy natural events, such as meteorite impacts and lightning strikes. Our study presents a large major-element compositional dataset for glasses within aerodynamic fallout from the Trinity nuclear test (“trinitite”) and a geochemically based analysis of the glass compositional trends. Silica-rich and alkali-rich trinitite glasses show compositions and textures consistent with formation through melting of individual mineral grains—quartz and alkali feldspar, respectively—from the test-site sediment. Furthermore, the volumetrically dominant glass phase—called the CaMgFe glass—shows extreme major-element compositional variability. Compositional trends in the CaMgFe glass are most consistent with formation through volatility-controlled condensation from compositionally heterogeneous plasma. Radioactivity occurs only in CaMgFe glass, indicating that co-condensation of evaporated bulk ground material and trace device material was the main mechanism of radioisotope incorporation into trinitite. CaMgFe trinitite glasses overlap compositionally with basalts, rhyolites, fulgurites, tektites, and microtektites but display greater compositional diversity than all of these naturally formed glasses. Indeed, the most refractory CaMgFe glasses compositionally resemble early solar system condensates—specifically, CAIs.

58 GEOSCIENCES↗

Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes

Abstract Mafic volcanic activity is dominated by effusive to mildly explosive eruptions. Plinian and ignimbrite-forming mafic eruptions, while rare, are also possible; however, the conditions that promote such explosivity are still being explored. Eruption style is determined by the ability of gas to escape as magma ascends, which tends to be easier in low-viscosity, mafic magmas. If magma permeability is sufficiently high to reduce bubble overpressure during ascent, volatiles may escape from the magma, inhibiting violent explosive activity. In contrast, if the permeability is sufficiently low to retain the gas phase within the magma during ascent, bubble overpressure may drive magma fragmentation. Rapid ascent may induce disequilibrium crystallization, increasing viscosity and affecting the bubble network with consequences for permeability, and hence, explosivity. To explore the conditions that promote strongly explosive mafic volcanism, we combine microlite textural analyses with synchrotron x-ray computed microtomography of 10 pyroclasts from the 12.6 ka mafic Curacautín Ignimbrite (Llaima Volcano, Chile). We quantify microlite crystal size distributions (CSD), microlite number densities, porosity, bubble interconnectivity, bubble number density, and geometrical properties of the porous media to investigate the role of magma degassing processes at mafic explosive eruptions. We use an analytical technique to estimate permeability and tortuosity by combing the Kozeny-Carman relationship, tortuosity factor, and pyroclast vesicle textures. The groundmass of our samples is composed of up to 44% plagioclase microlites, > 85% of which are < 10 µm in length. In addition, we identify two populations of vesicles in our samples: (1) a convoluted interconnected vesicle network produced by extensive coalescence of smaller vesicles (> 99% of pore volume), and (2) a population of very small and completely isolated vesicles (< 1% of porosity). Computed permeability ranges from 3.0 × 10 −13 to 6.3 × 10 −12 m 2 , which are lower than the similarly explosive mafic eruptions of Tarawera (1886; New Zealand) and Etna (112 BC; Italy). The combination of our CSDs, microlite number densities, and 3D vesicle textures evidence rapid ascent that induced high disequilibrium conditions, promoting rapid syn-eruptive crystallization of microlites within the shallow conduit. We interpret that microlite crystallization increased viscosity while simultaneously forcing bubbles to deform as they grew together, resulting in the permeable by highly tortuous network of vesicles. Using the bubble number densities for the isolated vesicles (0.1-3 −3 × 10 4 bubbles per mm 3 ), we obtain a minimum average decompression rate of 1.4 MPa/s. Despite the textural evidence that the Curacautín magma reached the percolation threshold, we propose that rapid ascent suppressed outgassing and increased bubble overpressures, leading to explosive fragmentation. Further, using the porosity and permeability of our samples, we estimated that a bubble overpressure > 5 MPa could have been sufficient to fragment the Curacautín magma. Other mafic explosive eruptions report similar disequilibrium conditions induced by rapid ascent rate, implying that syn-eruptive disequilibrium conditions may control the explosivity of mafic eruptions more generally.

Valdivia, Pedro↗

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗

Subexponentially Growing Hilbert Space and Nonconcentrating Distributions in a Constrained Spin Model

Motivated by recent experiments with two-component Bose–Einstein condensates, we study fully-connected spin models subject to an additional constraint. The constraint is responsible for the Hilbert space dimension to scale only linearly with the system size. We discuss the unconventional statistical physical and thermodynamic properties of such a system, in particular the absence of concentration of the underlying probability distributions. As a consequence, expectation values are less suitable to characterize such systems, and full distribution functions are required instead. Sharp signatures of phase transitions do not occur in such a setting, but transitions from singly peaked to doubly peaked distribution functions of an “order parameter” may be present.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Consequences and benefits of utilizing continuous vibro-impact representations in constrained pipeline conveying fluid systems

Here, the effectiveness of continuous vibro-impact forcing representations for the cantilevered pipe that conveys fluid is explored and analyzed. The previously accepted forcing model utilizing a smoothened trilinear spring is estimated using three continuous forcing representations, namely, polynomial, rational polynomial, and hyperbolic tangent. The accuracy of the estimated forcing functions is investigated and analyzed by calculating the root mean square error, and bifurcation diagrams are generated and compared to the nominal system. Additionally, the dynamic response of the system is further characterized using Poincare maps, power spectra, and basins of attraction. Once all continuous forcing representations are analyzed and compared to the nominal system, the computational cost of each method is examined, and further limitations of the hyperbolic tangent method are discovered. It is proved that the hyperbolic tangent forcing representation most accurately captures the dynamic response of the pipeline, and the least accurate representation is the rational polynomial representation. Additionally, considerable computational cost is saved when employing the hyperbolic tangent representation compared to the discontinuous representation.

42 ENGINEERING↗

Constraining self-interacting dark matter with the full dataset of PandaX-II

Self-interacting dark matter (SIDM) is a leading candidate proposed to solve discrepancies between predictions of the prevailing cold dark matter theory and observations of galaxies. Many SIDM models predict the existence of a light force carrier that mediates strong dark matter self-interactions. If the mediator couples to the standard model particles, it could produce characteristic signals in dark matter direct detection experiments. We report searches for signals of SIDM models with a light mediator using the full dataset of the PandaX-II experiment, basing on a total exposure of 132 tonne-days. No significant excess over background is found, and our likelihood analysis leads to a strong upper limit on the dark matter-nucleon coupling strength. We further combine the PandaX-II constraints and those from observations of the light element abundances in the early universe, and show that direct detection and cosmological probes can provide complementary constraints on dark matter models with a light mediator.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementing a unified solver for nonlinearly constrained optimization

SQP and interior-point methods (also referred to as Lagrange-Newton methods) typically share key algorithmic components, such as strategies for computing descent directions and mechanisms that promote global convergence. Building on this insight, we introduce a unifying framework with eight building blocks that abstracts the workflows of Lagrange-Newton methods. We then present Uno, a modular C++ solver that implements our unifying framework and allows the automatic combination of a wide range of strategies with no programming effort from the user. Uno is meant to (1) organize mathematical optimization strategies into a coherent hierarchy; (2) offer a wide range of efficient and robust methods that can be compared for a given instance; (3) enable researchers to experiment with novel optimization strategies; and (4) reduce the cost of development and maintenance of multiple optimization solvers. Uno’s software design allows user to compose new customized solvers for emerging optimization areas such as robust optimization or optimization problems with complementarity constraints, while building on reliable nonlinear optimization techniques. We demonstrate that Uno is highly competitive against state-of-the-art solvers filterSQP, IPOPT, SNOPT, MINOS, LANCELOT, LOQO, and CONOPT on a subset of 429 small problems from the CUTE collection. Uno is available as open-source software under the MIT license at https://github.com/cvanaret/Uno and via its C, Julia, Python, Fortran, and AMPL interfaces.

97 MATHEMATICS AND COMPUTING↗

Sommerfeld effect in a constrained electromechanical system

In this paper, the dynamics of an electromechanical system is investigated. In this model, the phase of the system is a dynamical variable. The existence of four periodic orbits is proved. Two of them are asymptotically stable and the others are unstable. In the absence of viscous damping, when the parameters associated with the dimensionless voltages are adequately changed, each orbit is subject to a change of stability of the kind stable → unstable → stable with repetition of this pattern. This phenomenon is known as Sommerfeld effect. Moreover, the stability of these orbits depends on the zeros of the Bessel function J{sub 1}.

97 MATHEMATICS AND COMPUTING↗

The use of Jacobi wavelets for constrained approximation of rational Bézier curves

This paper presents an efficient method to solve the approximation problem of the rational Bézier curve by continuous piecewise polynomial curve in L{sub 2}-norm. For this purpose, extended Jacobi wavelets together with the Gauss–Jacobi quadrature rules are employed. The proposed technique has some advantages such as: simplicity, high accuracy and fast computation. Also, the method in the paper performs multi-degree reduction at one time and does not need the stepwise computing. Several examples are given to illustrate the effectiveness of the algorithm.

97 MATHEMATICS AND COMPUTING↗