skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 1, 2026

Title: Residual-based data-driven variational multiscale reduced order models for parameter-dependent problems
Abstract In this paper, we propose a novel residual-based data-driven closure strategy for reduced order models (ROMs) of under-resolved, convection-dominated problems. The new ROM closure model is constructed in a variational multiscale (VMS) framework by using the available full order model data and a model form ansatz that depends on the ROM residual. We emphasize that this closure modeling strategy is fundamentally different from the current data-driven ROM closures, which generally depend on the ROM coefficients. We investigate the new residual-based data-driven VMS ROM closure strategy in the numerical simulation of three test problems: (i) a one-dimensional parameter-dependent advection-diffusion problem; (ii) a two-dimensional time-dependent advection-diffusion-reaction problem with a small diffusion coefficient ($$\varepsilon = 1e-4$$ ε = 1 e - 4 ); and (iii) a two-dimensional flow past a cylinder at Reynolds number$$Re=1000$$ R e = 1000 . Our numerical investigation shows that the new residual-based data-driven VMS-ROM is more accurate than the standard coefficient-based data-driven VMS-ROM.  more » « less
Award ID(s):
2012253
PAR ID:
10628758
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer
Date Published:
Journal Name:
Computational and Applied Mathematics
Volume:
44
Issue:
6
ISSN:
2238-3603
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Let$$\mathbb {F}_q^d$$ F q d be thed-dimensional vector space over the finite field withqelements. For a subset$$E\subseteq \mathbb {F}_q^d$$ E F q d and a fixed nonzero$$t\in \mathbb {F}_q$$ t F q , let$$\mathcal {H}_t(E)=\{h_y: y\in E\}$$ H t ( E ) = { h y : y E } , where$$h_y:E\rightarrow \{0,1\}$$ h y : E { 0 , 1 } is the indicator function of the set$$\{x\in E: x\cdot y=t\}$$ { x E : x · y = t } . Two of the authors, with Maxwell Sun, showed in the case$$d=3$$ d = 3 that if$$|E|\ge Cq^{\frac{11}{4}}$$ | E | C q 11 4 andqis sufficiently large, then the VC-dimension of$$\mathcal {H}_t(E)$$ H t ( E ) is 3. In this paper, we generalize the result to arbitrary dimension by showing that the VC-dimension of$$\mathcal {H}_t(E)$$ H t ( E ) isdwhenever$$E\subseteq \mathbb {F}_q^d$$ E F q d with$$|E|\ge C_d q^{d-\frac{1}{d-1}}$$ | E | C d q d - 1 d - 1
    more » « less
  2. Abstract We explore the decay of bound neutrons in the JUNO liquid scintillator detector into invisible particles (e.g.,$$n\rightarrow 3 \nu $$ n 3 ν or$$nn \rightarrow 2 \nu $$ n n 2 ν ), which do not produce an observable signal. The invisible decay includes two decay modes:$$ n \rightarrow { inv} $$ n inv and$$ nn \rightarrow { inv} $$ n n inv . The invisible decays ofs-shell neutrons in$$^{12}\textrm{C}$$ 12 C will leave a highly excited residual nucleus. Subsequently, some de-excitation modes of the excited residual nuclei can produce a time- and space-correlated triple coincidence signal in the JUNO detector. Based on a full Monte Carlo simulation informed with the latest available data, we estimate all backgrounds, including inverse beta decay events of the reactor antineutrino$${\bar{\nu }}_e$$ ν ¯ e , natural radioactivity, cosmogenic isotopes and neutral current interactions of atmospheric neutrinos. Pulse shape discrimination and multivariate analysis techniques are employed to further suppress backgrounds. With two years of exposure, JUNO is expected to give an order of magnitude improvement compared to the current best limits. After 10 years of data taking, the JUNO expected sensitivities at a 90% confidence level are$$\tau /B( n \rightarrow { inv} ) > 5.0 \times 10^{31} \, \textrm{years}$$ τ / B ( n inv ) > 5.0 × 10 31 years and$$\tau /B( nn \rightarrow { inv} ) > 1.4 \times 10^{32} \, \textrm{years}$$ τ / B ( n n inv ) > 1.4 × 10 32 years
    more » « less
  3. Abstract Given a compact doubling metric measure spaceXthat supports a 2-Poincaré inequality, we construct a Dirichlet form on$$N^{1,2}(X)$$ N 1 , 2 ( X ) that is comparable to the upper gradient energy form on$$N^{1,2}(X)$$ N 1 , 2 ( X ) . Our approach is based on the approximation ofXby a family of graphs that is doubling and supports a 2-Poincaré inequality (see [20]). We construct a bilinear form on$$N^{1,2}(X)$$ N 1 , 2 ( X ) using the Dirichlet form on the graph. We show that the$$\Gamma $$ Γ -limit$$\mathcal {E}$$ E of this family of bilinear forms (by taking a subsequence) exists and that$$\mathcal {E}$$ E is a Dirichlet form onX. Properties of$$\mathcal {E}$$ E are established. Moreover, we prove that$$\mathcal {E}$$ E has the property of matching boundary values on a domain$$\Omega \subseteq X$$ Ω X . This construction makes it possible to approximate harmonic functions (with respect to the Dirichlet form$$\mathcal {E}$$ E ) on a domain inXwith a prescribed Lipschitz boundary data via a numerical scheme dictated by the approximating Dirichlet forms, which are discrete objects. 
    more » « less
  4. Abstract We investigate the low moments$$\mathbb {E}[|A_N|^{2q}],\, 0 E [ | A N | 2 q ] , 0 < q 1 of secular coefficients$$A_N$$ A N of the critical non-Gaussian holomorphic multiplicative chaos, i.e. coefficients of$$z^N$$ z N in the power series expansion of$$\exp (\sum _{k=1}^\infty X_kz^k/\sqrt{k})$$ exp ( k = 1 X k z k / k ) , where$$\{X_k\}_{k\geqslant 1}$$ { X k } k 1 are i.i.d. rotationally invariant unit variance complex random variables. Inspired by Harper’s remarkable result on random multiplicative functions, Soundararajan and Zaman recently showed that if each$$X_k$$ X k is standard complex Gaussian,$$A_N$$ A N features better-than-square-root cancellation:$$\mathbb {E}[|A_N|^2]=1$$ E [ | A N | 2 ] = 1 and$$\mathbb {E}[|A_N|^{2q}]\asymp (\log N)^{-q/2}$$ E [ | A N | 2 q ] ( log N ) - q / 2 for fixed$$q\in (0,1)$$ q ( 0 , 1 ) as$$N\rightarrow \infty $$ N . We show that this asymptotics holds universally if$$\mathbb {E}[e^{\gamma |X_k|}]<\infty $$ E [ e γ | X k | ] < for some$$\gamma >2q$$ γ > 2 q . As a consequence, we establish the universality for the tightness of the normalized secular coefficients$$A_N(\log (1+N))^{1/4}$$ A N ( log ( 1 + N ) ) 1 / 4 , generalizing a result of Najnudel, Paquette, and Simm. Another corollary is the almost sure regularity of some critical non-Gaussian holomorphic chaos in appropriate Sobolev spaces. Moreover, we characterize the asymptotics of$$\mathbb {E}[|A_N|^{2q}]$$ E [ | A N | 2 q ] for$$|X_k|$$ | X k | following a stretched exponential distribution with an arbitrary scale parameter, which exhibits a completely different behavior and underlying mechanism from the Gaussian universality regime. As a result, we unveil a double-layer phase transition around the critical case of exponential tails. Our proofs combine Harper’s robust approach with a careful analysis of the (possibly random) leading terms in the monomial decomposition of$$A_N$$ A N
    more » « less
  5. Noble element time projection chambers are a leading technology for rare event detection in physics, such as for dark matter and neutrinoless double beta decay searches. Time projection chambers typically assign event position in the drift direction using the relative timing of prompt scintillation and delayed charge collection signals, allowing for reconstruction of an absolute position in the drift direction. In this paper, alternate methods for assigning event drift distance via quantification of electron diffusion in a pure high pressure xenon gas time projection chamber are explored. Data from the NEXT-White detector demonstrate the ability to achieve good position assignment accuracy for both high- and low-energy events. Using point-like energy deposits from$$^{83\textrm{m}}$$ 83 m Kr calibration electron captures ($$E\sim 45$$ E 45  keV), the position of origin of low-energy events is determined to 2 cm precision with bias$$< 1~$$ < 1 mm. A convolutional neural network approach is then used to quantify diffusion for longer tracks ($$E\ge ~1.5$$ E 1.5  MeV), from radiogenic electrons, yielding a precision of 3 cm on the event barycenter. The precision achieved with these methods indicates the feasibility energy calibrations of better than 1% FWHM at Q$$_{\beta \beta }$$ β β in pure xenon, as well as the potential for event fiducialization in large future detectors using an alternate method that does not rely on primary scintillation. 
    more » « less