Kangan Li , Antonio Rodríguez-Ferran , Guglielmo Scovazzi
{"title":"Crack branching and merging simulations with the shifted fracture method","authors":"Kangan Li , Antonio Rodríguez-Ferran , Guglielmo Scovazzi","doi":"10.1016/j.cma.2024.117528","DOIUrl":"10.1016/j.cma.2024.117528","url":null,"abstract":"<div><div>We propose a relatively simple and mesh-independent approach to model crack branching and merging using the Shifted Fracture Method (SFM), within the class of Shifted Boundary Methods. The proposed method achieves mesh independency by accurately accounting for the area of the fracture surface, in contrast to traditional element-deletion/node-release techniques. In the SFM, the <em>true</em> fracture is embedded into the computational grid, but the fracture interface conditions are modified (shifted) by means of Taylor expansions to the <em>surrogate</em> fracture composed of full edges/faces in two/three dimensions. This avoids numerical integration on cut elements, so that the data structures and geometrical treatment of cut elements are simple, while mesh-independent results and accurate fracture approximations are still maintained. We demonstrate the capabilities of the proposed approach in a number of prototypical numerical experiments.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117528"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637437","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Direct data-driven algorithms for multiscale mechanics","authors":"E. Prume , C. Gierden , M. Ortiz , S. Reese","doi":"10.1016/j.cma.2024.117525","DOIUrl":"10.1016/j.cma.2024.117525","url":null,"abstract":"<div><div>We propose a randomized data-driven solver for multiscale mechanics problems which improves accuracy by escaping local minima and reducing dependency on metric parameters, while requiring minimal changes relative to non-randomized solvers. We additionally develop an adaptive data-generation scheme to enrich data sets in an effective manner. This enrichment is achieved by utilizing material tangent information and an error-weighted k-means clustering algorithm. The proposed algorithms are assessed by means of three-dimensional test cases with data from a representative volume element model.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117525"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637439","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rutger A. Biezemans , Claude Le Bris , Frédéric Legoll , Alexei Lozinski
{"title":"MsFEM for advection-dominated problems in heterogeneous media: Stabilization via nonconforming variants","authors":"Rutger A. Biezemans , Claude Le Bris , Frédéric Legoll , Alexei Lozinski","doi":"10.1016/j.cma.2024.117496","DOIUrl":"10.1016/j.cma.2024.117496","url":null,"abstract":"<div><div>We study the numerical approximation of advection–diffusion equations with highly oscillatory coefficients and possibly dominant advection terms by means of the Multiscale Finite Element Method (MsFEM). The latter method is a now classical, finite element type method that performs a Galerkin approximation on a problem-dependent basis set, itself precomputed in an offline stage. The approach is implemented here using basis functions that locally resolve both the diffusion and the advection terms. Variants with additional bubble functions and possibly weak inter-element continuity are proposed. Some theoretical arguments and a comprehensive set of numerical experiments allow to investigate and compare the stability and the accuracy of the approaches. The best approach constructed is shown to be adequate for both the diffusion- and advection-dominated regimes, and does not rely on an auxiliary stabilization parameter that would have to be properly adjusted.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117496"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637441","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jixiang Zhang , Zhenzhong Chen , Ge Chen , Xiaoke Li , Pengcheng Zhao , Qianghua Pan
{"title":"A hyperspherical area integral method based on a quasi-Newton approximation for reliability analysis","authors":"Jixiang Zhang , Zhenzhong Chen , Ge Chen , Xiaoke Li , Pengcheng Zhao , Qianghua Pan","doi":"10.1016/j.cma.2024.117533","DOIUrl":"10.1016/j.cma.2024.117533","url":null,"abstract":"<div><div>The First-Order Reliability Method (FORM) is renowned for its high computational efficiency, but its accuracy declines when addressing the nNar Limit-State Function (LSF). The Second-Order Reliability Method (SORM) offers greater accuracy; however, its approximation formula can sometimes introduce errors. Additionally, SORM requires extra calculations involving the Hessian matrix, which can reduce its efficiency. To balance efficiency and accuracy, a Hyperspherical Area Integral Method based on a Quasi-Newton Approximation (HAI-QNAM) for reliability analysis is proposed. This method initially employs a quasi-Newton method to determine the Most Probable Point (MPP) of failure, calculate the reliability index, and obtain the approximate Hessian matrix. Then, based on the Curved Surface Integral (CSI) method, the area of the approximate failure domain and the area of the hypersphere are obtained. Using the proportionality of their areas, the failure probability is then calculated. Finally, the proposed method's accuracy and efficiency are validated through examples.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117533"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637436","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Isaac Bermúdez , Jessika Camaño , Ricardo Oyarzúa , Manuel Solano
{"title":"A conforming mixed finite element method for a coupled Navier–Stokes/transport system modeling reverse osmosis processes","authors":"Isaac Bermúdez , Jessika Camaño , Ricardo Oyarzúa , Manuel Solano","doi":"10.1016/j.cma.2024.117527","DOIUrl":"10.1016/j.cma.2024.117527","url":null,"abstract":"<div><div>We consider the coupled Navier–Stokes/transport equations with nonlinear transmission conditions, which constitute one of the most common models utilized to simulate a reverse osmosis effect in water desalination processes when considering feed and permeate channels coupled through a semi-permeate membrane. The variational formulation consists of a set of equations where the velocities, the concentrations, along with tensors and vector fields introduced as auxiliary unknowns and two Lagrange multipliers are the main unknowns of the system. The latter are introduced to deal with the trace of functions that do not have enough regularity to be restricted to the boundary. In addition, the pressures can be recovered afterwards by a postprocessing formula. As a consequence, we obtain a nonlinear Banach spaces-based mixed formulation, which has a perturbed saddle point structure. We analyze the continuous and discrete solvability of this problem by linearizing the perturbation and applying the classical Banach fixed point theorem along with the Banach–Nečas–Babuška result. Regarding the discrete scheme, feasible choices of finite element subspaces that can be used include Raviart–Thomas spaces for the auxiliary tensor and vector unknowns, piecewise polynomials for the velocities and the concentrations, and continuous polynomial space of lowest order for the traces, yielding stable discrete schemes. An optimal <em>a priori</em> error estimate is derived, and numerical results illustrating both, the performance of the scheme confirming the theoretical rates of convergence, and its applicability, are reported.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117527"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637438","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"A multigrid two-scale modeling approach for nonlinear multiphysical systems","authors":"Alaa Armiti-Juber, Tim Ricken","doi":"10.1016/j.cma.2024.117523","DOIUrl":"10.1016/j.cma.2024.117523","url":null,"abstract":"<div><div>High fidelity modeling of multiphysical systems is typically achieved using nonlinear coupled differential equations, often with multiscale model coefficients. These simulations are performed using finite-element methods with implicit time stepping. Within each time step, nonlinearities are numerically linearized using Newton-like iterative solvers, which increases the computational complexity. For multiscale systems fulfilling a scale-separation criterion, this complexity can be mitigated by upscaling the high fidelity models to describe the effective behavior of the system alone. To extend the validity of these models to systems with partial scale-separation, we propose a Multi-Grid Two-Scale (MGTS) modeling approach. This approach consists of the nonlinear upscaled models as a coarse-scale model and a linear fine-scale corrector. The derivation of the corrector is based on linearizing the high-fidelity models about their upscaled versions. This approach has the ability to capture most of the fine-scale perturbations, while maintaining a reduced computational complexity as a consequence of restricting the iterative solvers to coarse-scale grids in the upscaled domain. Additionally, the coarse- and fine-scale models are only weakly coupled, enabling a parallelization-in-time feature of the MGTS model. We apply the MGTS approach to a nonlinear model for fluid-saturated poro-elastic materials in thin domains. The performance of the MGTS approach is demonstrated by executing several numerical experiments based on the finite-element method.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117523"},"PeriodicalIF":6.9,"publicationDate":"2024-11-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637440","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
S. Rezaeiravesh , C. Gscheidle , A. Peplinski , J. Garcke , P. Schlatter
{"title":"In-situ estimation of time-averaging uncertainties in turbulent flow simulations","authors":"S. Rezaeiravesh , C. Gscheidle , A. Peplinski , J. Garcke , P. Schlatter","doi":"10.1016/j.cma.2024.117511","DOIUrl":"10.1016/j.cma.2024.117511","url":null,"abstract":"<div><div>The statistics obtained from turbulent flow simulations are generally uncertain due to finite time averaging. Most techniques available in the literature to accurately estimate these uncertainties typically only work in an offline mode, that is, they require access to all available samples of a time series at once. In addition to the impossibility of online monitoring of uncertainties during the course of simulations, such an offline approach can lead to input/output (I/O) deficiencies and large storage/memory requirements, which can be problematic for large-scale simulations of turbulent flows. Here, we designed, implemented and tested a framework for estimating time-averaging uncertainties in turbulence statistics in an in-situ (online/streaming/updating) manner. The proposed algorithm relies on a novel low-memory update formula for computing the sample-estimated autocorrelation functions (ACFs). Based on this, smooth modeled ACFs of turbulence quantities can be generated to accurately estimate the time-averaging uncertainties in the corresponding sample mean estimators. The resulting uncertainty estimates are highly robust, accurate, and quantitatively the same as those obtained by standard offline estimators. Moreover, the computational overhead added by the in-situ algorithm is found to be negligible allowing for online estimation of uncertainties for multiple points and quantities. The framework is general and can be used with any flow solver and also integrated into the simulations over conformal and complex meshes created by adopting adaptive mesh refinement techniques. The results of the study are encouraging for the further development of the in-situ framework for other uncertainty quantification and data-driven analyses relevant not only to large-scale turbulent flow simulations, but also to the simulation of other dynamical systems leading to time-varying quantities with autocorrelated samples.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117511"},"PeriodicalIF":6.9,"publicationDate":"2024-11-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637442","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Oscar Lopez-Pamies , John E. Dolbow , Gilles A. Francfort , Christopher J. Larsen
{"title":"Classical variational phase-field models cannot predict fracture nucleation","authors":"Oscar Lopez-Pamies , John E. Dolbow , Gilles A. Francfort , Christopher J. Larsen","doi":"10.1016/j.cma.2024.117520","DOIUrl":"10.1016/j.cma.2024.117520","url":null,"abstract":"<div><div>Notwithstanding the evidence against them, classical variational phase-field models continue to be used and pursued in an attempt to describe fracture nucleation in elastic brittle materials. In this context, the main objective of this paper is to provide a comprehensive review of the existing evidence against such a class of models as descriptors of fracture nucleation. To that end, a review is first given of the plethora of experimental observations of fracture nucleation in nominally elastic brittle materials under quasi-static loading conditions, as well as of classical variational phase-field models, without and with energy splits. These models are then confronted with the experimental observations. The conclusion is that they cannot possibly describe fracture nucleation in general. This because classical variational phase-field models cannot account for material strength as an independent macroscopic material property. The last part of the paper includes a brief summary of a class of phase-field models that can describe fracture nucleation. It also provides a discussion of how pervasively material strength has been overlooked in the analysis of fracture at large, as well as an outlook into the modeling of fracture nucleation beyond the basic setting of elastic brittle materials.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117520"},"PeriodicalIF":6.9,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660133","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dongyu Liu , Jian Zhang , Alejandro M. Aragón , Angelo Simone
{"title":"The discontinuity-enriched finite element method for multiple intersecting discontinuities","authors":"Dongyu Liu , Jian Zhang , Alejandro M. Aragón , Angelo Simone","doi":"10.1016/j.cma.2024.117432","DOIUrl":"10.1016/j.cma.2024.117432","url":null,"abstract":"<div><div>We extend the Discontinuity-Enriched Finite Element Method (DE-FEM) to simulate intersecting discontinuities, such as those encountered in polycrystalline materials, multi-material wedge problems, and branched cracks. The proposed hierarchical enrichment functions capture weak and strong discontinuities at junctions within a single formulation. Several numerical applications to branched cracks and polycrystalline microstructures under both thermal and mechanical loads are presented to demonstrate the proposed method. Results indicate that DE-FEM can accurately capture complex discontinuous primal and gradient fields and attain convergence rates comparable to those of standard FEM using fitted meshes. The main advantages of DE-FEM equipped with the proposed junction enrichment functions lie in the method’s ability to model intersecting discontinuities using meshes that are completely decoupled from them and its robustness in reproducing correct displacement and strain jumps across them, as demonstrated by a patch test. This work thus highlights the potential of DE-FEM for applications to problems characterized by the presence of multiple intersecting discontinuities, posing a valid alternative to traditional FEM and eXtended/Generalized Finite Element (X/GFEM) Methods.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117432"},"PeriodicalIF":6.9,"publicationDate":"2024-11-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142660132","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jiang Mo , Wang-Ji Yan , Ka-Veng Yuen , Michael Beer
{"title":"Efficient non-probabilistic parallel model updating based on analytical correlation propagation formula and derivative-aware deep neural network metamodel","authors":"Jiang Mo , Wang-Ji Yan , Ka-Veng Yuen , Michael Beer","doi":"10.1016/j.cma.2024.117490","DOIUrl":"10.1016/j.cma.2024.117490","url":null,"abstract":"<div><div>Non-probabilistic convex models are powerful tools for structural model updating with uncertain‑but-bounded parameters. However, existing non-probabilistic model updating (NPMU) methods often struggle with detecting parameter correlation due to limited prior information. Worth still, the unique core steps of NPMU, involving nested inner layer forward uncertainty propagation and outer layer inverse parameter updating, present challenges in efficiency and convergence. In response to these challenges, a novel and flexible NPMU scheme is introduced, integrating analytical correlation propagation and parallel interval bounds prediction to enable automatic detection of parameter correlations. In the forward uncertainty propagation phase, a linear coordinate transformation is applied to map the original parameter space to a standard hypercube space, simplifying correlation-involved bounds prediction into conventional interval bounds prediction. Moreover, an analytical correlation propagation formula is derived using a second-order response approximation to sidestep the complexities of geometry-based correlation calculations. To expedite forward propagation, a derivative-aware neural network model is employed to replace the physical solver, facilitating improved fitting capabilities and automatic differentiation, including the calculation of Jacobian and Hessian matrices essential for correlation propagation. The neural network's inherent parallelism accelerates interval bounds prediction through parallel computation of samples. In the inverse parameter updating phase, the block coordinate descent algorithm is embraced to narrow the search space and boost convergence capabilities, while the perturbation method is utilized to determine the optimal starting point for optimization. Two numerical examples illustrate the efficacy of the proposed method in updating structural models while considering correlations.</div></div>","PeriodicalId":55222,"journal":{"name":"Computer Methods in Applied Mechanics and Engineering","volume":"433 ","pages":"Article 117490"},"PeriodicalIF":6.9,"publicationDate":"2024-11-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142637364","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}