scholarly journals Planform evolution in convection–an embedded centre manifold

Author(s):  
A. J. Roberts

AbstractThe new motion of embedding a centre manifold in some higher-dimensional manifold leads to a practical approach to the rational low-dimensional approximation of a wide class of dynamical systems; it also provides a simple geometric picture for these approximations. In particular, I consider the problem of finding an approximate, but accurate, description of the evolution of a two-dimensional planform of convection. Inspired by a simple example, the straightforward adiabatic iteration is proposed to estimate an embedding manifold and arguments are presented for its effectiveness. Upon applying the procedure to a model convective planform problem I find that the resulting approximations perform remarkably well–much better than the traditional Swift-Hohenberg approximation for planform evolution.

Author(s):  
S. M. Cox ◽  
A. J. Roberts

AbstractCentre manifolds arise in a rational approach to the problem of forming low-dimensional models of dynamical systems with many degrees of freedom. When a dynamical system with a centre manifold is subject to a small forcing, F, there are two effects: to displace the centre manifold; and to alter the evolution thereon. We propose a formal scheme for calculating the centre manifold of such a forced dynamical system. Our formalism permits the calculation of these effects, with errors of order |F|2. We find that the displacement of the manifold allows a reparameterisation of its description, and we describe two “natural” ways in which this can be carried out. We give three examples: an introductory example; a five-mode model of the atmosphere to display the quasi-geostrophic approximation; and the forced Kuramoto-Sivashinsky equation.


2003 ◽  
Vol 15 (8) ◽  
pp. 1715-1749 ◽  
Author(s):  
Blaise Agüera y Arcas ◽  
Adrienne L. Fairhall ◽  
William Bialek

A spiking neuron “computes” by transforming a complex dynamical input into a train of action potentials, or spikes. The computation performed by the neuron can be formulated as dimensional reduction, or feature detection, followed by a nonlinear decision function over the low-dimensional space. Generalizations of the reverse correlation technique with white noise input provide a numerical strategy for extracting the relevant low-dimensional features from experimental data, and information theory can be used to evaluate the quality of the low-dimensional approximation. We apply these methods to analyze the simplest biophysically realistic model neuron, the Hodgkin-Huxley (HH) model, using this system to illustrate the general methodological issues. We focus on the features in the stimulus that trigger a spike, explicitly eliminating the effects of interactions between spikes. One can approximate this triggering “feature space” as a two-dimensional linear subspace in the high-dimensional space of input histories, capturing in this way a substantial fraction of the mutual information between inputs and spike time. We find that an even better approximation, however, is to describe the relevant subspace as two dimensional but curved; in this way, we can capture 90% of the mutual information even at high time resolution. Our analysis provides a new understanding of the computational properties of the HH model. While it is common to approximate neural behavior as “integrate and fire,” the HH model is not an integrator nor is it well described by a single threshold.


2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Wenchong Tian ◽  
Hao Wu

Abstract Transfer operators such as Perron–Frobenius and Koopman operator play a key role in modeling and analysis of complex dynamical systems, which allow linear representations of nonlinear dynamics by transforming the original state variables to feature spaces. However, it remains challenging to identify the optimal low-dimensional feature mappings from data. The variational approach for Markov processes (VAMP) provides a comprehensive framework for the evaluation and optimization of feature mappings based on the variational estimation of modeling errors, but it still suffers from a flawed assumption on the transfer operator and therefore sometimes fails to capture the essential structure of system dynamics. In this paper, we develop a powerful alternative to VAMP, called kernel embedding based variational approach for dynamical systems (KVAD). By using the distance measure of functions in the kernel embedding space, KVAD effectively overcomes theoretical and practical limitations of VAMP. In addition, we develop a data-driven KVAD algorithm for seeking the ideal feature mapping within a subspace spanned by given basis functions, and numerical experiments show that the proposed algorithm can significantly improve the modeling accuracy compared to VAMP.


2013 ◽  
Vol 1 ◽  
Author(s):  
MARK BRAVERMAN ◽  
KONSTANTIN MAKARYCHEV ◽  
YURY MAKARYCHEV ◽  
ASSAF NAOR

AbstractThe (real) Grothendieck constant ${K}_{G} $ is the infimum over those $K\in (0, \infty )$ such that for every $m, n\in \mathbb{N} $ and every $m\times n$ real matrix $({a}_{ij} )$ we have $$\begin{eqnarray*}\displaystyle \max _{\{ x_{i}\} _{i= 1}^{m} , \{ {y}_{j} \} _{j= 1}^{n} \subseteq {S}^{n+ m- 1} }\sum _{i= 1}^{m} \sum _{j= 1}^{n} {a}_{ij} \langle {x}_{i} , {y}_{j} \rangle \leqslant K\max _{\{ \varepsilon _{i}\} _{i= 1}^{m} , \{ {\delta }_{j} \} _{j= 1}^{n} \subseteq \{ - 1, 1\} }\sum _{i= 1}^{m} \sum _{j= 1}^{n} {a}_{ij} {\varepsilon }_{i} {\delta }_{j} . &&\displaystyle\end{eqnarray*}$$ The classical Grothendieck inequality asserts the nonobvious fact that the above inequality does hold true for some $K\in (0, \infty )$ that is independent of $m, n$ and $({a}_{ij} )$. Since Grothendieck’s 1953 discovery of this powerful theorem, it has found numerous applications in a variety of areas, but, despite attracting a lot of attention, the exact value of the Grothendieck constant ${K}_{G} $ remains a mystery. The last progress on this problem was in 1977, when Krivine proved that ${K}_{G} \leqslant \pi / 2\log (1+ \sqrt{2} )$ and conjectured that his bound is optimal. Krivine’s conjecture has been restated repeatedly since 1977, resulting in focusing the subsequent research on the search for examples of matrices $({a}_{ij} )$ which exhibit (asymptotically, as $m, n\rightarrow \infty $) a lower bound on ${K}_{G} $ that matches Krivine’s bound. Here, we obtain an improved Grothendieck inequality that holds for all matrices $({a}_{ij} )$ and yields a bound ${K}_{G} \lt \pi / 2\log (1+ \sqrt{2} )- {\varepsilon }_{0} $ for some effective constant ${\varepsilon }_{0} \gt 0$. Other than disproving Krivine’s conjecture, and along the way also disproving an intermediate conjecture of König that was made in 2000 as a step towards Krivine’s conjecture, our main contribution is conceptual: despite dealing with a binary rounding problem, random two-dimensional projections, when combined with a careful partition of ${ \mathbb{R} }^{2} $ in order to round the projected vectors to values in $\{ - 1, 1\} $, perform better than the ubiquitous random hyperplane technique. By establishing the usefulness of higher-dimensional rounding schemes, this fact has consequences in approximation algorithms. Specifically, it yields the best known polynomial-time approximation algorithm for the Frieze–Kannan Cut Norm problem, a generic and well-studied optimization problem with many applications.


2022 ◽  
Vol 0 (0) ◽  
pp. 0
Author(s):  
Héctor Barge ◽  
José M. R. Sanjurjo

<p style='text-indent:20px;'>In this paper we study generalized Poincaré-Andronov-Hopf bifurcations of discrete dynamical systems. We prove a general result for attractors in <inline-formula><tex-math id="M1">\begin{document}$ n $\end{document}</tex-math></inline-formula>-dimensional manifolds satisfying some suitable conditions. This result allows us to obtain sharper Hopf bifurcation theorems for fixed points in the general case and other attractors in low dimensional manifolds. Topological techniques based on the notion of concentricity of manifolds play a substantial role in the paper.</p>


Author(s):  
H.A. Cohen ◽  
W. Chiu ◽  
J. Hosoda

GP 32 (molecular weight 35000) is a T4 bacteriophage protein that destabilizes the DNA helix. The fragment GP32*I (77% of the total weight), which destabilizes helices better than does the parent molecule, crystallizes as platelets thin enough for electron diffraction and electron imaging. In this paper we discuss the structure of this protein as revealed in images reconstructed from stained and unstained crystals.Crystals were prepared as previously described. Crystals for electron microscopy were pelleted from the buffer suspension, washed in distilled water, and resuspended in 1% glucose. Two lambda droplets were placed on grids over freshly evaporated carbon, allowed to sit for five minutes, and then were drained. Stained crystals were prepared the same way, except that prior to draining the droplet, two lambda of aqueous 1% uranyl acetate solution were applied for 20 seconds. Micrographs were produced using less than 2 e/Å2 for unstained crystals or less than 8 e/Å2 for stained crystals.


Doklady BGUIR ◽  
2020 ◽  
Vol 18 (7) ◽  
pp. 87-95
Author(s):  
M. S. Baranava ◽  
P. A. Praskurava

The search for fundamental physical laws which lead to stable high-temperature ferromagnetism is an urgent task. In addition to the already synthesized two-dimensional materials, there remains a wide list of possible structures, the stability of which is predicted theoretically. The article suggests the results of studying the electronic properties of MAX3 (M = Cr, Fe, A = Ge, Si, X = S, Se, Te) transition metals based compounds with nanostructured magnetism. The research was carried out using quantum mechanical simulation in specialized VASP software and calculations within the Heisenberg model. The ground magnetic states of twodimensional MAX3 and the corresponding energy band structures are determined. We found that among the systems under study, CrGeTe3 is a semiconductor nanosized ferromagnet. In addition, one is a semiconductor with a bandgap of 0.35 eV. Other materials are antiferromagnetic. The magnetic moment in MAX3 is localized on the transition metal atoms: in particular, the main one on the d-orbital of the transition metal atom (and only a small part on the p-orbital of the chalcogen). For CrGeTe3, the exchange interaction integral is calculated. The mechanisms of the formation of magnetic order was established. According to the obtained exchange interaction integrals, a strong ferromagnetic order is formed in the semiconductor plane. The distribution of the projection density of electronic states indicates hybridization between the d-orbital of the transition metal atom and the p-orbital of the chalcogen. The study revealed that the exchange interaction by the mechanism of superexchange is more probabilistic.


2007 ◽  
Vol 5 ◽  
pp. 195-200
Author(s):  
A.V. Zhiber ◽  
O.S. Kostrigina

In the paper it is shown that the two-dimensional dynamical system of equations is Darboux integrable if and only if its characteristic Lie algebra is finite-dimensional. The class of systems having a full set of fist and second order integrals is described.


2004 ◽  
Vol 126 (5) ◽  
pp. 861-870 ◽  
Author(s):  
A. Thakur ◽  
X. Liu ◽  
J. S. Marshall

An experimental and computational study is performed of the wake flow behind a single yawed cylinder and a pair of parallel yawed cylinders placed in tandem. The experiments are performed for a yawed cylinder and a pair of yawed cylinders towed in a tank. Laser-induced fluorescence is used for flow visualization and particle-image velocimetry is used for quantitative velocity and vorticity measurement. Computations are performed using a second-order accurate block-structured finite-volume method with periodic boundary conditions along the cylinder axis. Results are applied to assess the applicability of a quasi-two-dimensional approximation, which assumes that the flow field is the same for any slice of the flow over the cylinder cross section. For a single cylinder, it is found that the cylinder wake vortices approach a quasi-two-dimensional state away from the cylinder upstream end for all cases examined (in which the cylinder yaw angle covers the range 0⩽ϕ⩽60°). Within the upstream region, the vortex orientation is found to be influenced by the tank side-wall boundary condition relative to the cylinder. For the case of two parallel yawed cylinders, vortices shed from the upstream cylinder are found to remain nearly quasi-two-dimensional as they are advected back and reach within about a cylinder diameter from the face of the downstream cylinder. As the vortices advect closer to the cylinder, the vortex cores become highly deformed and wrap around the downstream cylinder face. Three-dimensional perturbations of the upstream vortices are amplified as the vortices impact upon the downstream cylinder, such that during the final stages of vortex impact the quasi-two-dimensional nature of the flow breaks down and the vorticity field for the impacting vortices acquire significant three-dimensional perturbations. Quasi-two-dimensional and fully three-dimensional computational results are compared to assess the accuracy of the quasi-two-dimensional approximation in prediction of drag and lift coefficients of the cylinders.


2021 ◽  
Vol 183 (2) ◽  
Author(s):  
D. Benedetto ◽  
E. Caglioti ◽  
S. Caracciolo ◽  
M. D’Achille ◽  
G. Sicuro ◽  
...  

AbstractWe consider the assignment problem between two sets of N random points on a smooth, two-dimensional manifold $$\Omega $$ Ω of unit area. It is known that the average cost scales as $$E_{\Omega }(N)\sim {1}/{2\pi }\ln N$$ E Ω ( N ) ∼ 1 / 2 π ln N with a correction that is at most of order $$\sqrt{\ln N\ln \ln N}$$ ln N ln ln N . In this paper, we show that, within the linearization approximation of the field-theoretical formulation of the problem, the first $$\Omega $$ Ω -dependent correction is on the constant term, and can be exactly computed from the spectrum of the Laplace–Beltrami operator on $$\Omega $$ Ω . We perform the explicit calculation of this constant for various families of surfaces, and compare our predictions with extensive numerics.


Sign in / Sign up

Export Citation Format

Share Document