Modeling the contributions to acoustic nonlinearity from complex dislocation networks using 3D dislocation dynamics

C C. L. Williams (Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,) N N. Bertin (Lawrence Livermore National Laboratory 2 , Livermore, CA 94550) F F. Costanzo (Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,) M M. H. Lear (Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,) P P. Shokouhi (Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,)

Abstract

Nonlinear ultrasonic parameters are highly sensitive to microstructural features that affect macroscale material behavior, providing a nondestructive means to characterize their evolution. Although dislocations are known to be a strong source of acoustic nonlinearity, establishing quantitative links between the acoustic nonlinearity parameter (β), measured via Second Harmonic Generation, and dislocation morphology—such as dislocation length and density—remains an open challenge. This work advances the numerical modeling of dislocation–β relationships using 3D dislocation dynamics (DD) simulations in two approaches: a “static” method computing strain and stress fields from dislocation configurations in the absence of external loading, and a “quasi-static” method to estimate β from the curvature of dislocation lines under applied load. First, the static method is combined with finite element analysis to investigate a recent assertion that heterogeneous initial strain fields can induce higher harmonic generation in a linear elastic medium; the present results do not corroborate this outcome. Then, the quasi-static method is applied to multiple-dislocation scenarios through parametric studies, revealing behaviors not predicted by analytical models, such as the competing interactions of edge and screw dislocations and the significant influence of applied stress on β. Finally, the simulations are used to model SHG experimental results and validate the hypothesis that β can decrease during plastic deformation, despite increasing dislocation density. As the DD code used here is open-source, it provides a practical platform for future investigation into microstructure–β relationships important to the interpretation of SHG results.

Article Details

Volume / Issue Vol. 139, Issue 12
Published March 28, 2026
ISSN 0021-8979
Publisher American Institute of Physics

Journal Info

Journal of Applied Physics

American Institute of Physics

ISSN: 0021-8979 Physical Sciences

Authors (5)

C

C. L. Williams

Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,

N

N. Bertin

Lawrence Livermore National Laboratory 2 , Livermore, CA 94550

F

F. Costanzo

Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,

M

M. H. Lear

Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,

P

P. Shokouhi

Department of Engineering Science and Mechanics, The Pennsylvania State University 1 , University Park, Pennsylvania 16802,