Modeling Method of C/C-ZrC Composites and Prediction of Equivalent Thermal Conductivity Tensor Based on Asymptotic Homogenization

This article proposes a modeling method for C/C-ZrC composite materials.According to the superposition of Gaussian random field,the original gray model is obtained,and the threshold segmentation method is used to generate the C-ZrC inclusion model.Finally,the fiber structure is added to construct the microstructure of the three-phase plain weave composite.The reconstructed inclusions can meet the randomness of the shape and have a uniform distribution.Using an algorithm based on asymptotic homogenization and finite element method,the equivalent thermal conductivity prediction of the microstructure finite element model was carried out,and the influence of component volume fraction on material thermal properties was explored.The sensitivity of model parameters was studied,including the size,mesh sensitivity,Gaussian complexity,and correlation length of the RVE model,and the optimal calculation model was selected.The results indicate that the volume fraction of the inclusion phase has a significant impact on the equivalent thermal conductivity of the material.As the volume fraction of carbon fiber and ZrC increases,the equivalent thermal conductivity tensor gradually decreases.This model can be used to explore the impact of materialmicrostructure on the results,and numerical simulations have studied the relationship between structure and performance,providing the possibility of designing microstructure based on performance.

Theoretical and experimental investigation of the resonance responses and chaotic dynamics of a bistable laminated composite shell in the dynamic snap-through mode

The dynamic model of a bistable laminated composite shell simply supported by four corners is further developed to investigate the resonance responses and chaotic behaviors.The existence of the 1:1 resonance relationship between two order vibration modes of the system is verified.The resonance response of this class of bistable structures in the dynamic snap-through mode is investigated,and the four-dimensional(4D)nonlinear modulation equations are derived based on the 1:1 internal resonance relationship by means of the multiple scales method.The Hopf bifurcation and instability interval of the amplitude frequency and force amplitude curves are analyzed.The discussion focuses on investigating the effects of key parameters,e.g.,excitation amplitude,damping coefficient,and detuning parameters,on the resonance responses.The numerical simulations show that the foundation excitation and the degree of coupling between the vibration modes exert a substantial effect on the chaotic dynamics of the system.Furthermore,the significant motions under particular excitation conditions are visualized by bifurcation diagrams,time histories,phase portraits,three-dimensional(3D)phase portraits,and Poincare maps.Finally,the vibration experiment is carried out to study the amplitude frequency responses and bifurcation characteristics for the bistable laminated composite shell,yielding results that are qualitatively consistent with the theoretical results.

The First Five Years of a Phase Theory for Complex Systems and Networks

In this paper,we review the development of a phase theory for systems and networks in its first five years,represented by a trilogy:Matrix phases and their properties;The MIMO LTI system phase response,its physical interpretations,the small phase theorem,and the sectored real lemma;The synchronization of a multi-agent network using phase alignment.Towards the end,we also summarize a list of ongoing research on the phase theory and speculate what will happen in the next five years.

Active tuning of anisotropic phonon polaritons in natural van der Waals crystals with negative permittivity substrates and its application in energy transport

Phonon polaritons(PhPs)exhibit directional in-plane propagation and ultralow losses in van der Waals(vdW)crystals,offering new possibilities for controlling the flow of light at the nanoscale.However,these PhPs,including their directional propagation,are inherently determined by the anisotropic crystal structure of the host materials.Although in-plane anisotropic PhPs can be manipulated by twisting engineering,such as twisting individual vdW slabs,dynamically adjusting their propagation presents a significant challenge.The limited application of the twisted bilayer structure in bare films further restricts its usage.In this study,we present a technique in which anisotropic PhPs supported by bare biaxial vdW slabs can be actively tuned by modifying their local dielectric environment.Excitingly,we predict that the iso-frequency contour of PhPs can be reoriented to enable propagation along forbidden directions when the crystal is placed on a substrate with a moderate negative permittivity.Besides,we systematically investigate the impact of polaritonic coupling on near-field radiative heat transfer(NFRHT)between heterostructures integrated with different substrates that have negative permittivity.Our main findings reveal that through the analysis of dispersion contour and photon transmission coefficient,the excitation and reorientation of the fundamental mode facilitate increased photon tunneling,thereby enhancing heat transfer between heterostructures.Conversely,the annihilation of the fundamental mode hinders heat transfer.Furthermore,we find the enhancement or suppression of radiative energy transport depends on the relative magnitude of the slab thickness and the vacuum gap width.Finally,the effect of negative permittivity substrates on NFRHT along the[001]crystalline direction ofα-MoO3 is considered.The spectral band where the excited fundamental mode resulting from the negative permittivity substrates is shifted to the first Reststrahlen Band(RB 1)ofα-MoO_(3) and is widened,resulting in more significant enhancement of heat flux from RB 1.We anticipate our results will motivate new direction for dynamical tunability of the PhPs in photonic devices.

Theoretical analysis of surface waves in piezoelectric medium with periodic shunting circuits

The investigations of surface waves in the piezoelectric medium bring out great possibility in designing smart surface acoustic wave(SAW)devices.It is important to study the dispersion properties and manipulation mechanism of surface waves in the semi-infinite piezoelectric medium connected with periodic arrangement of shunting circuits.In this study,the extended Stroh formalism is developed to theoretically analyze the dispersion relations of surface waves under different external circuits.The band structures of both the Rayleigh wave and the Bleustein-Gulyaev(BG)wave can be determined and manipulated with proper electrical boundary conditions.Furthermore,the electromechanical coupling effects on the band structures of surface waves are discussed to figure out the manipulation mechanism of adjusting electric circuit.The results indicate that the proposed method can explain the propagation behaviors of surface waves under the periodic electrical boundary conditions,and can provide an important theoretical guidance for designing novel SAW devices and exploring extensive applications in practice.

Numerical simulation of the mechanical behavior of superconducting tape in conductor on round core cable using the cohesive zone model

Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.

Controllable rectification on the thermal conductivity of porous YBa_(2)Cu_(3)O_(7−x) superconductors from 3D-printing

Superconducting YBa_(2)Cu_(3)O_(7−x)(YBCO)bulks have promising applications in quasi-permanent magnets,levitation,etc.Recently,a new way of fabricating porous YBCO bulks,named direct-ink-writing(DIW)3D-printing method,has been reported.In this method,the customized precursor paste and programmable shape are two main advantages.Here,we have put forward a new way to customize the YBCO 3D-printing precursor paste which is doped with Al_(2)O_(3)nanoparticles to obtain YBCO with higher thermal conductivity.The great rheological properties of precursor paste after being doped with Al_(2)O_(3)nanoparticles can help the macroscopic YBCO samples with high thermal conductivity fabricated stably with high crystalline and lightweight properties.Test results show that the peak thermal conductivity of Al_(2)O_(3)-doped YBCO can reach twice as much as pure YBCO,which makes a great effort to reduce the quench propagation speed.Based on the microstructure analysis,one can find that the thermal conductivity of Al_(2)O_(3)-doped YBCO has been determined by its components and microstructures.In addition,a macroscopic theoretical model has been proposed to assess the thermal conductivity of different microstructures,whose calculated results take good agreement with the experimental results.Meanwhile,a microstructure with high thermal conductivity has been found.Finally,a macroscopic YBCO bulk with the presented high thermal conductivity microstructure has been fabricated by the Al_(2)O_(3)-doped method.Compared with YBCO fabricated by the traditional 3D-printed,the Al_(2)O_(3)-doped structural YBCO bulks present excellent heat transfer performances.Our customized design of 3D-printing precursor pastes and novel concept of structural design for enhancing the thermal conductivity of YBCO superconducting material can be widely used in other DIW 3D-printing materials.

Editorial: Computational mechanics of granular materials

Most of granular materials are highly heteroge- neous, composed of voids and particles with different sizes and shapes. Geological matter, soil and clay in nature, geo-structure, concrete, etc. are practical ex- amples among them. From the microscopic view, a lo- cal region in the medium is occupied by particles with small but finite sizes and granular material is naturally modeled as an assembly of discrete particles in contacts On the other hand, the local region is identified with a material point in the overall structure and this discon- tinuous medium can then be represented by an effective continuum on the macroscopic level

Editorial:Mechanics of intelligent materials and structures

Recently, intelligent or smart materials and structures have been received more and more attention due to their distinguished multi-field coupling properties and wide applications in aerospace, automobiles, civil structures, medical devices, information storage, energy harvesting and so on. It is of academic challenge to fully understand the complex multi-field coupling behaviors of various smart materials and structures, and of engineering sig- nificance to enhance the performance and reliability of these materials and structures in industrial applications. The papers in the special topic of Mechanics of Intelligent Materials and Structures focus on the understanding of the electromechanical, magneto-elastic, and magneto-rheological coupling behav- iors and properties of smart materials and structures for applications in vibration control, resonators, and various functional devices.

Slip boundary effect on the critical Reynolds number of subcritical transition in channel flow

In this letter,the effect of slip boundary on the origin of subcritical transition in two-dimensional chan-nel flows is studied numerically and theoretically.It is shown that both the positive and the negative slip lengths will increase the critical Reynolds number of localized wave packet and hence postpone the transition.By applying a variable transformation and expanding the variables about a small slip length,it is illustrated that the slip boundary effect only exists in the second and higher order modulations of the no-slip solution,and hence explains the power law found in simulations,i.e.the relative increment of the critical Reynolds number due to the slip boundary is proportional to the square of the slip length.

Mathematical modeling of mixed convective MHD Falkner-Skan squeezed Sutterby multiphase flow with non-Fourier heat flux theory and porosity

In a wide variety of mechanical and industrial applications,e.g.,space cooling,nuclear reactor cooling,medicinal utilizations(magnetic drug targeting),energy generation,and heat conduction in tissues,the heat transfer phenomenon is involved.Fourier’s law of heat conduction has been used as the foundation for predicting the heat transfer behavior in a variety of real-world contexts.This model’s production of a parabolic energy expression,which means that an initial disturbance would immediately affect the system under investigation,is one of its main drawbacks.Therefore,numerous researchers worked on such problem to resolve this issue.At last,this problem was resolved by Cattaneo by adding relaxation time for heat flux in Fourier’s law,which was defined as the time required to establish steady heat conduction once a temperature gradient is imposed.Christov offered a material invariant version of Cattaneo’s model by taking into account the upper-connected derivative of the Oldroyd model.Nowadays,both models are combinedly known as the Cattaneo-Christov(CC)model.In this attempt,the mixed convective MHD Falkner-Skan Sutterby nanofluid flow is addressed towards a wedge surface in the presence of the variable external magnetic field.The CC model is incorporated instead of Fourier’s law for the examination of heat transfer features in the energy expression.A two-phase nanofluid model is utilized for the implementation of nano-concept.The nonlinear system of equations is tackled through the bvp4c technique in the MATLAB software 2016.The influence of pertinent flow parameters is discussed and displayed through different sketches.Major and important results are summarized in the conclusion section.Furthermore,in both cases of wall-through flow(i.e.,suction and injection effects),the porosity parameters increase the flow speed,and decrease the heat transport and the influence of drag forces.

Critical current degradation in an epoxy-impregnated rare-earth Ba_(2)Cu_(3)O_(7-x)coated conductor caused by damage during a quench

High-temperature superconducting(HTS)rare-earth Ba_(2)Cu_(3)O_(7-x)(REBCO)coated conductors(CCs)have significant potential in high-current and high-field applications.However,owing to the weak interface strength of the laminated composite REBCO CCs,the damage induced by the thermal mismatch stress under a combination of epoxy impregnation,cooling,and quenching can cause premature degradation of the critical current.In this study,a three-dimensional(3D)electromagnetic-thermal-mechanical model based on the H-formulation and cohesive zone model(CZM)is developed to study the critical current degradation characteristics in an epoxy-impregnated REBCO CC caused by the damage during a quench.The temperature variation,critical current degradation of the REBCO CC,and its degradation onset temperature calculated by the numerical model are in agreement with the experimental data taken from the literature.The delamination of the REBCO CC predicted by the numerical model is consistent with the experimental result.The numerical results also indicate that the shear stress is the main contributor to the damage propagation inside the REBCO CC.The premature degradation of the critical current during a quench is closely related to the interface shear strength inside the REBCO CC.Finally,the effects of the coefficient of thermal expansion(CTE)of the epoxy resin,thickness of the substrate,and substrate material on the critical current degradation characteristics of the epoxy-impregnated REBCO CC during a quench are also discussed.These results help us understand the relationship between the current-carrying degradation and damage in the HTS applications.

The Boundary Element Method for Ordinary State-Based Peridynamics

The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.

Nonlinear Filtering With Sample-Based Approximation Under Constrained Communication:Progress, Insights and Trends

The nonlinear filtering problem has enduringly been an active research topic in both academia and industry due to its ever-growing theoretical importance and practical significance.The main objective of nonlinear filtering is to infer the states of a nonlinear dynamical system of interest based on the available noisy measurements. In recent years, the advance of network communication technology has not only popularized the networked systems with apparent advantages in terms of installation,cost and maintenance, but also brought about a series of challenges to the design of nonlinear filtering algorithms, among which the communication constraint has been recognized as a dominating concern. In this context, a great number of investigations have been launched towards the networked nonlinear filtering problem with communication constraints, and many samplebased nonlinear filters have been developed to deal with the highly nonlinear and/or non-Gaussian scenarios. The aim of this paper is to provide a timely survey about the recent advances on the sample-based networked nonlinear filtering problem from the perspective of communication constraints. More specifically, we first review three important families of sample-based filtering methods known as the unscented Kalman filter, particle filter,and maximum correntropy filter. Then, the latest developments are surveyed with stress on the topics regarding incomplete/imperfect information, limited resources and cyber security.Finally, several challenges and open problems are highlighted to shed some lights on the possible trends of future research in this realm.

INFLUENCE OF NOSE PERTURBATIONS ON BEHAVIORS OF ASYMMETRIC VORTICES OVER SLENDER BODY

The influence of nose perturbations on the behaviors of asymmetric vortices over a slender body with a three-caliber ogive nose is studied in this paper. The tests of a nose-disturbed slender body with surface pressure measurement were conducted at a low speed wind tunnel with subcritical Reynolds number of 1×105 at angle of attack α=50°. The experiment results show that the behaviors and structure of asymmetric vortices over the slender body are mainly controlled by manual perturbation on the nose of body as compared with geometrical minute irregularities on the test model from the machining tolerances. The effect of the perturbation axial location on asymmetric vortices is the strongest if its location is near the model apex. There are four sensitive circumferential locations of manual perturbation at which bistable vortices over the slender body are switched by the perturbation. The flowfield near the reattachment line of lee side is more sensitive to the perturbation, because the saddle point to saddle point topological structure in this reattachment flowfield is unstable. Various types of perturbation do not change the perturbation effect on the behaviors of bistable asymmetric vortices.

Space-time evolution rules of acoustic emission location of unloaded coal sample at different loading rates

By using MTS815 rock mechanics test system,a series of acoustic emission(AE) location experiments were performed under unloading confining pressure,increasing the axial stress.The AE space-time evolution regularities and energy releasing characteristics during deformation and failure process of coal of different loading rates are compared,the influence mechanism of loading rates on the microscopic crack evolution were studied,combining the AE characteristics and the macroscopic failure modes of the specimens,and the precursory characteristics of coal failure were also analyzed quantitatively.The results indicate that as the loading rate is higher,the AE activity and the main fracture will begin earlier.The destruction of coal body is mainly the function of shear strain at lower loading rate and tension strain at higher rate,and will transform from brittleness to ductility at critical velocities.When the deformation of the coal is mainly plasticity,the amplitude of the AE ringing counting rate increases largely and the AE energy curves appear an obvious ''step'',which can be defined as the first failure precursor point.Statics of AE information shows that the strongest AE activity begins when the axial stress level was 92-98%,which can be defined as the other failure precursor point.As the loading rate is smaller,the coal more easily reaches the latter precursor point after the first one,so attention should be aroused to prevent dynamic disaster in coal mining when the AE activity reaches the first precursor point.

UNSTEADY FLOWS OF A GENERALIZED SECOND GRADE FLUID WITH THE FRACTIONAL DERIVATIVE MODEL BETWEEN TWO PARALLEL PLATES

The fractional calculus approach in the constitutive relationship model of a generalized second grade fluid is introduced.Exact analytical solutions are obtained for a class of unsteady flows for the generalized second grade fluid with the fractional derivative model between two parallel plates by using the Laplace transform and Fourier transform for fractional calculus.The unsteady flows are generated by the impulsive motion or periodic oscillation of one of the plates.In addition,the solutions of the shear stresses at the plates are also determined.

Effects of constitutive parameters on adiabatic shear localization for ductile metal based on JOHNSON-COOK and gradient plasticity models

By using the widely used JOHNSON-COOK model and the gradient-dependent plasticity to consider microstructural effect beyond the occurrence of shear strain localization,the distributions of local plastic shear strain and deformation in adiabatic shear band(ASB)were analyzed.The peak local plastic shear strain is proportional to the average plastic shear strain,while it is inversely proportional to the critical plastic shear strain corresponding to the peak flow shear stress.The relative plastic shear deformation between the top and base of ASB depends on the thickness of ASB and the average plastic shear strain.A parametric study was carried out to study the influence of constitutive parameters on shear strain localization.Higher values of static shear strength and work to heat conversion factor lead to lower critical plastic shear strain so that the shear localization is more apparent at the same average plastic shear strain.Higher values of strain-hardening exponent,strain rate sensitive coefficient,melting point,thermal capacity and mass density result in higher critical plastic shear strain,leading to less apparent shear localization at the same average plastic shear strain.The strain rate sensitive coefficient has a minor influence on the critical plastic shear strain,the distributions of local plastic shear strain and deformation in ASB.The effect of strain-hardening modulus on the critical plastic shear strain is not monotonous.When the maximum critical plastic shear strain is reached,the least apparent shear localization occurs.

Analysis of localized shear deformation of ductile metal based on gradient-dependent plasticity

Shear localization in linear strain softening heterogeneous material under simple shear was investigated analytically. The closed-form solutions obtained based on gradient plasticity theory considering interactions and interplaying among microstructures due to heterogeneity of metal material show that in the normal direction of shear band, elastic shear displacement is linear; while plastic and total shear displacement are non-linear. Elastic shear strain in the band is uniform and the non-uniformity of total shear displacement stems from localized plastic shear displacement. In the center of the band, plastic and total shear displacement all reach their maximum values. In strain-softening process, elastic displacement decreases as flow shear stress decreases. Contrarily, plastic and total shear displacement increase and manifest shear localization occurs progressively. Under the same shear stress level, plastic and total shear displacement increase as strain softening modulus and elastic shear modulus decrease. The present analytical solutions were compared with many experimental results and the agreement is good.

Level set method for numerical simulation of a cavitation bubble,its growth, collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.