Bubbles rising through otherwise quiescent liquid present a challenge in fluid mechanics: A gas volume is released, buoyancy drives it upward, and one might expect the path to remain nearly vertical unless the surrounding liquid is disturbed. In practice, freely rising bubbles often develop lateral motion, moving along zigzag or spiral trajectories while their shapes deform and their wakes become unsteady. Bubble trajectory affects residence time, interfacial renewal, local mixing, momentum exchange, and the distribution of gas within a liquid phase. The scientific difficulty is that the observed motion is produced by several coupled processes occurring at once: the bubble deforms, liquid flows around the interface, the wake evolves, and hydrodynamic forces shift in direction as the bubble changes orientation. A long-standing interpretation has treated wake instability and vortex shedding as central to the onset of non-vertical bubble paths.  The vortex shed behind the bubble can coexist with path oscillation, but that does not by itself explain how the bubble surface is steered at the instant when its lateral motion changes direction. For a deformable bubble, the interface is both a boundary enclosing the gas phase and part of the dynamics. Local flow along the bubble surface can redistribute pressure and shear, reshape the lower interface, and shift the direction of the net hydrodynamic action.

In a recent research paper published in International Journal of Heat and Fluid Flow, Dr. He Liu, Dr. Yajing Yang, and Professor Yanju Wei from Xi’an Jiaotong University,  examined a surface-flow-based interpretation of bubble path instability in quiescent water. They identified alternating dominance between clockwise and counter-clockwise interfacial streams, formed after bypass-flow separation near the stagnation point, as the driver of lateral motion and rolling torque. The technically distinct contribution is the relocation of the causal mechanism from downstream vortex shedding to local reflux and counter-flow interaction along the bubble’s lower surface. Briefly, the research team combined high-speed experimental observation with interface-resolved numerical simulation. In their experiments, air bubbles were released from a needle into distilled water and recorded from two orthogonal directions, allowing the bubble centroid and shape evolution to be tracked in three dimensions. The numerical model, implemented for laminar incompressible two-phase flow with surface tension, was used to resolve the surrounding velocity field and the interfacial motion in greater detail than the optical measurements alone could provide. This pairing mattered because the path itself shows only the global consequence; the proposed mechanism depends on what happens locally along the bubble surface.

After detachment, the bubble initially rose almost vertically. Once it had ascended a finite distance, the trajectory departed from a straight path and developed into zigzag, oblique zigzag, spiral, or transitional forms. The reported time-averaged rising velocity was about 310 mm/s, and the measured oscillation frequencies showed a clear relation among the motion components: the frequency associated with vertical oscillation was nearly twice that of the lateral components. That relationship is consistent with a bubble whose lateral motion reverses over a repeated deformation cycle rather than one undergoing a purely random sideward drift.

The authors found that during a zigzag cycle, the lower bubble surface did not retain a fixed geometry. It evolved from a backslash-like profile, into a V-shaped lower surface, and then into a forward-slash-like profile. These shape changes were not described as passive distortions. They were tied directly to lateral steering. When the lower surface assumed the V-like form, its two arms acted as steering structures that redirected the bubble’s motion. A specific design choice in the analysis, subtracting the bubble centroid velocity to examine the relative velocity field, changed the physical interpretation: it separated translational motion from local rotational and surface-flow behavior, making the competing interfacial streams visible.

In this relative frame, the bypass flow separated at a stagnation point and produced two counter-rotating surface streams. One travelled counter-clockwise along one side of the interface; the other travelled clockwise along the opposite side. They alternately strengthened and weakened over the cycle. When the counter-clockwise component dominated, it promoted rightward translation and leftward rolling; when the clockwise component dominated, it promoted leftward translation and rightward rolling. The bubble’s quasi-sinusoidal lateral motion therefore arose from alternating dominance between these two surface-flow components.

The team also carried out simulation studies to investigate the relation between surface flow and vortex shedding and noted that counter-rotating surface flows converged near the lower part of the bubble, and the stronger stream could push past the lower stagnation region, impinge on the weaker one, and generate a reflux zone. This local interaction produced sharp changes in surface curvature and inflection points on the bubble. It also displaced portions of the weaker stream away from the interface, contributing to vortex detachment into the wake.   They appear as a consequence of the interfacial competition and flow reversal, rather than as the primary origin of the side-to-side motion. The same logic was extended to spiral motion. When the competition between the clockwise and counter-clockwise surface streams remained in the vertical plane, the bubble followed a zigzag path. When the competing motion was redirected into the horizontal plane, the interaction became a chasing-like motion around the bubble and generated a spiral trajectory. Experiments across bubbles of different initial diameters showed zigzag, spiral, zigzag-to-spiral transition, oblique zigzag, and steady vertical ascent, with no strict one-to-one relation between bubble size and trajectory class. The distinction between zigzag and spiral motion was therefore interpreted as a difference in the orientation of the same surface-flow competition, not as evidence for fundamentally separate mechanisms.

The findings of Professor Yanju Wei and colleagues have practical relevance for engineering systems in which bubbles are not simply dispersed gas volumes but moving, deforming hydrodynamic objects that influence transport performance. In gas–liquid contactors, bubble columns, chemical reactors, flotation devices, and thermal-fluid equipment, designers often rely on empirical descriptions of bubble rise velocity, residence time, interfacial area, and mixing intensity.   If the lateral motion of a bubble is driven by alternating clockwise and counter-clockwise surface streams, then trajectory control should not be approached solely through wake suppression or bulk turbulence management. It also requires attention to conditions that modify interfacial mobility, lower-surface reflux, stagnation-point behavior, and bubble deformation.

One implication we believe is especially important for reactor and heat-transfer design: bubble path instability can enhance lateral displacement and local liquid agitation even in otherwise quiescent liquid. A bubble that zigzags or spirals sweeps a larger volume than one rising vertically, which may improve local mixing, gas–liquid contact, and renewal of the liquid near the interface. At the same time, because the mechanism depends on interfacial flow competition, small changes in surface condition, fluid cleanliness, bubble size distribution, or confinement may alter the degree of lateral wandering.  The work also has implications for numerical modelling of bubbly flows. Many engineering-scale simulations cannot resolve each deforming interface, so they depend on closure relations for lift, drag, path oscillation, and dispersion. A wake-based interpretation may miss the timing and origin of the lateral force if the relevant event begins at the lower bubble surface. The authors’ analysis gives model developers a more physically specific target: the competition of counter-rotating surface flows and the resulting rolling torque.  For process control, the findings suggest that trajectory type—zigzag or spiral—should be regarded less as a fixed bubble-size category and more as a geometric expression of the same interfacial instability. The study by Liu, Yang, and Wei therefore provides a useful design insight: controlling bubble motion may require controlling how bypass flow is redirected along the interface, not just adjusting gas injection rate or relying on average bubble diameter.

Spacecraft operating in coordinated groups depend on more than the stability of individual vehicles. In formation-based missions, each spacecraft must regulate its attitude in relation to a shared reference or leader while exchanging only limited information with neighboring spacecraft. This requirement makes attitude consensus control a central problem in distributed spacecraft systems, particularly when coordinated pointing, formation maintenance, cooperative observation, or communication tasks depend on synchronized orientation. A distributed strategy is especially valuable because it avoids reliance on a central information source and allows each spacecraft to act using local communication links, which is more consistent with the practical constraints of multi-spacecraft operation.  The problem becomes more difficult when the spacecraft are not idealized as rigid bodies. Real spacecraft often carry flexible appendages such as solar arrays or antennas, and these structures can vibrate during attitude maneuvers. Such couples with the rotational motion of the spacecraft body and can degrade pointing accuracy and closed-loop performance. When this flexibility is represented directly, the spacecraft is more appropriately described as a distributed-parameter system governed by partial differential equations rather than by a finite set of rigid-body ordinary differential equations. This distinction is important because simplified modal truncation can neglect higher-frequency modes and may introduce spillover-related stability concerns. A control law designed for rigid spacecraft therefore cannot be transferred directly to flexible spacecraft without accounting for the interaction between hub attitude motion and appendage deformation.  Time performance adds another layer to the control challenge. Many attitude consensus strategies guarantee only asymptotic convergence, meaning that the tracking error approaches a neighborhood of the desired state as time tends to infinity. For mission planning, this is often not sufficiently explicit. Finite-time and fixed-time control methods provide stronger guarantees, but prescribed-time control is more directly useful because the convergence time can be selected in advance. In a multi-spacecraft system, such a capability allows the designer to specify not only how accurately the followers should track the leader, but also when that tracking accuracy must be reached. Extending this idea to flexible spacecraft, however, is mathematically demanding because the controller must impose prescribed timing while also keeping structural vibration bounded.

In a recent research paper published in Advances in Space Research, Professor Shilei Cao and Professor Jian Liu from Harbin Institute of Technology, together with Man Yang from HIT Satellite Technology Co., Ltd., developed a distributed hierarchical control framework for multiple flexible spacecraft. The strategy combines a prescribed-time distributed estimator for leader attitude and angular velocity with a local adaptive prescribed-performance controller based on state transformation. Its distinct technical feature is the joint enforcement of predefined convergence time and tracking accuracy while proving bounded appendage vibration under flexible rigid-body coupling. The controller uses adaptive laws and a Nussbaum-type function to handle unknown disturbances and time-varying actuator efficiency and bias drift.

The researchers formulated each flexible spacecraft as a hub-and-appendage system, with the attitude angle describing rigid-body orientation and the appendage deformation represented over its spatial length. Hamilton’s principle was used to derive the governing PDE-based model, including the bending stiffness, linear density, central-body inertia, structural damping, control input, and external disturbance. The network was described by a directed graph with a virtual leader supplying the desired attitude reference. The graph assumption required that at least one follower receive leader information and that every spacecraft have a directed path from the leader. This is a modest but important communication condition: it allows leader information to propagate through local links rather than requiring direct leader access for every follower.

The authors organized control design in two layers. In the estimator layer, each follower spacecraft constructed prescribed-time estimates of the leader’s attitude and angular velocity from neighbor information. They used graph-theoretic properties of the communication matrix and prescribed-time Lyapunov arguments to show that the estimation errors vanish within a predefined time. This part of the design is essential because the local controller is built on the estimated leader states rather than assuming that every follower directly measures the leader. By separating leader-state reconstruction from local tracking, the strategy preserves the distributed structure while still giving each spacecraft the reference information needed for tracking.

The local control layer then addressed the prescribed-performance tracking problem. The team introduced a performance boundary so that the tracking error would remain within a user-defined envelope and reach the final allowable accuracy within a specified settling time. Rather than enforcing this constrained error directly, the researchers transformed the constrained attitude error into an unconstrained variable. That design choice has a clear scientific consequence: it converts a prescribed transient-and-steady-state performance requirement into a smooth nonlinear control problem, making it possible to prove boundedness and performance preservation through Lyapunov analysis.

The adaptive controller also had to deal with time-varying actuator faults and lumped disturbances. To do this, the design incorporated a Nussbaum-type function and adaptive laws for unknown bounds, while additional terms compensated for rigid-flexible coupling between the hub attitude and appendage vibration. The stability proof combined the transformed tracking-error Lyapunov function with energy terms for the flexible appendage. A cross-energy term was introduced to handle the coupling that arises naturally in the derivative of the appendage and attitude energy. This is one of the more technically revealing parts of the research work, because it treats the vibration dynamics as part of the control proof rather than as an afterthought.

The authors performed numerical simulations and used a four-follower spacecraft network and a virtual leader under time-varying actuator faults and external disturbances. The estimator allowed the followers to recover the leader’s attitude and angular velocity within the predefined estimation time. Under the proposed control scheme, the followers tracked the leader with errors entering the prescribed bounds within the assigned 15-second settling time. Comparisons with distributed asymptotic control, distributed adaptive fault-tolerant control, and a predefined-time rigid-spacecraft control method clarified the contribution. The proposed method maintained prescribed tracking accuracy, kept appendage vibrations bounded and convergent, and did not require additional tip control forces, while the adaptive variables and control commands remained bounded in the simulations.

The engineering value of the research team findings is clearest for spacecraft formations in which orientation must be coordinated within a known operational time window. The control strategy developed by Cao, Liu, and Yang is relevant to these systems because it does not just drive the followers toward consensus eventually; it allows the convergence time and tracking accuracy to be specified in advance.   A controller with prescribed-time behavior can be aligned with observation windows, coordinated pointing sequences, or formation reconfiguration periods where attitude synchronization must be completed before the next mission phase begins. The study’s simulations show the follower spacecraft tracking the leader within a predefined 15-second settling time, while keeping the tracking errors inside prescribed performance bounds.  A second application lies in the control of spacecraft with large or lightly damped flexible appendages, including solar panels, antennas, booms, and other deployable structures.   The authors’ PDE-based treatment is therefore useful for systems where flexible motion cannot be safely reduced to a simple rigid-body approximation.

The findings are also applicable to fault-tolerant spacecraft operation. Actuator degradation, efficiency loss, and bias drift can compromise attitude control, especially in multi-spacecraft systems where one vehicle’s tracking error may affect coordinated behavior. The proposed controller explicitly considers time-varying actuator efficiency and bias drift rather than assuming fixed fault parameters.   The use of adaptive techniques and a Nussbaum-type function give the controller a mechanism to maintain prescribed-performance tracking despite uncertain actuator effectiveness.  The distributed nature of the strategy gives it further engineering relevance for larger formations. Each follower estimates the leader’s attitude and angular velocity using local neighbor information, reducing reliance on global communication. The authors’ layered estimator-controller design therefore supports a practical architecture: local information exchange first reconstructs the leader state, then each spacecraft performs prescribed-performance tracking while accounting for flexibility and actuator faults.

Small hydropower offers a practical route for generating electricity from local water resources, but its success depends strongly on turbine reliability under real operating conditions. Small-scale facilities are different from large installations, because usually they rely on simpler supporting infrastructure and have less room for costly debris-removal equipment. In such settings, the hydraulic machine is expected to extract energy efficiently and also to tolerate water conditions that may include foreign matter. This requirement places unusual pressure on runner design, because the same blade passages that guide the flow and transfer angular momentum can also become sites of blockage. When debris accumulates inside a small turbine, the problem is not just a local mechanical inconvenience; it directly threatens one of the principal advantages of hydropower, namely the ability to provide stable power from a continuous hydraulic resource.

A single blade centrifugal sewage pump offers an interesting route around this constraint because its wide internal passage is naturally suited to the transport of foreign matter. When operated in reverse, such a pump can function as a turbine, creating the possibility of a non-blocking water turbine based on an existing pump architecture. That idea, however, brings its own fluid-dynamic complications. A single blade runner is strongly asymmetric, and earlier studies have treated radial thrust and whirling as central reliability concerns. In turbine mode, another issue becomes equally important: the internal flow path was originally designed for pump operation. The blade angles, volute geometry, runner passage, clearances, and outlet flow arrangement therefore encounter a reversed operating condition in which efficiency and performance stability become strongly dependent on the internal flow.

The scientific challenge is to move beyond average turbine performance and identify how hydraulic losses are generated inside this strongly asymmetric turbomachine. For conventional multi-blade pump-as-turbine systems, much prior work has dealt with predicting turbine-mode performance from pump-mode behavior and with improving performance through blade or runner modifications. For the single blade reverse-running case, the analysis is more demanding because the flow structure changes markedly with runner rotation, so losses cannot be understood as steady, circumferentially uniform quantities. The relationship between turbine-mode internal flow and hydraulic loss generation had not been clarified, which leave an important gap between the practical appeal of the non-blocking runner and the hydraulic knowledge needed to improve it.

In a recent research paper published in International Journal of Heat and Fluid Flow, Professor Yasuyuki Nishi and Ms. Natsumi Itoh from Ibaraki University, developed a turbine-mode hydraulic loss analysis method for a single blade reverse running pump turbine. The method separates effective head into theoretical head and individual hydraulic losses, including runner friction loss, runner loss, casing loss, inlet pipe loss, and outlet pipe loss. Its distinct contribution is that it uses unsteady CFD data over one runner rotation to evaluate both time-averaged and instantaneous loss components. This allowed the authors to connect phase-dependent hydraulic losses directly to suction-surface separation, outlet vortices, casing pressure behavior, and blade-inlet backflow.

Nishi and Itoh in their study combined performance measurements, particle image velocimetry at the blade inlet, unsteady three-dimensional CFD, and a turbine-mode loss analysis method. The test machine they used was a closed single blade centrifugal runner installed in a volute casing and operated in reverse at 900 min⁻¹. The experimental program measured head, torque, output, and efficiency while also resolving the circumferential and radial velocity components near the blade inlet. The numerical model reproduced the turbine geometry, including the volute, runner, inlet and outlet pipes, and clearances around the shrouds, so that the computation could be used as a performance predictor and also a spatially resolved diagnostic of loss formation.

The authors’ CFD results reproduced the measured performance with reasonable agreement, especially near the maximum efficiency flow coefficient of φ = 0.051, where the experimental efficiency was about 0.560. The circumferential component of the absolute velocity at the blade inlet, which dominated the blade inlet flow, was also captured well. Differences appeared in the radial component at some circumferential positions and phases, but the main swirl-dominated character of the inlet flow was represented sufficiently for the authors’ loss analysis.

The team found turbine did not behave as a quasi-steady axisymmetric machine and for instance at the maximum efficiency flow rate, the head coefficient, output coefficient, and efficiency changed strongly with blade phase angle. The head and output coefficients reached their largest values near θ* = 40° and then decreased toward a minimum near θ* = 325°. Efficiency followed a related but not identical pattern, reaching its maximum near θ* = 206° because the output coefficient had a local maximum around θ* = 185°. This phase dependence is central to the study and shows that the single blade runner creates an uneven flow field and also continuously reorganizes the balance between theoretical head, hydraulic loss, and power extraction as it rotates.

The loss analysis separated the total hydraulic loss into runner friction loss, runner loss other than friction, casing loss, inlet pipe loss, and outlet pipe loss. At the maximum efficiency flow rate, runner loss dominated the total hydraulic loss, accounting for 55.9%, followed by outlet pipe loss at 26.7% and casing loss at 10.4%. Runner friction loss contributed 6.7%, and inlet pipe loss was only 0.3%. This distribution gives the study its interpretive focus: the major performance penalties were not distributed evenly through the machine, nor were they controlled primarily by simple wall friction.

A decisive physical link emerged between runner loss and separation at the blade inlet. Flow entered the runner locally from the suction-surface side of the blade inlet end rather than uniformly across the inlet. This local inflow produced separation on the blade suction surface, and the region of high total pressure loss expanded or contracted with blade phase.  When the separated region extended over much of the suction surface, the runner loss increased. The design consequence is direct: the compatibility between blade inlet angle and turbine-mode inflow controls the scale of inlet separation, and that separation controls much of the runner loss.

The outlet pipe loss had a different origin and the downstream of the runner, the calculations identified a large central vortex and additional vortices generated near the blade outlet on the shroud side and near the runner outlet. Regions of high total pressure loss aligned with these vortex structures. The large central vortex changed relatively little with runner rotation, while the vortices near the blade outlet and runner outlet expanded at phases associated with higher outlet pipe loss. As the blade outlet and runner outlet deliver different vortex structures into the outlet pipe, the local vortex-induced total pressure loss rises or falls.

The findings of Professor Yasuyuki Nishi and Ms. Natsumi Itoh have direct engineering value for the design of non-blocking small hydropower turbines based on reverse-running sewage pumps. The wide passage of a single blade centrifugal pump is attractive where foreign matter can enter the flow, but the study shows that hydraulic performance depends strongly on how the reverse turbine flow interacts with the blade inlet, runner passage, casing, and outlet pipe. One application is runner redesign for turbine-mode operation. Since the dominant runner loss is caused by separation on the suction surface at the blade inlet, the blade inlet angle and blade angle distribution can be reconsidered for reverse-flow operation rather than being evaluated only from the standpoint of pump-mode design. This could guide practical modifications such as reshaping the blade inlet, adjusting the incidence condition, or developing a runner geometry that reduces separated flow while preserving a wide flow path.

A second application is loss-targeted optimization of pump-as-turbine systems. The study quantified the loss distribution, showing that runner loss, outlet pipe loss, and casing loss are the major contributors at the maximum efficiency condition. This helps engineers prioritize design effort. Instead of treating low efficiency as a general problem of single blade machines, the work identifies where improvement is most likely to matter: suppressing blade-inlet separation, reducing vortex-related outlet pipe losses, and moderating phase-dependent casing losses. The findings are also useful for improving operational stability as the head coefficient, output coefficient, efficiency, and hydraulic losses changed with runner rotation, meaning that unsteady performance is inherent to this configuration. Understanding which flow structures cause these variations can support designs with smaller performance fluctuations, reduced vibration risk, and more stable output. For small hydropower installations, this matters because reliability and continuous operation are as important as peak efficiency. In practical terms, the study by Nishi and Itoh provides a hydraulic map for turning single blade sewage-pump geometry into a more viable reverse-running turbine technology.

Ceramic–matrix composite laminates are important structural materials, especially for parts that must work under high temperature and mechanical loading. When these materials contain a circular hole, blunt notch, or similar strain-concentrating feature, their behavior is more complicated than standard elastic stress analysis can explain. In these laminates, they can also activate matrix cracking, interfacial sliding, fiber bridging, fiber fracture, and other inelastic displacements that reshape how load is carried across the net-section.   In a conventional elastic analysis, the highest stress is expected at the notch edge, and design logic follows directly from that peak. For ceramic–matrix composites, however, matrix cracking can relax the stress concentration while increasing local strain concentration. The material is not simply weaker because it contains a notch; its apparent notch sensitivity depends on how damage-mediated inelastic deformation redistributes load. A proper description must therefore capture both the reduction of local stress and the growth of local strain, instead of treating one as a substitute for the other. The difficulty is partly one of scale. The macroscopic stress–strain field depends on damage processes occurring across individual layers, fiber orientations, matrix cracks, and fiber–matrix interfaces. Transverse layers may crack first near the notch, cracks may extend into longitudinal layers, and broken fibers may undergo sliding and displacement governed by interfacial shear. These events are not just microscopic details. They leave measurable signatures that influence, and reflect, the stress state in the concentration region. However, a direct connection between these microstructural features and the macroscopic mechanical response has remained difficult to quantify.

In a recent research paper published in International Journal of Mechanical Sciences, Dr. Xiaoyi Guan and Professor Zhengmao Yang from the Institute of Mechanics at the Chinese Academy of Sciences, working with Dr. Yana Wang and Dr. Jian Jiao from the National Key Laboratory of Advanced Composites at AECC Beijing Institute of Aeronautical Materials, examined how stress and strain redistribute ahead of notches in ceramic–matrix composite laminates and how those macroscopic fields relate to damage features observed by X-ray tomography.  The researchers built their analysis around three linked modeling levels. At the laminate scale, they used an inelastic constitutive model implemented in finite element analysis to describe the nonlinear stress–strain response of the composite. Alongside this, they developed an analytical model based on Neuber’s theory to predict stress and strain concentration behavior associated with matrix-cracking-induced inelasticity. At the microscale, they used a shear-lag-based model to connect fiber break displacement with local distal stress. The finite element model provides a more detailed field calculation, the analytical model offers a faster route to stress–strain prediction, and the micromechanical model gives physical meaning to damage features extracted from tomography.

The material system was a SiC fiber/SiC matrix laminate with a thin BN coating between fiber and matrix, fabricated through a prepreg–melt infiltration route. The laminate had a [0°/90°] arrangement, with specimens containing circular holes across several diameter-to-width ratios. Monotonic tensile tests, supported by two-dimensional digital image correlation, provided load–displacement curves and local strain fields. The use of DIC was important because the central question depended not only on global strength or stiffness, but on the spatial distribution of strain along the net-section near the notch. The load–displacement responses changed systematically with notch size. Specimens with smaller diameter-to-width ratios showed a clearer sequence of elastic deformation, matrix-cracking-related stiffness reduction, and a later fiber-dominated load-bearing stage. Larger holes shortened or suppressed this staged behavior, indicating that the notch geometry changed how much damage tolerance remained before final failure. The finite element model reproduced the elastic response closely and captured the strain fields under representative loading conditions. Its prediction of matrix cracking stress was also close to the experimental value, while the later stages of failure were treated more cautiously because additional damage mechanisms become increasingly important near ultimate failure. The comparison between the analytical model and finite element calculations is one of the more useful parts of the study. For stress concentration, both models gave closely aligned predictions over important ranges, especially before and around the onset of broader inelastic deformation. As inelastic displacements developed near the notch, the stress concentration factor decreased rapidly; around the matrix-cracking stress level, the reduction relative to the elastic state reached roughly 30–35%. The strain behavior moved in the opposite direction. Strain concentration increased when inelasticity occurred, with a peak near the condition where stress relief was strongest. In physical terms, local inelastic strain is the mechanism through which peak stress is moderated.

The spatial stress and strain distributions along the net-section sharpen this interpretation. Near the notch edge, elastic analysis overestimates stress once inelastic deformation begins, because matrix cracking and related deformation relax the local stress. Farther from the notch, however, equilibrium requires load redistribution, and the stress can be higher than an elastic calculation would suggest. This means that the region of greatest design concern cannot be inferred from elastic peak stress alone. For strain, both the finite element and analytical approaches gave predictions that generally lay within the DIC-measured bands, with the analytical model often providing the safer estimate under moderate strain conditions. A further contribution is the elastic–inelastic domain map. By plotting the transition behavior as a function of applied stress and notch size, the researchers separated regions where the analytical model is efficient and adequate from regions where finite element prediction using the inelastic constitutive model is preferable. The same boundary applies to both stress and strain analyses, which makes the map particularly useful as a practical modeling guide rather than a purely descriptive result.

The tomography analysis then connected these macroscopic fields to local damage. In the fractured specimen with a larger notch, crack distributions differed between 0° and 90° layers, with the 90° layers showing more transverse cracking. Fracture fraction generally decreased with distance from the notch edge. Fiber break displacement in the 0° layers followed a similar spatial trend, becoming smaller farther from the notch. When these measured displacements were inserted into the micromechanical relation, the estimated stress distribution correlated strongly with the stresses predicted by both the finite element and analytical models. The design choice to use fiber break displacement rather than crack opening displacement was scientifically consequential, because fiber breaks remained identifiable in the tomographic images even when matrix cracks could appear closed after fracture.

The research work of Professor Zhengmao Yang and colleagues has several engineering applications, most notably in the design, modeling, and damage assessment of ceramic–matrix composite components that contain notches, holes, or other geometric features, particularly in structures exposed to high temperature and mechanical loading. The clearest application is in aero-engine and other high-temperature structural parts made from SiC/SiC ceramic–matrix composites. These components often include holes, cut-outs, joints, cooling passages, attachment points, or similar features that disturb the local stress field. The study gives engineers a more realistic way to understand what happens around these features. Instead of assuming that the highest elastic stress alone controls the response, it shows how stress and strain can redistribute once matrix cracking and other inelastic mechanisms begin. In notched CMC laminates, matrix cracking and local inelastic deformation can reduce the peak stress near the notch while increasing the local strain concentration. That means engineers should not assess these components only by the elastic stress concentration factor. For practical design, both stress redistribution and strain concentration need to be considered when selecting notch sizes, hole diameters, ligament widths, and allowable load levels. The study also supports damage-tolerant design of CMC laminates. Smaller notches may still allow a more gradual damage process, beginning with elastic response, followed by matrix cracking and later fiber-dominated load bearing. Larger notches, by contrast, can shorten this staged response and reduce the remaining damage tolerance. This distinction is important for components that must continue to carry load even after local matrix cracking has started. There is also a clear application in finite element modeling and structural simulation. The hierarchical framework gives engineers a way to decide when a faster analytical model is sufficient and when a more detailed inelastic finite element analysis is needed. The domain map is useful in this respect because it links model choice to applied stress, notch size, and the elastic–inelastic state of the laminate.  By linking fiber break displacement observed through X-ray tomography with the local stress distribution, the study offers a way to interpret damage patterns after loading or fracture. In practical terms, tomography of damaged CMC parts could help engineers identify where stress concentration was most severe and how damage developed around notches.

Artificial structural color is produced when light interacts with deliberately shaped micro- and nano-scale surface features rather than with molecular pigments or coatings. In such systems, optical response is inseparable from geometry. Period, orientation, sectional profile, blaze angle, and height all influence how incident light is redirected, separated, or concentrated. Blazed nano-gratings are an important class of optical surfaces because their asymmetric sawtooth profile can direct diffracted energy more efficiently than a symmetric grating. Their usefulness, however, depends on whether the intended three-dimensional topography can be manufactured with sufficient precision and repeatability. This requirement creates a demanding problem for ultra-precision manufacturing. It is not enough to generate a periodic nanoscale texture; the grating must have a controlled spatial orientation, a prescribed spacing, a defined blaze angle, a consistent height, and a stable profile across the machined area. Traditional approaches such as mechanical ruling and lithographic fabrication can provide high precision in selected geometric parameters, especially grating period and blaze angle, but flexible control of grating orientation and full 3D topography remains difficult. Other methods can introduce orientation variation or nanoscale patterning, yet they may be constrained by topographical consistency, process complexity, cost, or material adaptability.

Vibration-assisted diamond cutting offers a more direct and potentially flexible route for producing micro/nano structures on metal surfaces. Still, many existing vibration-assisted or ruling-based methods form gratings through tool flank interference, local ploughing, or direct replication of tool geometry. These mechanisms can limit the independent control of grating geometry, especially when the desired structure is a sharp blazed nano-grating with a specified 3D form. The issue becomes more demanding on harder metallic materials, where elastic recovery, chip flow, cutting force, and tool wear can alter the final surface from the intended design. In the proposed oblique triangle vibration-assisted diamond cutting process, a single crystalline diamond tool is mounted on a two-degree-of-freedom non-resonant vibration generator, while a rotational stage sets the tool azimuth angle. This changes the relationship between the tool coordinate system and the workpiece coordinate system: the vibration plane remains tied to the rake-face geometry of the tool, but it is no longer coincident with the nominal cutting plane when the tool is rotated. Grating orientation can be set by tool azimuth angle, while spacing depends on equivalent cutting speed and vibration frequency; height and blaze angle can then be tuned through the triangle vibration trajectory.

In a recent research paper published in Journal of Materials Processing Technology, Dr. Changcheng Lin, Dr. Qinghou Cheng, Dr. Yingxue Yao, and Professor Yang Yang from Harbin Institute of Technology (Shenzhen) developed an oblique triangle vibration-assisted diamond cutting process for deterministic fabrication of blazed nano-gratings with controlled 3D topography. They also developed a kinematic surface prediction model that maps tool motion and tool geometry to the final grating topography. The process was experimentally demonstrated on aluminum alloy, brass, 304 stainless steel, and nickel-plated stainless steel with controlled profiles and limited tool wear.

The researchers used this geometric relationship to establish a bidirectional mapping between process parameters and grating topography. They expressed full grating shape through five independent parameters: width, spacing, blaze angle, height, and orientation. Afterwards, they built a numerical prediction model by discretizing the tool and workpiece, mapping the tool motion into the workpiece coordinate system, and updating the machined surface height as the tool swept through the material. The simulation served as the analytical bridge between a desired 3D grating and the process settings needed to produce it. Setting the equivalent blaze angle below the tool clearance angle also had a clear scientific consequence: it reduced flank-face ploughing and helped preserve the intended nanoscale profile rather than allowing elastic recovery to distort the intended final shape.

Machining experiments first examined whether grating orientation and spacing could be separated in practice. At constant nominal cutting speed, changing the tool azimuth angle rotated the grating orientation as expected, while the equivalent grating spacing changed because the effective cutting speed in the grating normal plane changed. When nominal cutting speed was adjusted to compensate for that projection effect, the equivalent spacing remained close to the target value while the orientation varied from negative to positive azimuth angles. SEM observations and simulation results showed close agreement, giving the process model practical weight rather than leaving it as a purely geometrical exercise.

Professor Yang Yang and colleagues then performed AFM measurements to clarify the 3D quality of the produced gratings and found that, at an equivalent grating spacing of 764 nm, the structures retained a regular asymmetric triangular blazed profile over different orientations, with measured grating heights around 115 ± 5 nm and an equivalent blaze angle near 12.5 degrees. At a fixed azimuth angle of −30 degrees, the process produced submicron-period gratings over a range from 300 to 1000 nm with generally good agreement between experiment and simulation. At the smallest tested spacing of 200 nm, edge-radius effects and nanoscale micro-ploughing became more pronounced, which altered profile consistency. The authors used this observation constructively, defining the practical submicron range for the process under their tested conditions as 300 to 1000 nm.

The process dynamics were examined through real-time trajectory measurement, three-directional cutting forces, and chip morphology. The triangle vibration trajectory remained stable during machining, with reported maximum displacement errors of 0.146 μm in the X_φ direction and 0.055 μm in the Z_φ direction. Cutting force analysis showed intermittent cutting behavior, while the main and feed force components changed strongly with tool azimuth angle. At an azimuth angle of 30 degrees, the resultant cutting force was about 25% lower than in orthogonal cutting, and its direction changed with the tool angle, corresponding to altered chip flow. This force-and-chip evidence gave the formation mechanism a physical basis: topography formation was not only a matter of programmed motion, but also of how the oblique tool posture redistributed cutting load and material flow.

The findings of Harbin Institute of Technology scientists have direct relevance to the manufacture of optical functional surfaces where the performance of a component depends on nanoscale geometry rather than on bulk material alone. Blazed nano-gratings are used to direct diffracted light efficiently, and the ability to control their spacing, orientation, blaze angle, height, and width gives engineers a more reliable route for designing surfaces whose optical response is prescribed before machining. This is especially useful for diffractive optical elements, spectrometer components, optical variable devices, and structural-color surfaces, where small deviations in grating profile can alter brightness, hue, angular response, or diffraction efficiency. The paper specifically connects blazed nano-gratings with diffractive waveguide augmented-reality displays and high-performance spectrometers, where precise 3D topography is critical to optical behavior. A major engineering application is in deterministic surface texturing for metals. Many optical micro/nano-structures are fabricated on limited material classes, but this process was demonstrated on aluminum alloy, brass, 304 stainless steel, and nickel-plated stainless steel. That matters because real devices often require different substrates for stiffness, corrosion resistance, reflectivity, durability, or integration with mechanical assemblies. The ability to texture both non-ferrous and ferrous metals with controlled grating geometry broadens the manufacturing window for optical components that must also meet structural or environmental requirements.

The new method also supports design of polarization- and angle-dependent optical surfaces. Because grating orientation can be controlled through the tool azimuth angle, the same machining principle could be used to fabricate surfaces with spatially varied orientations. Such surfaces are relevant to optical variable devices, anti-counterfeiting elements, polarization-sensitive imaging components, and information-encoded textures. The paper’s structural-color experiments show that the machined gratings can generate orientation- and incident-angle-dependent color responses, supporting their use in engineered visual and optical signatures. From a manufacturing perspective, the most important application may be process planning. The bidirectional mapping between process parameters and grating topography allows engineers to move from a target geometry to machining conditions more systematically. This reduces trial-and-error fabrication and makes nano-grating production more compatible with precision manufacturing workflows. The intermittent cutting behavior and reduced tool wear are also practically important, especially when extending diamond cutting to harder metallic substrates where force, heat, and tool degradation often limit scalability.

Phonon transport in low-dimensional carbon materials can be described in particle terms, with phonons drifting, scattering, and reflecting at boundaries. But the same transport can also be treated in wave terms, where coherence, localization, and resonance enter directly into the problem. Scattering-based strategies have been used to suppress thermal conductivity by introducing rough edges, interfaces, pores, defects, and dopants, whereas wave-based strategies have focused on periodic architectures that reshape vibrational spectra and induce resonant disruption of heat-carrying modes. Graphene nanoribbons provide a particularly sensitive platform for this distinction because their thermal transport is strongly shaped by reduced dimensionality, edge morphology, and geometric confinement, so even small structural modifications can alter transport through physically distinct mechanisms. Pillared graphene nanoribbons introduce an additional level of complexity. The pillars can act as scattering sites that shorten phonon mean free paths, yet they can also support localized vibrational modes that hybridize with propagating modes in the ribbon body. Once that happens, the transport problem is no longer just one of extra boundary scattering. It becomes a question of how much of the conductivity reduction comes from particle-like interruption of trajectories and how much comes from wave-mediated modification of the phonon spectrum itself.

Earlier work on related nanostructures had already shown that phonon wave effects can be important in periodic or resonant systems, and prior studies on pillared graphene nanoribbons had largely concentrated on local resonant hybridization. That emphasis still leaves an interpretive gap in the analysis. In graphene-based nanostructures, intrinsic phonon mean free paths are long and transport is unusually sensitive to boundaries. A pillared architecture therefore invites two intertwined explanations for reduced conductivity, one rooted in enhanced boundary scattering and the other in resonance-driven wave effects. In a recent research paper published in International Journal of Thermal Sciences, Dr. Shixian Liu, Dr. Fei Yin and Professor V.I. Khvesyuk from the Bauman Moscow State Technical University working together with Dr. Zhicheng Zong and Professor Nuo Yang from the National University of Defense Technology, developed a quantitative method for separating phonon particle and wave contributions to thermal conductivity reduction in pillared graphene nanoribbons. They combined calibrated Monte Carlo simulations with molecular dynamics and defined wave and particle ratios from the difference between the two approaches. They also introduced the concept of resonance hybridization depth to describe how far pillar-induced hybridization extends into the ribbon and showed that this depth grows with pillar height.

The researchers approached the problem by pairing two simulation frameworks that encode different physical emphases. Monte Carlo calculations were used in a phonon Boltzmann transport setting to capture particle-like behavior, especially phonon-phonon scattering and phonon-boundary scattering. Molecular dynamics, by contrast, retained the atomic vibrational picture and therefore can naturally include both scattering and wave-related effects such as resonance hybridization and localization. This comparison turns the difference between the two methods into a physically interpretable quantity rather than treating it as a methodological inconvenience. Before moving to pillared systems, they calibrated the Monte Carlo boundary treatment so that ordinary graphene nanoribbons without pillars gave thermal conductivities consistent with molecular dynamics. That step matters because it allowed later discrepancies in the pillared structures to be interpreted as wave effects rather than as artifacts of mismatched methods.

Once the benchmark was established, the pillared graphene nanoribbons showed a clear separation between the two transport contributions. In every case, Monte Carlo predicted a reduction in thermal conductivity because the pillars introduced additional boundary scattering and shortened the phonon mean free path. Molecular dynamics predicted a stronger reduction, indicating that scattering alone was not the whole story. The authors defined relative thermal conductivities for pillared structures and then extracted a wave ratio from the difference between Monte Carlo and molecular dynamics results. That formulation let them decompose the total conductivity reduction into particle and wave contributions in a way that remained directly tied to computed transport data.

The collaborative team reported that particle effects dominated the suppression of thermal conductivity in these graphene-based structures, while wave effects made a smaller but still measurable contribution. This is physically important because it distinguishes pillared graphene nanoribbons from systems such as pillared silicon nanowires, where wave effects can be stronger. The difference, as the paper argues, is tied to graphene’s long phonon mean free path and strong sensitivity to boundary modification, along with the reduced set of available phonon modes in a quasi-one-dimensional ribbon confined within a two-dimensional plane. In other words, the same architectural idea does not produce the same transport balance across materials or geometries. Here, the pillars act first as strong scatterers and only second as resonant wave modifiers.

Geometry then reshaped that balance in revealing ways. Increasing pillar height reduced thermal conductivity further. The scientific consequence is twofold: taller pillars increase the fraction of phonons that undergo boundary scattering, and they also introduce more resonant modes, including lower-frequency modes that hybridize more effectively with heat-carrying vibrations. Width behaved less simply. As ribbon width decreased, the wave contribution first increased and then declined, which the authors interpreted as a saturation effect. Once the wave contribution saturates, continued narrowing mainly strengthens particle-related suppression through boundary and pillar scattering. Temperature added another layer. Rising temperature slightly reduced the wave ratio, consistent with shorter phonon wavelengths and more frequent scattering events that weaken the coherence needed for wave-like transport.

The study also goes beyond quantifying of conductivity reduction and introduces a structural descriptor for resonance influence. By analyzing dispersion relations, group velocity reduction, and the frequency-integrated product of group velocity and density of states, the researchers extracted a resonance hybridization depth from the width-dependent decay of this spectral transport indicator. This quantity represents the spatial extent over which local pillar resonances significantly modify phonon transport in the ribbon. Its value increased with pillar height, which indicates that taller pillars do not just create stronger local resonances but broaden the region over which those resonances alter transport. That idea gives the analysis a clear conceptual endpoint: wave effects are not only present or absent, they occupy a quantifiable spatial domain within the nanoribbon architecture. Nanostructured thermal transport is full of architectures described as resonant, coherent, metamaterial-like, or scattering-dominated, but many of those labels remain loosely assigned when different mechanisms lead to the same macroscopic outcome. The authors show that in pillared graphene nanoribbons, reduced thermal conductivity cannot be read as automatic evidence of wave-dominated transport. The pillars do generate resonance effects, and those effects are quantifiable, but the larger share of suppression comes from particle-like boundary scattering. That point matters because it changes how one interprets structural design in graphene-based thermal materials.

There is also a methodological advancement that reaches past this particular geometry. Indeed, the comparative use of Monte Carlo and molecular dynamics gives a practical route for separating transport contributions that are otherwise entangled. It is a careful strategy because it does not ask either method to do what it cannot naturally do. The particle-based formalism emphasizes scattering-driven transport; the atomistic formalism can retain wave -related behavior; the comparison between them becomes the source of physical insight. This reasoning also suggests a useful template for other nanoscale systems in which dual phonon character is difficult to isolate. Moreover, the authors’ introduction of resonance hybridization depth adds another useful layer. It turns a fairly abstract notion, the spatial reach of resonant phonon modification, into something that can be estimated and compared across structures. Within the paper’s own scope, that supports a more concrete design logic for tuning width and pillar height depending on whether one wants stronger scattering, stronger resonance hybridization, or a particular balance between the two. The broader implication is not that all pillared nanostructures behave the same way, but that their particle-like and wave-like effects can be separated, quantified, and connected to structural design.

Air-breathing propulsion at hypersonic speed depends on the ability to sustain combustion within short aerodynamic timescale. In a scramjet combustor, the incoming flow remains supersonic, and the fuel–air mixture must ignite, release heat, and remain anchored before it is convected out of the combustor. This creates a demanding coupling between fluid residence time, chemical ignition delay, turbulent mixing, and flame stabilization. Even when the inflow temperature is high enough to make reaction possible, stable combustion is not guaranteed, because the chemical induction time may still exceed the residence time available within the cavity shear layer. For this reason, forced ignition and flame-holding strategies have become main issues in scramjet combustion research. Cavity flame holders are often used because the recirculation region and shear layer can increase the effective residence time of the reacting mixture and provide a region where heat and radicals can support flame stabilization. However, cavity-based stabilization is not a purely geometric problem. The ignition process depends strongly on the thermochemical state of the mixture within the shear layer, where fuel, main airflow, recirculated products, and torch gases interact. A cavity that is aerodynamically suitable still needs to provide a shear-layer residence time compatible with the chemical induction time. In a recent research paper published in Journal of Engineering for Gas Turbines and Power Professor Shinichiro Ogawa from the Department of Aerospace and Marine-System Engineering at Osaka Metropolitan University, developed a numerical evaluation method that combines the forced ignition model in the shear layer with a plug flow reactor calculation for cavity-held scramjet ignition. The technically distinct feature is that the model links local shear-layer residence time, torch-modified gas composition, and detailed hydrocarbon chemistry through a Damköhler-number ignition criterion. It was then used to separate forced-ignition behavior under H₂O-vitiated and nonvitiated airflow conditions. The approach also enabled direct comparison of Type A and Type B cavity geometries under varied micro-rocket torch input energies. Instead of treating the combustor as a fully resolved turbulent reacting flow, the analysis extracted the shear-layer velocity, temperature, residence time, and gas composition through the forced ignition model and then used a one-dimensional chemical reaction analysis to determine whether the ignition delay was short enough for ignition to occur within the available residence time.

When the residence time in the shear layer exceeded the calculated ignition delay, the condition corresponded to forced ignition; when the ignition delay was longer, ignition failed within that region. This is a useful reduction because it ties the chemistry directly to the cavity’s flow time rather than treating ignition as a temperature threshold alone. The comparison with prior combustion experiments gave the model its practical anchor: cases that ignited experimentally fell into the region where the calculated Damköhler number exceeded unity, while nonignition cases corresponded to values at or below unity. The validation covered two cavity geometries and several torch operating conditions, including cases with different helium dilution levels in the torch gas.

A methodological comparison between two plug-flow implementations helped establish the calculation route. The Lagrangian particle simulation and the chain-of-reactors approach produced nearly the same ignition-delay behavior for representative Type A and Type B cavity cases, with only small differences near the temperature-rise front. Because the chain-of-reactors method required lower computational cost and gave stable solutions during rapid reaction progress, it was used for the forced-ignition limit calculations and by enabling multiple parameter variations, the chain reactor strategy allowed the study to map how torch energy, cavity residence time, and H₂O vitiation interact.

Professor Shinichiro Ogawa used vitiation analysis and compared shear-layer ignition with and without water vapor contamination under otherwise matched Type A cavity conditions. In the vitiated case, the onset of temperature and CO₂ changes occurred farther upstream than in the nonvitiated case. The temperature-rise gradient was initially gentler with H₂O present, but the reacting flow eventually reached comparable high-temperature levels near the cavity ramp. The CO₂ behavior sharpened the interpretation: H₂O vitiation increased CO₂ production and shifted its onset, which indicates a change in oxidation progress rather than a simple passive dilution effect. The author performed species-based analysis which gave the chemical picture more texture. He found that fuel consumption of methane and ethylene became stronger downstream of the early shear-layer region under H₂O vitiation, while OH production increased markedly over the region where ignition developed. Formaldehyde and ketene formation also rose, consistent with activation of oxygen-containing intermediate pathways. The presence of H₂O therefore did not just delay or weaken the reaction through heat-capacity effects. Under the modeled conditions, it also promoted radical and intermediate formation in ways that supported oxidation and moved the ignition process upstream.

Torch energy also changed the sensitivity of the ignition process. As the net input energy increased, the forced-ignition limit temperature fell for both cavity types, indicating that stronger torch input partly compensated for conditions in which ignition chemistry required a longer induction time. This effect was clearest when the net input energy increased from 20 to 25 kW, where the required airflow temperature dropped markedly in both cavity configurations. The influence of H₂O was strongest in the lower-to-moderate torch-energy range, where chemical delay still governed whether ignition could occur within the shear layer residence time. At higher torch energy, the thermal and radical content supplied by the torch became more dominant, reducing the relative importance of vitiation. The cavity comparison clarifies why residence time matters. Type A, with its shorter residence time, showed greater sensitivity to H₂O vitiation because a small change in ignition delay had a larger consequence when the available flow time was limited. Type B, with a longer shear-layer residence time, was less affected. The design choice of cavity length therefore had a direct scientific consequence: shorter residence time amplified the chemical influence of vitiated-air composition, whereas longer residence time reduced the dependence of forced ignition on the presence of H₂O.

The findings of Professor Shinichiro Ogawa have direct engineering relevance for the design and interpretation of scramjet combustor ignition systems, especially when ground-test data are used to support flight-oriented combustor development. In many high-enthalpy ground tests, the main airflow contains H₂O from the combustion heater, while flight air would not contain the same level of water vapor. Ogawa’s analysis shows that this difference can shift the forced-ignition behavior by changing radical chemistry, ignition delay, and the location where reaction begins in the cavity shear layer. For engineers, the implication is that ignition limits measured in vitiated facilities need to be interpreted with explicit attention to H₂O effects before being used for flight-oriented assessment. A combustor that appears to ignite reliably in a ground facility may require different torch energy, airflow temperature, or cavity residence time under nonvitiated conditions.

For cavity design, the practical message from the study is that residence time controls how strongly facility chemistry enters the ignition margin. Compact cavities require greater attention to torch input energy and facility-to-flight chemical differences, whereas longer-residence-time configurations reduce that dependence. If a compact cavity is preferred for aerodynamic or structural reasons, the ignition system must be designed with greater attention to torch input energy and facility-to-flight chemical differences. If the design allows a longer residence-time cavity, the combustor may be less dependent on chemically favorable vitiated-air conditions. Increasing the net input energy lowered the forced-ignition limit temperature and reduced the relative influence of H₂O vitiation. This suggests that torch power can be used as a design parameter to preserve ignition margin under conditions with longer ignition delay, especially in short-residence-time configurations. However, the paper also notes that higher torch energy increases thermal loading on the torch body, especially during prolonged uncooled operation. Therefore, the results can help define a practical operating envelope: enough torch energy to secure ignition, but not so much that durability or permissible operating time becomes limiting.

A further application of Professor Shinichiro Ogawa study is in reduced-order combustor modeling and by linking shear-layer residence time, ignition delay, torch-gas composition, and Damköhler number, the model provides a useful engineering tool for screening ignition limits before more expensive full reacting-flow simulations or combustion tests are performed. This is particularly valuable during early-stage design, where many combinations of cavity geometry, torch condition, and inflow temperature must be screened before full reacting-flow simulations or combustion tests. The paper therefore improved our understanding of H₂O vitiation and also in developing a more rational workflow for designing and interpreting forced-ignition systems in scramjet combustors.

Industrial robots are indispensable in manufacturing tasks where programmed motion must be repeated with speed, stability, and geometric consistency. However, in precision applications, repeatability alone is not sufficient and a robot may return to nearly the same position again and again, while still failing to reach the exact pose expected from its nominal kinematic model. This difference between repeatable motion and absolute positioning accuracy is especially important in operations such as machining, assembly, inspection, and other tasks where small spatial errors can accumulate into measurable process deviations.  A considerable part of the positioning error in a serial industrial robot originates from kinematic deviations in the links, joints, assembly relationships, and transmission system. These errors may be introduced during manufacturing and installation, but they can also change gradually as the robot continues to operate. Conventional kinematic calibration addresses this problem by measuring robot poses, identifying deviations in the model parameters, and compensating the command or model accordingly. When accurate measuring instruments are used, such as a laser tracker, the absolute positioning performance can be improved substantially. The difficulty, however, is that this type of calibration is not naturally compatible with continuous industrial operation. It often requires specialized equipment, manual intervention, and measurement configurations that are separate from the robot’s normal work trajectory. In many cases, calibration also occupies the end-effector interface or physically restricts the robot’s motion, which means that the robot must stop working while calibration is performed.

This interruption is not a minor practical inconvenience. It becomes a deeper engineering limitation when calibration has to be repeated after long-term operation, especially because joint transmission wear can gradually reduce positioning accuracy again. A method that depends on repeated use of expensive metrology equipment and work suspension can be accurate, but its practical use becomes more difficult in manufacturing environments where uptime and process continuity matter.  In a recent research paper published in Precision Engineering, Professor Zhouxiang Jiang, Dr. Ruoheng Ding, Dr. Yuxuan Liu, Dr. Zhongjie Long from the Beijing Information Science & Technology University working together with Professor Bao Song from Huazhong University of Science & Technology developed a real-time kinematic calibration method that uses cameras and visual markers mounted on different robot joints instead of relying continuously on a laser tracker. They introduced two dimension-reduced kinematic-error models that split the parameter identification problem into early-joint and later-joint components under work-trajectory and marker-visibility constraints. They also designed configuration-specific backpropagation neural networks to convert inaccurate marker poses measured by cameras into accurate joint pose estimates. A further technical feature is the training-data strategy that combines workspace and jointspace regularity to improve prediction accuracy in the experimental robot.

The researchers chose measurement configurations as in a traditional full-workspace calibration. Joints 1–3 were constrained by the work trajectory because they define the major geometry of the robot motion and cannot be displaced arbitrarily without changing the task path. Joints 4–6 retained broader freedom because they chiefly determine orientation. Marker visibility added another practical filter: a configuration remained useful only if the marker plane could be seen by the camera within an acceptable angular range. This design choice directly shaped the scientific consequence of the method: the calibration data became compatible with robot work, but the model had to be restructured so that identifiability would not collapse under the restricted pose set.

The authors used simulation to separate the effects of model structure and pose prediction before moving to the experimental robot. They compared the dimension-reduced models with a conventional model under the same visibility and trajectory constraints, and also with a conventional full-workspace model. The reduced models produced identification and calibration behavior much closer to the unconstrained traditional case than to the constrained high-dimensional case. Along the work trajectory, the reduced-model strategy even gave smaller residuals than the traditional full-workspace model in the simulated comparison, which is consistent with the idea that calibration is most effective in the region where measurement configurations are generated. The neural-network component addressed the measurement side of the problem. The authors modeled systematic vision errors arising from lens distortion, camera calibration error, and illumination-dependent marker appearance, then used backpropagation neural networks to learn mappings from measured marker poses to joint poses. In simulation, the trained networks predicted joint poses with very small errors, and the authors observed better accuracy near the center of the testing sample, where actual pose errors would normally be concentrated if the robot drift remains modest.

The team used for the experimental validation a six-degree-of-freedom robot, first calibrated with a laser tracker to establish the baseline parameters. A camera-marker system then supported real-time calibration along a designed work trajectory. One important experimental detail is the refinement of training data for joint 6. Workspace-regular training alone did not give sufficiently accurate orientation prediction, because the corresponding joint angles were irregularly distributed. Adding data with regularity in joint space and combining it with the workspace-based samples improved the learned mapping. With these networks and the reduced models, the method gave lower maximum residuals than the constrained conventional model both in part of the visible workspace and along the work trajectory. Along the trajectory, the reported maximum residual for the reduced-model approach was 0.235 mm, compared with 0.326 mm for the constrained conventional model and 0.288 mm for the conventional full-workspace model. In the visible workspace comparison, the reduced-model approach was also close to the unconstrained traditional calibration, with maximum residuals of 0.375 mm and 0.398 mm, respectively.

The findings of Professor Zhouxiang Jiang  and colleagues have direct relevance for manufacturing environments where industrial robots are expected to maintain high absolute positioning accuracy without repeated interruption for conventional calibration. In robotic machining, drilling, trimming, grinding, polishing, and precision assembly, even a robot with good repeatability can gradually lose geometric accuracy because of kinematic parameter deviations and joint transmission wear. The method developed by Jiang, Ding, Liu, Long, and Song offers a route to monitor and restore positioning accuracy while keeping the calibration process closely tied to the robot’s actual work trajectory rather than a separate metrology routine. That is especially useful for production cells where stopping the robot, installing laser tracker targets, removing tools, or manually adjusting measurement hardware would reduce throughput and increase operating cost.

A practical application is real-time or near-real-time accuracy maintenance in robotic workstations. By attaching small visual markers to selected joints and using cameras positioned around the workspace, a robot cell could track changes in joint pose during operation and detect when positioning error exceeds an acceptable threshold. The original trajectory could then be slightly adjusted to pass through selected measurement configurations, allowing calibration data to be collected without a full shutdown. This is important for long production runs, where accuracy degradation may not occur suddenly but accumulates gradually as the robot continues working. The dimension-reduced calibration strategy is also valuable for constrained industrial tasks. Many robots cannot move freely through an ideal calibration workspace once they are installed near fixtures, machine tools, workpieces, guarding, or other equipment. The paper shows that calibration can be redesigned around trajectory and visibility constraints, rather than treating those constraints as obstacles. For engineering practice, this means calibration can be localized to the region where the robot actually works, which is often more relevant than improving accuracy uniformly across the full theoretical workspace. Another application lies in lower-cost robot deployment. Laser trackers and similar instruments remain highly accurate, but they are expensive and not always practical for frequent recalibration. A vision-based system trained to map marker measurements to joint poses could reduce dependence on repeated use of high-end metrology equipment after the first calibration. This would be useful for small and medium manufacturers, flexible production lines, and robotic cells that require periodic accuracy recovery but cannot justify frequent manual metrology intervention. The findings suggest a practical path toward robot workcells that can maintain accuracy without repeatedly stopping for conventional metrology-based calibration. By linking camera measurement, learned pose correction, and reduced kinematic modeling, the method brings calibration closer to the conditions under which the robot actually operates.

Prefabricated concrete construction places unusual demands on its connection systems. Unlike cast-in-place members, where reinforcement continuity is formed directly within a continuous concrete body, prefabricated structures achieve structural continuity through connection regions that transfer force between separately manufactured components. The fully-grouted sleeve connection is one of the principal solutions for this purpose and its performance depends on the coordinated action of reinforcing steel, grout, sleeve confinement, and interfacial bond, so its mechanical response cannot be understood only as the tensile behaviour of a steel bar or the compressive behaviour of grout.  Under earthquake–fire coupling conditions, the problem becomes more demanding: cyclic loading introduces repeated tension and compression, while thermal exposure changes the mechanical properties of the grout and steel and may alter the bond conditions within the sleeve. Under these conditions, conventional external measurements, such as axial load, displacement, and final failure mode, remain essential, but they do not fully describe how tensile damage, compressive damage, rebar ductile damage, and interfacial force transfer evolve inside the sleeve before fracture or pull-out occurs.

Previous studies have provided important experimental knowledge on fully-grouted sleeve connections under monotonic tension, cyclic loading, and elevated-temperature or post-fire conditions. They have clarified major failure modes, including steel bar fracture, pull-out, and bending, and have shown that thermal exposure can affect ultimate strength, displacement capacity, and bond-slip behaviour. However, few challenges still exist: first, a reliable quantitative indicator for damage in fully-grouted sleeve connections under mixed cyclic loading–thermal coupling is still not well established and secondly, the internal stress mechanism has often been treated mainly through the axial force state, leaving the radial force and third-axis bending response less fully explained. To address these gaps, a recent research paper published in Engineering Fracture Mechanics, Dr. Yitong Wang and Professor Guoxin Wang from Dalian University of Technology working with Associate Professor Fujian Yang from Changzhou University developed an optimized finite element model for fully-grouted sleeve connections under mixed-mode cyclic loading–thermal coupling, including both force-controlled and displacement-controlled stages. They introduced a quantitative damage pathway for post-fire grout under tensile–compressive cyclic loading and combined it with ductile damage evaluation of the rebar. They also developed a triaxial force interpretation using radial force, axial force, and third-axis bending moment, together with ratio-based measures linking bending response to radial and axial force components.

The researchers built their analysis around an optimized ABAQUS model of the fully-grouted sleeve connection, using a two-dimensional axisymmetric representation to reflect the geometry while keeping the calculation efficient. The model incorporated post-fire constitutive relations for grout and HRB400 rebar, tensile and compressive damage parameters for the grout, and ductile damage for the rebar. Two contact interfaces were central to the calculation: the rebar–grout surface and the grout–sleeve surface.  They validated against reference post-fire cyclic loading experiments on grouted sleeve specimens heated to target temperatures and then loaded under high-stress cyclic conditions. The authors found the model reproduced the main load–displacement behaviour and captured the observed distinction between rebar fracture outside the sleeve and bond-slip pull-out. It was also checked against time–strain behaviour at the sleeve midpoint and against large-displacement cyclic loading data at room temperature and 600 °C. These comparisons gave the later damage analysis a stronger basis, because the model was not used only as a qualitative visualization tool; it was tied back to measurable cyclic response. Afterward, the damage maps gave a more detailed interpretation of failure than the experiments alone could provide. Both tensile and compressive grout damage were concentrated more severely at the rebar–grout interface than at the grout–sleeve interface, and damage decreased from the inner surface toward the outer surface.   The rebar–grout interface is where force transfer is most directly imposed, while the sleeve modifies the stress field through confinement rather than acting as the primary bond surface. When rebar fracture occurred, grout compressive damage was more severe than in bond-slip failure, which the authors linked to the role of mechanical interlock. Once bond-slip develops, the interlocking contribution is reduced, and the associated grout damage is correspondingly lower.

The team also performed rebar damage analysis which added another useful distinction. Inside the sleeve, the reinforcement showed almost no ductile damage, confirming the protective influence of the sleeve over that embedded length. Damage began near the sleeve end and increased toward the exposed rebar end. Under bond-slip failure, the ductile damage outside the sleeve was more uniform, rather than sharply developing into the fracture pattern and this is an important separation between failure appearance and internal damage state. The authors did performance degradation analysis which showed that yield behaviour was comparatively insensitive to grout strength and thermal damage within the investigated range. Yield displacement varied only modestly, and yield force decreased slightly with temperature. Ultimate behaviour was more temperature-sensitive. Ultimate displacement remained relatively stable up to around 400 °C and then declined between 400 and 600 °C, with different reduction levels for the two grouts. Ultimate force also stayed broadly stable within 400 °C, followed by a decrease above that level. The choice to combine simulation data with existing high-stress cyclic experimental results allowed the authors to treat 400 °C not as an isolated observation from one specimen set, but as a practical turning point in the post-fire cyclic performance pattern.

The force-mechanism analysis moved beyond the usual axial interpretation and by extracting radial force, axial force, and the third-axis bending moment, the authors showed that axial loading induces a coupled triaxial response through sleeve restraint and grout-mediated force transfer. High-stress and large-displacement cyclic schemes produced different histories of damage accumulation. Although their final tensile stages were comparable, the earlier cyclic phase changed the internal condition of the connection. High-stress cyclic loading produced greater accumulated damage and lower capacity than large-displacement cyclic loading, while the radial force and bending moment followed closely related trends. The analytical strategy therefore linked the loading scheme directly to a mechanical consequence: the cyclic path altered triaxial force development and damage accumulation, not only the final axial capacity.

The engineering implications of the research work reported by Yitong Wang, Guoxin Wang, Fujian Yang are mainly in post-fire safety assessment, seismic qualification, and damage-informed design of prefabricated concrete connections. Fully-grouted sleeve connections are therefore critical local regions in precast construction.  A connection may still carry load in the early part of testing, while damage has already developed internally at the rebar–grout interface, where bond transfer, confinement, and local cracking interact most strongly. The new study is vital because it gives engineers a more resolved way to examine that hidden damage process. We think the most immediate application is post-fire evaluation of prefabricated concrete structures. After a fire, an engineer must decide whether a precast column, wall, beam, or joint region can remain in service, whether it requires strengthening, or whether the connection should be treated as unsafe. The study shows that yield behaviour may not change dramatically across the investigated thermal range, but ultimate performance becomes more vulnerable once the temperature exceeds about 400 °C.  A connection that appears acceptable at lower load levels may still have reduced deformation capacity, weaker bond resistance, or greater susceptibility to pull-out under later cyclic demand. In structural safety decisions, residual ductility and ultimate resistance are often just as important as initial stiffness or yield strength. The study also has implications for seismic qualification of precast concrete connections. Earthquake-resistant design depends on how a connection behaves under repeated loading, not only on its monotonic tensile capacity. By comparing high-stress and large-displacement cyclic schemes, the work shows that the loading history influences damage accumulation and triaxial force development. This can help researchers and engineers interpret cyclic qualification tests more carefully, especially when post-fire damage is involved.

A further application is connection detailing and repair design. The finding that grout damage is more severe at the rebar–grout surface than at the grout–sleeve surface points to the inner bond-transfer region as a critical zone. This could guide future improvements in grout formulation, sleeve geometry, anchorage design, inspection priorities, or strengthening strategies. By quantifying grout damage, rebar ductile damage, radial force, axial force, and third-axis bending moment, the model provides a basis for locating high-risk regions inside the sleeve and for understanding why axial response alone does not fully represent the mechanical state of the connection.

Piezoelectric shells are coupled electromechanical structures in which mechanical deformation, electric polarization, and electric potential must be considered together during vibration analysis. Such coupling manifests particularly complex behaviors for spherical configurations. The curvature of the shell, the radial direction of polarization, and the vibration modes of the spherical structure all influence how mechanical deformation and electric response develop together. The shell motion drives fluid flow, and in turn the surrounding fluid modifies the natural frequencies and damping characteristics of the shell. For an inviscid compressible fluid, the main effect is associated with fluid inertia and pressure coupling, whereas for a viscous fluid, energy dissipation dominates the fluid-structure interaction problem. When the enclosed fluid is non-Newtonian, viscous dissipation and elastic relaxation can act together, changing both the vibration frequency and the rate of energy loss. This distinction is important especially for small spherical shells and high-frequency vibration modes because, under such conditions, the time scale of structural oscillation may approach the relaxation time of the fluid, so a fluid that appears Newtonian at low frequency may show measurable viscoelastic behavior during vibration. A complete description must therefore account for shear relaxation, compressional relaxation, fluid compressibility, and the three-dimensional deformation of the piezoelectric shell. A formulation limited to radial motion, uncoupled elasticity, or purely Newtonian fluid behavior would describe only part of the coupled dynamics.  The problem was technically demanding because the piezoelectric shell supported both torsional and spheroidal vibration modes, and the enclosed fluid involved dilatational and equivoluminal motions whose behavior depended on shear and compressional relaxation. These two fields had to be solved in a consistent way and then coupled through interface conditions that enforced compatibility of motion and stress at the inner shell surface.

In a recently published research paper in the Journal of Sound and Vibration, Dr. Yuze Cao, Professor Bin Wu, and Professor Weiqiu Chen from Zhejiang University developed a three-dimensional analytical formulation for the free vibrations of a spherically isotropic piezoelectric shell filled with a compressible non-Newtonian fluid. The new model was distinct because it included both shear and compressional relaxation effects in the enclosed fluid and still retained the coupled electroelastic response of the shell. They separated torsional and spheroidal vibration modes, solved the piezoelectric shell equations using displacement functions and the generalized Frobenius power series method, and represented the fluid motion through velocity potentials. The final characteristic equations, successfully solved with the Müller iteration algorithm, provided complex vibration frequencies and quality factors for different modes, shell sizes, fluid viscosities, and fluid relaxation conditions.

Briefly, the researchers first formulated the piezoelectric shell and the enclosed fluid separately, then coupled them through interface conditions and used the resulting frequency equations to examine the effects of material and geometric parameters. For the piezoelectric shell, they began from linear piezoelectricity in spherical coordinates, with the radial direction aligned with the polarization axis. By introducing three displacement functions, they managed to separate the vibration problem into two independent classes. The first class described torsional modes and reduced to an uncoupled second-order differential equation, whereas the second class described spheroidal modes and remained coupled through radial displacement, tangential displacement, and electric potential. This separation allowed the torsional vibration to be handled independently, while the more complicated spheroidal motion was solved through a matrix form of the Frobenius power series method. In practical terms, the analytical strategy converted the three-dimensional electroelastic field into radial functions associated with spherical harmonics, so that each angular mode could be studied through ordinary differential equations rather than through the direct solution of the original three-dimensional field problem for each mode. The breathing mode received separate treatment because its displacement was purely radial and the tangential component did not contribute to the piezoelastic field.

Cao, Wu, and Chen used a compressible linear viscoelastic fluid model that included both the deviatoric shear relaxation effect and the spherical compressional relaxation effect. Small-amplitude harmonic motion permitted linearization, and the generalized Navier–Stokes equations were solved by introducing velocity potential functions. The Helmholtz decomposition separated dilatational and equivoluminal contributions to the fluid motion, and the finite-velocity condition at the center of the filled shell excluded singular solutions. The coupled frequency equations arose when the shell and fluid fields were matched at the inner interface. Continuity of displacement and stress linked the radial and tangential motion of the piezoelectric shell with the fluid velocity and fluid stress, while electrically open-circuited boundary conditions completed the electroelastic problem. They first validated the approach against available theoretical predictions for a PZT-4 shell filled with a non-viscous compressible fluid, where the natural frequencies were real because the fluid model contained no damping mechanism. The close agreement with prior exact results supported the analytical formulation and the numerical root-finding procedure.

The authors presented numerical examples and showed how the model behaved when dissipation and relaxation were present. For torsional modes, fluid viscoelasticity had little influence on the vibration frequency. Spheroidal modes responded more strongly: fluid-induced added mass effects and fluid viscosity reduced the vibration frequency, while fluid viscoelasticity could increase it by introducing an energy-storage contribution. The quality factor followed a more mode-dependent pattern. In many cases, the viscoelastic fluid model predicted higher quality factors than a purely viscous fluid model, but some lower-frequency torsional modes showed the opposite trend because the elastic relaxation contribution remained weak while dissipation was still enhanced.

The breathing mode was especially sensitive to the compressional relaxation effect of the fluid. They found that when the shell vibrated radially, the enclosed fluid underwent volumetric deformation, so the predicted quality factor of a model that neglected compressional relaxation could differ significantly from that of the compressible non-Newtonian fluid model. By contrast, for non-breathing spheroidal modes, shear-associated viscous dissipation dominated attenuation, and the distinction between the Maxwell-type treatment and the compressible non-Newtonian model became smaller. Shell size and vibration order sharpened these effects. As the radius decreased or the radial mode order increased, the vibration frequency rose, making the fluid relaxation time more relevant to the coupled dynamics. In glycerol–water mixtures, the Newtonian model predicted a steady decline in frequency and quality factor with increasing glycerol fraction, but the non-Newtonian model could reverse this trend beyond a critical concentration because elastic energy storage began to compete effectively with viscous loss.

The findings of Cao, Wu, and Chen are directly relevant to the design of piezoelectric spherical containers, resonators, and sensing elements that operate with enclosed complex fluids. In such systems, the fluid cannot always be treated as a simple added mass or as a purely viscous damping medium. The study showed that fluid viscosity, viscoelastic relaxation, shell radius, and vibration mode could each change the natural frequency and quality factor of the coupled system. This is important for engineers who need to predict resonance accurately, especially when Newtonian assumptions may not capture the full frequency and damping response. One practical application is in piezoelectric resonators used for fluid characterization. Because the vibration response depends on viscosity and relaxation time, a spherical piezoelectric shell could, in principle, be used to infer the properties of non-Newtonian liquids from shifts in frequency and quality factor. The breathing mode is especially relevant when the compressional relaxation effect is important, since radial vibration produces volumetric deformation of the enclosed fluid. Other spheroidal modes may be more sensitive to shear-related dissipation. This modal selectivity gives designers a way to choose vibration modes according to the fluid properties they want to probe.

The results are also useful for small-scale acoustic or electromechanical devices containing polymer solutions, glycerol-based mixtures, biological fluids, or other viscoelastic media. At small shell radii or higher vibration orders, the vibration frequency increases, making relaxation effects more visible. The study therefore helps identify when a Newtonian model may be acceptable and when a non-Newtonian description is required. This distinction matters in miniaturized resonators, precision detectors, and fluid-filled piezoelectric components where frequency drift or a loss of quality factor can affect device performance. Another engineering implication is damping control. The analysis showed that viscoelasticity did not simply increase dissipation; under some conditions, elastic energy storage in the viscoelasticity fluid could increase the quality factor relative to a purely viscous model. This means that the choice of enclosed fluid could be used deliberately to tune the dynamic response of a piezoelectric shell. For resonant sensing, a higher quality factor may improve frequency resolution, whereas stronger damping may be useful where vibration suppression is desired.

Free-form surfaces are important in precision manufacturing because their functional value often depends on geometry that is difficult to describe, machine, and verify using simple planar or rotational assumptions. In aerospace components, mold processing, and related high-value manufacturing contexts, surface form is actually part of the dimensional information needed for machining program generation and for later verification of surface quality. Once such parts enter a machining workflow, metrology must provide geometric information with sufficient accuracy, but also with enough practical robustness to operate under realistic surface and environmental conditions. Non-contact measurement methods provide dense and rapid surface acquisition in many cases, however their usefulness can be limited by material properties and measurement conditions. Highly reflective metal workpieces, blackened workpieces after heat treatment, and non-unloaded workpieces requiring in-site measurement can all create difficulties for optical methods. Chromatic confocal probes address some of these limitations, but their sensitivity to surface slope remains a concern for high-gradient free-form geometries. For this reason, contact measurement, particularly coordinate measuring machine and contact-probe-based on-machine measurement, remains central in free-form surface machining metrology. Its strength lies in robustness and measurement accuracy, but its efficiency is tied directly to the number and placement of probing points. In a recent research paper published in Precision Engineering, Dr. Ma Zhifu, Professor Lu Yong, Dr. Ma Shoudong, and Professor Deng Kenan from the Harbin Institute of Technology developed an adaptive sampling method in which a dense triangular mesh of a known free-form surface is simplified and the remaining vertices are used directly as contact measurement points. The method modifies edge-contraction mesh simplification by mapping contraction targets to existing vertices, defining boundary and corner-vertex constraints, and storing simplification cost in a vertex-indexed form.

The researchers approached the sampling problem through the geometry of a triangular mesh. A dense mesh, obtained from the surface parametric equation or CAD model, served as the representation of the free-form surface. Rather than placing measurement points independently, the method simplifies this mesh and then uses the retained vertices as sampling points. This is a meaningful design choice: if the simplified mesh preserves both geometric information and reasonable spatial coverage, then its vertices become a natural measurement plan rather than an arbitrary set of probe locations. The simplification procedure was built around edge contraction, but with an important modification. Conventional quadratic error metric simplification determines a new contraction target vertex by solving for the position that minimizes the accumulated squared distance to adjacent triangle planes. That calculation can require an invertible coefficient matrix, and when such conditions are not satisfied, additional selection rules are needed. To reduce this source of complexity, the proposed method restricts contraction targets to existing vertices. An edge is contracted by mapping one vertex to another existing vertex, and the simplification cost is stored by vertex sequence rather than by all edge sequences. This changes the computational character of the simplification step while keeping the retained vertices on the original surface, which matters because those vertices later become the measurement points.

The authors gave also careful attention to boundary behavior. Boundary edges and internal edges are distinguished, with boundary vertices constrained to pair only with other boundary vertices. Corner vertices, defined by boundary-edge angles below 90 degrees, are excluded from vertex-pair contraction. These rules are not incidental implementation details; they protect the mesh from degradation during simplification and preserve the surface boundary structure that would otherwise be vulnerable during repeated contractions. Quadratic error alone was not sufficient for the authors’ purpose, because it can assign the same simplification cost to local configurations with different geometric features. Curvature was therefore incorporated through quadratic mean curvature, calculated from discrete Gaussian and mean curvature. Yet curvature weighting by itself can drive too many sampling points into high-curvature zones and leave low-curvature regions sparsely represented. The central technical addition is the area preservation factor, defined from the ratio between area lost during contraction and newly generated area. By linking contraction cost to triangle-area change, the method connects a local mesh operation to the spatial region effectively represented by each sampling point. This analytical strategy has a direct scientific consequence: it prevents feature-sensitive sampling from becoming excessively clustered, because each retained vertex is encouraged to carry a more balanced surface patch. Simulation on a wave surface provided the first evaluation. The researchers defined the surface over a square domain and discretized it into a dense triangular mesh before simplification. When mesh simplification used only QEM, elongated triangles appeared. Adding QMC increased concentration differences between high- and low-curvature regions, whereas the proposed multi-constraint method reduced the area ratio markedly, confirming the role of the area preservation factor in improving distribution uniformity. Accuracy followed the same pattern. Using 64 sampling points, the proposed method produced the smallest reconstruction errors, with a PVE of 0.13 mm and an RMSE of 0.033 mm, clearly lower than the comparison strategies. When the number of sampling points was varied, the proposed method retained a steadier accuracy response. Even after a substantial reduction in sampling points, the increase in PVE remained comparatively modest. The authors used a robot-milled wave surface on a polyoxymethylene block, followed by CMM measurement. A high-density uniform point set provided the reference surface parameters against which the reduced sampling strategies were compared. With 64 points, the proposed strategy reconstructed the surface with maximum and mean milling errors close to the high-density reference values. These deviations met the stated accuracy criteria and were closer to the reference values than those obtained by uniform, QEM, or edge-constrained sampling.

The engineering applications and implications of the work of Harbin Institute of Technology scientists are in improving contact-based inspection of free-form surfaces, especially when measurement accuracy must be maintained without collecting an unnecessarily large number of probing points. It is relevant to CMM inspection and contact-probe-based on-machine measurement, where point distribution controls inspection time, reconstruction accuracy, and surface-quality evaluation.   We can think of a direct engineering implication is that inspection programs for known free-form surfaces can be made more efficient by deriving sampling points from a simplified mesh rather than using uniform or blind sampling. Uniform sampling is simple and widely usable, although it assigns the same sampling logic to regions whose geometry may vary substantially. By simplifying the CAD- or parametric-surface mesh under multiple constraints and using the retained vertices as measurement points, the sampling plan remains tied to the actual geometric character of the part rather than to a fixed grid. For industrial inspection, that is important because many free-form surfaces contain regions where curvature, local shape variation, and surface area coverage must be balanced rather than considered separately.

Optical measurement methods can be difficult when surfaces are reflective, blackened after heat treatment, or measured under less controlled in-site conditions, while chromatic confocal measurement can be limited by surface slope. The new study therefore positions contact measurement as a robust option for such cases. A better adaptive sampling strategy could help make contact-based inspection more practical in machining environments where dense probing would consume too much time but sparse probing could miss important geometric deviation. Another application is in reducing inspection time while preserving useful accuracy. In the simulation, the proposed method achieved lower reconstruction error than uniform sampling, QEM-based sampling, and edge-constrained sampling when using the same number of sampling points. The sampling-point quantity analysis also indicates that the method retains comparatively stable accuracy even when the point count is reduced. A robot-milled wave surface was measured using a CMM, and the proposed sampling strategy gave maximum and mean milling errors close to the high-density reference measurement. The reported deviations from the high-density reference remained within the stated accuracy criterion, supporting practical inspection of machined free-form surfaces rather than only numerical reconstruction. By combining quadratic error, quadratic mean curvature, and an area preservation factor, the method balances surface-feature sensitivity with sampling uniformity. The area preservation factor makes the sampling distribution more balanced, helping each measurement point represent a reasonable surface patch.

High-manganese austenitic steels are important engineering materials because their mechanical response reflects more than conventional dislocation motion; it also depends on deformation mechanisms that vary with alloy chemistry, crystallographic texture, prior deformation, and temperature. Among these materials, twinning-induced plasticity steels occupy a particularly important place. Their combination of strength and ductility arises from the interaction between dislocation slip and deformation twinning, where twins formed during plastic deformation can subdivide grains, obstruct further dislocation movement, and sustain strain hardening over a large deformation range. For engineering use, however, monotonic strength and ductility are not sufficient. Components made from rolled sheet steels may experience repeated loading during service, and their reliability depends on how cyclic plasticity accumulates damage over many loading and unloading cycles. This is especially important in the low-cycle fatigue regime, where plastic strain is not negligible and the material response can evolve substantially during the early stages of cyclic loading. In high-Mn TWIP steels, low-cycle fatigue is complicated by the fact that the same microstructural features that improve tensile deformation may also alter crack initiation, cyclic hardening, hysteresis behavior, and fatigue life. In Fe-Mn-C-based TWIP steels, stacking fault energy is strongly linked to the active deformation mechanism. At room temperature, the material may deform through a coupled contribution of twinning and dislocation slip. At elevated temperature, increased stacking fault energy can suppress twinning and shift the response toward slip-dominated plasticity. A pre-strain that is beneficial at room temperature, for example, may not have the same effect when twinning is reduced and cyclic deformation is governed mainly by slip.

In a recently published paper in International Journal of Fatigue, Professor Di Song, Professor Michael Vormwald, and Professor Heinz Thomas Beier from the University of Electronic Science and Technology of China and Technical University of Darmstadt, developed a fatigue life-prediction model for Fe-19Mn-0.4C-1.7Cr-1.2Al high-manganese steel that accounts for temperature-dependent changes between TWIP-assisted deformation and slip-dominated deformation. The model separates damage from pre-strain and cyclic loading, introduces distinct parameters for slip and twinning effects, and uses the plastic strain range at the twentieth cycle as the key damage descriptor. The research team examined Fe-19Mn-0.4C-1.7Cr-1.2Al high-manganese steel, cut from rolled sheet into specimens loaded at different orientations relative to the rolling direction. At room temperature, the calculated stacking fault energy places the material in a regime where twinning-induced plasticity and dislocation slip coexist. At 449 K, the calculated stacking fault energy rises well above the TWIP threshold, so the elevated-temperature tests probe a condition in which dislocation slip is expected to dominate.

The authors performed tensile testing at 449 K which showed a clear shift in the strength-ductility balance. Yield strength and ultimate tensile strength decreased compared with 295 K, while elongation to failure increased substantially. The orientation dependence remained visible: the 90° direction showed the highest strength, whereas the 0° direction gave the greatest ductility. The elevated temperature therefore did not erase anisotropy, but it changed the mechanical scale on which anisotropy was expressed. The authors also observed that Young’s modulus at 449 K was slightly higher than at 295 K in their measurements, while treating this trend cautiously because more temperature levels and more precise measurements would be needed to define a general temperature dependence. The cyclic tests were strain-controlled and fully reversed, using several total strain amplitudes at 449 K. The researchers compared non-pre-deformed specimens with specimens subjected to tension-unloading pre-strain before fatigue loading. At the same strain amplitude, the 90° specimens generally developed the highest stress amplitudes, while the 45° specimens showed narrower hysteresis loops and lower plastic strain ranges. Lower plastic strain range in the 45° orientation corresponded to longer fatigue life, giving the early-cycle plastic strain response stronger interpretive value than stress amplitude alone. They found at 449 K, stress amplitudes were reduced relative to room temperature, but plastic strain ranges increased. The material also showed stronger cyclic hardening during the initial tens of cycles under elevated-temperature loading. Pre-strain at 449 K increased the initial stress amplitude in cases above 1% strain amplitude and made the subsequent stress-amplitude evolution smoother by reducing the hardening rate. Additionally, when they increased the temperature from 295 K to 449 K shortened fatigue life across the tested orientations, but not in a uniform manner. The 0° orientation showed a stronger life reduction at higher strain amplitudes, while the 90° orientation exhibited a more even reduction across strain amplitudes. At 295 K, prior work of the authors showed that pre-strain could improve fatigue life, especially under lower strain amplitudes. At 449 K, pre-strain markedly reduced fatigue life, with the sharpest reduction in the 90° direction and under higher strain amplitudes.

The team built the life-prediction model around that mechanistic distinction. Classical Basquin, Coffin-Manson, and Basquin-Coffin-Morrow descriptions could not unify the data across orientation, loading history, and temperature. Instead of relying on half-life quantities alone, the proposed model uses the plastic strain range at the twentieth cycle as a damage-related parameter. This choice discriminates different loading modes with no extra variables incorporated, captures early stabilized cyclic plasticity before later convergence obscures orientation and pre-strain effects. The model divides damage into pre-strain-induced and cyclic components, introduces coefficients for slip and twinning contributions, and incorporates temperature through the transition in stacking fault energy and deformation mechanism.

The findings of Professors Song, Vormwald, and Beier have direct engineering relevance for the use of high-manganese TWIP steels in structural components that experience cyclic loading under non-room-temperature conditions. These steels are attractive because they can combine high strength with large ductility, but the research work shows that their fatigue performance depends strongly on how temperature changes the active deformation mechanism. For engineers, this means that fatigue design cannot rely only on room-temperature strain-life data when the component may operate at elevated temperature. A part that benefits from twinning-assisted deformation at 295 K may behave differently at 449 K, where dislocation slip becomes dominant and fatigue life is reduced. This is especially important for rolled sheet components, where the loading direction relative to the rolling direction can influence cyclic stress response, plastic strain range, and fatigue life. The study shows that the 45° orientation can exhibit lower plastic strain range and longer fatigue life under comparable cyclic conditions, while the 90° orientation often develops higher stress amplitudes and can be more sensitive to damaging pre-strain at elevated temperature. In practical forming and component layout, this suggests that the orientation of critical load paths should be considered during design, especially for parts expected to undergo repeated plastic strain.

The results are also relevant to manufacturing routes that involve pre-deformation, forming, or prior plastic strain before service. At room temperature, pre-strain may improve fatigue life in high-Mn TWIP steel because deformation twinning can hinder later dislocation slip. At elevated temperature, however, the same pre-strain can sharply reduce fatigue life, especially in the 90° orientation and under higher strain amplitudes. This has clear implications for formed automotive, aerospace, or civil engineering components made from high-Mn steel: the service temperature and loading direction must be considered together with the forming history. The life-prediction model developed in the study is also useful from an engineering assessment perspective. By incorporating deformation mechanism, temperature, loading orientation, pre-strain, and early-cycle plastic strain range, it provides a more realistic basis for fatigue-life estimation than classical strain-life models alone. Such a model can support safer component design, material qualification, and fatigue assessment when high-Mn TWIP steels are used in structures exposed to variable thermal and cyclic loading conditions.

Dynamic load identification is a demanding inverse problem in structural vibration analysis because applied forces often cannot be measured directly and must be determined from the structural responses they produce. In continuous systems, accelerations recorded at selected points may be the only available evidence of loads acting at unknown positions and changing rapidly with time. The problem becomes more demanding when several loads act at the same time. For a single excitation, the problem is often more manageable, especially when the load position is known or the event is a short impact. Under multi-point excitation, however, the responses produced by different loads are superposed through the structure’s modal behavior. When every possible loading-point combination must be tested, the computational burden can become a practical barrier, especially in continuous systems where engineers need to recover both the force position and the force history. In a recently research paper published in International Journal of Mechanical Sciences Dr. Zhengshu Wang, Professor Jinhui Jiang, Dr.  Wen Jing, and Dr. Rutong Chen from the Nanjing University of Aeronautics and Astronautics, developed a rapid dynamic load identification method that combines numerical integration, FastICA, and modal shape comparison. The technically distinct element is the formulation of modal load histories as mixed observations of statistically independent physical load signals, allowing multiple load contributions to be separated before localization. The new method then uses proportional relationships among modal shape values to determine load locations and recover corresponding magnitudes. This differs from prior search-based localization strategies by decoupling the multi-load problem and reducing the growth of candidate-location screening from combinational to load-wise.

The research team begins by establishing a mapping from measured structural response to modal loads. The researchers formulate the continuous system in modal coordinates and use a numerical integration procedure, illustrated through the Newmark explicit method, to express the measured response as a convolution-like relation involving modal load histories. Solving this inverse relation yields modal loads of different orders. This step is important because it moves the problem from physical response space, where several load effects are entangled, into a modal-load representation where the mixing relationship can be treated more directly.  The key design choice is the treatment of the modal load vector as a set of observed mixed signals. Since the contribution of each physical load to each modal load is proportional to the modal shape evaluated at the load location, the transformation between physical loads and modal loads has the same mathematical form as the standard Independent Component Analysis model. FastICA is then used to separate statistically independent load components from the identified modal loads. The scientific consequence is clear: multi-load localization no longer requires the simultaneous testing of all possible loading-point combinations. Each separated component can be associated with its own modal-shape signature and localized individually. After separation, modal shape comparison supplies the spatial interpretation. Because ICA cannot determine the absolute order or scale of the separated sources, the method does not rely on direct equality between the separated mixing matrix and the physical modal transformation matrix. Instead, it uses proportional relationships among modal shape values. For each separated load component, the method constructs a modal-shape deviation measure and searches for the location where that deviation is minimized. Once the position is identified, the load magnitude can be recovered by accounting for the modal shape relationship.

Afterwards, the researchers examined pairs of sinusoidal loads with different frequencies, combinations of sinusoidal and impact loading, broadband random loading combined with impact loading, and broadband random loading combined with sinusoidal excitation. Across these cases, the method localizes loads with either exact agreement or small errors relative to the beam length, even when response signals include 5% measurement noise. The reconstructed load histories also remain close to the applied histories, with signal-to-noise ratios reported for continuous loads and peak relative error metrics used for impact loads. One interesting point is that the identified load signals sometimes have better signal-to-noise measures than the original response signals, which the authors attribute to regularization in the inverse step and the noise-reduction behavior of signal separation. Computational efficiency is part of the method’s design logic. The authors compared the method against a modal load reconstruction approach using optimization. Even when the proposed method searches all possible loading positions, they found it requires less time than the comparison method, while approaching the computational time needed for magnitude reconstruction when locations are already known. The reason is structural: FastICA decouples the multi-load localization problem, reducing the required search from a combinational form to a load-by-load search, and modal shape comparison keeps each screening operation simple. The complex wing model extends the evaluation beyond the one-dimensional beam. With two sinusoidal loads applied to a finite element wing structure, the method identifies both load coordinates accurately and reconstructs the corresponding load histories with high reported signal-to-noise ratios. A separate model-error analysis introduces a 10% thickness error into the beam model. The team finds that localization remains relatively stable, while magnitude reconstruction is more affected. This distinction is physically reasonable within the proposed method because spatial localization depends mainly on the modal shape distribution, whereas magnitude reconstruction depends more directly on quantitative model parameters. Experimental verification on a simply supported beam adds a practical test: a sinusoidal excitation and a double-impact load are applied, accelerations are measured at five points, and the method identifies locations within the scale of the force sensor while reconstructing both load types with reported accuracy and faster computation than the comparison approach.

The findings of the Nanjing University of Aeronautics and Astronautics scientists are most directly relevant to engineering situations where structures experience unknown, time-varying loads and where placing force sensors at the actual loading positions is difficult or impossible. The new method uses measured vibration responses to infer both the location and time history of multiple dynamic loads, which makes it useful for structural vibration analysis, load monitoring, fault diagnosis, and structural health assessment in continuous systems. A major application is in aerospace structures, especially because the method was evaluated on a complex wing model as well as a simply supported beam. This extension matters because wing structures, including skins, spars, and ribs, may experience several simultaneous excitations during service. A rapid method that can separate and locate multiple loads from response signals could help engineers identify where dynamic excitation is entering the structure and how those loads evolve over time. This is especially useful for vibration control, fatigue assessment, and post-event interpretation after impacts or abnormal loading events.

The new method also has value for mechanical systems exposed to mixed loading conditions. The simulations include sinusoidal loads, random broadband loads, impact loads, and combined load cases. These conditions resemble practical situations where rotating machinery, aerodynamic excitation, operational vibration, and accidental impacts may overlap. Because the method can separate statistically independent load contributions, it could help distinguish different sources of excitation rather than treating the measured response as one inseparable vibration record. Another application is rapid localization for monitoring systems with limited sensor access. The investigators emphasized that conventional multi-load localization may require searching many possible loading-point combinations, while the FastICA-based strategy decouples the problem and improves computational efficiency. This makes the approach relevant to online or near-real-time diagnostic workflows, where localization must be completed within practical computational limits. The experimental beam validation gives the method practical weight beyond numerical simulation. Using measured responses, together with the reported noise and model-error tests, the paper shows that the method can support response-based identification of multiple dynamic loads when direct force measurement is difficult. Its most natural role is as a response-based tool for rapid localization and load-history reconstruction in continuous systems.

Elastic beam systems are important in mechanical, civil, marine, aerospace, and manufacturing engineering because many slender load-bearing components can be understood through their bending and vibration behavior. A beam may carry load, transmit motion, support rotating or reciprocating machinery, and receive dynamic input from more than one location at the same time. Under these conditions, vibration depends not only on isolated resonance, but also on how dynamic energy moves among connected substructures. Traditional vibration reduction strategies for beam-like systems rely on devices tuned to specific frequency ranges and these devices can be effective when the excitation environment is well defined, but practical structures do not always operate within a narrow or fixed frequency band. The problem becomes more involved when the system contains primary and secondary beams that are mechanically connected, because the secondary components are not just passive appendages. If their connection to the main structure is properly designed, they can participate in vibration redistribution and may serve as part of the vibration control mechanism itself.

In a recent research paper published in Mechanical Systems and Signal Processing, Professor Yuhao Zhao, Dr. Yang Cao, Dr. Mingfei Chen, Dr.  Cunhong Yin from Guizhou University working together with Dr. Zongfang Wu from Pipe China Southwest Pipeline Company Guiyang Maintenance and Emergency Repair Center developed a theoretical and experimental nonlinear coupling beam system subjected to two excitation sources. They introduced magnetic coupling nonlinearities between two excited primary beams and an internal secondary beam, creating nonlinear vibration transfer pathways within the structure. A technically important feature is that the modeling and experiments treat phase-difference effects within the same coupled beam system. The work also demonstrated that annular magnetic elements can realize the cubic nonlinear stiffness needed for the proposed coupling mechanism.

The researchers investigated a three-beam system in which Beam-A and Beam-B served as the excited primary sub-beams, while Beam-C was used as an internal secondary beam.  The two excitation sources had the same angular frequency but different phases, allowing the coupling effect between inputs to be examined. They developed a theoretical model using a Lagrangian formulation and represented the beam displacements through modal expansions, included kinetic and potential energy contributions from the three beams, accounted for the linear coupling springs, and also introduced nonlinear stiffness terms associated with the coupling nonlinearities. The magnetic nonlinear couplings were modeled with both linear and cubic nonlinear stiffness components. To check the reliability of the theoretical procedure, they compared the Lagrangian predictions with results obtained through a Galerkin-based approach and the close agreement between the two modeling routes supported the subsequent parameter studies. The team found under one parameter set, the nonlinear coupling system remained in what the authors describe as a linear vibration control state while under another parameter set, nonlinear response features appeared at specific excitation frequencies, including multiple local peaks at a single excitation frequency. The nonlinear state was further identified as quasi-periodic through phase-diagram analysis. This distinction matters because it shows that vibration control is not tied to a single response regime; the same coupling architecture can suppress vibration while the dynamic response remains linear-like or becomes distinctly nonlinear. The role of Beam-C is especially important because Beam-A and Beam-B receive vibration energy directly from the two excitation sources, while the nonlinear connections create routes through which part of that energy can transfer into Beam-C. In that sense, Beam-C and the coupling nonlinearities operate together as a nonlinear vibration absorber embedded within the coupled beam system. The design choice of placing nonlinear transfer paths between the excited primary beams and the internal beam directly changes the energy redistribution mechanism, allowing the secondary beam to participate in suppressing the resonance response of the primary substructures.

The nonlinear stiffness parameter also changed the vibration control behavior. Increasing the nonlinear stiffness improved the vibration suppression effect, but it also increased the likelihood of nonlinear phenomena appearing in the response. The study frames this as a design balance: stronger nonlinear coupling may provide greater suppression, while also requiring attention to the resulting nonlinear dynamic states. The phase difference between the two excitation sources produced another major finding. As the phase difference changed, response peaks shifted in magnitude, and the effect was not identical across frequency bands. When the vibrations generated by the two sources were effectively in phase at coupling points, the response could be amplified; when they were closer to opposite phase, the response could be suppressed. The researchers examined this behavior across several nonlinear stiffness values and showed that the response envelope still remained lower than that of the original beam system when coupling nonlinearities were introduced. Beam-A and Beam-B displayed similar peak-value patterns, but the phases at which their maxima occurred were not the same, meaning that the two primary sub-beams did not reach their largest responses simultaneously. The authors performed experimental work  using stretching and compression tests which confirmed that the magnetic coupling elements have suitable cubic nonlinear stiffness characteristics. The full testbed included Beam-A, Beam-B, Beam-C, springs, magnetic nonlinear couplings, two vibration sources, acceleration sensors, rigid supports, NI equipment, power amplifiers, and computer control. Ten independent experiments were carried out because the initial phase difference between the two excitation sources could not be directly controlled. The experimental acceleration data supported the theoretical prediction, showing reduced peak responses in the tested resonance regions and different stabilized states under different initial phase relationships. The time-domain acceleration responses also changed from sinusoidal behavior after the nonlinear couplings were introduced, consistent with the predicted nonlinear dynamics.

The findings of Professor Yuhao Zhao and colleagues are directly relevant to engineering systems in which beam-like structures are exposed to vibration from more than one source. Many real structures do not receive excitation from a single idealized point. Machinery is often mounted through several supports, rotating equipment can transmit vibration through multiple bases, and coupled frame or beam assemblies may experience simultaneous dynamic inputs with different phase relationships. The study is therefore useful for systems such as pipeline supports, shafting structures, robotic arms, frame-type offshore platforms, shipboard mechanical assemblies, and high-speed industrial equipment where vibration travels through connected structural members. A major engineering application is the design of built-in vibration control systems using secondary beams or internal structural members. Instead of adding a separate absorber, the study shows that an existing secondary beam can be connected to the primary beams through nonlinear elements and used as part of the vibration suppression mechanism. This is valuable in compact engineering structures where extra space, added mass, or external vibration-control devices may be undesirable. The secondary beam becomes dynamically useful because the nonlinear connections create pathways for vibration energy to move away from the excited primary beams. The magnetic coupling nonlinearities are also important from a design perspective and by using annular magnets to create cubic nonlinear stiffness, engineers can introduce nonlinear vibration transfer without relying only on mechanical springs or contact-based elements. Such a concept could be adapted for adjustable or replaceable vibration-control modules, especially where tuning the nonlinear stiffness is needed to match the operating conditions of a structure. The study also indicates that nonlinear stiffness selection is an important design consideration, since increasing stiffness can improve vibration suppression while making nonlinear response states more likely.

Another important application concerns multi-source vibration diagnosis and control. The study demonstrates that phase differences between two excitation sources can strongly affect resonance peaks, and that the effect varies across frequency bands. This means that engineers evaluating vibration in machinery-supported beams, pipelines, or frame structures should consider not only excitation frequency and amplitude, but also phase relationships between input sources. In practice, controlling or adjusting phase differences between excitation sources may become another route for reducing vibration. The findings therefore support a design approach in which secondary substructures are deliberately integrated into the vibration-control function of coupled beam systems.

Acoustic vortex beams are a class of structured sound fields distinguished by a helical phase distribution around the propagation axis. Unlike ordinary acoustic waves, which usually propagate with a comparatively uniform phase front, these beams possess an azimuthally varying phase term of the form e^ilθ, where θ is the angular coordinate and l is the topological charge. This topological charge defines the number of phase windings around the axis and is directly associated with the orbital angular momentum carried by the beam. Because of that helical structure, acoustic vortex beams have attracted interest in areas where controlled wavefront shaping and angular momentum transfer are important. A central challenge in this area has been the practical generation and detection of acoustic vortices without relying on cumbersome auxiliary systems. Most existing approaches generate them either through active transducer arrays or through passive structures designed to reshape incident sound fields. Rotating aerodynamic sources suggest a different possibility, since their tonal radiation is already shaped by periodic motion and strong spectral components tied to rotation frequency and blade passing frequency. In a recent research paper published in Applied Acoustics, Associate Professor Lianyun Liu and Professor Zhigang Chu from Chongqing University examined whether the tonal sound produced by rotating rotors can be understood as acoustic vortex fields. They combined this idea with a virtual rotating receiver measurement strategy that extracts Doppler-shifted spectra and adjacent-phase differences from a static circular microphone array. They also demonstrated single-rotor and coaxial twin-rotor generation of vortex beams, including multiplexing at shared or distinct operating frequencies and contactless inference of rotor speed and blade number.

The authors began from the tonal sound emitted by rotating rotors. Instead of treating rotor noise as something tied to one particular source mechanism, they considered the rotor more broadly as a periodically rotating acoustic source whose tonal content appears at the rotation frequency, the blade-passing frequency, and higher harmonics. From there, Liu and Chu showed that the phase lag between observers placed at different azimuthal positions around the rotor naturally gives rise to a helical phase structure. This means that the blade-passing component behaves as an acoustic vortex, while higher harmonics and other tonal components can also be understood within the same theoretical description. More importantly, the full tonal field can be described as a superposition of acoustic vortex beams, which becomes the central theoretical claim of the paper.

Liu and Chu designed the experiments to test that idea under realistic rotor operating conditions. A coaxial microphone array with inner and outer rings recorded both near-field and far-field sound from eight different rotors, including disks, axial propellers, axial fans, a centrifugal impeller, and a fiber impeller. They then applied the virtual rotating receiver method to reconstruct Doppler-shifted spectra while suppressing plane-wave contamination by removing the azimuthal average of the static microphone signals. That choice was important because once the spatially uniform part of the signal was removed, the remaining spectra preferentially retained sound with azimuthal structure, which is where the helical signature of acoustic vortices is expressed.

The near-field measurements revealed a consistent pattern. For all tested rotors, the blade-passing components were strongest on the inner array, and second harmonics were often clearly present as well. When the virtual rotating receiver spectra were examined, the dominant spectral ridges followed the Doppler behavior expected for vortex-carrying sound fields, and their slope with observer rotation matched blade number. At the same time, the phase difference between adjacent virtual receivers changed sign at the point corresponding to zero Doppler shift, indicating reversal of orbital angular momentum handedness. Those two observations reinforced one another: the frequency shift established the expected rotational relation, while the phase reversal showed that the field behaved as a true acoustic vortex rather than as an ordinary tonal signal. The authors found that axial propellers retained strong blade-passing vortex signatures, whereas several non-axial rotors were instead dominated by lower-order vortex components or showed much weaker high-charge content at distance. The authors related this to the weakening or disappearance of blade-passing components in far-field spectra and to the stronger decay associated with vortices of higher topological charge. That contrast gave the results a useful internal coherence: the same theory still applied, but what survived into the far field depended on rotor type and on which tonal components propagated strongly enough to remain distinguishable from other aerodynamic noise.

Liu and Chu examined two representative cases in more detail by filtering the sound at selected harmonic frequencies. For the centrifugal impeller, strong near-field vortices appeared across several orders, and comparison between the total spectral amplitude and the vortex-associated amplitude showed that most of the near-field acoustic energy at those frequencies was carried in vortex form. In the far field, that correspondence became much weaker at higher orders, and the paper used a partial-wave interpretation to explain why measurable energy could still appear there without indicating strong ideal vortex propagation. The large axial propeller showed a different pattern: even after propagation to the far field, several vortex orders remained prominent, and the blade-passing component was clearly dominant. The coaxial twin-rotor experiments then extended the argument from generation to multiplexing. Different rotor pairings produced simultaneous vortex beams with different topological charges at the same operating frequency, or at different operating frequencies for frequency-separated transmission. The virtual rotating receiver analysis resolved these as distinct channels rather than as a single ordinary tonal field. A simple coding demonstration then showed that the vortex-associated signal levels could distinguish whether one or both rotors were active, giving a proof-of-principle multiplexing scheme based on rotor-generated acoustic vortex beams.

The research work of Liu and Chu shows that acoustic vortex generation can arise intrinsically from rotor tonal radiation, with the vortex operating frequency and topological charge determined directly by rotation frequency and blade number. The methods used in the study are also important. The paper shows that acoustic vortex generation can be identified using a static microphone array processed with the virtual rotating receiver method, without needing a specially built vortex source. That gives the theory a practical experimental basis and clarifies why the authors emphasize contactless measurement of rotor speed and blade number. In their industrial fan example, ordinary sound pressure level peaks were not sufficient on their own, but the virtual rotating receiver spectra provided the information needed to determine both the rotation frequency and the blade number.  The authors’ findings have practical implications and rotor-generated sound can be exploited as a structured acoustic signal with measurable and usable vortex properties. This is important because in practical terms this supports contactless measurement of rotor speed and blade number, enables multiplexed acoustic communication using distinct vortex channels, and points toward future strategies for rotor noise control and acoustic energy harvesting through manipulation of acoustic orbital angular momentum.

Ultra-precision diamond turning has gradually expanded from the production of smooth optical surfaces toward the controlled generation of micro- and nano-scale topographies whose function depends as much on geometry as on surface finish. Within that shift, microchannels on cylindrical substrates have become particularly relevant, not simply as geometric features but as structures that influence wettability, fluid transport, and thermal behavior. Their fabrication, however, imposes a combination of constraints that are not easily reconciled within existing machining paradigms. The tool must follow controlled trajectories at small scales while maintaining stability over comparatively large displacements, and the resulting surface must remain free of discontinuities or residual artifacts that degrade performance. Fast tool servo systems have long been used to introduce controlled, high-frequency motion into ultra-precision cutting. When extended to two degrees of freedom, such systems offer a route toward generating more complex surface morphologies, including structures with coupled depth and orientation features but the transition from conceptual capability to practical implementation has proven difficult. A persistent issue lies in the inability of conventional designs to simultaneously deliver large stroke, high operational frequency, and cleanly decoupled motion along orthogonal axes. Serial configurations tend to provide good decoupling but suffer from inertia accumulation and limited bandwidth, whereas parallel architectures can provide dynamic advantages but introduce coupling between axes that can increase control demands and limit positioning precision. These limitations become especially pronounced in cylindrical microchannel fabrication. Existing approaches based on slow tool servo motion allow large displacements but operate at low speed, while methods relying on fast tool servo motion often reconstruct channels through overlapping micro-dimples, which can leave residual material at intersection boundaries and affect surface quality. The technical gap is therefore in the structural design: a mechanism is needed that can amplify the inherently small displacement of piezoelectric actuators, maintain high dynamic responsiveness, and preserve independent motion along each axis without adding substantial structural or control complexity.  In a recent research paper published in International Journal of Mechanical Sciences, Mr. Qinghou Cheng, Professor Yangkun Zhang, Professor Yingxue Yao, and Professor Yang Yang from the Harbin Institute of Technology (Shenzhen) developed a parallel two-degree-of-freedom piezoelectric tool holder that integrates displacement amplification and mechanical motion decoupling within a compact flexural structure. The design achieves large stroke and high resonant frequency without relying on control-based decoupling strategies. They also established a coordinated turning method that uses independent axis control to fabricate microchannels with tunable depth and orientation.

Briefly, the proposed tool holder is organized around a parallel configuration in which two input stages are independently driven by piezoelectric stacks and connected to a common end-effector through a system of actuated and passive flexural joints. Each axis is guided by a pair of prismatic mechanisms that enforce motion along a single direction while suppressing parasitic displacement. The structural novelty lies in the integration of two distinct parallelogram mechanisms within each actuated joint: one provides displacement amplification through controlled flexural deformation, while the other enforces kinematic guidance to maintain linear motion. This arrangement produces an interesting consequence. Because displacement amplification arises from geometric transformation rather than additional structural modules, the mechanism avoids the increase in inertia typically associated with extended amplification stages. At the same time, the use of passive flexural joints aligned with each axis ensures that motion transmitted to the end-effector remains directionally constrained, effectively achieving mechanical decoupling without reliance on control algorithms. To understand how this configuration behaves under load and excitation, the authors construct a dynamic stiffness model in which the mechanism is discretized into rigid bodies and flexural elements. Each flexible component contributes a frequency-dependent stiffness matrix, and the global behavior emerges from the equilibrium of forces and moments across all nodes. This modeling approach allows both static characteristics, such as stroke and stiffness, and dynamic properties, including resonant frequencies, to be evaluated within a unified framework.

Systematic variation of key geometric parameters reveals a fundamental trade-off between stroke and resonant frequency. Increasing the tilt angle of the amplification mechanism, for example, enhances frequency response but eventually reduces displacement, while changes in beam thickness alter stiffness and frequency in opposite directions. The optimization strategy therefore does not pursue a single objective but instead balances competing requirements by imposing constraints on minimum stroke and stiffness while maximizing resonant frequency. The resulting design achieves a calculated stroke of approximately 136 μm and a first resonant frequency near 924 Hz, indicating that the structural configuration can sustain both large displacement and high-speed operation. The authors performed finite element simulations to provide a secondary check on these predictions and found that the simulated stroke, stiffness, and modal frequencies remain within roughly ten percent of the analytical values, with discrepancies attributed mainly to the idealization of rigid bodies in the theoretical model. Stress analysis shows that peak stresses remain below material yield limits, suggesting that the compliant elements operate within a stable deformation regime. The team found that under low-frequency sinusoidal excitation, the tool holder produces strokes of 131 μm along the X-axis and 141 μm along the Y-axis, with cross-coupling ratios below one percent in both directions. This low level of coupling is particularly significant because it reflects structural rather than control-based decoupling which means that each axis can be driven independently without compensation strategies. Moreover, their dynamic testing further indicated that the system maintains resonant frequencies above 700 Hz in both directions, with stable amplitude response up to several hundred hertz. Positioning resolution, measured using high-precision displacement sensors, reaches approximately 18–19 nm, demonstrating that large stroke capability does not come at the expense of fine positional control. The study then moves from characterization to application by integrating the tool holder into a cylindrical turning setup. Here, the two axes are assigned distinct functional roles: motion along the cutting depth direction controls hierarchical features of the microchannels, while motion along the feed direction governs their orientation relative to the rotating workpiece.

Fabrication experiments illustrate how different signal inputs translate into surface morphology. Single-axis excitation along the depth direction generates periodic variations in channel depth, while excitation along the feed direction produces oriented patterns governed by spindle rotation. When both axes are actuated simultaneously, the resulting structures exhibit controlled variation in both depth and orientation, with measured dimensions closely matching theoretical predictions. The consistency of wavelength and amplitude across different operating conditions indicates that the mechanical decoupling remains effective during actual cutting, not only in isolated tests. The most immediate engineering value of the work of Harbin Institute of Technology researchers is in the fabrication of cylindrical microchannels whose geometry can be precisely controlled especially that the tool holder allows the channel depth, or hierarchical feature, and the channel orientation to be adjusted independently or in combination. The study links such microchannel structures to wettability control, directional droplet transport, and heat dissipation, all of which depend strongly on how fluids interact with textured surfaces. In that sense, the cylindrical surface is not just patterned; it is given a programmable microgeometry that can influence how liquids spread, adhere, or move along the surface.

A second important application is high-efficiency machining of complex microstructured surfaces. Because the tool holder combines large stroke, high operational bandwidth, nanoscale positioning resolution, and very low cross-axis coupling, it can generate controlled two-dimensional tool paths during turning rather than relying on slower servo motion or overlapping micro-dimples. This is important because the investigators are addressing both productivity and surface quality at the same time: the process can produce orientation-controllable hierarchical microchannels continuously on curved surfaces while avoiding the residual boundary defects associated with dimple-overlap methods. The engineering implication is that complex cylindrical surface morphologies can be formed more directly, with the final geometry governed by the coordinated input signals rather than by mechanical limitations of the actuator system.

Strain recovery becomes difficult to interpret when stress and temperature change at the same time, because the alloy is no longer moving along the simple actuation or superelastic routes that usually supply the data for constitutive calibration. High-temperature shape memory alloys, where non-proportional thermo-mechanical histories matter for practical use and where standard calibration practice still relies mainly on isobaric thermal cycling or isothermal stress cycling. Once temperature gradients, simultaneous loading, and heterogeneous local response enter the picture, a model built only from simple paths has to carry much more than a smooth macroscopic hysteresis loop; it has to represent how transformation, transformation-induced plasticity, and crystallographic slip accumulate inside a polycrystal whose grains do not all respond in the same way.

In a recent research paper published in International Journal of Plasticity, PhD candidate Adrien Cassagne, Professor Dimitris Lagoudas, and Professor Jean-Briac le Graverend from the Texas A & M University, the authors developed a crystal-plasticity mean-field framework for high-temperature shape memory alloys subjected to complex thermo-mechanical loading. The formulation combines thermo-elasticity, martensitic transformation, TRIP, and plastic or viscoplastic slip with an implicit β-transition rule that transfers macroscopic loading to grain-scale response. They also introduced a grain-size-dependent transformation-threshold law and a stress-dependent saturating transformation-strain variable tied to local von Mises stress. Using parameters calibrated from isobaric experiments, they applied the model to out-of-phase and in-phase loading paths in a Ti-rich NiTiHf polycrystal.

The authors build the constitutive framework by decomposing total strain into thermo-elastic strain, transformation strain, TRIP strain, and plastic or viscoplastic strain, then letting each inelastic part evolve from its own driving force at the grain scale. That structure matters because complex thermo-mechanical paths do not redistribute deformation through a single mechanism. A model that lumps the response too early would miss the fact that martensitic transformation changes the local state that then drives TRIP and slip. Cassagne and colleagues keep those channels separate, then couple them through dislocation-dependent transformation resistance and an irrecoverable martensite fraction. In practical terms, that means the model can let the same loading path generate recoverable strain, dislocation-assisted residual strain, and altered transformation thresholds within one common framework.

Two ingredients carry much of the paper’s original modeling work. One is the grain-size-dependent activation law for forward and reverse transformation thresholds. The other is the stress-sensitive coefficient δc that saturates the transformation-strain magnitude as a function of local von Mises stress. Those additions are not decorative. The grain-size criterion gives the model a microstructural basis for smoothing hysteresis across the aggregate, which is especially useful in a mean-field setting where each grain contributes statistically through its volume fraction. The saturation term changes the logic of how stress feeds transformation strain: instead of letting strain continue to scale too aggressively from the moment stress rises, the formulation lets the response build toward a ceiling. The paper shows why that choice is physically useful for the NiTiHf alloy under discussion, where higher stress does not simply translate into unlimited growth of transformation strain once preferred variants have already been activated.

They implement the polycrystal as a 250-grain aggregate with random orientations and a Gaussian grain-size distribution, then connect local grain stresses to the macroscopic state through the implicit β-transition rule. Here the modeling choice is quite telling. The β variable gives each grain access to a non-linear accommodation of internal stress through the aggregate average, so the model can generate harder local response without forcing the sharp domino-like propagation that full-field compatibility can produce. That is exactly the kind of adjustment complex thermo-mechanical paths need, because transformation onset in one part of the aggregate should influence the rest of the material, though not in a mechanically abrupt way. The calibration proceeds on isobaric tests across several stress levels, after which the same parameter set is used for out-of-phase and in-phase paths.

The computed isobaric responses recover the experimentally observed strain evolution with smooth hardening, and the transformation strain rises strongly with stress up to about 470 MPa before approaching saturation, while total strain continues to grow through TRIP. That separation is important: it tells the reader that the model does not treat all extra deformation at high stress as new transformation strain. When the authors move to out-of-phase loading, the simulated loops retain the expected hysteretic form and reproduce the experimentally observed trends in transformation and actuation strain. For in-phase loading, the framework reproduces the unusual first-cycle character associated with an initially self-accommodated martensitic state, where the reverse-transformation peak is largely absent, and it also produces the stronger hysteresis seen on subsequent cycles. The paper treats that first-cycle feature as a real mechanistic signature of variant selection during initial reverse transformation, not as a curiosity at the edge of the data.

The work of Texas A & M University researchers built constitutive model for shape memory alloys as well as reorganizes the way complex loading paths are handled in HTSMAs by tying the description of path dependence to grain-scale activation, stress-dependent saturation of transformation strain, and a self-consistent mean-field transfer of internal stress. That is a different modeling logic from one that begins with separate phenomenological rules for each loading route. Here, isobaric calibration remains the experimental foundation, yet the model attempts to carry that information into non-proportional loading through internal variables whose meaning remains local and mechanistic. For a field that often has to choose between macroscopic tractability and microstructural fidelity, that is a substantive shift in emphasis.

The grain-size dependence is especially meaningful and by embedding transformation thresholds in the statistical grain structure of the aggregate, the authors move hysteresis control away from a purely path-tagged rule and toward a microstructure-linked activation picture. In this paper, that step serves two purposes at once. It gives the mean-field model a way to spread transformation more smoothly across grains, and it gives the constitutive law a basis for responding to complex loading without being written separately for every test type. The Hall–Petch-like form used for the thresholds is part of that logic: smaller relative grains carry different activation costs, which changes how the aggregate enters and exits transformation. That kind of treatment matters for HTSMA design because actuator behavior is often judged at the scale of the whole component, yet the component response is assembled from grain-level events whose sequence changes under mixed thermal and mechanical driving.

The β-transition rule also has broader implications. In this paper it functions as more than a numerical bridge. It gives the model room to represent intergranular hardening through an implicit internal variable, and that in turn lets the homogenized response keep contact with grain-scale stress redistribution during non-proportional loading. For alloys expected to operate under thermal gradients, changing loads, or cycling histories that do not follow textbook isothermal or isobaric paths, that is an important capability. It means the constitutive description can remain sensitive to local accommodation and variant activity without paying the full cost of a full-field calculation for every cycle. The paper does not treat mean-field averaging as a retreat from physics. It treats it as a controlled way to keep the physics that matter most for the problem at hand.

A second implication comes from the in-phase results. The model reproduces the shape of the first cycle that begins from a self-accommodated martensitic state, including the muted reverse-transformation peak. That point reaches beyond one loading protocol. It shows that initial microstructural state can decisively reorganize the first thermo-mechanical cycle, and that a constitutive framework can register that reorganization when variant-resolved transformation remains explicit. The same framework also separates the roles of transformation strain and TRIP across different paths, with the paper showing that stress level during forward or reverse transformation alters where TRIP accumulates most strongly. For people building predictive tools for HTSMA actuators, that matters because actuation, residual strain, and cycle history do not arise from one scalar hysteresis parameter. They emerge from coupled internal processes whose order changes with the loading path.