Search

This article presents an innovative approach for developing an efficient reduced-order model to study the dispersion of urban air pollutants. The need for real-time air quality monitoring has become increasingly important, given the rise in pollutant emissions due to urbanization and its adverse effects on human health. The proposed methodology involves solving the linear advection-diffusion problem, where the solution of the Reynolds-averaged Navier-Stokes equations gives the convective field. At the same time, the source term consists of an empirical time series. However, the computational requirements of this approach, including microscale spatial resolution, repeated evaluation, and low time scale, necessitate the use of high-performance computing facilities, which can be a bottleneck for real-time monitoring. To address this challenge, a problem-specific methodology was developed that leverages a data-driven approach based on Proper Orthogonal Decomposition with regression (POD-R) coupled with Galerkin projection (POD-G) endorsed with the discrete empirical interpolation method (DEIM). The proposed method employs a feedforward neural network to non-intrusively retrieve the reduced-order convective operator required for online evaluation. The numerical framework was validated on synthetic emissions and real wind measurements. The results demonstrate that the proposed approach significantly reduces the computational burden of the traditional approach and is suitable for real-time air quality monitoring. Overall, the study advances the field of reduced order modeling and highlights the potential of data-driven approaches in environmental modeling and large-scale simulations.

We derive an analytical expression that describes the complete three-dimensional carrier-envelope phase (CEP) distribution of in the focal volume of ultrashort pulsed Gaussian beams focused by spherical mirrors or lenses. The focal CEP map depends on the so-called factor $g$ specifying the frequency-dependence of the beam width of the source few-cycle pulse, on its chirp and on the small chromatic aberration introduced by a lens without appreciably distorting or broadening the few-cycle pulse. We show how to tailor the CEP map of mirror-focused and lens-focused few-cycle pulses in order to produce negligible transversal and axial CEP variations in specific regions of the focal volume for phase-sensitive interactions of light with matter taking place in a volume or on a surface. We propose a quasi-achromatic doublet lens that can implement in practice these tailored CEP distributions., Comment: 9 pages, 6 figures

Widely used complex code refactoring tools lack a solid reasoning about the correctness of the transformations they implement, whilst interest in proven correct refactoring is ever increasing as only formal verification can provide true confidence in applying tool-automated refactoring to industrial-scale code. By using our strategic rewriting based refactoring specification language, we present the decomposition of a complex transformation into smaller steps that can be expressed as instances of refactoring schemes, then we demonstrate the semi-automatic formal verification of the components based on a theoretical understanding of the semantics of the programming language. The extensible and verifiable refactoring definitions can be executed in our interpreter built on top of a static analyser framework., Comment: In Proceedings VPT 2017, arXiv:1708.06887

Recently, Horv\'ath, Song, and Terlaky [\emph{A novel unified approach to invariance condition of dynamical system, submitted to Applied Mathematics and Computation}] proposed a novel unified approach to study, i.e., invariance conditions, sufficient and necessary conditions, under which some convex sets are invariant sets for linear dynamical systems. In this paper, by utilizing analogous methodology, we generalize the results for nonlinear dynamical systems. First, the Theorems of Alternatives, i.e., the nonlinear Farkas lemma and the \emph{S}-lemma, together with Nagumo's Theorem are utilized to derive invariance conditions for discrete and continuous systems. Only standard assumptions are needed to establish invariance of broadly used convex sets, including polyhedral and ellipsoidal sets. Second, we establish an optimization framework to computationally verify the derived invariance conditions. Finally, we derive analogous invariance conditions without any conditions.

In this paper, we consider local and uniform invariance preserving steplength thresholds on a set when a discretization method is applied to a linear or nonlinear dynamical system. For the forward or backward Euler method, the existence of local and uniform invariance preserving steplength thresholds is proved when the invariant sets are polyhedra, ellipsoids, or Lorenz cones. Further, we also quantify the steplength thresholds of the backward Euler methods on these sets for linear dynamical systems. Finally, we present our main results on the existence of uniform invariance preserving steplength threshold of general discretization methods on general convex sets, compact sets, and proper cones both for linear and nonlinear dynamical systems.

Steplength thresholds for invariance preserving of three types of discretization methods on a polyhedron are considered. For Taylor approximation type discretization methods we prove that a valid steplength threshold can be obtained by finding the first positive zeros of a finite number of polynomial functions. Further, a simple and efficient algorithm is proposed to numerically compute the steplength threshold. For rational function type discretization methods we derive a valid steplength threshold for invariance preserving, which can be computed by using an analogous algorithm as in the first case. The relationship between the previous two types of discretization methods and the forward Euler method is studied. Finally, we show that, for the forward Euler method, the largest steplength threshold for invariance preserving can be computed by solving a finite number of linear optimization problems.

In this paper, we propose a novel, unified, general approach to investigate sufficient and necessary conditions under which four types of convex sets, polyhedra, polyhedral cones, ellipsoids and Lorenz cones, are invariant sets for a linear continuous or discrete dynamical system. In proving invariance of ellipsoids and Lorenz cones for discrete systems, instead of the traditional Lyapunov method, our novel proofs are based on the S-lemma, which enables us to extend invariance conditions to any set represented by a quadratic inequality. Such sets include nonconvex and unbounded sets. Finally, according to the framework of our novel method, sufficient and necessary conditions for continuous systems are derived from the sufficient and necessary conditions for the corresponding discrete systems that are obtained by Euler methods.