Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Bounds”

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 109 records · Page 6

A Survey on Error-Bounded Lossy Compression for Scientific Datasets

Error-bounded lossy compression has been effective in significantly reducing the data storage/transfer burden while preserving the reconstructed data fidelity very well. Many error-bounded lossy compressors have been developed for a wide range of parallel and distributed use cases for years. They are designed with distinct compression models and principles, such that each of them features particular pros and cons. In this article, we provide a comprehensive survey of emerging error-bounded lossy compression techniques. The key contribution is fourfold. (1) We summarize a novel taxonomy of lossy compression into six classic models. (2) We provide a comprehensive survey of 10 commonly used compression components/modules. (3) We summarized pros and cons of 47 state-of-the-art lossy compressors and present how state-of-the-art compressors are designed based on different compression techniques. (4) We discuss how customized compressors are designed for specific scientific applications and use-cases. We believe this survey is useful to multiple communities including scientific applications, high-performance computing, lossy compression, and big data.

Error-Bounded Lossy Compression↗

Empirical Bounding Analysis and User Recommendations for a Neutron Multiplicity Detector

Neutron multiplicity detectors are useful for a variety of applications including nuclear emergency response, nuclear nonproliferation, safeguards, and criticality safety. When measuring black-box problems (i.e., when the system being measured is completely unknown) with systems that have relatively low detection efficiency, expert analysis is frequently required to determine appropriate bounds for system parameters such as neutron multiplication and neutron leakage. This is because the detection efficiency can vary wildly, which also means that the resulting system parameters can have large variations. This work applies an empirical approach to provide bounds on system parameters (such as system multiplication, neutron leakage, and detector efficiency). The results of this bounding analysis are then used as inputs to provide recommendations to users regarding criticality safety. In addition to describing the method used, this work provides sample results from a measurement campaign performed at the National Criticality Experiments Research Center (NCERC).

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗

A Scalable Interior‐Point Gauss–Newton Method for PDE‐Constrained Optimization With Bound Constraints

Here, we present a scalable approach to solve a class of partial differential equation (PDE)‐constrained optimization problems with bound constraints. This approach utilizes a robust full‐space interior‐point (IP)‐Gauss–Newton optimization method. To cope with the poorly‐conditioned IP‐Gauss–Newton saddle‐point linear systems that need to be solved approximately, once per optimization step, we propose two spectrally related preconditioners. These preconditioners leverage the limited informativeness of data in regularized PDE‐constrained optimization problems. A block Gauss–Seidel preconditioner is proposed for the GMRES‐based solution of the IP‐Gauss–Newton linear systems. It is shown, for a large‐class of PDE‐ and bound‐constrained optimization problems, that the spectrum of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix is asymptotically independent of discretization and is not impacted by the ill‐conditioning that notoriously plagues interior‐point methods. We exploit symmetry of the IP‐Gauss–Newton linear systems and propose a regularization and log‐barrier Hessian preconditioner for the preconditioned conjugate gradient (PCG)‐based solution of the equivalent IP‐Gauss–Newton–Schur complement linear systems. The eigenvalues of the block Gauss–Seidel preconditioned IP‐Gauss–Newton matrix, that are not equal to one, are identical to the eigenvalues of the regularization and log‐barrier Hessian preconditioned Schur complement matrix. The scalability of the approach is demonstrated on two example problems. The numerical solution of these optimization problems is shown to require a discretization independent number of IP‐Gauss–Newton linear solves. Furthermore, the linear systems are solved in a discretization and IP ill‐conditioning independent number of preconditioned Krylov subspace iterations. The parallel scalability of the preconditioner, achieved via algebraic multigrid component solvers when applicable, and the aforementioned algorithmic scalability permits a parallel scalable means to compute solutions of a large class of PDE‐ and bound‐constrained problems.

PDE-constrained optimization↗

Analysis of sparse recovery for Legendre expansions using envelope bound

We provide novel sufficient conditions for the uniform recovery of sparse Legendre expansions using ℓ 1 minimization, where the sampling points are drawn according to orthogonalization (uniform) measure. So far, conditions of the form m ≳ Θ 2 s x log factors have been relied on to determine the minimum number of samples m that guarantees successful reconstruction of s-sparse vectors when the measurement matrix is associated to an orthonormal system. However, in case of sparse Legendre expansions, the uniform bound Θ of Legendre systems is so high that these conditions are unable to provide meaningful guarantees. Here, in this paper, we present an analysis which employs the envelop bound of all Legendre polynomials instead, and prove a new recovery guarantee for s-sparse Legendre expansions, m ≳ Θs 2 x log factors, which is independent of Θ. Arguably, this is the first recovery condition established for orthonormal systems without assuming the uniform boundedness of the sampling matrix. The key ingredient of our analysis is an extension of chaining arguments, recently developed in Bourgain and Chkifa et al., to handle the envelope bound. Furthermore, our recovery condition is proved via restricted eigenvalue property, a less demanding replacement of restricted isometry property which is perfectly suited to the considered scenario. Along the way, we derive simple criteria to detect good sample sets. Our numerical tests show that sets of uniformly sampled points that meet these criteria will perform better recovery on average.

97 MATHEMATICS AND COMPUTING↗

Five-point superluminality bounds

We investigate how the speed of propagation of physical excitations is encoded in the coefficients of five-point interactions. This leads to a superluminality bound on scalar five-point interactions, which we present here for the first time. To substantiate our result, we also consider the case of four-point interactions for which bounds from S-matrix sum rules exist and show that these are parametrically equivalent to the bounds obtained within our analysis. Finally, we extend the discussion to a class of higher-point interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improving modular bootstrap bounds with integrality

We propose methods that efficiently impose integrality — i.e., the condition that the coefficients of characters in the partition function must be integers — into numerical modular bootstrap. We demonstrate the method with a number of examples where it can be used to strengthen modular bootstrap results. First, we show that, with a mild extra assumption, imposing integrality improves the bound on the maximal allowed gap in dimensions of operators in theories with a U(1) c symmetry at c = 3, and reduces it to the value saturated by the SU(4) 1 WZW model point of c = 3 Narain lattices moduli space. Second, we show that our method can be used to eliminate all but a discrete set of points saturating the bound from previous Virasoro modular bootstrap results. Finally, when central charge is close to 1, we can slightly improve the upper bound on the scaling dimension gap.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analytic bounds on late-time axion-scalar cosmologies

The cosmological dynamics of multiple scalar/pseudoscalar fields are difficult to solve, especially when the field-space metric is curved. This presents a challenge in determining whether a given model can support cosmic acceleration, without solving for the on-shell solution. In this work, we present bounds on late-time FLRW-cosmologies in classes of theories that involve arbitrary numbers of scalar and pseudoscalar fields coupled both kinetically (leading to a curved field space metric) and through scalar potentials. Such bounds are proven analytically, independently of initial conditions, with no approximation in the field equations and without referring to explicit solutions. Besides their broad applications to cosmological model building, our bounds can be applied to studying asymptotic cosmologies of certain classes of string compactifications.

79 ASTRONOMY AND ASTROPHYSICS↗

AdS bulk locality from sharp CFT bounds

It is a long-standing conjecture that any CFT with a large central charge and a large gap Δ gap in the spectrum of higher-spin single-trace operators must be dual to a local effective field theory in AdS. We prove a sharp form of this conjecture by deriving numerical bounds on bulk Wilson coefficients in terms of Δ gap using the conformal bootstrap. Our bounds exhibit the scaling in Δ gap expected from dimensional analysis in the bulk. Our main tools are dispersive sum rules that provide a dictionary between CFT dispersion relations and S-matrix dispersion relations in appropriate limits. This dictionary allows us to apply recently-developed flat-space methods to construct positive CFT functionals. We show how AdS 4 naturally resolves the infrared divergences present in 4D flat-space bounds. Our results imply the validity of twice-subtracted dispersion relations for any S-matrix arising from the flat-space limit of AdS/CFT.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Universal bounds on CFT Distance Conjecture

For any unitary conformal field theory in two dimensions with the central charge c, we prove that, if there is a nontrivial primary operator whose conformal dimension ∆ vanishes in some limit on the conformal manifold, the Zamolodchikov distance t to the limit is infinite, the approach to this limit is exponential ∆ = exp(−αt + O(1)), and the decay rate obeys the universal bounds c−1/2 ≤ α ≤ 1. In the limit, we also find that an infinite tower of primary operators emerges without a gap above the vacuum and that the conformal field theory becomes locally a tensor product of a sigma-model in the large radius limit and a compact theory. As a corollary, we establish a part of the Distance Conjecture about gravitational theories in three-dimensional anti-de Sitter space. In particular, our bounds on α indicate that the emergence of exponentially light states is inevitable as the moduli field corresponding to t rolls beyond the Planck scale along the steepest path and that this phenomenon can begin already at the curvature scale of the bulk geometry. We also comment on implications of our bounds for gravity in asymptotically flat spacetime by taking the flat space limit and compare with the Sharpened Distance Conjecture.

AdS-CFT Correspondence↗

A Converse to Lieb–Robinson Bounds in One Dimension Using Index Theory

Abstract Unitary dynamics with a strict causal cone (or “light cone”) have been studied extensively, under the name of quantum cellular automata (QCAs). In particular, QCAs in one dimension have been completely classified by an index theory. Physical systems often exhibit only approximate causal cones; Hamiltonian evolutions on the lattice satisfy Lieb–Robinson bounds rather than strict locality. This motivates us to study approximately locality preserving unitaries (ALPUs). We show that the index theory is robust and completely extends to one-dimensional ALPUs. As a consequence, we achieve a converse to the Lieb–Robinson bounds: any ALPU of index zero can be exactly generated by some time-dependent, quasi-local Hamiltonian in constant time. For the special case of finite chains with open boundaries, any unitary satisfying the Lieb–Robinson bound may be generated by such a Hamiltonian. We also discuss some results on the stability of operator algebras which may be of independent interest.

Physics↗

Tighter reformulations using classical Dawson and Sankoff bounds for approximating two-stage chance-constrained programs

We extend and improve recent results given by Singh and Watson on using classical bounds on the union of sets in a chance-constrained optimization problem. Specifically, we revisit the so-called Dawson and Sankoff bound that provided one of the best approximations of a chance constraint in the previous analysis. Next, we show that our work is a generalization of the previous work, and in fact the inequality employed previously is a very relaxed approximation with assumptions that do not generally hold. Computational results demonstrate on average over a 43% improvement in the bounds. As a byproduct, we provide an exact reformulation of the floor function in optimization models.

97 MATHEMATICS AND COMPUTING↗

Modeling approaches for addressing unrelaxable bound constraints with unconstrained optimization methods

We explore novel approaches for solving nonlinear optimization problems with unrelaxable bound constraints, which must be satisfied before the objective function can be evaluated. Our method reformulates the unrelaxable bound-constrained problem as an unconstrained optimization problem that is amenable to existing unconstrained optimization methods. The reformulation relies on a domain warping to form a merit function; the choice of the warping determines the level of exactness with which the unconstrained problem can be used to find solutions to the bound-constrained problem, as well as key properties of the unconstrained formulation such as smoothness. We develop theory when the domain warping is a multioutput sigmoidal warping, and we explore the practical elements of applying unconstrained optimization methods to the formulation. Here, we develop an algorithm that exploits the structure of the sigmoidal warping to guarantee that unconstrained optimization algorithms applied to the merit function will find a stationary point to the desired tolerance.

97 MATHEMATICS AND COMPUTING↗

Sequential membrane- and protein-bound organelles compartmentalize genomes during phage infection

Many eukaryotic viruses require membrane-bound compartments for replication, but no such organelles are known to be formed by prokaryotic viruses. Bacteriophages of the Chimalliviridae family sequester their genomes within a phage-generated organelle, the phage nucleus, which is enclosed by a lattice of the viral protein ChmA. We show that inhibiting phage nucleus formation arrests infections at an early stage in which the injected phage genome is enclosed within a membrane-bound early phage infection (EPI) vesicle. Early phage genes are expressed from the EPI vesicle, demonstrating its functionality as a prokaryotic, transcriptionally active, membrane-bound organelle. We also show that the phage nucleus is essential, with genome replication beginning after the injected DNA is transferred from the EPI vesicle to the phage nucleus. Our results show that Chimalliviridae require two sophisticated subcellular compartments of distinct compositions and functions that facilitate successive stages of the viral life cycle.

59 BASIC BIOLOGICAL SCIENCES↗

Theoretical estimates of flammability bounds for thin condensed fuel diffusion flames in microgravity using detailed models of chemistry and radiation

Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Isotope Exchange between Mercuric [Hg(II)] Chloride and Hg(II) Bound to Minerals and Thiolate Ligands: Implications for Enriched Isotope Tracer Studies

Enriched mercury (Hg) stable isotopes have been widely used as tracers in field and laboratory investigations of Hg(II) biogeochemical transformations such as methylation and demethylation. Few studies, however, have considered concurrent isotope exchange reactions between newly spiked and pre-existing Hg(II) in environmental matrices, which may alter redistribution and thus transformation of the spiked and pre-existing Hg(II). Using enriched 198 Hg [as mercuric Hg(II) or HgCl n species], this study investigated isotope exchange between 198 Hg and pre-existing Hg(II) bound to metacinnabar (β-HgS), sediments, low-molecular-weight (LMW) thiols, and dissolved organic matter (DOM). The impact of isotope exchange on methylmercury production in the presence of organic ligands was also evaluated with an iron-reducing bacterium Geobacter sulfurreducens PCA in a phosphate buffered solution (pH 7.4). Here, we found that spiked 198 Hg readily exchanged with mineral-bound ambient Hg(II) despite concurrent Hg(II) adsorption and immobilization on the solids. Rapid exchange (<3 min) was also observed between spiked 198 Hg and 200 Hg pre-equilibrated with LMW thiols and DOM in solution. While the exchange did not cause net changes in Hg(II) chemical speciation, it resulted in redistribution of Hg(II) isotopes bound to the ligands and thus an apparently similar methylation rate and magnitude of the spiked 198 Hg and pre-existing 200 Hg by PCA cells when 198 Hg and 200 Hg were added at 1:1 ratio. These observations underscore the importance of isotope exchange when an enriched Hg isotope is applied in environmental matrices, as the exchange could potentially lead to biased rate calculations of Hg(II) transformation and bioaccumulation and thus risk assessments of new Hg(II) input to the natural ecosystems.

58 GEOSCIENCES↗

Search for a bound di-neutron by comparing 3 He(e,e'p)d and 3 H(e,e'p)X measurements

We report on a search for a bound di-neutron by comparing electron-induced proton-knockout $(e,e'p)$ measurements from Helium-3 ($^3$He) and Tritium ($^3$H). The measurements were performed at Jefferson Lab Hall A with a 4.326 GeV electron beam, and kinematics of large momentum transfer $Q^2 \approx 1.9$ (GeV/$c$)$^2$ and $x_B>1$, to minimize contributions from non quasi-elastic (QE) reaction mechanisms. Analyzing the measured $^3$He missing mass ($M_{miss}$) and missing energy ($E_{miss}$) distributions, we can distinguish the two-body break-up reaction, in which the residual proton-neutron system remains bound as a deuteron. In the $^3$H mirror case, under the exact same kinematic conditions, we do not identify a signature for a bound di-neutron with similar binding energy to that of the deuteron. We calculate exclusion limits as a function of the di-neutron binding energy and find that, for binding equivalent to the deuteron, the two-body break-up cross section on $^3$H is less than 0.9% of that on $^3$He in the measured kinematics at the 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Particle removal from a flat surface using a translating bounded vortex flow

A bounded vortex flow is a hydrodynamic approach for removal of particles from a surface without scattering the particles onto nearby surfaces. The bounded vortex flow field is generated by a nozzle that combines azimuthally tilted jets arranged in a circular pattern and a central suction port. When the nozzle face is directed toward an ‘impingement surface’, the flow develops a wall-normal intake vortex below the suction outlet, which causes high shear stress on the impingement surface. When particles are present on the impingement surface, the high shear stress causes particles to roll along the surface and to be lifted off the surface and transported up the core of the wall-normal vortex into the suction outlet. In typical applications, the nozzle would be translated along the impingement surface to clean particles from the surface. The current paper reports on an experimental study of the effect of nozzle translation on the effectiveness of the bounded vortex flow field for particle mitigation. The effectiveness of particle mitigation was examined as a function of flow rate through the nozzle, particle size, and nozzle translation velocity relative to the impingement surface. As a result, numerical computations are used to relate the flow rate to the maximum shear stress on the impingement surface, which is then used to theoretically predict onset of particle motion.

42 ENGINEERING↗