skip to main content


Title: Back-calculation of soil parameters from displacement-controlled cavity expansion under geostatic stress by FEM and machine learning
Abstract

Estimating soil properties from the mechanical reaction to a displacement is a common strategy, used not only in in situ soil characterization (e.g., pressuremeter and dilatometer tests) but also by biological organisms (e.g., roots, earthworms, razor clams), which sense stresses to explore the subsurface. Still, the absence of analytical solutions to predict the stress and deformation fields around cavities subject to geostatic stress, has prevented the development of characterization methods that resemble the strategies adopted by nature. We use the finite element method (FEM) to model the displacement-controlled expansion of cavities under a wide range of stress conditions and soil properties. The radial stress distribution at the cavity wall during expansion is extracted. Then, methods are proposed to prepare, transform and use such stress distributions to back-calculate the far field stresses and the mechanical parameters of the material around the cavity (Mohr-Coulomb friction angle$$\phi $$ϕ, Young’s modulusE). Results show that: (i) The initial stress distribution around the cavity can be fitted to a sum of cosines to estimate the far field stresses; (ii) By encoding the stress distribution as intensity images, in addition to certain scalar parameters, convolutional neural networks can consistently and accurately back-calculate the friction angle and Young’s modulus of the soil.

 
more » « less
NSF-PAR ID:
10372030
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Acta Geotechnica
Volume:
18
Issue:
4
ISSN:
1861-1125
Page Range / eLocation ID:
p. 1755-1768
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The crystal structure and bonding environment of K2Ca(CO3)2bütschliite were probed under isothermal compression via Raman spectroscopy to 95 GPa and single crystal and powder X-ray diffraction to 12 and 68 GPa, respectively. A second order Birch-Murnaghan equation of state fit to the X-ray data yields a bulk modulus,$${K}_{0}=46.9$$K0=46.9GPa with an imposed value of$${K}_{0}^{\prime}= 4$$K0=4for the ambient pressure phase. Compression of bütschliite is highly anisotropic, with contraction along thec-axis accounting for most of the volume change. Bütschliite undergoes a phase transition to a monoclinicC2/mstructure at around 6 GPa, mirroring polymorphism within isostructural borates. A fit to the compression data of the monoclinic phase yields$${V}_{0}=322.2$$V0=322.2 Å3$$,$$,$${K}_{0}=24.8$$K0=24.8GPa and$${K}_{0}^{\prime}=4.0$$K0=4.0using a third order fit; the ability to access different compression mechanisms gives rise to a more compressible material than the low-pressure phase. In particular, compression of theC2/mphase involves interlayer displacement and twisting of the [CO3] units, and an increase in coordination number of the K+ion. Three more phase transitions, at ~ 28, 34, and 37 GPa occur based on the Raman spectra and powder diffraction data: these give rise to new [CO3] bonding environments within the structure.

     
    more » « less
  2. Abstract

    As the use of spectral/hpelement methods, and high-order finite element methods in general, continues to spread, community efforts to create efficient, optimized algorithms associated with fundamental high-order operations have grown. Core tasks such as solution expansion evaluation at quadrature points, stiffness and mass matrix generation, and matrix assembly have received tremendous attention. With the expansion of the types of problems to which high-order methods are applied, and correspondingly the growth in types of numerical tasks accomplished through high-order methods, the number and types of these core operations broaden. This work focuses on solution expansion evaluation at arbitrary points within an element. This operation is core to many postprocessing applications such as evaluation of streamlines and pathlines, as well as to field projection techniques such as mortaring. We expand barycentric interpolation techniques developed on an interval to 2D (triangles and quadrilaterals) and 3D (tetrahedra, prisms, pyramids, and hexahedra) spectral/hpelement methods. We provide efficient algorithms for their implementations, and demonstrate their effectiveness using the spectral/hpelement libraryNektar++by running a series of baseline evaluations against the ‘standard’ Lagrangian method, where an interpolation matrix is generated and matrix-multiplication applied to evaluate a point at a given location. We present results from a rigorous series of benchmarking tests for a variety of element shapes, polynomial orders and dimensions. We show that when the point of interest is to be repeatedly evaluated, the barycentric method performs at worst$$50\%$$50%slower, when compared to a cached matrix evaluation. However, when the point of interest changes repeatedly so that the interpolation matrix must be regenerated in the ‘standard’ approach, the barycentric method yields far greater performance, with a minimum speedup factor of$$7\times $$7×. Furthermore, when derivatives of the solution evaluation are also required, the barycentric method in general slightly outperforms the cached interpolation matrix method across all elements and orders, with an up to$$30\%$$30%speedup. Finally we investigate a real-world example of scalar transport using a non-conformal discontinuous Galerkin simulation, in which we observe around$$6\times $$6×speedup in computational time for the barycentric method compared to the matrix-based approach. We also explore the complexity of both interpolation methods and show that the barycentric interpolation method requires$${\mathcal {O}}(k)$$O(k)storage compared to a best case space complexity of$${\mathcal {O}}(k^2)$$O(k2)for the Lagrangian interpolation matrix method.

     
    more » « less
  3. Abstract

    We propose a new observable for the measurement of the forward–backward asymmetry$$(A_{FB})$$(AFB)in Drell–Yan lepton production. At hadron colliders, the$$A_{FB}$$AFBdistribution is sensitive to both the electroweak (EW) fundamental parameter$$\sin ^{2} \theta _{W}$$sin2θW, the weak mixing angle, and the parton distribution functions (PDFs). Hence, the determination of$$\sin ^{2} \theta _{W}$$sin2θWand the updating of PDFs by directly using the same$$A_{FB}$$AFBspectrum are strongly correlated. This correlation would introduce large bias or uncertainty into both precise measurements of EW and PDF sectors. In this article, we show that the sensitivity of$$A_{FB}$$AFBon$$\sin ^{2} \theta _{W}$$sin2θWis dominated by its average value around theZpole region, while the shape (or gradient) of the$$A_{FB}$$AFBspectrum is insensitive to$$\sin ^{2} \theta _{W}$$sin2θWand contains important information on the PDF modeling. Accordingly, a new observable related to the gradient of the spectrum is introduced, and demonstrated to be able to significantly reduce the potential bias on the determination of$$\sin ^{2} \theta _{W}$$sin2θWwhen updating the PDFs using the same$$A_{FB}$$AFBdata.

     
    more » « less
  4. Abstract

    We investigate electroabsorption (EA) in organic semiconductor microcavities to understand whether strong light-matter coupling non-trivially alters their nonlinear optical [$${\chi }^{(3)}\left(\omega,{{{{\mathrm{0,0}}}}}\right)$$χ(3)ω,0, 0] response. Focusing on strongly-absorbing squaraine (SQ) molecules dispersed in a wide-gap host matrix, we find that classical transfer matrix modeling accurately captures the EA response of low concentration SQ microcavities with a vacuum Rabi splitting of$$\hslash \Omega \approx 200$$Ω200meV, but fails for high concentration cavities with$$\hslash \Omega \approx 420$$Ω420meV. Rather than new physics in the ultrastrong coupling regime, however, we attribute the discrepancy at high SQ concentration to a nearly dark H-aggregate state below the SQ exciton transition, which goes undetected in the optical constant dispersion on which the transfer matrix model is based, but nonetheless interacts with and enhances the EA response of the lower polariton mode. These results indicate that strong coupling can be used to manipulate EA (and presumably other optical nonlinearities) from organic microcavities by controlling the energy of polariton modes relative to other states in the system, but it does not alter the intrinsic optical nonlinearity of the organic semiconductor inside the cavity.

     
    more » « less
  5. Abstract

    Cavities concentrate light and enhance its interaction with matter. Confining to microscopic volumes is necessary for many applications but space constraints in such cavities limit the design freedom. Here we demonstrate stable optical microcavities by counteracting the phase evolution of the cavity modes using an amorphous Silicon metasurface as cavity end mirror. Careful design allows us to limit the metasurface scattering losses at telecom wavelengths to less than 2% and using a distributed Bragg reflector as metasurface substrate ensures high reflectivity. Our demonstration experimentally achieves telecom-wavelength microcavities with quality factors of up to 4600, spectral resonance linewidths below 0.4 nm, and mode volumes below$$2.7{\lambda }^{3}$$2.7λ3. The method introduces freedom to stabilize modes with arbitrary transverse intensity profiles and to design cavity-enhanced hologram modes. Our approach introduces the nanoscopic light control capabilities of dielectric metasurfaces to cavity electrodynamics and is industrially scalable using semiconductor manufacturing processes.

     
    more » « less