{"title":"Unsupervised neural-network solvers for multi-material Riemann problems","authors":"Liang Xu , Ziyan Liu , Yiwei Feng , Tiegang Liu","doi":"10.1016/j.cpc.2024.109470","DOIUrl":null,"url":null,"abstract":"<div><div>Machine learning has the potential to provide a non-traditional and feasible approach for solving Riemann problems to model the coupling effects of multi-material flows. However, most recent research on predicting Riemann solutions with neural networks is limited to addressing single-material flows and featured as the supervised learning, or is limited to solving specific problems and difficult to apply to a wide range of initial conditions. In this work, we explore physics-constrained neural networks, termed PCNN-RS, as multi-material Riemann solvers without any labeled data. Based on the frame of a general neural network, physics-constrained functions that conform to the shock/rarefaction relationships between initial states and interfacial states are constructed after the output layer, transforming the unlabeled output into a theoretically zero-valued functional form. This allows training learning models with standard loss functions solely using input data. The interfacial pressure of multi-material Riemann problem is predicted using the surrogate model, and other interfacial states can be directly derived through simple calculations. In addition, the basic principle of scaling of initial conditions and Riemann solutions with general equations of state is established theoretically. Based on this property, a transformation of input and output data is proposed to enhance the wide applicability of the Riemann-solver surrogate model. Furthermore, an optimization of samples is presented to reduce the training dataset and shorten the training time. The PCNN-RS is able to make accurate predictions, even when utilizing a compact neural network architecture with fewer neurons, and it is easily applied to the ghost-fluid-based sharp interface methods. It possesses the ability to simulate various interface evolutions for the interaction between two materials.</div></div>","PeriodicalId":285,"journal":{"name":"Computer Physics Communications","volume":"308 ","pages":"Article 109470"},"PeriodicalIF":7.2000,"publicationDate":"2024-12-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer Physics Communications","FirstCategoryId":"101","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S001046552400393X","RegionNum":2,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Machine learning has the potential to provide a non-traditional and feasible approach for solving Riemann problems to model the coupling effects of multi-material flows. However, most recent research on predicting Riemann solutions with neural networks is limited to addressing single-material flows and featured as the supervised learning, or is limited to solving specific problems and difficult to apply to a wide range of initial conditions. In this work, we explore physics-constrained neural networks, termed PCNN-RS, as multi-material Riemann solvers without any labeled data. Based on the frame of a general neural network, physics-constrained functions that conform to the shock/rarefaction relationships between initial states and interfacial states are constructed after the output layer, transforming the unlabeled output into a theoretically zero-valued functional form. This allows training learning models with standard loss functions solely using input data. The interfacial pressure of multi-material Riemann problem is predicted using the surrogate model, and other interfacial states can be directly derived through simple calculations. In addition, the basic principle of scaling of initial conditions and Riemann solutions with general equations of state is established theoretically. Based on this property, a transformation of input and output data is proposed to enhance the wide applicability of the Riemann-solver surrogate model. Furthermore, an optimization of samples is presented to reduce the training dataset and shorten the training time. The PCNN-RS is able to make accurate predictions, even when utilizing a compact neural network architecture with fewer neurons, and it is easily applied to the ghost-fluid-based sharp interface methods. It possesses the ability to simulate various interface evolutions for the interaction between two materials.
期刊介绍:
The focus of CPC is on contemporary computational methods and techniques and their implementation, the effectiveness of which will normally be evidenced by the author(s) within the context of a substantive problem in physics. Within this setting CPC publishes two types of paper.
Computer Programs in Physics (CPiP)
These papers describe significant computer programs to be archived in the CPC Program Library which is held in the Mendeley Data repository. The submitted software must be covered by an approved open source licence. Papers and associated computer programs that address a problem of contemporary interest in physics that cannot be solved by current software are particularly encouraged.
Computational Physics Papers (CP)
These are research papers in, but are not limited to, the following themes across computational physics and related disciplines.
mathematical and numerical methods and algorithms;
computational models including those associated with the design, control and analysis of experiments; and
algebraic computation.
Each will normally include software implementation and performance details. The software implementation should, ideally, be available via GitHub, Zenodo or an institutional repository.In addition, research papers on the impact of advanced computer architecture and special purpose computers on computing in the physical sciences and software topics related to, and of importance in, the physical sciences may be considered.