An Interface-enriched Generalized Finite Element Method (IGFEM) is proposed for the coupled hydro-mechanical analysis of deformable, saturated porous media consisting of distinct, perfectly bonded material phases. The spatial discretization of the momentum balance equation and storage equation is derived using IGFEM, followed by the time discretization of these equations via the generalized Newmark method. This leads to a fully coupled system of nonlinear equations, which is solved iteratively using a monolithic update scheme. The IGFEM formulation is proficient in accurately capturing weak discontinuities in both the solid phase displacement field and the fluid phase pressure field at material interfaces, by placing enriched nodes directly on these interfaces. Several numerical examples demonstrate that the proposed IGFEM formulation not only matches the accuracy of standard FEM with a conformal mesh, but also outperforms the eXtended/Generalized Finite Element Method (XFEM/GFEM). Moreover, it can accurately capture complex, non-planar interfaces without requiring mesh alignment, highlighting the method's flexibility and robustness for practical hydro-mechanical analyses of porous media with geometrically intricate boundaries. Overall, IGFEM provides a highly accurate and efficient approach for solving transient coupled problems involving material interfaces.

We uncover a chain of nonlinear modal interactions in softly clamped nanostring resonators. The process involves the sequential coupling of five mechanical modes, during frequency sweeps, yielding a broad nonlinear response with nearly constant amplitude. We demonstrate that soft clamping enables this cascaded energy transfer and amplifies the effective geometric nonlinearity of the driven mode by an order of magnitude. Analytical and finite element-based reduced-order models capture the key features of the coupling cascade and clarify its underlying mechanism. The phenomenon is generic in nonlinear vibrational systems and can be tailored through soft-clamping design strategies.

This paper addresses the topology optimization of thermocouples for cooling applications, considering stress constraints to enhance reliability under service loads. We provide a first approach to derive sensitivities using SIMP (solid isotropic material with penalization) for thermo-electro-mechanical systems with temperature-dependent material properties. The proposed formulation decouples the thermoelectrical system from the mechanical degrees of freedom reducing computational memory usage from a fully coupled approach. The study focuses on the formulation of thermocouples for cooling applications using the Peltier effect, which considers electrical power limits, electrical working points, and material stress thresholds. Furthermore, while the thermoelectrical problem does not show the need for filtering techniques, including the mechanical degrees of freedom, we show that we recover undesirable porous optimized designs. We provide 2D thermocouple example optimizations with geometries and boundary conditions based on a practical case for the implementation of thermoelectric coolers in the Minimum Ionizing Particle Timing Detector (MTD) at CERN. The optimizations are performed with increased complexity, including the unfiltered thermoelectrical and thermo-electro-mechanical problems and a Helmholtz-filtered examples. The optimizations are compared with constant and nonlinear material properties with temperature and with respect to the consideration of air-conductance losses within the devices. Although more efficient topologies can be achieved without the need for volume constraints, we include an example with a constraint of 60% volume to understand its effect on the design and provide a methodology to reduce semiconductor-associated costs at lower efficiency costs. Finally, we explore the same formulation in 3D. The results provide guidelines for manufacturing compliant thermocouples, increasing their reliability without decreasing efficiency.

Continuous fiber fused filament fabrication (CF4) is a layer-by-layer additive manufacturing technique that deposits continuous fiber fused filaments (CFFFs) with a significant in-plane variation of the fiber trajectory, thereby offering great flexibility in fabricating variable-stiffness composite laminates (VSCLs). We introduce a topology optimization method for the design of additively manufactured VSCLs made of overlapping, fiber-reinforced bars. The proposed method is based on geometry projection (GP) techniques, whereby the bars are represented by high-level geometric primitives. As in other GP techniques, this high-level parameterization is mapped onto a fixed structured finite element mesh for conducting analysis, as in density-based topology optimization techniques. However, unlike previous GP techniques that have demonstrated their applicability in designing structures as assemblies of individual fiber-reinforced components, this work focuses on the design of composite structures that adhere to CF4 manufacturing processes. Therefore, we first formulate a material interpolation scheme that better captures the stiffness at the composite's joints obtained from bar overlaps as a stack. Second, the proposed material interpolation employs composite laminate theory to capture the in-plane and out-of-plane behavior of the structure. Third, to produce designs that conform to the CF4 process, we also proposed a novel length constraint formulation in the form of penalization on the projection scheme, which ensures a minimum length for all the bars. This minimum length limit does not require adding a constraint to the optimization problem. The efficacy and efficiency of the proposed method are demonstrated by a series of compliance minimization problems with in-plane and/or out-of-plane loading. The methodology is also applied to the design of a displacement inverter compliant mechanism.

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.

Nonlinear dynamic simulations of mechanical resonators have been facilitated by the advent of computational techniques that generate nonlinear reduced order models (ROMs) using the finite element (FE) method. However, designing devices with specific nonlinear characteristics remains inefficient since it requires manual adjustment of the design parameters and can result in suboptimal designs. Here, we integrate an FE-based nonlinear ROM technique with a derivative-free optimization algorithm to enable the design of nonlinear mechanical resonators. The resulting methodology is used to optimize the support design of high-stress nanomechanical Si ₃N ₄ string resonators, in the presence of conflicting objectives such as simultaneous enhancement of Q-factor and nonlinear Duffing constant. To that end, we generate Pareto frontiers that highlight the trade-offs between optimization objectives and validate the results both numerically and experimentally. To further demonstrate the capability of multi-objective optimization for practical design challenges, we simultaneously optimize the design of nanoresonators for three key figure-of-merits in resonant sensing: power consumption, sensitivity and response time. The presented methodology can facilitate and accelerate designing (nano) mechanical resonators with optimized performance for a wide variety of applications. (Figure presented.)

Topology optimization (TO) has found applications across a wide range of disciplines but remains underutilized in practice. Key barriers to broader adoption include the absence of versatile commercial software, the need for in-depth knowledge of the methodology from the user, and high computational demands. Additionally, challenges such as ensuring manufacturability, tuning hyper-parameters, and integrating subjective design elements like esthetics further hinder its widespread use. Emerging technologies like augmented reality and virtual reality offer transformative potential for TO. By enabling intuitive, gesture-based human–computer interactions, these immersive environments bridge the gap between human intuition and computational processes. They provide the means to integrate subjective human judgment into optimization workflows in real time, creating a paradigm shift toward interactive and immersive design. Here, we introduce the concept of immersive topology optimization (ITO) as a novel design paradigm that leverages augmented reality for TO. By incorporating real-time human interaction during the optimization process, and the subsequent visualization of the design in its intended target location, ITO has the potential to reduce lead times, enhance manufacturability, and improve design integration. To demonstrate this ITO design paradigm, we present ARCADE: Augmented Reality Computational Analysis and Design Environment. Developed in Swift for the Apple Vision Pro mixed reality headset, ARCADE enables users to define, manipulate, and solve structural compliance minimization problems within an augmented reality setting. We provide the source code in Swift for the optimization procedure. While initially demonstrated for minimizing compliance, our framework could also be extended to other disciplines, paving the way for a new era of interactive and immersive computational design.

Clamp-on ultrasonic flowmeters suffer from crosstalk—i.e., measurement errors due to the interference of signals generated in solid regions and solid–fluid interfaces with the required signal from the fluid. Although several approaches have been proposed to alleviate crosstalk, they only work in specific ranges of flow rates and pipe diameters, and some also introduce additional issues. We propose a novel clamp-on system design where the transmitting and receiving wedges are embedded with directional noise filtering mechanisms based on phononic crystals (PnCs) possessing directional band gaps (DBGs). PnCs are artificial materials consisting of periodic structures arrayed in a matrix medium exhibiting band gaps – i.e., frequency ranges where waves are attenuated – due to Bragg scattering. DBGs enable PnCs to propagate waves in specific directions while suppressing them in other directions. By guiding the input signal through the transmitting wedge to the wall, we minimize the generation of noise signals due to secondary reflections within the wedge. Similarly, by using the directionality of the DBG PnC in the receiver, we limit the effects of noise signals (that arrive in different directions) in the receiver. We numerically verify the DBG PnC embedded wedges’ performance by comparing wave propagation aspects of the PnC embedded clamp-on system with a standard clamp-on device. To that end, we develop accurate wave propagation models based on the Discontinuous Galerkin finite element method. By incorporating DBG PnCs into the wedges, we obtain about 20 dB increase in the signal-to-noise ratio compared to the clamp-on system with standard wedges.

Neural networks (NNs) hold great promise for advancing inverse design via topology optimization (TO), yet misconceptions about their application persist. This article focuses on neural topology optimization (neural TO), which leverages NNs to reparameterize the decision space and reshape the optimization landscape. While the method is still in its infancy, our analysis tools reveal critical insights into the NNs’ impact on the optimization process. We demonstrate that the choice of NN architecture significantly influences the objective landscape and the optimizer’s path to an optimum. Notably, NNs introduce non-convexities even in otherwise convex landscapes, potentially delaying convergence in convex problems but enhancing exploration for non-convex problems. This analysis lays the groundwork for future advancements by highlighting: (1) the potential of neural TO for non-convex problems and dedicated GPU hardware (the “good”), (2) the limitations in smooth landscapes (the “bad”), and (3) the complex challenge of selecting optimal NN architectures and hyperparameters for superior performance (the “ugly”).

Enriched finite element methods (e-FEMs) have become a popular choice for modeling problems containing material discontinuities (e.g., multi-phase materials and fracture). The main advantage as compared to the standard finite element method (FEM) remains the versatility in the choice of discretizations, since e-FEMs resolve discontinuities by completely decoupling them from the finite element mesh. However, modeling complex kinematics such as branching and merging of discrete cracks remains challenging. This article extends previous research on the Discontinuity-Enriched Finite Element Method (DE-FEM) for simulating quasi-static crack propagation in brittle materials. In DE-FEM enrichments are added to nodes created directly along discontinuities. Most notably, we demonstrate DE-FEM can resolve complex kinematics, namely the modeling of multiple cracks propagating and merging—and with a straightforward computer implementation. We validate the formulation with experimental results carried out on a compact tension specimen. Other numerical examples show the capability of DE-FEM in capturing crack paths similar to those observed in the literature.

The quality factor (Q factor) of nanomechanical resonators is influenced by geometry and stress, a phenomenon called dissipation dilution. Studies have explored maximizing this effect, leading to softly-clamped resonator designs. This paper proposes a topology optimization methodology to design two-dimensional nanomechanical resonators with high Q factors by maximizing dissipation dilution. A formulation based on the ratio of geometrically nonlinear to linear modal stiffnesses of a prestressed finite element model is used, with its corresponding adjoint sensitivity analysis formulation. Systematic design in square domains yields geometries with comparable Q factors to literature. We analyze the trade-offs between resonance frequency and quality factor, and how these are reflected in the geometry of resonators. We further apply the methodology to optimize a resonator on a full hexagonal domain. By using the entire mesh—i.e., without assuming any symmetries—we find that the optimizer converges to a two-axis symmetric design comprised of four tethers.

Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form p_k=p_σ>p_α=1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

The computational analysis of nanophotonic devices is usually carried out via the standard finite element method (FEM). However, FEM requires meshes that are fitted to the devices’ boundaries, so making changes to the geometry (and thus the mesh) results in an inefficient process at best. Such an approach is therefore at odds when conducting design, which requires the analysis of multiple device geometries until reaching a satisfactory solution. Computational design tools such as topology optimization are often used, but the use of density-based representations of geometry inevitably leads to other issues—e.g., pixelized fuzzy boundaries with “gray material” (that does not correspond to dielectric nor vacuum) have an adverse effect on the devices’ interaction with electromagnetic waves. In this paper we propose an interface-enriched generalized finite element method (IGFEM) for the analysis of two-dimensional electromagnetic scattering and eigenvalue problems. IGFEM enables the use of finite element meshes that are completely decoupled from the problem's geometry. The analysis procedure is further coupled to a level set description of topology, resulting in a versatile enriched approach to topology optimization; this level set-based interface-enriched topology optimization procedure is devoid of the issues mentioned above regarding density-based methods, and yields crisp “black-and-white” designs that are devoid of jagged fuzzy edges. We first demonstrate that the analysis procedure achieves the same convergence rate as that of standard FEM using geometry-fitted meshes. We then compare the convergence properties of IGFEM with Nitsche's method on a problem containing an embedded straight interface. Finally, we conduct topology optimization for designing both a 2-D metalens and a 2-D reflector, maximizing their ability to focus light onto a target point.

Ultrasonic flowmeters face unique challenges since, in addition to withstanding high fluid pressures, they have to avoid crosstalk, which is the interaction of the signals traveling through the fluid and the solid pipe. To avoid the crosstalk, which leads to poor accuracy or complete loss of the required signal, we develop a mounting mechanism based on phononic crystals (PnCs), which are artificial periodic materials possessing band gaps (BGs) due to Bragg scattering. These PnC structures should also possess high mechanical strength to sustain the fluid pressure. Designing PnCs for such applications is challenging as the BG width and the resistance to mechanical loading are conflicting objectives. To circumvent this, we propose a step-by-step design procedure to optimize both mechanical strength and wave attenuation performance of a single-phase 3D PnC waveguide using parametric sweeping and sensitivity analysis. We use finite element analysis (FEA) to characterize the behavior of the periodic unit cell and the waveguide. Since accurate dynamic FEA at high frequencies is computationally demanding, we develop surrogate models at different levels of the design process. We also consider additive manufacturing aspects in the design procedure, which we validate by 3D-printing the final design and measuring the parameters via computer tomography.

Current design methods for flexure (or compliant) mechanisms regard stress as a secondary, limiting factor. This is remarkable because stress is also known as a useful design parameter. In this paper we propose the Stress And Geometry (STAGE) method, to design the geometry of a flexure mechanism together with a desired stress field. From this design, the stress-free to-be-fabricated geometry is computed using the inverse finite element method. To demonstrate the potential of the method, the geometry of the well-known crossed-flexure pivot is taken as example. We first show how this mechanism can be redesigned for the same functional geometry with various internal stresses. This results for a specific choice of stress field in a design of a crossed-flexure pivot with 23% lower peak stresses during motion as compared to the known designs, for a ±45° rotation. We then present a second example, of a Folded Leaf Spring (FLS). With a parameter sweep the optimal stress field is calculated, showing a peak stress reduction of 28% during motion. This result was validated with an experiment, showing a normalized mean absolute error of 5.5% between experiment and theory. With a second experiment it was verified that the functional geometry of the FLS with internal stresses was equal to the one without internal stresses, with geometric deviations smaller than half the thickness of the flexures.

Underwater noise resulting from the monopile driving process can cause severe damage to marine wildlife, such as hearing injury, behavioral disturbance, or even death. Although current noise-attenuation techniques used in this process have shown a significant noise reduction at high frequency ranges, mitigating low-frequency noise is still extremely challenging. To address the problem, here we propose an elastic metamaterial-based structure composed of single-phase resonant structures. The proposed structure, which we call a meta-interface, is introduced between the monopile and the hammer and is used to remove energy from the input signal associated with high noise levels. To that end, we first identify the frequency ranges associated with high sound pressure levels, which were shown to be related to the monopile's eigenmodes. Then we design the meta-interface's periodic unit cells so that the elastic/acoustic waves at identified frequency ranges are attenuated. A meta-interface is then realized by replicating the unit cell along the monopile wall (matching the thickness) to form a ring-shaped layer, and then by stacking up these concentric layers. A frequency analysis of the pile driving system with the meta-interface shows that the new noise levels attain a significant attenuation in frequency ranges lower than 1000Hz. This demonstrates a novel solution for the low-frequency underwater noise issue during the hammering of offshore monopiles.

Inorganic scintillators often use exotic, expensive materials to increase their light yield. Although material chemistry is a valid way to increase the light collection, these methods are expensive and limited to the material properties. As such, alternative methods such as the use of specific reflective coatings and crystal optical shapes are critical for the scintillator crystal design procedure. In this paper, we explore the modeling of a scintillator and silicon-photomultiplier (SiPM) assembly detector using GEANT4. GEANT4, an open-source software for particle–matter interaction based on ray-tracing, allows the modeling of a scintillator-based detector while offering methods to simplify and study the computational requirements for a precise calculation of the light collection. These studies incorporate two different geometries compatible with the barrel timing layer (BTL) particle detector that is being built for the compact muon solenoid (CME) experiment at CERN. Furthermore, the geometry of our model is parameterized using splines for smoother results and meshed using GMSH to perform genetic numerical optimization of the crystal shape through genetic algorithms, in particular non-dominated sorting genetic algorithm II (NGSAII). Using NSGA-II, we provide a series of optimized scintillator geometries and study the trade-offs of multiple possible objective functions including the light output, light collection, light collection per energy deposited, and track path length. The converged Pareto results according to the hypervolume indicator are compared to the original simplified design, and a recommendation towards the use of the light collection per energy deposition and track path length is given based on the results. The results provide increases in this objective of up to 18% for a constant volume for a geometry compatible with the current design of the BTL detector.

Although strain engineering and soft-clamping techniques for attaining high Q-factors in nanoresonators have received much attention, their impact on nonlinear dynamics is not fully understood. In this study, we show that nonlinearity of high-Q Si₃N₄ nanomechanical string resonators can be substantially tuned by support design. Through careful engineering of support geometries, we control both stress and mechanical nonlinearities, effectively tuning nonlinear stiffness of two orders of magnitude. Our approach also allows control over the sign of the Duffing constant resulting in nonlinear softening of the mechanical mode that conventionally exhibits hardening behavior. We elucidate the influence of support design on the magnitude and trend of the nonlinearity using both analytical and finite element-based reduced-order models that validate our experimental findings. Our work provides evidence of the role of soft-clamping on the nonlinear dynamic response of nanoresonators, offering an alternative pathway for nullifying or enhancing nonlinearity in a reproducible and passive manner.

Fundamentals of Enriched Finite Element Methods provides an overview of the different enriched finite element methods, detailed instruction on their use, and their real-world applications, recommending in what situations they are best implemented. It starts with a concise background on the theory required to understand the underlying principles behind the methods before outlining detailed instruction on implementation of the techniques in standard displacement-based finite element codes. The strengths and weaknesses of each are discussed, as are computer implementation details, including a standalone generalized finite element package, written in Python. The applications of the methods to a range of scenarios, including multiphase, fracture, multiscale, and immersed boundary (fictitious domain) problems are covered, and readers can find ready-to-use code, simulation videos, and other useful resources on the companion website of the book.