Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SCALARS”

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

Search for physics beyond the standard model in multilepton final states in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A search for physics beyond the standard model in events with at least three charged leptons (electrons or muons) is presented. The data sample corresponds to an integrated luminosity of 137 fb$^{−1}$ of proton-proton collisions at $ \sqrt{s} $ = 13 TeV, collected with the CMS detector at the LHC in 2016–2018. The two targeted signal processes are pair production of type-III seesaw heavy fermions and production of a light scalar or pseudoscalar boson in association with a pair of top quarks. The heavy fermions may be manifested as an excess of events with large values of leptonic transverse momenta or missing transverse momentum. The light scalars or pseudoscalars may create a localized excess in the dilepton mass spectra. The results exclude heavy fermions of the type-III seesaw model for masses below 880 GeV at 95% confidence level in the scenario of equal branching fractions to each lepton flavor. This is the most restrictive limit on the flavor-democratic scenario of the type-III seesaw model to date. Assuming a Yukawa coupling of unit strength to top quarks, branching fractions of new scalar (pseudoscalar) bosons to dielectrons or dimuons above 0.004 (0.03) and 0.04 (0.03) are excluded at 95% confidence level for masses in the range 15–75 and 108–340 GeV, respectively. These are the first limits in these channels on an extension of the standard model with scalar or pseudoscalar particles.[graphic not available: see fulltext]

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Impact of new experimental data on the C2HDM: the strong interdependence between LHC Higgs data and the electron EDM

The complex two-Higgs doublet model (C2HDM) is one of the simplest extensions of the Standard Model with a source of CP-violation in the scalar sector. It has a $\mathbb{Z}$ 2 symmetry, softly broken by a complex coefficient. There are four ways to implement this symmetry in the fermion sector, leading to models known as Type-I, Type-II, Lepton Specific and Flipped. In the latter three models, there is a priori the surprising possibility that the 125 GeV Higgs boson couples mostly as a scalar to top quarks, while it couples mostly as a pseudoscalar to bottom quarks. This “maximal” scenario was still possible with the data available in 2017. Since then, there have been more data on the 125 GeV Higgs boson, direct searches for CP-violation in angular correlations of $τ$-leptons produced in Higgs boson decays, new results on the electron electric dipole moment, new constraints from LHC searches for additional Higgs bosons and new results on $b$ → $sγ$ transitions. Highlighting the crucial importance of the physics results of LHC’s Run 2, we combine all these experiments and show that the “maximal” scenario is now excluded in all models. Still, one can have a pseudoscalar component in $hτ\overline{τ}$ couplings in the Lepton-Specific case as large as 87% of the scalar component for all mass orderings of the neutral scalar bosons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Search for $t\bar{t}H/A$ → $t\bar{t}t\bar{t}$ production in the multilepton final state in proton–proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

A search for a new heavy scalar or pseudo-scalar Higgs boson ($H/A$) produced in association with a pair of top quarks, with the Higgs boson decaying into a pair of top quarks ($H/A$ → $t\bar{t}$) is reported. The search targets a final state with exactly two leptons with same-sign electric charges or at least three leptons. The analysed dataset corresponds to an integrated luminosity of 139 fb –1 of proton–proton collisions collected at a centre-of-mass energy of 13 TeV with the ATLAS detector at the LHC. Two multivariate classifiers are used to separate the signal from the background. No significant excess of events over the Standard Model expectation is observed. The results are interpreted in the context of a type-II two-Higgs-doublet model. The observed (expected) upper limits at 95% confidence level on the $t\bar{t}$$H/A$ production cross-section times the branching ratio of $H/A$ → $t\bar{t}$ range between 14 (10) fb and 6 (5) fb for a heavy Higgs boson with mass between 400 GeV and 1000 GeV, respectively. Assuming that only one particle, either the scalar $H$ or the pseudo-scalar $A$, contributes to the $t\bar{t}t\bar{t}$ final state, values of tan $β$ below 1.2 or 0.5 are excluded for a mass of 400 GeV or 1000 GeV, respectively. These exclusion ranges increase to tan $β$ below 1.6 or 0.6 when both particles are considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Higgs inflation, vacuum stability, and leptogenesis

We consider the introduction of a complex scalar field carrying a global lepton number charge to the Standard Model and the Higgs inflation framework. The conditions are investigated under which this model can simultaneously ensure Higgs vacuum stability up to the Planck scale, successful inflation, non-thermal Leptogenesis via the pendulum mechanism, and light neutrino masses. These can be simultaneously achieved when the scalar lepton is minimally coupled to gravity, that is, when standard Higgs inflation and reheating proceed without the interference of the additional scalar degrees of freedom. If the scalar lepton also has a non-minimal coupling to gravity, a multi-field inflation scenario is induced, with interesting interplay between the successful inflation constraints and those from vacuum stability and Leptogenesis. The parameter region that can simultaneously achieve the above goals is explored.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Effective Field Theory for the perturbations of a slowly rotating black hole

We develop the effective theory for perturbations around black holes with scalar hair, in two directions. First, we show that the scalar-Gauss-Bonnet theory, often used as an example exhibiting scalar black hole hair, can be deformed by galileon operators leading to order unity changes to its predictions. The effective theory for perturbations thus provides an efficient framework for describing and constraining broad classes of scalar-tensor theories, of which the addition of galileon operators is an example. Second, we extend the effective theory to perturbations around an axisymmetric, slowly rotating black hole, at linear order in the black hole spin. We also discuss the inclusion of parity-breaking operators in the effective theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

tapir: A tool for topologies, amplitudes, partial fraction decomposition and input for reductions

The demand for precision predictions in the field of high energy physics has dramatically increased over recent years. Experiments conducted at the LHC, as well as precision measurements at the intensity frontier such as Belle II require equally precise theoretical predictions to make full use of the acquired data. To match the experimental precision, second-, third- and, for certain quantities, even higher-order calculations in perturbative quantum field theory are required. To facilitate such calculations, computer software automating as many steps as possible is required. Yet, each calculation poses different challenges and thus, a high level of configurability is required. In this context we present tapir: a tool for identification, manipulation and minimization of Feynman integral families. It is designed to integrate in toolchains based on the computer algebra system FORM, the use of which is common practice in the field. tapir can be used to reduce the complexity of multi-loop problems with cut-filters, topology mapping, partial fraction decomposition and alike. Program Title:tapir CPC Library link to program files:https://doi.org/10.17632/ptc9t46xyn.1 Developer's repository link:https://gitlab.com/tapir-devs/tapir Licensing provisions: GPLv3 Programming language:python 3, C++ Nature of problem: Multi-loop computations require the automatization of a large number of different tasks related to Feynman integral topologies. Among them are the identification and minimization of integral topologies, partial fraction decomposition of topologies in the case of linearly dependent propagators as well as mapping scalar products of loop momenta to scalar functions. Solution method: The minimization of topologies is performed by comparison of their respective Nickel indices [1], even further minimization utilizes Pak's algorithm [2]. To efficiently map scalar products of loop momenta to scalar functions FORM [3] code is generated. Additional comments including restrictions and unusual features: Minimization based on Pak's algorithm slows down for many lines and scales. A coarser minimization using the Nickel indices, however, is still possible. [1]B. Nickel, D. Meiron, G.A.J. Baker, Compilation of 2-pt and 4-pt graphs for continuous spin model, Report, University of Guelph, 1977.[2]A. Pak, J. Phys. Conf. Ser. 368 (2012) 012049, https://doi.org/10.1088/1742-6596/368/1/012049, arXiv:1111.0868.[3]B. Ruijl, T. Ueda, J. Vermaseren, FORM version 4.2, arXiv:1707.06453, 7 2017.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ordinary and exotic mesons in the extended Linear Sigma Model

The extended Linear Sigma Model (eLSM) is a hadronic model based on the global symmetries of QCD and the corresponding explicit, anomalous, and spontaneous breaking patterns. In its basic three-flavor form, its mesonic part contains the dilaton/glueball as well as the nonets of pseudoscalar, scalar, vector, and axial–vector mesons, thus chiral symmetry is linearly realized. In the chiral limit and neglecting the chiral anomaly, only one term – within the dilaton potential – breaks dilatation invariance, and all terms are chirally symmetric. Spontaneous symmetry breaking is implemented by a generalization of the Mexican-hat potential, with explicit symmetry breaking responsible for its tilting. The overall mesonic phenomenology up to ~2 GeV is in agreement with the PDG compilation of masses and partial and total decay widths. The eLSM was enlarged in a straightforward way to include other conventional quark–antiquark nonets (pseudovector and orbitally excited vector mesons, tensor and axial-tensor mesons, radially excited (pseudo)scalar mesons, etc.), as well as two nonets of hybrid mesons, the lightest one with exotic quantum numbers J $\mathscr{PC}$ = 1 -+ not allowed for $\bar{q}q$ objects, such as the resonance π 1 (1600) and the recently discovered η 1 (1855) . In doing so, different types of chiral multiplets are introduced: heterochiral and homochiral multiplets, which differ in the way they transform under chiral transformations. Moreover, besides the scalar glueball that is present from the beginning as dilaton, other glueballs, the tensor, the pseudoscalar and the vector glueballs were coupled to the eLSM: the scalar resonance ƒ 0 (1710) turns out to be mostly gluonic, the tensor glueball couples strongly to vector mesons, and the pseudoscalar glueball couples sizably to $ππη'$ and can be assigned to X (2370) or X (2600). In all cases above, masses and decays can be analyzed allowing for a better understanding of both conventional and non-conventional mesons: whenever data are available, a comparison is performed and, when this is not the case, predictions of decay widths and decay ratios are outlined. The eLSM contains chiral partners on an equal footing and is therefore well suited for studies of chiral symmetry restoration at nonzero temperature and densities: this is done by coupling it to the Polyakov loop. Finally, the QCD phase diagram and the location of the critical endpoint were investigated within this framework.

Chiral symmetry↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2024

Present and future of $\mathcal{C}$osmo $\mathcal{L}$attice

We discuss the present state and planned updates of Cosmo Lattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review the current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of Cosmo Lattice that we plan to release publicly. We comment on new physics modules, which include axion-gauge interactions $\phi F\tilde{F}$, non-minimal gravitational couplings $\phi ^2R$, creation and evolution of cosmic-defect networks, and magnetohydrodynamics. Finally, we also discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher-accuracy spatial derivatives.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Minimally entangled state preparation of localized wave functions on quantum computers

Initializing a single site of a lattice scalar field theory into an arbitrary state with support throughout the quantum register requires $O(2^{nQ})$ entangling gates on a quantum computer with $n_Q$ qubits per site. It is conceivable that instead initializing to functions that are good approximations to states may have utility in reducing the number of required entangling gates. In the case of a single site of a noninteracting scalar field theory, initializing to a symmetric exponential wave function requires $n_Q – 1$ entangling gates, the minimal number necessary to create an entangled state of all $n_Q$ qubits. This is compared with the $2^{nQ – 1} + n_Q – 3 + δ_{nQ, 1}$ required for a symmetric Gaussian wave function. In this work, we explore the initialization of one-site $(n_Q = 4)$, two-site $(n_Q = 3)$, and three-site ( n Q = 3 ) noninteracting scalar field theories with symmetric exponential wave functions using IBM's quantum simulators and quantum devices (Poughkeepsie and Tokyo). With the digitizations attainable with $n_Q = 3,4$, these tensor-product wave functions are found to have large overlap with a Gaussian wave function and provide a suitable low-noise initialization for improvement and Somma inflation. Further, in performing these simulations, we have employed a workflow that interleaves calibrations to mitigate systematic errors in production. The calibrations allow tolerance cuts on gate performance including the fidelity of the symmetrizing Hadamard gate, both in vacuum $(|\textbf{0}\gg^{\bigotimes nQ})$ and in medium ($n_Q –1$ qubits initialized to an exponential function). The results obtained in this work are relevant to systems beyond scalar field theories, such as the deuteron radial wave function, two- and three-dimensional Cartesian-space wave functions, and nonrelativistic multinucleon systems built on a localized eigenbasis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Microphysics of early dark energy

Early dark energy (EDE) relies on scalar field dynamics to resolve the Hubble tension, by boosting the pre-recombination length scales and thereby raising the CMB-inferred value of the Hubble constant into agreement with late universe probes. However, the collateral effect of scalar field microphysics on the linear perturbation spectra appears to preclude a fully satisfactory solution. H 0 is not raised without the inclusion of a late universe prior, and the “S 8 -tension”, a discrepancy between early- and late-universe measurements of the structure growth parameter, is exacerbated. What if EDE is not a scalar field? Here, we investigate whether different microphysics, encoded in the constitutive relationships between pressure and energy density fluctuations, can relieve these tensions. We show that EDE with an anisotropic sound speed can soften both the H 0 and S 8 tensions while still providing a quality fit to CMB data. Future observations from the CMB-S4 experiment may be able to distinguish the underlying microphysics at the 4σ level, and thereby test whether a scalar field or some richer physics is at work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

How do axisymmetric black holes grow monopole and dipole hair?

Here, we study the dynamical formation of scalar monopole and dipole hair in scalar Gauss-Bonnet theory and dynamical Chern-Simons theory. We prove that the spherically symmetric mode of the dipole hair is completely determined by the product of the mass of the spacetime and the value of the monopole hair. We then show that the dynamics of the ℓ = 1 mode of the dipole hair is intimately tied to the appearance of the event horizon during axisymmetric collapse, which results in the radiation of certain modes that could have been divergent in the future of the collapse. We confirm these analytical predictions by simulating the gravitational collapse of a rapidly rotating neutron star in the decoupling limit, both in scalar Gauss-Bonnet and dynamical Chern-Simons theory. Our results, combined with those provide a clear physical picture of the dynamics of scalar monopole and dipole radiation in axisymmetric and spherical gravitational collapse in these theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Hidden conformal symmetry from the lattice

We analyze newly expanded and refined data from lattice studies of an SU(3) gauge theory with eight Dirac fermions in the fundamental representation. We focus on the light composite states emerging from these studies, consisting of a set of pseudoscalars and a single light scalar. We first consider the view that this theory is just outside the conformal window. In this case, the pseudoscalars arise from spontaneous breaking of chiral symmetry. Identifying the scalar in this case as an approximate dilaton, we fit the lattice data to a dilaton effective field theory, finding that it yields a good fit even at lowest order. For comparison, we then consider the possibility that the theory is inside the conformal window. The fermion mass provides a deformation, triggering confinement. We employ simple scaling laws to fit the lattice data, and find that it is of lesser quality.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Search for a μ + μ − resonance in four-muon final states at Belle II

We report on a search for a resonance X decaying to a pair of muons in e + e − → μ + μ − X events in the 0.212 – 9.000 GeV / c 2 mass range, using 178 fb − 1 of data collected by the Belle II experiment at the SuperKEKB collider at a center of mass energy of 10.58 GeV. The analysis probes two different models of X beyond the standard model: a Z ′ vector boson in the L μ − L τ model and a muonphilic scalar. We observe no evidence for a signal and set exclusion limits at the 90% confidence level on the products of cross section and branching fraction for these processes, ranging from 0.046 fb to 0.97 fb for the L μ − L τ model and from 0.055 fb to 1.3 fb for the muonphilic scalar model. For masses below 6 GeV / c 2 , the corresponding constraints on the couplings of these processes to the standard model range from 0.0008 to 0.039 for the L μ − L τ model and from 0.0018 to 0.040 for the muonphilic scalar model. These are the first constraints on the muonphilic scalar from a dedicated search. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Search for $t\bar{t}H/A \rightarrow t\bar{t}t\bar{t}$ production in proton–proton collisions at $\sqrt{s}=13$ $\text {TeV}$ with the ATLAS detector

A search is presented for a heavy scalar (H) or pseudo-scalar (A) predicted by the two-Higgs-doublet models, where the H/A is produced in association with a top-quark pair $(t\bar{t}H/A),$ and with the H/A decaying into a $t\bar{t}$ pair. The full LHC Run 2 proton–proton collision data collected by the ATLAS experiment is used, corresponding to an integrated luminosity of $139~\text {fb}^{-1}.$ Events are selected requiring exactly one or two opposite-charge electrons or muons. Data-driven corrections are applied to improve the modelling of the $t\bar{t}$ +jets background in the regime with high jet and b-jet multiplicities. These include a novel multi-dimensional kinematic reweighting based on a neural network trained using data and simulations. An H/A-mass parameterised graph neural network is trained to optimise the signal-to-background discrimination. In combination with the previous search performed by the ATLAS Collaboration in the multilepton final state, the observed upper limits on the $t\bar{t}H/A \rightarrow t\bar{t}t\bar{t}$ production cross-section at 95% confidence level range between 14 fb and 5.0 fb for an H/A with mass between 400 GeV and 1000 GeV , respectively. Assuming that both the H and A contribute to the $t\bar{t}t\bar{t}$ cross-section, tan β values below 1.7 or 0.7 are excluded for a mass of 400 GeV or 1000 GeV , respectively. The results are also used to constrain a model predicting the pair production of a colour-octet scalar, with the scalar decaying into a $t\bar{t}$ pair.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Probing neutrino mass models through resonances at neutrino telescopes

We study the detection prospects of relatively light charged scalars in radiative Majorana neutrino mass models, such as the Zee model and its variants using scalar leptoquarks, at current and future neutrino telescopes. In particular, we show that these scalar mediators can give rise to Glashow-like resonance features in the ultrahigh energy (UHE) neutrino event spectrum at the IceCube neutrino observatory and its high-energy upgrade IceCube-Gen2. The same scalars can also give rise to observable nonstandard neutrino interactions (NSIs), and we show that the UHE neutrinos provide a complementary probe of NSI. We also discuss an interesting possibility of producing such resonances by incoming sterile neutrino components in the case where neutrinos are pseudo-Dirac particles.

Physics↗
DOE OSTI2023

Lagrangian and Eulerian Supersaturation Statistics in Turbulent Cloudy Rayleigh–Bénard Convection: Applications for LES Subgrid Modeling

Abstract Turbulent fluctuations of scalar and velocity fields are critical for cloud microphysical processes, e.g., droplet activation and size distribution evolution, and can therefore influence cloud radiative forcing and precipitation formation. Lagrangian and Eulerian water vapor, temperature, and supersaturation statistics are investigated in direct numerical simulations (DNS) of turbulent Rayleigh–Bénard convection in the Pi Convection Cloud Chamber to provide a foundation for parameterizing subgrid-scale fluctuations in atmospheric models. A subgrid model for water vapor and temperature variances and covariance and supersaturation variance is proposed, valid for both clear and cloudy conditions. Evaluation of phase change contributions through an a priori test using DNS data shows good performance of the model. Supersaturation is a nonlinear function of temperature and water vapor, and relative external fluxes of water vapor and heat (e.g., during entrainment-mixing and phase change) influence turbulent supersaturation fluctuations. Although supersaturation has autocorrelation and structure functions similar to the independent scalars (temperature and water vapor), the autocorrelation time scale of supersaturation differs. Relative scalar fluxes in DNS without cloud make supersaturation PDFs less skewed than the adiabatic case, where they are highly negatively skewed. However, droplet condensation changes the PDF shape response: it becomes positively skewed for the adiabatic case and negatively skewed when the sidewall relative fluxes are large. Condensation also increases correlations between water vapor and temperature in the presence of relative scalar fluxes but decreases correlations for the adiabatic case. These changes in correlation suppress supersaturation variability for the nonadiabatic cases and increase it for the adiabatic case. Implications of this work for subgrid microphysics modeling using a Lagrangian stochastic scheme are also discussed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Concurrent two-way coupling of global and local models across internal boundaries with non-matching discretizations

Coupling local and global models enables efficient simulation of multiscale systems, where global models capture large-scale behavior and local models, with enhanced physics, resolve finer details over a smaller region. Here, this paper presents a mathematically consistent method for coupling physics-based models of varying fidelity across adjacent, non-overlapping subdomains, even when discretizations do not match at the immersed interdomain interfaces. Incompressible Navier-Stokes equations (NSE) constitute the global model while residual-based turbulence model serves as the local high-fidelity model. In addition, a scalar advection-diffusion equation that models the convection of an active scalar field is appended to the turbulence model in the local domain. This scalar field does not have its complement in the global model, giving rise to unequal number of equations at the immersed boundary between local and global models. Interdomain coupling terms are derived via the Variational Multiscale Discontinuous Galerkin (VMDG) method with new developments in scale representation and efficient fine-scale estimation. While transient laminar flows modeled with NSE in the global domain can be resolved with relatively coarse mesh, turbulent flow calculations in the local model require much finer spatial discretizations as well as smaller time-step for appropriately resolving the turbulent flow physics. The proposed framework also accommodates non-matching meshes at the immersed boundaries. Test problems in 2D and 3D numerically showcase the concurrent two-way coupling of unknown fields across the immersed boundaries. The 3D test presents a case with an unequal number of equations, where the scalar field represents the convection of contaminant concentration. This provides more detailed physics in the local region and highlights its application in climate modeling and atmospheric sciences.

Variational Multiscale Discontinuous Galerkin (VMD↗