Search

Actuator saturation is a common nonlinearity. In wave energy conversion, force saturation conveniently limits drivetrain size and cost with minimal impact on energy generation. However, such nonlinear dynamics typically demand numerical simulation, which increases computational cost and diminishes intuition. This paper instead uses describing functions to approximate a force saturation nonlinearity as a linear impedance mismatch. In the frequency domain, the impact of controller impedance mismatch (such as force limit, finite bandwidth, or parameter error) on electrical power production is shown analytically and graphically for a generic nondimensionalized single degree of freedom wave energy converter in regular waves. Results are visualized with Smith charts. Notably, systems with a specific ratio of reactive to real mechanical impedance are least sensitive to force limits, a criteria which conflicts with resonance and bandwidth considerations. The describing function method shows promise to enable future studies such as large-scale design optimization and co-design., Comment: 6 pages, 7 figures. For code, see https://github.com/symbiotic-engineering/IFAC_CAMS_2024/ . To be presented at IFAC CAMS 2024 conference and to appear in IFAC-PapersOnLine

The series expansion of $x^m (-\log x)^l$ in terms of the shifted Chebyshev Polynomials $T_n^*(x)$ requires evaluation of the integral family $\int_0^1 x^m (-\log x)^l dx / \sqrt{x-x^2}$. We demonstrate that these can be reduced by partial integration to sums over integrals with exponent $m=0$ which have known representations as finite sums over polygamma functions., Comment: 9 pages, no figures. Integrals table in the anc directory

Direct numerical simulations (DNS) of microscale fluid-structure interactions (mFSI) in multicomponent multiphase flows pose many challenges such as the thermodynamic consistency of multiphysics couplings, the tracking of moving interfaces, and the dynamics of moving triple-phase contact lines. We propose and validate a generic DNS approach: Diffuse-Resistance-Domain (DRD). It overcomes the above challenges by employing Onsager's principle to formulate dynamic models and combining traditional diffuse-interface models for fluid-fluid interfacial dynamics with a novel implementation of complex fluid-solid interfacial conditions via smooth interpolations of dynamic-resistance coefficients across interfaces. After careful validation by benchmark simulations, we simulated three challenging mFSI problems taken from different fields., Comment: 6 pages, 4 figures

In this work, we propose an algorithm for a filter based on the Fast Fourier Transform (FFT), which, due to its characteristics, allows for an efficient computational implementation, ease of use, and minimizes amplitude variation in the filtered signal. The algorithm was implemented using the programming languages Python, R, and MATLAB. Initial results led to the conclusion that there was less amplitude loss in the filtered signal compared to the FIR filter. Future work may address a more rigorous methodology and comparative assessment of computational cost., Comment: 15 pages, 12 figures, code examples, 1 table

This brief note aims at condensing some results on the 32-point approximate DFT and discussing its arithmetic complexity., Comment: Corrected a typo in (4). 8 pages, 2 tables

We consider the fractional influence maximization problem, i.e., identifying users on a social network to be incentivized with potentially partial discounts to maximize the influence on the network. The larger the discount given to a user, the higher the likelihood of its activation (adopting a new product or innovation), who then attempts to activate its neighboring users, causing a cascade effect of influence through the network. Our goal is to devise efficient algorithms that assign initial discounts to the network's users to maximize the total number of activated users at the end of the cascade, subject to a constraint on the total sum of discounts given. In general, the activation likelihood could be any non-decreasing function of the discount, whereas, our focus lies on the case when the activation likelihood is an affine function of the discount, potentially varying across different users. As this problem is shown to be NP-hard, we propose and analyze an efficient (1-1/e)-approximation algorithm. Furthermore, we run experiments on real-world social networks to show the performance and scalability of our method., Comment: 5 pages, 2 figures, and 1 table

Off-policy evaluation (OPE) is crucial for evaluating a target policy's impact offline before its deployment. However, achieving accurate OPE in large state spaces remains challenging.This paper studies state abstractions-originally designed for policy learning-in the context of OPE. Our contributions are three-fold: (i) We define a set of irrelevance conditions central to learning state abstractions for OPE. (ii) We derive sufficient conditions for achieving irrelevance in Q-functions and marginalized importance sampling ratios, the latter obtained by constructing a time-reversed Markov decision process (MDP) based on the observed MDP. (iii) We propose a novel two-step procedure that sequentially projects the original state space into a smaller space, which substantially simplify the sample complexity of OPE arising from high cardinality., Comment: 42 pages, 5 figures

In portfolio risk minimization, the inverse covariance matrix of returns is often unknown and has to be estimated in practice. This inverse covariance matrix also prescribes the hedge trades in which a stock is hedged by all the other stocks in the portfolio. In practice with finite samples, however, multicollinearity gives rise to considerable estimation errors, making the hedge trades too unstable and unreliable for use. By adopting ideas from current methodologies in the existing literature, we propose 2 new estimators of the inverse covariance matrix, one which relies only on the l2 norm while the other utilizes both the l1 and l2 norms. These 2 new estimators are classified as shrinkage estimators in the literature. Comparing favorably with other methods (sample-based estimation, equal-weighting, estimation based on Principal Component Analysis), a portfolio formed on the proposed estimators achieves substantial out-of-sample risk reduction and improves the out-of-sample risk-adjusted returns of the portfolio, particularly in high-dimensional settings. Furthermore, the proposed estimators can still be computed even in instances where the sample covariance matrix is ill-conditioned or singular, Comment: B.Sc. Dissertation

This paper introduces a novel method for the stability analysis of positive feedback systems with a class of fully connected feedforward neural networks (FFNN) controllers. By establishing sector bounds for fully connected FFNNs without biases, we present a stability theorem that demonstrates the global exponential stability of linear systems under fully connected FFNN control. Utilizing principles from positive Lur'e systems and the positive Aizerman conjecture, our approach effectively addresses the challenge of ensuring stability in highly nonlinear systems. The crux of our method lies in maintaining sector bounds that preserve the positivity and Hurwitz property of the overall Lur'e system. We showcase the practical applicability of our methodology through its implementation in a linear system managed by a FFNN trained on output feedback controller data, highlighting its potential for enhancing stability in dynamic systems., Comment: 6 pages, 7 figures, to be published in IEEE Control Systems Letters (L-CSS)

We calculate quasinormal modes of scalar and neutrino perturbations around the charged black hole in Asymptotically Safe Gravity. We show that the charge and coupling constant change the quasinormal spectrum considerably. We show that previous calculations of scalar quasinormal modes in this background [F. Javed, Phys. Dark Univ., 44, 101450 (2024)] suffer from a large numerical error exceeding the effect, that is, the deviations of the frequencies from their Schwarzschild limits. In the high frequency (eikonal) limit an explicit analytic formula for quasinormal modes is derived, which confirms the correspondence between the null circular geodesics and eikonal quasinormal frequencies.

In this work, the concept of Braced Fourier Continuation and Regression (BFCR) is introduced. BFCR is a novel and computationally efficient means of finding nonlinear regressions or trend lines in arbitrary one-dimensional data sets. The Braced Fourier Continuation (BFC) and BFCR algorithms are first outlined, followed by a discussion of the properties of BFCR as well as demonstrations of how BFCR trend lines may be used effectively for anomaly detection both within and at the edges of arbitrary one-dimensional data sets. Finally, potential issues which may arise while using BFCR for anomaly detection as well as possible mitigation techniques are outlined and discussed. All source code and example data sets are either referenced or available via GitHub, and all associated code is written entirely in Python., Comment: 16 pages, 9 figures, associated Github link: https://github.com/j4sabuda/Braced-Fourier-Continuation-and-Regression -6/30/2024 update corrected and reworded erroneous figure references, minor typos

Multifidelity surrogate modelling combines data of varying accuracy and cost from different sources. It strategically uses low-fidelity models for rapid evaluations, saving computational resources, and high-fidelity models for detailed refinement. It improves decision-making by addressing uncertainties and surpassing the limits of single-fidelity models, which either oversimplify or are computationally intensive. Blending high-fidelity data for detailed responses with frequent low-fidelity data for quick approximations facilitates design optimisation in various domains. Despite progress in interpolation, regression, enhanced sampling, error estimation, variable fidelity, and data fusion techniques, challenges persist in selecting fidelity levels and developing efficient data fusion methods. This study proposes a new fusion approach to construct multi-fidelity surrogate models by constructing gradient-only surrogates that use only gradients to construct regression surfaces. Results are demonstrated on foundational example problems that isolate and illustrate the fusion approach's efficacy, avoiding the need for complex examples that obfuscate the main concept., Comment: 8 pages, 4 figures, SACAM2024 Conference, 22-23 January 2024

In this work we provide theoretical estimates for the ranks of the power functions $f(k) = k^{-\alpha}$, $\alpha>1$ in the quantized tensor train (QTT) format for $k = 1, 2, 3, \ldots, 2^{d}$. Such functions and their several generalizations (e.~g. $f(k) = k^{-\alpha} \cdot e^{-\lambda k}, \lambda > 0$) play an important role in studies of the asymptotic solutions of the aggregation-fragmentation kinetic equations. In order to support the constructed theory we verify the values of QTT-ranks of these functions in practice with the use of the TTSVD procedure and show an agreement between the numerical and analytical results., Comment: 6 pages, 1 figure, 20 references

Birthday problem is a well-known classic problem in probability theory widely applied in cryptography. Although bubble sort is a popular algorithm leading to some interesting theoretical problems in computer science, its relation to birthday problem has not been found yet. This paper indicates how Rayleigh distribution naturally arises in bubble sort by relating it to birthday problem, which presents a novel direction for analysing bubble sort and birthday problem. Then asymptotic distributions and statistical characteristics of bubble sort and birthday problem with very small absolute errors are presented. Moreover, this paper proves that some common optimizations of bubble sort could lead to average performance degradation., Comment: 8 pages, 2 figures

We consider the 2D acoustic system with the Gaussian pulse as the initial data. This case was proposed at the first Workshop on benchmark problems in computational aeroacoustics, and it is commonly used for the verification of numerical methods. We construct an efficient algorithm to evaluate the exact solution for a given time t and distance r. For a precision eps, it takes c*ln(1/eps) operations (the evaluation of a Bessel function counts as one operation) where c does not depend on t and r. This becomes possible by using three different integral representations and an asymptotic series depending on t and r.

We propose an efficient implementation of the numerical tensor-train (TT) based algorithm solving the multicomponent coagulation equation preserving the nonnegativeness of solution. Unnatural negative elements in the constructed approximation arise due to the errors of the low-rank decomposition and discretization scheme. In this work, we propose to apply the rank-one corrections in the TT-format proportional to the minimal negative element. Such an element can be found via application of the global optimization methods that can be fully implemented within efficient operations in the tensor train format. We incorporate this trick into the time-integration scheme for the multicomponent coagulation equation and also use it for post-processing of the stationary solution for the problem with the source of particles., Comment: 12 pages, 4 figures, 4 tables

We develop a novel identification strategy as well as a new estimator for context-dependent causal inference in non-parametric triangular models with non-separable disturbances. Departing from the common practice, our analysis does not rely on the strict monotonicity assumption. Our key contribution lies in leveraging on diffusion models to formulate the structural equations as a system evolving from noise accumulation to account for the influence of the latent context (confounder) on the outcome. Our identifiability strategy involves a system of Fredholm integral equations expressing the distributional relationship between a latent context variable and a vector of observables. These integral equations involve an unknown kernel and are governed by a set of structural form functions, inducing a non-monotonic inverse problem. We prove that if the kernel density can be represented as an infinite mixture of Gaussians, then there exists a unique solution for the unknown function. This is a significant result, as it shows that it is possible to solve a non-monotonic inverse problem even when the kernel is unknown. On the methodological front we leverage on a novel and enriched Contaminated Generative Adversarial (Neural) Networks (CONGAN) which we provide as a solution to the non-monotonic inverse problem.

This note aims to present novel positive linear operators involving the Wright function. Furthermore, the present research established the moments of these newly defined operators and estimated the convergence rate using the classical modulus of continuity. Additionally, the convergence rate in the Lipschitz spaces, their A-statistical convergence property, and $({\lambda}, {\gamma})$-statistical convergence properties have been covered., Comment: 18

We introduce cppdlr, a C++ library implementing the discrete Lehmann representation (DLR) of functions in imaginary time and Matsubara frequency, such as Green's functions and self-energies. The DLR is based on a low-rank approximation of the analytic continuation kernel, and yields a compact and explicit basis consisting of exponentials in imaginary time and simple poles in Matsubara frequency. cppdlr constructs the DLR basis and associated interpolation grids, and implements standard operations. It provides a flexible yet high-level interface, facilitating the incorporation of the DLR into both small-scale applications and existing large-scale software projects.

In the classical best approximation pair (BAP) problem, one is given two nonempty, closed, convex and disjoint subsets in a finite- or an infinite-dimensional Hilbert space, and the goal is to find a pair of points, each from each subset, which realizes the distance between the subsets. This problem, which has a long history, has found applications in science and technology. We discuss the problem in more general normed spaces and with possibly non-convex subsets, and focus our attention on the issues of uniqueness and existence of the solution to the problem. To the best of our knowledge these fundamental issues have not received much attention. In particular, we present several sufficient geometric conditions for the (at most) uniqueness of a BAP relative to these subsets. These conditions are related to the structure of the boundaries of the subsets, their relative orientation, and the structure of the unit sphere of the space. In addition, we present many sufficient conditions for the existence of a BAP, possibly without convexity . Our results allow us to significantly extend the horizon of the recent alternating simultaneous Halpern-Lions-Wittmann-Bauschke (A-S-HLWB) algorithm [Censor, Mansour and Reem, The alternating simultaneous Halpern-Lions-Wittmann-Bauschke algorithm for finding the best approximation pair for two disjoint intersections of convex sets, arXiv:2304.09600 (2023)] for solving the BAP problem.

This work highlights an approach for incorporating realistic uncertainties into scientific computing workflows based on finite elements, focusing on applications in computational mechanics and design optimization. We leverage Mat\'ern-type Gaussian random fields (GRFs) generated using the SPDE method to model aleatoric uncertainties, including environmental influences, variating material properties, and geometric ambiguities. Our focus lies on delivering practical GRF realizations that accurately capture imperfections and variations and understanding how they impact the predictions of computational models and the topology of optimized designs. We describe a numerical algorithm based on solving a generalized SPDE to sample GRFs on arbitrary meshed domains. The algorithm leverages established techniques and integrates seamlessly with the open-source finite element library MFEM and associated scientific computing workflows, like those found in industrial and national laboratory settings. Our solver scales efficiently for large-scale problems and supports various domain types, including surfaces and embedded manifolds. We showcase its versatility through biomechanics and topology optimization applications. The flexibility and efficiency of SPDE-based GRF generation empower us to run large-scale optimization problems on 2D and 3D domains, including finding optimized designs on embedded surfaces, and to generate topologies beyond the reach of conventional techniques. Moreover, these capabilities allow us to model geometric uncertainties of reconstructed submanifolds, such as the surfaces of cerebral aneurysms. In addition to offering benefits in these specific domains, the proposed techniques transcend specific applications and generalize to arbitrary forward and backward problems in uncertainty quantification involving finite elements., Comment: 38 pages, 21 figures

This report elaborates on approximations for the discrete Fourier transform by means of replacing the exact Cooley-Tukey algorithm twiddle-factors by low-complexity integers, such as $0, \pm \frac{1}{2}, \pm 1$., Comment: Originally written in Oct 2016. Contains 55 pages, 17 figures, 4 tables

An approach to modeling the oligomer composition distribution function in the irreversible step growth homopolymerization process based on a mixture of oligomers of arbitrary composition is developed. The approach is based on consideration of probabilities of processes in the system, proceeding on the principle of P. Flory and obtaining on this basis an infinite system of differential equations linking mole fractions of each oligomer of a finite mixture with the degree of transformation., Comment: Libre Office.org 7.5, pdf, 19 pages, 3 figures

The square root velocity transformation provides a convenient and numerically efficient approach to functional and shape data analysis of curves. We study fundamental geometric properties of curves under this transformation. Moreover, utilizing natural geometric constructions, we employ the approach for intrinsic comparison within several classes of surfaces and augmented curves, which arise in the real world applications such as tubes, ruled surfaces, spherical strips, protein molecules and hurricane tracks.

In conventional finite element simulations, foil windings with a thin foil and many turns require many mesh elements. This renders models quickly computationally infeasible. With the use of homogenization approaches, the finite element mesh does not need to resolve the small-scale structure of the foil winding domain. Present homogenization approaches take resistive and inductive effects into account. With an increase of the operation frequency of foil windings, however, capacitive effects between adjacent turns in the foil winding become relevant. This paper presents an extension to the standard foil winding model that covers the capacitive behavior of foil windings., Comment: 8 pages, 12 figures