Free Energy Wells and Overlap Gap Property in Sparse PCA

We study a variant of the sparse PCA (principal component analysis) problem in the “hard” regime, where the inference task is possible yet no polynomial-time algorithm is known to exist. Prior work, based on the low-degree likelihood ratio, has conjectured a precise expression for the best possible (subexponential) runtime throughout the hard regime. Following instead a statistical physics-inspired point of view, we show bounds on the depth of free energy wells for various Gibbs measures naturally associated to the problem. These free energy wells imply hitting time lower bounds that corroborate the low-degree conjecture: we show that a class of natural MCMC (Markov chain Monte Carlo) methods (with worst-case initialization) cannot solve sparse PCA with less than the conjectured runtime. These lower bounds apply to a wide range of values for two tuning parameters: temperature and sparsity misparametrization. Finally, we prove that the overlap gap property (OGP), a structural property that implies failure of certain local search algorithms, holds in a significant part of the hard regime. © 2022 Wiley Periodicals, Inc.     

Time resolved measurements of the flow generated by suction feeding fish

The majority of aquatic vertebrates are suction feeders: by rapidly expanding the mouth cavity they generate a fluid flow outside of their head in order to draw prey into their mouth. In addition to the biological relevance, the generated flow field is interesting fluid mechanically as it incorporates high velocities, is localized in front of the mouth, and is unsteady, typically lasting between 10 and 50 ms. Using manometry and high-speed particle image velocimetry, this is the first study to quantify pressure within and outside the mouth of a feeding fish while simultaneously measuring the velocity field outside the mouth. Measurements with a high temporal (2 ms) and spatial (<1 mm) resolution were made for several feeding events of a single largemouth bass (Micropterus salmoides). General properties of the flow were evaluated, including the transient velocity field, its relationship to pressure within the mouth and pressure at the prey. We find that throughout the feeding event a relationship exists for the magnitude of fluid speed as a function of distance from the predator mouth that is based on scaling the velocity field according to the size of the mouth opening and the magnitude of fluid speed at the mouth. The velocity field is concentrated within an area extending approximately one mouth diameter from the fish and the generated pressure field is even more local to the mouth aperture. Although peak suction pressures measured inside the mouth were slightly larger than those that were predicted using the equations of motion, we find that these equations give a very accurate prediction of the timing of peak pressure, so long as the unsteady nature of the flow is included.     

Il'iushin's postulate and resulting thermodynamic conditions on elasto-plastic coupling

Il'iushin's postulate is restated within a general thermodynamic strain space formulation of rate independent plasticity by means of plastic internal variables. This yields a general expression in terms of appropriate thermodynamic potentials. A combination of a thermodynamic condition, derived from the general development, with the results of Il'iushin's postulate, furnishes explicit conditions on elasto-plastic coupling. A specific example is presented, with the plastic work being the only plastic internal variable. Necessary and sufficient consitbns on the elastic moduli and their change with plastic deformation are derived, for the thermodynamic condition to be satisfied.     

High-intensity turbulence measurements in a Taylor-Couette flow reactor

Turbulence-intensity measurements were made in a Taylor-Couette flow reactor consisting of two counter-rotating concentric cylinders designed for the purpose of studying turbulent premixed-flame propagation. In the annulus separating the two cylinders, a nearly homogeneous turbulent flow is generated. The intensities of turbulent velocity fluctuations in the annulus in both axial and circumferential directions were measured by using laser-Doppler velocimetry for a wide range of cylinder rotation rates, corresponding to low through high (120 cm/s) intensities relative to typical laminar flame speeds for lean methane-air mixtures. The experimental measurements indicate a linear relation between turbulence intensities and average cylinder surface speed and demonstrate the usefulness of the Taylor-Couette apparatus for studies of premixed-flame propagation in high-intensity turbulent flow.     

Elasto-plastic coupling within a thermodynamic strain space formulation of plasticity

Elasto-plastic coupling is studied within a general thermodynamic train space formulation of rate-independent plasticity by means of plastic internal variables. The strain space formulation offers a unified approach in treating both stable and unstable material behavior simultaneously. The conditions on elasto-plastic coupling are imposed by the second law of thermodynamics and the consistency equation in the strain space formulation, and are expressed by two inequalities. This is further illustrated by specific examples where explicit necessary and sufficient conditions on the elastic moduli and their change with plastic deformation are derived for the two inequalities to be satisfied.     

Aerosol formation in an isothermal stagnation flow

An experiment has been performed in a laminar stagnation point flow in which two non-premixed reactants produced an aerosol of sub-micron particles. The reactants were NH₃ and HCl. The rate of mixing of the reactants was determined by the velocity gradient or strain rate of the flow; the response of the aerosol dynamics to the flow field was measured with a laser light scattering technique. Laser Doppler Spectroscopy was used to measure the particle size. It was found that the particle size was independent of the strain rate of the flow. On the other hand, the particle number density decreased as the strain rate increased. It is argued that the intensity of light scattered from the aerosol is, therefore, a measure of the amount of product of the relatively slow NH₃-HCl reaction.     

The Forchheimer equation: A theoretical development

In this paper we illustrate how the method of volume averaging can be used to derive Darcy's law with the Forchheimer correction for homogeneous porous media. Beginning with the Navier-Stokes equations, we find the volume averaged momentum equation to be given by $$\langle v_\beta \rangle = - \frac{K}{{\mu _\beta }} \cdot (\nabla \langle p_\beta \rangle ^\beta - \rho _\beta g) - F\cdot \langle v_\beta \rangle .$$ The Darcy's law permeability tensor, K, and the Forchheimer correction tensor, F, are determined by closure problems that must be solved using a spatially periodic model of a porous medium. When the Reynolds number is small compared to one, the closure problem can be used to prove that F is a linear function of the velocity, and order of magnitude analysis suggests that this linear dependence may persist for a wide range of Reynolds numbers.     

Two-phase flow in heterogeneous porous media III: Laboratory experiments for flow parallel to a stratified system

Two-phase flow in stratified porous media is a problem of central importance in the study of oil recovery processes. In general, these flows are parallel to the stratifications, and it is this type of flow that we have investigated experimentally and theoretically in this study. The experiments were performed with a two-layer model of a stratified porous medium. The individual strata were composed of Aerolith-10, an artificial: sintered porous medium, and Berea sandstone, a natural porous medium reputed to be relatively homogeneous. Waterflooding experiments were performed in which the saturation field was measured by gamma-ray absorption. Data were obtained at 150 points distributed evenly over a flow domain of 0.1 × 0.6 m. The slabs of Aerolith-10 and Berea sandstone were of equal thickness, i.e. 5 centimeters thick. An intensive experimental study was carried out in order to accurately characterize the individual strata; however, this effort was hampered by both local heterogeneities and large-scale heterogeneities.The theoretical analysis of the waterflooding experiments was based on the method of large-scale averaging and the large-scale closure problem. The latter provides a precise method of discussing the crossflow phenomena, and it illustrates exactly how the crossflow influences the theoretical prediction of the large-scale permeability tensor. The theoretical analysis was restricted to the quasi-static theory of Quintard and Whitaker (1988), however, the dynamic effects described in Part I (Quintard and Whitaker 1990a) are discussed in terms of their influence on the crossflow.Roman Letters A interfacial area between the -region and the -region contained within V, m² - a vector that maps onto , m - b vector that maps onto , m - b vector that maps onto , m - B second order tensor that maps onto , m² - C second order tensor that maps onto , m² - E energy of the gamma emitter, keV - f fractional flow of the -phase - g gravitational vector, m/s² - h characteristic length of the large-scale averaging volume, m - H height of the stratified porous medium , m - i unit base vector in the x-direction - K local volume-averaged single-phase permeability, m² - K - {K}, large-scale spatial deviation permeability - { K} large-scale volume-averaged single-phase permeability, m² - K ^* large-scale single-phase permeability, m² - K ^** equivalent large-scale single-phase permeability, m² - K local volume-averaged -phase permeability in the -region, m² - K local volume-averaged -phase permeability in the -region, m² - K - {K } , large-scale spatial deviation for the -phase permeability, m² - K ^* large-scale permeability for the -phase, m² - l thickness of the porous medium, m - l characteristic length for the -region, m - l characteristic length for the -region, m - L length of the experimental porous medium, m - characteristic length for large-scale averaged quantities, m - n outward unit normal vector for the -region - n outward unit normal vector for the -region - n unit normal vector pointing from the -region toward the -region (n = - n ) - N number of photons - p pressure in the -phase, N/m² - p ₀ reference pressure in the -phase, N/m² - local volume-averaged intrinsic phase average pressure in the -phase, N/m² - large-scale volume-averaged pressure of the -phase, N/m² - large-scale intrinsic phase average pressure in the capillary region of the -phase, N/m² - - , large-scale spatial deviation for the -phase pressure, N/m² - p_c , capillary pressure, N/m² - p _c capillary pressure in the -region, N/m² - p capillary pressure in the -region, N/m² - {p _c } ^c large-scale capillary pressure, N/m² - q -phase velocity at the entrance of the porous medium, m/s - q -phase velocity at the entrance of the porous medium, m/s - S_wi irreducible water saturation - S /, local volume-averaged saturation for the -phase - S _i initial saturation for the -phase - S _r residual saturation for the -phase - S ^* { }^*/}^*, large-scale average saturation for the -phase - S saturation for the -phase in the -region - S saturation for the -phase in the -region - t time, s - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the -phase, m/s - {v } large-scale averaged velocity for the -phase, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v local volume-averaged phase average velocity for the -phase in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - v -{v } , large-scale spatial deviation for the -phase velocity in the -region, m/s - V large-scale averaging volume, m³ - y position vector relative to the centroid of the large-scale averaging volume, m - {y}^c large-scale average of y over the capillary region, m Greek Letters local porosity - local porosity in the -region - local porosity in the -region - local volume fraction for the -phase - local volume fraction for the -phase in the -region - local volume fraction for the -phase in the -region - {}^* { }^*+{ }^*, large-scale spatial average volume fraction - { }^* large-scale spatial average volume fraction for the -phase - mass density of the -phase, kg/m³ - mass density of the -phase, kg/m³ - viscosity of the -phase, N s/m² - viscosity of the -phase, Ns/m² - V /V , volume fraction of the -region ( + =1) - V /V , volume fraction of the -region ( + =1) - attenuation coefficient to gamma-rays, m^-1 - -      

The Interaction between Quasilinear Elastodynamics and the Navier-Stokes Equations

The interaction between a viscous fluid and an elastic solid is modeled by a system of parabolic and hyperbolic equations, coupled to one another along the moving material interface through the continuity of the velocity and traction vectors. We prove the existence and uniqueness (locally in time) of strong solutions in Sobolev spaces for quasilinear elastodynamics coupled to the incompressible Navier-Stokes equations. Unlike our approach in [5] for the case of linear elastodynamics, we cannot employ a fixed-point argument on the nonlinear system itself, and are instead forced to regularize it by a particular parabolic artificial viscosity term. We proceed to show that with this specific regularization, we obtain a time interval of existence which is independent of the artificial viscosity; together with a priori estimates, we identify the global solution (in both phases), as well as the interface motion, as a weak limit in strong norms of our sequence of regularized problems.