Search

Based on the control strategies, transmission mechanism and clinical progression of COVID-19, we propose a compartmental model with time delays. Two time delays are introduced into the model to describe the incubation period of the disease and the quarantine period of uninfected individuals who have contacts with infected people. In order to reveal the spread rule for COVID-19, we study the threshold dynamics for the model. The basic reproduction number \begin{document}$ \mathcal R_0 $\end{document} is obtained. When \begin{document}$ \mathcal R_0 , the disease-free equilibrium is locally asymptotically stable, when \begin{document}$ \mathcal R_0>1 $\end{document} , the disease-free equilibrium is unstable and the disease is uniformly persistent. As the model's applications, we study COVID-19 transmission in the United States. The parameters are chosen to fit public data in the US. The numerical results indicate that an outbreak peak time in the US will appear in the middle of March. Sensitive analysis results show that enhancing the control measures, such as keeping social distance, wearing masks, isolation etc., can significantly contribute to the prevention and control of COVID-19 infection.

This paper applies a powerful scheme, namely Bernoulli sub-equation function method, to some partial differential equations with high non-linearity. Many new travelling wave solutions, such as mixed dark-bright soliton, exponential and complex domain, are reported. Under a suitable choice of the values of parameters, wave behaviours of the results obtained in the paper – in terms of 2D, 3D and contour surfaces – are observed.

In this paper, we study the null controllability of a nonlinear age, space and two-sex structured population dynamics model. This model is such that the nonlinearity and the couplage are at birth level. We consider a population with males and females and we are dealing with two cases of null controllability problems. The first problem is related to the total extinction, which means that, we estimate a time \begin{document}$ T $\end{document} to bring the male and female subpopulation density to zero. The second case concerns null controllability of male or female subpopulation. Since the absence of males or females in the population stops births; so, if we have the total extinction of the females at time \begin{document}$ T, $\end{document} and if \begin{document}$ A $\end{document} is the life span of the individuals, at time \begin{document}$ T+A $\end{document} one will get certainly the total extinction of the population. Our method uses first an observability inequality related to the adjoint of an auxiliary system, a null controllability of the linear auxiliary system, and after the Schauder's fixed point theorem.

In this paper, we propose an efficient numerical method for a comprehensive infection model that is formulated by a system of nonlinear coupling advection-diffusion-reaction equations. Using some subtle mixed explicit-implicit treatments, we construct a linearized and decoupled discrete scheme. Moreover, the proposed scheme is capable of preserving the positivity of variables, which is an essential requirement of the model under consideration. The proposed scheme uses the cell-centered finite difference method for the spatial discretization, and thus, it is easy to implement. The diffusion terms are treated implicitly to improve the robustness of the scheme. A semi-implicit upwind approach is proposed to discretize the advection terms, and a distinctive feature of the resulting scheme is to preserve the positivity of variables without any restriction on the spatial mesh size and time step size. We rigorously prove the unique existence of discrete solutions and positivity-preserving property of the proposed scheme without requirements for the mesh size and time step size. It is worthwhile to note that these properties are proved using the discrete variational principles rather than the conventional approaches of matrix analysis. Numerical results are also provided to assess the performance of the proposed scheme.

Classical bearing capacity theory was developed mainly based on spatially uniform soil properties, which cannot account for the influence of inherent soil variability. If the soil strength is heterogeneous, then using the average strength may overestimate the bearing capacity of foundations, because the failure mechanism may preferentially mobilise the weaker soils. In this study the aim is to establish a theoretical model using upper-bound solutions applied to the bearing capacity analysis of shallow foundations on undrained clay considering spatial variability. The model is derived on the principle of least energy dissipation using a four-parameter variation on Prandtl's mechanism. The theoretical model developed is verified by the random finite-element (FE) method in spatially varying soil conditions. The results show that the model can accurately capture the effect of spatially varying strength on the shallow foundation failure mechanism. The difference of bearing capacity factor between the proposed model and the FE model is within 5%, which demonstrates that the four-parameter model has an accuracy that is comparable to FE analysis with many hundreds of degrees of freedom. Another advantage of the theoretical model is that the possible non-convergence in FE analysis can be avoided, and hence, the calculation efficiency is significantly enhanced. The model is therefore suitable for rapid quantification of bearing capacity in spatially varying soils.

Hyperbolic systems of conservation laws in multiple spatial dimensions display features absent in the one-dimensional case, such as involutions and non-trivial stationary states. These features need to be captured by numerical methods without excessive grid refinement. The active flux method is an extension of the finite volume scheme with additional point values distributed along the cell boundary. For the equations of linear acoustics, an exact evolution operator can be used for the update of these point values. It incorporates all multi-dimensional information. The active flux method is stationarity preserving, i.e., it discretizes all the stationary states of the PDE. This paper demonstrates the experimental evidence for the discrete stationary states of the active flux method and shows the evolution of setups towards a discrete stationary state.

We show that the solution (in the sense of distribution) to the Cauchy problem with the periodic boundary condition associated with the modified Benjamin–Ono equation is unique in $L^\infty _t(H^s(\mathbb{T} ))$ for $s>1/2$. The proof is based on the analysis of a normal form equation obtained by infinitely many reduction steps using integration by parts in time after a suitable gauge transform.

The Rosenau–Korteweg-deVries–Kawahara equation describes the dynamics of dense discrete systems or small-amplitude gravity capillary waves on water of a finite depth. In this paper, we prove the well-posedness of the classical solutions for the Cauchy problem.

We discuss the asymptotic stability of stationary solutions to the incompressible Navier–Stokes equations on the whole space in Besov spaces. A critical estimate for the semigroup generated by the Laplacian with a perturbation is the main ingredient of the argument.

We study the family of blowup solutions to semilinear elliptic equations in two-space dimensions with exponentially-dominated nonnegative nonlinearities. Such a family admits an exclusion of the boundary blowup, finiteness of blowup points, and pattern formation. Then, Hamiltonian control of the location of blowup points, residual vanishing, and mass quantization arise under the estimate from below of the nonlinearity. Finally, if the principal growth rate of nonlinearity is exactly exponential and the residual part has a gap relative to this term, there is a locally uniform estimate of the solution which ensures its asymptotic non-degeneracy.

This article concerns the multiplicity of solutions for a class of nonlocal and nonhomogeneous Neumann boundary value problems involving the p ( x )-Laplacian, in which both nonlinear terms assume critical growth. We use variational method, exploring an important truncation argument and properties of the genus.

We design a virtual element method for the numerical treatment of the two-dimensional parabolic variational inequality problem on unstructured polygonal meshes. Due to the expected low regularity of the exact solution, the virtual element method is based on the lowest-order virtual element space that contains the subspace of the linear polynomials defined on each element. The connection between the nonnegativity of the virtual element functions and the nonnegativity of the degrees of freedom, i.e., the values at the mesh vertices, is established by applying the Maximum and Minimum Principle Theorem. The mass matrix is computed through an approximate L 2 polynomial projection, whose properties are carefully investigated in the paper. We prove the well-posedness of the resulting scheme in two different ways that reveal the contractive nature of the VEM and its connection with the minimization of quadratic functionals. The convergence analysis requires the existence of a nonnegative quasi-interpolation operator, whose construction is also discussed in the paper. The variational crime introduced by the virtual element setting produces five error terms that we control by estimating a suitable upper bound. Numerical experiments confirm the theoretical convergence rate for the refinement in space and time on three different mesh families including distorted squares, nonconvex elements, and Voronoi tesselations.

The understanding of nonlinear, high dimensional flows, e.g, atmospheric and ocean flows, is critical to address the impacts of global climate change. Data Assimilation techniques combine physical models and observational data, often in a Bayesian framework, to predict the future state of the model and the uncertainty in this prediction. Inherent in these systems are noise (Gaussian and non-Gaussian), nonlinearity, and high dimensionality that pose challenges to making accurate predictions. To address these issues we investigate the use of both model and data dimension reduction based on techniques including Assimilation in Unstable Subspaces, Proper Orthogonal Decomposition, and Dynamic Mode Decomposition. Algorithms that take advantage of projected physical and data models may be combined with Data Analysis techniques such as Ensemble Kalman Filter and Particle Filter variants. The projected Data Assimilation techniques are developed for the optimal proposal particle filter and applied to the Lorenz'96 and Shallow Water Equations to test the efficacy of our techniques in high dimensional, nonlinear systems., 30 pages, 13 figures, 3 tables To appear in Computers & Mathematics with Applications, 2021,ISSN 0898-1221

An efficient strategy for the numerical solution of time-fractional diffusion-reaction problems is devised. A standard finite difference discretization of the space derivative is initially applied which results in a linear stiff term. Then a rectangular product-integration (PI) rule is implemented in an implicit-explicit (IMEX) framework in order to implicitly treat this linear stiff term and handle in an explicit way the non-linear, and usually non-stiff, term. The kernel compression scheme (KCS) is successively adopted to reduce the overload of computation and storage need for the persistent memory term. To reduce the computational effort the semidiscretized problem is described in a matrix-form, so as to require the solution of Sylvester equations only with small matrices. Theoretical results on the accuracy, together with strategies for the optimal selection of the main parameters of the whole method, are derived and verified by means of numerical experiments carried out in two-dimensional domains. The computational advantages with respect to other approaches are also shown and some applications to the detection of pattern formation are illustrated at the end of the paper.

We address the issue of designing robust stabilization terms for the nonconforming virtual element method. To this end, we transfer the problem of defining the stabilizing bilinear form from the elemental nonconforming virtual element space, whose functions are not known in closed form, to the dual space spanned by the known functionals providing the degrees of freedom. By this approach, we manage to construct different bilinear forms yielding optimal or quasi-optimal stability bounds and error estimates, under weaker assumptions on the tessellation than the ones usually considered in this framework. In particular, we prove optimality under geometrical assumptions allowing a mesh to have a very large number of arbitrarily small edges per element. Finally, we numerically assess the performance of the VEM for several different stabilizations fitting with our new framework on a set of representative test cases.

In this paper we develop a general computer-assisted proof method for periodic solutions to delay differential equations. The class of problems considered includes systems of delay differential equations with an arbitrary number of (forward and backward) delays. When the nonlinearities include nonpolynomial terms we introduce auxiliary variables to first rewrite the problem into an equivalent polynomial one. We then apply a flexible fixed point technique in a space of geometrically decaying Fourier coefficients. We showcase the efficacy of this method by proving periodic solutions in the well-known Mackey–Glass delay differential equation for the classical parameter values.

We study fluid-saturated porous materials that undergo poro-elastic deformations in thin domains. The mechanics in such materials are described using a biphasic model based on the theory of porous media (TPM) and consisting of a system of differential equations for material’s displacement and fluid’s pressure. These equations are in general strongly coupled and nonlinear, such that exact solutions are hard to obtain and numerical solutions are computationally expensive. This paper reduces the complexity of the biphasic model in thin domains with a scale separation between domain’s width and length. Based on standard asymptotic analysis, we derive a reduced model that combines two sub-models. Firstly, a limit model consists of averaged equations that describe the fluid pore pressure and displacement in the longitudinal direction of the domain. Secondly, a corrector model re-captures the mechanics in the transverse direction. The validity of the reduced model is finally tested using a set of numerical examples. These demonstrate the computational efficiency of the reduced model, while maintaining reliable solutions in comparison with original biphasic TPM model in thin domain., Deutsche Forschungsgemeinschaft, Projekt DEAL

Statistical solutions have recently been introduced as an alternative solution framework for hyperbolic systems of conservation laws. In this work, we derive a novel a posteriori error estimate in the Wasserstein distance between dissipative statistical solutions and numerical approximations obtained from the Runge-Kutta Discontinuous Galerkin method in one spatial dimension, which rely on so-called regularized empirical measures. The error estimator can be split into deterministic parts which correspond to spatio-temporal approximation errors and a stochastic part which reflects the stochastic error. We provide numerical experiments which examine the scaling properties of the residuals and verify their splitting., Baden-Württemberg Stiftung, Projekt DEAL, Deutsche Forschungsgemeinschaft

We propose a new three-dimensional formulation based on the mixed tangential-displacement normal-normal-stress method for elasticity. In elastic tangential-displacement normal-normal-stress elements, the tangential component of the displacement field and the normal component of the stress vector are degrees of freedom and continuous across inter-element interfaces. Tangential-displacement normal-normal-stress finite elements have been shown to be locking-free with respect to shear locking in thin elements, which makes them suitable for the discretization of laminates or macro-fiber composites. In the current paper, we extend the formulation to piezoelectric materials by adding the electric potential as degree of freedom.

Part of the work was done while the authors were visiting the South China University of Science and Technology (SUSTech), Shenzhen, in 2018, and we are grateful (in particular to Professor Cai Heng Li) for the hospitality that we received.The authors would like to thank the Isaac Newton Institute for Mathematical Sciences, Cambridge, for support and hospitality during the programme Groups, representations and applications: new perspectives (supported by EPSRC grant no.EP/R014604/1), where further work on this paper was undertaken. In particular we acknowledge a Simons Fellowship (Cameron) and a Kirk Distinguished Visiting Fellowship (Praeger) during this programme. Schneider thanks the Centre for the Mathematics of Symmetry and Computation of The University of Western Australia and Australian Research Council Discovery Grant DP160102323 for hosting his visit in 2017 and acknowledges the support of the CNPq projects Produtividade em Pesquisa (project no.: 308212/2019-3) and Universal (project no.:421624/2018-3). Diagonal groups are one of the classes of finite primitive permutation groups occurring in the conclusion of the O'Nan-Scott theorem. Several of the other classes have been described as the automorphism groups of geometric or combinatorial structures such as affine spaces or Cartesian decompositions, but such structures for diagonal groups have not been studied in general. The main purpose of this paper is to describe and characterise such structures, which we call diagonal semilattices. Unlike the diagonal groups in the O'Nan-Scott theorem, which are defined over finite characteristically simple groups, our construction works over arbitrary groups, finite or infinite. A diagonal semilattice depends on a dimension m and a group T. For m=2, it is a Latin square, the Cayley table of T, though in fact any Latin square satisfies our combinatorial axioms. However, for m≥3, the group T emerges naturally and uniquely from the axioms. (The situation somewhat resembles projective geometry, where projective planes exist in great profusion but higher-dimensional structures are coordinatised by an algebraic object, a division ring.) A diagonal semilattice is contained in the partition lattice on a set Ω, and we provide an introduction to the calculus of partitions. Many of the concepts and constructions come from experimental design in statistics. We also determine when a diagonal group can be primitive, or quasiprimitive (these conditions turn out to be equivalent for diagonal groups). Associated with the diagonal semilattice is a graph, the diagonal graph, which has the same automorphism group as the diagonal semilattice except in four small cases with m

This article examines the accuracy for large times of asymptotic expansions from periodic homogenization of wave equations. As usual, $\epsilon$ denotes the small period of the coefficients in the wave equation. We first prove that the standard two scale asymptotic expansion provides an accurate approximation of the exact solution for times $t$ of order $\epsilon^{-2+\delta}$ for any $\delta>0$. Second, for longer times, we show that a different algorithm, that is called criminal because it mixes different powers of $\epsilon$, yields an approximation of the exact solution with error $O(\epsilon^N)$ for times $\epsilon^{-N}$ with $N$ as large as one likes. The criminal algorithm involves high order homogenized equations that, in the context of the wave equation, were first proposed by Santosa and Symes and analyzed by Lamacz. The high order homogenized equations yield dispersive corrections for moderate wave numbers. We give a systematic analysis for all time scales and all high order corrective terms.

In this article we derive the expression of \textit{renormalized energies} for unit-valued harmonic maps defined on a smooth bounded domain in \(\mathbb{R}^2\) whose boundary has several connected components. The notion of renormalized energies was introduced by Bethuel-Brezis-H\'elein in order to describe the position of limiting Ginzburg-Landau vortices in simply connected domains. We show here, how a non-trivial topology of the domain modifies the expression of the renormalized energies. We treat the case of Dirichlet boundary conditions and Neumann boundary conditions as well.

The goal of this note is to present a simple proof of existence and uniqueness of the solution of the Leray problem for high viscosities and homogeneous or non homogeneous boundary conditions. Furthermore we address the issue of uniqueness of the Poiseuille flow in a pipe.

Sparse direct solvers using Block Low-Rank compression have been proven efficient to solve problems arising in many real-life applications. Improving those solvers is crucial for being able to 1) solve larger problems and 2) speed up computations. A main characteristic of a sparse direct solver using low-rank compression is when compression is performed. There are two distinct approaches: (1) all blocks are compressed before starting the factorization, which reduces the memory as much as possible, or (2) each block is compressed as late as possible, which usually leads to better speedup. The objective of this paper is to design a composite approach, to speedup computations while staying under a given memory limit. This should allow to solve large problems that cannot be solved with Approach 2 while reducing the execution time compared to Approach 1.We propose a memory-aware strategy where each block can be compressed either at the beginning or as late as possible. We first consider the problem of choosing when to compress each block, under the assumption that all information on blocks is perfectly known, i.e., memory requirement and execution time of a block when compressed or not. We show that this problem is a variant of the NP-complete Knapsack problem, and adapt an existing 2-approximation algorithm for our problem. Unfortunately, the required information on blocks depends on numerical properties and in practice cannot be known in advance. We thus introduce models to estimate those values. Experiments on the PaStiX solver demonstrate that our new approach can achieve an excellent trade-off between memory consumption and computational cost. For instance on matrix Geo1438, Approach 2 uses three times as much memory as Approach 1 while being three times faster. Our new approach leads to an execution time only 30% larger than Approach 2 when given a memory 30% larger than the one needed by Approach 1., Les solveurs directs creux utilisant de la compression low-rank ont montré leur efficacité pour résoudre de grands systèmes. Les améliorer permetde 1) résoudre de plus grands problèmes et 2) accélérer la résolution. Une caractéristique principale de ces solveurs est quand la compression est réalisée. Il existe deux approches: (1) tous les blocs sont compressés avant le début de la factorisation, ce qui minimise la consommation mémoire, ou (2) chaque bloc est compressé au plus tard, ce qui permet d’avoir de meilleures performances. L’objectif de cet article est de proposer une approche intermédiaire, qui accélère les calculs au maximum tout en restant en dessous d’une limite mémoire donnée. Cela devrait permettre de résoudre de grands problèmes que l’approche 2 ne peut pas résoudre ou de réduire le temps d’exécution de l’approche 1. Nous proposonsune stratégie memory-aware où chaque bloc peut ˆêtre compressé au début ou au plus tard. Nous commençons par considérer le problème où toutes les informations (consommation mémoire et temps d’exécution en version compressée ou non) sont parfaitement connues. Nous montrons que ce problème est une variante du problème NP-complet Knapsack et adaptons une 2-approximation pour notre problème. Malheureusement, toutes ces informations dépendent de propriétés numériques et ne sont pas connues `a l’avance. Des modèles sont donc introduits afin d’estimer ces valeurs. Des expérimentations avec le solveur PaStiX montrent que notre approche permet d’avoir un excellent compromis entre consommation mémoire et temps d’exécution. Par exemple, pour la matrice Geo1438, l’approche 2 utilise trois fois plus de mémoire que l’approche 1, qui est trois fois plus lente. Notre nouvelle méthode permet d’obtenir à la fois un temps d’exécution seulement 30% supérieur à celui de l’approche 2 tout en ayant une consommation mémoire seulement supérieure de 30% à celle de l’approche 1.

In this paper we consider singularly perturbed nonlinear Schrödinger equations with electromagnetic potentials and involving continuous nonlinearities with subcritical, critical or supercritical growth. By means of suitable variational techniques, truncation arguments and Lusternik–Schnirelman theory, we relate the number of nontrivial complex-valued solutions with the topology of the set where the electric potential attains its minimum value.

In this article, we are interested in the study of the well-posedness as well as of the long time behavior, in terms of finite-dimensional attractors, of a coupled Allen–Cahn/Cahn–Hilliard system associated with dynamic boundary conditions. In particular, we prove the existence of the global attractor with finite fractal dimension.

In the present paper, we are concerned with the semilinear viscoelastic wave equation in an inhomogeneous medium Ω subject to two localized dampings. The first one is of the type viscoelastic and is distributed around a neighborhood ω of the boundary according to the Geometric Control Condition. The second one is a frictional damping and we consider it hurting the geometric condition of control. We show that the energy of the wave equation goes uniformly and exponentially to zero for all initial data of finite energy taken in bounded sets of finite energy phase-space.

Problems with variable exponents have attracted a great deal of attention lately and various existence, nonexistence and stability results have been established. The importance of such problems has manifested due to the recent advancement of science and technology and to the wide application in areas such as electrorheological fluids (smart fluids) which have the property that the viscosity changes drastically when exposed to heat or electrical fields. To tackle and understand these models, new sophisticated mathematical functional spaces have been introduced, such as the Lebesgue and Sobolev spaces with variable exponents. In this work, we are concerned with a system of wave equations with variable-exponent nonlinearities. This system can be regarded as a model for interaction between two fields describing the motion of two “smart” materials. We, first, establish the existence of global solutions then show that solutions of enough regularities stabilize to the rest state ( 0 , 0 ) either exponentially or polynomially depending on the range of the variable exponents. We also present some numerical tests to illustrate our theoretical findings.

The aim of this paper is to study higher-order Caginalp phase-field systems based on the Maxwell–Cattaneo law, instead of the classical Fourier law. More precisely, one obtains well-posedness results, as well as the existence of finite-dimensional attractors.

We present a gentle approach to the justification of effective media approximations, for PDE's set outside the union of $n \gg 1$ spheres with low volume fraction. To illustrate our approach, we consider three classical examples: the derivation of the so-called strange term, made popular by Cioranescu and Murat, the derivation of the Brinkman term in the Stokes equation, and a scalar analogue of the effective viscosity problem. Under some separation assumption on the spheres, valid for periodic and random distributions of the centers, we recover effective models as $n$ goes to infinity by simple arguments., 20 pages. Final version, to appear in Asympt. Anal

Prostate cancer (PCa) is the second most frequently diagnosed male cancer worldwide. Early diagnosis of PCa, response to therapy and prognosis still represent a challenge. Nearly 60% of PCa patients undergo radiation therapy (RT) which might cause side effects. In spite of numerous researches in this field, predictive biomarkers for radiation toxicity are still not elucidated. MicroRNAs as posttranscriptional regulators of gene expression are shown to be changed during and after irradiation. Manipulation with miRNA levels might be used to modulate response to RT-to reverse radioresistance-to induce radiosensitivity, or if needed, to reduce sensitivity to treatment to avoid side effects. In this review we have listed and described miRNAs involved in response to RT in PCa, and highlighted potential candidates for future biological tests predicting radiation response to RT, with the special focus on side effects of RT. Individual radiation response is a result of the interactions between physical characteristics of radiation treatment and biological background of each patient, and miRNA expression changes among others. According to described literature we concluded that let-7, miR-21, miR-34a, miR-146a, miR-155, and members of miR-17/92 cluster might be promising candidates for biological tests predicting radiosensitivity of PCa patients undergoing radiation treatment, and as future agents for modulation of radiation response. Predictive miRNA panels, especially for acute and late side effects of RT can serve as a starting point for decisions for individualized RT planning. We believe that this review might be one step closer to understanding molecular mechanisms underlying individual radiation response of patients with PCa.

In this paper, we investigate the relationships between proper efficiency and the solutions of a general scalarization problem in multi-objective optimization. We provide some conditions under which the solutions of the dealt with scalar program are properly efficient and vice versa. We also show that, under some conditions, if the considered general scalar problem is unbounded, then the original multi-objective problem does not have any properly efficient solution. In another part of the work, we investigate a general transformation of the objective functions which preserves proper efficiency. We show that several important results existing in the literature are direct consequences of the results of the present paper.

With the climate change processes over times, all professions and trades in Three Gorges Reservoir Area will be influenced. One of the biggest challenges is the risk of rising sea level. In this situation, a large number of uncertainties for climate changes will be faced in Three Gorges Reservoir Area. Therefore, it is of importance to investigate the complexity of decision making on investing in the long term rising sea level risk related projects in Three Gorges Reservoir Area. This paper investigates the sea level and the temperature as the underlying assets in Three Gorges Reservoir Area. A real option model is constructed to evaluate potential sea level rising risk. We formulate European and American real option models into a linear parabolic variational inequalities and propose a power penalty approach to solve it. Then we obtain a nonlinear parabolic equation. It shows that the nonlinear parabolic equation is unique and solvable. Also, the solutions of the nonlinear parabolic equation converge to the solutions of the parabolic variational inequalities at the rate of order O(λ −k/2). Since the analytic solution of nonlinear parabolic equation is difficult to obtain, a fitted finite volume method is developed to solve it in case of European and American options, and the convergence of the nonlinear parabolic equation is obtained. An empirical analysis is presented to illustrate our theoretical results.

Beam-like compliant elements have found wide-ranging application in many fields of engineering and science often where 3D large deflections can be of concern such as soft robotics, DNA mechanics and helicopter/wind turbine rotor blades. The homotopy analysis method (HAM) is used for the first time to obtain a novel analytical solution in converged series form for the arbitrary large deflection of geometrically exact beams subject to both conservative and follower loading scenarios. The homotopy analysis method, which offers desirable characteristics such as being free from small or large parameters, coupled with auxiliary parameters controlling convergence, is applied directly to the intrinsic governing equations of a geometrically exact beam theory. The system of first-order differential governing equations of geometrically exact beams with intrinsic formulation is free from rotation and displacement variables, and offers a low degree of nonlinearity (quadratic at most) and compact mathematical form, making it suitable for analytical solutions. Due to the relatively poor convergence of the original HAM algorithm, the iterative HAM technique is employed which is known to accelerate convergence and to improve the computational efficiency of the homotopy analysis method. The obtained homotopy series offers a number of novel features in the context of the analytical solutions for the large deflection of beams, including (a) the direct calculation of internal forces and moments which is significant for engineering design purposes, (b) being able to capture 3D deflections, (c) considering transverse shear effects which can be important for thicker beams or when the Young’s modulus to shear modulus ratio is significant (such as composite materials) and (d) considering conservative and follower tip and distributed loads, in a unified framework. In order to investigate the efficacy, applicability and accuracy of the proposed method, a number of numerical examples are considered in which a cantilever beam subject to tip or distributed loads undergoes large deflection. Large deflection results for both conservative and follower loads are compared against those of less comprehensive analytical solutions as well as against numerical methods including finite element and Chebyshev collocation methods where good agreement is observed. These results demonstrate the applicability and effectiveness of HAM for the large deflection analysis of geometrically exact beams.

This paper is concerned with the recovery of (approximate) solutions to parabolic problems from incomplete and possibly inconsistent observational data, given on a time-space cylinder that is a strict subset of the computational domain under consideration. Unlike previous approaches to this and related problems our starting point is a regularized least squares formulation in a continuous infinite-dimensional setting that is based on stable variational time-space formulations of the parabolic partial differential equation. This allows us to derive a priori as well as a posteriori error bounds for the recovered states with respect to a certain reference solution. In these bounds the regularization parameter is disentangled from the underlying discretization. An important ingredient for the derivation of a posteriori bounds is the construction of suitable Fortin operators which allow us to control oscillation errors stemming from the discretization of dual norms. Moreover, the variational framework allows us to contrive preconditioners for the discrete problems whose application can be performed in linear time, and for which the condition numbers of the preconditioned systems are uniformly proportional to that of the regularized continuous problem. In particular, we provide suitable stopping criteria for the iterative solvers based on the a posteriori error bounds. The presented numerical experiments quantify the theoretical findings and demonstrate the performance of the numerical scheme in relation with the underlying discretization and regularization.

In a previous paper, we obtained several "compact versions" of Rubio de Francia's weighted extrapolation theorem, which allowed us to extrapolate the compactness of linear operators from just one space to the full range of weighted Lebesgue spaces, where these operators are bounded. In this paper, we study the extrapolation of compactness for bilinear operators in terms of bilinear Muckenhoupt weights. As applications, we easily recover and improve earlier results on the weighted compactness of commutators of bilinear Calder\'{o}n-Zygmund operators, bilinear fractional integrals and bilinear Fourier multipliers. More general versions of these results are recently due to Cao, Olivo and Yabuta (arXiv:2011.13191), whose approach depends on developing weighted versions of the Fr\'echet--Kolmogorov criterion of compactness, whereas we avoid this by relying on "softer" tools, which might have an independent interest in view of further extensions of the method., Comment: v3: final version, incorporated referee comments, to appear in Indagationes Mathematicae, 27 pages

We introduce and analyze a stress-based formulation for Zener's model in linear viscoelasticity. The method is aimed to tackle efficiently heterogeneous materials that admit purely elastic and viscoelastic parts in their composition. We write the mixed variational formulation of the problem in terms of a class of tensorial wave equation and obtain an energy estimate that guaranties the well-posedness of the problem through a standard Galerkin procedure. We propose and analyze mixed continuous and discontinuous Galerkin space discretizations of the problem and derive optimal error bounds for each semidiscrete solution in the corresponding energy norm., Comment: arXiv admin note: text overlap with arXiv:2005.01021

To identify determinants of labor productivity, correlation relationships of productivity of various indicators was evaluated, reflecting the influence of a wide range of socio-economic and innovative factors in the regions of Russia for 2015-2017. It has been established that many indicators characterizing socio-economic and innovative factors do not have a significant relationship with labor productivity and are multicollenarity (they have correlation relationships among themselves). For each year, according to statistics of 82 regions of Russia, econometric models in the form of a well-known standard internal linear function - an analog of the Cobb-Douglas production function are constructed. The obtained models indicate a positive impact on labor productivity, capital-labor ratio, foreign investment, wage levels, income inequality, inflation in industrial goods and export markets. The obtained simulation results showing a significant positive effect on labor productivity exerted by income inequality and the inflation rate on the industrial goods market, which indicates the action of economic mechanisms in the country in developed countries. When performing research, it was found that the positive impact of capital-labor ratio on productivity in 2015-2017 was reduced. This indicates a decrease in the efficiency of use of fixed capital in countries, about the crisis in the country's economy.

In this paper, we study the long-time behavior of the following plate equation ε ( t ) u t t + g ( u t ) + Δ 2 u + λ u + f ( u ) = h , where the coefficient ε depends explicitly on time, the nonlinear damping and the nonlinearity both have critical growths.

In this paper, we give a short proof of the Lorentz estimates for gradients of very weak solutions to the linear parabolic equations with the Muckenhoupt class A q -weights u t − div ( A ( x , t ) ∇ u ) = div ( F ) , in a bounded domain Ω × ( 0 , T ) ⊂ R N + 1 , where A has a small mean oscillation, and Ω is a Lipchistz domain with a small Lipschitz constant.

In this paper, we consider a class of nonlocal equations where the convolution kernel is given by a Bessel potential symbol of order $\alpha$ for $\alpha > 1$. Based on the properties of the convolution operator, we apply a global bifurcation technique to show the existence of a highest, even, $2\pi$-periodic traveling-wave solution. The regularity of this wave is proved to be exactly Lipschitz., Comment: 22 pages. arXiv admin note: text overlap with arXiv:1810.00248 by other authors

We present a construction of a multi-scale Gaussian beam parametrix for the Dirichlet boundary value problem associated with the wave equation, and study its convergence rate to the true solution in the highly oscillatory regime. The construction elaborates on the wave-atom parametrix of Bao, Qian, Ying, and Zhang and extends to a multi-scale setting the technique of Gaussian beam propagation from a boundary of Katchalov, Kurylev and Lassas., Comment: 64 pages, 7 figures, minor update

Geological disasters such as slope failure and landslides can cause loss of life and property. Therefore, reproducing their evolution process is of great importance for risk assessment and mitigation. The recently developed SIMSAND critical state sand model combined with the smoothed particle hydrodynamics (SPH) method is adopted in this work to study slope failure under large deformations. To illustrate the efficiency and accuracy of the SIMSAND-SPH approach, a series of slope collapse studies using the discrete element method (DEM) considering various particle shapes (i.e. spherical, tetrahedral and elongated) is adopted as benchmarks. The parameters of the SIMSAND model are calibrated using DEM triaxial tests. In comparison to the DEM simulations, the runout distance and final slope height are well characterized with the SIMSAND-SPH approach with less computational cost. All comparisons show that the SIMSAND-SPH approach is highly efficient and accurate, which can be an alternative numerical tool to simulate real scale granular flow.

The susceptible-infected-recovered (SIR) model is perhaps the most basic epidemiological model for the evolution of disease spread within a population. Because of its direct representation of fundamental physical quantities, a true solution to an SIR model possesses a number of qualitative properties, such as conservation of the total population or positivity or monotonicity of its constituent populations, that may only be guaranteed to hold numerically under step-size restrictions on the solver. Operator-splitting methods with order greater than two require backward sub-steps in each operator, and the effects of these backward sub-steps on the step-size restrictions for guarantees of qualitative correctness of numerical solutions are not well studied. In this study, we analyze the impact of backward steps on step-size restrictions for guaranteed qualitative properties by applying third- and fourth-order operator-splitting methods to the SIR epidemic model. We find that it is possible to provide step-size restrictions that guarantee qualitative property preservation of the numerical solution despite the negative sub-steps, but care must be taken in the choice of the method. Results such as this open the door for the design and application of high-order operator-splitting methods to other mathematical models in general for which qualitative property preservation is important.

The fractional Laplacian, ( − △ ) s , s ∈ ( 0 , 1 ) , appears in a wide range of physical systems, including Levy flights, some stochastic interfaces, and theoretical physics in connection to the problem of stability of the matter. In this paper, a matrix transfer technique (MTT) is employed combining with spectral/element method to solve fractional diffusion equations involving the fractional Laplacian. The convergence of the MTT method is analyzed by the abstract operator theory. Our method can be applied to solve various fractional equation involving fractional Laplacian on some complex domains. Numerical results indicate exponential convergence in the spatial discretization which is in good agreement with the theoretical analysis.

We bound the variance and other moments of a random vector based on the range of its realizations, thus generalizing inequalities of Popoviciu and of Bhatia and Davis concerning measures on the line to several dimensions. This is done using convex duality and (infinite-dimensional) linear programming. The following consequence of our bounds exhibits symmetry breaking, provides a new proof of Jung’s theorem, and turns out to have applications to the aggregation dynamics modelling attractive–repulsive interactions: among probability measures on [Formula: see text] whose support has diameter at most [Formula: see text], we show that the variance around the mean is maximized precisely by those measures that assign mass [Formula: see text] to each vertex of a standard simplex. For [Formula: see text], the [Formula: see text] th moment—optimally centered—is maximized by the same measures among those satisfying the diameter constraint.

The Hilbert L-matrix A s = [ a i j ( s ) ] , where a i j ( s ) = 1 / ( max ⁡ { i , j } + s ) with i , j ≥ 0 , was introduced in [3] . As a surprising property, we showed that its 2-norm is constant for s ≥ s 0 , where the critical point s 0 is unknown but relies in the interval ( 1 / 4 , 1 / 2 ) . In this note, using some delicate calculations we sharpen this result by improving the upper and lower bounds of the interval surrounding s 0 . Moreover, we establish that the same property persists for the p-norm of A s matrices.

We consider an Klein–Gordon relativistic equation with a boundary dissipation of fractional derivative type. We study of stability of the system using semigroups theory and classical theorems over asymptotic behavior.