Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak scaling”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First Detection of the GI-type of Intrinsic Alignments of Galaxies Using the Self-calibration Method in a Photometric Galaxy Survey

Weak gravitational lensing is one of the most promising cosmological probes to constrain dark matter, dark energy, and the nature of gravity at cosmic scales. Intrinsic alignments (IAs) of galaxies have been recognized as one of the most serious systematic effects facing gravitational lensing. Such alignments must be isolated and removed to obtain a pure lensing signal. Furthermore, the alignments are related to the processes of galaxy formation, so their extracted signal can help in understanding such formation processes and improving their theoretical modeling. We report in this Letter the first detection of the gravitational shear–intrinsic shape (GI) correlation and the intrinsic shape–galaxy density (Ig) correlation using the self-calibration method in a photometric redshift survey. These direct measurements are made from the KiDS-450 photometric galaxy survey with a significance of 3.65σ in the third bin for the Ig correlation, and 3.51σ for the GI cross-correlation between the third and fourth bins. The self-calibration method uses the information available from photometric surveys without needing to specify an IA model and will play an important role in validating IA models and IA mitigation in future surveys such as the Rubin Observatory Legacy Survey of Space and Time, Euclid, and WFIRST.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Interpreting Sunyaev–Zel’dovich observations with MillenniumTNG: mass and environment scaling relations

ABSTRACT Sunyaev–Zel’dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton-y parameter, Y, receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y–M and τ–M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (${\sim }10^{13} \, {\rm M}_{\odot }$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L* galaxies.

Hadzhiyska, Boryana (ORCID:0000000223123121)↗
DOE OSTI2020

Non-Gaussianity in the weak lensing correlation function likelihood – implications for cosmological parameter biases

ABSTRACT We study the significance of non-Gaussianity in the likelihood of weak lensing shear two-point correlation functions, detecting significantly non-zero skewness and kurtosis in 1D marginal distributions of shear two-point correlation functions in simulated weak lensing data. We examine the implications in the context of future surveys, in particular LSST, with derivations of how the non-Gaussianity scales with survey area. We show that there is no significant bias in 1D posteriors of Ωm and σ8 due to the non-Gaussian likelihood distributions of shear correlations functions using the mock data (100 deg2). We also present a systematic approach to constructing approximate multivariate likelihoods with 1D parametric functions by assuming independence or more flexible non-parametric multivariate methods after decorrelating the data points using principal component analysis (PCA). While the use of PCA does not modify the non-Gaussianity of the multivariate likelihood, we find empirically that the 1D marginal sampling distributions of the PCA components exhibit less skewness and kurtosis than the original shear correlation functions. Modelling the likelihood with marginal parametric functions based on the assumption of independence between PCA components thus gives a lower limit for the biases. We further demonstrate that the difference in cosmological parameter constraints between the multivariate Gaussian likelihood model and more complex non-Gaussian likelihood models would be even smaller for an LSST-like survey. In addition, the PCA approach automatically serves as a data compression method, enabling the retention of the majority of the cosmological information while reducing the dimensionality of the data vector by a factor of ∼5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Interacting dark sector within ETHOS: Cosmological constraints from SPT cluster abundance with DES and HST weak lensing data

We use galaxy cluster abundance measurements from the South Pole Telescope enhanced by multicomponent matched filter confirmation and complemented with mass information obtained using weak-lensing data from Dark Energy Survey Year 3 (DES Y3) and targeted Hubble Space Telescope observations for probing deviations from the cold dark matter paradigm. Concretely, we consider a class of dark sector models featuring interactions between dark matter (DM) and a dark radiation (DR) component within the framework of the effective theory of structure formation (ETHOS). We focus on scenarios that lead to power suppression over a wide range of scales, and thus can be tested with data sensitive to large scales, as realized, for example, for DM–DR interactions following from an unbroken non-Abelian 𝑆⁢𝑈⁡(𝑁) gauge theory (interaction rate with power-law index 𝑛 = 0 within the ETHOS parametrization). Cluster abundance measurements are mostly sensitive to the amount of DR interacting with DM, parametrized by the ratio of DR temperature to the cosmic microwave background (CMB) temperature, 𝜉 DR = 𝑇 DR /𝑇 CMB . We find an upper limit 𝜉 DR < 17% at 95% credibility. When the cluster data are combined with Planck 2018 CMB data along with baryon acoustic oscillation (BAO) measurements we find 𝜉 DR < 10%, corresponding to a limit on the abundance of interacting DR that is around 3 times tighter than that from CMB + BAO data alone. We also discuss the complementarity of weak lensing informed cluster abundance studies with probes sensitive to smaller scales, explore the impact on our analysis of massive neutrinos, and comment on a slight preference for the presence of a nonzero interacting DR abundance, which enables a physical solution to the 𝑆 8 tension.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Noncoordinating Secondary Sphere Ion Modulates Supramolecular Clustering of Lanthanides

The role of counterions in molecular recognition of lanthanides is underexplored, especially when they exhibit weak interactions with the metal cations. Here, we report a complementary and comprehensive investigation integrating theoretical calculations with X-ray absorption fine structure spectroscopy, dynamic light scattering, and small-angle X-ray scattering to reveal atomic-scale structural features beyond the immediate coordination sphere of a system used for rare-earth element separations. Here, our results indicate the formation of an unusual T-shaped outer-sphere lanthanide complex, containing two ligands and two nitrate ions in the first coordination sphere, whereas the third nitrate is weakly coordinated and resides in the second shell. This unique structural arrangement causes inhomogeneous charge distribution, leading to self-assembly of the complexes into larger nanoclusters through sterically directed electrostatic interactions in the nonpolar medium. Our findings point to the importance of "noncoordinating" anions in defining the degree of supramolecular aggregation and ion cluster assembly.

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

Systematic Study of Projection Biases in Weak Lensing Analysis

The nature of Dark Energy is possibly one of the most fundamental questions in Physics today. Comprising 70% of the energy density of the Universe, and being responsible for the accelerated expansion of the Universe, this mysterious form of energy with negative pressure is challenging our current understanding of the fundamental laws of nature. Cosmological parameters like total relative matter Density Ωm, the normalization of power spectrum σ8 and their combination S8 = σ8 p Ωm/0.3 are used to quantify our understanding of the evolution of the large-scale structure in the Universe and shed light on the accelerated expansion, caused by the Dark Energy, that started about halfway through the life of the Universe. Sky surveys, like Dark Energy Survey (DES), utilize weak lensing analysis to use the observations of galaxies and their shapes and extract these cosmological parameters. By doing an exhaustive study and quantifying the projection biases in the weak lensing statistics, the research work presented in this dissertation aims to contribute to the precise measurements of the cosmological parameters Ωm, σ8, and S8. For the analysis, we used the standard cosmological model that includes dark energy as the cosmological constant. And, we used the data taken by the DES experiment, which is one of the most complete, deep, and wide field surveys of the sky today.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Dynamic Pressure Drag on Rising Buoyant Thermals in a Neutrally Stable Environment

Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.

Meteorology & Atmospheric Sciences↗
DOE OSTI2024

Exciton–Phonon Coupling Induces a New Pathway for Ultrafast Intralayer-to-Interlayer Exciton Transition and Interlayer Charge Transfer in WS 2 –MoS 2 Heterostructure: A First-Principles Study

Despite the weak, van der Waals interlayer coupling, photoinduced charge transfer vertically across atomically thin interfaces can occur within surprisingly fast, sub-50 fs time scales. An early theoretical understanding of charge transfer is based on a noninteracting picture, neglecting excitonic effects that dominate optical properties of such materials. We employ an ab initio many-body perturbation theory approach, which explicitly accounts for the excitons and phonons in the heterostructure. Our large-scale first-principles calculations directly probe the role of exciton-phonon coupling in the charge dynamics of the WS 2 /MoS 2 heterobilayer. We find that the exciton-phonon interaction induced relaxation time of photoexcited excitons at the K valley of MoS 2 and WS 2 is 67 and 15 fs at 300 K, respectively, which sets a lower bound to the intralayer-to-interlayer exciton transfer time and is consistent with experiment reports. We further show that electron-hole correlations facilitate novel transfer pathways that are otherwise inaccessible to noninteracting electrons and holes.

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

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Vertical detachment energies of ammonia cluster anions using self-interaction-corrected methods

Systems with weakly bound extra electrons impose great challenges to semilocal density functional approximations (DFAs), which suffer from self-interaction errors. Small ammonia clusters are one such example of weakly bound anions where the extra electron is weakly bound. We applied two self-interaction correction (SIC) schemes, viz., the well-known Perdew–Zunger and the recently developed locally scaled SIC (LSIC) with the local spin density approximation (LSDA), Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation (GGA), and the SCAN meta-GGA functionals to calculate the vertical detachment energies (VDEs) of small ammonia cluster anions (NH3)n–. Our results show that the LSIC significantly reduces the errors in calculations of VDE with LSDA and PBE-GGA functionals leading to better agreement with the reference values calculated with coupled cluster singles and doubles with perturbative triples [CCSD(T)]. Accurate prediction of VDE as an absolute of the highest occupied molecular orbital (HOMO) is challenging for DFAs. Our results show that VDEs estimated from the negative of HOMO eigenvalues with the LSIC-LSDA and Perdew–Zunger SIC-PBE are within 11 meV of the reference CCSD(T) results. Here, the LSIC method performs consistently well for the VDE estimates, from both the total energy differences and the absolute HOMO eigenvalues.

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

Constraints on compact objects from the Dark Energy Survey 5-yr supernova sample

ABSTRACT Gravitational lensing magnification of Type Ia supernovae (SNe Ia) allows information to be obtained about the distribution of matter on small scales. In this paper, we derive limits on the fraction $\alpha$ of the total matter density in compact objects (which comprise stars, stellar remnants, small stellar groupings, and primordial black holes) of mass M > 0.03 ${\rm M}_{\odot }$ over cosmological distances. Using 1532 SNe Ia from the Dark Energy Survey Year 5 sample (DES-SN5YR) combined with a Bayesian prior for the absolute magnitude M, we obtain α < 0.12 at the 95 per cent confidence level after marginalization over cosmological parameters, lensing due to large-scale structure, and intrinsic non-Gaussianity. Similar results are obtained using priors from the cosmic microwave background, baryon acoustic oscillations, and galaxy weak lensing, indicating our results do not depend on the background cosmology. We argue our constraints are likely to be conservative (in the sense of the values we quote being higher than the truth), but discuss scenarios in which they could be weakened by systematics of the order of $\Delta \alpha \sim 0.04$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Two-neutrino double-𝛽 decay in pionless effective field theory from a Euclidean finite-volume correlation function

Two-neutrino double-β decay of certain nuclear isotopes is one of the rarest Standard Model processes observed in nature. Its neutrinoless counterpart is an exotic lepton-number nonconserving process that is widely searched to determine if the neutrinos are Majorana fermions. In order to connect the rate of these processes to the Standard Model and beyond the Standard Model interactions, it is essential that the corresponding nuclear matrix elements are constrained reliably from theory. Lattice quantum chromodynamics (LQCD) and low-energy effective field theories (EFTs) are expected to play an essential role in constraining the matrix element of the two-nucleon subprocess, which could in turn provide the input into ab initio nuclear-structure calculations in larger isotopes. Focusing on the two-neutrino process $nn\rightarrow pp(ee\bar{v}_{e}\bar{v}_e)$, the amplitude is constructed in this work in pionless EFT at next-to-leading order, demonstrating the emergence of a renormalization-scale independent amplitude and the absence of any new low-energy constant at this order beyond those present in the single-weak process. Most importantly, it is shown how a LQCD four-point correlation function in Euclidean and finite-volume spacetime can be used to constrain the Minkowski infinite-volume amplitude in the EFT. The same formalism is provided for the related single-weak process, which is an input to the double-β decay formalism. The LQCD-EFT matching procedure outlined for the double-weak amplitude paves the road toward constraining the two-nucleon matrix element entering the neutrinoless double-β decay amplitude with a light Majorana neutrino.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Evaluating the effects of structural disorder on the magnetic properties of Nd 2 Zr 2 O 7

Motivated by the variation in reported lattice parameters of floating-zone-grown Nd 2 Zr 2 O 7 crystals, we have performed a detailed study of the relationship between synthesis environment, structural disorder, and magnetic properties. Using a combination of polycrystalline standards, electron-probe microanalysis, and scattering techniques, we show that crystals grown under atmospheric conditions have a reduced lattice parameter relative to pristine polycrystalline powders due to occupation of the Nd site by excess Zr (i.e., negative stuffing). In contrast, crystals grown under high-pressure Ar are nearly stoichiometric with an average lattice parameter approaching the polycrystalline value. While minimal disorder of the oxygen sublattices is observed on the scale of the average structure, neutron pair-distribution function analysis indicates a highly local disorder of the oxygen coordination, which is only weakly dependent on growth environment. Most importantly, our magnetization, heat capacity, and single-crystal neutron scattering data show that the magnetic properties of crystals grown under high-pressure Ar match closely with those of stoichiometric powders. Neutron scattering measurements reveal that the signature of magnetic moment fragmentation—the coexistence of all-in-all-out (AIAO) magnetic Bragg peaks and diffuse pinch-point scattering due to spin-ice correlations–persists in these nearly stoichiometric crystals. However, in addition to an increased AIAO transition temperature, the diffuse signal is seemingly stabilized and remains nearly unchanged upon warming to 800 mK. This behavior indicates that both the AIAO magnetic order and spin-ice correlations are sensitive to deviations of the Nd stoichiometry.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Spatial–Temporal Patterns of Historical, Near-Term, and Projected Drought in the Conterminous United States

Major droughts in the United States have heavily impacted the hydrologic system, negatively effecting energy and food production. Improved understanding of historical drought is critical for accurate forecasts. Data from global climate models (GCMs), commonly used to assess drought, cannot effectively evaluate local patterns because of their low spatial scale. This research leverages downscaled (~4 km grid spacing) temperature and precipitation estimates from nine GCMs’ data under the business-as-usual scenario (Representative Concentration Pathway 8.5) to examine drought patterns. Drought severity is estimated using the Palmer Drought Severity Index (PDSI) with the Thornthwaite evapotranspiration method. The specific objectives were (1) To reproduce historical (1966–2005) drought and calculate near-term to future (2011–2050) drought patterns over the conterminous USA. (2) To uncover the local variability of spatial drought patterns in California between 2012 and 2018 using a network-based approach. Our estimates of land proportions affected by drought agree with the known historical drought events of the mid-1960s, late 1970s to early 1980s, early 2000s, and between 2012 and 2015. Network analysis showed heterogeneity in spatial drought patterns in California, indicating local variability of drought occurrence. The high spatial scale at which the analysis was performed allowed us to uncover significant local differences in drought patterns. This is critical for highlighting possible weak systems that could inform adaptation strategies such as in the energy and agricultural sectors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Code for the manuscript "Lagrangian Large Eddy Simulations vis Physics-informed Machine Learning"

We would like to disclose a code based on Python/Pytorch, in which we implement the data-driven physics-informed machine learning of Lagrangian large eddy simulation that is described in the manuscript "Lagrangian Large Eddy Simulations via Physics-Informed Machine Learning ". High-Reynolds number Homogeneous Isotropic Turbulence is fully described within the Navier-Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as Large Eddy Simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about the effects of sub-grid scales on the resolved scales. Here, we take an alternative approach and design novel LES heuristics stated in terms of Lagrangian particles moving with the flow. Our \emph{Lagrangian LES}, thus L-LES, is described by equations generalizing the weakly compressible Smoothed Particle Hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from Direct Numerical Simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The sub-grid scale contributions are modeled separately with physical constraints to account for the effects from un-resolved scales. We build the resulting model under the Differentiable Programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our Lagrangian LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

Tian, Yifeng↗
DOE OSTI2020

CODEX weak lensing mass catalogue and implications on the mass–richness relation

The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at 0.35 < z < 0.65. It was selected from ROSAT data in the 10 000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX programme on cluster mass measurement, including a reanalysis of CFHTLS Wide data, with 25 individual lensing-constrained cluster masses. We employ lensfit shape measurement and perform a conservative colour–space selection and weighting of background galaxies. Using the combination of shape noise and an analytic covariance for intrinsic variations of cluster profiles at fixed mass due to large-scale structure, miscentring, and variations in concentration and ellipticity, we determine the likelihood of the observed shear signal as a function of true mass for each cluster. We combine 25 individual cluster mass likelihoods in a Bayesian hierarchical scheme with the inclusion of optical and X-ray selection functions to derive constraints on the slope α, normalization β, and scatter σln λ|μ of our richness–mass scaling relation model in log-space: ${\langle {\rm In}\,\, \lambda\!\!\mid\!\!\mu\rangle = \alpha\mu + \beta,} $ with μ = ln (M 200c /M piv ), and M piv = 10 14.81 M ⊙ . We find a slope $\alpha = 0.49^{+0.20}_{-0.15}$, normalization $\exp (\beta) = 84.0^{+9.2}_{-14.8}$, and $\sigma _{\ln \lambda | \mu } = 0.17^{+0.13}_{-0.09}$ using CFHT richness estimates. In comparison to other weak lensing richness–mass relations, we find the normalization of the richness statistically agreeing with the normalization of other scaling relations from a broad redshift range (0.0 < z < 0.65) and with different cluster selection (X-ray, Sunyaev–Zeldovich, and optical).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗