Journal of Scientific Computing最新文献_第9页

Structure Preserving Polytopal Discontinuous Galerkin Methods for the Numerical Modeling of Neurodegenerative Diseases 用于神经退行性疾病数值建模的结构保持型多顶点非连续伽勒金方法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-20 DOI: 10.1007/s10915-024-02581-7

Mattia Corti, Francesca Bonizzoni, Paola F. Antonietti

{"title":"Structure Preserving Polytopal Discontinuous Galerkin Methods for the Numerical Modeling of Neurodegenerative Diseases","authors":"Mattia Corti, Francesca Bonizzoni, Paola F. Antonietti","doi":"10.1007/s10915-024-02581-7","DOIUrl":"https://doi.org/10.1007/s10915-024-02581-7","url":null,"abstract":"Many neurodegenerative diseases are connected to the spreading of misfolded prionic proteins. In this paper, we analyse the process of misfolding and spreading of both (alpha )-synuclein and Amyloid-(beta ), related to Parkinson’s and Alzheimer’s diseases, respectively. We introduce and analyze a positivity-preserving numerical method for the discretization of the Fisher-Kolmogorov equation, modelling accumulation and spreading of prionic proteins. The proposed approximation method is based on the discontinuous Galerkin method on polygonal and polyhedral grids for space discretization and on (vartheta -)method time integration scheme. We prove the existence of the discrete solution and a convergence result where the Implicit Euler scheme is employed for time integration. We show that the proposed approach is structure-preserving, in the sense that it guarantees that the discrete solution is non-negative, a feature that is of paramount importance in practical application. The numerical verification of our numerical model is performed both using a manufactured solution and considering wavefront propagation in two-dimensional polygonal grids. Next, we present a simulation of (alpha )-synuclein spreading in a two-dimensional brain slice in the sagittal plane. The polygonal mesh for this simulation is agglomerated maintaining the distinction of white and grey matter, taking advantage of the flexibility of PolyDG methods in the mesh construction. Finally, we simulate the spreading of Amyloid-(beta ) in a patient-specific setting by using a three-dimensional geometry reconstructed from magnetic resonance images and an initial condition reconstructed from positron emission tomography. Our numerical simulations confirm that the proposed method is able to capture the evolution of Parkinson’s and Alzheimer’s diseases.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141529168","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A Hybrid SBP-SAT/Fourier Pseudo-spectral Method for the Transient Wigner Equation Involving Inflow Boundary Conditions 涉及流入边界条件的瞬态维格纳方程的 SBP-SAT/Fourier 伪谱混合方法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-20 DOI: 10.1007/s10915-024-02582-6

Zhangpeng Sun, Wenqi Yao, Qiuping Yu

{"title":"A Hybrid SBP-SAT/Fourier Pseudo-spectral Method for the Transient Wigner Equation Involving Inflow Boundary Conditions","authors":"Zhangpeng Sun, Wenqi Yao, Qiuping Yu","doi":"10.1007/s10915-024-02582-6","DOIUrl":"https://doi.org/10.1007/s10915-024-02582-6","url":null,"abstract":"In this paper, a hybrid SBP-SAT/pseudo-spectral method is proposed for solving the time-dependent Wigner equation. High-order summation-by-parts (SBP) operators are utilized to discretize the Wigner equation spatially, where the inflow boundary conditions are weakly imposed by adding simultaneous approximation terms (SATs) to the semi-discretized Wigner equation. The pseudo-differential term, governing the quantum effect, is discretized in a pseudo-spectral manner with spectral accuracy. (L^2)-stabilities of both the semi-discretized (excluding time) and fully discretized systems are thoroughly discussed, with the inclusion of an arbitrary-stage explicit Runge–Kutta scheme for time integration. Numerical experiments are conducted, including simulations of a harmonic oscillator, a Gaussian wave packet, and a typical RTD with its I–V characteristic curves. The numerical results demonstrate: (1) the accuracy order of the numerical scheme in discretizing the Wigner equation in phase space matches the theoretical value; (2) observation of typical quantum effects, including tunneling and negative resistance; and (3) rapid convergence of numerical solutions relative to the accuracy order of SBP operators.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512490","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Randomized Methods for Computing Optimal Transport Without Regularization and Their Convergence Analysis 计算无正规化最优传输的随机方法及其收敛性分析

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-20 DOI: 10.1007/s10915-024-02570-w

Yue Xie, Zhongjian Wang, Zhiwen Zhang

{"title":"Randomized Methods for Computing Optimal Transport Without Regularization and Their Convergence Analysis","authors":"Yue Xie, Zhongjian Wang, Zhiwen Zhang","doi":"10.1007/s10915-024-02570-w","DOIUrl":"https://doi.org/10.1007/s10915-024-02570-w","url":null,"abstract":"The optimal transport (OT) problem can be reduced to a linear programming (LP) problem through discretization. In this paper, we introduced the random block coordinate descent (RBCD) methods to directly solve this LP problem. Our approach involves restricting the potentially large-scale optimization problem to small LP subproblems constructed via randomly chosen working sets. By using a random Gauss-Southwell-q rule to select these working sets, we equip the vanilla version of (({textbf {RBCD}}_0)) with almost sure convergence and a linear convergence rate to solve general standard LP problems. To further improve the efficiency of the (({textbf {RBCD}}_0)) method, we explore the special structure of constraints in the OT problems and leverage the theory of linear systems to propose several approaches for refining the random working set selection and accelerating the vanilla method. Inexact versions of the RBCD methods are also discussed. Our preliminary numerical experiments demonstrate that the accelerated random block coordinate descent (ARBCD) method compares well with other solvers including stabilized Sinkhorn’s algorithm when seeking solutions with relatively high accuracy, and offers the advantage of saving memory.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512489","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

On Approximate Matrix Factorization and TASE W-Methods for the Time Integration of Parabolic Partial Differential Equations 关于抛物型偏微分方程时间积分的近似矩阵因式分解和 TASE W 方法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-19 DOI: 10.1007/s10915-024-02579-1

Dajana Conte, Severiano González-Pinto, Domingo Hernández-Abreu, Giovanni Pagano

{"title":"On Approximate Matrix Factorization and TASE W-Methods for the Time Integration of Parabolic Partial Differential Equations","authors":"Dajana Conte, Severiano González-Pinto, Domingo Hernández-Abreu, Giovanni Pagano","doi":"10.1007/s10915-024-02579-1","DOIUrl":"https://doi.org/10.1007/s10915-024-02579-1","url":null,"abstract":"Linearly implicit methods for Ordinary Differential Equations combined with the application of Approximate Matrix Factorization (AMF) provide efficient numerical methods for the solution of large semi-discrete parabolic Partial Differential Equations in several spatial dimensions. Interesting particular subclasses of such linearly implicit methods are the so-called W-methods and the TASE W-methods recently introduced in González-Pinto et al. (Appl Numer Math, 188:129–145, 2023) with the aim of reducing the computational cost of the TASE Runge–Kutta methods in Bassenne et al. (J Comput Phys 424:109847, 2021) and Calvo et al. (J Comp Phys 436:110316, 2021). In this paper, we study the application of the AMF approach in combination with TASE W-methods. While for AMF W-methods the temporal order of consistency is immediately obtained from that of the underlying W-method, this property needs a more thorough analysis for the newly introduced AMF-TASE W-methods. For these latter methods it is described which are the additional order conditions to be fulfilled and it is shown that the parallel structure of the methods is crucial to retain the order of consistency of the underlying TASE W-method. Numerical experiments are presented in three spatial dimensions to assess the consistency result and to show that the proposed schemes are competitive with other well-known good performing AMF W-methods.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141532492","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A $$theta $$ -L Approach for the Simulation of Solid-State Dewetting Problems with Strongly Anisotropic Surface Energies 模拟具有强各向异性表面能量的固态润湿问题的 $$theta $$ -L 方法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-19 DOI: 10.1007/s10915-024-02589-z

Weijie Huang, Wei Jiang, Yan Wang

引用次数: 0

Second-Order Decoupled Linear Energy-Law Preserving gPAV Numerical Schemes for Two-Phase Flows in Superposed Free Flow and Porous Media 自由流与多孔介质叠加两相流体的二阶解耦线性能量守恒 gPAV 数值方案

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-19 DOI: 10.1007/s10915-024-02576-4

Yali Gao, Daozhi Han

引用次数: 0

A Compact Coupling Interface Method with Second-Order Gradient Approximation for Elliptic Interface Problems 针对椭圆界面问题的二阶梯度逼近紧凑耦合界面方法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-18 DOI: 10.1007/s10915-024-02587-1

Ray Zirui Zhang, Li-Tien Cheng

{"title":"A Compact Coupling Interface Method with Second-Order Gradient Approximation for Elliptic Interface Problems","authors":"Ray Zirui Zhang, Li-Tien Cheng","doi":"10.1007/s10915-024-02587-1","DOIUrl":"https://doi.org/10.1007/s10915-024-02587-1","url":null,"abstract":"We propose the Compact Coupling Interface Method, a finite difference method capable of obtaining second-order accurate approximations of not only solution values but their gradients, for elliptic complex interface problems with interfacial jump conditions. Such elliptic interface boundary value problems with interfacial jump conditions are a critical part of numerous applications in fields such as heat conduction, fluid flow, materials science, and protein docking, to name a few. A typical example involves the construction of biomolecular shapes, where such elliptic interface problems are in the form of linearized Poisson–Boltzmann equations, involving discontinuous dielectric constants across the interface, that govern electrostatic contributions. Additionally, when interface dynamics are involved, the normal velocity of the interface might be comprised of the normal derivatives of solution, which can be approximated to second-order by our method, resulting in accurate interface dynamics. Our method, which can be formulated in arbitrary spatial dimensions, combines elements of the highly-regarded Coupling Interface Method, for such elliptic interface problems, and Smereka’s second-order accurate discrete delta function. The result is a variation and hybrid with a more compact stencil than that found in the Coupling Interface Method, and with advantages, borne out in numerical experiments involving both geometric model problems and complex biomolecular surfaces, in more robust error profiles.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141529166","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Maximum-Norm Error Estimates of Fourth-Order Compact and ADI Compact Finite Difference Methods for Nonlinear Coupled Bacterial Systems 非线性耦合细菌系统的四阶紧凑型和 ADI 紧凑型有限差分法的最大正则误差估算

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-18 DOI: 10.1007/s10915-024-02588-0

Jie Xu, Shusen Xie, Hongfei Fu

引用次数: 0

Fourth-Order Accurate Compact Scheme for First-Order Maxwell’s Equations 一阶麦克斯韦方程的四阶精确紧凑方案

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-15 DOI: 10.1007/s10915-024-02583-5

I. Versano, E. Turkel, S. Tsynkov

引用次数: 0

Supervised Low-Rank Semi-nonnegative Matrix Factorization with Frequency Regularization for Forecasting Spatio-temporal Data 用于时空数据预测的带频率正则化的有监督低库半负矩阵因式分解法

IF 2.5 2区数学

Journal of Scientific Computing Pub Date : 2024-06-14 DOI: 10.1007/s10915-024-02565-7

Keunsu Kim, Hanbaek Lyu, Jinsu Kim, Jae-Hun Jung

{"title":"Supervised Low-Rank Semi-nonnegative Matrix Factorization with Frequency Regularization for Forecasting Spatio-temporal Data","authors":"Keunsu Kim, Hanbaek Lyu, Jinsu Kim, Jae-Hun Jung","doi":"10.1007/s10915-024-02565-7","DOIUrl":"https://doi.org/10.1007/s10915-024-02565-7","url":null,"abstract":"We propose a novel methodology for forecasting spatio-temporal data using supervised semi-nonnegative matrix factorization (SSNMF) with frequency regularization. Matrix factorization is employed to decompose spatio-temporal data into spatial and temporal components. To improve clarity in the temporal patterns, we introduce a nonnegativity constraint on the time domain along with regularization in the frequency domain. Specifically, regularization in the frequency domain involves selecting features in the frequency space, making an interpretation in the frequency domain more convenient. We propose two methods in the frequency domain: soft and hard regularizations, and provide convergence guarantees to first-order stationary points of the corresponding constrained optimization problem. While our primary motivation stems from geophysical data analysis based on GRACE (Gravity Recovery and Climate Experiment) data, our methodology has the potential for wider application. Consequently, when applying our methodology to GRACE data, we find that the results with the proposed methodology are comparable to previous research in the field of geophysical sciences but offer clearer interpretability.","PeriodicalId":50055,"journal":{"name":"Journal of Scientific Computing","volume":null,"pages":null},"PeriodicalIF":2.5,"publicationDate":"2024-06-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141512495","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0