Engineering with Computers最新文献

{"title":"Parameterized shape optimization of a bi-leaflet heart valved conduit for pediatric applications.","authors":"Chuan Luo, Kewei Li, Abigail R Herschman, Mingze Sun, Haim Waisman, Vijay Vedula, Jeffrey W Kysar, David Kalfa","doi":"10.1007/s00366-026-02311-7","DOIUrl":"10.1007/s00366-026-02311-7","url":null,"abstract":"<p><p>Congenital heart defects could affect the right ventricular outflow tract in pediatric patients. As a result, pediatric pulmonary valve replacements are often needed to effectively mediate unidirectional blood flow from the right ventricle to the pulmonary artery. The present work aims to optimize a parameterized shape of a bi-leaflet heart-valved conduit using multi-objective optimization algorithms. We developed an integrated framework that automatically facilitates the design, structural mechanics simulation, and design optimization of the prosthetic valve. Bezier curves are employed to represent the geometric profile of both the free edge and the attachment edge of the leaflets, while a genetic algorithm updates the parameterized design variables during optimization. A quasi-static finite element analysis (FEA) model simulates valve opening and closure mechanics under a prescribed hemodynamic pressure profile. The objective is to find the optimal leaflet shape, enhancing durability and functionality by minimizing the leaflets' maximum principal stress and the orifice area at valve closure. An optimized design is selected from the final Pareto fronts for prototyping. Numerical results demonstrate an improved performance of the optimized valve over the initial design, indicating the complex impact of the valve geometry on valve performance metrics and underscoring the imperative for design optimization. The performance of optimal valve design is further analyzed using fluid-structure interaction (FSI) modeling to evaluate its performance under dynamic loading conditions.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":""},"PeriodicalIF":4.9,"publicationDate":"2026-04-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13059700/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147644578","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Coupling finite volume-lattice Boltzmann methods for advanced heat transfer simulations.","authors":"Yang Zhou, Alessandro De Rosis, Alistair Revell","doi":"10.1007/s00366-026-02288-3","DOIUrl":"https://doi.org/10.1007/s00366-026-02288-3","url":null,"abstract":"<p><p>We present a high-performance coupled framework that advances the integration of the finite volume method (FVM) and the lattice Boltzmann method (LBM) for multi-physics thermal flow simulations, including heat conduction, conjugated heat transfer, natural and forced convection, and phase change. The proposed scheme employs a central-moments-based collision operator for both velocity and temperature fields, substantially improving numerical stability and accuracy over traditional approaches within the LBM community. The reconstruction strategy, combining regularised and high-order truncated equilibrium methods, ensures smooth and accurate data exchange at FVM-LBM coupling interfaces. The implementation employs the Parallel Location and Exchange coupling library, enabling efficient and scalable communication between the FVM and LBM. Validation against standard benchmark problems and complex melting scenarios demonstrates excellent numerical accuracy and convergence. These algorithmic advances establish the proposed framework as a significant step forward in coupled FVM-LBM methods for multiscale thermal flow problems.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"56"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12935858/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147325045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Near real-time adaptive isotropic and anisotropic image-to-mesh conversion for cerebral aneurysm simulations.","authors":"Kevin Garner, Chander Sadasivan, Nikos Chrisochoides","doi":"10.1007/s00366-026-02287-4","DOIUrl":"10.1007/s00366-026-02287-4","url":null,"abstract":"<p><p>This paper presents two performance optimization techniques for a mesh adaptation method that is designed to help streamline the discretization of complex vascular geometries within the numerical modeling process. This method is integrated into a pipeline with an image-to-mesh conversion tool to generate adaptive anisotropic meshes from segmented medical images. The pipeline is shown to satisfy quality, fidelity, smoothness, and robustness requirements while providing near real-time performance for medical image-to-mesh conversion. Tested with two brain aneurysm cases and utilizing up to 96 CPU cores within a single, multicore node on Purdue University's Anvil supercomputer, the parallel adaptive anisotropic meshing method utilizes a hierarchical load balancing model (designed for large, cc-NUMA shared memory architectures) and contains an optimized local reconnection operation that performs three times faster than its original implementation from previous studies. While utilizing a new user-defined sizing function, we also show an adaptive isotropic method that generates meshes with good quality and fidelity of up to approximately 50 million elements in less than a minute while the adaptive anisotropic method is shown to generate approximately the same number of elements in about 5 min.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"45"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12913359/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146225822","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Implicit sub-stepping scheme for critical state soil models.","authors":"Hoang-Giang Bui, Jelena Niníc, Günther Meschke","doi":"10.1007/s00366-026-02309-1","DOIUrl":"https://doi.org/10.1007/s00366-026-02309-1","url":null,"abstract":"<p><p>The stress integration of critical soil model is usually based on implicit Euler algorithm, where the stress predictor is corrected by employing a return mapping algorithm. In the case of a large load step, the solution of the local nonlinear system to compute the plastic multiplier may not be attained. To overcome this problem, a sub-stepping scheme is generally used to improve the convergence of the local nonlinear system solution strategy. Nevertheless, the complexity of the tangent operator of the sub-stepping scheme is high. This complicates the use of Newton-Raphson algorithm to obtain global quadratic convergence. In this paper, a formulation for consistent tangent operator is developed for implicit sub-stepping integration for the modified Cam-Clay model and unified Clay and Sand model. This formulation is highly efficient and can be used with problem involving arbitrary large load step, such as tunnel simulation.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"76"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13038472/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147608358","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A differential equation-driven update strategy for density-based topology optimization: implementation with MATLAB codes.","authors":"Yang Liu, Wei Tan","doi":"10.1007/s00366-025-02237-6","DOIUrl":"https://doi.org/10.1007/s00366-025-02237-6","url":null,"abstract":"<p><p>Differential equation-driven evolution strategies are often associated with boundary-driven topology optimization methods, such as the level set method. However, differential equations can also be utilized effectively in density-based approaches. This paper presents a design update scheme formulated using differential equations to evolve elemental densities in topology optimization. The proposed scheme transforms the differential equation into an absolute increment format, closely related to the optimality criteria (OC) method, which is traditionally implemented in a relative increment format in density-based methods. The relative increment format of the OC method typically ensures an efficient and stable optimization process, whereas the absolute increment format tends to enable a more active and responsive optimization process, potentially leading to optimized results with improved performance. Furthermore, the absolute increment format can be converted into a relative one if needed. This study explores compliance minimization problems for both isotropic composite and single-material cases. Detailed MATLAB implementations for these cases are presented and thoroughly explained. Numerical examples demonstrate that the differential equation-driven update scheme effectively addresses density distribution optimization problems, offering an alternative to classical density methods.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 1","pages":"34"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12862039/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146112616","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Isogeometric suitable coupling methods for partitioned multiphysics simulation with application to fluid-structure interaction.","authors":"Jing-Ya Li, Hugo M Verhelst, Henk den Besten, Matthias Möller","doi":"10.1007/s00366-026-02299-0","DOIUrl":"10.1007/s00366-026-02299-0","url":null,"abstract":"<p><p>This paper presents spline-based coupling methods for partitioned multiphysics simulations, specifically designed for isogeometric analysis (IGA) based solvers. Traditional vertex-based coupling approaches face significant challenges when applied to IGA solvers, including geometric accuracy issues, interpolation errors, and substantial communication overhead. The methodology draws on the IGA mathematical framework to deliver coupling solutions that preserve the high-order continuity and exact geometric representation of splines. We develop two complementary strategies: (1) a spline-vertex coupling method that enables efficient interaction between IGA and conventional solvers, and (2) a fully isogeometric coupling approach that maximizes accuracy for IGA-to-IGA communication. Both theoretical analysis and extensive numerical experiments demonstrate that our spline-based methods significantly reduce communication overhead compared to traditional approaches while simultaneously enhancing geometric accuracy through exact boundary representation and maintaining higher-order solution continuity across the coupled interfaces. We quantitatively confirm the communication efficiency benefits through systematic measurements of both transfer times and data volumes across various mesh refinement levels, with experimental results closely aligning with our theoretical predictions. Our benchmark studies further demonstrate the geometric fidelity advantages through exact boundary representation, while also highlighting how the inherent mathematical structure of splines naturally preserves solution derivatives across interfaces without requiring additional computation or specialized transfer algorithms. This work not only provides efficient coupling strategies tailored to IGA-based solvers but also establishes a practical bridge between IGA and traditional discretization methods in partitioned multiphysics simulations. By offering viable options for coupling conventional solvers with IGA-based components, our approach enables broader adoption of IGA in established simulation workflows while ensuring accurate and high-performance interface communications.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"69"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13005864/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147503431","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A computational framework to predict the spreading of Alzheimer's disease.","authors":"Ana Vazquez-Palomo, Covadonga Betegón, Johannes Weickenmeier, Emilio Martínez-Pañeda","doi":"10.1007/s00366-026-02313-5","DOIUrl":"10.1007/s00366-026-02313-5","url":null,"abstract":"<p><p>Alzheimer's disease is characterised by the spreading of misfolded proteins and progressive structural changes in the brain. Despite significant clinical research, understanding how microscopic protein dynamics translate into macroscopic tissue degeneration remains a major challenge. In this work, we present a three-dimensional, finite element-based computational framework to model disease progression by combining multi-protein transport and brain tissue deformation within anatomically realistic geometries. The propagation of toxic tau and amyloid-[Formula: see text] proteins is described using reaction-diffusion equations of the Fisher-Kolmogorov type, incorporating prion-like growth mechanisms and anisotropic transport along white matter fibre tracts. Brain atrophy is represented through a hyperelastic constitutive model driven by protein-dependent volume loss. Subject-specific simulations are achieved through an automated preprocessing pipeline that generates finite element meshes and reconstructs axonal orientation fields from medical imaging data. The model reproduces key morphological patterns observed in Alzheimer's disease and shows good quantitative agreement with longitudinal imaging measurements. Overall, the proposed framework offers an extensible computational platform for studying Alzheimer's disease progression across subject-specific brain geometries. The models developed, including the image processing framework (BrainImage2Mesh) and the coupled bio-chemo-mechanical COMSOL finite element implementation, are made freely available to download at https://mechmat.web.ox.ac.uk/codes.</p><p><strong>Supplementary information: </strong>The online version contains supplementary material available at 10.1007/s00366-026-02313-5.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"78"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC13043544/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147622106","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A differentiable variational model for structural self-contact and fracture.","authors":"Mirko Ciceri, Charlie Aveline, Dilaksan Thillaithevan, Robert Hewson, Matthew Santer","doi":"10.1007/s00366-026-02285-6","DOIUrl":"https://doi.org/10.1007/s00366-026-02285-6","url":null,"abstract":"<p><p>Numerical modelling of structural self-contact and crack propagation presents significant challenges due to the inherently discontinuous and non-differentiable nature of the underlying physical phenomena. Traditional contact models demand explicit definition and tracking of contact points, while fracture models often rely on predefined crack initiation sites, sharp interfaces, and re-meshing. This study introduces a novel framework that overcomes these limitations within a unified and numerically stable variational formulation. The contact phenomenon is described through the hyperelastic third medium contact model and fracture is represented by a phase field. Structures are embedded in a third medium that stiffens under compression, enabling the transfer of forces between structural members. Crack propagation occurs in regions in which it is energetically favourable for the system to evolve toward a fully damaged state, specifically where the critical energy release rate is exceeded. Careful treatment is required when coupling the two phenomena, particularly concerning the void material behaviour. This work presents an efficient and differentiable numerical model that captures both nonlinear phenomena within a unified framework. This framework will allow designers and engineers to efficiently analyse complex nonlinear structural behaviours, previously requiring separate models that involved pre-defined crack initiation sites and contact points. Lastly, the differentiable nature of the model facilitates straightforward future integration into topology optimisation pipelines, providing designers the ability to intentionally design for and leverage self-contact interactions and material failure as functional, performance-enhancing features.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 2","pages":"52"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12920418/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147270075","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Weak wall boundary conditions for compressible flows.","authors":"Monu Jaiswal, Manoj R Rajanna, Md Rhyhanul Islam, Ming-Chen Hsu, Yuri Bazilevs","doi":"10.1007/s00366-025-02232-x","DOIUrl":"10.1007/s00366-025-02232-x","url":null,"abstract":"<p><p>Weak imposition of essential boundary conditions (i.e., weak BCs) for the Navier-Stokes equations of incompressible flows allows a certain amount of controlled numerical flow slip on the solid surface. Numerical flow slip mimics the presence of a thin boundary layer that would otherwise need to be captured using a fine mesh resolution. As a result, weak BCs enable the use of coarser meshes near solid walls without sacrificing numerical solution accuracy, which significantly reduces the computational costs, especially for 3D, wall-bounded turbulent flows. However, weak BCs for compressible flows are not as well understood as those for the incompressible-flow case. In particular, numerical instabilities were observed in some cases where the weak BCs were simultaneously imposed for the velocity and temperature fields. In the present effort, to address these stability issues, we develop a methodology for the design of compressible-flow weak BC operators and demonstrate the improved performance of the resulting weak BC formulations using challenging 2D and 3D test cases.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 1","pages":"16"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12804276/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145997632","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Bridging experiments and defects' mechanics: a data-driven toolbox for configurational force analysis.","authors":"Abdalrhaman Koko, Alya Abdelnour, Thorsten H Becker, T James Marrow","doi":"10.1007/s00366-025-02262-5","DOIUrl":"10.1007/s00366-025-02262-5","url":null,"abstract":"<p><p>Understanding the mechanical behaviour of defective materials is key to predicting failure and enhancing performance. Traditional fracture mechanics often requires assumptions about geometry and loading that are unavailable in experimental systems. We present a MATLAB-based computational toolbox that extracts configurational forces and mixed-mode SIFs directly from experimentally measured displacement or deformation gradient fields, like digital image/volume correlation and high (angular) resolution electron backscatter diffraction. The toolbox implements path-independent energy integrals, including the <i>J</i>- and <i>M</i>-integrals, and introduces a novel mode decomposition formulation that isolates mode I-III SIFs contributions without predefined specimen geometries, applied loads, or boundary conditions. Applications to microcracks, dislocations, and fatigue cracks demonstrate its robust, geometry-independent characterisation, which can enable data-driven analysis of defect behaviour in anisotropic and complex materials. The framework is material-agnostic in principle and operates directly on experimental fields; however, its current implementation assumes small-strain kinematics, making it most applicable to linear and anisotropic elastic and elastoplastic materials such as metals and ceramics.</p><p><strong>Supplementary information: </strong>The online version contains supplementary material available at 10.1007/s00366-025-02262-5.</p>","PeriodicalId":11696,"journal":{"name":"Engineering with Computers","volume":"42 1","pages":"21"},"PeriodicalIF":4.9,"publicationDate":"2026-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12804342/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145997599","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}