BANDEAU4_ESMC_copie_1.jpg

Non-reciprocity for the time-modulated wave equation and diffusion equation through the lens of high-order homogenization
Marie Touboul  1, 2, 3, *@  , Bruno Lombard  4@  , Raphaël Assier  5@  , Sébastien Guenneau  3, 6@  , Richard V. Craster  7, 8, 9@  
1 : Department of Mathematics [Imperial College London]
2 : Propagation des Ondes : Étude Mathématique et Simulation
CNRS, ENSTA Paris, L'Institut National de Recherche en Informatique et e n Automatique (INRIA), Institut Polytechnique de Paris
3 : UMI 2004 Abraham de Moivre-CNRS, Imperial College London
4 : Laboratoire de Mécanique et d'Acoustique [Marseille]  (LMA)
Aix Marseille Université, Ecole Centrale de Marseille, Centre National de la Recherche Scientifique, Aix Marseille Université : UMR7031, Ecole Centrale de Marseille : UMR7031, Centre National de la Recherche Scientifique : UMR7031
4 impasse Nikola TeslaCS 4000613453 Marseille Cedex 13 -  France
5 : University of Manchester [Manchester]
6 : Blackett Laboratory
7 : Department of Mathematics [Imperial College London]
Imperial College London, 180 Queen's Gate, London SW7 2AZ, United Kingdom -  United Kingdom
8 : UMI 2004 Abraham de Moivre-CNRS [Imperial College London]
London SW7 2AZ, UK -  United Kingdom
9 : Department of Mechanical Engineering [Imperial College London]
Imperial College London, London SW7 2AZ -  United Kingdom
* : Corresponding author

Time-modulated laminate for the wave equation

We first consider a 1D periodic elastic medium in the antiplane elasticity framework. Moreover, the physical parameters, i.e. the mass density and the shear modulus, are modulated in space and time in a wave-like fashion with modulation speed c.We further assume that the constant modulation speed is such that the configuration is in the subsonic or supersonic cases. The range of modulation speed considered here excludes the sonic case, where an effect similar to a sonic blow-up occurs (Cassedy, 1963).

Homogenized equation at the leading order

Assuming that the characteristic wavelength is much larger than the periodicity, we introduce a small parameter for the low-frequency setting. Using the two-scale asymptotic technique, homogenization can be applied. At the leading order, the effective equation is obtained, where three effective parameters appear: an effective mass density, an effective shear modulus, and an effective Willis coupling coefficient W0 due to the modulation. The associated dispersion relation describes the behavior of wave propagation in the modulated medium.

This effective equation implies that reciprocity is broken as soon as W0 is non-zero, as there is no longer symmetry between the wavevector k and -k. However, if only one physical parameter (either the shear modulus or the density) is modulated, the Willis coupling coefficient W0 vanishes and reciprocity is recovered. This is consistent with the literature (Nassar, 2017; Huidobro, 2021), but inconsistent with Bloch-Floquet analysis (Cassedy, 1963), which highlights a limitation of the leading-order homogenized model.

Homogenized equation at the second order

To correct the limitations of the leading-order model, we push the homogenization to the second order (Touboul, 2024a). When the density is constant, the effective dispersion relation for the second-order model developed in this work shows that non-reciprocity is ensured due to a non-zero term breaking the symmetry between k and -k. A similar result holds when the other parameter is constant.

The accuracy of this model is illustrated for a bilayered medium. The dispersion relation computed through a Bloch-Floquet analysis in a moving frame is compared to the one derived from the homogenized models. Time-domain simulations in the microstructured medium are also performed to illustrate non-reciprocity.

The case of the diffusion equation

This same methodology is applied to the diffusion equation. In this case, a time-modulation in a wave-like fashion leads to non-reciprocal behavior and propagation of the field, while it would be purely attenuated in a non-modulated medium. However, leading-order models (Torrent, 2018) or first-order models (Xu, 2022) fail to describe this phenomenon when only one parameter is modulated. As with the wave equation, this issue is corrected by considering second-order homogenization (Touboul, 2024b). The same methodology is also applied to the case when the density is modulated in the heat equation, leading, therefore, to a corrective advective term (Li, 2022) that cancels out non-reciprocity at the leading order but not at the second order.

References

  • Nassar H., Xu X., Norris A., and Huang G., "Modulated phononic crystals: Non-reciprocal wave propagation and Willis materials," J. Mech. Phys. Solids, vol. 101, pp. 10–29, 2017.
  • Huidobro P., Silveirinha M., Galiffi E., and Pendry J., "Homogenization theory of space-time metamaterials," Phys. Rev. Appl., vol. 16, pp. 014044, 2021.
  • Xu L., Xu G., Li J., Li Y., Huang J., and Qiu C.W., "Thermal Willis Coupling in Spatiotemporal Diffusive Metamaterials," Phys. Rev. Lett., vol. 129, 2022.
  • Torrent D., Poncelet O., and Batsale J.C., "Nonreciprocal Thermal Material by Spatiotemporal Modulation," Phys. Rev. Lett., vol. 120, pp. 125501, 2018.
  • Cassedy E., and Oliner A., "Dispersion relations in time-space periodic media: Part I—Stable interactions," Proceedings of the IEEE, vol. 51, pp. 1342–1359, 1963.
  • Touboul M., Lombard B., Assier R., Guenneau S., and Craster R.V., "High-order homogenisation of the time-modulated wave equation: non-reciprocity for a single varying parameter," Proc. R. Soc. A, vol. 480, No. 2289, 2024a.
  • Touboul M., Lombard B., Assier R., Guenneau S., and Craster R.V., " Propagation and non-reciprocity in time-modulated diffusion through the lens of high-order homogenization," Proc. R. Soc. A, vol. 480, No. 2301, 2024b.
  • Li J., Li Y., Cao P.-C., Qi M., Zheng X., Peng Y.-G., Li B., Zhu X.-F., Alù A., Chen H., Qiu C.-W., " Reciprocity of thermal diffusion in time-modulated systems," Nature Communications, Vol. 13, No. 1, 2022.

 


Loading... Loading...