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.


Title: On the reconstruction of unknown driving forces from low-mode observations in the 2D Navier–Stokes equations
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
Award ID(s):
2213363
PAR ID:
10523365
Author(s) / Creator(s):
Publisher / Repository:
Proc. R. Soc. Edinb. A: Math.
Date Published:
Journal Name:
Proceedings of the Royal Society of Edinburgh: Section A Mathematics
ISSN:
0308-2105
Page Range / eLocation ID:
1 to 24
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Experimental design is a classical statistics problem, and its aim is to estimate an unknown vector from linear measurements where a Gaussian noise is introduced in each measurement. For the combinatorial experimental design problem, the goal is to pick a subset of experiments so as to make the most accurate estimate of the unknown parameters. In this paper, we will study one of the most robust measures of error estimation—the D-optimality criterion, which corresponds to minimizing the volume of the confidence ellipsoid for the estimation error. The problem gives rise to two natural variants depending on whether repetitions of experiments are allowed or not. We first propose an approximation algorithm with a 1/e-approximation for the D-optimal design problem with and without repetitions, giving the first constant-factor approximation for the problem. We then analyze another sampling approximation algorithm and prove that it is asymptotically optimal. Finally, for D-optimal design with repetitions, we study a different algorithm proposed by the literature and show that it can improve this asymptotic approximation ratio. All the sampling algorithms studied in this paper are shown to admit polynomial-time deterministic implementations. 
    more » « less
  2. We consider a large-scale parallel-server loss system with an unknown arrival rate, where each server is able to adjust its processing speed. The objective is to minimize the system cost, which consists of a power cost to maintain the servers' processing speeds and a quality of service cost depending on the tasks' processing times, among others. We draw on ideas from stochastic approximation to design a novel speed scaling algorithm and prove that the servers' processing speeds converge to the globally asymptotically optimum value. Curiously, the algorithm is fully distributed and does not require any communication between servers. Apart from the algorithm design, a key contribution of our approach lies in demonstrating how concepts from the stochastic approximation literature can be leveraged to effectively tackle learning problems in large-scale, distributed systems. En route, we also analyze the performance of a fully heterogeneous parallel-server loss system, where each server has a distinct processing speed, which might be of independent interest. 
    more » « less
  3. Polymeric materials that couple deformation and electrostatics have the potential for use in soft sensors and actuators with potential applications ranging from robotic, biomedical, energy, aerospace and automotive technologies. In contrast to the mechanics of polymers that has been studied using statistical mechanics approaches for decades, the coupled response under deformation and electrical field has largely been modeled only phenomenologically at the continuum scale. In this work, we examine the physics of the coupled deformation and electrical response of an electrically-responsive polymer chain using statistical mechanics. We begin with a simple anisotropic model for the electrostatic dipole response to electric field of a single monomer, and use a separation of energy scales between the electrostatic field energy and the induced dipole field energy to reduce the nonlocal and infinite-dimensional statistical averaging to a simpler local finite-dimensional averaging. In this simplified setting, we derive the equations of the most likely monomer orientation density using the maximum term approximation, and a chain free energy is derived using this approximation. These equations are investigated numerically and the results provide insight into the physics of electro-mechanically coupled elastomer chains. Closed-form approximations are also developed in the limit of small electrical energy with respect to thermal energy; in the limit of small mechanical tension force acting on the chain; and using asymptotic matching for general chain conditions. 
    more » « less
  4. There is a wide variety of applications that require sorting and separation of micro- particles from a large cluster of similar objects. Existing methods can distinguish micro-particles by their bulk properties, such as their size, density, and electric polarizability. These methods, however, are not selective with respect to the individual geometry of the particles. In this work, we focus on the use of a resonance effect between a microparticle and an evanescent light field known as the Whispering Gallery Mode (WGM) force. The WGM force is highly sensitive to the radius of the particle and is both controllable and tunable. In this paper, we explore through simulation the design of a WGM-based device for micro-particle separation. In this device, particles flow in through an inlet and are carried over two actuation regions given by waveguides carrying laser light to generate the evanescent field. Particles are observed by a camera, allowing for feedback control on the power of the lasers. While the basic control structure is simple, there are several challenges, including unknown disturbances to the fluid flow, limited laser power, and uni-directional control over each actuation region. We combine Expectation Maximization with Kalman filtering to both estimate the unknown disturbance and filter the measurements into a position estimate. We then develop simple hybrid controllers and compare them to the ideal setting (without any constraints) based on a Linear–Quadratic–Gaussian (LQG) control approach. 
    more » « less
  5. We investigate the motion of a thin vortex filament in the presence of buoyancy. The asymptotic model of Moore & Saffman ( Phil. Trans. R. Soc. Lond. A, vol. 272, 1972, pp. 403–429) is extended to take account of buoyancy forces in the force balance on a vortex element. The motion of a buoyant vortex is given by the transverse component of force balance, while the tangential component governs the dynamics of the structure in the core. We show that the local acceleration of axial flow is generated by the external pressure gradient due to gravity. The equations are then solved for vortex rings. An analytic solution for a buoyant vortex ring at a small initial inclination is obtained and asymptotically agrees with the literature. 
    more » « less