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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

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.

@en

Ultrasonic flowmeters that use transit-time ultrasonic transducers face measurement errors due to 'crosstalk,' whereby the working signal travels through the pipe wall and couplings, interfering with the signal from the fluid. Although various procedures have been proposed to solve the issue of crosstalk, they're limited to low-frequency ranges, or they are not effective in high-pressure environments. We propose a mounting mechanism based on a single-phase 3-D phononic crystal (PnC) waveguide that can mitigate crosstalk at high frequencies (megahertz range) and thus improve the flowmeters' measurement accuracy. PnCs are artificial materials consisting of periodically arranged scatterers thereby showing bandgaps (BGs) - ranges of frequencies where elastic/acoustic waves are attenuated - due to Bragg scattering. We design PnC wave filters by engineering the BG frequency range to the working signal of the ultrasonic flowmeter. We then fabricate the waveguide using additive manufacturing and connect it between the transducer and the pipe wall. Transient ultrasonic experiments show that transducers with PnC mountings attain a 40 dB crosstalk reduction in comparison with a standard transducer mounting configuration.

@en

This broad review summarizes recent advances and "hot"research topics in nanophononics and elastic, acoustic, and mechanical metamaterials based on results presented by the authors at the EUROMECH 610 Colloquium held on April 25-27, 2022 in Benicássim, Spain. The key goal of the colloquium was to highlight important developments in these areas, particularly new results that emerged during the last two years. This work thus presents a "snapshot"of the state-of-the-art of different nanophononics- and metamaterial-related topics rather than a historical view on these subjects, in contrast to a conventional review article. The introduction of basic definitions for each topic is followed by an outline of design strategies for the media under consideration, recently developed analysis and implementation techniques, and discussions of current challenges and promising applications. This review, while not comprehensive, will be helpful especially for early-career researchers, among others, as it offers a broad view of the current state-of-the-art and highlights some unique and flourishing research in the mentioned fields, providing insight into multiple exciting research directions.

@en

The published article contains a mistake in the notation of the strain–displacement matrix (Formula presented.) in Voigt notation concerning the isoparametric mapping and in the notation of the corresponding derivatives. As the mistakes relate only to the notation, and the implementation in the code was done correctly, they bear no consequences for the remainder of the manuscript and the results presented therein. The (Formula presented.) matrix defined inline after Eq. (11) on page 5 (repeated on pages 8 and 17) should be written as (Formula presented.) for an enriched element with n original shape functions and m enrichment functions. The differential operator (Formula presented.) in Eq. 12 should then have been defined in Voigt notation and in global coordinates as (Formula presented.) for elastostatics in 2-D and 3-D, respectively, and (Formula presented.) for heat conductivity in 2-D and 3-D, respectively. The derivatives in global coordinates are computed from the derivatives in local coordinates as (Formula presented.) for standard and enriched shape functions, respectively, where (Formula presented.) is the Jacobian of the intersected original element and (Formula presented.) is the Jacobian of the integration element. Furthermore, the derivative of (Formula presented.) with respect to the enriched node location (Formula presented.) in Equation 28 should be written as (Formula presented.) Note that the terms in (Formula presented.) corresponding to the standard shape functions are zero (this equation replaces Eq. 29 of the manuscript): (Formula presented.) as the Jacobian of the parent element is not influenced by the enriched node location, and the derivatives of the standard shape functions with respect to local coordinates of the parent element are constant. The components in (Formula presented.) corresponding to the enrichment functions are computed using (this equation replaces Eq. 30 in the manuscript) (Formula presented.) where the second term is zero because the derivatives of the enrichment function with respect to local coordinates are constant and independent of the enriched node location in global coordinates.

@en

In recent years, the Q-factor of Si 3 N 4 nanomechanical resonators has significantly been increased by soft-clamping techniques using large and complex support structures. To date, however, obtaining similar performance with smaller supports has remained a challenge. Here, we make use of torsion beam supports to tune the Q-factor of Si 3 N 4 string resonators. By design optimization of the supports, we obtain a 50% Q-factor enhancement compared to the standard clamped-clamped string resonators. By performing experimental and numerical studies, we show that further improvement of the Q-factor is limited by a trade-off between maximizing stress and minimizing torsional support stiffness. Thus, our study also provides insight into dissipation limits of high-stress string resonators and outlines how advanced designs can be realized for reaching ultimate f 0 × Q product while maintaining a small footprint.

@en

Mechanical nonlinearities dominate the motion of nanoresonators already at relatively small oscillation amplitudes. Although single and coupled two-degree-of-freedom models have been used to account for experimentally observed nonlinear effects, it is shown that these models quickly deviate from experimental findings when multiple modes influence the nonlinear response. Here, we present a nonlinear reduced-order modeling methodology based on finite-element method simulations for capturing the global nonlinear dynamics of nanomechanical resonators. Our physics-based approach obtains the quadratic and cubic nonlinearities of resonators over a wide frequency range that spans 70 MHz. To qualitatively validate our approach, we perform experiments on a graphene nanodrum driven optothermally and show that the model can replicate diverse ranges of nonlinear phenomena, including multistability, parametric resonance, and different internal resonances without considering any empirical nonlinear fitting parameters. By providing a direct link between microscopic geometry, material parameters, and nonlinear dynamic response, we clarify the physical significance of nonlinear parameters that are obtained from fitting the dynamics of nanomechanical systems, and provide a route for designing devices with desired nonlinear behavior.

@en

Phononic crystals can be designed to show band gaps—ranges of frequencies whose propagation is strongly attenuated in the material. In essence, their working principle is based on destructive interference of waves reflecting from the periodic arrangement of material interfaces (i.e., Bragg scattering). Consequently, capturing accurately the behavior at material interfaces requires appropriate numerical modeling and computational design techniques. However, the commonly used density-based representation in popular topology optimization methods results in a diffuse staircased boundary. The heavily refined finite element meshes required to compensate for this boundary description results in exceedingly large and expensive optimization problems. In this paper, we demonstrate the adverse effect of the density-based boundary description. Furthermore, we propose a level set-based topology optimization procedure with an enriched finite element method that shows improved performance when compared to the density-based approach.

@en

Phononic crystal band gaps (BGs), which are realized by Bragg scattering, have a central frequency and width related to the unit cell's size and the impedance mismatch between material phases. BG tuning has generally been performed by either trial and error or by computational tools such as topology optimization. In either case, understanding how to systematically change the design for a particular band structure is missing. This paper addresses this by closely studying the displacement modes within the wavebands that are responsible for the BG. We look at the variation in different displacement modes due to the changes in the geometry and correlate these changes to their corresponding band structures. We then use this insight to design the unit cell for a particular application, for instance, for generating partial BGs.

@en

The band-gap frequencies of elastic metamaterials are ideally determined by a metamaterial architecture; yet, in practical situations, are often dependent on the material damping in their constituent(s). The analysis of viscoelastic metamaterials requires however substantial computational resources and, except for oversimplified cases, is solely done numerically. Here, we propose an analytical procedure based on the spectral element method (SEM) to analyze bulk metamaterials with viscoelastic damping as continuous systems. Due to intrinsic limitations of the SEM to deal with complex geometries, we develop a procedure to build an approximate model based on SEM frame elements. The viscoelastic behavior is included by means of complex viscoelasticity moduli expressed by the generalized Maxwell mechanical model. We validate this approach by analyzing metamaterial plates and verify the findings experimentally. We demonstrate that our SEM-based analytical model can accurately capture wave transmission around the first band-gap frequencies. Therefore, our extension of the SEM approach to analyze three-dimensional meta-structures is promising to characterize wave propagation in realistic viscoelastic structures (with any type of linear viscoelastic behavior) in an accurate and computationally efficient way.

@en

Elastic metamaterials have been mostly designed to control structural vibration for a specific frequency range. However, the stresses developed due to the wave propagation are not commonly consider in their design phase. This work proposes a multi-objective design optimization of a metamaterial-based interface to tune its dynamic characteristics, while keeping stresses below the allowed value.@en

The performance improvement of renewable energy sources has become nowadays a crucial topic to boost the energy transition. With respect to offshore wind farms, larger wind turbines are currently being designed to increase their power capacity. To support such turbines into the seabed, monopile foundations are mainly used, due to their low-cost manufacturing and transportation. Because of the increased wind turbines’ dimensions, larger monopiles are needed, which have shown to shift high sound waves to lower-frequency ranges during their installation. A precise prediction of the high noise level is necessary to properly design noise-cancelling devices, which prevent negative impacts on marine life. Because the evaluation of sound levels at a full-scale monopile installation is time-consuming and expensive, small-scale tests are being explored. Here we propose a small-scale test for investigating the high noise levels source. To that end, proper scaling laws and instrumentation have been adopted to assure the similarity between the small- and full-scale cases. The experimental results were also compared to a simplified numerical model. The findings of this work showed that investigating a priori the monopile’s features provides an underwater noise prediction, which can be used to design the upcoming noise-cancelling devices.@en

Quantifying the nanomechanical properties of soft-matter using multi-frequency atomic force microscopy (AFM) is crucial for studying the performance of polymers, ultra-thin coatings, and biological systems. Such characterization processes often make use of cantilever's spectral components to discern nanomechanical properties within a multi-parameter optimization problem. This could inadvertently lead to an over-determined parameter estimation with no clear relation between the identified parameters and their influence on the experimental data. In this work, we explore the sensitivity of viscoelastic characterization in polymeric samples to the experimental observables of multi-frequency intermodulation AFM. By performing simulations and experiments we show that surface viscoelasticity has negligible effect on the experimental data and can lead to inconsistent and often non-physical identified parameters. Our analysis reveals that this lack of influence of the surface parameters relates to a vanishing gradient and non-convexity while minimizing the objective function. By removing the surface dependency from the model, we show that the characterization of bulk properties can be achieved with ease and without any ambiguity. Our work sheds light on the sensitivity issues that can be faced when optimizing for a large number of parameters and observables in AFM operation, and calls for the development of new viscoelastic models at the nanoscale and improved computational methodologies for nanoscale mapping of viscoelasticity using AFM.

@en

In this work, an object-oriented geometric engine is proposed to solve problems with discontinuities, for instance, material interfaces and cracks, by means of unfitted, immersed, or enriched finite element methods (FEMs). Both explicit and implicit representations, such as geometric entities and level sets, are introduced to describe configurations of discontinuities. The geometric engine is designed in an object-oriented way and consists of several modules. For efficiency, a (Formula presented.) -d tree data structure that partitions the background mesh is constructed for detecting cut elements whose neighbors are found by means of a dual graph structure. Moreover, the implementation for creating enriched nodes, integration elements, and physical groups is described in detail, and the corresponding pseudo-code is also provided. The complexity and efficiency of the geometric engine are investigated by solving 2-D and 3-D discontinuous models. The capability of the geometric engine is demonstrated on several numerical examples. Topology optimization and problems with intersecting discontinuities are handled with enriched FEMs, where enriched discretizations obtained from the geometric engine are used for the analysis. Furthermore, polycrystalline structures that overlap with an unfitted mesh are considered, where integration elements are created so they align with grain boundaries. Another example shows that the Stanford bunny, which is discretized by a surface mesh with triangular elements, can be fully immersed into a 3-D background mesh. Finally, we share a list of main findings and conclude that the proposed geometric engine is general, robust, and efficient.

@en