Non-standard coupled extensional and bending bias tests for planar pantographic lattices. Part I: numerical simulations

In dell’Isola et al. (Zeitschrift für Angewandte Math und Physik 66(6):3473–3498, 2015, Proc R Soc Lond A Math Phys Eng Sci 472(2185), 2016), the concept of pantographic sheet is proposed. The aim is to design a metamaterial showing: (i) a large range of elastic response; (ii) an extreme toughness in extensional deformation; (iii) a convenient ratio between toughness and weight. However, these required properties must coexist with non-detrimental mechanical characteristics in the presence of other kinds of imposed displacements. The aim of this paper is to prove via numerical simulations that pantographic sheets may effectively resist to coupled bending and extensional deformations. The four-parameter model introduced shows its versatility as it is able to encompass all the considered types of (large) deformations. The numerical integration scheme which we use is based on the same concepts exploited in Turco et al. (Zeitschrift für Angewandte Math und Physik 67(4):1–28, 2016): They prove that the Hencky-type discretization is very efficient also in nonlinear large deformations and large displacements regimes. In Part II of this paper, we will show that the used models are very effective to describe experimental evidence.     

Non-standard coupled extensional and bending bias tests for planar pantographic lattices. Part II: comparison with experimental evidence

In dell’Isola et al. (Zeitschrift für Angewandte Math und Physik 66(6):3473–3498, 2015, Proc R Soc Lond A Math Phys Eng Sci 472(2185):1–23, 2016) pantographic sheets are proposed as a basic constituent for a novel metamaterial. In Part I, see Turco et al. (Zeitschrift für Angewandte Math und Physik, doi:10.1007/s00033-016-0713-4, 2016), two different numerical models are applied in order to design an experimental setup aimed to prove the effectiveness of introduced concept. The aim of this paper is to prove that the Hencky-type model introduced for planar pantographic sheets allows for the correct prediction, in a large range of imposed displacements, of the experimental measurements concerning specimens undergoing coupled bending and extensional deformations. The four-parameter numerical model introduced is shown to have a large range of applicability: Indeed without changing the values of the material parameters previously attributed in simple extensional tests to a specific specimen by a best-fit procedure, it is possible to forecast its behavior in all the considered type of imposed deformations. The measurements performed include the determination of reactive forces exerted by used hard devices, and the numerical modeling is able to predict very carefully quantitatively and qualitatively also this complex aspect of phenomenology, where previously attempted models seem to have failed.     

Numerical identification procedure between a micro-Cauchy model and a macro-second gradient model for planar pantographic structures

In order to design the microstructure of metamaterials showing high toughness in extension (property to be shared with muscles), it has been recently proposed (Dell’Isola et al. in Z Angew Math Phys 66(6):3473–3498, 2015) to consider pantographic structures. It is possible to model such structures at a suitably small length scale (resolving in detail the interconnecting pivots/cylinders) using a standard Cauchy first gradient theory. However, the computational costs for such modelling choice are not allowing for the study of more complex mechanical systems including for instance many pantographic substructures. The microscopic model considered here is a quadratic isotropic Saint-Venant first gradient continuum including geometric nonlinearities and characterized by two Lamé parameters. The introduced macroscopic two-dimensional model for pantographic sheets is characterized by a deformation energy quadratic both in the first and second gradient of placement. However, as underlined in Dell’Isola et al. (Proc R Soc Lond A 472(2185):20150790, 2016), it is needed that the second gradient stiffness depends on the first gradient of placement if large deformations and large displacements configurations must be described. The numerical identification procedure presented in this paper consists in fitting the macro-constitutive parameters using several numerical simulations performed with the micro-model. The parameters obtained by the best fit identification in few deformation problems fit very well also in many others, showing that the reduced proposed model is suitable to get an effective model at relevantly lower computational effort. The presented numerical evidences suggest that a rigorous mathematical homogenization result most likely holds.     

Finite-Codimensional Compressions of Symmetric and Self-Adjoint Linear Relations in Krein Spaces

Theorems due to Stenger (Bull Am Math Soc 74:369–372, 1968) and Nudelman (Int Equ Oper Theory 70:301–305, 2011) in Hilbert spaces and their generalizations to Krein spaces in Azizov and Dijksma (Int Equ Oper Theory 74(2):259–269, 2012) and Azizov et al. (Linear Algebra Appl 439:771–792, 2013) generate additional questions about properties a finite-codimensional compression \({T_0}\) of a symmetric or self-adjoint linear relation \({T}\) may or may not inherit from \({T}\). These questions concern existence of invariant maximal nonnegative subspaces, definitizability, singular critical points and defect indices.     

Agler-Commutant Lifting on an Annulus

This note presents a commutant lifting theorem (CLT) of Agler type for the annulus \({\mathbb A}\) . Here the relevant set of test functions are the minimal inner functions on \({\mathbb A}\) —those analytic functions on \({\mathbb A}\) which are unimodular on the boundary and have exactly two zeros in \({\mathbb A}\) —and the model space is determined by a distinguished member of the Sarason family of kernels over \({\mathbb A}\) . The ideas and constructions borrow freely from the CLT of Ball et al. (Indiana Univ Math J 48(2):653–675, 1999) and Archer (Unitary dilations of commuting contractions. PhD thesis, University of Newcastle, 2004) for the polydisc, and Ambrozie and Eschmeier (A commutant lifting theorem on analytic polyhedra. Topological algebras, their applications, and related topics, 83108, Banach Center Publications, vol 67. Polish Academy of Sciences, Warsaw, 2005) for the ball in \({\mathbb C^n}\) , as well as generalizations of the de Branges–Rovnyak construction like found in Agler (On the representation of certain holomorphic functions defined on a polydisc. Topics in operator theory: Ernst D. Hellinger memorial volume, operator theory: advances and applications, vol 48. Birkhäuser, Basel, pp 47–66, 1990) and Ambrozie et al. (J Oper Theory 47(2):287–302, 2002). It offers a template for extending the result in McCullough and Sultanic (Complex Anal Oper Theory 1(4):581–620, 2007) to infinitely many test functions. Among the needed new ingredients is the formulation of the factorization implicit in the statement of the results in Ball et al. (Indiana Univ Math J 48(2):653–675, 1999) and Archer (Unitary dilations of commuting contractions. PhD thesis, University of Newcastle, 2004) and McCullough and Sultanic (Complex Anal Oper Theory 1(4):581–620, 2007) in terms of certain functional Hilbert spaces of Hilbert space valued functions.     

Vanishing of Rabinowitz Floer homology on negative line bundles

Following Frauenfelder (Rabinowitz action functional on very negative line bundles, Habilitationsschrift, Munich/München, 2008), Albers and Frauenfelder (Bubbles and onis, 2014. arXiv:1412.4360) we construct Rabinowitz Floer homology for negative line bundles over symplectic manifolds and prove a vanishing result. Ritter (Adv Math 262:1035–1106, 2014) showed that symplectic homology of these spaces does not vanish, in general. Thus, the theorem \(\mathrm {SH}=0\Leftrightarrow \mathrm {RFH}=0\) (Ritter in J Topol 6(2):391–489, 2013), does not extend beyond the symplectically aspherical situation. We give a conjectural explanation in terms of the Cieliebak–Frauenfelder–Oancea long exact sequence Cieliebak et al. (Ann Sci Éc Norm Supér (4) 43(6):957–1015, 2010).     

Minimising movements for the motion of discrete screw dislocations along glide directions

In Alicandro et al. (J Mech Phys Solids 92:87–104, 2016) a simple discrete scheme for the motion of screw dislocations toward low energy configurations has been proposed. There, a formal limit of such a scheme, as the lattice spacing and the time step tend to zero, has been described. The limiting dynamics agrees with the maximal dissipation criterion introduced in Cermelli and Gurtin (Arch Ration Mech Anal 148, 1999) and predicts motion along the glide directions of the crystal. In this paper, we provide rigorous proofs of the results in [3], and in particular of the passage from the discrete to the continuous dynamics. The proofs are based on \(\Gamma \)-convergence techniques.     

Hedging with temporary price impact

We consider the problem of hedging a European contingent claim in a Bachelier model with temporary price impact as proposed by Almgren and Chriss (J Risk 3:5–39, 2001). Following the approach of Rogers and Singh (Math Financ 20:597–615, 2010) and Naujokat and Westray (Math Financ Econ 4(4):299–335, 2011), the hedging problem can be regarded as a cost optimal tracking problem of the frictionless hedging strategy. We solve this problem explicitly for general predictable target hedging strategies. It turns out that, rather than towards the current target position, the optimal policy trades towards a weighted average of expected future target positions. This generalizes an observation of Gârleanu and Pedersen (Dynamic portfolio choice with frictions. Preprint, 2013b) from their homogenous Markovian optimal investment problem to a general hedging problem. Our findings complement a number of previous studies in the literature on optimal strategies in illiquid markets as, e.g., Gârleanu and Pedersen (Dynamic portfolio choice with frictions. Preprint, 2013b), Naujokat and Westray (Math Financ Econ 4(4):299–335, 2011), Rogers and Singh (Math Financ 20:597–615, 2010), Almgren and Li (Option hedging with smooth market impact. Preprint, 2015), Moreau et al. (Math Financ. doi: 10.1111/mafi.12098, 2015), Kallsen and Muhle-Karbe (High-resilience limits of block-shaped order books. Preprint, 2014), Guasoni and Weber (Mathematical Financ. doi: 10.1111/mafi.12099, 2015a; Nonlinear price impact and portfolio choice. Preprint, 2015b), where the frictionless hedging strategy is confined to diffusions. The consideration of general predictable reference strategies is made possible by the use of a convex analysis approach instead of the more common dynamic programming methods.     

Reformulation of Hensel’s Lemma and extension of a theorem of Ore

In this paper, we extend the theorem of Ore regarding factorization of polynomials over p-adic numbers to henselian valued fields of arbitrary rank thereby generalizing the main results of Khanduja and Kumar (J Pure Appl Algebra 216:2648–2656, 2012) and Cohen et al. (Mathematika 47:173–196, 2000). As an application, we derive the analogue of Dedekind’s Theorem regarding splitting of rational primes in algebraic number fields as well as of its converse for general valued fields extending similar results proved for discrete valued fields in Khanduja and Kumar (Int J Number Theory 4:1019–1025, 2008). The generalized version of Ore’s Theorem leads to an extension of a result of Weintraub dealing with a generalization of Eisenstein Irreducibility Criterion (cf. Weintraub in Proc Am Math Soc 141:1159–1160, 2013). We also give a reformulation of Hensel’s Lemma for polynomials with coefficients in henselian valued fields which is used in the proof of the extended Ore’s Theorem and was proved in Khanduja and Kumar (J Algebra Appl 12:1250125, 2013) in the particular case of complete rank one valued fields.     

Non-Stationary Abstract Friedrichs Systems

Inspired by the results of Ern et al. (Commun Partial Differ Equ 32:317–341, 2007) on the abstract theory for Friedrichs symmetric positive systems, we give the existence and uniqueness result for the initial- (boundary) value problem for the non-stationary abstract Friedrichs system. Despite the absence of the well-posedness result for such systems, there were already attempts for their numerical treatment by Burman et al. (SIAM J Numer Anal 48:2019–2042, 2010) and Bui-Thanh et al. (SIAM J Numer Anal 51:1933–1958, 2013). We use the semigroup theory approach and prove that the operator involved satisfies the conditions of the Hille–Yosida generation theorem. We also address the semilinear problem and apply the new results to a number of examples, such as the symmetric hyperbolic system, the unsteady div–grad problem, and the wave equation. Special attention was paid to the (generalised) unsteady Maxwell system.     

New characterizations of the Takagi function via functional equations

We provide two new characterizations of the Takagi function as the unique bounded solution of some systems of two functional equations. The results are independent of those obtained by Kairies (Wy? Szko? Ped Krakow Rocznik Nauk Dydakt Prace Mat 196:73–82, 1998), Kairies (Aequ Math 53:207–241, 1997), Kairies (Aequ Math 58:183–191, 1999) and Kairies et al. (Rad Mat 4:361–374, 1989; Errata, Rad Mat 5:179–180, 1989).     

A remark on energy estimates concerning extremals for Trudinger–Moser inequalities on a disc

In this short note, we generalized an energy estimate due to Malchiodi–Martinazzi (J Eur Math Soc 16:893–908, 2014) and Mancini–Martinazzi (Calc Var 56:94, 2017). As an application, we used it to reprove existence of extremals for Trudinger–Moser inequalities of Adimurthi–Druet type on the unit disc. Such existence problems in general cases had been considered by Yang  (Trans Am Math Soc 359:5761–5776, 2007; J Differ Equ 258:3161–3193, 2015) and Lu–Yang (Discrete Contin Dyn Syst 25:963–979, 2009) by using another method.     

Absolute Cesàro Series Spaces and Matrix Operators

In this paper we derive a series space \(\vert C_{\lambda,\mu} \vert _{k}\) using the well known absolute Cesàro summability \(\vert C_{\lambda,\mu} \vert _{k}\) of Das (Proc. Camb. Philol. Soc. 67:321–326, 1970), compute its \(\beta\)-dual, give some algebraic and topological properties, and characterize some matrix operators defined on that space. So we generalize some results of Bosanquet (J. Lond. Math. Soc. 20:39–48, 1945), Flett (Proc. Lond. Math. Soc. 7:113–141, 1957), Mehdi (Proc. Lond. Math. Soc. (3)10:180–199, 1960), Mazhar (Tohoku Math. J. 23:433–451, 1971), Orhan and Sar?göl (Rocky Mt. J. Math. 23(3):1091–1097, 1993) and Sar?göl (Commun. Math. Appl. 7(1):11–22, 2016; Math. Comput. Model. 55:1763–1769, 2012).     

On a rigidity result of singularities

The aim of this note is to prove, in the spirit of a rigidity result for isolated singularities of Schlessinger see Schlessinger (Invent Math 14:17–26, 1971) or also Kleiman and Landolfi (Compositio Math 23:407–434, 1971), a variant of a rigidity criterion for arbitrary singularities (Theorem 2.1 below). The proof of this result does not use Schlessinger’s Deformation Theory [Schlessinger (Trans Am Math Soc 130:208–222, 1968) and Schlessinger (Invent Math 14:17–26, 1971)]. Instead it makes use of Local Grothendieck-Lefschetz Theory, see (Grothendieck 1968, Éxposé 9, Proposition 1.4, page 106) and a Lemma of Zariski, see (Zariski, Am J Math 87:507–536, 1965, Lemma 4, page 526). I hope that this proof, although works only in characteristic zero, might also have some interest in its own.     

On the maximum TSP with <Emphasis Type="Italic">γ</Emphasis>-parameterized triangle inequality

The maximum TSP with γ-parameterized triangle inequality is defined as follows. Given a complete graph G = (V, E, w) in which the edge weights satisfy w(uv) ≤ γ · (w(ux) + w(xv)) for all distinct nodes \({u,x,v \in V}\), find a tour with maximum weight that visits each node exactly once. Recently, Zhang et al. (Theor Comput Sci 411(26–28):2537–2541, 2010) proposed a \({\frac{\gamma+1}{3\gamma}}\)-approximation algorithm for \({\gamma\in\left[\frac{1}{2},1\right)}\). In this paper, we show that the approximation ratio of Kostochka and Serdyukov’s algorithm (Upravlyaemye Sistemy 26:55–59, 1985) is \({\frac{4\gamma+1}{6\gamma}}\), and the expected approximation ratio of Hassin and Rubinstein’s randomized algorithm (Inf Process Lett 81(5):247–251, 2002) is \({\frac{3\gamma+\frac{1}{2}}{4\gamma}-O\left(\frac{1}{\sqrt{n}}\right)}\), for \({\gamma\in\left[\frac{1}{2},+\infty\right)}\). These improve the result in Zhang et al. (Theor Comput Sci 411(26–28):2537–2541, 2010) and generalize the results in Hassin and Rubinstein and Kostochka and Serdyukov (Inf Process Lett 81(5):247–251, 2002; Upravlyaemye Sistemy 26:55–59, 1985).     

Data dependence for the amplitude equation of surface waves

We consider the amplitude equation for nonlinear surface wave solutions of hyperbolic conservation laws. This is an asymptotic nonlocal, Hamiltonian evolution equation with quadratic nonlinearity. For example, this equation describes the propagation of nonlinear Rayleigh waves (Hamilton et al. in J Acoust Soc Am 97:891–897, 1995), surface waves on current-vortex sheets in incompressible MHD (Alì and Hunter in Q Appl Math 61(3):451–474, 2003; Alì et al. in Stud Appl Math 108(3):305–321, 2002) and on the incompressible plasma–vacuum interface (Secchi in Q Appl Math 73(4):711–737, 2015). The local-in-time existence of smooth solutions to the Cauchy problem for the amplitude equation in noncanonical variables was shown in Hunter (J Hyperbolic Differ Equ 3(2):247–267, 2006), Secchi (Q Appl Math 73(4):711–737, 2015). In the present paper we prove the continuous dependence in strong norm of solutions on the initial data. This completes the proof of the well-posedness of the problem in the classical sense of Hadamard.     

Improved local convergence analysis of the Gauss–Newton method under a majorant condition

We present a local convergence analysis of Gauss-Newton method for solving nonlinear least square problems. Using more precise majorant conditions than in earlier studies such as Chen (Comput Optim Appl 40:97–118, 2008), Chen and Li (Appl Math Comput 170:686–705, 2005), Chen and Li (Appl Math Comput 324:1381–1394, 2006), Ferreira (J Comput Appl Math 235:1515–1522, 2011), Ferreira and Gonçalves (Comput Optim Appl 48:1–21, 2011), Ferreira and Gonçalves (J Complex 27(1):111–125, 2011), Li et al. (J Complex 26:268–295, 2010), Li et al. (Comput Optim Appl 47:1057–1067, 2004), Proinov (J Complex 25:38–62, 2009), Ewing, Gross, Martin (eds.) (The merging of disciplines: new directions in pure, applied and computational mathematics 185–196, 1986), Traup (Iterative methods for the solution of equations, 1964), Wang (J Numer Anal 20:123–134, 2000), we provide a larger radius of convergence; tighter error estimates on the distances involved and a clearer relationship between the majorant function and the associated least squares problem. Moreover, these advantages are obtained under the same computational cost.     

Stability of Blowup for a 1D Model of Axisymmetric 3D Euler Equation

The question of the global regularity versus finite- time blowup in solutions of the 3D incompressible Euler equation is a major open problem of modern applied analysis. In this paper, we study a class of one-dimensional models of the axisymmetric hyperbolic boundary blow-up scenario for the 3D Euler equation proposed by Hou and Luo (Multiscale Model Simul 12:1722–1776, 2014) based on extensive numerical simulations. These models generalize the 1D Hou–Luo model suggested in Hou and Luo Luo and Hou (2014), for which finite-time blowup has been established in Choi et al. (arXiv preprint. arXiv:1407.4776, 2014). The main new aspects of this work are twofold. First, we establish finite-time blowup for a model that is a closer approximation of the three-dimensional case than the original Hou–Luo model, in the sense that it contains relevant lower-order terms in the Biot–Savart law that have been discarded in Hou and Luo Choi et al. (2014). Secondly, we show that the blow-up mechanism is quite robust, by considering a broader family of models with the same main term as in the Hou–Luo model. Such blow-up stability result may be useful in further work on understanding the 3D hyperbolic blow-up scenario.     

Global Existence and Asymptotic Behavior of Solutions to the Hyperbolic Keller-Segel Equation with a Logistic Source

In this paper we consider a hyperbolic Keller-Segel system with a logistic source in two dimension. We show the system has a global smooth solution upon small perturbation around a constant equilibrium and the solution satisfies a dissipative energy inequality. To do this we find a convex entropy functional and a compensating matrix, which transforms the partially dissipative system into a uniformly dissipative one. Those two ingredients were crucial for the study of a partially dissipative hyperbolic system (Hanouzet and Natalini in Arch. Ration. Mech. Anal. 169(2):89–117, 2003; Kawashima in Ph.D. Thesis, Kyoto University, 1983; Yong in Arch. Ration. Mech. Anal. 172(2):247–266, 2004).     

Variational Integrators for Interconnected Lagrange–Dirac Systems

Interconnected systems are an important class of mathematical models, as they allow for the construction of complex, hierarchical, multiphysics, and multiscale models by the interconnection of simpler subsystems. Lagrange–Dirac mechanical systems provide a broad category of mathematical models that are closed under interconnection, and in this paper, we develop a framework for the interconnection of discrete Lagrange–Dirac mechanical systems, with a view toward constructing geometric structure-preserving discretizations of interconnected systems. This work builds on previous work on the interconnection of continuous Lagrange–Dirac systems (Jacobs and Yoshimura in J Geom Mech 6(1):67–98, 2014) and discrete Dirac variational integrators (Leok and Ohsawa in Found Comput Math 11(5), 529–562, 2011). We test our results by simulating some of the continuous examples given in Jacobs and Yoshimura (2014).