Low cycle fatigue life prediction remains a challenge in the structural assessment of high-strength alloys used under demanding cyclic loading. For nickel-based superalloys, the difficulty arises because fatigue failure depends on multiple interacting factors, namely, material state, welding history, local deformation, cyclic strain evolution, and accumulated damage, which make life prediction particularly complex. GH4169 is a precipitation-hardened nickel-based alloy used where high strength and resistance to elevated-temperature degradation are required. It is commonly joined by inertia friction welding, but this introduces additional complexity: the weld joint must be evaluated not only as a nominal material but also as a structural region whose cyclic response is shaped by processing and local deformation. Traditional low cycle fatigue prediction has usually relied on strain-life relations, empirical equations, strain-energy approaches, and damage or critical-plane concepts. These methods remain important because they incorporate physically meaningful variables. However, their practical limitation lies in the dependence on empirically fitted parameters and the challenge of capturing fatigue life when multiple interacting factors are involved. In welded GH4169 joints, the prediction problem is made sharper by the need to connect measurable deformation during cyclic loading with final fatigue life in a way that is both data-efficient and physically interpretable.

Deep learning offers an alternative route by extracting patterns from complex datasets without requiring every feature to be predefined. A model trained directly on deformation images may receive information-rich input, but much of that information—such as speckle motion, contrast, and local texture—can be optically complex rather than mechanically decisive for fatigue life. The key question is: what form of experimental information allows neural network to learn a meaningful relationship between cyclic deformation and fatigue life? To address this, Professor Ming-Liang Zhu and Professor Fu-Zhen Xuan from East China University of Science and Technology developed convolutional neural network models in a recent paper published in Fatigue & Fracture of Engineering Materials & Structures for predicting the low cycle fatigue life of GH4169 inertia friction welded joints using three input formats: deformation images, average strain values, and combined average strain with damage fraction. The best-performing model used the combined strain and damage-fraction dataset and achieved the highest reported prediction accuracy in the study, whose framework is shown in Fig.1.

The research team generated a controlled low cycle fatigue dataset from GH4169 homogeneous inertia friction welded joint specimens. The material composition, welding parameters, and tensile properties were specified. Cyclic testing was performed at room temperature under stress control with a stress ratio of -1 and sinusoidal loading. Stress amplitudes ranged from 930 to 1160 MPa, producing fatigue lives from approximately 1,170 to over 25,000 cycles.

In their experimental design, the team simultaneously collected fatigue response data and optical deformation information. Speckle patterns were prepared on specimen surfaces, and a camera-based observation system was used during cyclic loading. Instead of retaining all images, they divided each specimen’s life into five equal periods and selected deformation images associated with peak strain within sampled cycles. This is important because peak strain is a mechanically relevant state within the cyclic response and by using peak-strain-associated images, the dataset was guided toward deformation states more directly related to fatigue accumulation. They generated three data streams from the same experimental foundation. The first consisted of cropped deformation images, used as input to a three-dimensional convolutional neural network (CNN). The second replaced raw images with average axial strain values extracted by digital image correlation (DIC). The third added a damage fraction calculated from the ratio of sampled cycle count to the corresponding fatigue life, following a linear cumulative damage representation. Notably, comparing these inputs was the main methodological strength of the work, because it isolated the effect of image richness from that of mechanically processed information.

The image-based network used convolutional layers to process deformation images from multiple life periods before passing concatenated features to fully connected layers for life prediction. The strain-based and combined-input models achieved higher predictive accuracy. With the smaller image dataset, the test coefficient of determination(R²) was 0.4652; with the larger dataset, it fell to 0.2089. The authors interpreted this as evidence that more image data did not necessarily provide more fatigue-relevant information. Additional images may have introduced optical variation that was less directly connected to fatigue life, increasing the complexity of the learning task when the network had to infer the connection from surface texture to deformation state and then to fatigue life. When average strain values replaced deformation images as network input, prediction accuracy improved markedly. The corresponding one-dimensional convolutional networks reached test R² of 0.8159 for the smaller dataset and 0.9371 for the larger dataset. This improvement is important because the strain values were derived from the same image source that gave poorer results when used directly. It clarifies the role of image processing: digital image correlation acted as a physics-guided feature extraction step, translating optical deformation into a compact variable with direct fatigue relevance.

The combined strain and damage fraction model achieved the strongest prediction. Under the smaller dataset, the test R² increased to 0.8478; and under the larger dataset it reached 0.9560 (Fig.2). The addition of damage fraction allowed the network to receive not only a deformation feature but also a normalized indication of where the sampled state lay within the specimen’s fatigue process. The larger combined dataset produced the most reliable predictions, with test points distributed within the narrower error band. The work therefore supports a clear hierarchy: raw deformation images were less effective, strain values were substantially more informative, and strain values combined with a physics-based damage descriptor yielded the best fatigue life prediction among the tested models.

The findings of East China University of Science and Technology researchers have direct engineering relevance for fatigue assessment of nickel-based welded components, especially where inspection must move beyond visual observation and toward measurable indicators of remaining life. GH4169 inertia friction welded joints are used in demanding mechanical systems, and their low cycle fatigue response is strongly tied to local deformation under repeated loading. By showing that average strain extracted from deformation images provides a much stronger prediction basis than raw images alone, the study points toward a practical monitoring strategy: optical measurements can be useful, but their engineering value increases when converted into mechanically meaningful strain features. For components operating under cyclic loading, this distinction matters. Surface texture, contrast, and image noise may complicate life prediction if they are treated as direct model input. In an engineering setting, the more useful route is to process deformation images through digital image correlation, extract peak strain-related information, and use those values as compact descriptors of the fatigue state. This makes the approach more compatible with inspection systems that must provide interpretable and repeatable indicators rather than opaque image-based judgments.

The integration of damage fraction adds another practical layer. By combining strain response with a measure of accumulated fatigue damage, the model connects what is observed at a given stage of loading with where the component lies in its fatigue life. This is especially relevant for non-destructive evaluation of in-service equipment, where maintenance decisions depend not only on whether deformation is occurring, but on how that deformation relates to life consumption. The demonstrated improvement in prediction accuracy suggests that data-driven fatigue assessment can benefit from physics-based descriptors when they are chosen carefully.

In design and maintenance workflows, the approach could support more informed evaluation of welded joints, fatigue-critical regions, and components subjected to controlled cyclic loading. It may help engineers compare fatigue states across specimens or service intervals using strain-based features rather than relying only on final failure data. Within the tested range, the strongest engineering message is the value of a focused data-physics framework that estimates low cycle fatigue life from experimentally accessible deformation and damage information.

Cavitation erosion results from the intense interaction between a collapsing bubble and a solid surface. A vapor-filled bubble can expand and collapse very quickly, and the resulting liquid motion can break the protective surface film, deform the material, and create localized pits. Near a rigid wall, the bubble usually collapses asymmetrically. The pressure difference caused by the wall pushes the bubble toward the surface, generates a high-velocity micro-jet, and changes the rebounding bubble into a ring-shaped structure. When this structure collapses, it can produce fragmented flow directed toward the wall. Separating the damage caused by these different stages is important for understanding single-bubble erosion. The main difficulty is both temporal and physical. High-speed imaging can capture changes in bubble shape, jet penetration, movement, rebound, and fragmentation, but it does not directly track how the material responds during these events. Surface microscopy can reveal pits and depressed regions, but these features are usually examined after the collapse sequence or after repeated impacts. A depressed area may be linked to a jet, while a ring of pits may result from toroidal bubble collapse. Surface morphology alone, however, cannot always determine the contribution of each event. Chronoamperometry offers a way to connect bubble motion with surface damage on the relevant timescale. Pure aluminum carries a thin passive oxide film. When collapse-induced loading ruptures that film and exposes the underlying metal to an oxygen-containing electrolyte, oxidation produces a transient current. The height of the resulting current peak reflects the area and depth of the disturbed surface, while repassivation returns the current toward its baseline. This electrochemical response can therefore provide a real-time measure of erosion that complements optical observations rather than replacing them.

In a recently published research paper in International Journal of Multiphase Flow,  Professor Qingmiao Ding, Professor Yanyu Cui, Dr. Yunlong Shan, Dr. Jiakui Yu, and Dr. Futai Guo from the Civil Aviation University of China developed an experimental method that combines high-speed bubble imaging with real-time chronoamperometric measurement of passivation-layer damage on pure aluminum. They introduced a controlled eccentric bubble position that diverts the first-collapse micro-jet away from the working electrode while allowing ring-vortex crushing flow to reach it.

The researchers generated cavitation bubbles with a pulsed Nd:YAG laser and positioned them at controlled distances from a polished, high-purity aluminum working electrode immersed in a saline solution. They defined the dimensionless stand-off distance, γ, as the distance between the bubble centre and the wall divided by the maximum bubble radius. Repeated measurements were used to determine the mean peak current at each position.

Bubble motion changed markedly across the examined distances. Far from the wall, the cavity remained nearly spherical through much of its first contraction, moved gradually toward the surface, and formed a penetrating micro-jet during collapse. Rebound produced a ring vortex that continued toward the boundary and fragmented during later collapse, but the resulting flows had to travel a relatively long distance before reaching the aluminum.

The authors found a different sequence appeared at γ = 1.3 and the wall distorted the contracting bubble from a spherical to an ellipsoidal form. Its first collapse generated a hemispherical jet directed toward the surface, after which the bubble became a ring vortex. During the second collapse, that toroidal structure collapsed asymmetrically and produced a crushing flow accompanied by numerous small bubbles spreading close to the wall. At γ = 0.1, the bubble was tightly adjacent to the boundary and became strongly flattened. The first-collapse jet moved radially along the surface, while the subsequent ring vortex expanded, broke apart, and drove fragmented flow against the near-wall region.

They also performed chronoamperometric measurements which showed a non-monotonic dependence of erosion on stand-off distance. The largest peak current occurred at γ = 0.1, identifying this condition as the most damaging among those tested. Damage then decreased, varied at intermediate positions, and rose again near γ = 1.3 before declining with increasing separation. At sufficiently large distances, current transients became barely detectable. Distance therefore did more than scale impact strength; it altered the collapse mode and, consequently, the mechanism by which energy reached the surface.

To separate micro-jet damage from ring-vortex damage, the team introduced a dimensionless eccentricity of ε = 0.8. This design choice displaced the bubble vertically so that the first-collapse jet struck the surrounding insulating substrate rather than the aluminum electrode, while the more spatially distributed crushing flow from the toroidal collapse could still reach the metal. At γ = 1.3, the measured peak current changed by less than 1% when eccentricity was introduced. The micro-jet thus contributed little to erosion under that condition, whereas the second-collapse ring-vortex flow accounted for most of the electrochemical response. At γ = 0.1, eccentricity reduced the peak current to less than half its centred value, showing that the micro-jet made a substantial contribution there.

The team observed repeated-bubble exposure to produce surface morphologies consistent with this interpretation. At γ = 0.1, the aluminum developed a central depressed region surrounded by an asymmetric annulus of concentrated pits. The central deformation was associated with the jet, while the annular pitting corresponded to fragmented crushing flow from the later toroidal collapse. At γ = 1.3, the centre remained comparatively flat and the damage consisted mainly of densely distributed, shallower pits in an asymmetric ring. The agreement among bubble imaging, current response, and surface morphology supported the distinction between the two erosion mechanisms.

The findings of Professor Yanyu Cui and colleagues provide a useful basis for interpreting and managing cavitation erosion in hydraulic machinery, marine equipment, piping systems, valves, pumps, propellers, and other components exposed to repeated bubble collapse near solid boundaries. A key engineering implication is that erosion severity cannot be estimated from bubble proximity alone. The dimensionless stand-off distance changes the collapse mode and determines whether surface damage is driven mainly by a direct micro-jet, by the fragmented flow produced during ring-vortex collapse, or by both mechanisms acting together. A centrally depressed region surrounded by annular pitting is consistent with the combined action of a wall-directed jet and the later collapse of a toroidal bubble. In contrast, dense ring-shaped pitting with a comparatively flat centre indicates that ring-vortex crushing flow is likely to be the dominant source of erosion. Such morphological signatures may help engineers relate observed surface damage to the underlying bubble dynamics rather than treating all cavitation pits as products of a single impact mechanism.

The non-monotonic relation between stand-off distance and damage is also relevant to equipment design. The most severe erosion occurred when the bubble was very close to the wall, but pronounced damage reappeared at an intermediate distance where the direct jet contributed little. Geometric modifications intended to redirect bubbles, alter near-wall flow, or move collapse zones away from vulnerable surfaces should account for the possibility that toroidal collapse and fragmented crushing flow can remain damaging even when direct jet impact is limited. The combined measurement approach proposed in the new study has practical value for laboratory evaluation of cavitation-resistant materials and protective surfaces. Chronoamperometry records passivation-layer disruption in real time, while high-speed imaging identifies the corresponding collapse stage and surface microscopy confirms the accumulated damage pattern.

Extrusion-based additive manufacturing gives polymer processing a form of geometric freedom that conventional molding cannot easily match and can build customized parts directly from digital toolpaths, produce internal features without dedicated molds, and adapt quickly from one geometry to another. However, the same layer-by-layer and road-by-road deposition that gives the process its flexibility also limits performance when dense, load-bearing polymer parts are needed. A printed component made by fused deposition or fused filament fabrication is not simply a shaped polymer body; it is an assembly of deposited roads whose interfaces record the thermal and mechanical history of printing. The weakest direction is often the build direction, where tensile loading must be carried through interlayer bonds rather than through a fully consolidated material volume. A second challenge is production rate and increasing extrusion speed alone cannot fully resolve it because the melt must still be generated, delivered, deposited, and bonded along a sequential path. Higher flow rates also affect nozzle pressure, material residence time, and deposition stability. The trade-off between resolution and throughput therefore remains deeply connected to the physics of material extrusion, not only to machine settings.

In a recently published research paper in Journal of Manufacturing Processes, Zexin Wang, Shilin Wang, Jiawei Zhu and Professor Ranran Jian from Qingdao University of Science and Technology working together with Dr. Xiang Cheng from China Aerospace Science and Technology Corporation developed an injection-based additive manufacturing process in which an extrusion-printed thermoplastic shell serves as an intrinsic mold for a subsequently injected core. The work is unique in the unified reciprocating screw mechanism, which performs shell printing by rotation and core filling by axial injection within the same material system. They also developed a coupled simulation strategy to analyze shell-core interfacial melting, pressure-driven deformation, cooling shrinkage, and debonding.

During shell fabrication, screw rotation controls the extrusion and deposition of successive layers. The same unit then shifts to axial motion and injects molten polypropylene into the cavity enclosed by the printed shell. This dual operation replaces the slowly deposited internal infill with a dense core formed through rapid volumetric filling. The researchers treated the new mechanism as a coupled thermal, fluid, and structural problem. They examined flow simulations and how nozzle diameter and injection speed influence shear heating during core delivery. The authors found for polypropylene, small nozzles sharply restricted the safe injection rate because the melt temperature rose too quickly. Larger nozzle diameters provided a better balance between extrusion stability and injection capacity, while very small nozzles generated higher pressure and stronger velocity fluctuations during shell printing. The nozzle used for validation was therefore selected to perform reliably in both phases of the process.

The research team paid good attention to the shell-core interface  and found that molten polypropylene entering the cavity transferred heat to the colder printed shell and partially re-melted its inner surface. This local re-melting promoted fusion across the interface rather than leaving the core in mechanical contact with the shell alone. At the same time, injection pressure could deform the shell, and uneven cooling could generate shrinkage stresses between the two regions. The co-simulation therefore considered interfacial melting, shell deformation, structural collapse, and debonding. Shell thickness and injection flow rate proved especially important. Thicker shells resisted deformation and retained greater thermal stability, whereas slower filling prolonged heat transfer into the shell.

The investigators printed tensile specimens with the build direction aligned with the loading direction. This placed the comparison in the mechanically weakest orientation for conventional extrusion-based manufacturing and made the test particularly sensitive to interlayer bonding. The reference specimens contained fully printed grid infill, whereas the hybrid specimens combined a printed shell with an injected core. The injection-based parts showed an approximately 80% increase in ultimate tensile strength and an approximately 89% increase in elongation at break. These changes indicate that the hybrid process improved both resistance to failure and the capacity for deformation before fracture. They also examined fracture morphology and noticed the extrusion-based specimens retained the layered structure and local voids associated with deposited roads, and failure occurred mainly along weak interlayer boundaries. Their fracture surfaces also showed stress whitening and localized crazing, consistent with greater plastic deformation. The improvement arose from two related effects: injection replaced much of the road-by-road internal structure with a consolidated core, and the heat of the incoming melt strengthened bonding where the core met the printed shell.

The findings of Professor Ranran Jian and colleagues point to practical engineering use wherever polymer components need the geometric flexibility of additive manufacturing but also require better consolidation, higher build-direction strength, and shorter production time than conventional extrusion-based printing can usually provide. The shell-mold printing and core-filling strategy is especially relevant for customized structural polymer parts in which the outer geometry must be digitally defined, but the interior does not need to be built slowly as a printed grid. By replacing printed infill with an injected polymeric core, the process creates a denser internal structure and reduces the number of weak interlayer interfaces that normally control failure in the build direction. This makes the approach useful for load-bearing or semi-load-bearing thermoplastic components where conventional FDM or FFF parts may be limited by anisotropy. Housings, brackets, fixtures, tooling aids, protective casings, and customized mechanical supports could benefit from the combination of printed shape control and injection-like core consolidation, provided the part can be designed with cavities that can be filled through suitable gates and venting paths. The work also suggests value for small-batch manufacturing, where traditional injection molding may be too costly or slow because a dedicated mold is required. Here, the printed shell becomes the mold, allowing mold fabrication and part fabrication to merge into one additive route.

In the authors’ analysis, the time saved becomes more substantial as the volume of printed infill that would otherwise be deposited line by line is replaced by rapid cavity injection. This is important for industrial additive manufacturing, where production time often determines whether a process is practical beyond prototyping. From a design standpoint, we can think of the new study gives engineers useful guidance  and shell thickness, nozzle diameter, injection speed, cavity height, and filling pressure must be chosen together because they control thermal penetration, shell deformation, interfacial fusion, and shrinkage behavior. Controlled shell deformation may even help compensate for cooling shrinkage, but excessive pressure or heat must be avoided.

Copper is an important engineering material because it combines efficient electrical and thermal transport with useful structural integrity. Its high electrical and thermal conductivity make it suitable for components that carry current, remove heat, or operate in thermally demanding assemblies. At the same time, many such components increasingly require geometries that are difficult to produce by conventional subtractive or forming methods. Laser powder bed fusion offers an attractive manufacturing approach in this setting because it can build dense metal parts directly from powder while allowing considerable geometric freedom. The main challenge comes from the interaction between copper and the near-infrared lasers commonly used in commercial LPBF systems. Copper reflects much of the incident radiation at these wavelengths, so only a limited fraction of the applied energy is absorbed by the powder bed. Also the high thermal conductivity of copper can bring another challenge because heat spreads rapidly away from the irradiated region, which make it harder to maintain a stable melt pool and reliable fusion between adjacent scan tracks and layers. The result is a narrow processing range in which relatively small changes in laser power or scan speed can shift the material toward lack-of-fusion defects, porosity, spheroidization, or unstable keyhole behaviour.

Several approaches have been explored to improve the LPBF processing of copper. Higher-power near-infrared lasers can increase energy delivery, while shorter-wavelength laser sources can improve absorption. These solutions, however, may require costly equipment or introduce practical limitations related to system capability and laser–material interaction. A different route is to modify the powder itself. Alloying additions can improve printability, but elements that dissolve substantially in copper may also impair electrical conductivity through electron scattering. For applications in which copper is selected primarily for its transport properties, this trade-off is difficult to ignore.

Particle-based modification provides another possibility. A dispersed additive can alter the optical and thermal response of the powder bed while remaining distinct from the copper matrix to a meaningful extent. LaB₆ is of particular interest because of its high absorptivity under near-infrared irradiation and its potential to remain as a second phase in copper. Earlier studies have considered nanoscale LaB₆ additions, but nanoparticles are prone to agglomeration and may reduce powder flowability, complicating powder spreading during LPBF. Their preparation can also be costly.

Micrometre-scale LaB₆ particles offer a potentially more practical alternative, yet their behaviour during copper LPBF has remained insufficiently understood. Whether such particles can improve laser absorption and densification while also influencing microstructural development, mechanical response, and electrical or thermal transport requires direct examination. The present investigation addresses that question by introducing LaB₆ microparticles into copper powder and examining how they affect the printing process and the resulting material.

In a recently published research paper in Journal of Manufacturing Processes Dr. Jinchao Zhao, Professor  Jiabin Liu, Dr.  Liuyi Huang, Professor  Kai Ren, and Professor Yanlong Cao from Zhejiang University developed a 1 wt% LaB₆ microparticle-modified copper feedstock for near-infrared laser powder bed fusion. The powder preparation combined ultrasonic vibration mixing with three-dimensional mixing to distribute LaB₆ across copper particles. During printing, the material produced a dual-scale LaB₆ structure composed of retained larger particles and in situ precipitated nanoscale LaB₆ particles. This combination increased near-infrared absorptivity, broadened the high-density processing window, and strengthened the printed copper without eliminating its high electrical and thermal conductivity.

The researchers prepared the composite feedstock by combining gas-atomized copper powder with irregular LaB₆ microparticles using ultrasonic vibration mixing followed by three-dimensional mixing under argon protection. Elemental mapping showed that the LaB6 particles were distributed across the copper-powder surfaces rather than concentrated in isolated regions. This dispersion was important because the intended effect depended on LaB6 being available throughout the powder bed to interact with incident laser energy. They found that at 1080 nm, adding LaB₆ microparticles increased powder absorbance by 50.4%, changing the densification behaviour during LPBF. Across a substantially broader range of processing conditions than pure copper, the composite achieved relative densities above 99.9%, whereas unmodified copper remained more susceptible to unfused regions and porosity when laser input was insufficient.

The benefit was therefore not confined to a single optimized condition. The authors found that the high-density processing window of the composite was approximately twice as wide as that of pure copper, which indicates greater tolerance to parameter variation during printing. The composite also developed larger melt-pool widths and depths under comparable conditions. By improving laser-energy coupling within the powder bed, LaB6 promoted more stable melting and supported fusion between neighbouring tracks and successive layers.

The team performed microstructural analysis which showed that the initial microparticle population did not behave uniformly during LPBF. Larger LaB₆ particles remained partially or fully unmelted in some regions and were redistributed within the liquid copper by melt-pool flow before becoming embedded on solidification. Their number decreased as volumetric energy density increased, consistent with greater particle melting at higher energy input. At the same time, transmission electron microscopy identified a second LaB6 population: regularly shaped nanoparticles with an average size of approximately 60 nm dispersed in the copper matrix.

This nanoscale population was interpreted as the result of complete melting of smaller LaB₆ particles, dissolution into the liquid copper, and subsequent precipitation during rapid solidification. La and B became supersaturated in the melt and reprecipitated as LaB₆ nanoparticles within grains and at grain boundaries. They found using X-ray diffraction no measurable shift in copper lattice parameters or the formation of new phases, supporting the interpretation that LaB6 remained a second phase rather than producing a substantially altered copper solid solution. The team found that both pure copper and composite samples retained a columnar grain structure aligned with the build direction and the composite showed weaker texture, a higher fraction of low-angle grain boundaries, and greater kernel average misorientation, consistent with increased lattice distortion and dislocation density. These changes were reflected in the tensile response: the LaB₆/Cu composite reached a yield strength of 260 MPa and an ultimate tensile strength of 385 MPa, compared with 127 MPa and 246.5 MPa for LPBF pure copper. Its elongation decreased from 34% to 25.6%, although the composite remained appreciably ductile.

Overall, the findings of Zhejiang University researchers are relevant to copper components that must retain high electrical and thermal transport while being manufactured in geometries that are difficult to obtain through conventional routes. The composite retained 87.4% IACS electrical conductivity and a calculated thermal conductivity of 359 W·m⁻¹·K⁻¹ while providing substantially higher strength than the LPBF pure-copper material.  Electrical components with complex internal or external geometries are one clear area of relevance. Copper is widely used where current must be carried efficiently, but the present results indicate that a modest LaB₆ addition can preserve much of that electrical capability while improving the mechanical response of LPBF-built material.

Thermal-management hardware is another direct application area identified by the paper’s starting context. The retained thermal conductivity supports the use of the material where heat must be transferred through a printed copper body, while LPBF provides access to geometrically intricate forms. Such flexibility matters for designs that depend on internal passages, local thickness variations, integrated mounting features, or compact heat-transfer pathways that would be difficult to machine from bulk copper. The present study does not demonstrate a specific thermal-management device, but it establishes a materials-processing route compatible with that class of component.

The manufacturing implications are equally important and by introducing 1 wt% LaB₆ increased powder absorptivity by 50.4% and expanded the high-density parameter range by more than 1.8 times relative to pure copper. This broader process window gives engineers greater latitude when selecting laser power and scan speed, and it can reduce the sensitivity of production to small parameter variations.

Microparticle-based feedstock modification may also be attractive where powder preparation must remain practical. The paper emphasizes that micrometre-scale particles distribute more readily than nanoparticles and can reduce preparation complexity and cost. Their behaviour during printing is especially useful: larger LaB₆ particles can remain dispersed in the copper matrix, whereas smaller particles melt and reprecipitate as nanoscale LaB₆ during solidification. The resulting dual-scale structure links processability with strengthening, while retaining the conductive character required for electrically and thermally functional copper parts.

Compact Strip Production rolling mills belong to a class of large industrial systems in which mechanical transmission, hydraulic actuation, electrical drive control, and plastic deformation of metal operate as one physically coupled process. The roll system must sustain large deformation resistance from the strip, the drive train must deliver torque under changing load conditions, and the hydraulic screw-down system must maintain the roll gap while responding to force fluctuations generated during rolling. When these actions occur during the production of high-strength thin strip, even small dynamic disturbances can be amplified through the mill structure and through the deformation zone itself. A recurring difficulty in CSP rolling is the appearance of strong vibration during rolling operations, particularly when demanding strip grades and thin specifications are produced. Such vibration is not only a structural response of the mill stand or a torsional response of the drive shaft. It can be accompanied by changes in rolling force, motor torque current, hydraulic pressure, roll displacement, interface friction, and strip surface quality. The scientific problem is therefore not simply to identify a vibrating component, but to understand how different subsystems exchange dynamic influence during rolling.   The deformation zone, where the roll and strip interact under pressure, is the location where torque, vertical force, friction stress, roll gap variation, and material deformation are converted into one another. For modeling clarity, the main drive transmission system and the vertical hydraulic screw-down system are often treated as separate dynamic subsystems. That separation is useful and, in many cases, necessary for establishing solvable mathematical models, but the present rolling condition requires attention to what occurs where both subsystems act on the strip. Torsional fluctuation in the roll can change the circumferential velocity at the interface, thereby modifying friction stress and rolling force. Vertical vibration can change the roll gap and contact arc length, thereby altering rolling load torque and feeding disturbance back into the transmission system. The rolling deformation zone is therefore the site of strip reduction and also a dynamic coupling region.

In a recently published research paper in International Journal of Precision Engineering and Manufacturing Associate Professor Yifang Zhang, Dr. Yiwei Wan, Dr. Huiwei Yan,  Professor Cheng He, Professor Li Cui, Dr. Xu Ding, Dr. Tianyi Chen & Dr. Pingye Wan from Shanghai Polytechnic University developed a vertical-torsional-deformation zone coupling closed-loop dynamics model for a Compact Strip Production rolling mill. The model links the main drive transmission system, hydraulic screw-down system, and rolling deformation zone through parameter pathways involving roll torsion angle, circumferential velocity, interface friction stress, rolling force, contact arc length, and rolling load torque. They also developed and applied a synchronized industrial monitoring approach that allowed the proposed closed-loop mechanism to be checked against field vibration signals and lubrication-adjustment experiments.

The research team began with industrial monitoring on the F3 mill during rolling operations. They developed a measurement system capable of synchronously collecting torsional and vertical vibration information, while also drawing motor current and hydraulic rolling force data from the plant’s online process data acquisition system. Torque and bending moment in the drive system were measured through resistance strain gauges mounted on the rotating shaft, while vertical roll motion was captured using acceleration and displacement sensors installed near the roll bearing seat. This synchronized acquisition was important because the coupling mechanism depends on frequency relationships among signals measured under the same rolling condition.

The field signals revealed a structured frequency relationship between the drive and vertical systems. Torsional vibration in the main drive showed a dominant component near 41 Hz. The vertical vibration contained an 82 Hz dominant component as well as a 41 Hz component, linking the vertical response to the same fundamental frequency present in the torsional response. The hydraulic rolling force and motor torque current signals also contained harmonic components near the same 41 Hz range observed in the torsional and vertical vibration responses and it is this recurrence across drive, hydraulic, and vertical measurements gave the frequency-domain evidence a clear physical coherence, showing harmonically related behavior across subsystems. The team then calculated the inherent torsional characteristics of the F3 rolling mill using a finite element model built from the mill structure. The second-order torsional modal frequency was 42.4 Hz, close to the approximately 41 Hz component observed in the monitored vibration. External disturbances and harmonic combinations generated during rolling could approach the natural torsional frequency of the system, allowing strong vibration to develop.  For the main drive transmission, the researchers used an equivalent two-inertia nonlinear torsional model with excitation from both the electrical drive end and the roll load end. Nonlinear stiffness and damping terms were retained, and the authors used regular perturbation method to obtain the angular response under multi-source harmonic excitation. The solution showed that the torsional response contains frequency components generated by combinations of the excitation frequencies. When one of these combined components approaches the inherent modal frequency, strong torsional vibration can be induced.

That torsional motion was then connected to the rolling deformation zone. Fluctuation in the roller torsion angle changes the circumferential velocity of the roll. Under mixed lubrication conditions, this velocity change affects the tangential friction stress at the rolling interface. Since the total rolling force depends on both deformation resistance and friction stress, torsional vibration can produce rolling force fluctuation, which then excites the vertical roll system. The authors modeled the reverse pathway through the hydraulic screw-down system, represented as an equivalent two-mass system involving the roller and hydraulic cylinder under hydraulic pressure and rolling force fluctuations. Solving the vertical vibration response showed that vertical roll displacement changes the roll gap and contact arc length. Changes in contact arc length alter the rolling load torque, feeding excitation back into the transmission system. In mechanical terms, vertical vibration changes the torque demand placed on the drive; torsional vibration changes the frictional and force conditions imposed on the vertical system.

The authors’ lubrication experiment provided a practical validation of this closed-loop interpretation and by adjusting the lubrication oil supply at the roll gap, the researchers altered a parameter located directly in the deformation-zone coupling path. When lubrication oil content was within an appropriate range, vertical and torsional vibration amplitudes were lower than at either insufficient or excessive lubrication levels. In the comparison before and after lubrication adjustment, the torsional vibration amplitude decreased from 39.68 to 8.35, while the vertical vibration amplitude decreased from 40.01 to 18.52. The dominant vibration shifted away from the strong coupled 41 Hz condition, and the vertical response no longer showed excitation at that critical frequency. A change at the roll-gap interface therefore changed the vibration state of the coupled mill system.

The findings of Associate Professor Yifang Zhang  et al have direct engineering relevance for the diagnosis and suppression of vibration in Compact Strip Production rolling mills, especially during the rolling of high-strength thin strip. The main practical value is the shift from treating torsional vibration, vertical vibration, hydraulic force fluctuation, and motor torque current as separate symptoms toward treating them as linked expressions of one coupled dynamic system. For mill engineers, this means that vibration control should not begin only with the drive train, the hydraulic screw-down system, or the mill stand in isolation.

One important application is in industrial monitoring and the new study shows that meaningful diagnosis requires synchronized measurement of signals from the electrical drive, mechanical transmission, hydraulic screw-down system, and roll system. When related frequency components appear across these signals, especially near the inherent modal frequency of the mill, they can indicate a coupled vibration condition rather than a local disturbance. This provides a practical basis for condition monitoring systems that do more than record vibration amplitude. They can track frequency relationships among motor current, rolling force, torque, roll displacement, and acceleration, allowing operators to identify when the mill is approaching a strongly coupled vibration state.

The results also support more targeted vibration suppression strategies. Since the rolling deformation zone transmits disturbances between torsional and vertical motion, engineering adjustments at the roll gap can influence the behavior of the entire system. Adjusting the lubrication oil supply within an appropriate range reduced the coupling degree of the system and substantially lowered both torsional and vertical vibration amplitudes. This suggests that lubrication control is not only a tribological or surface-quality measure, but also a dynamic control parameter for the rolling mill. The new model can also guide process optimization for high-strength thin strip production and by linking roll torsion angle, circumferential velocity, interface friction stress, rolling force, contact arc length, and rolling load torque, it provides engineers with a structured way to understand how changes in operating conditions may feed back through the mill. Such a model can help engineers define more stable operating windows, avoid excitation near critical modal frequencies, and improve rolling stability through process-level control as well as equipment-level measures.

Liquid-loaded membranes represent a familiar yet mechanically demanding class of fluid–structure interaction problems. A thin membrane held under initial tension may appear, at first sight, to respond to liquid self-weight in a straightforward way: the liquid applies pressure, the membrane deflects, and a larger liquid volume should produce a deeper deformation. This expectation is natural because conventional membrane models often treat the liquid load as a prescribed pressure acting on the undeformed configuration. Under that view, the liquid level is effectively fixed by the initial volume, and the membrane’s response follows from the balance between pressure and tension. The difficulty arises when the membrane deformation is no longer negligible relative to the available liquid depth. As the membrane deflects, it creates additional space that must be filled by part of the liquid. The liquid level therefore changes after deformation, and the hydrostatic pressure acting on the membrane is no longer determined only by the initial height. The load follows the evolving configuration. This seemingly small distinction becomes central, because the membrane shape, the liquid level, and the region of liquid contact are all coupled through volume conservation. A model that ignores this coupling may still perform well for small deflections, high membrane tension, short spans, or low-density liquids, but it cannot explain the full range of observed membrane–liquid configurations. In a recent research paper published in International Journal of Engineering Science, Dr. Weiting Chen and Professor Quanzi Yuan from the University of Chinese Academy of Sciences developed an analytical model for liquid-loaded initially stressed membranes in which the hydrostatic pressure is determined by the deformed liquid configuration rather than by the initial liquid height. They derived closed-form solutions for one-dimensional and two-dimensional axisymmetric membrane deflections, including membrane-prevails, equipoise, and liquid-prevails regimes.  This allowed them to identify a dimensionless control parameter that determines the membrane–liquid configuration independently of liquid volume.

The researchers treated the membrane as initially tensioned, linearly elastic, and thin enough that bending stiffness could be neglected. Rather than prescribing the liquid pressure from the initial liquid height, they enforced liquid-volume conservation after deformation.   For the one-dimensional membrane, the classical model gives a parabolic deflection controlled by the dimensionless parameter λ = √(ρgL²/T). In the present model, the same parameter acquires direct physical meaning as the measure of competition between gravity-driven loading and membrane tension. When 0 < λ < π, the membrane prevails, and the liquid remains in contact with the full membrane span. The deformed liquid level decreases as λ increases, following an analytical expression derived from the volume constraint. At λ = π, the system reaches the equipoise configuration, where the deformed liquid level is exactly flush with the bottom. For λ > π, the liquid prevails, and the interaction region shrinks. In that regime, the boundary of contact is determined analytically, and its location depends on λ but not on the original liquid volume.

Chen and Yuan used a PET membrane clamped in a tension apparatus with a transparent tank, varying the membrane tension to control λ. The experiments reproduced the three predicted configurations. The membrane prevailed for λ below π, equipoise occurred at λ = π, and liquid prevailed when λ exceeded π. The measured liquid levels and interaction regions agreed closely with the analytical predictions, while the classical model did not reproduce the observed regime selection with the same accuracy. The authors also compared their formulation with a nonlinear classical model for an inextensible membrane under uniform pressure. Including geometric nonlinearity changed the predicted deflection, especially as λ and the initial liquid level increased, but it did not resolve the central discrepancy. The reason is instructive: curvature and tension variation matter, yet they do not replace the need to update the liquid loading according to the deformed configuration.

The team extended same reasoning to two-dimensional axisymmetric membranes. With ξ = √(ρgR²/T), the axisymmetric formulation yields analytical solutions involving Bessel functions. The transition occurs at the first nontrivial zero of J0, approximately 2.4. Below this value, the membrane prevails; at the critical value, equipoise is reached; above it, the liquid prevails and the liquid-contact radius decreases according to the dimensionless parameter. This extension shows that the volume-independent regime selection is not a peculiarity of the one-dimensional geometry.

The findings are directly useful for engineering systems in which a thin, tensioned membrane supports or confines a liquid but the liquid level is free to adjust during deformation. A practical example is flexible covering films exposed to rainwater. In such systems, designers often estimate deformation from the applied water depth or total accumulated volume. Chen and Yuan’s analysis indicates that this can be misleading when the membrane deflection significantly alters the water configuration. The relevant design question becomes whether membrane tension, span, and liquid density place the system in a membrane-prevailing, balanced, or liquid-prevailing regime. This gives engineers a clearer criterion for deciding when ponding remains shallow and distributed, and when liquid contact may shrink into a central region with much larger local deflection.

In these cases, the membrane is not simply a passive surface under a fixed pressure; it reshapes the liquid domain. The dimensionless parameter identified in the paper provides a compact way to tune design variables. Increasing initial membrane tension or reducing the characteristic span lowers the competition parameter, keeping the system closer to the classical small-deflection regime. Larger spans, lower tensions, or denser liquids push the system toward stronger configuration-dependent behavior. This is useful because the engineer can adjust geometry or prestress before relying on more complex numerical simulations. The work shows when classical uniform-pressure membrane theory is likely sufficient and when it is not. When this dimensionless parameter is small, the classical model and the configuration-dependent model give nearly identical predictions. For larger values, especially beyond the critical regime, classical predictions may miss the magnitude of deflection as well as the actual liquid–membrane contact area. That distinction matters for load paths, seal design, drainage planning, support spacing, and failure-risk assessment. The findings of Chen and Yuan may also guide experimental design for soft membranes, flexible electronics, biological-mimetic membranes, and liquid blister systems where researchers need to distinguish material effects from loading-configuration effects. By separating the role of liquid volume from the role of density, length scale, and tension, the model offers a more reliable way to interpret membrane deformation tests and its value in that it identifies the mechanical control parameter that must be respected before such models can be meaningfully applied.

Frictional mechanical structures present a demanding design problem because their performance depends on stiffness and load transfer as well as on surface degradation, heat generation, and irreversible energy dissipation at moving interfaces. Components exposed to sliding contact must often satisfy competing requirements: they must remain mechanically efficient, use material economically, and resist wear in localized surface regions where contact pressure, frictional work, and temperature rise are concentrated. Conventional topology optimization has become a mature tool for distributing material within a design domain, yet the problems most commonly addressed are compliance, frequency response, heat transfer, and stress.  A structural layout that is optimal for stiffness may not be optimal for a frictional surface whose degradation is driven by local energy dissipation. At the same time, direct treatment of wear in topology optimization is complicated by the need for a physically meaningful scalar measure that can be differentiated with respect to design variables and used to guide the optimization. In a recently published research paper in Thin-Walled Structures Dr. Guikai Guo, Jie Lei, Jichang Yang, and Professor Huanhuan Gao from Jilin University developed a SIMP-based thermoelastic topology optimization method that minimizes structural compliance while limiting material use and wear-related entropy generation at the sliding surface. Its technically distinct element is the use of the degradation entropy generation theorem to represent mechanical wear rate through frictional entropy production. They derived finite element formulations for steady-state heat conduction and static thermoelastic deformation, together with adjoint sensitivities for compliance and entropy generation rate. The new method was tested on cantilever, half MBB, and L-shaped structures under frictional heating and wear-surface loading conditions. The research team used steady-state heat conduction to compute the temperature field produced by frictional heat input, while static thermoelastic analysis accounts for the deformation and stiffness response under mechanical loading and thermally induced loads. Within the finite element framework, the global thermal conductivity matrix, thermoelastic coupling matrix, and structural stiffness matrix are assembled from element contributions. This construction matters because the design variables alter material density, and through density interpolation they also alter thermal conductivity, Young’s modulus, and the thermal stress coefficient.

The authors introduced wear through entropy generation during sliding. For the frictional contact process considered in their study, the dominant entropy contribution is taken to arise from plastic deformation work converted into heat. Under the stated assumptions of steady conditions, localized entropy generation, negligible energy transport by material loss, and complete conversion of frictional work into thermal energy inside the control volume, the entropy generation rate at the wear surface becomes a function of friction coefficient, normal force, sliding velocity, and contact temperature. With friction coefficient, normal pressure, and sliding velocity fixed in the instantaneous optimization problem, the design sensitivity enters mainly through the temperature field. A change in material distribution changes heat conduction, contact temperature, and ultimately the entropy generation rate. The authors used the SIMP material interpolation scheme and solve the resulting gradient-based problem using the method of moving asymptotes. Adjoint sensitivity analysis is derived for both the compliance objective and the entropy generation constraint, avoiding direct computation of displacement and temperature derivatives for every design variable and by formulating the entropy constraint in adjoint form, the thermodynamic wear measure becomes computationally compatible with density-based topology optimization.

The team used three benchmark structures to examine the behavior of the method: a cantilever beam, a half MBB beam, and an L-shaped beam. In the cantilever case, the wear region is placed near the upper right surface, with frictional force and heat generated from the imposed normal force and sliding velocity. When the entropy generation constraint is tightened, the optimized layout changes its material distribution near the wear region and alters the temperature field. The study reports that enforcing the constraint reduces the entropy generation rate associated with wear, while compliance increases only modestly. The design response is not a simple addition of material at the hot surface; rather, the optimizer redistributes material in a way that manages thermal transfer and local stiffness simultaneously. Additionally, they noticed as the entropy generation rate limit decreases, the topology becomes more elaborate, with added supporting members and modified internal branching. Material distribution is refined in regions where the thermal and mechanical fields interact, and the lower right region uses less material to suppress excessive heat transfer and entropy production.  There, material tends to concentrate more in vertical members, while stricter entropy constraints increase low-temperature regions and reduce high-temperature regions near the wear-related thermal field. The optimization introduces additional voids near the wear surface, adjusts supporting members near the fixed boundary, and creates load-bearing paths that reduce deformation in the wear region. They also found that increasing the filtering radius produces simpler and more uniform structures, lowers entropy production in the examined cantilever model, and raises compliance. Density penalization factors influence the clarity and continuity of the optimized layout, especially because thermal conductivity, elastic modulus, and thermal stress coefficient are interpolated separately. The friction coefficient has a direct role in both frictional force and heat source intensity; as it rises, entropy production and compliance increase, and the topology shifts toward the friction surface. Wear region size also changes the design response: when the wear surface expands in the half MBB beam, the heat-generation and load-application regions expand as well, internal trusses decrease, and the branch structure shifts toward the enlarged wear area.

The findings of Professor Huanhuan Gao  and colleagues are directly relevant to the design of mechanical components whose service performance is limited by sliding contact, frictional heating, and progressive wear. In such components, structural stiffness alone is not a sufficient design target, because the same material layout that provides efficient load transfer may also intensify thermal accumulation or entropy generation near the contact surface. The entropy-constrained topology optimization approach developed in this work offers engineers a way to treat wear resistance as an active design requirement rather than as a later material-selection or surface-treatment issue. One immediate application lies in frictional machine elements such as gears, brake pads, cams, pistons, and sliding bearings, where contact surfaces experience repeated mechanical loading and heat generation. The proposed method can guide material distribution around wear-prone regions so that stiffness is retained while entropy generation at the sliding interface is reduced. This is particularly valuable for components in which local surface degradation can shorten service life, alter contact geometry, or reduce mechanical reliability over time.

The approach also has practical value for lightweight structural design. Conventional topology optimization often removes material to satisfy a volume constraint while maintaining stiffness, but frictional components require a more careful balance. By incorporating entropy generation rate as a wear-related constraint, designers can identify layouts that do not simply minimize compliance, but also regulate temperature fields and thermal pathways associated with frictional work. The numerical examples show that this can be achieved with only a modest increase in compliance, suggesting that wear-aware stiffness optimization can be integrated into early-stage engineering design without abandoning structural efficiency. Another important application is in simulation-driven design of structures exposed to both frictional heating and mechanical loading. The method provides a finite-element framework in which heat conduction, thermoelastic deformation, frictional loading, and degradation-related entropy production are evaluated together. This allows engineers to compare design alternatives under controlled assumptions before manufacturing or testing physical prototypes. Parameter studies on filtering radius, interpolation factors, friction coefficient, and wear-region size further show how operating conditions and modeling choices influence optimized layouts. Finally, for engineering practice, the main value is that wear resistance can be brought into the topology optimization stage itself. Components exposed to sliding contact can therefore be shaped for stiffness, material economy, and reduced entropy generation at critical wear surfaces.

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.