Analysis of Lie symmetries with conservation laws and solutions for the generalized (3+1)-dimensional time fractional Camassa–Holm–Kadomtsev–Petviashvili equation

In this paper, under investigated is a generalized $(3+1)$-dimensional Camassa–Holm–Kadomtsev–Petviashvili (gCH-KP) equation, which describes the role of dispersion in the formation of patterns in liquid drops. With the help of the semi-inverse method, the Euler–Lagrange equation and Agrawal’s method, the time fractional gCH-KP equation is derived in the sense of Riemann–Liouville fractional derivatives. Further, the symmetry of the $(3+1)$-dimensional time fractional gCH-KP equation is studied by fractional order symmetry. Meanwhile, based on the new conservation theorem, the conservation laws of $(3+1)$-dimensional time fractional gCH-KP equation are constructed. Finally, the solutions to the equation are given via a bilinear method and the radial basis functions (RBFs) meshless approach.     

Global well-posedness of the d-D Oldroyd-B type models with fractional Laplacian dissipation

This paper focuses on the Cauchy problem of the d-dimensional incompressible Oldroyd-B type models for viscoelastic flow with fractional Laplacian dissipation, namely, with ${(?\Delta )}^{{\eta }_1}u$ and ${(?\Delta )}^{{\eta }_2}\tau$. For ${\eta }_1\ge \frac{1}{2}+\frac{d}{4}$, ${\eta }_2>0$ and ${\eta }_1+{\eta }_2\ge 1+\frac{d}{2}$, we obtain the global regularity of strong solutions when the initial data $(u_0,{\tau }_0)$ are sufficiently smooth.     

A remark concerning divergence accuracy order for H(div)-conforming finite element flux approximations

The construction of finite element approximations in $\mathbf{H}(\text{div},\Omega )$ usually requires the Piola transformation to map vector polynomials from a master element to vector fields in the elements of a partition of the region $\Omega$. It is known that degradation may occur in convergence order if non affine geometric mappings are used. On this point, we revisit a general procedure for the improvement of two-dimensional flux approximations discussed in a recent paper of this journal (Comput. Math. Appl. 74 (2017) 3283–3295). The starting point is an approximation scheme, which is known to provide $L^2$-errors with accuracy of order $k+1$ for sufficiently smooth flux functions, and of order $r+1$ for flux divergence. An example is $R{T}_k$ spaces on quadrilateral meshes, where $r=k$ or $k?1$ if linear or bilinear geometric isomorphisms are applied. Furthermore, the original space is required to be expressed by a factorization in terms of edge and internal shape flux functions. The goal is to define a hierarchy of enriched flux approximations to reach arbitrary higher orders of divergence accuracy $r+n+1$ as desired, for any $n\ge 1$. The enriched versions are defined by adding higher degree internal shape functions of the original family of spaces at level $k+n$, while keeping the original border fluxes at level $k$. The case $n=1$ has been discussed in the mentioned publication for two particular examples. General stronger enrichment $n>1$ shall be analyzed and applied to Darcy’s flow simulations, the global condensed systems to be solved having same dimension and structure of the original scheme.     

Note on nonparametric spectral analysis of wideband spectrum with missing data via sample-and-hold interpolation and deconvolution

In [1] a procedure for bias-free estimation of the autocorrelation function is introduced for equidistantly sampled data with randomly occurring samples being invalid. The method incorporates sample-and-hold interpolation of the missing data points. The occurring dynamic error of the primary estimate of the correlation function is treated by a deconvolution procedure with two parameters $c_0$ and $c_1$ with $c_0+2c_1=1$, which are the on-diagonal and the aside-diagonal parameters of a specific correction matrix (at all lag times except zero). The parameters $c_0$ and $c_1$ were obtained as a function of the probability α of a sample to be valid by numerical simulation. However, explicit expressions for the parameters $c_0(\alpha )=1-\frac{2}{\alpha }+\frac{2}{{\alpha }^2}$ and $c_1(\alpha )=\frac{1}{\alpha }-\frac{1}{{\alpha }^2}$ can be derived, which might improve the usability of the deconvolution procedure in [1].     

Unsteady micropolar fluid flow in a thin domain with Tresca fluid–solid interface law

We consider a micropolar fluid flow in a two-dimensional domain. We assume that the velocity field satisfies a non-linear slip boundary condition of friction type on a part of the boundary while the micro-rotation field satisfies non-homogeneous Dirichlet boundary conditions. We prove the existence and uniqueness of a solution. Then motivated by lubrication problems we assume that the thickness and the roughness of the domain are of order $0<\epsilon <<1$ and we study the asymptotic behaviour of the flow as $\epsilon$ tends to zero. By using the two-scale convergence technique we derive the limit problem which is totally decoupled for the limit velocity and pressure $(v^0,p^0)$ on one hand and the limit micro-rotation $Z^0$ on the other hand. Moreover we prove that $v^0$, $p^0$ and $Z^0$ are uniquely determined via auxiliary well-posed problems.     

GPU implementation of Borůvka’s algorithm to Euclidean minimum spanning tree based on Elias method

We present both sequential and data parallel approaches to build hierarchical minimum spanning forest (MSF) or tree (MST) in Euclidean space (EMSF/EMST) for applications whose input $N$ points are uniformly or boundedly distributed in Euclidean space. Each iteration of the sequential approach takes $O\left(N\right)$ time complexity through combining Borůvka’s algorithm with an improved component-based neighborhood search algorithm, namely sliced spiral search, which is a newly proposed improvement to Bentley’s spiral search for finding a component graph’s closest outgoing point on the plane. It works based on the uniqueness property in Euclidean space, and allows $O\left(1\right)$ time complexity for one search from a query point to find the component’s closest outgoing point at different iterations of Borůvka’s algorithm. The data parallel approach includes a newly proposed two-direction breadth-first search (BFS) implementation on graphics processing unit (GPU) platform, which is specialized for selecting a spanning tree’s shortest outgoing edge. This GPU two-direction parallel BFS enables a tree traversal operation to start from any one of its vertex acting as root. These GPU parallel implementations work by assigning $N$ threads with one thread associated to one input point, one thread occupies $O\left(1\right)$ local memory and the whole algorithm occupies $O\left(N\right)$ global memory. Experiments are conducted on both uniformly distributed data sets and TSPLIB database. We evaluate computation time of the proposed approaches on more than 80 benchmarks with size $N$ growing up to $10^6$ points on personal laptop.     

Dynamic behavior of droplet through a confining orifice:A lattice Boltzmann study

The motion of gravity-driven deformable droplets passing through a confining orifice in two-dimensional ($2D$) space is numerically studied by the phase-field-based multiple-relaxation-time (MRT) lattice Boltzmann (LB) model, and the ratio of orifice-to-droplet diameter is less than 1. Droplets are placed just above a sink with an orifice in the middle, accelerate under gravity and encounter the orifice plate. In this work, we mainly consider the effects of the Bond number ($Bo$), orifice-to-droplet diameter ratio ($r=d∕D$), plate thickness ($H_t$), wettability (or contact angle) and the diameter ratio of two droplets ($r_d=D_1∕D_2$) on the dynamic behavior of droplet through the orifice. The results show that these issues have great influences on the typical flow patterns (i.e., release and capture). With the decrease of contact angle, the droplet is more easily captured, and there exists a critical equilibrium contact angle ${\theta }^{eq}$ when the Bond number and the orifice-to-droplet diameter ratio as well as the thickness of the plate are specified. For the case with $\theta >{\theta }^{eq}$, the droplet can finally pass through the orifice, otherwise, the droplet cannot pass through the orifice. In addition, the droplet is more likely to pass through the orifice as the thickness of the obstacle increases. Actually, when the obstacle thickness is large enough, droplet breaks into three segments and a liquid slug is formed in a hydrophilic orifice. Finally, for the evolution of two droplets with a larger diameter ratio ($r_d=1.0$), the combined droplet finally passes through the orifice due to greater inertia than the cases with $r_d=0$ and $r_d=0.43$. Besides, we also establish the relation $r=0.57^{\frac{2}{3}}B{o}^{?\frac{1}{3}}$ which can be used to separate droplet release from capture at $H_t=1.2mm$.     

Ground state solutions of Pohoz̆aev type for fractional Choquard equations with general nonlinearities

In this paper, we study the fractional Choquard equation

${(?\Delta )}^{s}u+u=\left(\right|x{|}^{?\mu }1F\left(u\right))f\left(u\right),\text{in}{\mathbb{R}}^N,$

where $N\ge 3$, $0<s<1$, $0<\mu <min\{N,4s\}$, and $f\in C(\mathbb{R},\mathbb{R})$ satisfies the general Berestycki–Lions conditions. Combining constrained variational method with deformation lemma, we obtain a ground state solution of Pohoz?aev type for the above equation. The result improves some ones in Shen et al. (2016).     

Divide and capture: An improved cryptanalysis of the encryption standard algorithm RSA

RSA is a well known standard algorithm used by modern computers to encrypt and decrypt messages. In some applications, to save the decryption time, it is desirable to have a short secret key d compared to the modulus N. The first significant attack that breaks RSA with short secret key given by Wiener in 1990 is based on the continued fraction technique and it works with $d<\frac{1}{\sqrt[4]{18}}{N}^{.25}$. A decade later, in 2000, Boneh and Durfee presented an improved attack based on lattice technique which works with d < N^.292. Until this day, Boneh–Durfee attack remain as the best attack on RSA with short secret key. In this paper, we revisit the continued fraction technique and propose a new attack on RSA. Our main result shows that when $d<\sqrt{t(2\sqrt{2}+8/3)}{N}^{.75}/\sqrt{e},$ where e is the public exponent and t is a chosen parameter, our attack can break the RSA with the running time of O(tlog (N)). Our attack is especially well suited for the case where e is much smaller than N. When e ≈ N, the Boneh–Durfee attack outperforms ours. As a result, we could simultaneously run both attacks, our new attack and the classical Boneh–Durfee attack as a backup.     

An accelerated cyclic-reduction-based solvent method for solving quadratic eigenvalue problem of gyroscopic systems

The quadratic eigenvalue problem (QEP) $({\lambda }^{2}M+\lambda G+K)x=0$, with $M^T=M$ being positive definite, $K^T=K$ being negative definite and $G^T=?G$, is associated with gyroscopic systems. In Guo (2004), a cyclic-reduction-based solvent (CRS) method was proposed to compute all eigenvalues of the above mentioned QEP. Firstly, the problem is converted to find a suitable solvent of the quadratic matrix equation (QME) $M{X}^2+GX+K=0$. Then using a Cayley transformation and a proper substitution, the QME is transformed into the nonlinear matrix equation (NME) $Z+A^T{Z}^{?1}A=Q$ with $A=M+K+G$ and $Q=2(M?K)$. The problem finally can be solved by applying the CR method to obtain the maximal symmetric positive definite solution of the NME as long as the QEP has no eigenvalues on the imaginary axis or for some cases where the QEP has eigenvalues on the imaginary axis. However, when all eigenvalues of the QEP are far away from or near the origin, the Cayley transformation seems not to be the best one and the convergence rate of the CRS method proposed in Guo (2004) might be further improved. In this paper, inspired by using a doubling algorithm to solve the QME, we use a Möbius transformation instead of the Cayley transformation to present an accelerated CRS (ACRS) method for solving the QEP of gyroscopic systems. In addition, we discuss the selection strategies of optimal parameter for the ACRS method. Numerical results demonstrate the efficiency of our method.     

Existence and asymptotic behavior of positive ground state solutions for coupled nonlinear fractional Kirchhoff-type systems

This paper is concerned with the following linearly coupled fractional Kirchhoff-type system

$\left\{\begin{array}{c}\left(a+b{\int }_{{\mathbb{R}}^3}|{(?△)}^{\frac{\alpha }{2}}u{|}^2\mathrm{d}x\right){(?△)}^{\alpha }u+\lambda u=f\left(u\right)+\gamma v,\text{in}{\mathbb{R}}^3,\hfill \\ \left(c+d{\int }_{{\mathbb{R}}^3}|{(?△)}^{\frac{\alpha }{2}}v{|}^2\mathrm{d}x\right){(?△)}^{\alpha }v+\mu v=g\left(v\right)+\gamma u,\text{in}{\mathbb{R}}^3,\hfill \\ u,v\in H^{\alpha }\left({\mathbb{R}}^3\right),\hfill \end{array}\right.$

where $a,c,\lambda ,\mu >0$, $b,d\ge 0$ are constants, $\alpha \in [\frac{3}{4},1)$ and $\gamma >0$ is a coupling parameter. Under the general Berestycki–Lions conditions on the nonlinear terms $f$ and $g$, we prove the existence of positive vector ground state solutions of Poho?aev type for the above system via variational methods. Moreover, the asymptotic behavior of these solutions as $\gamma \to 0^{+}$ is explored as well. Recent results from the literature are generally improved and extended.     

Effect of Porosity on free and forced vibration characteristics of the GPL reinforcement composite nanostructures

This article investigates the influence of porosity on free and forced vibration characteristics of a nanoshell reinforced by graphene platelets (GPL). The material properties of piece-wise graphene-reinforced composites (GPLRCs) are assumed to be graded in the thickness direction of a cylindrical nanoshell and estimated using a nanomechanical model. In addition, because of imperfection of the current structure, three kinds of porosity distributions are considered. The nanostructure is modeled using modified strain gradient theory (MSGT) which is a size-dependent theory with three length scale parameters. The novelty of the current study is to consider the effects of porosity, GPLRC and MSGT on dynamic and static behaviors of the nanostructure. Considering three length scale parameters ( $l_0=5h$, $l_1=3h$, $l_2=5h$ ) in MSGT leads to a better agreement with MD simulation in comparison by other theories. Finally, effects of different factors on static and dynamic behaviors of the porous nanostructure are examined in detail.