à l’étape a1), les modes de résonance sont déterminés par simulation numérique, la modélisation de la pièce de référence prenant en compte les dimensions et la géométrie de la pièce de référence, la masse volumique apparente du matériau constituant la pièce de référence, le module d’élasticité dudit matériau, et le coefficient de Poisson dudit matériau ; à l’étape a2), on compare des modèles numériques tridimensionnels de la pièce de référence représentant chacun un champ desdites contraintes ou desdites déformations mécaniques généré lorsque la pièce de référence est en résonance à une fréquence de résonance respective ; le volume de la région d'intérêt est inférieur à 0,2 fois et supérieur à 0,01 fois le volume de la pièce de référence ; les étapes a) et b) sont réalisées simultanément, la recherche du mode de sollicitation optimal étant réalisée de la manière suivante : - simulation numérique d’une pluralité de modes de sollicitation, la simulation numérique déterminant pour chaque mode de sollicitation, un champ de contraintes mécaniques dans la pièce de référence et une onde de sortie théorique (c'est-à-dire une simulation de l’onde de sortie dans ce mode de sollicitation) ; - analyse des champs de contraintes mécaniques et des ondes de sortie théoriques de manière à sélectionner, comme mode de sollicitation optimal, le mode de sollicitation générant, dans la région d’intérêt, une contrainte mécanique ou une déformation maximale et mettant en résonance la pièce de référence ; à l’étape c), le paramètre de non-linéarité est la pente d’une droite représentative de l’évolution d’un décalage de fréquence en fonction de l’évolution de l’amplitude de l’onde de sortie lorsque ladite amplitude de l’onde de sortie est modifiée, de préférence augmentée, depuis le mode de sollicitation optimal, le décalage de fréquence, pour une amplitude de l’onde de sortie, étant le rapport de la valeur absolue de la différence entre la fréquence de résonance optimale (f 0 ) dans le mode de résonance optimal et la fréquence de résonance (f) déterminée pour ladite amplitude de l’onde de sortie, divisé par la fréquence de résonance optimale (f 0 ) ; à l’étape d), on compare le paramètre de non-linéarité de la pièce à tester à un seuil déterminé à partir du paramètre de non-linéarité de la pièce de référence, puis on classe la pièce à tester en fonction de la différence entre le paramètre de non-linéarité de la pièce à tester et le seuil, de préférence en fonction du signe de ladite différence ; le défaut est un espace vide au sein de la pièce à tester, ou un espace rempli par un matériau différent du reste de la pièce à tester ; la pièce est en un matériau inorganique ; la pièce est en un métal, de préférence fritté, un matériau céramique, un matériau vitrocéramique ou en un mélange de ces matériaux, en particulier un composite ; dans le mode de sollicitation optimal, le pic principal d’un spectre de fréquence de l’onde de sortie est à une fréquence comprise entre 1 Hz et 200 KHz, de préférence inférieure à 100 KHz, de préférence supérieure à 20 KHz ; l’onde de sortie est à une onde acoustique ; ladite pièce à tester est choisie parmi : - un bloc ou linteau de gorge, - un bloc de cuve, - une brique ou un bloc de mur ou de côté, - un bloc de coin ou d’angle, - une brique à nez, - un bloc ou une dalle de sole, - une brique ou un sommier de voute, - une brique ou un bloc d’entourage de tuyère, - une brique de trou ou chenal de coulée, - un bloc porte électrode, - une pièce de bec réfractaire de four verrier, - un bloc d’injecteur, - une gorge de four de verrerie, - une pièce d’un échangeur thermique du four, - une plaque ou une tuile réfractaire d’un revêtement de chaudière, - une coquille de protection de tube de réchauffeur d’un incinérateur, - une tuile d’un incinérateur, - une pièce céramique d’un absorbeur solaire, - une tuile ou une pièce de protection d’une chambre de turbine de combustion.