An official website of the United States government
ABSTRACT Numerical stabilization techniques are often employed in under‐resolved simulations of convection‐dominated flows to improve accuracy and mitigate spurious oscillations. Specifically, the evolve–filter–relax (EFR) algorithm is a framework that consists of evolving the solution, applying a filtering step to remove high‐frequency noise, and relaxing through a convex combination of filtered and original solutions. The stability and accuracy of the EFR solution strongly depend on two parameters, the filter radius and the relaxation parameter . Standard choices for these parameters are usually fixed in time, and related to the full order model setting, that is, the grid size for and the time step for . The key novelties with respect to the standard EFR approach are: (i) time‐dependent parameters and , and (ii) data‐driven adaptive optimization of the parameters in time, considering a fully‐resolved simulation as reference. In particular, we propose three different classes of optimized‐EFR (Opt‐EFR) strategies, aiming to optimize one or both parameters. The new Opt‐EFR strategies are tested in the under‐resolved simulation of a turbulent flow past a cylinder at . The Opt‐EFR proved to be more accurate than standard approaches by up to 99, while maintaining a similar computational time. In particular, the key new finding of our analysis is that such accuracy can be obtained only if the optimized objective function includes: (i) aglobalmetric (as the kinetic energy), and (ii)spatial gradients' information.
more »
« less
This content will become publicly available on April 11, 2026
Relaxation-based schemes for on-the-fly parameter estimation in dissipative dynamical systems
Abstract This article studies two particular algorithms, a relaxation least squares algorithm and a relaxation Newton iteration scheme, for reconstructing unknown parameters in dissipative dynamical systems. Both algorithms are based on a continuous data assimilation (CDA) algorithm for state reconstruction of Azouaniet al(2014J. Nonlinear Sci.24277–304). Due to the CDA origins of these parameter recovery algorithms, these schemes provide on-the-fly reconstruction, that is, as data is collected, of unknown state and parameters simultaneously. It is shown how both algorithms give way to a robust general framework for simultaneous state and parameter estimation. In particular, we develop a general theory, applicable to a large class of dissipative dynamical systems, which identifies structural and algorithmic conditions under which the proposed algorithms achieve reconstruction of the true parameters. The algorithms are implemented on a high-dimensional two-layer Lorenz 96 model, where the theoretical conditions of the general framework are explicitly verifiable. They are also implemented on the two-dimensional Rayleigh–Bénard convection system to demonstrate the applicability of the algorithms beyond the finite-dimensional setting. In each case, systematic numerical experiments are carried out probing the efficacy of the proposed algorithms, in addition to the apparent benefits and drawbacks between them.
more »
« less
- PAR ID:
- 10599411
- Publisher / Repository:
- IOP Science
- Date Published:
- Journal Name:
- Inverse Problems
- Volume:
- 41
- Issue:
- 5
- ISSN:
- 0266-5611
- Page Range / eLocation ID:
- 055001
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
This paper develops a co-simulation framework based on the use of the package LiveLinkTMfor Matlab to perform parameters optimization of dynamical systems implemented in COMSOL Multiphysics. The identification problem is recast as an optimization problem which is solved in Matlab. Code for the key steps of the approach is described in detail, and an implementation based on the particle swarm optimization (PSO) algorithm is proposed. The effectiveness and general applicability of the framework are shown for two energy systems: lithium-ion battery (LIB) and gasoline particulate filter (GPF). Matlab codes and COMSOL models for both case studies are made publicly available and can be used as a starting point to solve parameter identification problems for systems beyond the case studies presented here.more » « less
-
Motivated by recent progress in data assimilation, we develop an algorithm to dynamically learn the parameters of a chaotic system from partial observations. Under reasonable assumptions, we supply a rigorous analytical proof that guarantees the convergence of this algorithm to the true parameter values when the system in question is the classic three-dimensional Lorenz system. Such a result appears to be the first of its kind for dynamical parameter estimation of nonlinear systems. Computationally, we demonstrate the efficacy of this algorithm on the Lorenz system by recovering any proper subset of the three non-dimensional parameters of the system, so long as a corresponding subset of the state is observable. We moreover probe the limitations of the algorithm by identifying dynamical regimes under which certain parameters cannot be effectively inferred having only observed certain state variables. In such cases, modifications to the algorithm are proposed that ultimately result in recovery of the parameter. Lastly, computational evidence is provided that supports the efficacy of the algorithm well beyond the hypotheses specified by the theorem, including in the presence of noisy observations, stochastic forcing, and the case where the observations are discrete and sparse in time.more » « less
-
This article is concerned with the problem of determining an unknown source of non-potential, external time-dependent perturbations of an incompressible fluid from large-scale observations on the flow field. A relaxation-based approach is proposed for accomplishing this, which makes use of a nonlinear property of the equations of motions to asymptotically enslave small scales to large scales. In particular, an algorithm is introduced that systematically produces approximations of the flow field on the unobserved scales in order to generate an approximation to the unknown force; the process is then repeated to generate an improved approximation of the unobserved scales, and so on. A mathematical proof of convergence of this algorithm is established in the context of the two-dimensional Navier–Stokes equations with periodic boundary conditions under the assumption that the force belongs to the observational subspace of phase space; at each stage in the algorithm, it is shown that the model error, represented as the difference between the approximating and true force, asymptotically decreases to zero in a geometric fashion provided that sufficiently many scales are observed and certain parameters of the algorithm are appropriately tuned.more » « less
-
Estimating and quantifying uncertainty in unknown system parameters from limited data remains a challenging inverse problem in a variety of real-world applications. While many approaches focus on estimating constant parameters, a subset of these problems includes time-varying parameters with unknown evolution models that often cannot be directly observed. This work develops a systematic particle filtering approach that reframes the idea behind artificial parameter evolution to estimate time-varying parameters in nonstationary inverse problems arising from deterministic dynamical systems. Focusing on systems modeled by ordinary differential equations, we present two particle filter algorithms for time-varying parameter estimation: one that relies on a fixed value for the noise variance of a parameter random walk; another that employs online estimation of the parameter evolution noise variance along with the time-varying parameter of interest. Several computed examples demonstrate the capability of the proposed algorithms in estimating time-varying parameters with different underlying functional forms and different relationships with the system states (i.e. additive vs. multiplicative).more » « less
