Search

In a recent paper referred to in the title, the author used the concept of quasi self-adjointness to obtain conservation laws for a system of the ( 2 + 1 ) -dimensional Bogoyavlensky–Konopelchenko equation. Apart from the adjoint system of equations, all the results on the quasi self-adjointness and the subsequent conservation laws obtained are inaccurate. In this comment, we clarify these inaccuracies and also generate conservation laws for a potential form of the underlying equation through Noether theorem. [ABSTRACT FROM AUTHOR]

This paper derives well-posedness and asymptotic results that provide qualitative information about the behavior, mechanism and strategies used by living organisms to navigate their biological networks. Chemical driven swimming is a captivating phenomenon that is observed in various living organisms like bacteria and protozoa but the problem in weighted networks is more complex, since the equations of parabolic-parabolic Keller-Segel model coupled with incompressible Navier-Stokes equations must be reformulated in a discrete setting. The starting point is to transpose the coupled system from the Euclidean case to the connected networks with certain network-theoretic simplified approaches, which yields fruitful key results. These results not only enable the construction of global solutions but also serve as a foundation for determining information about the stability and large time behavior of the system. Then, decay rates are well established to predict important features, such as how quickly weak solutions at a given point decrease over time due to dissipative processes. Additionally, the L 1 - convergence of cell densities towards the self-similar Gaussian solution of the heat equation is well proved by time dependent scaling, which shows that the solution maintains its shape and only scales in time and space as time evolves. Finally, this paper includes many numerical tests through a recent robust numerical scheme to illustrate the theoretical results and to develop computational control and prediction of living organisms' trajectories around central nodes in networked flows. [ABSTRACT FROM AUTHOR]

In this paper, distributed strain nucleus is presented to compute the modes I/II stress intensity factors (SIFs) and electric displacement intensity factors (EDIFs) for multiple unequal cracks placed between a functionally graded piezoelectric materials (FGPMs) and a piezoelectric half-plane. Employing the Fourier transform, the governing electro-elastic equations are solved in terms of the Burgers vectors b x , b y and b ϕ. These equations converted into a system of singular integral equations with Cauchy type then are solved numerically with the help of collocation methods. Crack problems can be solved by establishing the stress and displacement fields in the crack's absence and distributing strain nuclei along the line of the crack. The distributed strain nucleus is a powerful method to calculate accurate solutions to plane crack problems. In fact, the technique is excellent for solutions of complex crack patterns and it is suitable to analyze multiple cracks. Hence this technique allows considering multiple cracks with arbitrary lengths and positions at the interface. The validation of results has been checked by available analytical, finite element and experimental papers. The primary objective of this paper is to study the interacting effects among multiple unequal interface cracks, the non-homogeneity parameter, geometry and material properties and loading parameters on the field intensity factors. Numerical results show that the normalized modes I/II SIFs and EDIFs decreasing with increasing stiffness and thickness of coating. [ABSTRACT FROM AUTHOR]

This paper introduces a novel approach for the efficient determination of confined electron states in quantum semiconductor structures through the introduction of a generalized energy eigenvalue problem formulated within the framework of boundary integral analysis. The proposed method enables the direct determination of the energy eigenvalues and normalized wavefunctions for bound quantum states. The novel technique aims to address the challenges associated with efficiently modeling the behavior of electrons within confined regions, offering insights into optimizing the performance of a wide range of quantum semiconductor structures. By employing boundary integral techniques, the paper establishes an efficient numerical framework that effectively addresses the complexity of quantum confinement effects. The methodology is demonstrated through numerical simulations, showcasing its effectiveness and accuracy in predicting electron states within different quantum wire structures. [ABSTRACT FROM AUTHOR]

This paper presents a numerical comparison of finite-element methods resulting in local mass conservation at the element level for Poisson's problem, namely the primal hybrid and mixed methods. These formulations result in an indefinite system. Alternative formulations yielding a positive-definite system are obtained after hybridizing each method. The choice of approximation spaces yields methods with enhanced accuracy for the pressure variable, and results in systems with identical size and structure after static condensation. A regular pressure precision mixed formulation is also considered based on the classical RT k space. The simulations are accelerated using Open multi-processing (OMP) and Thread-Building Blocks (TBB) multithreading paradigms alongside either a coloring strategy or atomic operations ensuring a thread-safe execution. An additional parallel strategy is developed using C++ threads, which is based on the producer-consumer paradigm, and uses locks and semaphores as synchronization primitives. Numerical tests show the optimal parallel strategy for these finite-element formulations, and the computational performance of the methods are compared in terms of simulation time and approximation errors. Additional results are developed during the process. Numerical solvers often fail to find an accurate solution to the highly indefinite systems arising from finite-element formulations, and this paper documents a matrix ordering strategy to stabilize the resolution. A procedure to enable static condensation based on the introduction of piecewise constant functions that also fulfills Neumann's compatibility condition, and yet computes an average pressure per element is presented. • Numerical comparison of the mixed, primal hybrid, and other locally conservatives finite element methods. • Introduction of a twice hybridized primal hybrid method, yielding a positive definite matrix after static condensation. • Comparison of multi-threaded strategies and thread-safety mechanisms. • A numerical procedure to stabilize the linear system resolution for indefinite matrices is described. • A general procedure that enables static condensation and fulfills Neumann's compatibility condition. [ABSTRACT FROM AUTHOR]

This paper considers the discretization of the time-dependent Navier–Stokes equations with the family of inf-sup stabilized Scott–Vogelius pairs recently introduced in [John/Li/Merdon/Rui, Math. Models Methods Appl. Sci., 2024] for the Stokes problem. Therein, the velocity space is obtained by enriching the H 1 -conforming Lagrange element space with some H (div) -conforming Raviart–Thomas functions, such that the divergence constraint is satisfied exactly. In these methods arbitrary shape-regular simplicial grids can be used. In the present paper two alternatives for discretizing the convective terms are considered. One variant leads to a scheme that still only involves volume integrals, and the other variant employs upwinding known from DG schemes. Both variants ensure the conservation of linear momentum and angular momentum in some suitable sense. In addition, a pressure-robust and convection-robust velocity error estimate is derived, i.e., the velocity error bound does not depend on the pressure and the constant in the error bound for the kinetic energy does not blow up for small viscosity. After condensation of the enrichment unknowns and all non-constant pressure unknowns, the method can be reduced to a P k − P 0 -like system for arbitrary velocity polynomial degree k. Numerical studies verify the theoretical findings. [ABSTRACT FROM AUTHOR]

The interior inverse elastic scattering problem is considered in this paper. The scattering in the cavity is illuminated by an time-harmonic elastic point source at fixed frequency. The boundary of the cavity is assumed to be a small and smooth perturbation of some circle. In this paper, we are interested in the reconstruction of the boundary according to the measured data on a closed curve inside the cavity. In implementation, we split the displacement field into compressional and shear waves by the Helmholtz decomposition. Then we use coordinate transformation, power series expansion based on perturbed parameters and Fourier expansion of corresponding power series coefficients to simplify the original boundary value problem to a successive sequence of one-dimensional boundary value problem. We linearize the inverse problem by eliminating the higher-order terms in the power series expansion; and design a reconstruction algorithm by the explicit relationship between the Fourier coefficients of the boundary function and the total field. In order to obtain better resolution, we also use a nonlinear correction method in the numerical experiments. The numerical results have shown the effectiveness and the stableness of the proposed method. [ABSTRACT FROM AUTHOR]

In this paper, we present explicit expressions for conforming finite element function spaces, basis functions, and degrees of freedom on the pentatope (the 4-simplex) and tetrahedral prism elements. More generally, our objective is to construct high-order finite element function spaces that maintain conformity with infinite-dimensional spaces of a carefully chosen de Rham complex in four dimensions. This paper is a natural extension of the companion paper entitled "Conforming finite element function spaces in four dimensions, part I: Foundational principles and the tesseract" by Nigam and Williams (2024). In contrast to Part I, in this paper we focus on two of the most popular elements which do not possess a full tensor-product structure in all four coordinate directions. We note that these elements appear frequently in existing space-time finite element methods. In order to build our finite element spaces, we utilize powerful techniques from the recently developed 'Finite Element Exterior Calculus'. Subsequently, we translate our results into the well-known language of linear algebra (vectors and matrices) in order to facilitate implementation by scientists and engineers. [ABSTRACT FROM AUTHOR]

In recent years, the heat kernel has proven to be a valuable tool for graph characterization and graph-based object clustering. It serves as the fundamental solution to the heat diffusion equation associated with the discrete graph Laplacian, describing how information flows across graph edges over time by exponentiating the Laplacian eigensystem. This paper focuses on the novel concept of the generalized heat kernel of a network, recently introduced by Kalala Mutombo et al. These authors build upon the k-path Laplacian operator for graphs, pioneered by Estrada et al., to incorporate long-range interactions (LRI) into information transmission across nodes and edges of a network/graph. LRI are handled through the use of Mellin and Laplace transforms. This paper makes a notable contribution by highlighting the practical application of the generalized heat kernel. We achieve this by demonstrating the effectiveness of the generalised heat kernel invariants in object clustering with real data using principal component analysis (PCA). Moreover, we investigate how incorporating long-range interactions (LRI) impacts object characterization and clustering, revealing superior results compared to conventional diffusion methods. Through experimentation, we show that object clustering remains achievable even with small values of the Mellin and Laplace parameters, contrasting with the requirement of an infinite value in the absence of LRI. • Application of graph generalised heat kernel invariants in object clustering via PCA. • Long-range interactions enhance graph-based object clustering. • Object clustering achieved with small Mellin/Laplace parameters. [ABSTRACT FROM AUTHOR]

In this paper the construction of C k basis functions on meshes with hanging nodes for hp -finite element methods is proposed. The construction is limited to axis-parallel, rectangular meshes, which can, however, be of any dimension. The basis functions are composed by using hierarchical tensor product shape functions which are defined on patches of so-called constraining nodes and combine univariate Hermite interpolation polynomials and Gegenbauer polynomials. The paper discusses how to integrate the constructed C k basis functions into the assembling process of stiffness matrices and load vectors by using connectivity matrices. Numerical examples illustrate the applicability of the proposed concept. [ABSTRACT FROM AUTHOR]

The stability, robustness, accuracy, and efficiency of space-time finite element methods crucially depend on the choice of approximation spaces for test and trial functions. This is especially true for high-order, mixed finite element methods which often must satisfy an inf-sup condition in order to ensure stability. With this in mind, the primary objective of this paper and a companion paper is to provide a wide range of explicitly stated, conforming, finite element spaces in four dimensions. In this paper, we construct explicit high-order conforming finite elements on 4-cubes (tesseracts); our construction uses tools from the recently developed 'Finite Element Exterior Calculus'. With a focus on practical implementation, we provide details including Piola-type transformations, and explicit expressions for the volumetric, facet, face, edge, and vertex degrees of freedom. In addition, we establish important theoretical properties, such as the exactness of the finite element sequences, and the unisolvence of the degrees of freedom. [ABSTRACT FROM AUTHOR]

This paper employs a numerical approach to investigate the impact of liquidity shocks on European options in modeling markets. To accurately capture the behavior of European options under liquidity shocks, a coupled system of differential equations is employed, consisting of a degenerate parabolic equation and a diffusion-free equation. The primary focus is on developing and analyzing implicit-explicit methods for numerically simulating European option pricing, specifically considering the presence of liquidity shocks while ensuring the positivity of the solution. The paper also includes convergence analysis and establishes the discrete comparison principle for the developed methods. Numerical experiments are conducted using both uniform and nonuniform meshes to validate the theoretical findings, demonstrating the efficiency and accuracy of the proposed methods. [ABSTRACT FROM AUTHOR]

A novel higher-order multi-scale computational approach for efficient nonlocal gradient elasticity simulation of composite materials is proposed in this paper. Based on the use of decoupling technique, raw fourth-order nonlocal gradient equations with periodically oscillatory coefficients are decoupled as the new second-order multi-scale equations for reducing the continuity requirements when implementing multi-scale simulation. Next, a novel macro-micro coupled multi-scale computational model with higher-order corrected terms is rigorously devised for high-resolution multi-scale and nonlocal analysis of composite materials via multi-scale asymptotic analysis. Then, the local error analysis of higher-order multi-scale solutions in the point-wise sense illustrates that establishing higher-order asymptotic corrected terms plays a vital role in investigating the micromechanical mechanisms. Moreover, a rigorous global error estimation with an explicit rate of higher-order multi-scale solutions is first derived in the energy norm sense. In addition, an efficient multi-scale numerical algorithm is presented to effectively simulate nonlocal gradient elasticity problems of composite materials based on finite element method (FEM). And we also derive the convergence estimation of the proposed numerical algorithm. Finally, numerical experiments are conducted to gauge the efficiency and accuracy of the proposed higher-order multi-scale method. This paper offers a unified higher-order two-scale computational framework for enabling nonlocal gradient elasticity behavior simulation of composite materials. [ABSTRACT FROM AUTHOR]

The aim of this paper is to propose an Eulerian type of numerical finite difference approximation on structured grids for an advection dominated level set equation enhanced with a small curvature term. The proposed numerical scheme is based on upwind discretization in both the advection and the curvature term. This class of approximations is used in algebraic solvers, such as fast sweeping or fast marching methods, which respect the causality of solutions for which any information is propagating along characteristics, in order to solve the algebraic system in a finite number of computationally simple iterations or updates. An implicit numerical scheme for a time relaxed level set equation with non-zero right-hand side is used in this paper and solved by a nonlinear fast sweeping method. The results of the von Neumann stability analysis of the proposed scheme are presented, and numerical examples are used to verify the analysis. [ABSTRACT FROM AUTHOR]