Communications in Computational Physics
type: channel
description: Communications in Computational Physics (CiCP) is an international, peer reviewed, interdisciplinary journal edited by a group of world-wide leading scientists in diverse fields to guarantee its highest standard. It publishes original research and survey papers of high scientific value in computational modeling of physical problems. Results in multi-physics and multi-scale innovative computational methods and modeling in all physical sciences will be featured. The journal is published by Global Science Press in printed and electronic forms.
publisher: Global Science Press
language: en
publication name: Communications in Computational Physics
issn: 1815-2406
items: Structured First-Layer Initialization Pre-Training Techniques to Accelerate Training Process Based on $\varepsilon$-Rank
items: A Novel High-Resolution TENO Scheme with Adaptive Artificial Anti-Diffusion
items: The First- and Second-Order Energy Stable, Mass Conservative and Bounds-Preserving Schemes for Two-Phase Incompressible Flow with Rock Compressibility
items: Local Discontinuous Galerkin Method for Nonlinear BSPDEs of Neumann Boundary Conditions with Deep Backward Dynamic Programming Time-Marching
items: Scale-Decomposed Physics-Informed Neural Networks for Singular Perturbation Problems
items: Solving Multi-Group Neutron Diffusion Eigenvalue Problem with Decoupling Residual Loss Function
items: A Locally Mass-Conservative Enriched Petrov-Galerkin Method Without Penalty for the Darcy Flow in Porous Media
items: Maximum-Principle-Preserving and Positivity-Preserving Central WENO Schemes on Overlapping Meshes
items: Long-Time Integration of Nonlinear Wave Equations with Neural Operators
items: Steady-State Simulation of Euler Equations by the Discontinuous Galerkin Method with the Jump Filter
items: High-Order Structure-Preserving Spectral Difference Methods for Hyperbolic Conservation Laws
items: Do We Need Decay-Preserving Error Estimate for Solving Parabolic Equations with Initial Singularity?
items: Least-Squares Recovery of Dual Variables in ROMs for Saddle-Point Problems
items: Quantum Simulation of the Fokker-Planck Equation via Schrödingerization
items: Reconstruction of Dynamical Systems Without Time Label
items: Derivation of Asymptotic Limiting Equations for Multigroup Anisotropic-Scattering Transport Processes in Optically Thick Regimes
items: Towards the Efficient Calculation of Quantity of Interest from Steady Euler Equations II: A CNNs-Based Automatic Implementation
items: Neural-Network-Augmented Empirical Interpolation Method for Field Reconstruction with Noise and Vibration Tolerance
items: Thermodynamically Consistent Modeling and Energy Stable Numerical Simulation of Multicomponent Compressible Flow in Poroelastic Media
items: Enhanced Second-Order Gauss-Seidel Projection Methods for the Landau-Lifshitz Equation
items: A Highly Efficient CN-Like Scheme for the Transient Wigner Equation
items: United Filter for Jointly Estimating State and Parameters of Stochastic Dynamical Systems
items: An Improved Ghost-Point Technique for the Electrostatic Equation with Curved Boundaries and Discontinuous Permittivity in ECR Ion Source PIC Simulations
items: Weak Minimization-Sampling Network for Solving PDEs
items: Structure-Preserving SAV-RK$s$ Spectral Methods for Nonlinear Hamiltonian System
items: Deep Learning-Based Reduced-Order Methods for Fast Transient Dynamics
items: A New Approximate Method to Mean Field Stochastic Differential Equation with One-Sided Lipschitz Drift Coefficient
items: Uncertainty Quantification of Phase Transition Problems with an Injection Boundary
items: Adaptive Neural Network Basis Methods for Partial Differential Equations with Low-Regular Solutions
items: A Unifying Moving Mesh Method for Curves, Surfaces, and Domains Based on Mesh Equidistribution and Alignment

Structured First-Layer Initialization Pre-Training Techniques to Accelerate Training Process Based on $\varepsilon$-Rank
type: item
description: Training deep neural networks for scientific computing remains computationally expensive due to the slow formation of diverse feature representations in early training stages. Recent studies [37] identify a staircase phenomenon in training dynamics, where loss decreases are closely correlated with increases in $\varepsilon$-rank, reflecting the effective number of linearly independent neuron functions. Motivated by this observation, this work proposes a structured first-layer initialization (SFLI) pre-training technique to enhance the diversity of neural features at initialization by constructing $\varepsilon$-linearly independent neurons in the input layer. We present systematic initialization schemes compatible with various activation functions and integrate the strategy into multiple neural architectures, including modified multi-layer perceptrons and physics-informed residual adaptive networks. Only needing to add one line of code to conventional stochastic gradient descent algorithms, extensive numerical experiments on function approximation and PDE benchmarks, demonstrate that SFLI significantly improves the initial $\varepsilon$-rank, accelerates convergence, mitigates spectral bias, and enhances prediction accuracy.
subject: Articles
  topic: section
subject: Staircase phenomenon
  topic: keywords
subject: ε-rank
  topic: keywords
subject: structured first-layer initialization
  topic: keywords
subject: deep neural network
  topic: keywords
subject: training acceleration
  topic: keywords
creator: Tao Tang
creator: Jiang Yang
creator: Yuxiang Zhao
creator: Quanhui Zhu
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-03-16
publication date: 2026-03-16
doi: 10.4208/cicp.OA-2025-0185

A Novel High-Resolution TENO Scheme with Adaptive Artificial Anti-Diffusion
type: item
description: Low-dissipation higher-order approaches are crucial for flows characterized by shocks and turbulence. The broad applicability of conventional WENO/TENO systems for accurately capturing shocks makes them a preferred option for simulating complex flow phenomena. Nonetheless, the inherent numerical diffusion in these techniques poses a significant challenge in simulating multiscale turbulent flow. This study mitigates the inherent numerical damping by using locally adaptive higher-order artificial anti-diffusion. A switching function quantifies the extent of anti-dissipation in the smooth zone and assesses the local variation of flow field variables. To assure stability, the anti-diffusion coefficient is meticulously assessed using two fixed free parameters. Next strategy utilizes a fourth-order central difference scheme in the highly smooth region, incorporating an additional coefficient modified in a binary fashion. The adaptive regulating parameter $C_T$ in the fifth-order TENO scheme is integrated with the anti-diffusive function. It is recognized that the new scheme exhibits less numerical diffusion than the original TENO scheme. This is also reflected in the numerical experiments, encompassing a series of benchmark test scenarios featuring shock and fine-scale structures. It is shown that the proposed scheme could resolve the turbulent scales with much higher resolution and capture the shock sharply. It is worth noting that these strategies can easily be implemented in the other variations of the TENO scheme.
subject: Articles
  topic: section
subject: ENO-like stencil selection
  topic: keywords
subject: adaptive artificial anti-diffusion
  topic: keywords
subject: WENO
  topic: keywords
subject: TENO
  topic: keywords
subject: low-dissipation schemes
  topic: keywords
subject: fourth-order central difference
  topic: keywords
creator: Velmani Murugesan
creator: Manoj T. Nair
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-03-16
publication date: 2026-03-16
doi: 10.4208/cicp.OA-2025-0113

The First- and Second-Order Energy Stable, Mass Conservative and Bounds-Preserving Schemes for Two-Phase Incompressible Flow with Rock Compressibility
type: item
description: In this paper, we construct a unified framework of first- and second-order schemes in time for thermodynamically consistent modeling of two-phase incompressible and immiscible flow in porous media with rock compressibility. We rigorously prove that the proposed schemes can preserve the energy dissipation law, conserve the mass of each phase as well as pore volumes and are bounds-preserving for both phases without any restrictions on the time step size. Moreover, the proposed schemes can achieve local mass conservation for both phases and preserve the pore volumes when the saturation of the whole region is within the bounds, and we only need to solve one linear system and several linear algebraic equations, whereas for points where saturation is outside the bounds, we just need to solve an additional nonlinear algebraic equation and hence the proposed schemes can achieve global mass conservation. The developed schemes exhibit high efficiency owing to a substantial reduction in the scale of nonlinear computations and their straightforward implementation. The key point is that we propose a new Lagrange multiplier method that is based on the judgment of saturation bound. This addresses the limitation of the classical Lagrange multiplier method, which is only capable of ensuring global mass conservation. Finally, a variety of illustrative numerical examples including several benchmark problems are provided to validate the accuracy and efficiency of the proposed schemes.
subject: Articles
  topic: section
subject: Two-phase flow
  topic: keywords
subject: energy evolution equation
  topic: keywords
subject: Lagrange multiplier
  topic: keywords
subject: rock compressibility
  topic: keywords
creator: Yujing Yan
creator: Huangxin Chen
creator: Xiaoli Li
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1332
ending page: 1362
doi: 10.4208/cicp.OA-2025-0089

Local Discontinuous Galerkin Method for Nonlinear BSPDEs of Neumann Boundary Conditions with Deep Backward Dynamic Programming Time-Marching
type: item
description: This paper aims to present a local discontinuous Galerkin (LDG) method for solving nonlinear backward stochastic partial differential equations (BSPDEs) with Neumann boundary conditions. We establish the $L^2$-stability and optimal error estimates of the proposed numerical scheme. Two numerical examples are provided to demonstrate the performance of the LDG method, where we incorporate a deep learning algorithm to address the challenge of the curse of dimensionality in backward stochastic differential equations (BSDEs). The results show the effectiveness and accuracy of the LDG method in tackling BSPDEs with Neumann boundary conditions.
subject: Articles
  topic: section
subject: Local discontinuous Galerkin method
  topic: keywords
subject: backward stochastic partial differential equations
  topic: keywords
subject: Neumann boundary problems
  topic: keywords
subject: stability analysis
  topic: keywords
subject: error estimates
  topic: keywords
subject: deep learning algorithm
  topic: keywords
creator: Yixiang Dai
creator: Yunzhang Li
creator: Jing Zhang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1299
ending page: 1331
doi: 10.4208/cicp.OA-2024-0296

Scale-Decomposed Physics-Informed Neural Networks for Singular Perturbation Problems
type: item
description: A novel physics-informed neural network (PINN), scale-decomposed PINN (SD-PINN), is proposed for singular perturbation problems, such as boundary-layer flows. Singular perturbation problems exhibit sharp changes in solutions at different scales, owing to large gradients in equations, which presents a challenge to a conventional PINN. In SD-PINN, the solutions at different scales are represented by different neural networks, and their matching is achieved by imposed matching conditions. In comparison with the conventional procedure for singular perturbation problems, no domain decomposition is used. A nonlinear stretching transformation is introduced to prevent the occurrence of semi-infinite regions in the neural networks. Six benchmark singular perturbation problems are used to evaluate SD-PINN, including linear and nonlinear second-order ordinary differential equations, two-dimensional hyperbolic problems, and a flat-plate boundary-layer flow. It is demonstrated that SD-PINN can reproduce asymptotic and non-asymptotic solutions for small and finite perturbation parameters, respectively. Application of SD-PINN to the inference of perturbation parameters from available data is also discussed. Codes are available at https://github.com/LeiZhang-code/SD-PINN.
subject: Articles
  topic: section
subject: Singular perturbation
  topic: keywords
subject: physics-informed neural network
  topic: keywords
subject: boundary layer
  topic: keywords
subject: machine learning
  topic: keywords
creator: Lei Zhang
creator: Guowei He
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1267
ending page: 1298
doi: 10.4208/cicp.OA-2024-0149

Solving Multi-Group Neutron Diffusion Eigenvalue Problem with Decoupling Residual Loss Function
type: item
description: In the midst of the neural network’s success in solving partial differential equations, tackling eigenvalue problems using neural networks remains a challenging task. However, the Physics Constrained-General Inverse Power Method Neural Network (PC-GIPMNN) approach was proposed and successfully applied to solve the single-group critical problems in reactor physics. This paper aims to solve critical problems in multi-group scenarios and in more complex geometries. Hence, inspired by the merits of traditional source iterative method, which can overcome the ill-condition of the right side of the equations effectively and solve the multi-group problem effectively, we propose two residual loss function called Decoupling Residual loss function and Direct Iterative loss function. Our loss function can deal with multi-group eigenvalue problem, and also single-group eigenvalue problem. Using the new residual loss functions, our study solves one-dimensional, two-dimensional, and three-dimensional multi-group problems in nuclear reactor physics without prior data. In numerical experiments, our approach demonstrates superior generalization capabilities compared to previous work.
subject: Articles
  topic: section
subject: Deep learning
  topic: keywords
subject: eigenvalue problem
  topic: keywords
subject: nuclear reactors
  topic: keywords
subject: multi-group problem
  topic: keywords
creator: Shupei Yu
creator: Qiaolin He
creator: Shiquan Zhang
creator: Qihong Yang
creator: Yu Yang
creator: Helin Gong
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1231
ending page: 1266
doi: 10.4208/cicp.OA-2024-0176

A Locally Mass-Conservative Enriched Petrov-Galerkin Method Without Penalty for the Darcy Flow in Porous Media
type: item
description: In this work we present an enriched Petrov-Galerkin (EPG) method for the simulation of the Darcy flow in porous media. The new method enriches the approximation trial space of the conforming continuous Galerkin (CG) method with bubble functions and enriches the approximation test space of the CG method with piecewise constant functions, and it does not require any penalty term in the weak formulation. Moreover, we propose a framework for constructing the bubble functions and consider a decoupled algorithm for the EPG method based on this framework, which enables the process of solving pressure to be decoupled into two steps. The first step is to solve the pressure by the standard CG method, and the second step is a post-processing correction of the first step. Compared with the CG method, the proposed EPG method is locally mass-conservative, while keeping fewer degrees of freedom than the discontinuous Galerkin (DG) method. In addition, this method is more concise in the error analysis than the enriched Galerkin (EG) method. The coupled flow and transport in porous media is considered to illustrate the advantages of locally mass-conservative properties of the EPG method. We establish the optimal convergence of numerical solutions and present several numerical examples to illustrate the performance of the proposed method.
subject: Articles
  topic: section
subject: Enriched Petrov-Galerkin method
  topic: keywords
subject: local mass conservation
  topic: keywords
subject: post-processing
  topic: keywords
subject: error analysis
  topic: keywords
creator: Huangxin Chen
creator: Piaopiao Dong
creator: Shuyu Sun
creator: Zixuan Wang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1199
ending page: 1230
doi: 10.4208/cicp.OA-2023-0321

Maximum-Principle-Preserving and Positivity-Preserving Central WENO Schemes on Overlapping Meshes
type: item
description: In this paper, we propose a class of maximum-principle-preserving central WENO schemes for scalar conservation laws, and positivity-preserving central WENO schemes for compressible Euler equations. Formulated in a finite volume framework on overlapping meshes, the central schemes require neither flux splitting nor numerical fluxes that are often exact or approximate Riemann solvers. A new fifth-order WENO reconstruction is applied for the spatial discretization, and the linear weights of such reconstruction can be any positive number as long as their sum equals one, which leads to much simpler implementation. The sufficient conditions are provided for the cell average values to preserve maximum principle or positivity property with Euler forward time discretization. The method can be generalized to high-order strong stability preserving Runge-Kutta method without technical difficulties. Extensive numerical examples are presented to illustrate the accuracy and performance of the proposed methods.
subject: Articles
  topic: section
subject: Maximum-principle-preserving
  topic: keywords
subject: positivity-preserving
  topic: keywords
subject: central scheme
  topic: keywords
subject: WENO reconstruction
  topic: keywords
subject: overlapping mesh
  topic: keywords
creator: Zhanjing Tao
creator: Zhengfu Xu
creator: Jianxian Qiu
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1164
ending page: 1198
doi: 10.4208/cicp.OA-2024-0165

Long-Time Integration of Nonlinear Wave Equations with Neural Operators
type: item
description: Neural operators have shown promise in solving many types of Partial Differential Equations (PDEs). They are significantly faster compared to traditional numerical solvers once they have been trained with a certain amount of observed data. However, their numerical performance in solving time-dependent PDEs, particularly in long-time prediction of dynamic systems, still needs improvement. In this paper, we focus on solving the long-time integration of nonlinear wave equations via neural operators by replacing the initial condition with the prediction in a recurrent manner. Given limited observed temporal trajectory data, we utilize some intrinsic features of these nonlinear wave equations, such as conservation laws and well-posedness, to improve the algorithm design and reduce accumulated error. Our numerical experiments examine these improvements in the Korteweg-de Vries (KdV) equation, the sine-Gordon equation, and the Klein-Gordon wave equation on the irregular domain.
subject: Articles
  topic: section
subject: Neural operator
  topic: keywords
subject: nonlinear wave equations
  topic: keywords
subject: long-time integration
  topic: keywords
creator: Guanhang Lei
creator: Zhen Lei
creator: Lei Shi
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1137
ending page: 1163
doi: 10.4208/cicp.OA-2024-0274

Steady-State Simulation of Euler Equations by the Discontinuous Galerkin Method with the Jump Filter
type: item
description: In this study, we introduce a jump filter for discontinuous Galerkin (DG) methods on both rectangular and triangular meshes, specifically designed for steady-state solutions of the Euler equations. Traditional limiters have occasionally been observed to produce slight post-shock oscillations, but the integration of the jump filter effectively reduces these fluctuations. This improvement allows the residual to stabilize at a level indistinguishable from machine zero in steady-state simulations. The jump filter operates by applying distinct damping factors to the different moments of the solution, based on the jump information at cell interfaces. This strategy creates transition zones from smooth regions to shock waves, which is critical for achieving steady-state solutions with residuals approaching the threshold of machine zero. The jump filter on unstructured meshes using Dubiner polynomials associated with the Sturm-Liouville operator is also investigated. Furthermore, the DG method with the jump filter retains many of the key advantages of the classical DG approach, including compactness, conservation, stability, and high-order accuracy. Extensive numerical experiments, conducted on both rectangular and triangular meshes, confirm the effectiveness of the jump filter in not only suppressing numerical oscillations but also significantly driving the residual towards machine zero.
subject: Articles
  topic: section
subject: Euler equations
  topic: keywords
subject: discontinuous Galerkin method
  topic: keywords
subject: local viscosity
  topic: keywords
subject: jump filter
  topic: keywords
subject: steady-state convergence
  topic: keywords
subject: rectangular and triangular meshes
  topic: keywords
creator: Lei Wei
creator: Yinhua Xia
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1103
ending page: 1136
doi: 10.4208/cicp.OA-2025-0091

High-Order Structure-Preserving Spectral Difference Methods for Hyperbolic Conservation Laws
type: item
description: This paper presents a study focused on solving hyperbolic conservation laws using arbitrary-order spectral difference (SD) methods. The study is structured around several crucial aspects. Firstly, in order to ensure the maximum principle for scalar equations and positivity preservation for Euler systems, we adopt the concept of flux limiters. This adoption leads to the development of the structure-preserving SD scheme with a flux limiter (SDFL), which is proven to preserve the original high-order accuracy. However, the SDFL scheme with lower order might lack conservational properties, despite its strong performance in short-term simulations. Consequently, we have developed a specific variant of the SDFL scheme with conservational properties, referred to as the CSDFL scheme. Secondly, we introduce a modified WENO-ZQ (MWENO-ZQ) reconstruction to suppress spurious oscillations when simulating problems with strong discontinuities. Finally, we conduct extensive numerical experiments to validate the effectiveness of the proposed high-order SD (SDFL, CSDFL) methods with MWENO-ZQ reconstruction. The results demonstrate the robustness and efficiency of these techniques in solving problems involving strong discontinuities, low pressure, and low density.
subject: Articles
  topic: section
subject: Spectral difference method
  topic: keywords
subject: conservation law
  topic: keywords
subject: maximum principle
  topic: keywords
subject: positivity-preserving
  topic: keywords
creator: Yupeng Ren
creator: Zhimin Zhang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1059
ending page: 1102
doi: 10.4208/cicp.OA-2024-0058

Do We Need Decay-Preserving Error Estimate for Solving Parabolic Equations with Initial Singularity?
type: item
description: Solutions exhibiting weak initial singularities arise in various equations, including diffusion and subdiffusion equations. When employing the well-known L1 scheme to solve subdiffusion equations with weak singularities, numerical simulations reveal that this scheme exhibits varying convergence rates for different choices of model parameters (i.e., domain size, final time $T$, and reaction coefficient $\kappa$). This elusive phenomenon is not unique to the L1 scheme but is also observed in other numerical methods for reaction-diffusion equations such as the backward Euler (IE) scheme, Crank-Nicolson (C-N) scheme, and two-step backward differentiation formula (BDF2) scheme. The existing literature lacks an explanation for the existence of two different convergence regimes, which has puzzled us for a long while and motivated us to study this inconsistency between the standard convergence theory and numerical experiences. In this paper, we provide a general methodology to systematically obtain error estimates that incorporate the exponential decaying feature of the solution. We term this novel error estimate the ‘decay-preserving error estimate’ and apply it to the aforementioned IE, C-N, and BDF2 schemes. Our decay-preserving error estimate consists of a low-order term with an exponential coefficient and a high-order term with an algebraic coefficient, both of which depend on the model parameters. Our estimates reveal that the varying convergence rates are caused by a trade-off between these two components in different model parameter regimes. By considering the model parameters, we capture different states of the convergence rate that traditional error estimates fail to explain. This approach retains more properties of the continuous solution. We validate our analysis with numerical results.
subject: Articles
  topic: section
subject: Diffusion equation
  topic: keywords
subject: subdiffusion equation
  topic: keywords
subject: weak initial singularity
  topic: keywords
subject: decay-preserving error estimate
  topic: keywords
subject: convergence rate
  topic: keywords
creator: Jiwei Zhang
creator: Zhimin Zhang
creator: Chengchao Zhao
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1027
ending page: 1058
doi: 10.4208/cicp.OA-2025-0130

Least-Squares Recovery of Dual Variables in ROMs for Saddle-Point Problems
type: item
description: In this work, we introduce a methodology for recovering the reduced Lagrange multiplier in Reduced Order Models (ROMs) of general saddle-point problems. Specifically, the multiplier is determined as the least-squares minimizer of the residual of the reduced primal variable in a dual norm. We prove that this procedure yields a unique solution under the condition that the full-order primal and dual spaces satisfy the inf-sup stability condition. This is an extension to general saddle-point problems of the method recently introduced in [6] to recover the reduced pressure in ROMs of incompressible flows. We further show that the proposed approach is equivalent to solving the reduced mixed problem with an enriched reduced primal basis, augmented by the supremizers of the reduced multiplier. We establish optimal error estimates for the reduced multiplier applicable to general saddle-point problems. Numerical experiments for heat diffusion in composite media validate our theoretical findings and demonstrate the suitability of the method to solve engineering problems.
subject: Articles
  topic: section
subject: Reduced order models (ROM)
  topic: keywords
subject: saddle point problems
  topic: keywords
subject: dual variable recovery
  topic: keywords
creator: M. Azaïez
creator: T. Chacόn Rebollo
creator: M. Oulghelou
creator: I. Sánchez Muñoz
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 1002
ending page: 1026
doi: 10.4208/cicp.OA-2025-0040

Quantum Simulation of the Fokker-Planck Equation via Schrödingerization
type: item
description: This paper studies a quantum simulation technique for solving the Fokker-Planck equation. Traditional semi-discretization methods often fail to preserve the underlying Hamiltonian dynamics and may even modify the Hamiltonian structure, particularly when incorporating boundary conditions. We address this challenge by employing the Schrödingerization method – it converts any linear partial and ordinary differential equation with non-Hermitian dynamics into systems of Schrödinger-type equations. It does so via the so-called warped phase transformation that maps the equation into one higher dimension. We explore the application in two distinct forms of the Fokker-Planck equation. For the conservation form, we show that the semi-discretization-based Schrödingerization is preferable, especially when dealing with non-periodic boundary conditions. Additionally, we analyze the Schrödingerization approach for unstable systems that possess positive eigenvalues in the real part of the coefficient matrix or differential operator. Our analysis reveals that the direct use of Schrödingerization has the same effect as a stabilization procedure. For the heat equation form, we propose a quantum simulation procedure based on the time-splitting technique, and give explicitly its corresponding quantum circuit. We discuss the relationship between operator splitting in the Schrödingerization method and its application directly to the original problem, illustrating how the Schrödingerization method accurately reproduces the time-splitting solutions at each step. Furthermore, we explore finite difference discretizations of the heat equation form using shift operators. Utilizing Fourier bases, we diagonalize the shift operators, enabling efficient simulation in the frequency space. Providing additional guidance on implementing the diagonal unitary operators, we conduct a comparative analysis between diagonalizations in the Bell and the Fourier bases, and show that the former generally exhibit…
subject: Articles
  topic: section
subject: Fokker-Planck equation
  topic: keywords
subject: Schrödingerization
  topic: keywords
subject: quantum simulation
  topic: keywords
subject: shift operators
  topic: keywords
subject: diagonal unitary operators
  topic: keywords
creator: Shi Jin
creator: Nana Liu
creator: Yue Yu
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-18
publication date: 2026-01-18
starting page: 969
ending page: 1001
doi: 10.4208/cicp.OA-2024-0097

Reconstruction of Dynamical Systems Without Time Label
type: item
description: In this paper, we study the method to reconstruct dynamical systems from data without time labels. Data without time labels appear in many applications, such as molecular dynamics, single-cell RNA sequencing, etc. Reconstruction of dynamical system from time sequence data has been studied extensively. However, these methods do not apply if time labels are unknown. Without time labels, sequence data become distribution data. Based on this observation, we propose to treat the data as samples from a probability distribution and try to reconstruct the underlying dynamical system by minimizing the distribution loss, sliced Wasserstein distance more specifically. Extensive experiment results demonstrate the effectiveness of the proposed method.
subject: Articles
  topic: section
subject: Dynamical system recovery
  topic: keywords
subject: data without time label
  topic: keywords
subject: Wasserstein distance
  topic: keywords
creator: Zhijun Zeng
creator: Chenglong Bao
creator: Pipi Hu
creator: Yi Zhu
creator: Zuoqiang Shi
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 941
ending page: 968
doi: 10.4208/cicp.OA-2025-0042

Derivation of Asymptotic Limiting Equations for Multigroup Anisotropic-Scattering Transport Processes in Optically Thick Regimes
type: item
description: Asymptotic limit methods play a crucial role in deriving macroscopic fluid equations from equations governing mesoscopic particle transport mechanics, constituting a significant research area within kinetic theory. In this paper, we present asymptotic expansions applied to the multigroup linear transport equation with a general anisotropic scattering kernel. We derive the diffusion limit for such transport equations under three different scalings. For the conventional scaling, the standard multigroup diffusion equations are derived as an asymptotic approximation to the multigroup transport equations. For the other two scalings where the scattering effect between groups is strong, new group-collapsed (monoenergetic) diffusion approximations with different coefficients depending on the form of the scaling are derived. Numerical experiments validate the prediction of the theoretical models. This paper provides a deeper understanding of the diffusion physics for the multigroup anisotropic particle transport in the optically thick regime.
subject: Articles
  topic: section
subject: Transport equation
  topic: keywords
subject: multigroup approximation
  topic: keywords
subject: asymptotic analysis
  topic: keywords
subject: diffusion limit
  topic: keywords
creator: Weiming Li
creator: Chang Liu
creator: Peng Song
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 919
ending page: 940
doi: 10.4208/cicp.OA-2024-0181

Towards the Efficient Calculation of Quantity of Interest from Steady Euler Equations II: A CNNs-Based Automatic Implementation
type: item
description: In [J. Wang, G. Hu, arxiv: 2302.14262], a dual-consistent dual-weighted residual-based $h$-adaptive method has been proposed based on a Newton-GMG framework, toward the accurate calculation of a given quantity of interest from Euler equations. The performance of such a numerical method is satisfactory, i.e., the stable convergence of the quantity of interest can be observed. In this paper, we will focus on the efficiency issue to further develop this method. Three approaches are studied for addressing the efficiency issue, i.e., i). using convolutional neural networks as a solver for dual equations, ii). designing an automatic adjustment strategy for the tolerance in the $h$-adaptive process to conduct the local refinement and/or coarsening of mesh grids, and iii). introducing OpenMP, a shared memory parallelization technique, to accelerate the module such as the solution reconstruction in the method. The feasibility of each approach and numerical issues are discussed in depth, and significant acceleration from those approaches in simulations can be observed clearly from a number of numerical experiments. In convolutional neural networks, it is worth mentioning that the dual consistency plays an important role in guaranteeing the efficiency of the whole method and that unstructured meshes are employed in all simulations.
subject: Articles
  topic: section
subject: DWR-based $h$-adaptivity
  topic: keywords
subject: dual consistency
  topic: keywords
subject: convolutional neural networks
  topic: keywords
subject: automatic adjustment of tolerance
  topic: keywords
subject: Newton-GMG method for Euler equations
  topic: keywords
creator: Jingfeng Wang
creator: Guanghui Hu
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 884
ending page: 918
doi: 10.4208/cicp.OA-2023-0236

Neural-Network-Augmented Empirical Interpolation Method for Field Reconstruction with Noise and Vibration Tolerance
type: item
description: Field reconstruction methods based on machine learning often face challenges when handling unstructured data and designing complex model structures. On the other hand, Reduced Order Models (ROM) struggle with robustness in the presence of observation noise and sensor vibrations. Striking a balance between lightweight design, the ability to handle unstructured data, and robustness presents a significant challenge. In this paper, we introduce the EIM-NN algorithm, which leverages neural networks to determine the coefficients of the subspace found by the EIM algorithm. Furthermore, we present the EIM-TNN algorithm, which enhances robustness by designing a loss function incorporating Tikhonov regularization. Our neural network consists of only two fully connected layers, allowing it to handle unstructured data while maintaining a lightweight profile. Experimental results demonstrate the algorithm’s ability to significantly enhance robustness against noise and sensor vibrations without compromising its accuracy in fitting the original data. Additionally, the algorithm’s lightweight nature ensures that the added training time and memory requirements over EIM remain acceptable, making it adaptable to a range of industrial applications.
subject: Articles
  topic: section
subject: Field reconstruction
  topic: keywords
subject: EIM-TNN
  topic: keywords
subject: lightweight
  topic: keywords
subject: observation noise
  topic: keywords
subject: sensor vibration
  topic: keywords
creator: Han Li
creator: Helin Gong
creator: Chuanju Xu
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 855
ending page: 883
doi: 10.4208/cicp.OA-2024-0099

Thermodynamically Consistent Modeling and Energy Stable Numerical Simulation of Multicomponent Compressible Flow in Poroelastic Media
type: item
description: Modeling and numerical simulation of coupled poromechanical problems with multicomponent compressible flow are of particular importance in many fields including shale and natural gas engineering, carbon dioxide sequestration and geotechnical engineering. In this paper, using the second law of thermodynamics, we rigorously derive a thermodynamically consistent model for multicomponent compressible flow in deformable porous media coupled with poroelasticity. The model herein takes molar densities as the primary unknowns rather than pressure and molar fractions as well as introduces fluid and solid free energies, so that it naturally follows an energy dissipation law. Additionally, the Maxwell-Stefan model of multicomponent diffusion is generalized as a multicomponent fluid-solid coupling model accounting for the solid deformation, which not only satisfies Onsager’s reciprocal principle but also yields a thermodynamically consistent poro-visco-elastic equation. For numerical simulation, we propose a novel energy stable and mass conservative numerical scheme for the model. We first design the semi-discrete time scheme using the stabilized energy factorization approach to deal with the multicomponent Helmholtz free energy as well as subtle semi-implicit treatments for the coupling between multicomponent fluids and solids. A nontrivial treatment is the use of the discrete Gibbs-Duhem equation substituting for the pressure gradient contributed by the multicomponent fluids in the solid mechanical balance equation, which establishes the thermodynamically consistent relation between fluids and solids at the discrete level. Based on the cell-centered finite volume method on staggered grids, the fully discrete scheme is constructed using the upwind strategy for both molar densities and porosity. The scheme is proved to preserve the discrete energy dissipation law and Onsager’s reciprocal principle as well as to conserve the mass of fluid components and solids. Numerical experim…
subject: Articles
  topic: section
subject: Multicomponent flow
  topic: keywords
subject: poroelasticity
  topic: keywords
subject: Maxwell-Stefan model
  topic: keywords
subject: porous media
  topic: keywords
subject: thermodynamic consistency
  topic: keywords
subject: energy stability
  topic: keywords
creator: Jisheng Kou
creator: Amgad Salama
creator: Huangxin Chen
creator: Shuyu Sun
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 815
ending page: 854
doi: 10.4208/cicp.OA-2023-0254

Enhanced Second-Order Gauss-Seidel Projection Methods for the Landau-Lifshitz Equation
type: item
description: The dynamics of magnetization in ferromagnetic materials are modeled by the Landau-Lifshitz equation, which presents significant challenges due to its inherent nonlinearity and non-convex constraint. These complexities necessitate efficient numerical methods for micromagnetics simulations. The Gauss-Seidel Projection Method (GSPM), first introduced in 2001, is among the most efficient techniques currently available. However, existing GSPMs are limited to first-order accuracy. This paper introduces two novel second-order accurate GSPMs based on a combination of the biharmonic equation and the second-order backward differentiation formula, achieving computational complexity comparable to that of solving the scalar biharmonic equation implicitly. The first proposed method achieves unconditional stability through Gauss-Seidel updates, while the second method exhibits conditional stability with a Courant-Friedrichs-Lewy constant of 0.25. Through consistency analysis and numerical experiments, we demonstrate the efficacy and reliability of these methods. Notably, the first method displays unconditional stability in micromagnetics simulations, even when the stray field is updated only once per time step.
subject: Articles
  topic: section
subject: Gauss-Seidel projection methods
  topic: keywords
subject: Landau-Lifshitz equation
  topic: keywords
subject: backward differentiation formula
  topic: keywords
subject: second-order accuracy
  topic: keywords
creator: Panchi Li
creator: Xiao-Ping Wang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 799
ending page: 814
doi: 10.4208/cicp.OA-2025-0018

A Highly Efficient CN-Like Scheme for the Transient Wigner Equation
type: item
description: The transient Wigner equation (TWE) plays an important role in modeling quantum effects of nano-scale semiconductor devices. The pseudo-differential term contained within the TWE provides quantum descriptions for the model, however, it is costly in numerical simulations. We develop a CN-like scheme for time integration of the TWE, based on the two-step method. Additionally, the spatial discretization is the hybrid finite-difference/sinc-Galerkin scheme [H. Jiang, T. Lu, W. Yao and W. Zhang, SIAM J. Sci. Comput. 45 (2023)]. Rigorous proofs are provided to show that the CN-like scheme is unconditionally $L^2$-stable and has second-order accuracy in time. More importantly, the computational efficiency of the new CN-like scheme could be considered to be higher than that of any explicit multi-stage one-step time integration scheme. Numerical experiments are also carried out to verify the accuracy, stability and efficiency of the new CN-like scheme. In addition to the verification experiments, resonant tunneling diodes (RTDs) with various sets of parameters are simulated. I-V characteristics and associated curves are obtained to demonstrate the resonant tunneling effect and the variations in I-V characteristic curves in relation to changes in structure parameters.
subject: Articles
  topic: section
subject: Transient Wigner equation
  topic: keywords
subject: resonant tunneling diode
  topic: keywords
subject: two-step method
  topic: keywords
subject: unconditionally $L^2$-stable
  topic: keywords
subject: high efficiency
  topic: keywords
creator: Yidan Wang
creator: Haiyan Jiang
creator: Tiao Lu
creator: Wenqi Yao
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 775
ending page: 798
doi: 10.4208/cicp.OA-2024-0105

United Filter for Jointly Estimating State and Parameters of Stochastic Dynamical Systems
type: item
description: This paper tackles the intricate task of jointly estimating state and parameters in data assimilation for stochastic dynamical systems that are affected by noise and observed only partially. While the concept of “optimal filtering” serves as the customary approach to estimate the state of the target dynamical system, traditional methods such as Kalman filters and particle filters encounter significant challenges when dealing with high-dimensional and nonlinear problems. When we also consider the scenario where the model parameters are unknown, the problem transforms into a joint state-parameter estimation problem. Presently, the leading-edge technique known as the Augmented Ensemble Kalman Filter (AugEnKF) addresses this issue by treating unknown parameters as additional state variables and employing the Ensemble Kalman Filter to estimate the augmented state-parameter vector. Despite its considerable progress, AugEnKF does exhibit certain limitations in terms of accuracy and stability. To address these challenges, we introduce an innovative approach, referred to as the United Filter. This method combines a remarkably stable and efficient ensemble score filter (EnSF) for state estimation with a precise direct filter dedicated to online parameter estimation. Utilizing the EnSF’s generative capabilities grounded in diffusion models, the United Filter iteratively fine-tunes both state and parameter estimates within a single temporal data assimilation step. Thanks to the robustness of the EnSF, the proposed United Filter method offers a promising solution for enhancing our understanding and modeling of dynamical systems, as demonstrated by results from numerical experiments.
subject: Articles
  topic: section
subject: Optimal filtering
  topic: keywords
subject: state estimation
  topic: keywords
subject: parameter estimation
  topic: keywords
subject: diffusion model
  topic: keywords
subject: ensemble score filter
  topic: keywords
subject: Bayesian inference
  topic: keywords
creator: Feng Bao
creator: Zezhong Zhang
creator: Guannan Zhang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 747
ending page: 774
doi: 10.4208/cicp.OA-2024-0009

An Improved Ghost-Point Technique for the Electrostatic Equation with Curved Boundaries and Discontinuous Permittivity in ECR Ion Source PIC Simulations
type: item
description: In this paper, we propose an efficient Poisson solver for a two-dimensional axisymmetric domain to simulate plasma in a 2.45 GHz Electron Cyclotron Resonance (ECR) ion source. The solver employs a uniform Cartesian mesh for a curved domain, utilizing an unfitted boundary method. The presence of two Boron Nitride disks inside the plasma chamber creates a discontinuity in permittivity, leading to a Poisson problem with discontinuous coefficients. We discretize the Laplace operator by standard central difference formulas, and apply a ghost-point technique at nodes outside the domain and near the interfaces between the disks and the chamber, enforcing appropriate boundary and interface conditions for the ghost values. This method significantly simplifies and accelerates computations compared to boundary-conforming mesh approaches. The primary objective of this work is to upgrade our Poisson solver (employed in the Particle-In-Cell (PIC) code) from a fitted boundary method, such as COMSOL, to a custom implementation based on an unfitted boundary approach. The performance improvements are substantial, reducing computation time from a minimum of 5 seconds per solution with COMSOL to just 0.062 seconds with our custom solver.
subject: Articles
  topic: section
subject: Ghost-point method
  topic: keywords
subject: elliptic interface problems
  topic: keywords
subject: numerical treatment of material discontinuities
  topic: keywords
subject: finite-difference discretization
  topic: keywords
subject: embedded boundary techniques
  topic: keywords
subject: charge-particle simulation
  topic: keywords
subject: plasma source modeling
  topic: keywords
subject: computational electromagnetics
  topic: keywords
creator: Armando Coco
creator: Lorenzo Neri
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 713
ending page: 746
doi: 10.4208/cicp.OA-2025-0100

Weak Minimization-Sampling Network for Solving PDEs
type: item
description: In many engineering applications, people concentrate on deriving the weak solution to a given partial differential equation (PDE) model. Traditionally, both the trial space and the test space are constructed as the span of a set of basis functions with compact supports and the approximate solution is achieved by solving a (linear) system of equations. A recent advancement in this field is the Weak Adversarial Network (WAN), a machine learning-based method that employs deep neural networks (DNNs) to approximate trial and test spaces. WAN addresses PDEs through a min-max optimization framework, showcasing its effectiveness in solving high dimensional PDEs and surpassing the limitations of traditional numerical methods. In this paper, motivated by WAN, we propose a Weak Minimization-Sampling Network (WMSN) for solving PDEs which also leverages DNNs to approximate the trial and test spaces. Instead of the min-max loss function, our key idea is to reformulate the weak formulation into a loss function defined as an expectation over all possible DNNs in the test space with only the network parameters in the trial space to be trained. This leads to a minimization problem in the trial space and a sampling problem over the test space, which can be solved by stochastic gradient descent type methods with neat implementation. Through a series of numerical tests, we find that WMSN exhibits comparable or slightly better results with smaller computational cost when compared to WAN.
subject: Articles
  topic: section
subject: Partial differential equations
  topic: keywords
subject: weak formulation
  topic: keywords
subject: deep neural network
  topic: keywords
subject: random sampling
  topic: keywords
creator: Zhenqian E
creator: Weiguo Gao
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 696
ending page: 712
doi: 10.4208/cicp.OA-2024-0153

Structure-Preserving SAV-RK$s$ Spectral Methods for Nonlinear Hamiltonian System
type: item
description: It is well-established that high-order symplectic Runge–Kutta methods cannot preserve polynomial energy of Hamiltonian systems beyond degree two. However, the scalar auxiliary variable (SAV) approach, which was originally proposed for gradient flow systems, offers an effective strategy to overcome this limitation, albeit with a modified energy. In this paper, we study high-order Runge–Kutta methods for solving general Hamiltonian systems without imposing a lower bound on the high-order part of the energy functional. By integrating a symplectic Runge–Kutta method with a new SAV strategy, we develop a class of high-order nonlinear structure-preserving $s$-stage Runge–Kutta SAV (SAV-RK$s$) spectral methods, and prove that the nonlinear semi-discrete scheme can preserve the Hamiltonian energy and other conservative quantities such as mass/momentum (if they exist). In addition, we design a tailored fast Newton iteration algorithm to efficiently solve the resulting nonlinear algebraic system. Finally, we carry out extensive numerical simulations on several benchmark problems, including cases where the energy functional is either bounded or unbounded, to validate the accuracy and robustness of the proposed algorithms.
subject: Articles
  topic: section
subject: Hamiltonian system
  topic: keywords
subject: symplectic Runge–Kutta method
  topic: keywords
subject: SAV approach
  topic: keywords
subject: spectral method
  topic: keywords
creator: Jie Shen
creator: Lianpeng Xue
creator: Qifeng Zhang
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 661
ending page: 695
doi: 10.4208/cicp.OA-2025-0108

Deep Learning-Based Reduced-Order Methods for Fast Transient Dynamics
type: item
description: In recent years, large-scale numerical simulations played an essential role in estimating the effects of explosion events in urban environments, for the purpose of ensuring the security and safety of cities. Such simulations are computationally expensive and, often, the time taken for one single computation is large and does not permit parametric studies. The aim of this work is therefore to facilitate real-time and multi-query calculations by employing a non-intrusive Reduced Order Method (ROM). We propose a deep learning-based (DL) ROM scheme able to deal with fast transient dynamics. In the case of blast waves, the parametrised PDEs are time-dependent and non-linear. For such problems, the Proper Orthogonal Decomposition (POD), which relies on a linear superposition of modes, cannot approximate the solutions efficiently. The piecewise POD-DL scheme developed here is a local ROM based on time-domain partitioning and a first dimensionality reduction obtained through the POD. Autoencoders are used as a second and non-linear dimensionality reduction. The latent space obtained is then reconstructed from the time and parameter space through deep forward neural networks. The proposed scheme is applied to an example consisting of a blast wave propagating in air and impacting on the outside of a building. The efficiency of the deep learning-based ROM in approximating the time-dependent pressure field is shown.
subject: Articles
  topic: section
subject: Deep learning
  topic: keywords
subject: autoencoders
  topic: keywords
subject: artificial neural network
  topic: keywords
subject: reduced-order modelling
  topic: keywords
subject: nonlinear PDEs
  topic: keywords
subject: blast wave
  topic: keywords
subject: explosion event
  topic: keywords
creator: Martina Cracco
creator: Giovanni Stabile
creator: Andrea Lario
creator: Armin Sheidani
creator: Martin Larcher
creator: Folco Casadei
creator: Georgios Valsamos
creator: Gianluigi Rozza
rights: Copyright (c) 2026 Communications in Computational Physics
date: 2026-01-03
publication date: 2026-01-03
starting page: 635
ending page: 660
doi: 10.4208/cicp.OA-2024-0168

A New Approximate Method to Mean Field Stochastic Differential Equation with One-Sided Lipschitz Drift Coefficient
type: item
description: In this paper, we construct a numerical method to approximate a class of mean field stochastic differential equations whose drift coefficient satisfies a one-sided Lipschitz condition. The key idea is to independently approximate the distribution via solving a nonlinear Fokker-Planck equation which governs the density function of the solution. Meanwhile, based on the approximate density function, we construct an approximate stochastic differential equation to approach the mean field one. Then, we also apply a truncated Euler-Maruyama method to achieve discretization and derive error estimates. Finally, we present several numerical experiments to illustrate our theoretical analysis.
subject: Articles
  topic: section
subject: Mean field stochastic differential equation
  topic: keywords
subject: nonlinear Fokker-Planck equation
  topic: keywords
subject: truncated Euler-Maruyama method
  topic: keywords
subject: error estimate
  topic: keywords
creator: Jinhui Zhou
creator: Yongkui Zou
creator: Boyu Wang
rights: Copyright (c) 2025 Communications in Computational Physics
date: 2026-02-15
publication date: 2026-02-15
starting page: 615
ending page: 634
doi: 10.4208/cicp.OA-2025-0126

Uncertainty Quantification of Phase Transition Problems with an Injection Boundary
type: item
description: We develop an enthalpy-based modeling and computational framework to quantify uncertainty in Stefan problems with an injection boundary. Inspired by airfoil icing studies, we consider a system featuring an injection boundary inducing domain changes and a free boundary separating phases, resulting in two types of moving boundaries. Our proposed enthalpy-based formulation seamlessly integrates thermal diffusion across the domain with energy fluxes at the boundaries, addressing a modified injection condition for boundary movement. Uncertainty then stems from random variations in the injection boundary. The primary focus of our Uncertainty Quantification (UQ) centers on investigating the effects of uncertainty on free boundary propagation. Through mapping to a reference domain, we derive an enthalpy-based numerical scheme tailored to the transformed coordinate system, facilitating a simple and efficient simulation. Numerical and UQ studies in one and two dimensions validate the proposed model and the extended enthalpy method. They offer intriguing insights into ice accretion and other multiphysics processes involving phase transitions.
subject: Articles
  topic: section
subject: Uncertainty quantification
  topic: keywords
subject: phase transition
  topic: keywords
subject: injection boundary
  topic: keywords
subject: enthalpy method
  topic: keywords
creator: Zhenyi Zhang
creator: Shengbo Ma
creator: Zhennan Zhou
rights: Copyright (c) 2025 Communications in Computational Physics
date: 2026-02-15
publication date: 2026-02-15
starting page: 578
ending page: 614
doi: 10.4208/cicp.OA-2024-0036

Adaptive Neural Network Basis Methods for Partial Differential Equations with Low-Regular Solutions
type: item
description: This paper aims to devise an adaptive neural network basis method for numerically solving a second-order semi-linear partial differential equation with low-regular solutions in two/three dimensions. The method is obtained by combining basis functions from a class of shallow neural networks and the resulting multi-scale analogues, a residual strategy in adaptive methods and the non-overlapping domain decomposition method. Firstly, based on the solution residual, the domain $\Omega$ is iteratively decomposed and eventually partitioned into $K+1$ non-overlapping subdomains, denoted respectively as $\{\Omega_k\}_{k=0}^K$, where the exact solution is smooth on subdomain $\Omega_0$ and low-regular on subdomain $\Omega_k$ $(1\leq k\leq K)$. Secondly, the low-regular solutions on different subdomains $\Omega_k$ $(1\leq k\leq K)$ are approximated by neural networks with different scales, while the smooth solution on subdomain $\Omega_0$ is approximated by the initialized neural network. Thirdly, we determine the undetermined coefficients by solving the linear least squares problems directly or the nonlinear least squares problem via the Gauss-Newton method. The proposed method can be extended to multi-level case naturally. Finally, we use this adaptive method for several peak problems in two/three dimensions to show its high-efficient computational performance.
subject: Articles
  topic: section
subject: Neural network basis functions
  topic: keywords
subject: domain decomposition
  topic: keywords
subject: low-regular solutions
  topic: keywords
subject: least squares problem
  topic: keywords
subject: adaptive method
  topic: keywords
creator: Jianguo Huang
creator: Haohao Wu
creator: Tao Zhou
rights: Copyright (c) 2025 Communications in Computational Physics
date: 2026-02-15
publication date: 2026-02-15
starting page: 553
ending page: 577
doi: 10.4208/cicp.OA-2024-0310

A Unifying Moving Mesh Method for Curves, Surfaces, and Domains Based on Mesh Equidistribution and Alignment
type: item
description: A unifying moving mesh method is developed for general $m$-dimensional geometric objects in $d$-dimensions ($d \geq 1$ and $1 \leq m \leq d$) including curves, surfaces, and domains. The method is based on mesh equidistribution and alignment and does not require the availability of an analytical parametric representation of the underlying geometric object. Mathematical characterizations of shape and size of $m$-simplexes and properties of corresponding edge matrices and affine mappings are derived. The equidistribution and alignment conditions are presented in a unifying form for $m$-simplicial meshes. The equation for mesh movement is defined based on the moving mesh PDE approach, and suitable projection of the nodal mesh velocities is employed to ensure the mesh points are not moved out of the underlying geometric object. The analytical expression for the mesh velocities is obtained in a compact matrix form. The nonsingularity of moving meshes is proved. Numerical results for curves $(m=1)$ and surfaces $(m=2)$ in two and three dimensions are presented to demonstrate the ability of the developed method to move mesh points without causing singularity and control their concentration.
subject: Articles
  topic: section
subject: Unifying method for mesh movement
  topic: keywords
subject: moving mesh PDE
  topic: keywords
subject: mesh nonsingularity
  topic: keywords
subject: equidistribution
  topic: keywords
subject: alignment
  topic: keywords
creator: Min Zhang
creator: Weizhang Huang
rights: Copyright (c) 2025 Communications in Computational Physics
date: 2026-02-15
publication date: 2026-02-15
starting page: 513
ending page: 552
doi: 10.4208/cicp.OA-2025-0004