Computational and Mathematical Methods最新文献

Mikhlin’s integral equation is a classical integral equation for solving boundary value problems for Laplace’s equation. The kernel of the integral equation is known as the Neumann kernel. Recently, an integral equation for solving the Riemann–Hilbert problem was derived. The kernel of the new integral equation is a generalization of the Neumann kernel, and hence, it is called the generalized Neumann kernel. The objective of this paper is to present a detailed comparison between these two integral equations with emphasis on their similarities and differences. This comparison is done through applying both equations to solve Laplace’s equation with Dirichlet boundary conditions in simply connected domains with smooth and piecewise smooth boundaries.

Applied complex network theory has become an interesting research field in the last years. Many papers have appeared on this subject dealing with the topological description of real transport systems, from small networks like the Italian airport network to the worldwide air transportation network. A comprehensive topological description of those critical structures is relevant in order to understand their dynamics, capacities, and vulnerabilities. In this work, for the first time, we describe the Spanish airport network (SAN) as a complex network. Nodes are airports, and links are flight connections weighed by traffic flow. We study its topological features and traffic dynamics. Our analysis shows that SAN has complex dynamics similar to small-size air transportation networks of other developed economies. It shares properties of small-world and scale-free networks, and it is highly connected and efficient and has a disassortative pattern for high-degree nodes.

The purpose of this article is to explore different types of solutions for the Kadomtsev-Petviashvili-modified Kadomtsev-Petviashvili (KP-mKP) equation which is termed as KP-Gardner equation, extensively used to model strong nonlinear internal waves in (1 + 2)-dimensions on the stratified ocean shelf. This evolution equation is also used to describe weakly nonlinear shallow-water wave and dispersive interracial waves traveling in a mildly rotating channel with slowly varying topography. Introducing Liu’s approach regarding the complete discrimination system for polynomial and the trial equation technique, a set of new solutions to the KP-mKP equation containing Jacobi elliptic function have been derived. It is found that these analytical solutions numerically exhibit different nonlinear structures such as solitary waves, shock waves, and periodic wave profiles. The reliability and effectiveness are confirmed from the numerical graphs of the solutions. Finally, the existence and validity of the various topological structures of the solutions are confirmed from the phase portrait of the dynamical system. Based on this investigation, it is confirmed that the method is not only suited for obtaining the classification of the solutions but also for qualitative analysis, which means that it can also be extended to other fields of application.

Dengue fever is a viral mosquito-borne disease, a significant global health concern, with more than one third of the world population at risk of acquiring the disease. Caused by 4 antigenically distinct but related virus serotypes, named DENV-1, DENV-2, DENV-3, and DENV-4, infection by one serotype confers lifelong immunity to that serotype and a short period of temporary cross immunity to other related serotypes. Severe dengue is epidemiologically associated with a secondary infection caused by a heterologous serotype via the so-called antibody-dependent enhancement (ADE), an immunological process enhancing a new infection. Within-host dengue modeling is restricted to a small number of studies so far. With many open questions, the understanding of immunopathogenesis of severe disease during recurrent infections is important to evaluate the impact of newly licensed vaccines. In this paper, we revisit the modeling framework proposed by Sebayang et al. and perform a detailed sensitivity analysis of the well-known biological parameters and its possible combinations to understand the existing data sets. Using numerical simulations, we investigate features of viral replication, antibody production, and infection clearance over time for three possible scenarios: primary infection, secondary infection caused by homologous serotype, and secondary infection caused by heterologous serotype. Besides, describing well the infection dynamics as reported in the immunology literature, our results provide information on parameter combinations to best describe the differences on the immunological dynamics of secondary infections with homologous and heterologous viruses. The results presented here will be used as baseline to investigate a more complex within-host dengue model.

The Box-Cox transformations are used to make the data more suitable for statistical analysis. We know from the literature that this transformation of the data can increase the rate of convergence of the tail of the distribution to the generalized extreme value distribution, and as a byproduct, the bias of the estimation procedure is reduced. The reduction of bias of the Hill estimator has been widely addressed in the literature of extreme value theory. Several techniques have been used to achieve such reduction of bias, either by removing the main component of the bias of the Hill estimator of the extreme value index (EVI) or by constructing new estimators based on generalized means or norms that generalize the Hill estimator. We are going to study the Box-Cox Hill estimator introduced by Teugels and Vanroelen, in 2004, proving the consistency and asymptotic normality of the estimator and addressing the choice and estimation of the power and shift parameters of the Box-Cox transformation for the EVI estimation. The performance of the estimators under study will be illustrated for finite samples through small-scale Monte Carlo simulation studies.

Summary. We investigate models to describe respiratory diseases with fast mutating virus pathogens such that after some years the aquired resistance is lost and hosts can be infected with new variants of the pathogen. Such models were initially suggested for respiartory diseases like influenza, showing complex dynamics in reasonable parameter regions when comparing to historic empirical influenza like illness data, e.g., from Ille de France. The seasonal forcing typical for respiratory diseases gives rise to the different rich dynamical scenarios with even small parameter changes. Especially the seasonality of the infection leads for small values already to period doubling bifurcations into chaos, besides additional coexisting attractors. Such models could in the future also play a role in understanding the presently experienced COVID-19 pandemic, under emerging new variants and with only limited vaccine efficacies against newly upcoming variants. From first period doubling bifurcations, we can eventually infer at which close by parameter regions complex dynamics including deterministic chaos can arise.

In this paper, we derive upper bounds that characterize the rate of convergence of the SOR method for solving a linear system of the form Gx = b, where G is a real symmetric positive semidefinite n × n matrix. The bounds are given in terms of the condition number of G, which is the ratio κ = α/β, where α is the largest eigenvalue of G and β is the smallest nonzero eigenvalue of G. Let H denote the related iteration matrix. Then, since G has a zero eigenvalue, the spectral radius of H equals 1, and the rate of convergence is determined by the size of η, the largest eigenvalue of H whose modulus differs from 1. The bound has the form |η|² ≤ 1 − 1/(κc), where c = 2 + log₂n. The main consequence from this bound is that small condition number forces fast convergence while large condition number allows slow convergence.

A mathematical model to describe the dynamic of a multiserotype infectious disease at the population level is studied. Applied to dengue fever epidemiology, we analyse a mathematical model with time delay to describe the cross-immunity protection period, including a key parameter for the antibody-dependent enhancement (ADE) effect, the well-known features of dengue fever infection. Numerical experiments are performed to show the stability of the coexistence equilibrium, which is completely determined by the basic reproduction number and by the invasion reproduction number, as well as the bifurcation structures for different scenarios of dengue fever transmission in a population. The model shows a rich dynamical behavior, from fixed points and periodic oscillations up to chaotic behaviour with complex attractors.

Mathematical models play an important role in epidemiology. The inclusion of a spatial component in epidemiological models is especially important to understand and address many relevant ecological and public health questions, e.g., when wanting to differentiate transmission patterns across geographical regions or when considering spatially heterogeneous intervention measures. However, the introduction of spatial effects can have significant consequences on the observed model dynamics and hence must be carefully analyzed and interpreted. Cellular automata epidemiological models typically rely on simplified computational grids but can provide valuable insight into the spatial dynamics of transmission within a population by suitably accounting for the connections between individuals in the considered community. In this paper, we describe a stochastic cellular automata disease model based on an extension of the traditional Susceptible-Infected-Recovered (SIR) compartmentalization of the population, namely, the Susceptible-Hospitalized-Asymptomatic-Recovered (SHAR) formulation, in which infected individuals either present a severe form of the disease, thus requiring hospitalization, or belong to the so-called mild/asymptomatic class. The critical transmission threshold is derived analytically in the nonspatial SHAR formulation, and this generalizes previously obtained theoretical results for the SIR model. We present simulation results discussing the effect of key model parameters and of spatial correlations on model outputs and propose an algorithm for tracking the evolution of infection clusters within the considered population. Focusing on the role of import and criticality on the overall dynamics, we conclude that the current spatial setting increases the critical transmission threshold in comparison to the nonspatial model.

A cell-centered finite volume semi-Lagrangian method is presented for the numerical solution of two-dimensional coupled Burgers’ problems on unstructured triangular meshes. The method combines a modified method of characteristics for the time integration and a cell-centered finite volume for the space discretization. The new method belongs to fractional-step algorithms for which the convection and the viscous parts in the coupled Burgers’ problems are treated separately. The crucial step of interpolation in the convection step is performed using two local procedures accounting for the element where the departure point is located. The resulting semidiscretized system is then solved using a third-order explicit Runge-Kutta scheme. In contrast to the Eulerian-based methods, we apply the new method for each time step along the characteristic curves instead of the time direction. The performance of the current method is verified using different examples for coupled Burgers’ problems with known analytical solutions. We also apply the method for simulation of an example of coupled Burgers’ flows in a complex geometry. In these test problems, the new cell-centered finite volume semi-Lagrangian method demonstrates its ability to accurately resolve the two-dimensional coupled Burgers’ problems.