In Chaps. 14 and 15 , we applied Eshelby’s interior solutions and the equivalent inclusion approach. In this chapter, we investigate the development of flanking structures as an example to use Eshelby’s exterior solutions. Flanking structures [Fig. 16.1; Passchier (J Struct Geol 23(6–7):951–962, 2001)] are deflections of linear and planar fabric elements around a cross-cutting element such as a fault, vein, dyke, or any other heterogeneity. Many authors have investigated the development of flanking structures by a combination of field, numerical modeling, analogue modeling, and analytical approaches [Exner and Dabrowski (J Struct Geol 32(12):2009–2021, 2010); Exner et al. (J Struct Geol 26(12):2191–2201, 2004); Grasemann et al. (J Struct Geol 25(1):19–34, 2003), (J Struct Geol 33(11):1650–1661, 2011); Grasemann and Stuwe (J Struct Geol 23(4):715–724, 2001); Kocher and Mancktelow (J Struct Geol 27(8):1346–1354, 2005), (J Struct Geol 28(7):1139–1145, 2006); Wiesmayr and Brasemann (J Struct Geol 27(2):249–264, 2005)]. These studies have greatly improved the understanding of flanking structures. The main conclusions are that the final geometrical features of flanking structures depend on the initial conditions of the cutting element (shape, orientation, relative rheology) and the characteristics of the flow field. One must therefore be cautious when interpreting the kinematics of flanking structures (e.g., Grasemann et al. (J Struct Geol 25(1):19–34, 2003); Wiesmayr and Brasemann (J Struct Geol 27(2):249–264, 2005)]. Because flanking structures result from the progressive deformation around a cutting element heterogeneity, we can regard the cutting element as an Eshelby ellipsoid and apply the exterior solutions (Chap. 11 ) to the development of flanking structures. Compared to previous investigations, the application of Eshelby’s exterior solutions has the following advantages. First, there is no limit to the flow type and the initial geometrical conditions of the cutting element (any ellipsoidal shape, arbitrary orientation with respect to the host element and the bulk flow field). Second, there is no limit to the amount of strain that can be reached because the approach does not rely on meshing. Third, the approach can handle anisotropic viscosity readily, although at a more expensive computational cost (Jiang and Bhandari 2018; Bhandari 2021). Furthermore, although the approach is based on linear viscous materials so far (Chaps. 10 and 11 ), it can be extended to power-law viscous materials by a linearization approach [Jiang (J Struct Geol 68:247–272, 2014), (Tectonophysics 693:116–142, 2016); Chap. 17 ]. This chapter first summarizes the principal equations related to the exterior mechanical fields in the vicinity of a heterogeneous ellipsoid. An outline of the algorithm for the numerical modeling of flanking structures is then presented. The algorithm is implemented in MATLAB. The MATLAB programs are provided in the online resource of the book. The program is applicable to anisotropic viscous materials under any flow field. As examples of using the numerical approach and the MATLAB program, we present numerical results for the development of flanking structures in plane-strain general shearing flows. The cutting element and the host medium are assumed to be isotropic and linear viscous materials. Furthermore, for simplicity, we limit ourselves to cutting elements with one principal axis parallel to the flow vorticity so that the whole deforming system has monoclinic symmetry.
Read more