This content will become publicly available on May 1, 2023
 Award ID(s):
 2011519
 Publication Date:
 NSFPAR ID:
 10343111
 Journal Name:
 ESAIM: Mathematical Modelling and Numerical Analysis
 Volume:
 56
 Issue:
 3
 Page Range or eLocationID:
 1007 to 1025
 ISSN:
 28227840
 Sponsoring Org:
 National Science Foundation
More Like this

Embedding properties of network realizations of dissipative reduced order models Jörn Zimmerling, Mikhail Zaslavsky,Rob Remis, Shasri Moskow, Alexander Mamonov, Murthy Guddati, Vladimir Druskin, and Liliana Borcea Mathematical Sciences Department, Worcester Polytechnic Institute https://www.wpi.edu/people/vdruskin Abstract Realizations of reduced order models of passive SISO or MIMO LTI problems can be transformed to tridiagonal and blocktridiagonal forms, respectively, via dierent modications of the Lanczos algorithm. Generally, such realizations can be interpreted as ladder resistorcapacitorinductor (RCL) networks. They gave rise to network syntheses in the rst half of the 20th century that was at the base of modern electronics design and consecutively to MOR that tremendously impacted many areas of engineering (electrical, mechanical, aerospace, etc.) by enabling ecient compression of the underlining dynamical systems. In his seminal 1950s works Krein realized that in addition to their compressing properties, network realizations can be used to embed the data back into the state space of the underlying continuum problems. In more recent works of the authors Krein's ideas gave rise to socalled nitedierence Gaussian quadrature rules (FDGQR), allowing to approximately map the ROM statespace representation to its full order continuum counterpart on a judicially chosen grid. Thus, the state variables can be accessed directly from themore »

Various dispersion models can be expressed as special cases of the Generalized Dispersion Model (GDM), which is composed of a series of Padé polynomials. While important for its broad applicability, we found that some materials with Drude dispersive terms can be accurately modeled by mixing a 1^{st}order Padé polynomial with an extra conductivity term. This conductivity term can be separated from the auxiliary differential equation (ADE). Therefore, the proposed mixedorder model can achieve the same accuracy with fewer unknowns, thus realizing higher computational efficiency and lower memory consumption. For examples, we derive the model parameters and corresponding numerical errors for noble metals including Au, Ag, and Al in the optical regime. Finally, the proposed model’s efficiency improvements are validated through implementation within a Discontinuous Galerkin Time Domain (DGTD) framework. The proposed model can achieve up to 12.5% efficiency improvement in theory compared to the conventional GDM with the same accuracy. A numerical example validates that, in practice, 9% memory reduction and 11% acceleration can be realized.

Abstract Measurements made in the Columbia River Basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP 2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modelled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and MoninObukhov similarity and (ii) a modified PriestleyTaylor (PT) algorithm with physically based ecophysiological constraints which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified PT model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP 2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulkmore »

Abstract Harnessing data to discover the underlying governing laws or equations that describe the behavior of complex physical systems can significantly advance our modeling, simulation and understanding of such systems in various science and engineering disciplines. This work introduces a novel approach called physicsinformed neural network with sparse regression to discover governing partial differential equations from scarce and noisy data for nonlinear spatiotemporal systems. In particular, this discovery approach seamlessly integrates the strengths of deep neural networks for rich representation learning, physics embedding, automatic differentiation and sparse regression to approximate the solution of system variables, compute essential derivatives, as well as identify the key derivative terms and parameters that form the structure and explicit expression of the equations. The efficacy and robustness of this method are demonstrated, both numerically and experimentally, on discovering a variety of partial differential equation systems with different levels of data scarcity and noise accounting for different initial/boundary conditions. The resulting computational framework shows the potential for closedform model discovery in practical applications where large and accurate datasets are intractable to capture.

A bstract The BrownYork stress tensor provides a means for defining quasilocal gravitational charges in subregions bounded by a timelike hypersurface. We consider the generalization of this stress tensor to null hypersurfaces. Such a stress tensor can be derived from the onshell subregion action of general relativity associated with a Dirichlet variational principle, which fixes an induced Carroll structure on the null boundary. The formula for the mixedindex tensor T i j takes a remarkably simple form that is manifestly independent of the choice of auxiliary null vector at the null surface, and we compare this expression to previous proposals for null BrownYork stress tensors. The stress tensor we obtain satisfies a covariant conservation equation with respect to any connection induced from a rigging vector at the hypersurface, as a result of the null constraint equations. For transformations that act covariantly on the boundary structures, the BrownYork charges coincide with canonical charges constructed from a version of the WaldZoupas procedure. For anomalous transformations, the charges differ by an intrinsic functional of the boundary geometry, which we explicity verify for a set of symmetries associated with finite null hypersurfaces. Applications of the null BrownYork stress tensor to symmetries of asymptotically flatmore »