Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0233
Wenqiang Xiao, Bo Gong, Junshan Lin null, Jiguang Sun
{"title":"Band Structure Calculations of Dispersive Photonic Crystals in 3D using Holomorphic Operator Functions","authors":"Wenqiang Xiao, Bo Gong, Junshan Lin null, Jiguang Sun","doi":"10.4208/cicp.oa-2022-0233","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0233","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49652728","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}
. Constitutive modeling of heterogeneous hyperelastic materials is still a challenge due to their complex and variable microstructures. We propose a multiscale data-driven approach with a hierarchical learning strategy for the discovery of a generic physics-constrained anisotropic constitutive model for the heterogeneous hyperelastic materials. Based on the sparse multiscale experimental data, the constitutive artificial neural networks for hyperelastic component phases containing composite interfaces are established by the particle swarm optimization algorithm. A microscopic finite element coupled constitutive artificial neural networks solver is introduced to obtain the homogenized stress-stretch relation of heterogeneous materials with different microstructures. And a dense stress-stretch relation dataset is generated by training a neural network through the FE results. Further, a generic invariant representation of strain energy function (SEF) is proposed with a parameter set being implicitly expressed by artificial neural networks (SANN), which describes the hyperelastic properties of heterogeneous materials with different microstructures. A convexity constraint is imposed on the SEF to ensure that the multiscale constitutive model is physically relevant
{"title":"Learning Invariant Representation of Multiscale Hyperelastic Constitutive Law from Sparse Experimental Data","authors":"Rui He, Junzhi Cui, Zihao Yang, Jieqiong Zhang null, Xiaofei Guan","doi":"10.4208/cicp.oa-2023-0098","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0098","url":null,"abstract":". Constitutive modeling of heterogeneous hyperelastic materials is still a challenge due to their complex and variable microstructures. We propose a multiscale data-driven approach with a hierarchical learning strategy for the discovery of a generic physics-constrained anisotropic constitutive model for the heterogeneous hyperelastic materials. Based on the sparse multiscale experimental data, the constitutive artificial neural networks for hyperelastic component phases containing composite interfaces are established by the particle swarm optimization algorithm. A microscopic finite element coupled constitutive artificial neural networks solver is introduced to obtain the homogenized stress-stretch relation of heterogeneous materials with different microstructures. And a dense stress-stretch relation dataset is generated by training a neural network through the FE results. Further, a generic invariant representation of strain energy function (SEF) is proposed with a parameter set being implicitly expressed by artificial neural networks (SANN), which describes the hyperelastic properties of heterogeneous materials with different microstructures. A convexity constraint is imposed on the SEF to ensure that the multiscale constitutive model is physically relevant","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41985774","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0270
Yifei Wan null, Yinhua Xia
. For steady Euler equations in complex boundary domains, high-order shock-capturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical boundaries on Cartesian grids to achieve uniform high-order accuracy. In this paper, we utilize a fifth-order finite difference hybrid WENO scheme to simulate steady Euler equations, and the same fifth-order WENO extrapolation methods are developed to handle the curved boundary. The values of the ghost points outside the physical boundary can be obtained by applying WENO extrapolation near the boundary, involving normal derivatives acquired by the simplified inverse Lax-Wendroff procedure. Both equivalent expressions involving curvature and numerical differentiation are utilized to transform the tangential derivatives along the curved solid wall boundary. This hybrid WENO scheme is robust for steady-state convergence and maintains high-order accuracy in the smooth region even with the solid wall boundary condition. Besides, the essentially non-oscillation property is achieved. The numerical spectral analysis also shows that this hybrid WENO scheme has low dispersion and dissipation errors. Numerical examples are presented to validate the high-order accuracy and robust performance of the hybrid scheme for steady Euler equations in curved domains with Cartesian grids.
{"title":"A Hybrid WENO Scheme for Steady Euler Equations in Curved Geometries on Cartesian Grids","authors":"Yifei Wan null, Yinhua Xia","doi":"10.4208/cicp.oa-2022-0270","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0270","url":null,"abstract":". For steady Euler equations in complex boundary domains, high-order shock-capturing schemes usually suffer not only from the difficulty of steady-state convergence but also from the problem of dealing with physical boundaries on Cartesian grids to achieve uniform high-order accuracy. In this paper, we utilize a fifth-order finite difference hybrid WENO scheme to simulate steady Euler equations, and the same fifth-order WENO extrapolation methods are developed to handle the curved boundary. The values of the ghost points outside the physical boundary can be obtained by applying WENO extrapolation near the boundary, involving normal derivatives acquired by the simplified inverse Lax-Wendroff procedure. Both equivalent expressions involving curvature and numerical differentiation are utilized to transform the tangential derivatives along the curved solid wall boundary. This hybrid WENO scheme is robust for steady-state convergence and maintains high-order accuracy in the smooth region even with the solid wall boundary condition. Besides, the essentially non-oscillation property is achieved. The numerical spectral analysis also shows that this hybrid WENO scheme has low dispersion and dissipation errors. Numerical examples are presented to validate the high-order accuracy and robust performance of the hybrid scheme for steady Euler equations in curved domains with Cartesian grids.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"45229733","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2023-0068
Jianan Zeng, W. S. Null, L. Wu
{"title":"General Synthetic Iterative Scheme for Unsteady Rarefied Gas Flows","authors":"Jianan Zeng, W. S. Null, L. Wu","doi":"10.4208/cicp.oa-2023-0068","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0068","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48220126","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2023-0132
Yangyang Cao, Alexander Kurganov null, Yongle Liu
{"title":"Flux Globalization Based Well-Balanced Path-Conservative Central-Upwind Scheme for the Thermal Rotating Shallow Water Equations","authors":"Yangyang Cao, Alexander Kurganov null, Yongle Liu","doi":"10.4208/cicp.oa-2023-0132","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0132","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"70 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135194270","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 traditional stochastic homogenization method can obtain homogenized solutions of elliptic problems with stationary random coefficients. However, many random composite materials in scientific and engineering computing do not satisfy the stationary assumption. To overcome the difficulty, we propose a normalizing field flow induced two-stage stochastic homogenization method to efficiently solve the random elliptic problem with non-stationary coefficients. By applying the two-stage stochastic homogenization method, the original elliptic equation with random and fast oscillatory coefficients is approximated as an equivalent elliptic equation, where the equivalent coefficients are obtained by solving a set of cell problems. Without the stationary assumption, the number of cell problems is large and the corresponding computational cost is high. To improve the efficiency, we apply the normalizing field flow model to learn a reference Gaussian field for the random equivalent coefficients based on a small amount of data, which is obtained by solving the cell problems with the finite element method. Numerical results demonstrate that the newly proposed method is efficient and accurate in tackling high dimensional partial differential equations in composite materials with complex random microstructures
{"title":"A Normalizing Field Flow Induced Two-Stage Stochastic Homogenization Method for Random Composite Materials","authors":"Zihao Yang, Xintong Wang, Xiaofei Guan, Jizu Huang null, Xixin Wu","doi":"10.4208/cicp.oa-2023-0007","DOIUrl":"https://doi.org/10.4208/cicp.oa-2023-0007","url":null,"abstract":". The traditional stochastic homogenization method can obtain homogenized solutions of elliptic problems with stationary random coefficients. However, many random composite materials in scientific and engineering computing do not satisfy the stationary assumption. To overcome the difficulty, we propose a normalizing field flow induced two-stage stochastic homogenization method to efficiently solve the random elliptic problem with non-stationary coefficients. By applying the two-stage stochastic homogenization method, the original elliptic equation with random and fast oscillatory coefficients is approximated as an equivalent elliptic equation, where the equivalent coefficients are obtained by solving a set of cell problems. Without the stationary assumption, the number of cell problems is large and the corresponding computational cost is high. To improve the efficiency, we apply the normalizing field flow model to learn a reference Gaussian field for the random equivalent coefficients based on a small amount of data, which is obtained by solving the cell problems with the finite element method. Numerical results demonstrate that the newly proposed method is efficient and accurate in tackling high dimensional partial differential equations in composite materials with complex random microstructures","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"21 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135145211","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0236
Yaning Xie, Shuwang Li null, W. Ying
{"title":"A Fourth-Order Kernel-Free Boundary Integral Method for Interface Problems","authors":"Yaning Xie, Shuwang Li null, W. Ying","doi":"10.4208/cicp.oa-2022-0236","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0236","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42983991","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0253
Hojun You, Juhyun Kim null, Chongam Kim
{"title":"Implicit Quadrature-Free Direct Reconstruction Method for Efficient Scale-Resolving Simulations","authors":"Hojun You, Juhyun Kim null, Chongam Kim","doi":"10.4208/cicp.oa-2022-0253","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0253","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"1 1","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41409425","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0319
Xinjiang Chen null, A. Kurganov
{"title":"A Well-Balanced Partial Relaxation Scheme for the Two-Dimensional Saint-Venant System","authors":"Xinjiang Chen null, A. Kurganov","doi":"10.4208/cicp.oa-2022-0319","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0319","url":null,"abstract":"","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":" ","pages":""},"PeriodicalIF":3.7,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"48867170","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}
Pub Date : 2023-06-01DOI: 10.4208/cicp.oa-2022-0104
Meiling Zhao, Jiahui He null, Liqun Wang
. In this paper, the electromagnetic scattering from overfilled cavities with inhomogeneous anisotropic media is investigated. To solve the scattering problem, a Petrov-Galerkin finite element interface method on non-body-fitted grids is presented. We reduce the infinite domain of scattering to a bounded domain problem by introducing a transparent boundary condition. The level set function is used to capture complex boundary and interface geometry that is not aligned with the mesh. Non-body-fitted grids allow us to save computational costs during mesh generation and significantly reduce the amount of computer memory required. The solution is built by connecting two linear polynomials across the interfaces to satisfy the jump conditions. The proposed method can handle matrix coefficients produced by permittivity and permeability tensors of anisotropic media. The final linear system is sparse, making it more suitable for most iterative methods. Numerical experiments show that the proposed method has good convergence and realizability. Meanwhile, we discover that the absorbing properties of anisotropic media clearly and positively influence the reduction of radar cross section. It has also been demonstrated that the method can achieve second-order accuracy.
{"title":"Numerical Solutions of the Electromagnetic Scattering by Overfilled Cavities with Inhomogeneous Anisotropic Media","authors":"Meiling Zhao, Jiahui He null, Liqun Wang","doi":"10.4208/cicp.oa-2022-0104","DOIUrl":"https://doi.org/10.4208/cicp.oa-2022-0104","url":null,"abstract":". In this paper, the electromagnetic scattering from overfilled cavities with inhomogeneous anisotropic media is investigated. To solve the scattering problem, a Petrov-Galerkin finite element interface method on non-body-fitted grids is presented. We reduce the infinite domain of scattering to a bounded domain problem by introducing a transparent boundary condition. The level set function is used to capture complex boundary and interface geometry that is not aligned with the mesh. Non-body-fitted grids allow us to save computational costs during mesh generation and significantly reduce the amount of computer memory required. The solution is built by connecting two linear polynomials across the interfaces to satisfy the jump conditions. The proposed method can handle matrix coefficients produced by permittivity and permeability tensors of anisotropic media. The final linear system is sparse, making it more suitable for most iterative methods. Numerical experiments show that the proposed method has good convergence and realizability. Meanwhile, we discover that the absorbing properties of anisotropic media clearly and positively influence the reduction of radar cross section. It has also been demonstrated that the method can achieve second-order accuracy.","PeriodicalId":50661,"journal":{"name":"Communications in Computational Physics","volume":"20 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"2023-06-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"135144002","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}