Search papers, labs, and topics across Lattice.
This paper introduces a structure-oriented randomized neural network framework (SO-RaNN) designed to solve the Poisson-Nernst-Planck (PNP) and Poisson-Nernst-Planck-Navier-Stokes (PNP-NS) systems. The approach employs iterative solutions to decoupled linearized subproblems, incorporating mechanisms for positivity enforcement and mass preservation, alongside a structure-preserving method for velocity fields that ensures incompressibility. Key results include theoretical convergence guarantees and numerical experiments demonstrating high accuracy and adherence to physical constraints, such as positivity and mass matching.
SO-RaNN achieves exact mass matching and positivity in PNP systems while ensuring incompressibility in velocity fields, setting a new standard for accuracy in neural network-based simulations.
We develop a structure-oriented randomized neural network framework, termed SO-RaNN, for the Poisson-Nernst-Planck (PNP) system and the Poisson-Nernst-Planck-Navier-Stokes (PNP-NS) system. The decoupled linearized subproblems are solved iteratively by randomized neural networks in a space-time framework. For the concentration variables, a pointwise cut-off is used to enforce positivity at the value level, and discrete mass-scaling factors are computed at selected correction instants and interpolated in time, so as to ensure exact mass matching at those instants and to promote approximate mass preservation between them. To introduce an auxiliary discrete dissipation mechanism, we further employ an SAV-type post-processing correction, which yields monotonicity of the SAV auxiliary variable under the ideal SAV update. For the PNP-NS system, a structure-preserving randomized neural network (SP-RaNN) is used for the velocity field, so that the velocity approximation satisfies the incompressibility constraint pointwise by construction. On the theoretical side, we derive residual-based estimates for the raw, uncorrected RaNN solvers of the linearized subproblems, formulate a conditional local-in-time convergence result for the raw outer Picard iteration of the PNP system, and analyze the value-level positivity correction together with the mass-correction and SAV post-processing steps. For the PNP-NS system, we establish an approximation result for the SP-RaNN space and provide a conditional error statement for the corresponding linearized Oseen-type problem. Numerical experiments demonstrate approximation accuracy in the source-driven manufactured tests and illustrate the intended value-level positivity correction, selected-time mass matching, computed free-energy curves based on the final gauge-fixed potential, and divergence-free approximation in benchmark tests.
