An Invariant-Preserving Scheme for the Viscous Burgers-Poisson System

We formulate and analyze a new finite difference scheme for a shallow water model in the form of viscous Burgers-Poisson system with periodic boundary conditions. The proposed scheme belongs to a family of three-level linearized finite difference methods. It is proved to preserve both momentum and energy in the discrete sense. In addition, we proved that the method converges uniformly and has second order of accuracy in space. The analysis given in this work is the first time a pointwise error estimation is done on a second-order finite difference operator applied to the Burgers-Poisson system. We validate our findings by performing various numerical simulations on both viscous and inviscous problems. These numerical examples show the efficacy of the proposed method and confirm the proven theoretical results.

Prediction of the Delamination at the Steel and CFRP Interface of Hybrid Composite Part

The purpose of this study was to predict the adhesive behavior of steel and carbon-fiber-reinforced plastic (CFRP) hybrid parts based on the cohesive zone model (CZM). In this study, the steel sheet and CFRP were joined by epoxy resin in the CFRP prepreg during the curing process, which could generate delamination at their interface because of the springback of steel or the thermal contraction of the CFRP. First, double cantilever beam (DCB) and end-notched flexure (ENF) tests were performed to obtain various adhesion properties such as the critical energy release rate of mode I, mode II (G_I, G_II), and critical stress (σ_max). A finite element (FE) simulation was performed to predict delamination using CZM, which was also used to describe the interfacial behavior between the steel sheet and the CFRP. Finally, a U-shape drawing test was performed for the steel/CFRP hybrid parts, and these results were compared with analytical results.

Spin-orbit interactions of transverse sound

Spin-orbit interactions (SOIs) endow light with intriguing properties and applications such as photonic spin-Hall effects and spin-dependent vortex generations. However, it is counterintuitive that SOIs can exist for sound, which is a longitudinal wave that carries no intrinsic spin. Here, we theoretically and experimentally demonstrate that airborne sound can possess artificial transversality in an acoustic micropolar metamaterial and thus carry both spin and orbital angular momentum. This enables the realization of acoustic SOIs with rich phenomena beyond those in conventional acoustic systems. We demonstrate that acoustic activity of the metamaterial can induce coupling between the spin and linear crystal momentum k, which leads to negative refraction of the transverse sound. In addition, we show that the scattering of the transverse sound by a dipole particle can generate spin-dependent acoustic vortices via the geometric phase effect. The acoustic SOIs can provide new perspectives and functionalities for sound manipulations beyond the conventional scalar degree of freedom and may open an avenue to the development of spin-orbit acoustics.

Probing the phase space of coupled oscillators with Koopman analysis

With the development of probing and computing technology, the study of complex systems has become a necessity in various science and engineering problems, which may be treated efficiently with Koopman operator theory based on observed time series. In the current paper, combined with a singular value decomposition (SVD) of the constructed Hankel matrix, Koopman analysis is applied to a system of coupled oscillators. The spectral properties of the operator and the Koopman modes of a typical orbit reveal interesting invariant structures with periodic, quasiperiodic, or chaotic motion. By checking the amplitude of the principal modes along a straight line in the phase space, cusps of different sizes on the magnitude profiles are identified whenever a qualitative change of dynamics takes place. There seems to be no obstacle to extend the current analysis to high-dimensional nonlinear systems with intricate orbit structures.

A Multiphysics Peridynamic Model for Simulation of Fracture in Si Thin Films during Lithiation/Delithiation Cycles

Material failure is the main obstacle in fulfilling the potential of electrodes in lithium batteries. To date, different failure phenomena observed experimentally in various structures have become challenging to model in numerical simulations. Moreover, their mechanisms are not well understood. To fill the gap, here we develop a coupled chemo-mechanical model based on peridynamics, a particle method that is suitable for simulating spontaneous crack growth, to solve the fracture problems in silicon thin films due to lithiation/delithiation. The model solves mechanical and lithium diffusion problems, respectively, and uses a coupling technique to deal with the interaction between them. The numerical examples of different types of Si films show the advantage of the model in this category and well reproduce the fracture patterns observed in the experiments, demonstrating that it is a promising tool in simulating material failure in electrodes.

On Fractional Geometry of Curves

Fractional Differential Geometry of curves is discussed, with the help of a new fractional derivative, the Λ-fractional derivative, with the corresponding Λ-fractional space. Λ-Fractional derivative completely conforms with the demands of Differential Topology, for the existence of a differential. Therefore Fractional Differential Geometry is established in that Λ-space. The results are pulled back to the initial space.

A Dimension Splitting Generalized Interpolating Element-Free Galerkin Method for the Singularly Perturbed Steady Convection–Diffusion–Reaction Problems

By introducing the dimension splitting method (DSM) into the generalized element-free Galerkin (GEFG) method, a dimension splitting generalized interpolating element-free Galerkin (DS-GIEFG) method is presented for analyzing the numerical solutions of the singularly perturbed steady convection–diffusion–reaction (CDR) problems. In the DS-GIEFG method, the DSM is used to divide the two-dimensional CDR problem into a series of lower-dimensional problems. The GEFG and the improved interpolated moving least squares (IIMLS) methods are used to obtain the discrete equations on the subdivision plane. Finally, the IIMLS method is applied to assemble the discrete equations of the entire problem. Some examples are solved to verify the effectiveness of the DS-GIEFG method. The numerical results show that the numerical solution converges to the analytical solution with the decrease in node spacing, and the DS-GIEFG method has high computational efficiency and accuracy.

Advanced Application of Electrospun Polycaprolactone Fibers for Seed Germination Activity

The increasing intensity of coronavirus (COVID-19) spreading emphasizes the significant development in home food production to reduce the incoming socioeconomic impact from soaring food prices, supply chain fragility, and severe economic crisis. This preliminary study was initiated to demonstrate the possibility of using electrospun fibers as a potential substrate in the application of seed germination activity. The drive of this preliminary study was to integrate the electrospun nanofiber-based material in exploring the current surge in home food production via seed germination in order to introduce cheap source of food without being distracted by the pandemic impact in general. Mung bean (Vigna radiata L. Wilczek) was chosen as it is easy and fast to sprout. Four samples of poly (ε-caprolactone)- (PCL-) based fibers were prepared by means of electrospinning technique, with the optimized flow rate between 0.05 and 0.20 ml/min at a fixed distance of 10 cm needle tip to collector. Mung bean seeds were allowed to germinate on the fabricated electrospun PCL fibers for 96 hours. Our observations include germination percentage, seedling weight, radicle length, and plumule growth. The highest radicle length and plumule length of seedlings were 27.8 mm and 6.7 mm, respectively. There were no inhibitory effects on seed germination and minimal structural fragmentation of smaller diameter electrospun fibers as revealed by FESEM. These results show that the seeds were able to germinate on electrospun PCL fiber substrate, owing to the properties of high surface area and excellent fluid water uptake of PCL fibers.

Sintering Modeling of Thermal Barrier Coatings at Elevated Temperatures: A Review of Recent Advances

To achieve a higher efficiency in gas turbine engine by increasing the inlet-temperature of burning gas is one of the primary goals in aviation industry. The development of thermal barrier coating system (TBCs) continuously raises the inlet-temperature of gas turbine engine in the past decades. Due to the complexity of TBCs and harsh operation environments, the degradation and failure mechanisms of hot section components have not been fully understood, and consequently limits the application of TBCs. It was identified that high-temperature sintering of the topcoat in a typical TBC could be one of the major sources of its failure since the microstructures of the constituent coating layers evolve dynamically during the service period, resulting in significant changes of mechanical and thermal physical properties of the coating system. This paper intends to review recent advances of analytical and numerical modeling of sintering of topcoat in TBCs including the modeling methodology and applications of the models, particularly the implementation of finite element combined with specific materials constitutive functions. Critical remarks on the future development and applications of these models are also discussed in the end.

Viscoplastic plates

An asymptotic model is constructed to describe the bending of thin sheets, or plates, of viscoplastic fluid described by the Herschel–Bulkley constitutive law, which incorporates the von Mises yield condition and a nonlinear viscous stress. The model reduces to a number of previous ones from plasticity theory and viscous fluid mechanics in various limits. It is characterized by a yield criterion proposed by Ilyushin which compactly combines the effect of the bending moment and in-plane stress tensors through three particular invariants. The model is used to explore the bending of loaded flat plates, the deflection of impulsively driven circular plates, and the tension-controlled deflection of loaded beams.

The Hydrodynamics of a Swirling Blood Flow in the Left Heart and Aorta

This paper proposes a new approach to the quantitative analysis of the hydrodynamic structure of a blood flow in the flow channel running from the left atrium to the end of the aorta. This approach is based on the concept of the structural organization of tornado-like swirling jets in channels with a given geometric configuration. Considering the large amount of experimental data in our possession, it was shown that along the entire length of the flow channel, conditions exist for the generation and maintenance of a swirling structure of the jet throughout the entire cardiac cycle. This study has given rise to a new direction in research in fundamental physiology and medicine, which is of great practical importance for diagnosing and treating circulatory disorders accompanied by changes in the geometric configuration and biomechanical characteristics of the heart and great vessels.

Modeling Brittle Fractures in Epoxy Nanocomposites Using Extended Finite Element and Cohesive Zone Surface Methods

Linear elastic fracture modeling coupled with empirical material tensile data result in good quantitative agreement with the experimental determination of mode I fracture for both brittle and toughened epoxy nanocomposites. The nanocomposites are comprised of diglycidyl ether of bisphenol A cured with Jeffamine D-230 and some were filled with core-shell rubber nanoparticles of varying concentrations. The quasi-static single-edge notched bending (SENB) test is modeled using both the surface-based cohesive zone (CZS) and extended finite element methods (XFEM) implemented in the Abaqus software. For each material considered, the critical load predicted by the simulated SENB test is used to calculate the mode I fracture toughness. Damage initiates in these models when nodes at the simulated crack tip attain the experimentally measured yield stress. Prediction of fracture processes using a generalized truncated linear traction–separation law (TSL) was significantly improved by considering the case of a linear softening function. There are no adjustable parameters in the XFEM model. The CZS model requires only optimization of the element displacement at the fracture parameter. Thus, these continuum methods describe these materials in mode I fracture with a minimum number of independent parameters.

Crash Analysis of Aluminum/CFRP Hybrid Adhesive Joint Parts Using Adhesive Modeling Technique Based on the Fracture Mechanics

This study describes the numerical simulation results of aluminum/carbon-fiber-reinforced plastic (CFRP) hybrid joint parts using the explicit finite-element solver LS-DYNA, with a focus on capturing the failure behavior of composite laminates as well as the adhesive capacity of the aluminum–composite interface. In this study, two types of adhesive modeling techniques were investigated: a tiebreak contact condition and a cohesive zone model. Adhesive modeling techniques have been adopted as a widely commercialized model of structural adhesives to simulate adhesive failure based on fracture mechanics. CFRP was studied with numerical simulations utilizing LS-DYNA MAT54 to analyze the crash capability of aluminum/CFRP. To evaluate the simulation model, the results were compared with the force–displacement curve from numerical analysis and experimental results. A parametric study was conducted to evaluate the effect of different fracture toughness values used by designers to predict crash capability and adhesive failure of aluminum/CFRP parts.

Toughening of Bioceramic Composites for Bone Regeneration

Bioceramics are widely considered as elective materials for the regeneration of bone tissue, due to their compositional mimicry with bone inorganic components. However, they are intrinsically brittle, which limits their capability to sustain multiple biomechanical loads, especially in the case of load-bearing bone districts. In the last decades, intense research has been dedicated to combining processes to enhance both the strength and toughness of bioceramics, leading to bioceramic composite scaffolds. This review summarizes the recent approaches to this purpose, particularly those addressed to limiting the propagation of cracks to prevent the sudden mechanical failure of bioceramic composites.

A Comparative Study of the Fractional-Order System of Burgers Equations

This paper is related to the fractional view analysis of coupled Burgers equations, using innovative analytical techniques. The fractional analysis of the proposed problems has been done in terms of the Caputo-operator sense. In the current methodologies, first, we applied the Elzaki transform to the targeted problem. The Adomian decomposition method and homotopy perturbation method are then implemented to obtain the series form solution. After applying the inverse transform, the desire analytical solution is achieved. The suggested procedures are verified through specific examples of the fractional Burgers couple systems. The current methods are found to be effective methods having a close resemblance with the actual solutions. The proposed techniques have less computational cost and a higher rate of convergence. The proposed techniques are, therefore, beneficial to solve other systems of fractional-order problems.

Computational Failure Analysis under Overloading

The aim of this research work is to shed more light on performance-based design through a computational framework that assesses the residual strength of damaged plate-type configurations under overloading. Novel expressions are generated to analyze the power of crack-like stress raisers coupled with retardation effects. Analytical outcomes show that careful consideration of the overload location and crack size can be quite effective in improving safety design and failure mode estimation.

Aerodynamics-assisted, efficient and scalable kirigami fog collectors

To address the global water shortage crisis, one of the promising solutions is to collect freshwater from the environmental resources such as fog. However, the efficiency of conventional fog collectors remains low due to the viscous drag of fog-laden wind deflected around the collecting surface. Here, we show that the three-dimensional and centimetric kirigami structures can control the wind flow, forming quasi-stable counter-rotating vortices. The vortices regulate the trajectories of incoming fog clusters and eject extensive droplets to the substrate. As the characteristic structural length is increased to the size of vortices, we greatly reduce the dependence of fog collection on the structural delicacy. Together with gravity-directed gathering by the folds, the kirigami fog collector yields a collection efficiency of 16.1% at a low wind speed of 0.8 m/s and is robust against surface characteristics. The collection efficiency is maintained even on a 1 m² collector in an outdoor setting.

Water shortage not only occurs in arid regions, but also in humid area with little precipitation, despite abundant fog. Authors develop robust and scalable 3D centimetric kirigami structures to control wind flow and regulate the trajectories of incoming fog, yielding high fog collection efficiency.

Stabilization of the Lower Jamuna River in Bangladesh—Hydraulic and Morphological Assessment

This paper presents a hydraulic and morphological analysis of the Lower Jamuna in Bangladesh with a focus on two key bifurcations that are important for stabilization of the Lower Jamuna reach. We used ground measurements, historical data, multispectral satellite images from various sources as well as numerical models. We carried out hydraulic analyses of the changes and their peculiarities, such as flow distributions at the bifurcation and hysteresis of the stage–discharge relationships. We supplemented our analysis by using numerical models to simulate discharge distribution at the bifurcations under various flow and riverbed conditions. We developed an advanced and automated satellite image processing application for the Lower Jamuna, referred to as Morphology Monitor (MoMo), using the Google Earth Engine. MoMo was found to be an effective tool for a rapid assessment and analysis of the changes in deep-channel and sandbar areas. It is also useful for monitoring and assessing riverbank and char erosion and accretion, which is important not only for morphological but also ecological impact assessment. The application can be adapted as an operational tool as well. Furthermore, we assessed the evolution of deep channels at the bifurcations based on regularly and extensively measured bathymetry data. The analysis was carried out in complement with morphological modeling, particularly for short-term prediction. In this paper we present the major findings of the analysis and discuss their implications for adaptive river management.

Atomic Force Microscope Nanoindentation Analysis of Diffuse Astrocytic Tumor Elasticity: Relation with Tumor Histopathology

This study aims to investigate the influence of isocitrate dehydrogenase gene family (IDH) mutations, World Health Organization (WHO) grade, and mechanical preconditioning on glioma and adjacent brain elasticity through standard monotonic and repetitive atomic force microscope (AFM) nanoindentation. The elastic modulus was measured ex vivo on fresh tissue specimens acquired during craniotomy from the tumor and the peritumoral white matter of 16 diffuse glioma patients. Linear mixed-effects models examined the impact of tumor traits and preconditioning on tissue elasticity. Tissues from IDH-mutant cases were stiffer than those from IDH-wildtype ones among anaplastic astrocytoma patients (p = 0.0496) but of similar elasticity to IDH-wildtype cases for diffuse astrocytoma patients (p = 0.480). The tumor was found to be non-significantly softer than white matter in anaplastic astrocytomas (p = 0.070), but of similar elasticity to adjacent brain in diffuse astrocytomas (p = 0.492) and glioblastomas (p = 0.593). During repetitive indentation, both tumor (p = 0.002) and white matter (p = 0.003) showed initial stiffening followed by softening. Stiffening was fully reversed in white matter (p = 0.942) and partially reversed in tumor (p = 0.015). Tissue elasticity comprises a phenotypic characteristic closely related to glioma histopathology. Heterogeneity between patients should be further explored.

3D pull-out finite element simulation of the pedicle screw-trabecular bone interface at strain rates

Biomedical experimental studies such as pull-out (PO), screw loosening experience variability mechanical properties of fresh bone, legal procedures of cadaver bone samples and time-consuming problems. Finite Element Method (FEM) could overcome experimental problems in biomechanics. However, material modelling of bone is quite difficult, which has viscoelastic and viscoplastic properties. The study presents a bone material model which is constructed at the strain rates with the Johnson-Cook (JC) material model, one of the robust constitutive material models. The JC material constants of trabecular bone are determined by the curve fitting method at strain rates for the 3D PO finite element simulation, which defines the screw-bone interface relationship. The PO simulation is performed using the Abaqus/CAE software program. Bone fracture mechanisms are simulated with dynamic/explicit solutions during the PO phenomenon. The paper exposes whether the strain rate has effects on the PO performance. Moreover, simulation reveals the relationship between pedicle screw diameter and PO performance. The results obtained that the maximum pull-out force (POF) improves as both the screw diameter and the strain rate increase. For 5.5 mm diameter pedicle screw POFs were 487, 517 and 1708 N at strain rate 0.00015, 0.015 and 0.015 s^-1, respectively. The FOFs obtained from the simulation of the other screw were 730, 802 and 2008 N at strain rates 0.00015, 0.0015 and 0.015, respectively. PO phenomenon was also simulated realistically in the finite element analysis (FEA).

On the equilibrium bifurcation of axially deformable holonomic systems: solution of a long-standing enigma

The stability of equilibrium is a fundamental topic in mechanics and applied sciences. Apart from its central role in most engineering fields, it also arises in many natural systems at any scale, from folding/unfolding processes of macromolecules and growth-induced wrinkling in biological tissues to meteorology and celestial mechanics. As such, a few key models represent essential benchmarks in order to gain significant insights into more complex physical phenomena. Among these models, a cornerstone is represented by a structure made of two straight axially deformable bars, connected by an elastic hinge and simply supported at the ends, which are capable of buckling under a compressive axial force. This classical example has been proposed and analysed in some depth by Feodosyev but the attention is here focused on an apparently paradoxical result given by this model, i.e. the existence of a lower bound for the axial-to-flexural stiffness ratio in order for the bifurcation to take place. This enigma is solved theoretically by showing that, differently from other classical stability problems, constitutive and geometric nonlinearities cannot be a priori disconnected and an ideal linearized axial constitutive law cannot be employed in this case. The theory is validated with an experiment, and post-buckling and energy extrema of the proposed solution are discussed as well, highlighting possible snap-back and snap-through phenomena. Finally, the results are extended to the complementary case of tensile buckling.

Universal deformations in inhomogeneous isotropic nonlinear elastic solids

Universal (controllable) deformations of an elastic solid are those deformations that can be maintained for all possible strain-energy density functions and suitable boundary tractions. Universal deformations have played a central role in nonlinear elasticity and anelasticity. However, their classification has been mostly established for homogeneous isotropic solids following the seminal works of Ericksen. In this article, we extend Ericksen’s analysis of universal deformations to inhomogeneous compressible and incompressible isotropic solids. We show that a necessary condition for the known universal deformations of homogeneous isotropic solids to be universal for inhomogeneous solids is that inhomogeneities respect the symmetries of the deformations. Symmetries of a deformation are encoded in the symmetries of its pulled-back metric (the right Cauchy–Green strain). We show that this necessary condition is sufficient as well for all the known families of universal deformations except for Family 5.

Efficient quantum algorithm for dissipative nonlinear differential equations

Nonlinear differential equations model diverse phenomena but are notoriously difficult to solve. While there has been extensive previous work on efficient quantum algorithms for linear differential equations, the linearity of quantum mechanics has limited analogous progress for the nonlinear case. Despite this obstacle, we develop a quantum algorithm for dissipative quadratic n-dimensional ordinary differential equations. Assuming [Formula: see text], where R is a parameter characterizing the ratio of the nonlinearity and forcing to the linear dissipation, this algorithm has complexity [Formula: see text], where T is the evolution time, ϵ is the allowed error, and q measures decay of the solution. This is an exponential improvement over the best previous quantum algorithms, whose complexity is exponential in T. While exponential decay precludes efficiency, driven equations can avoid this issue despite the presence of dissipation. Our algorithm uses the method of Carleman linearization, for which we give a convergence theorem. This method maps a system of nonlinear differential equations to an infinite-dimensional system of linear differential equations, which we discretize, truncate, and solve using the forward Euler method and the quantum linear system algorithm. We also provide a lower bound on the worst-case complexity of quantum algorithms for general quadratic differential equations, showing that the problem is intractable for [Formula: see text] Finally, we discuss potential applications, showing that the [Formula: see text] condition can be satisfied in realistic epidemiological models and giving numerical evidence that the method may describe a model of fluid dynamics even for larger values of R.

A Molecular Simulation Approach to Bond Reorganization in Epoxy Resins: From Curing to Deformation and Fracture

We model bond formation and dissociation processes in thermosetting polymer networks from molecular dynamics simulations. For this, a coarsened molecular mechanics model is derived from quantum calculations to provide effective interaction potentials that enable million-atoms scale simulations. The importance of bond (re)organization is demonstrated for (i) simulating epoxy resin formation-for which our approach leads to realistic network models which can now account for degrees of curing up to 98%. Moreover, (ii) we elucidate the competition of bond dissociation and bond reformation during plastic deformation and fracture. On this basis, we rationalize the molecular mechanisms that account for the irreversible nature of damaging epoxy polymers by mechanical load.

Research on optimization of perforation parameters for formation fractures based on response surface optimization method

For staged multi-cluster fracturing, methods for controlling perforation friction to adjust the flow distribution of each cluster can effectively promote the uniform extension of multiple fractures but lacks a fast and quantitative optimization method for different perforation parameters of each cluster. By establishing a numerical model of single-stage three-cluster flow-limited fracturing under stress-seepage coupling, and based on the response surface optimization method, fully considering the impact of perforation parameters interaction among three perforation clusters, according to the regression equation fitted under the global response, the rapid optimization of perforation parameters of segmented multi-cluster fracturing model is realized. The results show that: in determining the three factors of the study, it is found that there is an obvious interaction between the number of intermediate cluster perforations and the number of cluster perforations on both sides, the number of cluster perforations on both sides and the diameter of intermediate cluster perforations, the response surface optimization method gives the optimal perforation parameter combination of three clusters of fractures under global response; When the perforation parameters were combined before optimization, the fracture length difference was 32.550m, and the intermediate perforation cluster evolved into invalid perforation cluster, when the perforation parameters were combined after optimization, the fracture length difference was 0.528m, the three perforation clusters spread uniformly, and there are no invalid clusters. At the same time, the regression equation under the response is optimized before and after the comparison between the predicted value of the equation and the actual simulation value. It is found that the estimated deviation rate of the equation before optimization is 1.2%, and the estimated deviation rate after optimization is 0.4%. The estimated deviation rates are all less, and the response regression equation based on the response surface optimization method can quickly optimize the perforation parameters. The response surface optimization method is suitable for the multi parameter optimization research of formation fracturing which is often affected by many geological and engineering factors. Combining with the engineering practice and integrating more factors to optimize the hydraulic fracturing parameters, it is of great significance to improve the success rate of hydraulic fracturing application.

Mathematical Methodology and Metallurgical Application of Turbulence Modelling: A Review

This paper focusses on three main numerical methods, i.e., the Reynolds-Averaged Navier-Stokes (RANS), Large Eddy Simulation (LES), and Direct Numerical Simulation (DNS) methods. The formulation and variation of different RANS methods are evaluated. The advantage and disadvantage of RANS models to characterize turbulent flows are discussed. The progress of LES with different subgrid scale models is presented. Special attention is paid to the inflow boundary condition for LES modelling. Application and limitation of the DNS model are described. Different experimental techniques for model validation are given. The consistency between physical experimentation/modelling and industrial cases is discussed. An emphasis is placed on the model validation through physical experimentation. Subsequently, the application of a turbulence model for three specific flow problems commonly encountered in metallurgical process, i.e., bubble-induced turbulence, supersonic jet transport, and electromagnetic suppression of turbulence, is discussed. Some future perspectives for the simulation of turbulent flow are formulated.

Tunable stiffness enables fast and efficient swimming in fish-like robots

Fish maintain high swimming efficiencies over a wide range of speeds. A key to this achievement is their flexibility, yet even flexible robotic fish trail real fish in terms of performance. Here, we explore how fish leverage tunable flexibility by using their muscles to modulate the stiffness of their tails to achieve efficient swimming. We derived a model that explains how and why tuning stiffness affects performance. We show that to maximize efficiency, muscle tension should scale with swimming speed squared, offering a simple tuning strategy for fish-like robots. Tuning stiffness can double swimming efficiency at tuna-like frequencies and speeds (0 to 6 hertz; 0 to 2 body lengths per second). Energy savings increase with frequency, suggesting that high-frequency fish-like robots have the most to gain from tuning stiffness.

Enabling novel dispersion and topological characteristics in mechanical lattices via stable negative inertial coupling

Mechanical topological insulators have enabled a myriad of unprecedented characteristics that are otherwise not conceivable in traditional periodic structures. While rich in dynamics, new developments in the domain of mechanical topological systems are hindered by their inherent inability to exhibit negative elastic or inertial couplings owing to the inevitable loss of dynamical stability. The aim of this paper is, therefore, to remedy this challenge by introducing a class of architected inertial metamaterials (AIMs) as a platform for designing mechanical lattices with novel topological and dispersion traits. We show that carefully coupling elastically supported masses via moment-free rigid linkages invokes a dynamically stable negative inertial coupling, which is essential for topological classes in need of such negative interconnection. The potential of the proposed AIMs is demonstrated via three examples: (i) a mechanical analogue of Majorana edge states, (ii) a square diatomic AIM that can sustain the quantum valley Hall effect (classically arising in hexagonal lattices), and (iii) a square tetratomic AIM with topological corner modes. We envision that the presented framework will pave the way for a plethora of robust topological mechanical systems.

Electropolishing—A Practical Method for Accessing Voids in Metal Films for Analyses

In many applications, voids in metals are observed as early degradation features caused by fatigue. In this publication, electropolishing is presented in the context of a novel sample preparation method that is capable of accessing voids in the interior of metal thin films along their lateral direction by material removal. When performed at optimized process parameters, material removal can be well controlled and the surface becomes smooth at the micro scale, resulting in the voids being well distinguishable from the background in scanning electron microscopy images. Compared to conventional cross-sectional sample preparation (embedded mechanical cross-section or focused ion beam), the accessed surface is not constrained by the thickness of the investigated film and laterally resolved void analyses are possible. For demonstrational purposes of this method, the distribution of degradation voids along the metallization of thermo-mechanically stressed microelectronic chips has been quantified.

Evaluation and Validation of Thorax Model Responses: A Hierarchical Approach to Achieve High Biofidelity for Thoracic Musculoskeletal System

As one of the most frequently occurring injuries, thoracic trauma is a significant public health burden occurring in road traffic crashes, sports accidents, and military events. The biomechanics of the human thorax under impact loading can be investigated by computational finite element (FE) models, which are capable of predicting complex thoracic responses and injury outcomes quantitatively. One of the key challenges for developing a biofidelic FE model involves model evaluation and validation. In this work, the biofidelity of a mid-sized male thorax model has been evaluated and enhanced by a multi-level, hierarchical strategy of validation, focusing on injury characteristics, and model improvement of the thoracic musculoskeletal system. At the component level, the biomechanical responses of several major thoracic load-bearing structures were validated against different relevant experimental cases in the literature, including the thoracic intervertebral joints, costovertebral joints, clavicle, sternum, and costal cartilages. As an example, the thoracic spine was improved by accurate representation of the components, material properties, and ligament failure features at tissue level then validated based on the quasi-static response at the segment level, flexion bending response at the functional spinal unit level, and extension angle of the whole thoracic spine. At ribcage and full thorax levels, the thorax model with validated bony components was evaluated by a series of experimental testing cases. The validation responses were rated above 0.76, as assessed by the CORA evaluation system, indicating the model exhibited overall good biofidelity. At both component and full thorax levels, the model showed good computational stability, and reasonable agreement with the experimental data both qualitatively and quantitatively. It is expected that our validated thorax model can predict thorax behavior with high biofidelity to assess injury risk and investigate injury mechanisms of the thoracic musculoskeletal system in various impact scenarios. The relevant validation cases established in this study shall be directly used for future evaluation of other thorax models, and the validation approach and process presented here may provide an insightful framework toward multi-level validating of human body models.

The Rheological Analytical Solution and Parameter Inversion of Soft Soil Foundation

In soft soil engineering projects, the building loads are always required to be symmetrically distributed on the surface of the foundation to prevent uneven settlement. Even if the buildings and soft clay are controlled by engineers, it can still lead to the rheology of the foundation. The analytical solution based on the Laplace integral transformation method has positive significance for providing a simple and highly efficient way to solve engineering problems, especially in the long-term uneven settlement deformation prediction of buildings on soft soil foundations. This paper proposes an analytical solution to analyze the deformation of soft soil foundations. The methodology is based on calculus theory, Laplace integral transformation, and viscoelastic theory. It combines an analytical solution with finite theory to solve the construction sequences and loading processes. In addition, an improved quantum genetic algorithm is put forward to inverse the parameters of soft soil foundations. The analytical solution based on Laplace integral transformation is validated through an engineering case. The results clearly illustrate the accuracy of the method.

Reduced Order Multiscale Simulation of Diffuse Damage in Concrete

Damage in concrete structures initiates as the growth of diffuse microcracks that is followed by damage localisation and eventually leads to structural failure. Weak changes such as diffuse microcracking processes are failure precursors. Identification and characterisation of these failure precursors at an early stage of concrete degradation and application of suitable precautionary measures will considerably reduce the costs of repair and maintenance. To this end, a reduced order multiscale model for simulating microcracking-induced damage in concrete at the mesoscale level is proposed. The model simulates the propagation of microcracks in concrete using a two-scale computational methodology. First, a realistic concrete specimen that explicitly resolves the coarse aggregates in a mortar matrix was generated at the mesoscale. Microcrack growth in the mortar matrix is modelled using a synthesis of continuum micromechanics and fracture mechanics. Model order reduction of the two-scale model is achieved using a clustering technique. Model predictions are calibrated and validated using uniaxial compression tests performed in the laboratory.

Traction-separation laws of graphene grain boundaries

Molecular dynamics simulations are used to extract the traction-separation laws (TSLs) of symmetric grain boundaries of graphene. Grain boundaries with realistic atomic structures are constructed using different types of dislocations. The TSLs of grain boundaries are extracted by using cohesive zone volume elements (CZVEs) ahead of the crack tip. The traction and separation of each cohesive zone volume element are calculated during the crack growth. The traction and separation values obtained for the cohesive elements predict that the TSLs of grain boundaries have a bilinear form. The areas under the traction-separation curves are used to calculate the separation energy of the grain boundaries. The results show that as the grain boundary misorientation angle increases the separation energy of the grain boundaries decreases. The impact of temperature on the traction separation laws is studied. The results show that, with an increase of the temperature from 0.1 K to 300 K, the separation energy first increases to reach its peak at around 25 K and then slightly decreases.

Electrospun Polylactide/Natural Rubber Fibers: Effect Natural Rubber Content on Fiber Morphology and Properties

Non-woven polylactide-natural rubber fiber materials with a rubber content of 5, 10 and 15 wt.% were obtained by electrospinning. The thermal, dynamic, and mechanical properties of the fibers were determined. It was shown that the average fiber diameter increased with adding of the NR content, while the linear and surface densities changed slightly. Using the differential scanning calorimetry, the thermal characteristics were obtained. It was found that the glass transition temperature of polylactide increased by 2–5 °C, and the melting temperature increased by 2–4 °C in the presence of natural rubber in the samples. By the method of electronic paramagnetic resonance at T = 50 and 70 °C it was determined that the mobility of the amorphous phase in PLA/NR fibers increased with the addition of NR. The adding of NR at a content of 15 wt.% increased the value of elongation at break by 3.5 times compared to pure PLA.

Fracture and Size Effect of PFRC Specimens Simulated by Using a Trilinear Softening Diagram: A Predictive Approach

The size effect on plain concrete specimens is well known and can be correctly captured when performing numerical simulations by using a well characterised softening function. Nevertheless, in the case of polyolefin-fibre-reinforced concrete (PFRC), this is not directly applicable, since using only diagram cannot capture the material behaviour on elements with different sizes due to dependence of the orientation factor of the fibres with the size of the specimen. In previous works, the use of a trilinear softening diagram proved to be very convenient for reproducing fracture of polyolefin-fibre-reinforced concrete elements, but only if it is previously adapted for each specimen size. In this work, a predictive methodology is used to reproduce fracture of polyolefin-fibre-reinforced concrete specimens of different sizes under three-point bending. Fracture is reproduced by means of a well-known embedded cohesive model, with a trilinear softening function that is defined specifically for each specimen size. The fundamental points of these softening functions are defined a priori by using empirical expressions proposed in past works, based on an extensive experimental background. Therefore, the numerical results are obtained in a predictive manner and then compared with a previous experimental campaign in which PFRC notched specimens of different sizes were tested with a three-point bending test setup, showing that this approach properly captures the size effect, although some values of the fundamental points in the trilinear diagram could be defined more accurately.

Hydrogen-Assisted Brittle Fracture Behavior of Low Alloy 30CrMo Steel Based on the Combination of Experimental and Numerical Analyses

Compact-tension (CT) specimens made of low alloy 30CrMo steels were hydrogen-charged, and then subjected to the fracture toughness test. The experimental results revealed that the higher crack propagation and the lower crack growth resistance (CTOD-R curve) are significantly noticeable with increasing hydrogen embrittlement (HE) indexes. Moreover, the transition in the microstructural fracture mechanism from ductile (microvoid coalescence (MVC)) without hydrogen to a mixed quasi-cleavage (QC) fracture and QC + intergranular (IG) fracture with hydrogen was observed. The hydrogen-enhanced decohesion (HEDE) mechanism was characterized as the dominant HE mechanism. According to the experimental testing, the coupled problem of stress field and hydrogen diffusion field with cohesive zone stress analysis was employed to simulate hydrogen-assisted brittle fracture behavior by using ABAQUS software. The trapezoidal traction-separation law (TSL) was adopted, and the initial TSL parameters from the best fit to the load-displacement and J-integral experimental curves without hydrogen were calibrated for the critical separation of 0.0393 mm and the cohesive strength of 2100 MPa. The HEDE was implemented through hydrogen influence in the TSL, and to estimate the initial hydrogen concentration based on matching numerical and experimental load-line displacement curves with hydrogen. The simulation results show that the general trend of the computational CTOD-R curves corresponding to initial hydrogen concentration is almost the same as that obtained from the experimental data but in full agreement, the computational CTOD values being slightly higher. Comparative analysis of numerical and experimental results shows that the coupled model can provide design and prediction to calculate hydrogen-assisted fracture behavior prior to extensive laboratory testing, provided that the material properties and properly calibrated TSL parameters are known.

Indentation mapping of stretched adhesive membranes

Unilateral adhesive contact between a rigid indenter and a uniformly stretched membrane of arbitrary shape is considered. The generalized Johnson–Kendall–Roberts (JKR)-type and Derjaguin– Muller–Toporov (DMT)-type models of non-axisymmetric adhesive contact are presented for short- and long-range adhesion, respectively, and the JKR–DMT transition is established in the framework of the generalized Maugis–Dugdale model. A refined method of matched asymptotic expansions is applied to construct the leading-order asymptotic model for indentation mapping of freestanding two-dimensional materials with an axisymmetric probe, using the approximate analytical solution obtained in explicit form for an infinite membrane in the limit of short-range adhesive contact with an indenter in the form of an elliptic paraboloid. The cases of a spherical indenter and a rectangular membrane are studied in detail.

Contrasting the Role of Pores on the Stress State Dependent Fracture Behavior of Additively Manufactured Low and High Ductility Metals

This study investigates the disparate impact of internal pores on the fracture behavior of two metal alloys fabricated via laser powder bed fusion (L-PBF) additive manufacturing (AM)—316L stainless steel and Ti-6Al-4V. Data from mechanical tests over a range of stress states for dense samples and those with intentionally introduced penny-shaped pores of various diameters were used to contrast the combined impact of pore size and stress state on the fracture behavior of these two materials. The fracture data were used to calibrate and compare multiple fracture models (Mohr-Coulomb, Hosford-Coulomb, and maximum stress criteria), with results compared in equivalent stress (versus stress triaxiality and Lode angle) space, as well as in their conversions to equivalent strain space. For L-PBF 316L, the strain-based fracture models captured the stress state dependent failure behavior up to the largest pore size studied (2400 µm diameter, 16% cross-sectional area of gauge region), while for L-PBF Ti-6Al-4V, the stress-based fracture models better captured the change in failure behavior with pore size up to the largest pore size studied. This difference can be attributed to the relatively high ductility of 316L stainless steel, for which all samples underwent significant plastic deformation prior to failure, contrasted with the relatively low ductility of Ti-6Al-4V, for which, with increasing pore size, the displacement to failure was dominated by elastic deformation.

Curvature-controlled delamination patterns of thin films on spherical substrates

Periodic delamination patterns in multilayer structures have exhibited extensive applications in microelectronics and optics devices. However, delamination behaviors of a closed thin shell on spherical substrates are still elusive. Herein, a unique instability mechanism of buckle delamination in a closed thin film weakly bonded to spherical substrates is studied by experiments, simulations, and theoretical analyses. The system of an Al film depositing on polystyrene spheres subjected to thermal mismatch strain is used for demonstration. Unlike traditional phenomena of wrinkling and wrinkle-induced delamination under increasing misfit strain, the weak adhesion between the core and shell results in a periodic pattern of delaminated hexagonal dimples that emerges directly from the smooth sphere configuration, before which no wrinkling occurs. Both substrate curvature and interfacial adhesion are revealed to control the dimple size and delamination width. These findings open a new venue for manifesting new controllable features for surface microfabrication.

Lignin Nanoparticles and Their Nanocomposites

Lignin nanomaterials have emerged as a promising alternative to fossil-based chemicals and products for some potential added-value applications, which benefits from their structural diversity and biodegradability. This review elucidates a perspective in recent research on nanolignins and their nanocomposites. It summarizes the different nanolignin preparation methods, emphasizing anti-solvent precipitation, self-assembly and interfacial crosslinking. Also described are the preparation of various nanocomposites by the chemical modification of nanolignin and compounds with inorganic materials or polymers. Additionally, advances in numerous potential high-value applications, such as use in food packaging, biomedical, chemical engineering and biorefineries, are described.