Data-driven reduced-order models via regularised Operator Inference for a single-injector combustion process
Not Available
97 MATHEMATICS AND COMPUTING↗
SEARCH · Search NASA
Results for “regularisation”
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.
Quark and gluon two-loop beam functions for leading-jet p T and slicing at NNLO
We compute the complete set of two-loop beam functions for the transverse momentum distribution of the leading jet produced in association with an arbitrary colour-singlet system. Our results constitute the last missing ingredient for the calculation of the jet-vetoed cross section at small veto scales at the next-to-next-to-leading order, as well as an important ingredient for its resummation to next-to-next-to-next-to-leading logarithmic order. Our calculation is performed in the soft-collinear effective theory framework with a suitable regularisation of the rapidity divergences occurring in the phase-space integrals. We discuss the occurrence of soft-collinear mixing terms that might violate the factorisation theorem, and demonstrate that they are naturally absorbed into the beam functions at two loops in the exponential rapidity regularisation scheme when performing a multipole expansion of the measurement function. As in our recent computation of the two-loop soft function, we present the results as a Laurent expansion in the jet radius R. We provide analytic expressions for all flavour channels in x space with the exception of a set of R-independent non-logarithmic terms that are given as numerical grids. We also perform a fully numerical calculation with exact R dependence, and find that it agrees with our analytic expansion at the permyriad level or better. Our calculation allows us to define a next-to-next-to-leading order slicing method using the leading-jet p T as a slicing variable. As a check of our results, we carry out a calculation of the Higgs and Z boson total production cross sections at the next-to-next-to-leading order in QCD.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
The collisional particle-in-cell method for the Vlasov–Maxwell–Landau equations
We introduce an extension of the particle-in-cell method that captures the Landau collisional effects in the Vlasov–Maxwell–Landau equations. The method arises from a regularisation of the variational formulation of the Landau equation, leading to a discretisation of the collision operator that conserves mass, charge, momentum and energy, while increasing the (regularised) entropy. The collisional effects appear as a fully deterministic effective force, thus the method does not require any transport–collision splitting. The scheme can be used in arbitrary dimension, and for a general interaction, including the Coulomb case. We validate the scheme on scenarios such as the Landau damping, the two-stream instability and the Weibel instability, demonstrating its effectiveness in the numerical simulation of plasma.
Bailo, Rafael (ORCID:0000000180183799)↗
Revealing the structure of light pseudoscalar mesons at the electron–ion collider
The questions of how the bulk of the Universe's visible mass emerges and how it is manifest in the existence and properties of hadrons are profound, and probe the heart of strongly interacting matter. Paradoxically, the lightest pseudoscalar mesons appear to be key to a further understanding of the emergent mass and structure mechanisms. These mesons, namely, the pion and kaon, are the Nambu-Goldstone boson modes of quantum chromodynamics (QCD). Unravelling their partonic structure and the interplay between emergent and Higgs-boson mass mechanisms is a common goal of three interdependent approaches - continuum QCD phenomenology, lattice-regularised QCD, and the global analysis of parton distributions - linked to experimental measurements of hadron structure. Experimentally, the anticipated electron-ion collider will enable a revolution in our ability to study pion and kaon structures, accessed by scattering from the 'meson cloud' of the proton through the Sullivan process. With the goal of enabling a suite of measurements that can address these questions, we examine key reactions that identify the critical detector-system requirements needed to map tagged pion and kaon cross-sections over a wide range of kinematics. The excellent prospects for extracting pion structural, functional, and form-factor data are outlined, and similar prospects for kaon structures are discussed in the context of a worldwide programme. Finally, the successful completion of the programme outlined herein will deliver deep, far-reaching insights into the emergence of pions and kaons, their properties, and their role as QCD's Goldstone boson modes.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Insights into the emergence of mass from studies of pion and kaon structure
There are two mass generating mechanisms in the standard model of particle physics (SM). One is related to the Higgs boson and fairly well understood. The other is embedded in quantum chromodynamics (QCD), the SM’s strong interaction piece; and although responsible for emergence of the roughly 1 GeV mass scale that characterises the proton and hence all observable matter, the source and impacts of this emergent hadronic mass (EHM) remain puzzling. As bound states seeded by a valence-quark and -antiquark, pseudoscalar mesons present a simpler problem in quantum field theory than that associated with the nucleon. Consequently, there is a large array of robust predictions for pion and kaon properties whose empirical validation will provide a clear window onto many effects of both mass generating mechanisms and the constructive interference between them. This has now become significant because new-era experimental facilities, in operation, construction, or planning, are capable of conducting such tests and thereby contributing greatly to resolving the puzzles of EHM. These aspects of experiment, phenomenology, and theory, along with contemporary successes and challenges, are reviewed herein. In addition to providing an overview of the experimental status, we focus on recent progress made using continuum Schwinger function methods and lattice-regularised QCD. Further, advances made using other theoretical tools are also sketched. Our primary goal is to highlight the potential gains that can accrue from a coherent effort aimed at finally reaching an understanding of the character and structure of Nature’s Nambu–Goldstone modes.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Dirac traces and the Tutte polynomial
Perturbative calculations involving fermion loops in quantum field theories require tracing over Dirac matrices. A simple way to regulate the divergences that generically appear in these calculations is dimensional regularisation, which has the consequence of replacing 4-dimensional Dirac matrices with d-dimensional counterparts for arbitrary complex values of d. In this work, a connection between traces of d-dimensional Dirac matrices and computations of the Tutte polynomial of associated graphs is proven. The time complexity of computing Dirac traces is analysed by this connection, and improvements to algorithms for computing Dirac traces are proposed.
Renormalization and Regularization↗
The analytic two-loop soft function for leading-jet p T
We present the calculation of the two-loop soft function for the transverse momentum distribution of the leading jet produced in association with any colour-singlet system (e.g. a Higgs or a Z boson). This constitutes a central ingredient for the resummation of the above distribution as well as the jet-vetoed cross section at the next-to-next-to-next-to-leading logarithmic order, both of which play an important role in the precision physics programme at the Large Hadron Collider. The calculation is performed in soft-collinear effective theory with an appropriate regularisation of the rapidity divergences that occur in the phase-space integrals. We obtain analytic results by employing an exponential regulator and by taking a Laurent expansion in the jet radius R . All expressions are presented as supplementary material attached to this article.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach
This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.
direct numerical simulation↗
Diquark correlations in hadron physics: Origin, impact and evidence
The last decade has seen a marked shift in how the internal structure of hadrons is understood. Modern experimental facilities, new theoretical techniques for the continuum bound-state problem and progress with lattice-regularised QCD have provided strong indications that soft quark+quark (diquark) correlations play a crucial role in hadron physics. For example, theory indicates that the appearance of such correlations is a necessary consequence of dynamical chiral symmetry breaking, viz. a corollary of emergent hadronic mass that is responsible for almost all visible mass in the universe; experiment has uncovered signals for such correlations in the flavour-separation of the proton's electromagnetic form factors; and phenomenology suggests that diquark correlations might be critical to the formation of exotic tetra- and penta-quark hadrons. Therefore, a broad spectrum of such information is evaluated herein, with a view to consolidating the facts and therefrom moving toward a coherent, unified picture of hadron structure and the role that diquark correlations might play.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions
The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.
17 WIND ENERGY↗
The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from the anisotropic power spectrum of the quasar sample between redshift 0.8 and 2.2
ABSTRACT We measure the clustering of quasars of the final data release (DR16) of eBOSS. The sample contains $343\, 708$ quasars between redshifts 0.8 ≤ z ≤ 2.2 over $4699\, \mathrm{deg}^2$. We calculate the Legendre multipoles (0,2,4) of the anisotropic power spectrum and perform a BAO and a Full-Shape (FS) analysis at the effective redshift zeff = 1.480. The errors include systematic errors that amount to 1/3 of the statistical error. The systematic errors comprise a modelling part studied using a blind N-body mock challenge and observational effects studied with approximate mocks to account for various types of redshift smearing and fibre collisions. For the BAO analysis, we measure the transverse comoving distance DM(zeff)/rdrag = 30.60 ± 0.90 and the Hubble distance DH(zeff)/rdrag = 13.34 ± 0.60. This agrees with the configuration space analysis, and the consensus yields: DM(zeff)/rdrag = 30.69 ± 0.80 and DH(zeff)/rdrag = 13.26 ± 0.55. In the FS analysis, we fit the power spectrum using a model based on Regularised Perturbation Theory, which includes redshift space distortions and the Alcock–Paczynski effect. The results are DM(zeff)/rdrag = 30.68 ± 0.90 and DH(zeff)/rdrag = 13.52 ± 0.51 and we constrain the linear growth rate of structure f(zeff)σ8(zeff) = 0.476 ± 0.047. Our results agree with the configuration space analysis. The consensus analysis of the eBOSS quasar sample yields: DM(zeff)/rdrag = 30.21 ± 0.79, DH(zeff)/rdrag = 3.23 ± 0.47, and f(zeff)σ8(zeff) = 0.462 ± 0.045 and is consistent with a flat ΛCDM cosmological model using Planck results.
Neveux, Richard↗
Optimal Management of Grid-Interactive Efficient Buildings via Safe Reinforcement Learning
Reinforcement learning (RL)-based methods have achieved significant success in managing grid-interactive efficient buildings (GEBs). However, RL does not carry intrinsic guarantees of constraint satisfaction, which may lead to severe safety consequences. Besides, in GEB control applications, most existing safe RL approaches rely only on the regularisation parameters in neural networks or penalty of rewards, which often encounter challenges with parameter tuning and lead to catastrophic constraint violations. To provide enforced safety guarantees in controlling GEBs, this paper designs a physics-inspired safe RL method whose decision-making is enhanced through safe interaction with the environment. Different energy resources in GEBs are optimally managed to minimize energy costs and maximize customer comfort. The proposed approach can achieve strict constraint guarantees based on prior knowledge of a set of developed hard steady-state rules. Simulations on the optimal management of GEBs, including heating, ventilation, and air conditioning (HVAC), solar photovoltaics, and energy storage systems, demonstrate the effectiveness of the proposed approach.
Huo, Xiang↗
Shock structure for the seven-equation, two-phase continuum-mixture model
Not Available
5-equation reduced model↗