Abstract The general family of Galerkin variational integrators are analyzed and a complete classification of such methods is proposed. This classification is based upon the type of basis function chosen to approximate the trajectories of material points and the numerical quadrature formula used in time. This approach leads to the definition of arbitrarily high order method in time. The proposed methodology is applied to the simulation of brownout phenomena occurring in helicopter take-off and landing.
{"title":"High-Order Variational Time Integrators for Particle Dynamics","authors":"E. Miglio, N. Parolini, M. Penati, R. Porcù","doi":"10.2478/caim-2018-0015","DOIUrl":"https://doi.org/10.2478/caim-2018-0015","url":null,"abstract":"Abstract The general family of Galerkin variational integrators are analyzed and a complete classification of such methods is proposed. This classification is based upon the type of basis function chosen to approximate the trajectories of material points and the numerical quadrature formula used in time. This approach leads to the definition of arbitrarily high order method in time. The proposed methodology is applied to the simulation of brownout phenomena occurring in helicopter take-off and landing.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"34 - 49"},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43817588","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract The political replacement effect is an interesting socio-political hypothesis introduced by Acemoglu and Robinson and statistically tested. It may determine, under some conditions, the phenomenon of innovation blocking, possibly leading to economic backwardness in a society. In a previous paper, we have introduced a kinetic model with stochastic evolutive game-type interactions, analyzing the relationship between the level of political competition in a society and the degree of economic liberalization. In the present paper we model we model the possibility of having a sort of phase transition occurring in the system when the phenomenon of blocking of the introduction of technological innovation, intended in a broad sense, appears. Crossing a critical point, the rules of interactions change by means of slightly different transition probabilities nevertheless determining very significant differences in the resulting long-term solutions.
{"title":"The Political Replacement Effect in a Kinetic Model of Social Dynamics with Phase Transition","authors":"M. Dolfin","doi":"10.2478/caim-2018-0013","DOIUrl":"https://doi.org/10.2478/caim-2018-0013","url":null,"abstract":"Abstract The political replacement effect is an interesting socio-political hypothesis introduced by Acemoglu and Robinson and statistically tested. It may determine, under some conditions, the phenomenon of innovation blocking, possibly leading to economic backwardness in a society. In a previous paper, we have introduced a kinetic model with stochastic evolutive game-type interactions, analyzing the relationship between the level of political competition in a society and the degree of economic liberalization. In the present paper we model we model the possibility of having a sort of phase transition occurring in the system when the phenomenon of blocking of the introduction of technological innovation, intended in a broad sense, appears. Crossing a critical point, the rules of interactions change by means of slightly different transition probabilities nevertheless determining very significant differences in the resulting long-term solutions.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"50 - 60"},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41682466","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract This paper aims at bridging Mathematical Modelling and Mathematics Education by studying the opinion dynamics of students who work in small groups during mathematics classrooms. In particular, we propose a model which hinges upon the pioneering work of Hegselmann and Krause on opinion dynamics and integrates recent results of interactionist research in Mathematical Education. More precisely, the proposed model incorporates the following features: 1) the feelings of each student towards the classmates (building upon the so-called I can" -you can" framework); 2) the different levels of preparation of the students; 3) the presence of the teacher, who may or may not intervene to drive the students towards the correct solution of the problem. Several numerical experiments are presented to assess the capability of the model in reproducing typical realistic scenarios.
{"title":"Student interactions during class activities: a mathematical model","authors":"D. Brunetto, C. Andrà, N. Parolini, M. Verani","doi":"10.2478/caim-2018-0011","DOIUrl":"https://doi.org/10.2478/caim-2018-0011","url":null,"abstract":"Abstract This paper aims at bridging Mathematical Modelling and Mathematics Education by studying the opinion dynamics of students who work in small groups during mathematics classrooms. In particular, we propose a model which hinges upon the pioneering work of Hegselmann and Krause on opinion dynamics and integrates recent results of interactionist research in Mathematical Education. More precisely, the proposed model incorporates the following features: 1) the feelings of each student towards the classmates (building upon the so-called I can\" -you can\" framework); 2) the different levels of preparation of the students; 3) the presence of the teacher, who may or may not intervene to drive the students towards the correct solution of the problem. Several numerical experiments are presented to assess the capability of the model in reproducing typical realistic scenarios.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"105 - 91"},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42360248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract Kinetic models have so far been used to model wealth distribution in a society. In particular, within the framework of the kinetic theory for active particles, some important models have been developed and proposed. They involve nonlinear interactions among individuals that are modeled according to game theoretical tools by introducing parameters governing the temporal dynamics of the system. In this present paper we propose an approach based on optimal control tools that aims to optimize this evolving dynamics from a social point of view. Namely, we look for time dependent control variables concerning the distribution of wealth that can be managed, for instance, by the government, in order to obtain a given social profile.
{"title":"On an optimal control strategy in a kinetic social dynamics model","authors":"D. Knopoff, G. Torres","doi":"10.2478/caim-2018-0014","DOIUrl":"https://doi.org/10.2478/caim-2018-0014","url":null,"abstract":"Abstract Kinetic models have so far been used to model wealth distribution in a society. In particular, within the framework of the kinetic theory for active particles, some important models have been developed and proposed. They involve nonlinear interactions among individuals that are modeled according to game theoretical tools by introducing parameters governing the temporal dynamics of the system. In this present paper we propose an approach based on optimal control tools that aims to optimize this evolving dynamics from a social point of view. Namely, we look for time dependent control variables concerning the distribution of wealth that can be managed, for instance, by the government, in order to obtain a given social profile.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"22 - 33"},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49035441","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract This paper deals with the modeling of interactions between the immune system and cancer cells, in the framework of the mathematical kinetic theory for active particles. The work deepens a previous paper of Belloquid et al. that assumes spatial homogeneity and discrete values of the activity of cancer and immune cells. A number of simulations are made with the aim to investigate how the state of the various cell populations evolves in time depending on the choice of the free parameters.
{"title":"A contribution to the mathematical modeling of immune-cancer competition","authors":"N. Dabnoun, M. Mongiovì","doi":"10.2478/caim-2018-0012","DOIUrl":"https://doi.org/10.2478/caim-2018-0012","url":null,"abstract":"Abstract This paper deals with the modeling of interactions between the immune system and cancer cells, in the framework of the mathematical kinetic theory for active particles. The work deepens a previous paper of Belloquid et al. that assumes spatial homogeneity and discrete values of the activity of cancer and immune cells. A number of simulations are made with the aim to investigate how the state of the various cell populations evolves in time depending on the choice of the free parameters.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"76 - 90"},"PeriodicalIF":1.3,"publicationDate":"2018-12-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44932155","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
S. Carillo, M. Chipot, V. Valente, G. V. Caffarelli
Abstract The existence and uniqueness of solution to a one-dimensional hyperbolic integro-differential problem arising in viscoelasticity is here considered. The kernel, in the linear viscoelasticity equation, represents the relaxation function which is characteristic of the considered material. Specifically, the case of a kernel, which does not satisfy the classical regularity requirements is analysed. This choice is suggested by applications according to the literature to model a wider variety of materials. A notable example of kernel, not satisfying the classical regularity requirements, is represented by a wedge continuous function. Indeed, the linear integro-differential viscoelasticity equation, characterised by a suitable wedge continuous relaxation function, is shown to give the classical linear wave equation via a limit procedure.
{"title":"On weak regularity requirements of the relaxation modulus in viscoelasticity","authors":"S. Carillo, M. Chipot, V. Valente, G. V. Caffarelli","doi":"10.2478/caim-2019-0014","DOIUrl":"https://doi.org/10.2478/caim-2019-0014","url":null,"abstract":"Abstract The existence and uniqueness of solution to a one-dimensional hyperbolic integro-differential problem arising in viscoelasticity is here considered. The kernel, in the linear viscoelasticity equation, represents the relaxation function which is characteristic of the considered material. Specifically, the case of a kernel, which does not satisfy the classical regularity requirements is analysed. This choice is suggested by applications according to the literature to model a wider variety of materials. A notable example of kernel, not satisfying the classical regularity requirements, is represented by a wedge continuous function. Indeed, the linear integro-differential viscoelasticity equation, characterised by a suitable wedge continuous relaxation function, is shown to give the classical linear wave equation via a limit procedure.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"10 1","pages":"78 - 87"},"PeriodicalIF":1.3,"publicationDate":"2018-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42837066","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We consider heat rectification in radial flows of turbulent helium II, where heat flux is not described by Fourier's law, but by a more general law. This is different from previous analyses of heat rectification, based on such law. In our simplified analysis we show that the coupling between heat flux and the gradient of vortex line density plays a decisive role in such rectification. Such rectification will be low at low and high values of the heat rate, but it may exhibit a very high value at an intermediate value of the heat rate. In particular, for a given range of values for the incoming heat ow, the outgoing heat flow corresponding to the exchange of internal and external temperatures would be very small. This would imply difficulties in heat removal in a given range of temperature gradients.
{"title":"Heat rectification in He II counterflow in radial geometries","authors":"L. Saluto, D. Jou, M. Mongiovì","doi":"10.2478/caim-2018-0017","DOIUrl":"https://doi.org/10.2478/caim-2018-0017","url":null,"abstract":"Abstract We consider heat rectification in radial flows of turbulent helium II, where heat flux is not described by Fourier's law, but by a more general law. This is different from previous analyses of heat rectification, based on such law. In our simplified analysis we show that the coupling between heat flux and the gradient of vortex line density plays a decisive role in such rectification. Such rectification will be low at low and high values of the heat rate, but it may exhibit a very high value at an intermediate value of the heat rate. In particular, for a given range of values for the incoming heat ow, the outgoing heat flow corresponding to the exchange of internal and external temperatures would be very small. This would imply difficulties in heat removal in a given range of temperature gradients.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"141 - 148"},"PeriodicalIF":1.3,"publicationDate":"2018-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43588964","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We identify sufficient conditions for the stability of some well-known finite difference schemes for the solution of the multivariable reaction-diffusion equations that model chemical reaction networks. Since the equations are mainly nonlinear, these conditions are obtained through local linearization. A recurrent condition is that the Jacobian matrix of the reaction part evaluated at some positive unknown solution is either D-semi-stable or semi-stable. We demonstrate that for a single reversible chemical reaction whose kinetics are monotone, the Jacobian matrix is D-semi-stable and therefore such schemes are guaranteed to work well.
{"title":"On the linear stability of some finite difference schemes for nonlinear reaction-diffusion models of chemical reaction networks","authors":"N. Muyinda, B. De Baets, Shodhan Rao","doi":"10.2478/caim-2018-0016","DOIUrl":"https://doi.org/10.2478/caim-2018-0016","url":null,"abstract":"Abstract We identify sufficient conditions for the stability of some well-known finite difference schemes for the solution of the multivariable reaction-diffusion equations that model chemical reaction networks. Since the equations are mainly nonlinear, these conditions are obtained through local linearization. A recurrent condition is that the Jacobian matrix of the reaction part evaluated at some positive unknown solution is either D-semi-stable or semi-stable. We demonstrate that for a single reversible chemical reaction whose kinetics are monotone, the Jacobian matrix is D-semi-stable and therefore such schemes are guaranteed to work well.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"121 - 140"},"PeriodicalIF":1.3,"publicationDate":"2018-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49146761","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract We study a stochastic Nonlinear Schrödinger Equation (NLSE), with additive white Gaussian noise, by means of the Nonlinear Fourier Transform (NFT). In particular, we focus on the propagation of discrete eigenvalues along a focusing fiber. Since the stochastic NLSE is not exactly integrable by means of the NFT, then we use a perturbation approach, where we assume that the signal-to-noise ratio is high. The zeroth-order perturbation leads to the deterministic NLSE while the first-order perturbation allows to describe the statistics of the discrete eigenvalues. This is important to understand the properties of the channel for recently devised optical transmission techniques, where the information is encoded in the nonlinear Fourier spectrum.
{"title":"Some results on discrete eigenvalues for the Stochastic Nonlinear Schrödinger Equation in fiber optics","authors":"L. Prati, L. Barletti","doi":"10.1515/caim-2015-0006","DOIUrl":"https://doi.org/10.1515/caim-2015-0006","url":null,"abstract":"Abstract We study a stochastic Nonlinear Schrödinger Equation (NLSE), with additive white Gaussian noise, by means of the Nonlinear Fourier Transform (NFT). In particular, we focus on the propagation of discrete eigenvalues along a focusing fiber. Since the stochastic NLSE is not exactly integrable by means of the NFT, then we use a perturbation approach, where we assume that the signal-to-noise ratio is high. The zeroth-order perturbation leads to the deterministic NLSE while the first-order perturbation allows to describe the statistics of the discrete eigenvalues. This is important to understand the properties of the channel for recently devised optical transmission techniques, where the information is encoded in the nonlinear Fourier spectrum.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"103 - 87"},"PeriodicalIF":1.3,"publicationDate":"2018-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48496010","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Abstract Segmentation is a typical task in image processing having as main goal the partitioning of the image into multiple segments in order to simplify its interpretation and analysis. One of the more popular segmentation model, formulated by Chan-Vese, is the piecewise constant Mumford-Shah model restricted to the case of two-phase segmentation. We consider a convex relaxation formulation of the segmentation model, that can be regarded as a nonsmooth optimization problem, because the presence of the l1-term. Two basic approaches in optimization can be distinguished to deal with its non differentiability: the smoothing methods and the nonsmoothing methods. In this work, a numerical comparison of some first order methods belongs of both approaches are presented. The relationships among the different methods are shown, and accuracy and efficiency tests are also performed on several images.
{"title":"Comparison of minimization methods for nonsmooth image segmentation","authors":"L. Antonelli, V. De Simone","doi":"10.1515/caim-2018-0005","DOIUrl":"https://doi.org/10.1515/caim-2018-0005","url":null,"abstract":"Abstract Segmentation is a typical task in image processing having as main goal the partitioning of the image into multiple segments in order to simplify its interpretation and analysis. One of the more popular segmentation model, formulated by Chan-Vese, is the piecewise constant Mumford-Shah model restricted to the case of two-phase segmentation. We consider a convex relaxation formulation of the segmentation model, that can be regarded as a nonsmooth optimization problem, because the presence of the l1-term. Two basic approaches in optimization can be distinguished to deal with its non differentiability: the smoothing methods and the nonsmoothing methods. In this work, a numerical comparison of some first order methods belongs of both approaches are presented. The relationships among the different methods are shown, and accuracy and efficiency tests are also performed on several images.","PeriodicalId":37903,"journal":{"name":"Communications in Applied and Industrial Mathematics","volume":"9 1","pages":"68 - 86"},"PeriodicalIF":1.3,"publicationDate":"2018-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"46495186","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}