Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion theory”

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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Size distribution and scattering phase function of aerosol particles retrieved from sky brightness measurements

Ground-based measurements of the solar transmission and sky radiance in a horizontal plane through the Sun are taken in several geographical regions and aerosol types: dust in a desert transition zone in Israel, sulfate particles in Eastern and Western Europe, tropical aerosol in Brazil, and mixed continental/maritime aerosol in California. Stratospheric aerosol was introduced after the eruption of Mount Pinatubo in June 1991. Therefore measurements taken before the eruption are used to analyze the properties of tropospheric aerosol; measurements from 1992 are also used to detect the particle size and concentration of stratospheric aerosol. The measurements are used to retrieve the size distribution and the scattering phase function at large scattering angles of the undisturbed aerosol particles. The retrieved properties represent an average on the entire atmospheric column. A comparison between the retrieved phase function for a scattering angle of 120 deg, with phase function predicted from the retrieved size distribution, is used to test the assumption of particle homogeneity and sphericity in radiative transfer models (Mie theory). The effect was found to be small (20% +/- 15%). For the stratospheric aerosol (sulfates), as expected, the phase function was very well predicted using the Mie theory. A model with a power law distribution, based on the spectral dependence of the optical thickness, alpha, cannot estimate accurately the phase function (up to 50% error for lambda = 0.87 microns). Before the Pinatubo eruption the ratio between the volumes of sulfate and coarse particles was very well correlated with alpha. The Pinatubo stratospheric aerosol destroyed this correlation. The aerosol optical properties are compared with analysis of the size, shape, and composition of the individual particles by electron microscopy of in situ samples. The measured volume size distribution before the injection of stratospheric aerosol consistently show two modes, sulfate particles with r(sub m) less than 0.2 microns and coarse paritcles with r(sub m) greater than 0.7 microns. The 'window' in the tropospheric aerosol in this radius range was used to observe a stable stratospheric aerosol in 1992, with r(sub m) approximately 0.5 microns. A combination of such optical thickness and sky measurements can be used to assess the direct forcing and the climatic impact of aerosol. Systematic inversion for the key aerosol types (sulfates, smoke, dust, and maritime aerosol) of the size distribution and phase function can give the relationship between the aerosol physical and optical properties that can be used to compute the radiative forcing. This forcing can be validated in dedicated field experiments.

Kaufman, Y. J.↗
NASA NTRS2011

Prominences: The Key to Understanding Solar Activity

Prominences are spectacular manifestations of both quiescent and eruptive solar activity. The largest examples can be seen with the naked eye during eclipses, making prominences among the first solar features to be described and catalogued. Steady improvements in temporal and spatial resolution from both ground- and space-based instruments have led us to recognize how complex and dynamic these majestic structures really are. Their distinguishing characteristics - cool knots and threads suspended in the hot corona, alignment along inversion lines in the photospheric magnetic field within highly sheared filament channels, and a tendency to disappear through eruption - offer vital clues as to their origin and dynamic evolution. Interpreting these clues has proven to be contentious, however, leading to fundamentally different models that address the basic questions: What is the magnetic structure supporting prominences, and how does so much cool, dense plasma appear in the corona? Despite centuries of increasingly detailed observations, the magnetic and plasma structures in prominences are poorly known. Routine measurements of the vector magnetic field in and around prominences have become possible only recently, while long-term monitoring of the underlying filament-channel formation process also remains scarce. The process responsible for prominence mass is equally difficult to establish, although we have long known that the chromosphere is the only plausible source. As I will discuss, however, the motions and locations of prominence material can be used to trace the coronal field, thus defining the magnetic origins of solar eruptions. A combination of observations, theory, and numerical modeling must be used to determine whether any of the competing theories accurately represents the physics of prominences. I will discuss the criteria for a successful prominence model, compare the leading models, and present in detail one promising, comprehensive scenario for prominence formation and evolution that could answer the two questions posed above.

Karpen, Judy T.↗
DOE OSTI2023

Uncertainty quantification for neutrino opacities in core-collapse supernovae and neutron star mergers

We perform an extensive study of the correlations between the neutrino-nucleon inverse mean free paths (IMFPs) and the underlying equation of states (EoSs). Strong interaction uncertainties in the neutrino mean free path are investigated in different density regimes. The nucleon effective mass, the nucleon chemical potentials, and the residual interactions in the medium play an important role in determining neutrino-nucleon interactions in a density-dependent manner. Here we study how the above quantities are constrained by an EoS consistent with (i) nuclear mass measurements, (ii) proton-proton scattering phase shifts, and (iii) neutron star observations. We then study the uncertainties of both the charged current and the neutral current neutrino-nucleon inverse mean free paths due to the variation of these quantities, using the Hartree-Fock+random phase approximation method. Finally, we calculate the Pearson correlation coefficients between (i) the EoS-based quantities and the EoS-based quantities; (ii) the EoS-based quantities and the IMFPs; (iii) the IMFPs and the IMFPs. We find a strong impact of residual interactions on neutrino opacity in the spin and spin-isospin channels, which are not well constrained by current nuclear modelings.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

The Use of Vegetation For Estimating Broken-Cloud Optical Properties From Surface Measurements

It is known that because of complex three-dimensional (3D) radiative effects of broken clouds, the retrieval of cloud optical properties from upward measurements based on a one-dimensional (1D) inversion technique almost surely fails. To remove radiative effects of 3D cloud structure, we have developed a new technique that retrieves cloud optical thickness for broken clouds above green vegetation from simultaneous surface measurements in the VIS and Near Infrared (NIR) spectral regions. The theoretical basis of the method is the very different spectral behavior of cloud liquid water drops and green vegetation. For example, cloud optical properties, and hence cloud reflectivities, change little between 650 and 860 nm, while the vegetated surface albedo changes from 0.05 to 0.5 between the same two wavelengths. This spectral contrast in surface albedo suggests using ground measurements at both wavelengths not independently, but as an algebraic combination (a spectral index). For a spectral band in the NIR region, the green vegetation acts as a powerful reflector that "illuminates" horizontally inhomogeneous clouds from below. This provides the extra information needed to largely remove the 3D radiative effects, especially in the case of broken clouds; this in turn allows the retrieval of cloud optical depth using traditional 1D radiative transfer theory. This approach is similar to the so-called Green's function problem for radiative transfer where a laser beam illuminates clouds and the resulting "spot-size" of the reflected light around the beam characterizes cloud properties. We generalize Green's function theory to surf ace-cloud interaction and develop new spectral indices from which broken-cloud optical depth can be retrieved.

Marshak, A.↗
DOE OSTI2021

Semiclassical shell-structure micro-macroscopic approach for the level density

Level density ρ(E,A) is derived for a one-component nucleon system with a given energy E and particle number A within the mean-field semiclassical periodic-orbit theory beyond the saddle-point method of the Fermi gas model. We obtain ρ∝I ν (S)/S ν , with I ν (S) being the modified Bessel function of the entropy S. Within the micro-macro-canonical approximation (MMA), for a small thermal excitation energy U, with respect to rotational excitations E rot , one obtains ν = 3/2 for ρ(E,A). In the case of excitation energy U larger than E rot but smaller than the neutron separation energy, one finds a larger value of ν = 5/2. A role of the fixed spin variables for rotating nuclei is discussed. The MMA level density ρ reaches the well-known grand-canonical ensemble limit (Fermi gas asymptote) for large S related to large excitation energies, and also reaches the finite micro-canonical limit for small combinatorial entropy S at low excitation energies (the constant “temperature” model). Fitting the ρ(E,A) of the MMA to the experimental data for low excitation energies, taking into account shell and, qualitatively, pairing effects, one obtains for the inverse level density parameter K a value which differs essentially from that parameter derived from data on neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2010

Thermal Performance of Cryogenic Multilayer Insulation at Various Layer Spacings

Multilayer insulation (MLI) has been shown to be the best performing cryogenic insulation system at high vacuum (less that 10 (exp 3) torr), and is widely used on spaceflight vehicles. Over the past 50 years, many investigations into MLI have yielded a general understanding of the many variables that are associated with MLI. MLI has been shown to be a function of variables such as warm boundary temperature, the number of reflector layers, and the spacer material in between reflectors, the interstitial gas pressure and the interstitial gas. Since the conduction between reflectors increases with the thickness of the spacer material, yet the radiation heat transfer is inversely proportional to the number of layers, it stands to reason that the thermal performance of MLI is a function of the number of layers per thickness, or layer density. Empirical equations that were derived based on some of the early tests showed that the conduction term was proportional to the layer density to a power. This power depended on the material combination and was determined by empirical test data. Many authors have graphically shown such optimal layer density, but none have provided any data at such low densities, or any method of determining this density. Keller, Cunnington, and Glassford showed MLI thermal performance as a function of layer density of high layer densities, but they didn't show a minimal layer density or any data below the supposed optimal layer density. However, it was recently discovered that by manipulating the derived empirical equations and taking a derivative with respect to layer density yields a solution for on optimal layer density. Various manufacturers have begun manufacturing MLI at densities below the optimal density. They began this based on the theory that increasing the distance between layers lowered the conductive heat transfer and they had no limitations on volume. By modifying the circumference of these blankets, the layer density can easily be varied. The simplest method of determining the thermal performance of MLI at cryogenic temperature is by boil-off calorimetry. Several blankets were procured and tested at various layer densities at the Cryogenics Test Laboratory at Kennedy Space Center. The densities that the blankets were tested over covered a wide range of layer densities including the analytical minimum. Several of the blankets were tested at the same insulation thickness while changing the layer density (thus a different number of reflector layers). Optimizing the layer density of multilayer insulation systems for heat transfer would remove a layer density from the complex method of designing such insulation systems. Additional testing was performed at various warm boundary temperatures and pressures. The testing and analysis was performed to simplify the analysis of cryogenic thermal insulation systems. This research was funded by the National Aeronautics and Space Administration's Exploration Technology Development Program's Cryogenic Fluid Management Project

Johnson, Wesley Louis↗
NASA NTRS1999

A Spherical Earth Solution for TOA Lightning Location Retrieval

With the advent of high-speed digital computers, the method of chi square minimization is a highly practical means for analyzing a wide variety of (otherwise intractable) nonlinear inversion problems in applied mathematical physics. Little thought or effort is required to apply the chi square method to obtain quick and reasonable estimates of a solution, and the method offers a means to assess retrieval errors. Because the method is simple and practical it is sometimes hastily applied to problems that can be solved by formal analytic or quasi-analytic means. Presently, Global Atmospherics Inc. (GAI) finds the minimum of a chi square function to analyze time-of-arrival (TOA) and magnetic bearing data derived from the National Lightning Detection Network (NLDN); ellipsoidal Earth geometry is assumed. An analytic solution to this problem has not yet been attained, but the consideration and solving of less general problem statements might eventually lead to a final solution. In the present study, the problem of retrieving lightning ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered. It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (such as the chi square method mentioned above). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A quasi-analytic extension of the spherical Earth solution for an oblate spheroidal Earth geometry is considered in a related study. The incorporation of magnetic bearing information into these analytic solutions would lead to a general and elegant analytic retrieval scheme that would most likely replace the chi square estimation theory currently employed by Global Atmospherics Inc. (GAI).

Koshak, William J,↗
DOE OSTI2021

An adaptive approach to machine learning for compact particle accelerators

Abstract Machine learning (ML) tools are able to learn relationships between the inputs and outputs of large complex systems directly from data. However, for time-varying systems, the predictive capabilities of ML tools degrade if the systems are no longer accurately represented by the data with which the ML models were trained. For complex systems, re-training is only possible if the changes are slow relative to the rate at which large numbers of new input-output training data can be non-invasively recorded. In this work, we present an approach to deep learning for time-varying systems that does not require re-training, but uses instead an adaptive feedback in the architecture of deep convolutional neural networks (CNN). The feedback is based only on available system output measurements and is applied in the encoded low-dimensional dense layers of the encoder-decoder CNNs. First, we develop an inverse model of a complex accelerator system to map output beam measurements to input beam distributions, while both the accelerator components and the unknown input beam distribution vary rapidly with time. We then demonstrate our method on experimental measurements of the input and output beam distributions of the HiRES ultra-fast electron diffraction (UED) beam line at Lawrence Berkeley National Laboratory, and showcase its ability for automatic tracking of the time varying photocathode quantum efficiency map. Our method can be successfully used to aid both physics and ML-based surrogate online models to provide non-invasive beam diagnostics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Scientific Visualization and Computational Science: Natural Partners

Scientific visualization is developing rapidly, stimulated by computational science, which is gaining acceptance as a third alternative to theory and experiment. Computational science is based on numerical simulations of mathematical models derived from theory. But each individual simulation is like a hypothetical experiment; initial conditions are specified, and the result is a record of the observed conditions. Experiments can be simulated for situations that can not really be created or controlled. Results impossible to measure can be computed.. Even for observable values, computed samples are typically much denser. Numerical simulations also extend scientific exploration where the mathematics is analytically intractable. Numerical simulations are used to study phenomena from subatomic to intergalactic scales and from abstract mathematical structures to pragmatic engineering of everyday objects. But computational science methods would be almost useless without visualization. The obvious reason is that the huge amounts of data produced require the high bandwidth of the human visual system, and interactivity adds to the power. Visualization systems also provide a single context for all the activities involved from debugging the simulations, to exploring the data, to communicating the results. Most of the presentations today have their roots in image processing, where the fundamental task is: Given an image, extract information about the scene. Visualization has developed from computer graphics, and the inverse task: Given a scene description, make an image. Visualization extends the graphics paradigm by expanding the possible input. The goal is still to produce images; the difficulty is that the input is not a scene description displayable by standard graphics methods. Visualization techniques must either transform the data into a scene description or extend graphics techniques to display this odd input. Computational science is a fertile field for visualization research because the results vary so widely and include things that have no known appearance. The amount of data creates additional challenges for both hardware and software systems. Evaluations of visualization should ultimately reflect the insight gained into the scientific phenomena. So making good visualizations requires consideration of characteristics of the user and the purpose of the visualization. Knowledge about human perception and graphic design is also relevant. It is this breadth of knowledge that stimulates proposals for multidisciplinary visualization teams and intelligent visualization assistant software. Visualization is an immature field, but computational science is stimulating research on a broad front.

Uselton, Samuel P.↗
NASA NTRS2000

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗
DOE OSTI2026

Quantum state preparation with resolution refinement

We introduce a method called resolution refinement that allows one to bootstrap eigenstate preparation on a quantum computer. We first prepare an eigenstate of a low-resolution Hamiltonian using any method of choice. The eigenstate is then lifted to higher resolution and adiabatically evolved to produce the corresponding eigenstate of a higher-fidelity Hamiltonian. We give examples of resolution refinement applied to both single-particle basis states as well as a spatial lattice grid. For basis refinement, we compute few-body ground states of the Busch model for interacting particles in a harmonic trap in one dimension. For lattice refinement, we compute Hartree-Fock nuclear states for a central Woods-Saxon potential in three dimensions, and we compute bound states and continuum states in a multi-species Hubbard model of fermions in one dimension. In all cases, the method is efficient and requires an adiabatic evolution time that scales with the inverse of the energy gap times the square root of the system size. We show that this very favorable scaling arises from the fact that resolution refinement does not make large changes to the structure or energies of the low-energy eigenstates.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2000

Freja Studies of the Current-Voltage Relation in Substorm-Related Events

Field-aligned currents and electrostatic potentials play important roles in the coupling between the magnetosphere and the ionosphere. If one assumes that the ionosphere-magnetosphere potential difference is mainly due to the mirror force, one can use the single particle adiabatic kinetic theory to describe the system. From this theory, a linear relationship j(sub II) = KV between field-aligned current density j(sub II) and potential drop V along the same field line can be derived, provided that the potential drop is not too large and not too small. With rare exceptions, observational tests of this relation have mainly concentrated on quiet magnetospheric situations, with acceleration voltages V approx. less than 5 kV. Here we use observations from the Freja satellite of precipitating auroral electrons at 1.700 km altitude to study substorm related events, with acceleration voltages up to 20 keV. The observations are found to be consistent with a linear current-voltage relation even i n these conditions, although with values of the field aligned K lower than previously reported (1-5 x 10(exp 11 S/sq m). This can be explained by lower densities and higher characteristic electron energies in the magnetospheric source region of the precipitating electrons. We analyze the data by three different methods, which are all found to be in general agreement. The results are in agreement with a previous study, where the spectra of precipitating electrons --were indirectly inferred by inversion of data from the EISCAT incoherent scatter radar, thereby validating the use of radar data for studies of auroral electrons. Comparisons with previous studies are made, emphasizing the dependence of the results on the type of auroral structure and magnetospheric conditions.

Olsson, A.↗
DOE OSTI2021

When and why PINNs fail to train: A neural tangent kernel perspective

Physics-informed neural networks (PINNs) have lately received great attention thanks to their flexibility in tackling a wide range of forward and inverse problems involving partial differential equations. However, despite their noticeable empirical success, little is known about how such constrained neural networks behave during their training via gradient descent. More importantly, even less is known about why such models sometimes fail to train at all. Here in this work, we aim to investigate these questions through the lens of the Neural Tangent Kernel (NTK); a kernel that captures the behavior of fully-connected neural networks in the infinite width limit during training via gradient descent. Specifically, we derive the NTK of PINNs and prove that, under appropriate conditions, it converges to a deterministic kernel that stays constant during training in the infinite-width limit. This allows us to analyze the training dynamics of PINNs through the lens of their limiting NTK and find a remarkable discrepancy in the convergence rate of the different loss components contributing to the total training error. To address this fundamental pathology, we propose a novel gradient descent algorithm that utilizes the eigenvalues of the NTK to adaptively calibrate the convergence rate of the total training error. Finally, we perform a series of numerical experiments to verify the correctness of our theory and the practical effectiveness of the proposed algorithms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗
DOE OSTI2022

High-resolution CH stretch spectroscopy of jet-cooled cyclopentyl radical: First insights into equilibrium structure, out-of-plane puckering, and IVR dynamics

First, high-resolution sub-Doppler infrared spectroscopic results for cyclopentyl radical (C 5 H 9 ) are reported on the α-CH stretch fundamental with suppression of spectral congestion achieved by adiabatic cooling to T rot ≈ 19(4) K in a slit jet expansion. Surprisingly, cyclopentyl radical exhibits a rotationally assignable infrared spectrum, despite 3N – 6 = 36 vibrational modes and an upper vibrational state density (ρ ≈ 40–90 #/cm –1 ) in the critical regime (ρ ≈ 100 #/cm –1 ) necessary for onset of intramolecular vibrational relaxation (IVR) dynamics. Such high-resolution data for cyclopentyl radical permit detailed fits to a rigid-rotor asymmetric top Hamiltonian, initial structural information for ground and vibrationally excited states, and opportunities for detailed comparison with theoretical predictions. Specifically, high level ab initio calculations at the coupled-cluster singles, doubles, and perturbative triples (CCSD(T))/ANO0, 1 level are used to calculate an out-of-plane bending potential, which reveals a C 2 symmetry double minimum 1D energy surface over a C 2v transition state. Here, the inversion barrier [V barrier ≈ 3.7(1) kcal/mol] is much larger than the effective moment of inertia for out-of-plane bending, resulting in localization of the cyclopentyl wavefunction near its C 2 symmetry equilibrium geometry and tunneling splittings for the ground state too small (<1 MHz) to be resolved under sub-Doppler slit jet conditions. The persistence of fully resolved high-resolution infrared spectroscopy for such large cyclic polyatomic radicals at high vibrational state densities suggests a “deceleration” of IVR for a cycloalkane ring topology, much as low frequency torsion/methyl rotation degrees of freedom have demonstrated a corresponding “acceleration” of IVR processes in linear hydrocarbons.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Putting all the X in one basket: Updated X-ray constraints on sub-GeV Dark Matter

Sub-GeV dark matter particles can annihilate or decay producing e ± pairs which upscatter the low-energy photon fields in the Galaxy and generate an X-ray emission (via the Inverse Compton effect). Here, using X-ray data from X MM -N EWTON , I NTEGRAL , N U S TAR and S UZAKU , we derive new constraints on this class of dark matter (DM). In the annihilation case, our new bounds are the strongest available for DM masses above 180 MeV, reaching $\langle$σv$\rangle$ ≲ 10 -28 cm 3 /s for m DM ≃ 1 GeV. In the decay case, our bounds are the strongest to date over a large fraction of the considered mass range, constraining τ ≳ 10 28 s for m DM ≃ 1 GeV and improving by up to 3 orders of magnitude upon existing limits.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Analyzing rotational bands in odd-mass nuclei using effective field theory and Bayesian methods

We recently developed an Effective Field Theory (EFT) for rotational bands in odd-mass nuclei. Here we use EFT expressions to perform a Bayesian analysis of data on the rotational energy levels of 99 Tc, 155,157 Gd, 159 Dy, 167,169 Er, 167,169 Tm, 183 W, 235 U and 239 Pu. The error model in our Bayesian analysis includes both experimental and EFT truncation uncertainties. It also accounts for the fact that low-energy constants (LECs) at even and odd orders are expected to have different sizes. We use Markov Chain Monte Carlo (MCMC) sampling to explore the joint posterior of the EFT and error-model parameters and show both the LECs and the breakdown scale can be reliably determined. We extract the LECs up to fourth order in the EFT and find that, provided we correctly account for EFT truncation errors in our likelihood, results for lower-order LECs are stable as we go to higher orders. LEC results are also stable with respect to the addition of higher-energy data. We extract the expansion parameter for all the nuclei listed above and find a clear correlation between the extracted and the expected value of the inverse breakdown scale, W, based on the single-particle and vibrational energy scales. However, the W that actually determines the convergence of the EFT expansion is markedly smaller than would be naively expected based on those scales.

EFT↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗