SEARCH · Engineering Papers
Results for “VIOLATIONS”
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.
Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays
The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.
Lattice dynamics of high-pressure hydrides studied by inelastic neutron scattering
Due to the small mass and anomalously large neutron scattering cross-section of proton (about 80 barns compared to a few barns for other nuclei), inelastic neutron scattering is considered as one of the most effective tools in studying optical vibrations of hydrogen atoms in metal hydrides. The current review is focused on the binary hydrides of 3d- and 4d-metals of groups VI–VIII, which were produced at high hydrogen pressures of several gigapascals in relatively large quantities of hundreds of mg, quenched to low temperature and studied by INS ex situ at ambient pressure with high statistical accuracy. One of the unusual effects revealed by INS is a strong increase in the strength of the metal-hydrogen interactions with decreasing atomic number of the d-metal accompanied by an increase in the Me-H distance. Based on the available experimental results, the spectra g(E) of the phonon density of states and temperature dependencies CV(T) of the heat capacity at constant volume at T up to 1000 K have been derived in this paper and presented both in the figures and in digital form. This provides the reference data for the theoretical investigations of the crystal structures and compositions of new practically important hydrides giving the opportunity to validate calculation methods by comparing the calculated g(E) and CV(T) with the accurate experimental dependencies for the binary hydrides. Recent INS studies showed [R.A. Klein et al., J. Alloy. Compd. 894 (2022) 162381] that the fingerprints of anomalously short H-H separations of 1.6 Å violating the “2 Å rule” can be easily and unambiguously identified in the complex INS spectra of quaternary hydrides (La,Ce)NiInH 1+x . This makes neutron spectroscopy an attractive means for obtaining valuable data in the search for novel hydrides with a record high hydrogen capacity.
Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations
A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.
A novel energy balance approach for a verifiable and accurate solution of radiation extinction in purely absorbing particle clouds
Here we consider the problem of radiation transport through purely absorbing particle clouds. The gold-standard solution of particle-resolved Monte Carlo ray-tracing method to this problem is computationally expensive and therefore solving the radiation transport equation (specifically, the Beer-Bouguer law (BB-law)) on a Eulerian mesh is often preferred. While the absorption coefficient in the real problem is infinite, the BB-law approximates it to be a finite number in the form of number density through a set of assumptions. For particle clouds that do not obey these assumptions, the BB-law predicts an incorrect exponential decay. Also, when the number density is computed using the nearest-neighbor approach, the BB-law solution diverges when the Eulerian mesh size becomes closer to or smaller than the particle size. This numerical divergence is due to the homogenization error. Although the filtering strategy for number density minimizes the homogenization error, it still converges to an incorrect exponential decay for particle clouds that break the BB-law constraints and the cost of the filtering is equivalent to that of the gold-standard solution. In this study, we develop a novel, highly accurate, verifiable and cost-effective solution to the radiation transport equation on a Eulerian domain using an energy balance approach where we derive an expression for the absorption coefficient as a function of particle and Eulerian mesh sizes. We apply our new method to Poisson and turbulent particle clouds that violate all the BB-law constraints and show that the solution is converging upon the mesh refinement, and eventually, we recover the same gold-standard solution for a much cheaper computational cost.
Polynomial range estimation as a troubled-cell indicator for high-order methods
Two troubled-cell indicators based on polynomial range estimation methods are used to flag cells that may violate positivity constraints. One method uses interval extension, and the second uses the range enclosure property of the Bernstein polynomial basis. Furthermore, both methods reduce compute time for the positivity preserver by limiting its application to a subset of cells. The Bernstein polynomial method remains effective as the problem dimensionality increases. Interval extension applied to the internal energy equation permits the use of the troubled-cell indicators for rational functions, though performance suffers compared to directly applying the indicators to polynomial functions.
An atomistic theory of nucleation: Self-organization via non-equilibrium work and fluctuations
For this work, insights from non-equilibrium statistical mechanics, highlighting the role of work and fluctuations at the microscale, are applied toward the development of a fundamental, rigorous and purely atomistic theory of nucleation. Nanoscale fluctuations in order, density and heat influence the local nucleation rate by orders of magnitude, necessitating their inclusion through a modern approach. Coarse-graining over the underlying Hamiltonian dynamics allows derivation of a microscale expression for the nucleation rate in terms of a classical path integral over far from equilibrium trajectories and their associated work. Second law violating states at the microscale, as found from the dynamics of small critical nucleation clusters, contribute exponentially to the observable macroscale nucleation rate.
Comparison of near-background concentrations of Argon-37 and Xenon-133 in the Atmosphere
Radioisotopes of noble gases xenon and argon can be important indicators of underground nuclear explosions. The Comprehensive Nuclear-Test-Ban Treaty (CTBT) includes monitoring capabilities to identify potential nuclear explosions conducted in violation of the CTBT. This monitoring currently focuses on measurement of the xenon isotopes 133Xe, 135Xe, 131mXe and 133mXe. However, it is predicted that within 100 days of an underground nuclear explosion (UNE) 37Ar would be released to the atmosphere at higher concentrations than xenon isotopes (Haas et al. 2010) and with a higher signal to background ratio, depending on the radioxenon background levels. Therefore, inclusion of 37Ar measurement capabilities at atmospheric International Monitoring System (IMS) stations may represent an improvement in the capability to detect a nuclear explosion. At an IMS station, it could be difficult to determine what constitutes an elevated concentration as a result of an UNE without first understanding the expected range of background concentrations. This work describes our analysis of atmospheric samples for 37Ar to evaluate the range of background concentrations. Samples were collected at multiple locations, with approximately half coming from a sampler co-located with an IMS xenon monitoring station (RN75). The range of 37Ar concentrations measured in atmospheric air samples was relatively narrow; for samples considered detectable, the minimum and maximum measured concentrations were 0.56 and 2.3 mBq/m3, respectively. Comparison of 37Ar and 133Xe concentrations measured at the IMS station indicated some correlation between the measured concentrations. The results presented here demonstrate the capability to detect background concentrations of 37Ar in atmospheric air and provide a basis for potential implementation of 37Ar monitoring at IMS stations.
Time-Domain Vortex Induced Vibration Modeling of Reference Dynamic Power Cable for the Gulf of Maine
Vortex-induced vibration (VIV) is a phenomenon known to decrease the fatigue life of dynamic power cables used in floating offshore wind systems through increased bending loading cycles. A variety of modeling solutions have been proposed to study this, but none enable fully-coupled simulations leveraging open-source tools. To address this, a time-domain VIV model has been successfully implemented into the open-source mooring dynamics model MoorDyn. The semi-empirical VIV model describes a lift force acting on segments of a flexible cylinder, making it well suited for implementation into MoorDyn. The new capability successfully predicted the frequencies and magnitudes of strain in both steady and oscillating flows when compared to the original validation results and matched peak spectral response when compared to experimental results, verifying successful implementation. MoorDyn with VIV was then leveraged to simulate a 15 MW floating turbine with a dynamic power cable in the Gulf of Maine using OpenFAST for six return periods of combined wind, wave, and current conditions. Increased curvature and tension fluctuations in higher flow speeds were observed when simulating VIV. Maximum tensions and curvatures also increased, with the 500 year conditions violating the curvature factor of safety of 2.0. These results highlight the importance of considering VIV when designing dynamic power cables for the Gulf of Maine, and indicate that cost effective mitigation strategies should be explored for the area. They also demonstrate the utility of this new modeling capability, which provides the first open-source tool for fully coupled time-domain floating offshore wind simulations with dynamic power cable VIV.
Ab initio study and assignment of electronic states in molecular RaCl
Radium compounds have attracted recently considerable attention due to both development of experimental techniques for high-precision laser spectroscopy of molecules with short-lived nuclei and amenability of certain radium compounds for direct cooling with lasers. Currently, radium monofluoride (RaF) is one of the most studied molecules among the radium compounds, both theoretically and recently also experimentally. Complementary studies of further diatomic radium derivatives are highly desired to assess the influence of chemical substitution on diverse molecular parameters, especially on those connected with laser cooling, such as vibronic transition probabilities, and those related to violations of fundamental symmetries. Here in this article high-precision ab initio studies of electronic and vibronic levels of diatomic radium monochloride (RaCl) are presented. Recently developed approaches for treating electronic correlation with Fock-space coupled cluster methods are applied for this purpose. Theoretical results are compared to an early experimental investigation by Lagerqvist and used to partially reassign the experimentally observed transitions and molecular electronic levels of RaCl. Effective constants of $\mathscr{P}$-odd hyperfine interaction W a and $\mathscr{P, T}$-odd scalar-pseudoscalar nucleus-electron interaction Ws in the ground electronic state of RaCl are estimated within the framework of a quasirelativistic Zeroth-Order Regular Approximation approach and compared to parameters in RaF and RaOH.
Demonstration of o-Ps detection with a cylindrical array of NaI detectors
Ortho-positronium (o-Ps) is a bound state of an electron and positron that can provide a means of searching for new physics. We searched and identified o-Ps in our cylindrical array of NaI(Tl) scintillators, developing a technique for tagging on o-Ps with any one of the 24 bars in the array. With high angular efficiency and this new technique, such an experimental setup could be used in a fundamental symmetries search. In particular, we explore the possibility of using such a setup with a search for CP- or CPT-violation in o-Ps.
Simulation of the spatial shift in detector response for polarized protons within a calorimeter
Measurement of the helicity dependent elastic electron-proton scattering cross section provides a key means of investigating parity violation within the proton. However, such measurements exhibit potential instrumental effects associated with the detection of polarized recoiled protons. In particular, spin-orbit interactions within a massive detector induce a systematic spatial shift in the detector signal. In this study, we determine the size of this shift using the Geant4 simulation toolkit. For a typical hadron calorimeter, we found a polarization dependent shift on the order of 0.01-0.1 mm, multiple orders of magnitude smaller than the typical spatial resolution seen in hadronic calorimeters. Additionally, we provide the custom modifications required of the Geant4 source code to implement the quasi-elastic scattering of polarized protons incident on nuclei in the detector. The modifications are readily extendable to generic matter sources, and can be used for the study of additional spin dependent observables in Geant4.
Precision Møller polarimetry for PREX-2 and CREX
The PREX-2 and CREX experiments in Hall A at Jefferson Lab are precision measurements of parity violating elastic electron scattering from complex nuclei. One requirement was that the incident electron beam polarization, typically , be known with 1% precision. We commissioned and operated a Møller polarimeter on the beam line that exceeds this requirement, achieving a precision of 0.89% for PREX-2, and 0.85% for CREX. The uncertainty is purely systematic, accumulated from several different sources, but dominated by our knowledge of the target polarization. Our analysis also demonstrates the need for accurate atomic wave functions in order to correct for the Levchuk Effect. We describe the details of the polarimeter operation and analysis, as well as (for CREX) a comparison to results from a different polarimeter based on Compton scattering.
Measured and projected beam backgrounds in the Belle II experiment at the SuperKEKB collider
The Belle II experiment at the SuperKEKB electron–positron collider aims to collect an unprecedented data set of 50 ab -1 to study CP -violation in the B -meson system and to search for Physics beyond the Standard Model. SuperKEKB is already the world’s highest-luminosity collider. In order to collect the planned data set within approximately one decade, the target is to reach a peak luminosity of 6 x 10 35 cm −2 s −1 by further increasing the beam currents and reducing the beam size at the interaction point by squeezing the betatron function down to β $^*_y$ = 0.3 mm. To ensure detector longevity and maintain good reconstruction performance, beam backgrounds must remain well controlled. We report on current background rates in Belle II and compare these against simulation. We find that a number of recent refinements have significantly improved the background simulation accuracy. Finally, we estimate the safety margins going forward. We predict that backgrounds should remain high but acceptable until a luminosity of at least 2.8 x 10 35 cm −2 s −1 is reached for β $^*_y$ = 0.6 mm. At this point, the most vulnerable Belle II detectors, the Time-of-Propagation (TOP) particle identification system and the Central Drift Chamber (CDC), have predicted background hit rates from single-beam and luminosity backgrounds that add up to approximately half of the maximum acceptable rates.
Characterization of the optical model of the T2K 3D segmented plastic scintillator detector unit cube
The magnetized near detector (ND280) of the T2K long-baseline neutrino oscillation experiment has been recently upgraded aiming to satisfy the requirement of reducing the systematic uncertainty from measuring the neutrino–nucleus interaction cross section, which is the largest systematic uncertainty in the search for leptonic charge-parity symmetry violation. A key component of the upgrade is SuperFGD, a 3D segmented plastic scintillator detector made of approximately 2,000,000 optically-isolated 1 cm 3 cubes. The SuperFGD cube unit shows promising optical performance, including a high light yield of about 40 photoelectrons (p.e.) per channel, a low cube-to-cube crosstalk rate below 3%, and a sub-nanosecond time resolution of 0.96 ns. By combining tracking and stopping power measurements of final state particles, this novel detector enables precise 3D-imaging of GeV neutrino interactions with reduced systematic uncertainties. A detailed Geant4 based optical simulation of the SuperFGD building block, i.e. a plastic scintillating cube read out by three wavelength shifting fibers, has been developed and validated with the different datasets collected in various beam tests. In this manuscript the description of the optical model as well as the comparison with data are reported.
Optimizing a detection system for fissile material in nuclear disarmament verification
In arms control treaties, verification plays a crucial role in detecting non-compliance, deterring future violations, and building trust between state parties. Neutron interrogation that induces fission reactions in fissile isotopes and measures the resulting fission neutrons, could be employed for this purpose. This study aims to develop a system which can determine the presence of fissile material while intrinsically protecting information. In this paper, we focus on optimizing the system for discriminating between an enriched uranium block from a depleted uranium (DU) block. The system was built and we report on benchmark measurements with DU and 16% enriched uranium blocks. Furthermore, the Excalibur (Experiment for Calibration with Uranium) neutron source, a neutron spectrometer (redBubble Technology Industries (BTI) N-Probe), and superheated droplet detectors were used for these measurements. MCNP simulations provided insights into detector responses to fissile materials with varying isotopic compositions, confirming that the system functioned as designed.
Isotope dependence of muon-to-electron conversion
The lepton-flavor-violating conversion of a muon into an electron in the field of a nucleus is one of the most sensitive probes of physics beyond the Standard Model and the experiments Mu2e, COMET, and DeeMe will explore uncharted terrain in the near future. An observation of this μ - → e - conversion process opens up the possibility to distinguish the underlying operator or new-physics model by exploiting the target-nucleus dependence of the conversion rate. To facilitate the choice of optimal targets we provide μ - → e - conversion rates for all stable isotopes and estimate nuclear-structure uncertainties. Our results enable studies of mixed or enriched target materials that are particularly promising for distinguishing scenarios in which the muon converts either on protons or neutrons.
Mass distributions of dijet resonances from excited quarks at proton-proton colliders
We study the expected experimental mass distributions of dijet resonances from excited quarks in proton-proton collisions at energies $\sqrt{s}$ = 13, 14, 27, 100, 300, and 500 TeV. We explore in detail the expected shapes at both the generator and experimental levels, and identify within the distributions the effects of the excited quark natural width, parton momentum distributions of the proton, radiation, and experimental resolution. We present both differential and cumulative probability distributions as a function of dijet mass, and the signal acceptance of a window in dijet mass centered on each resonance. We find that for a range of resonance masses, between 10% and 50% of $\sqrt{s}$, the dijet mass distributions and window acceptance are practically universal, approximately invariant under changes in resonance mass and $\sqrt{s}$. This universality is violated when the resonance mass reaches 60% of $\sqrt{s}$, because the steepness of the parton momentum distributions of the proton produces a significant tail at low dijet mass. This work supports our Snowmass 2021 study on the sensitivity to dijet resonances at proton-proton colliders.