The conductive interconnect structure that connects the electrical functional devices is an important micro-nano structure in stretchable electronics. Given the reliance of numerous devices on steady electrical currents for operation, stretchable electronics would benefit from interconnects with minimal resistance variation during deformation. This paper proposes a topology optimization method for the design of stretchable interconnect structures with stable resistance under large deformation. In the proposed method, an equal material method considering geometrically nonlinear and electromechanical coupling effects is developed to evaluate the resistance of a deformed structure. Besides, a new connectivity control method is proposed to ensure the connectivity between the inlet and outlet by making full use of the electrical problem itself. To achieve the design goal of connected interconnect structures with negligible resistance fluctuation during stretching, a topology optimization formulation is established, and the corresponding sensitivity is also analytically derived. Several numerical examples show that the proposed method is capable of computationally and intelligently generating stretchable structures with extremely small variations in resistance during stretching.
{"title":"Nonlinear electromechanical topology optimization method for stretchable electronic interconnect structures","authors":"Yunfeng Luo, Shiyuan Qu, Shutian Liu, YongAn Huang","doi":"10.1007/s00366-024-01996-y","DOIUrl":"https://doi.org/10.1007/s00366-024-01996-y","url":null,"abstract":"<p>The conductive interconnect structure that connects the electrical functional devices is an important micro-nano structure in stretchable electronics. Given the reliance of numerous devices on steady electrical currents for operation, stretchable electronics would benefit from interconnects with minimal resistance variation during deformation. This paper proposes a topology optimization method for the design of stretchable interconnect structures with stable resistance under large deformation. In the proposed method, an equal material method considering geometrically nonlinear and electromechanical coupling effects is developed to evaluate the resistance of a deformed structure. Besides, a new connectivity control method is proposed to ensure the connectivity between the inlet and outlet by making full use of the electrical problem itself. To achieve the design goal of connected interconnect structures with negligible resistance fluctuation during stretching, a topology optimization formulation is established, and the corresponding sensitivity is also analytically derived. Several numerical examples show that the proposed method is capable of computationally and intelligently generating stretchable structures with extremely small variations in resistance during stretching.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"105 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140932114","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}
Pub Date : 2024-05-08DOI: 10.1007/s00366-024-01987-z
Amirhossein Arzani, Lingxiao Yuan, Pania Newell, Bei Wang
Although deep learning has achieved remarkable success in various scientific machine learning applications, its opaque nature poses concerns regarding interpretability and generalization capabilities beyond the training data. Interpretability is crucial and often desired in modeling physical systems. Moreover, acquiring extensive datasets that encompass the entire range of input features is challenging in many physics-based learning tasks, leading to increased errors when encountering out-of-distribution (OOD) data. In this work, motivated by the field of functional data analysis (FDA), we propose generalized functional linear models as an interpretable surrogate for a trained deep learning model. We demonstrate that our model could be trained either based on a trained neural network (post-hoc interpretation) or directly from training data (interpretable operator learning). A library of generalized functional linear models with different kernel functions is considered and sparse regression is used to discover an interpretable surrogate model that could be analytically presented. We present test cases in solid mechanics, fluid mechanics, and transport. Our results demonstrate that our model can achieve comparable accuracy to deep learning and can improve OOD generalization while providing more transparency and interpretability. Our study underscores the significance of interpretable representation in scientific machine learning and showcases the potential of functional linear models as a tool for interpreting and generalizing deep learning.
{"title":"Interpreting and generalizing deep learning in physics-based problems with functional linear models","authors":"Amirhossein Arzani, Lingxiao Yuan, Pania Newell, Bei Wang","doi":"10.1007/s00366-024-01987-z","DOIUrl":"https://doi.org/10.1007/s00366-024-01987-z","url":null,"abstract":"<p>Although deep learning has achieved remarkable success in various scientific machine learning applications, its opaque nature poses concerns regarding interpretability and generalization capabilities beyond the training data. Interpretability is crucial and often desired in modeling physical systems. Moreover, acquiring extensive datasets that encompass the entire range of input features is challenging in many physics-based learning tasks, leading to increased errors when encountering out-of-distribution (OOD) data. In this work, motivated by the field of functional data analysis (FDA), we propose generalized functional linear models as an interpretable surrogate for a trained deep learning model. We demonstrate that our model could be trained either based on a trained neural network (post-hoc interpretation) or directly from training data (interpretable operator learning). A library of generalized functional linear models with different kernel functions is considered and sparse regression is used to discover an interpretable surrogate model that could be analytically presented. We present test cases in solid mechanics, fluid mechanics, and transport. Our results demonstrate that our model can achieve comparable accuracy to deep learning and can improve OOD generalization while providing more transparency and interpretability. Our study underscores the significance of interpretable representation in scientific machine learning and showcases the potential of functional linear models as a tool for interpreting and generalizing deep learning.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"43 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140932020","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}
Pub Date : 2024-05-08DOI: 10.1007/s00366-024-01986-0
Puyang Gao, Xiaolin Hu
In this paper, we develop a semi-implicit partitioned finite element and discontinuous Galerkin method for the non-isothermal non-Newtonian fluid structure interaction (NNFSI) problem within the arbitrary Lagrangian–Eulerian (ALE) framework. The structure is composed of the elastic solid material. The entire mathematical model consists of the governing equations of the non-Newtonian fluid and the structure, as well as the boundary conditions on the contacting interface. The rheological behavior of non-Newtonian fluid is described according to the power law constitutive equation. The whole system is split into several sub-equations and then appropriate finite element method or discontinuous Galerkin method is employed for the spatial discretizations of them. As for the deformation of the structure and the change of the fluid area and computational mesh, we employ the moving mesh technique to handle them. The problem involving a hot flexible rod fixed on the hot bottom of an irregular pipe is fully investigated. The influences of the fluid inlet velocity and the behavior of the fluid on the deformation of the rod and the temperature distribution are all analyzed.
{"title":"The development of an ALE finite element and discontinuous Galerkin method for the non-isothermal non-Newtonian FSI problem","authors":"Puyang Gao, Xiaolin Hu","doi":"10.1007/s00366-024-01986-0","DOIUrl":"https://doi.org/10.1007/s00366-024-01986-0","url":null,"abstract":"<p>In this paper, we develop a semi-implicit partitioned finite element and discontinuous Galerkin method for the non-isothermal non-Newtonian fluid structure interaction (NNFSI) problem within the arbitrary Lagrangian–Eulerian (ALE) framework. The structure is composed of the elastic solid material. The entire mathematical model consists of the governing equations of the non-Newtonian fluid and the structure, as well as the boundary conditions on the contacting interface. The rheological behavior of non-Newtonian fluid is described according to the power law constitutive equation. The whole system is split into several sub-equations and then appropriate finite element method or discontinuous Galerkin method is employed for the spatial discretizations of them. As for the deformation of the structure and the change of the fluid area and computational mesh, we employ the moving mesh technique to handle them. The problem involving a hot flexible rod fixed on the hot bottom of an irregular pipe is fully investigated. The influences of the fluid inlet velocity and the behavior of the fluid on the deformation of the rod and the temperature distribution are all analyzed.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"62 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140932058","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}
This paper introduces the application of Physics-Informed Neural Network (PINN), a novel class of scientific machine learning techniques, for analyzing the static bending response of nanobeams as essential structural elements in micro/nanoelectromechanical systems, including nanoprobes, atomic force microscope sensors, nanoswitches, nanoactuators, and nanoscale biosensors on a three-parameter nonlinear elastic foundation. The study combines Euler–Bernoulli beam theory and Eringen’s nonlocal continuum theory to derive the governing differential equation using the minimum total potential energy principle. PINN is utilized for approximating the differential equation solution and identifying the nanobeam’s nonlocal parameter through an inverse problem with measurement data. The loss function incorporates terms representing the initial and boundary conditions, along with the differential equation residual at specific points in the domain and boundary. The research demonstrates PINN’s efficacy in analyzing nanobeam behavior on nonlinear elastic foundations, providing valuable insights into responses under different loading and boundary conditions. The proposed approach's accuracy and efficiency are validated through comparisons with existing literature. Additionally, the study investigates the effects of activation functions, collocation points’ number and distribution, nonlocal parameter, foundation stiffness coefficients, loading types, and various boundary conditions on nanobeam bending behavior.
{"title":"PINN-based forward and inverse bending analysis of nanobeams on a three-parameter nonlinear elastic foundation including hardening and softening effect using nonlocal elasticity theory","authors":"Omid Kianian, Saeid Sarrami, Bashir Movahedian, Mojtaba Azhari","doi":"10.1007/s00366-024-01985-1","DOIUrl":"https://doi.org/10.1007/s00366-024-01985-1","url":null,"abstract":"<p>This paper introduces the application of Physics-Informed Neural Network (PINN), a novel class of scientific machine learning techniques, for analyzing the static bending response of nanobeams as essential structural elements in micro/nanoelectromechanical systems, including nanoprobes, atomic force microscope sensors, nanoswitches, nanoactuators, and nanoscale biosensors on a three-parameter nonlinear elastic foundation. The study combines Euler–Bernoulli beam theory and Eringen’s nonlocal continuum theory to derive the governing differential equation using the minimum total potential energy principle. PINN is utilized for approximating the differential equation solution and identifying the nanobeam’s nonlocal parameter through an inverse problem with measurement data. The loss function incorporates terms representing the initial and boundary conditions, along with the differential equation residual at specific points in the domain and boundary. The research demonstrates PINN’s efficacy in analyzing nanobeam behavior on nonlinear elastic foundations, providing valuable insights into responses under different loading and boundary conditions. The proposed approach's accuracy and efficiency are validated through comparisons with existing literature. Additionally, the study investigates the effects of activation functions, collocation points’ number and distribution, nonlocal parameter, foundation stiffness coefficients, loading types, and various boundary conditions on nanobeam bending behavior.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"47 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140887994","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}
Pub Date : 2024-05-07DOI: 10.1007/s00366-024-01978-0
Mohammad Naqib Rahimi, Georgios Moutsanidis
We introduce a novel immersed-like numerical framework that combines isogeometric analysis with smoothed particle hydrodynamics for simulating air-blast–structure interaction. The solid domain is represented by a Lagrangian point cloud, which is immersed into a background Eulerian fluid domain. The smoothed particle hydrodynamics framework is employed to solve the equations of motion of the solid point cloud, whereas isogeometric analysis is used for the fluid mechanics equations on the background domain. The coupling strategy relies on a penalty-based volumetric coupling scheme that penalizes the velocity difference between the two domains, and involves a minimal amount of modification to existing codes, resulting in a straightforward implementation. The immersed nature of the proposed approach, combined with volumetric coupling, eliminates the need for explicit tracking of fluid–structure interfaces and imposes no limitations on solid domain motion and topology. Ample mathematical details are provided, and the proposed method is verified and validated against established numerical tools and experimental studies. The results affirm the method’s accuracy, robustness, and ease with which it seamlessly integrates two distinct computational techniques.
{"title":"IGA-SPH: coupling isogeometric analysis with smoothed particle hydrodynamics for air-blast–structure interaction","authors":"Mohammad Naqib Rahimi, Georgios Moutsanidis","doi":"10.1007/s00366-024-01978-0","DOIUrl":"https://doi.org/10.1007/s00366-024-01978-0","url":null,"abstract":"<p>We introduce a novel immersed-like numerical framework that combines isogeometric analysis with smoothed particle hydrodynamics for simulating air-blast–structure interaction. The solid domain is represented by a Lagrangian point cloud, which is immersed into a background Eulerian fluid domain. The smoothed particle hydrodynamics framework is employed to solve the equations of motion of the solid point cloud, whereas isogeometric analysis is used for the fluid mechanics equations on the background domain. The coupling strategy relies on a penalty-based volumetric coupling scheme that penalizes the velocity difference between the two domains, and involves a minimal amount of modification to existing codes, resulting in a straightforward implementation. The immersed nature of the proposed approach, combined with volumetric coupling, eliminates the need for explicit tracking of fluid–structure interfaces and imposes no limitations on solid domain motion and topology. Ample mathematical details are provided, and the proposed method is verified and validated against established numerical tools and experimental studies. The results affirm the method’s accuracy, robustness, and ease with which it seamlessly integrates two distinct computational techniques.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"30 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140888063","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}
Pub Date : 2024-05-06DOI: 10.1007/s00366-024-01981-5
Tsung-Yeh Hsieh, Tsung-Hui Huang
Physics informed neural network (PINN) frameworks have been developed as a powerful technique for solving partial differential equations (PDEs) with potential data integration. However, when applied to advection based PDEs, PINNs confront challenges such as parameter sensitivity in boundary condition enforcement and diminished learning capability due to an ill-conditioned system resulting from the strong advection. In this study, we present a multiscale stabilized PINN formulation with a weakly imposed boundary condition (WBC) method coupled with transfer learning that can robustly model the advection diffusion equation. To address key challenges, we use an advection-flux-decoupling technique to prescribe the Dirichlet boundary conditions, which rectifies the imbalanced training observed in PINN with conventional penalty and strong enforcement methods. A multiscale approach under the least squares functional form of PINN is developed that introduces a controllable stabilization term, which can be regarded as a special form of Sobolev training that augments the learning capacity. The efficacy of the proposed method is demonstrated through the resolution of a series of benchmark problems of forward modeling, and the outcomes affirm the potency of the methodology proposed.
{"title":"A multiscale stabilized physics informed neural networks with weakly imposed boundary conditions transfer learning method for modeling advection dominated flow","authors":"Tsung-Yeh Hsieh, Tsung-Hui Huang","doi":"10.1007/s00366-024-01981-5","DOIUrl":"https://doi.org/10.1007/s00366-024-01981-5","url":null,"abstract":"<p>Physics informed neural network (PINN) frameworks have been developed as a powerful technique for solving partial differential equations (PDEs) with potential data integration. However, when applied to advection based PDEs, PINNs confront challenges such as parameter sensitivity in boundary condition enforcement and diminished learning capability due to an ill-conditioned system resulting from the strong advection. In this study, we present a multiscale stabilized PINN formulation with a weakly imposed boundary condition (WBC) method coupled with transfer learning that can robustly model the advection diffusion equation. To address key challenges, we use an advection-flux-decoupling technique to prescribe the Dirichlet boundary conditions, which rectifies the imbalanced training observed in PINN with conventional penalty and strong enforcement methods. A multiscale approach under the least squares functional form of PINN is developed that introduces a controllable stabilization term, which can be regarded as a special form of Sobolev training that augments the learning capacity. The efficacy of the proposed method is demonstrated through the resolution of a series of benchmark problems of forward modeling, and the outcomes affirm the potency of the methodology proposed.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"17 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-05-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140887984","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}
Pub Date : 2024-04-29DOI: 10.1007/s00366-024-01980-6
Margarita Chasapi, Pablo Antolin, Annalisa Buffa
This contribution presents a model order reduction framework for real-time efficient solution of trimmed, multi-patch isogeometric Kirchhoff-Love shells. In several scenarios, such as design and shape optimization, multiple simulations need to be performed for a given set of physical or geometrical parameters. This step can be computationally expensive in particular for real world, practical applications. We are interested in geometrical parameters and take advantage of the flexibility of splines in representing complex geometries. In this case, the operators are geometry-dependent and generally depend on the parameters in a non-affine way. Moreover, the solutions obtained from trimmed domains may vary highly with respect to different values of the parameters. Therefore, we employ a local reduced basis method based on clustering techniques and the Discrete Empirical Interpolation Method to construct affine approximations and efficient reduced order models. In addition, we discuss the application of the reduction strategy to parametric shape optimization. Finally, we demonstrate the performance of the proposed framework to parameterized Kirchhoff-Love shells through benchmark tests on trimmed, multi-patch meshes including a complex geometry. The proposed approach is accurate and achieves a significant reduction of the online computational cost in comparison to the standard reduced basis method.
{"title":"Fast parametric analysis of trimmed multi-patch isogeometric Kirchhoff-Love shells using a local reduced basis method","authors":"Margarita Chasapi, Pablo Antolin, Annalisa Buffa","doi":"10.1007/s00366-024-01980-6","DOIUrl":"https://doi.org/10.1007/s00366-024-01980-6","url":null,"abstract":"<p>This contribution presents a model order reduction framework for real-time efficient solution of trimmed, multi-patch isogeometric Kirchhoff-Love shells. In several scenarios, such as design and shape optimization, multiple simulations need to be performed for a given set of physical or geometrical parameters. This step can be computationally expensive in particular for real world, practical applications. We are interested in geometrical parameters and take advantage of the flexibility of splines in representing complex geometries. In this case, the operators are geometry-dependent and generally depend on the parameters in a non-affine way. Moreover, the solutions obtained from trimmed domains may vary highly with respect to different values of the parameters. Therefore, we employ a local reduced basis method based on clustering techniques and the Discrete Empirical Interpolation Method to construct affine approximations and efficient reduced order models. In addition, we discuss the application of the reduction strategy to parametric shape optimization. Finally, we demonstrate the performance of the proposed framework to parameterized Kirchhoff-Love shells through benchmark tests on trimmed, multi-patch meshes including a complex geometry. The proposed approach is accurate and achieves a significant reduction of the online computational cost in comparison to the standard reduced basis method.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"75 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140811184","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}
Pub Date : 2024-04-29DOI: 10.1007/s00366-024-01977-1
Chein-Shan Liu, Chung-Lun Kuo
Because the complete set of Trefftz functions for the 3D biharmonic equation is not yet well established, a multiple-direction Trefftz method (MDTM) and an in-plane biharmonic functions method (IPBFM) are deduced in the paper. Inspired by the Trefftz method for the 2D biharmonic equation, a novel MDTM incorporates planar directors into the 2D like Trefftz functions to solve the 3D biharmonic equation. These functions being a series of biharmonic polynomials of different degree, automatically satisfying the 3D biharmonic equation, are taken as the bases to expand the solution. Then, we derive a quite large class solution of the 3D biharmonic equation in terms of 3D harmonic functions, and 2D biharmonic functions in three sub-planes. The 2D biharmonic functions are formulated as the Trefftz functions in terms of the polar coordinates for each sub-plane. Introducing a projective variable, we can obtain the projective type general solution for the 3D Laplace equation, which is used to generate the 3D Trefftz type harmonic functions. Several numerical examples confirm the efficiency and accuracy of the proposed MDTM and IPBFM.
{"title":"The Trefftz methods for 3D biharmonic equation using directors and in-plane biharmonic functions","authors":"Chein-Shan Liu, Chung-Lun Kuo","doi":"10.1007/s00366-024-01977-1","DOIUrl":"https://doi.org/10.1007/s00366-024-01977-1","url":null,"abstract":"<p>Because the complete set of Trefftz functions for the 3D biharmonic equation is not yet well established, a multiple-direction Trefftz method (MDTM) and an in-plane biharmonic functions method (IPBFM) are deduced in the paper. Inspired by the Trefftz method for the 2D biharmonic equation, a novel MDTM incorporates planar directors into the 2D like Trefftz functions to solve the 3D biharmonic equation. These functions being a series of biharmonic polynomials of different degree, automatically satisfying the 3D biharmonic equation, are taken as the bases to expand the solution. Then, we derive a quite large class solution of the 3D biharmonic equation in terms of 3D harmonic functions, and 2D biharmonic functions in three sub-planes. The 2D biharmonic functions are formulated as the Trefftz functions in terms of the polar coordinates for each sub-plane. Introducing a projective variable, we can obtain the projective type general solution for the 3D Laplace equation, which is used to generate the 3D Trefftz type harmonic functions. Several numerical examples confirm the efficiency and accuracy of the proposed MDTM and IPBFM.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"1 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-04-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140811303","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}
Pub Date : 2024-04-28DOI: 10.1007/s00366-024-01976-2
Mark S. Shephard, Jacob Merson, Onkar Sahni, Angel E. Castillo, Aditya Y. Joshi, Dhyanjyoti D. Nath, Usman Riaz, E. Seegyoung Seol, Cameron W. Smith, Chonglin Zhang, Mark W. Beall, Ottmar Klaas, Rocco Nastasia, Saurabh Tendulkar
As fusion simulation codes increasingly account for the full geometric complexity of magnetically confined fusion systems, a need arose to provide tailored unstructured mesh technologies to address the specific needs of fusion plasma simulation codes and their coupling to other physics simulation codes. This paper presents a high level overview of a set of unstructured mesh developments that include; methods to effectively employ manufacturing CAD models in the construction of quality analysis model geometries, specialized mesh generation and adaptation tools; an infrastructure to support parallel particle-in-cell calculation on unstructured meshes; and an infrastructure for coupling of massively parallel mesh-based fusion codes.
{"title":"Unstructured mesh tools for magnetically confined fusion system simulations","authors":"Mark S. Shephard, Jacob Merson, Onkar Sahni, Angel E. Castillo, Aditya Y. Joshi, Dhyanjyoti D. Nath, Usman Riaz, E. Seegyoung Seol, Cameron W. Smith, Chonglin Zhang, Mark W. Beall, Ottmar Klaas, Rocco Nastasia, Saurabh Tendulkar","doi":"10.1007/s00366-024-01976-2","DOIUrl":"https://doi.org/10.1007/s00366-024-01976-2","url":null,"abstract":"<p>As fusion simulation codes increasingly account for the full geometric complexity of magnetically confined fusion systems, a need arose to provide tailored unstructured mesh technologies to address the specific needs of fusion plasma simulation codes and their coupling to other physics simulation codes. This paper presents a high level overview of a set of unstructured mesh developments that include; methods to effectively employ manufacturing CAD models in the construction of quality analysis model geometries, specialized mesh generation and adaptation tools; an infrastructure to support parallel particle-in-cell calculation on unstructured meshes; and an infrastructure for coupling of massively parallel mesh-based fusion codes.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"6 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140811005","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}
Pub Date : 2024-04-28DOI: 10.1007/s00366-024-01983-3
Zhixin Zeng, Xiong Zhang
A novel adaptive refinement method is proposed for ordinary state-based PD in both 2D and 3D simulation. We perform an efficient adaptive refinement by splitting a parent particle into several child particles directly which does not require any information from its adjacent particles. The state variables and neighbor list of the child particles can be obtained directly from their parent particles, so that the proposed adaptive refinement method is very efficient and can be easily applied in both 2D and 3D cases without complex adjacent particle list building. To maintain simulation accuracy, the deformation of the parent particles need to be taken into consideration in the refinement process. Therefore, the accumulated strain of a PD particle is defined to calculate the position vector of its child particles in current configuration. With the consideration of the deformation of the parent particle, the fake damage is avoid. And a new criterion is established based on the particle accumulated strain to efficiently determine the particles need to be refined. The numerical examples studied show that the proposed adaptive refinement correctly captures the complicated crack propagation process and eliminates the fake crack growth with slight extra simulation cost, compared to the coarse discretization.
{"title":"A novel peridynamics refinement method with dual-horizon peridynamics","authors":"Zhixin Zeng, Xiong Zhang","doi":"10.1007/s00366-024-01983-3","DOIUrl":"https://doi.org/10.1007/s00366-024-01983-3","url":null,"abstract":"<p>A novel adaptive refinement method is proposed for ordinary state-based PD in both 2D and 3D simulation. We perform an efficient adaptive refinement by splitting a parent particle into several child particles directly which does not require any information from its adjacent particles. The state variables and neighbor list of the child particles can be obtained directly from their parent particles, so that the proposed adaptive refinement method is very efficient and can be easily applied in both 2D and 3D cases without complex adjacent particle list building. To maintain simulation accuracy, the deformation of the parent particles need to be taken into consideration in the refinement process. Therefore, the accumulated strain of a PD particle is defined to calculate the position vector of its child particles in current configuration. With the consideration of the deformation of the parent particle, the fake damage is avoid. And a new criterion is established based on the particle accumulated strain to efficiently determine the particles need to be refined. The numerical examples studied show that the proposed adaptive refinement correctly captures the complicated crack propagation process and eliminates the fake crack growth with slight extra simulation cost, compared to the coarse discretization.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"16 1","pages":""},"PeriodicalIF":8.7,"publicationDate":"2024-04-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140810872","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}