Filipe A. dos Santos, A. Favata, Andrea Micheletti, Roberto Paroni, M. Picchi Scardaoni
The control of the shape of complex metasurfaces is a challenging task often addressed in the literature. This work presents a class of tessellated plates able to deform into surfaces of preprogrammed shape upon activation by any flexural load and that can be controlled by a single actuator. Quadric metasurfaces are obtained from infinitesimal origami maps of monohedral hexagonal tessellations of the plane, that is pavings in which all tiles are congruent to each other. Monohedral tessellated portions can be joined together to obtain more complex shapes, which can be locally synclastic or anticlastic and can have a certain roughness. We broaden previous work by providing a complete characterization of all the three known types of monohedral tessellations composed by irregular hexagons. The proposed two-dimensional structures may have applications in prosthetics, tissue engineering, wearable devices, energy harvesting devices, tunable focus mirrors and adaptive facades. The study is divided in two parts. In Part I, after introducing the discrete kinematics of tessellated plates, it is proved analytically that essentially each type of monohedral hexagonal tessellation possesses only one deformation mode. Afterwards, several numerical examples are provided to demonstrate the variety of achievable surface shapes. In Part II, first the metasurfaces corresponding to assigned tile geometries are given a continuum description, which establishes that the continuous interpolant is always a quadric. Then, experimental results on fabrication, assembly and surface accuracy are reported.
{"title":"Programming quadric metasurfaces via infinitesimal origami maps of monohedral hexagonal tessellations: Part I","authors":"Filipe A. dos Santos, A. Favata, Andrea Micheletti, Roberto Paroni, M. Picchi Scardaoni","doi":"10.1098/rspa.2023.0430","DOIUrl":"https://doi.org/10.1098/rspa.2023.0430","url":null,"abstract":"The control of the shape of complex metasurfaces is a challenging task often addressed in the literature. This work presents a class of tessellated plates able to deform into surfaces of preprogrammed shape upon activation by any flexural load and that can be controlled by a single actuator. Quadric metasurfaces are obtained from infinitesimal origami maps of monohedral hexagonal tessellations of the plane, that is pavings in which all tiles are congruent to each other. Monohedral tessellated portions can be joined together to obtain more complex shapes, which can be locally synclastic or anticlastic and can have a certain roughness. We broaden previous work by providing a complete characterization of all the three known types of monohedral tessellations composed by irregular hexagons. The proposed two-dimensional structures may have applications in prosthetics, tissue engineering, wearable devices, energy harvesting devices, tunable focus mirrors and adaptive facades. The study is divided in two parts. In Part I, after introducing the discrete kinematics of tessellated plates, it is proved analytically that essentially each type of monohedral hexagonal tessellation possesses only one deformation mode. Afterwards, several numerical examples are provided to demonstrate the variety of achievable surface shapes. In Part II, first the metasurfaces corresponding to assigned tile geometries are given a continuum description, which establishes that the continuous interpolant is always a quadric. Then, experimental results on fabrication, assembly and surface accuracy are reported.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139878410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ligand-based virtual screening aims to reduce the cost and duration of small molecule drug discovery campaigns. Shape similarity can be used to screen large databases, with the goal of predicting potential new hits by comparing with molecules with known favourable properties. This paper presents the theory underpinning RGMolSA, a new alignment-free and mesh-free surface-based molecular shape descriptor derived from the mathematical theory of Riemannian geometry. The treatment of a molecule as a series of intersecting spheres allows the description of its surface geometry using the Riemannian metric , obtained by considering the spectrum of the Laplacian. This gives a simple vector descriptor constructed of the weighted surface area and eight non-zero eigenvalues, which capture the surface shape. We demonstrate the potential of our method by considering a series of PDE5 inhibitors that are known to have similar shape as an initial test case. RGMolSA displays promise when compared with existing shape descriptors and in its capability to handle different molecular conformers.
{"title":"Riemannian geometry and molecular similarity I: spectrum of the Laplacian","authors":"Stuart J. Hall, Rachael Pirie, Daniel J. Cole","doi":"10.1098/rspa.2023.0343","DOIUrl":"https://doi.org/10.1098/rspa.2023.0343","url":null,"abstract":"\u0000 Ligand-based virtual screening aims to reduce the cost and duration of small molecule drug discovery campaigns. Shape similarity can be used to screen large databases, with the goal of predicting potential new hits by comparing with molecules with known favourable properties. This paper presents the theory underpinning RGMolSA, a new alignment-free and mesh-free surface-based molecular shape descriptor derived from the mathematical theory of Riemannian geometry. The treatment of a molecule as a series of intersecting spheres allows the description of its surface geometry using the\u0000 Riemannian metric\u0000 , obtained by considering the spectrum of the Laplacian. This gives a simple vector descriptor constructed of the weighted surface area and eight non-zero eigenvalues, which capture the surface shape. We demonstrate the potential of our method by considering a series of PDE5 inhibitors that are known to have similar shape as an initial test case. RGMolSA displays promise when compared with existing shape descriptors and in its capability to handle different molecular conformers.\u0000","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-02-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139966431","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Josep Fontana-McNally, Eva Miranda, Daniel Peralta-Salas
We prove that the correspondence between Reeb and Beltrami vector fields presented in Etnyre & Ghrist (Etnyre, Ghrist 2000 Nonlinearity13 , 441–458 ( doi:10.1088/0951-7715/13/2/306 )) can be made equivariant whenever additional symmetries of the underlying geometric structures are considered. As a corollary of this correspondence, we show that energy levels above the maximum of the potential energy of mechanical Hamiltonian systems can be viewed as stationary fluid flows, though the metric is not prescribed. In particular, we showcase the emblematic example of the n -body problem and focus on the Kepler problem. We explicitly construct a compatible Riemannian metric that makes the Kepler problem of celestial mechanics a stationary fluid flow (of Beltrami type) on a suitable manifold, the Kepler–Euler flow .
{"title":"An equivariant Reeb–Beltrami correspondence and the Kepler–Euler flow","authors":"Josep Fontana-McNally, Eva Miranda, Daniel Peralta-Salas","doi":"10.1098/rspa.2023.0499","DOIUrl":"https://doi.org/10.1098/rspa.2023.0499","url":null,"abstract":"We prove that the correspondence between Reeb and Beltrami vector fields presented in Etnyre & Ghrist (Etnyre, Ghrist 2000 <jats:italic>Nonlinearity</jats:italic> <jats:bold>13</jats:bold> , 441–458 ( <jats:ext-link xmlns:xlink=\"http://www.w3.org/1999/xlink\" ext-link-type=\"uri\" xlink:href=\"http://dx.doi.org/10.1088/0951-7715/13/2/306\">doi:10.1088/0951-7715/13/2/306</jats:ext-link> )) can be made equivariant whenever additional symmetries of the underlying geometric structures are considered. As a corollary of this correspondence, we show that energy levels above the maximum of the potential energy of mechanical Hamiltonian systems can be viewed as stationary fluid flows, though the metric is not prescribed. In particular, we showcase the emblematic example of the <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:mi>n</mml:mi> </mml:math> </jats:inline-formula> -body problem and focus on the Kepler problem. We explicitly construct a compatible Riemannian metric that makes the Kepler problem of celestial mechanics a stationary fluid flow (of Beltrami type) on a suitable manifold, the <jats:italic>Kepler–Euler flow</jats:italic> .","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025113","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Motivated by observations of loop formation and elimination phenomena in elastic filaments subjected to torsion and axial end displacement, we develop a tension–torsion tester to study the slack–extension responses of filaments with varied initial twists. The experiments are conducted by initially twisting the filament by a specific degree and subsequently adjusting the axial end displacement. By continuously monitoring the correlation between torque and filament configuration, we can accurately determine the critical points associated with buckling, loop formation and loop elimination. The interconversion of link, twist and writhe is also tracked in the experiment. We establish a relation between torque and end displacement of filament that provides an insight into the conditions leading to instabilities. Three buckling criteria are also evaluated in the context of measurement data. Simulations on the slack and extension of the twisted filaments based on the Cosserat rod theory are performed, providing reliable predictions of the configuration evolution. Our results highlight the importance of twist and slack in torsional buckling processes, providing guidance for the control of systems with twisted slender structures.
{"title":"Toward a further understanding of the loop formation and elimination in twisted filament: experiments and validation","authors":"Jiongjiong Hu, Jiahui Teng, Lei Liu, Dabiao Liu","doi":"10.1098/rspa.2023.0719","DOIUrl":"https://doi.org/10.1098/rspa.2023.0719","url":null,"abstract":"Motivated by observations of loop formation and elimination phenomena in elastic filaments subjected to torsion and axial end displacement, we develop a tension–torsion tester to study the slack–extension responses of filaments with varied initial twists. The experiments are conducted by initially twisting the filament by a specific degree and subsequently adjusting the axial end displacement. By continuously monitoring the correlation between torque and filament configuration, we can accurately determine the critical points associated with buckling, loop formation and loop elimination. The interconversion of link, twist and writhe is also tracked in the experiment. We establish a relation between torque and end displacement of filament that provides an insight into the conditions leading to instabilities. Three buckling criteria are also evaluated in the context of measurement data. Simulations on the slack and extension of the twisted filaments based on the Cosserat rod theory are performed, providing reliable predictions of the configuration evolution. Our results highlight the importance of twist and slack in torsional buckling processes, providing guidance for the control of systems with twisted slender structures.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025192","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
A novel methodology for determining the singular point of an enveloped surface is put forward. Unlike some existing methods, the presented method starts directly from the equation of the enveloped surface instead of that of the generating surface, and it is thus called a direct method. The calculation for the normal vector of the enveloped surface is well simplified with the help of the moving frame approach, which makes the presented method feasible. The singularity condition equation is extracted by using the theory of linear algebra. For singular points with different properties, proper solving techniques are established, including resultant elimination and simple elimination. Applying the developed method, the undercutting characteristics of the Archimedes worm wheel are investigated from the perspective of spatial meshing. The numerical results demonstrate that the worm wheel generally has one undercutting limit line, whose trend is along the tooth width of the wheel. Locating on one side of the tooth surface and near the tooth root is a dangerous part of the worm wheel undercutting. The proposed method is beneficial for the development of gear meshing science.
{"title":"Direct method to determine singular point of enveloped surface and its application to worm wheel tooth surface","authors":"Jian Cui, Yaping Zhao, Qingxiang Meng, Gongfa Li","doi":"10.1098/rspa.2023.0369","DOIUrl":"https://doi.org/10.1098/rspa.2023.0369","url":null,"abstract":"A novel methodology for determining the singular point of an enveloped surface is put forward. Unlike some existing methods, the presented method starts directly from the equation of the enveloped surface instead of that of the generating surface, and it is thus called a direct method. The calculation for the normal vector of the enveloped surface is well simplified with the help of the moving frame approach, which makes the presented method feasible. The singularity condition equation is extracted by using the theory of linear algebra. For singular points with different properties, proper solving techniques are established, including resultant elimination and simple elimination. Applying the developed method, the undercutting characteristics of the Archimedes worm wheel are investigated from the perspective of spatial meshing. The numerical results demonstrate that the worm wheel generally has one undercutting limit line, whose trend is along the tooth width of the wheel. Locating on one side of the tooth surface and near the tooth root is a dangerous part of the worm wheel undercutting. The proposed method is beneficial for the development of gear meshing science.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025397","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The buckling instability of cylindrical shells under axial compression has been one of the most renowned problems in structural engineering for several decades. Many pioneering works in the twentieth century have provided insights into understanding the shells’ infamous imperfection sensitivity and led to reliability-based designs. However, a recent surge in numerical studies of the snaking phenomenon explores the development of a localized stable post-buckling mode in axially compressed cylindrical shells. Hitherto, none of the experimental studies report on the evolution of azimuthal snaking. In this work, experimental studies are carried out with the objective of revealing the snaking phenomenon. The axial compression experiments are performed on 3D-printed shells made of thermoplastic polyurethane (TPU). The work’s novelty lies in the usage of TPU shells for slowing down the propagation of circumferential dimples and making it feasible to capture them using photogrammetry. Despite the match between the experimental and numerical mode shapes, the experiments reveal multiple routes for the snaking sequence. Furthermore, mode transitions such as reduction in circumferential wave number and transformation of symmetric mode into an asymmetric one are observed. These experimental results provide insights into the localized phenomenon of snaking and validate numerical solutions.
几十年来,圆柱形壳体在轴向压缩下的屈曲不稳定性一直是结构工程领域最著名的问题之一。二十世纪的许多开创性工作为了解壳体的缺陷敏感性提供了见解,并促成了基于可靠性的设计。然而,最近对蛇形现象的数值研究激增,探索了轴向压缩圆柱壳的局部稳定后屈曲模式的发展。迄今为止,还没有任何一项实验研究报告了方位蛇行的演变过程。在这项工作中,为了揭示蛇行现象,我们进行了实验研究。轴向压缩实验是在热塑性聚氨酯(TPU)制成的 3D 打印外壳上进行的。这项工作的新颖之处在于使用热塑性聚氨酯外壳来减缓圆周凹痕的传播,并使使用摄影测量捕捉凹痕成为可行。尽管实验和数值模式形状相吻合,但实验揭示了蛇形序列的多种路径。此外,还观察到模式转换,如圆周波数减少以及对称模式转变为不对称模式。这些实验结果提供了对局部蛇行现象的见解,并验证了数值解决方案。
{"title":"Experimental studies on snaking in 3D-printed cylindrical shells under axial compression using photogrammetry","authors":"V. Ravulapalli, G. Raju, V. Narayanamurthy","doi":"10.1098/rspa.2023.0631","DOIUrl":"https://doi.org/10.1098/rspa.2023.0631","url":null,"abstract":"The buckling instability of cylindrical shells under axial compression has been one of the most renowned problems in structural engineering for several decades. Many pioneering works in the twentieth century have provided insights into understanding the shells’ infamous imperfection sensitivity and led to reliability-based designs. However, a recent surge in numerical studies of the snaking phenomenon explores the development of a localized stable post-buckling mode in axially compressed cylindrical shells. Hitherto, none of the experimental studies report on the evolution of azimuthal snaking. In this work, experimental studies are carried out with the objective of revealing the snaking phenomenon. The axial compression experiments are performed on 3D-printed shells made of thermoplastic polyurethane (TPU). The work’s novelty lies in the usage of TPU shells for slowing down the propagation of circumferential dimples and making it feasible to capture them using photogrammetry. Despite the match between the experimental and numerical mode shapes, the experiments reveal multiple routes for the snaking sequence. Furthermore, mode transitions such as reduction in circumferential wave number and transformation of symmetric mode into an asymmetric one are observed. These experimental results provide insights into the localized phenomenon of snaking and validate numerical solutions.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025034","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
M. Harrison, M. Riva, M. Mousavi Nezhad, A. Guadagnini
We derive analytical formulations relating the spatial covariance ( CY ) of (log-transformed) hydraulic conductivities to auto- and cross-covariances of porosity ( ϕ ) and representative soil particle sizes within the framework of the classical Terzaghi model. The latter provides an empirical relationship which is widely used to obtain conductivity estimates. We frame the study within recent stochastic approaches and conceptualize appropriate transformations of ϕ and representative soil particle size as Generalized Sub-Gaussian (GSG) spatially cross-correlated random processes. Consistency of the theoretical framework against sample distributions of ϕ and particle size is assessed through the analysis of field data. A perturbation-based approach yields workable expressions of CY upon truncating the otherwise exact analytical solution at given orders of approximations. Our analytical (truncated) log-conductivity covariance is in agreement with its Monte Carlo-based counterpart. A Global Sensitivity Analysis relying on classical Sobol indices quantifies the relative importance of all parameters embedded in the formulation of CY . We show that parameters driving the GSG nature of the distribution of (transformed) porosity are key to the main features of CY . We also document the relevance of properly capturing emergences of possible cross-correlations between ϕ and representative particle size to reconstruct conductivity fields.
在经典特尔扎吉模型的框架内,我们推导出了(对数变换)水力传导性的空间协方差(C Y)与孔隙度(j)和代表性土壤颗粒尺寸的自协方差和交叉协方差的相关分析公式。后者提供了一种经验关系,被广泛用于估算导流系数。我们将这项研究纳入最新的随机方法,并将 ϕ 和代表性土壤粒径的适当变换概念化为广义子高斯(GSG)空间交叉相关随机过程。通过分析实地数据,评估了理论框架与 ϕ 和粒径样本分布的一致性。采用基于扰动的方法,在给定近似阶数下截断原本精确的分析解,即可得到可行的 C Y 表达式。我们的分析(截断)对数电导协方差与基于蒙特卡洛的协方差一致。全局敏感性分析依赖于经典的索布尔指数,量化了 C Y 公式中所有参数的相对重要性。我们表明,驱动(转换)孔隙度分布 GSG 性质的参数是 C Y 主要特征的关键。我们还证明了正确捕捉 ϕ 与代表性粒度之间可能存在的交叉相关性对重建电导场的重要性。
{"title":"Estimation of auto-covariance of log hydraulic conductivity from Generalized Sub-Gaussian porosity and particle size random fields","authors":"M. Harrison, M. Riva, M. Mousavi Nezhad, A. Guadagnini","doi":"10.1098/rspa.2023.0476","DOIUrl":"https://doi.org/10.1098/rspa.2023.0476","url":null,"abstract":"We derive analytical formulations relating the spatial covariance ( <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mi>C</mml:mi> <mml:mi>Y</mml:mi> </mml:msub> </mml:math> </jats:inline-formula> ) of (log-transformed) hydraulic conductivities to auto- and cross-covariances of porosity ( <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:mi>ϕ</mml:mi> </mml:math> </jats:inline-formula> ) and representative soil particle sizes within the framework of the classical Terzaghi model. The latter provides an empirical relationship which is widely used to obtain conductivity estimates. We frame the study within recent stochastic approaches and conceptualize appropriate transformations of <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:mi>ϕ</mml:mi> </mml:math> </jats:inline-formula> and representative soil particle size as Generalized Sub-Gaussian (GSG) spatially cross-correlated random processes. Consistency of the theoretical framework against sample distributions of <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:mi>ϕ</mml:mi> </mml:math> </jats:inline-formula> and particle size is assessed through the analysis of field data. A perturbation-based approach yields workable expressions of <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mi>C</mml:mi> <mml:mi>Y</mml:mi> </mml:msub> </mml:math> </jats:inline-formula> upon truncating the otherwise exact analytical solution at given orders of approximations. Our analytical (truncated) log-conductivity covariance is in agreement with its Monte Carlo-based counterpart. A Global Sensitivity Analysis relying on classical Sobol indices quantifies the relative importance of all parameters embedded in the formulation of <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mi>C</mml:mi> <mml:mi>Y</mml:mi> </mml:msub> </mml:math> </jats:inline-formula> . We show that parameters driving the GSG nature of the distribution of (transformed) porosity are key to the main features of <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:msub> <mml:mi>C</mml:mi> <mml:mi>Y</mml:mi> </mml:msub> </mml:math> </jats:inline-formula> . We also document the relevance of properly capturing emergences of possible cross-correlations between <jats:inline-formula> <mml:math xmlns:mml=\"http://www.w3.org/1998/Math/MathML\"> <mml:mi>ϕ</mml:mi> </mml:math> </jats:inline-formula> and representative particle size to reconstruct conductivity fields.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025114","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The problem of thermal convection in a layer of viscous incompressible fluid is analysed. The heat flux law is taken to be one of Cattaneo type. The time derivative of the heat flux is allowed to be a material derivative, or a general objective derivative. It is shown that only one objective derivative leads to results consistent with what one expects in real life. This objective derivative leads to a Cattaneo–Christov theory, and the results for linear instability theory are in agreement with those for a material derivative. It is further shown that none of the theories allow a standard nonlinear, energy stability analysis. A further heat flux due to P.M. Mariano is added and then an analysis is performed for stationary convection, oscillatory convection, and fully nonlinear theory. For the material derivative case, the analysis proceeds and global nonlinear stability is achieved. For Cattaneo–Christov theory, it appears necessary to add a regularization term in the equation for the heat flux, and even then the analysis only works in two space dimensions, and is conditional upon the size of the initial data. For the three-dimensional situation, it is shown how a nonlinear stability analysis may be achieved with a Navier–Stokes–Voigt fluid rather than a Navier–Stokes one.
{"title":"Thermal convection with a Cattaneo heat flux model","authors":"M. Gentile, B. Straughan","doi":"10.1098/rspa.2023.0771","DOIUrl":"https://doi.org/10.1098/rspa.2023.0771","url":null,"abstract":"The problem of thermal convection in a layer of viscous incompressible fluid is analysed. The heat flux law is taken to be one of Cattaneo type. The time derivative of the heat flux is allowed to be a material derivative, or a general objective derivative. It is shown that only one objective derivative leads to results consistent with what one expects in real life. This objective derivative leads to a Cattaneo–Christov theory, and the results for linear instability theory are in agreement with those for a material derivative. It is further shown that none of the theories allow a standard nonlinear, energy stability analysis. A further heat flux due to P.M. Mariano is added and then an analysis is performed for stationary convection, oscillatory convection, and fully nonlinear theory. For the material derivative case, the analysis proceeds and global nonlinear stability is achieved. For Cattaneo–Christov theory, it appears necessary to add a regularization term in the equation for the heat flux, and even then the analysis only works in two space dimensions, and is conditional upon the size of the initial data. For the three-dimensional situation, it is shown how a nonlinear stability analysis may be achieved with a Navier–Stokes–Voigt fluid rather than a Navier–Stokes one.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
I. I. Mullyadzhanov, A. S. Gudko, R. I. Mullyadzhanov, A. A. Gelash
We consider the model of the focusing one-dimensional nonlinear Schrödinger equation (fNLSE) in the presence of an unstable constant background, which exhibits coherent solitary wave structures—breathers. Within the inverse scattering transform (IST) method, we study the problem of the scattering data numerical computation for a broad class of breathers localized in space. Such a direct scattering transform (DST) procedure requires a numerical solution of the auxiliary Zakharov–Shabat system with boundary conditions corresponding to the background. To find the solution, we compute the transfer matrix using the second-order Boffetta–Osborne approach and recently developed high-order numerical schemes based on the Magnus expansion. To recover the scattering data of breathers, we derive analytical relations between the scattering coefficients and the transfer matrix elements. Then we construct localized single- and multi-breather solutions and verify the developed numerical approach by recovering the complete set of scattering data with the built-in accuracy providing the information about the amplitude, velocity, phase and position of each breather. To combine the conventional IST approach with the efficient dressing method for multi-breather solutions, we derive the exact relation between the parameters of breathers in these two frameworks.
我们考虑了存在不稳定恒定背景的聚焦一维非线性薛定谔方程(fNLSE)模型,该模型表现出相干孤波结构--呼吸器。在反向散射变换(IST)方法中,我们研究了在空间中定位的一大类呼吸器的散射数据数值计算问题。这种直接散射变换(DST)程序需要对带有与背景相对应的边界条件的辅助 Zakharov-Shabat 系统进行数值求解。为了求解,我们使用二阶 Boffetta-Osborne 方法和最近开发的基于马格努斯展开的高阶数值方案计算传递矩阵。为了恢复呼吸器的散射数据,我们推导出了散射系数与传递矩阵元素之间的分析关系。然后,我们构建了局部单呼吸器和多呼吸器解决方案,并通过提供每个呼吸器的振幅、速度、相位和位置信息的内置精度恢复整套散射数据来验证所开发的数值方法。为了将传统的 IST 方法与多呼吸器解的高效修整方法相结合,我们推导出了这两种框架中呼吸器参数之间的精确关系。
{"title":"Numerical direct scattering transform for breathers","authors":"I. I. Mullyadzhanov, A. S. Gudko, R. I. Mullyadzhanov, A. A. Gelash","doi":"10.1098/rspa.2023.0529","DOIUrl":"https://doi.org/10.1098/rspa.2023.0529","url":null,"abstract":"We consider the model of the focusing one-dimensional nonlinear Schrödinger equation (fNLSE) in the presence of an unstable constant background, which exhibits coherent solitary wave structures—breathers. Within the inverse scattering transform (IST) method, we study the problem of the scattering data numerical computation for a broad class of breathers localized in space. Such a direct scattering transform (DST) procedure requires a numerical solution of the auxiliary Zakharov–Shabat system with boundary conditions corresponding to the background. To find the solution, we compute the transfer matrix using the second-order Boffetta–Osborne approach and recently developed high-order numerical schemes based on the Magnus expansion. To recover the scattering data of breathers, we derive analytical relations between the scattering coefficients and the transfer matrix elements. Then we construct localized single- and multi-breather solutions and verify the developed numerical approach by recovering the complete set of scattering data with the built-in accuracy providing the information about the amplitude, velocity, phase and position of each breather. To combine the conventional IST approach with the efficient dressing method for multi-breather solutions, we derive the exact relation between the parameters of breathers in these two frameworks.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025130","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The number of non-negative integer matrices with given row and column sums features in a variety of problems in mathematics and statistics but no closed-form expression for it is known, so we rely on approximations. In this paper, we describe a new such approximation, motivated by consideration of the statistics of matrices with non-integer numbers of columns. This estimate can be evaluated in time linear in the size of the matrix and returns results of accuracy as good as or better than existing linear-time approximations across a wide range of settings. We show that the estimate is asymptotically exact in the regime of sparse tables, while empirically performing at least as well as other linear-time estimates in the regime of dense tables. We also use the new estimate as the starting point for an improved numerical method for either counting or sampling matrices with given margins using sequential importance sampling. Code implementing our methods is available.
{"title":"Improved estimates for the number of non-negative integer matrices with given row and column sums","authors":"Maximilian Jerdee, Alec Kirkley, M. E. J. Newman","doi":"10.1098/rspa.2023.0470","DOIUrl":"https://doi.org/10.1098/rspa.2023.0470","url":null,"abstract":"The number of non-negative integer matrices with given row and column sums features in a variety of problems in mathematics and statistics but no closed-form expression for it is known, so we rely on approximations. In this paper, we describe a new such approximation, motivated by consideration of the statistics of matrices with non-integer numbers of columns. This estimate can be evaluated in time linear in the size of the matrix and returns results of accuracy as good as or better than existing linear-time approximations across a wide range of settings. We show that the estimate is asymptotically exact in the regime of sparse tables, while empirically performing at least as well as other linear-time estimates in the regime of dense tables. We also use the new estimate as the starting point for an improved numerical method for either counting or sampling matrices with given margins using sequential importance sampling. Code implementing our methods is available.","PeriodicalId":20716,"journal":{"name":"Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences","volume":null,"pages":null},"PeriodicalIF":3.5,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140025229","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"综合性期刊","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}