Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inversion constraints”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Lateral variation in upper mantle temperature and composition beneath mid-ocean ridges inferred from shear-wave propagation, geoid, and bathymetry

Resolution of both the extent and mechanism of lateral heterogeneity in the upper mantle constraints the nature and scales of mantle convection. Oceanic regions are of particular interest as they are likely to provide the closest glimpse at the patterns of temperature anomalies and convective flow in the upper mantle because of their young age and simple crustal structure relative to continental regions. Lateral variations were determined in the seismic velocity and attenuation structure of the lithosphere and astenosphere beneath the oceans, and these seismological observations were combined with the data and theory of geoid and bathymetry anomalies in order to test and improve current models for seafloor spreading and mantle convection. Variations were determined in mantle properties on a scale of about 1000 km, comparable to the thickness of the upper mantle. Seismic velocity, geoid, and bathymetry anomalies are all sensitive to variations in upper mantle density, and inversions were formulated to combine quantitatively these different data and to search for a common origin. Variations in mantle density can be either of thermal or compositional origin and are related to mantle convection or differentiation.

Sheehan, Anne Francis↗
DOE OSTI2022

Physics Informed Neural Networks as Computational Physics Emulators

This report is a brief overview and evaluation of Physics Informed Neural Networks (PINNs). Karniadakis and co-workers, e.g., Karniadakis et al. (2021) assert that the PINNs approach integrates seamlessly both data and mathematical physics models, even in partially understood, uncertain and high-dimensional contexts. They further claim that PINNs are effective and efficient for ill-posed and inverse problems, and when combined with domain decomposition, are scalable to large problems and a tool to discover hidden physics. While they demonstrate the capabilities in specific academic instances, their overarching claims about PINNs seem to be an overstatement, at least at the current time. We have briefly considered a few of the limitations of PINNs in this investigation. It is not clear to us if the PINNs approach can ever be competitive with approaches that use specialized algorithms to achieve high-accuracy solutions of governing equations and other techniques that can combine observational data with such solutions. As an example of the latter, consistent with the principles of Bayesian inference, data assimilation, or more generally data-model fusion, is a process that fuses observational data typically with a computational model that respects certain constraints such as conservation laws. For example, improvements in observational network combined with data assimilation have been key in improving weather predictions over the past four decades Kalnay (2003).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Measurement of Thomson-scattering spectra with continuous angular resolution (invited)

Here, a novel Thomson-scattering diagnostic with continuous angular resolution over a span of 120° was developed for the characterization of plasmas produced at the Omega Laser Facility. Spectrally resolving light scattered from electron plasma wave features as a function of emission angle provides a means to efficiently probe a large range of plasma frequencies and k vectors. Together, these spectra contain critical constraints on the plasma-physics models used to interpret the data and enable experimental measurements of the electron-velocity distribution function over several orders of magnitude without assumptions about its mathematical form. Major components of the instrument include (1) a reflective collection objective that gathers light over a range of 120° × 12°; (2) a spatial-filter image relay for measurement localization; (3) cylindrical optics for producing a line image of the collection aperture; (4) a transmission grating spectrometer; and (5) a time-gated, image-intensified camera. Thomson-scattered light collected from an ~(50 – μM) 3 volume of plasma is recorded with 0.8-nm spectral and 1° angular resolution. Initial experiments examined the properties of the electron-velocity distribution in gas-jet-produced plasmas in the presence of heating via inverse bremsstrahlung absorption.

47 OTHER INSTRUMENTATION↗
NASA NTRS2008

A New Technique for Retrieval of Tropospheric and Stratospheric Ozone Profiles using Sky Radiance Measurements at Multiple View Angles: Application to a Brewer Spectrometer

This paper describes and applies a new technique for retrieving diurnal variability in tropospheric ozone vertical distribution using ground-based measurements of ultraviolet sky radiances. The measured radiances are obtained by a polarization-insensitive modified Brewer double spectrometer located at Goddard Space Flight Center, in Greenbelt, Maryland, USA. Results demonstrate that the Brewer angular (0-72deg viewing zenith angle) and spectral (303-320 nm) measurements of sky radiance in the solar principal plane provide sufficient information to derive tropospheric ozone diurnal variability. In addition, the Brewer measurements provide stratospheric ozone vertical distributions at least twice per day near sunrise and sunset. Frequent measurements of total column ozone amounts from direct-sun observations are used as constraints in the retrieval. The vertical ozone profile resolution is shown in terms of averaging kernels to yield at least four points in the troposphere-low stratosphere, including good information in Umkehr layer 0 (0-5 km). The focus of this paper is on the derivation of stratospheric and tropospheric ozone profiles using both simulated and measured radiances. We briefly discuss the necessary modifications of the Brewer spectrometer that were used to eliminate instrumental polarization sensitivity so that accurate sky radiances can be obtained in the presence of strong Rayleigh scattering and aerosols. The results demonstrate that including a site-specific and time-dependent aerosol correction, based on Brewer direct-sun observations of aerosol optical thickness, is critical to minimize the sky radiance residuals as a function of observing angle in the optimal estimation inversion algorithm and improve the accuracy of the retrieved ozone profile.

Tzortziou, Maria↗
NASA NTRS2004

Path following using visual odometry for Mars Rover in high-slip environments.

An architecture for autonomous operation of Mars rovers in high slip environments has been designed, implemented, and tested. This architecture is composed of several key technologies that enable the rover to accurately follow a designated path, compensate for slippage, and reach intended goals independent of the terrain over which it is traversing (within the mechanical constraints of the mobility system). These technologies include: visual odometry, full vehicle kinematics, a Kalman filter, and a slop compensation/path follower. Visual odometry tracks distinctive scene features in stereo imagery to estimate rover motion between successively acquired stereo image pairs using a maximum likelihood motion estimation algorithm. The full vehicle kinematics for a rocker-bogie suspension system estimates motion, with a no-slip assumption, by measuring wheel, rates, and rocker, bogie, and steering angles. The Kalman filter merges data from an Inertial Measurement Unit (IMU) and visual odometry. This merged estimate is then compared to the kinematic estimate to determining (tracking into account estimate uncertainties) if and how much slippage has occurred then a slip vector is calculated by differencing the current Kalman filter estimate from the kinematic estimate. This slip vector is then used, in conjunction with the inverse kinematics, to determine the necessary wheel velocities and steering angles to compensate for slip and follow the desired path.

Roumeliotis, Stergios I.↗
NASA NTRS2017

Map Projection Induced Variations in Locations of Polygon Geofence Edges

This Paper under-estimates answers to the following question under various constraints: If a geofencing algorithm uses a map projection to determine whether a position is inside/outside a polygon region, how far outside/inside the polygon can the point be and the algorithm determine that it is inside/outside (the opposite and therefore incorrect answer)? Geofencing systems for unmanned aircraft systems (UAS) often model stay-in and stay-out regions using 2D polygons with minimum and maximum altitudes. The vertices of the polygons are typically input as latitude-longitude pairs, and the edges as paths between adjacent vertices. There are numerous ways to generate these paths, resulting in numerous potential locations for the edges of stay-in and stay-out regions. These paths may be geodesics on a spherical model of the earth or geodesics on the WGS84 reference ellipsoid. In geofencing applications that use map projections, these paths are inverse images of straight lines in the projected plane. This projected plane may be a projection of a spherical earth model onto a tangent plane, called an orthographic projection. Alternatively, it may be a projection where the straight lines in the projected plane correspond to straight lines in the latitudelongitude coordinate system, also called a Plate Carr´ee projection. This paper estimates distances between different edge paths and an oracle path, which is a geodesic on either the spherical earth or the WGS84 ellipsoidal earth. This paper therefore estimates how far apart different edge paths can be rather than comparing their path lengths, which are not considered. Rather, the comparision is between the actual locations of the edges between vertices. For edges drawn using orthographic projections, this maximum distance increases as the distance from the polygon vertices to the projection point increases. For edges drawn using Plate Carr´ee projections, this maximum distance increases as the vertices become further from the equator. Distances between geodesics on a spherical earth and a WGS84 ellipsoidal earth are also analyzed, using the WGS84 ellipsoid as the oracle. Bounds on the 2D distance between a straight line and a great circle path, in an orthographically projected plane rather than on the surface of the earth, have been formally verified in the PVS theorem prover, meaning that they are mathematically correct in the absence of floating point errors.

Neeley, Paula↗
DOE OSTI2020

The Effects of Assumed Source Depth and Shear-Wave Velocity on Moment Tensors Estimated for Small, Contained Chemical Explosions in Granite

ABSTRACT The Source Phenomenology Experiment (SPE-Arizona) included of a series of chemical explosions detonated within a copper mine in Arizona. This study focuses on ground motions from detonations in the copper mine, which are analyzed to assess the uniqueness of the resulting source representation when the source region propagation characteristics have a range of possible models. P-wave velocities are well constrained by refraction data with less constraint of the S-wave velocities. The effects of explosion source depth and VS are assessed with Green’s functions for a range of models in which VP is held constant. Propagation models with a Poisson’s value of 0.25 and a source depth 30–60 m most accurately replicate the data. The explosion was detonated at a centroid depth of 30 m, so trade-offs in depth are demonstrated. The compensated linear vector dipole and explosion components of the Green’s functions convolved with a Mueller–Murphy source function are compared. Both produce significant energy in the 2–12 Hz band, due to surface-wave contributions with no clear depth dependencies above 20 Hz. The range of propagation models is used with the observational data to invert for the frequency-domain moment tensor. Fits to the data from these inversions have cross-correlation values of 0.64, demonstrating effectiveness in replicating the observations with the assumed propagation path effects and resulting source function. Inversions produce horizontal dipoles (Mxx and Myy), roughly half the maximum amplitude of Mzz, consistent with a compensated linear vector dipole source, which is frequency dependent. Denny and Johnson, Mueller–Murphy, Walter and Ford, and the revised Mueller–Murphy source models, parameterized for granite, are compared to the moment tensors. Despite a nonisotropic moment tensor source, the revised Mueller–Murphy isotropic source model best replicates the long-period moments, overshoot, and corner frequency.

Geochemistry & Geophysics↗
NASA NTRS2017

Variability and Quasi-Decadal Changes in the Methane Budget over the Period 2000-2012

Following the recent Global Carbon Project (GCP) synthesis of the decadal methane (CH4) budget over 2000- 2012, we analyse here the same dataset with a focus on quasi-decadal and inter-annual variability in CH4 emissions. The GCP dataset integrates results from top-down studies (exploiting atmospheric observations within an atmospheric inverse-modelling framework) and bottom-up models (including process-based models for estimating land surface emissions and atmospheric chemistry), inventories of anthropogenic emissions, and data-driven approaches. The annual global methane emissions from top-down studies, which by construction match the observed methane growth rate within their uncertainties, all show an increase in total methane emissions over the period 2000-2012, but this increase is not linear over the 13 years. Despite differences between individual studies, the mean emission anomaly of the top-down ensemble shows no significant trend in total methane emissions over the period 2000-2006, during the plateau of atmospheric methane mole fractions, and also over the period 2008-2012, during the renewed atmospheric methane increase. However, the top-down ensemble mean produces an emission shift between 2006 and 2008, leading to 22 [16-32] Tg CH4 yr(exp -1) higher methane emissions over the period 2008-2012 compared to 2002-2006. This emission increase mostly originated from the tropics, with a smaller contribution from mid-latitudes and no significant change from boreal regions. The regional contributions remain uncertain in top-down studies. Tropical South America and South and East Asia seem to contribute the most to the emission increase in the tropics. However, these two regions have only limited atmospheric measurements and remain therefore poorly constrained. The sectorial partitioning of this emission increase between the periods 2002-2006 and 2008-2012 differs from one atmospheric inversion study to another. However, all top-down studies suggest smaller changes in fossil fuel emissions (from oil, gas, and coal industries) compared to the mean of the bottom-up inventories included in this study. This difference is partly driven by a smaller emission change in China from the top-down studies compared to the estimate in the Emission Database for Global Atmospheric Research (EDGARv4.2) inventory, which should be revised to smaller values in a near future. We apply isotopic signatures to the emission changes estimated for individual studies based on five emission sectors and find that for six individual top-down studies (out of eight) the average isotopic signature of the emission changes is not consistent with the observed change in atmospheric 13CH4. However, the partitioning in emission change derived from the ensemble mean is consistent with this isotopic constraint. At the global scale, the top-down ensemble mean suggests that the dominant contribution to the resumed atmospheric CH4 growth after 2006 comes from microbial sources (more from agriculture and waste sectors than from natural wetlands), with an uncertain but smaller contribution from fossil CH4 emissions. In addition, a decrease in biomass burning emissions (in agreement with the biomass burning emission databases) makes the balance of sources consistent with atmospheric 13CH4 observations. In most of the top-down studies included here, OH concentrations are considered constant over the years (seasonal variations but without any inter-annual variability). As a result, the methane loss (in particular through OH oxidation) varies mainly through the change in methane concentrations and not its oxidants. For these reasons, changes in the methane loss could not be properly investigated in this study, although it may play a significant role in the recent atmospheric methane changes as briefly discussed at the end of the paper.

Saunois, Marielle↗
DOE OSTI2020

sMF-BO-2CoGP: A sequential multi-fidelity constrained Bayesian optimization framework for design applications

Bayesian optimization (BO) is an effective surrogate-based method that has been widely used to optimize simulation-based applications. While the traditional Bayesian optimization approach only applies to single-fidelity models, many realistic applications provide multiple levels of fidelity with various levels of computational complexity and predictive capability. In this work, we propose a multi-fidelity Bayesian optimization method for design applications with both known and unknown constraints. The proposed framework, called sMF-BO-2CoGP, is built on a multi-level CoKriging method to predict the objective function. An external binary classifier, which we approximate using a separate CoKriging model, is used to distinguish between feasible and infeasible regions. Finally, the sMF-BO-2CoGP method is demonstrated using a series of analytical examples and a flip-chip application for design optimization to minimize the deformation due to warping under thermal loading conditions.

42 ENGINEERING↗
DOE OSTI2023

The Magnificent Five Images of Supernova Refsdal: Time Delay and Magnification Measurements

In late 2014, four images of supernova (SN) "Refsdal," the first known example of a strongly lensed SN with multiple resolved images, were detected in the MACS J1149 galaxy-cluster field. Following the images' discovery, the SN was predicted to reappear within hundreds of days at a new position ~8'' away in the field. The observed reappearance in late 2015 makes it possible to carry out Refsdal's original proposal to use a multiply imaged SN to measure the Hubble constant H0, since the time delay between appearances should vary inversely with H0. Moreover, the position, brightness, and timing of the reappearance enable a novel test of the blind predictions of galaxy-cluster models, which are typically constrained only by the positions of multiply imaged galaxies. We have developed a new photometry pipeline that uses DOLPHOT to measure the fluxes of the five images of SN Refsdal from difference images. We apply four separate techniques to perform a blind measurement of the relative time delays and magnification ratios between the last image SX and the earlier images S1–S4. We measure the relative time delay of SX–S1 to be 376.0$^{+5.6}_{–5.5}$ days and the relative magnification to be 0.30$^{+0.5}_{–0.3}$. This corresponds to a 1.5% precision on the time delay and 17% precision for the magnification ratios and includes uncertainties due to millilensing and microlensing. In an accompanying paper, we place initial and blind constraints on the value of the Hubble constant.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1989

Manipulation strategies for massive space payloads

Control for the bracing strategy is being examined. It was concluded earlier that trajectory planning must be improved to best achieve the bracing motion. Very interesting results were achieved which enable the inverse dynamics of flexible arms to be calculated for linearized motion in a more efficient manner than previously published. The desired motion of the end point beginning at t=0 and ending at t=t sub f is used to calculate the required torque at the joint. The solution is separated into a causal function that is zero for t is less than 0 and an accusal function which is zero for t is greater than t sub f. A number of alternative end point trajectories were explored in terms of the peak torque required, the amount of anticipatory action, and other issues. The single link case is the immediate subject and an experimental verification of that case is being performed. Modeling with experimental verification of closed chain dynamics continues. Modeling effort has pointed out inaccuracies that result from the choice of numerical techniques used to incorporate the closed chain constraints when modeling our experimental prototype RALF (Robotic Arm Large and Flexible). Results were compared to TREETOPS, a multi body code. The experimental verification work is suggesting new ways to make comparisons with systems having structural linearity and joint and geometric nonlinearity. The generation of inertial forces was studied with a small arm that will damp the large arm's vibration.

Book, Wayne J.↗
NASA NTRSDate not supplied

Multiobjective Multidisciplinary Optimization of Low-Boom Supersonic Transports Using Multifidelity Models

A multidisciplinary optimization (MDO) method has been developed to design a computational fluid dynamics (CFD) based low-boom configuration that can be obtained from a Pareto solution of a low-fidelity multiobjective MDO problem with mission constraints. This paper refines the developed MDO method using multifidelity models for CFD-based multiobjective MDO. The refined MDO method can generate a low-boom configuration that satisfies the mission requirements, has the lowest takeoff gross weight and the longest range for the low-boom mission, trims the low-boom cruise flight with fuel redistributions, and has a reversed equivalent area distribution closely matching a low-boom target with ground noise level below 70 PLdB. The validity of the refined MDO method is demonstrated by a design study of a low-boom supersonic transport that carries 40 passengers, flies a low-boom mission with cruise Mach of 1.7 and range of 3500 nm, and cruises overwater at Mach 1.8 with range of 3882 nm. Moreover, the refined MDO method eliminates the difference between the assumed cruise weight for CFD-based low-boom inverse design optimization and the estimated cruise weight of the optimal inverse design solution with respect to the mission requirements.

multidisciplinary optimization↗
DOE OSTI2021

Peptide nucleic acid Hoogsteen strand linker design for major groove recognition of DNA thymine bases

Sequence-specific targeting of double-stranded DNA and non-coding RNA via triple-helix-forming peptide nucleic acids (PNAs) has attracted considerable attention in therapeutic, diagnostic and nanotechnological fields. An E-base (3-oxo-2,3-dihydropyridazine), attached to the polyamide backbone of a PNA Hoogsteen strand by a side-chain linker molecule, is typically used in the hydrogen bond recognition of the 4-oxo group of thymine and uracil nucleic acid bases in the major groove. We report on the application of quantum chemical computational methods, in conjunction with spatial constraints derived from the experimental structure of a homopyrimidine PNA·DNA-PNA hetero-triplex, to investigate the influence of linker flexibility on binding interactions of the E-base with thymine and uracil bases in geometry-optimised model systems. Hydrogen bond formation between the N2 E-base atom and target pyrimidine base 4-oxo groups in model systems containing a β-alanine linker (J Am Chem Soc 119:11116, 1997) was found to incur significant internal strain energy and the potential disruption of intra-stand aromatic base stacking interactions in an oligomeric context. Additionally, in geometry-optimised model systems containing a 3-trans olefin linker (Bioorg Med Chem Lett 14:1551, 2004) the E-base swung out away from the target pyrimidine bases into the solvent. These findings are in qualitative agreement with calorimetric measurements in hybridisation experiments at T–A and U–A inversion sites. In contrast, calculations on a novel 2-cis olefin linker design indicate that it could permit simultaneous E-base hydrogen bonding with the thymine 4-oxo group, circumvention and solvent screening of the thymine 5-methyl group, and maintenance of triplex intra-stand base stacking interactions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Novel application of neutrinos to evaluate U.S. nuclear weapons performance

There is a growing realization that neutrinos can be used as a diagnostic tool to better understand the inner workings of a nuclear weapon. Robust estimates demonstrate that an Inverse Beta Decay (IBD) neutrino scintillation detector built at the Nevada Test Site with a 1000-ton active target mass at a standoff distance of 500 m would detect thousands of antineutrino events per nuclear test. This would provide less than 4% statistical error on the measured antineutrino rate and 5% error on antineutrino energy. Extrapolating this to an error on the test device explosive yield requires knowledge from evaluated nuclear databases, non-equilibrium fission rates, and assumptions on internal neutron fluxes. Initial calculations demonstrate that the total number of neutrinos emitted per fission in the first 10 3 s after a short pulse of 239 Pu fission is about a factor of two less than that from Pu fissioning under steady state conditions. Furthermore, there are significant energy spectral differences as a function of time after the pulse that must be considered. These and other model dependencies will be discussed in the paper. In the absence of nuclear weapons testing, many of the technical and theoretical challenges of a full nuclear test could be mitigated with a low cost smaller scale 20 ton fiducial mass IBD demonstration detector placed near a pulsed reactor. Potential reactors include the Texas A&M University TRIGA 1 GW–10 ms pulsed facility or the Sandia Annular Core Research Reactor. The short duty cycle and repeatability of pulses would provide critical real environment testing and measurements, which would be valuable for planning a possible real test shot in the future. Furthermore, the antineutrino rate as a function of time data would provide unique constraints on fission databases and model assumptions. Finally, there are impactful science drivers such as sensitive searches for ∼1 eV 2 sterile neutrinos and ∼MeV scale axions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

A Framework for a Supervisory Expert System for Robotic Manipulators with Joint-Position Limits and Joint-Rate Limits

This report addresses the problem of path planning and control of robotic manipulators which have joint-position limits and joint-rate limits. The manipulators move autonomously and carry out variable tasks in a dynamic, unstructured and cluttered environment. The issue considered is whether the robotic manipulator can achieve all its tasks, and if it cannot, the objective is to identify the closest achievable goal. This problem is formalized and systematically solved for generic manipulators by using inverse kinematics and forward kinematics. Inverse kinematics are employed to define the subspace, workspace and constrained workspace, which are then used to identify when a task is not achievable. The closest achievable goal is obtained by determining weights for an optimal control redistribution scheme. These weights are quantified by using forward kinematics. Conditions leading to joint rate limits are identified, in particular it is established that all generic manipulators have singularities at the boundary of their workspace, while some have loci of singularities inside their workspace. Once the manipulator singularity is identified the command redistribution scheme is used to compute the closest achievable Cartesian velocities. Two examples are used to illustrate the use of the algorithm: A three link planar manipulator and the Unimation Puma 560. Implementation of the derived algorithm is effected by using a supervisory expert system to check whether the desired goal lies in the constrained workspace and if not, to evoke the redistribution scheme which determines the constraint relaxation between end effector position and orientation, and then computes optimal gains.

Mutambara, Arthur G. O.↗
DOE OSTI2025

Search for very-short-baseline oscillations of reactor antineutrinos with the SoLid detector

In this paper we report the first scientific result based on antineutrinos emitted from the BR2 reactor at SCK CEN. The SoLid experiment uses a novel type of highly granular detector whose basic detection unit combines two scintillators, polyvinyl toluene (PVT) and Li 6 F : ZnS ( Ag ) , to measure antineutrinos via their inverse- β -decay products. An advantage of PVT is its highly linear response as a function of deposited particle energy. The full-scale detector comprises 12 800 voxels and operates over a very short 6.3–8.9 m baseline from the reactor core. The detector segmentation and its three-dimensional imaging capabilities facilitate the extraction of the positron energy from the rest of the visible energy, allowing the latter to be utilized for signal-background discrimination. We present a result obtained from 280 reactor-on days (55 MW mean power) and 172 reactor-off days, respectively, of live data taking. A total of 29 479 ± 603 (stat) antineutrino candidates have been selected, corresponding to an average rate of 105 events per day and a signal-to-background ratio of 0.27. A search for disappearance of antineutrinos to a sterile state has been conducted using complementary model-dependent frequentist and Bayesian fits, providing constraints on the allowed region of the reactor antineutrino anomaly.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Context. As a result of the high-quality constraints available for the Sun, we are able to carry out detailed combined analyses using neutrino, spectroscopic, and helioseismic observations. These studies lay the ground for future improvements of the key physical components of solar and stellar models because ingredients such as the equation of state, the radiative opacities, or the prescriptions for macroscopic transport processes of chemicals are then used to study other stars in the Universe. Aims. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss the effect on their properties of changes in the micro and macro physical ingredients. Methods. We carried out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged 3D models that were claimed to resolve the solar modelling problem. We compared these models to helioseismic and neutrino constraints. Results. The properties of the solar models are significantly affected by the use of the recent OPLIB opacity tables and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those for which the OP opacities were used. We show that a modification of the temperature gradient just below the base of the convective zone is required to remove the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a localised increase in the opacity of a few percent. Conclusions. We conclude that the existing degeneracies and issues in solar modelling are not removed by using an increase in the solar metallicity, in contradiction to what has been suggested in the recent literature. Therefore, standard solar models cannot be used as an argument for a high-metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Baryon-Interacting Dark Matter: heating dark matter and the emergence of galaxy scaling relations

The empirical scaling relations observed in disk galaxies remain challenging for models of galaxy formation. The most striking among these is the Mass Discrepancy-Acceleration Relation (MDAR), which encodes both a tight baryonic Tully-Fisher relation (BTFR) and the observed diversity of galaxy rotation curves through the central surface density relation (CSDR). Building on our earlier work [1], we propose here that the MDAR is the result of interactions between baryons and ‘Baryon-Interacting Dark Matter’ (BIDM), which heat up the dark matter. Following a bottom-up, hydrodynamical approach, we find that the MDAR follows if: i) the BIDM equation of state approximates that of an ideal gas; ii) the BIDM relaxation time is order the Jeans time; iii) the heating rate is inversely proportional to the BIDM density. Remarkably, under these assumptions the set of hydrodynamical equations together with Poisson’s equation enjoy an anisotropic scaling symmetry. In the BIDM-dominated regime, this gives rise to an enhanced symmetry which fully captures the low-acceleration limit of the MDAR. We then show that, assuming a cored pseudo-isothermal profile at equilibrium, this set of equations gives rise to parameters reproducing the MDAR. Specifically, in the flat part of the rotation curve the asymptotic rotational velocity matches the parametric dependence of the BTFR. Moreover, in the central region of high-surface brightness galaxies, the profile reproduces the CSDR. Finally, by studying the time-dependent approach to equilibrium, we derive a global combination of the BTFR and CSDR, which matches the expectations in low surface-brightness galaxies. The form of the heating rate also makes model-independent predictions for various cosmological observables. We argue that our scenario satisfies existing observational constraints, and, intriguingly, offers a possible explanation to the EDGES anomaly.

79 ASTRONOMY AND ASTROPHYSICS↗