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.

ANALYTICAL SOLUTIONS FOR ELASTOSTATIC PROBLEMS OF PARTICLE-AND FIBER-REINFORCED COMPOSITES WITH INHOMOGENEOUS INTERPHASES

By transforming the governing equations for displacement components into Riccati equations, analytical solutions for displacements, strains and stresses for Representive Volume Elements (RVEs) of particle_ and fiber_reinforced composites containing inhomo geneous interphases were obtained. The analytical solutions derived here are new and general for power_law variations of the elastic moduli of the inhomogeneous interphases. Given a power exponent, analytical expressions for the bulk moduli of the composites with inho mogeneous interphases can be obtained. By changing the power exponent and the coefficients of the power terms, the solutions derived here can be applied to inhomogeneous interphases with many different property profiles. The results show that the modulus variation and the thickness of the inhomogeneous interphase have great effect on the bulk moduli of the composites. The particle will exhibit a sort of “size effect”, if there is an interphase.

Design and optimization of a SiC thermal emitter/absorber composed of periodic microstructures based on a non-linear method

Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.

Atomic diffusion across Ni_(50)Ti_(50)Cu explosive welding interface:Diffusion layer thickness and atomic concentration distribution

Molecular dynamics simulations are carried out to study atomic diffusion in the explosive welding process of NisoTis0-Cu （at.%）. By using a hybrid method which combines molecular dynamics simulation and classical diffusion the- ory, the thickness of the diffusion layer and the atomic concentration distribution across the welding interface are obtained. The results indicate that the concentration distribution curves at different times have a geometric similarity. According to the geometric similarity, the atomic concentration distribution at any time in explosive welding can be calculated. NisoTis0- Cu explosive welding and scanning electron microscope experiments are done to verify the results. The simulation results and the experimental results are in good agreement.

INTERFACE EFFECT ON THE EFFECTIVE BULK MODULUS OF A PARTICLE-REINFORCED COMPOSITE

Classical micromechanical methods for calculating the effective moduli of a heteroge- neous material are generalized to include the interface(surface)effect.By using Hashin's Composite Sphere Assemblage(CSA)model,a new expression of the bulk modulus for a particle-reinforced com- posite is derived.It is emphasized that the present study is within the finite-deformation framework such that the effective properties are not influenced by the interface stress itself solely,but influenced by the change of the interface stress due to changes of the shape and size of the interface.Hence some inadequacies in previous papers are pointed out.

A three-dimensional Eulerian method for the numerical simulation of high-velocity impact problems

In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.

Numerical study of dynamic phase transitions in shock tube

Shock tube problem of a van der Waals fluid with a relaxation model was investigated. In the limit of relaxation parameter tending towards zero, this model yields a specific Riemann solver. Relaxing and relaxed schemes were derived. For an incident shock in a fixed tube, numerical simulations show convergence toward the Riemann solution in one space dimension. Impact of parameters was studied theoretically and numerically. For certain initial shock profiles, nonclassical reflecting wave was observed. In two space dimensions, the effect of curved wave fronts was studied, and some interesting wave patterns were exposed.

A New Method for Separation of Waves in Improving the Conventional SHPB Technique

We present a new wave separation method to enable the split Hopkinson pressure bar （SHPB） technique to break through the limitation of the length of the incident bar and greatly to increase its measurable maximum strain. At the same time the dispersion effect of the elastic wave is significantly reduced. The fundamental principle of the new method is proven rigorously. The feasibility and credibility of the new method are also verified by experiments.

LIQUID-GAS COEXISTENCE EQUILIBRIUM IN A RELAXATION MODEL

Stability of liquid-gas coexistence equilibrium in a relaxation model for isothermal phase transition in a sealed one-dimensional tube was discussed. With matched asymptotic expansion, a linear system for first order perturbations was derived formally. By solving this system analytically, it is shown that small initial perturbations are damped out in general; yet they may maintain at certain level for special cases. Numerical evidence is presented. This manifests the regularization effects of relaxation.

Matching Boundary Conditions for a Heterogeneous Atomic Chain

We propose accurate boundary treatments for a heterogeneous atomic chain, in terms of matching boundary conditions （MBCs）. The main challenge lies in reproducing the physical reflection across the boundary to a correct amount. With reflection coefficients we demonstrate that the accuracy is improved when more atoms are used under the boundary condition. The inclusion of an atom in the embedded sublattice B may considerably enhance the performance. Numerical testing illustrates the effectiveness of the proposed MBCs.

Effect of geometric parameters and grazing incidence on magnetic polaritons excited in 1D multi-groove gratings

In this paper,effect of geometric parameters and grazing incidence on excitation of magnetic polaritons(MPs)in a 1D multigroove grating with different groove depths made of silver was studied.Numerical results reveal that when the distance between grooves is sufficiently small,the resonance wavelengths of MPs excited in the grooves of unequal depths exhibit red shift with decrease of the distance,contrary to the case with equal groove depths.The shift of the MP resonance wavelengths was explained with the LC circuit model.Furthermore,it was found that the resonance wavelengths of MPs depend linearly on the groove depths except when the difference between the groove depths is small.When the grooves have equal depths,a large drop of the absorptance can occur due to the interaction and cancellation of the electric field vectors in the region between the grooves.Finally,the results show that when a TM wave is at grazing incidence,MPs can be excited simultaneously in the grooves at a blue-shifted wavelength due to reduction of the effective capacitance,resulting in a dramatic enhancement of the absorptance.Therefore,the results in this work may provide useful guidance in the design of wavelength-selective absorbers based on MPs.

An energy-consistent fracture model for ferroelectrics

The fracture behavior of ferroelectrics has been intensively studied in recent decades, though currently a widely accepted fracture mechanism is still lacking. In this work, enlightened by previous experimental observations that crack propagation in ferroelectrics is always accompanied by domain switching, we propose a micromechanical model in which both crack propagation and domain switching are controlled by energy-based criteria. Both electric energy and mechanical energy can induce domain switching, while only mechanical energy can drive crack propagation. Furthermore, constrained domain switching is considered in this model, leading to the gradient domain switching zone near the crack tip. Analysis results show that stress-induced ferroelastic switching always has a toughening effect as the mechanical energy release rate serves as the driving force for both fracture and domain switching. In comparison, the electric-field-induced switching may have either a toughening or detoughening effect. The proposed model can qualitatively agree with the existing experimental results.

Development of an Omni-Directional Shear Horizontal Wave Transducer Based on a Radially Poled Piezoelec trie Ring

For the inspection of large plate-like structures,the omni-directional guided wave transducer-based system has been regarded as an effective tool since only a few transducers are required to cover the entire inspection area without blind zones.In comparison with Lamb waves,the shear horizontal(SH)wave is more promising because its fundamental mode is non-dispersive.In this work,we proposed an omni-directional SH wave piezoelec trie t ransducer(OSH-PT)based on a radially poled d24-mode PZT ring.Eoth the finite element simulations and experiments were carried out to demonstrate its performance in generating and receiving SH°wave.Results showed that the radially poled OSH-PT could generate single-mode SHo wave and receive SHo wave only over a wide frequency range from 70 to 200 kHz.The obtained signal-to-noise ratio can reach up to 26 dB in generation and 24 dB in reception.The omni-directivity of this OSH-PT is also very good with the deviation only about 6%in both generation and reception of SHo wave.Considering its easy fabrication,low cost and superior performances,this proposed OSH-PT may promote the applications of SHo wave-based inspection in structural health monitoring and nondestructive testing.

Mesh Deformation Method Based on Mean Value Coordinates Interpolation

Mesh deformation technique is widely used in many application fields, and has re- ceived a lot of attentions in recent years. This paper focuses on the methodology and algorithm of algebraic type mesh deformation for unstructured mesh in numerical discretization. To preserve mesh quality effectively, an algebraic approach for two and three dimensional unstructured mesh is developed based on mean value coordinates interpolation combined with node visibility analysis. The proposed approach firstly performs node visibility analysis to find out the visible boundary for each grid point to be moved, then evaluates the mean value coordinates of each grid point with respect to all vertices on its visible boundary. Thus the displacements of grid points can be calculated by interpolating the boundary movement by the mean value coordinates. Compared with other methods, the proposed method has good deformation capability and predictable com- putational cost, with no need to select parameters or functions. Applications of mesh deformation in different fields are presented to demonstrate the effectiveness of the proposed approach. The results of numerical experiments exhibit not only superior deformation capability of the method in traditional applications of fluid dynamic grid, but also great potential in modeling for large deformation analysis and inverse design problems.

Projection Improved SPAI Preconditioner for FGMRES

Krylov subspace methods are widely used for solving sparse linear algebraic equations,but they rely heavily on preconditioners,and it is difficult to find an effective preconditioner that is efficient and stable for all problems.In this paper,a novel projection strategy including the orthogonal and the oblique projection is proposed to improve the preconditioner,which can enhance the efficiency and stability of iteration.The proposed strategy can be considered as a minimization process,where the orthogonal projection minimizes the energy norm of error and the oblique projection minimizes the 2-norm of the residual,meanwhile they can be regarded as approaches to correct the approximation by solving low-rank inverse of the matrices.The strategy is a wide-ranging approach and provides a way to transform the constant preconditioner into a variable one.The paper discusses in detail the projection strategy for sparse approximate inverse(SPAI)preconditioner applied to flexible GMRES and conducts the numerical test for problems from different applications.The results show that the proposed projection strategy can significantly improve the solving process,especially for some non-converging and slowly convergence systems.