An official website of the United States government
Consider the scattering of a time-harmonic acoustic plane wave by a bounded elastic obstacle which is immersed in a homogeneous acoustic medium. This paper is concerned with an inverse acoustic-elastic interaction problem, which is to determine the location and shape of the elastic obstacle by using either the phased or phaseless far-field data. By introducing the Helmholtz decomposition, the model problem is reduced to a coupled boundary value problem of the Helmholtz equations. The jump relations are studied for the second derivatives of the single-layer potential in order to deduce the corresponding boundary integral equations. The well-posedness is discussed for the solution of the coupled boundary integral equations. An efficient and high order Nyström-type discretization method is proposed for the integral system. A numerical method of nonlinear integral equations is developed for the inverse problem. For the case of phaseless data, we show that the modulus of the far-field pattern is invariant under a translation of the obstacle. To break the translation invariance, an elastic reference ball technique is introduced. We prove that the inverse problem with phaseless far-field pattern has a unique solution under certain conditions. In addition, a numerical method of the reference ball technique based nonlinear integral equations is proposed for the phaseless inverse problem. Numerical experiments are presented to demonstrate the effectiveness and robustness of the proposed methods.
more »
« less
This content will become publicly available on December 1, 2025
Elastic Fields Around Multiple Stiff Prestressed Arcs Located on a Circle
The plane strain problem of an isotropic elastic matrix subjected to uniform far-field load and containing multiple stiff prestressed arcs located on the same circle is considered. The boundary conditions for the arcs are described by those of either Gurtin–Murdoch or Steigmann–Ogden theories in which the arcs are endowed with their own elastic energies. The material parameters for each arc can in general be different. The problem is reduced to the system of real variables hypersingular boundary integral equations in terms of two scalar unknowns expressed via the components of the stress tensors of the arcs. The unknowns are approximated by the series of trigonometric functions that are multiplied by the square root weight functions to allow for automatic incorporation of the tip conditions. The coefficients in series are found from the system of linear algebraic equations that are solved using the collocation method. The expressions for the stress intensity factors are derived and numerical examples are presented to illustrate the influence of governing dimensionless parameters.
more »
« less
- Award ID(s):
- 2112894
- PAR ID:
- 10631461
- Publisher / Repository:
- ASME
- Date Published:
- Journal Name:
- Journal of Applied Mechanics
- Volume:
- 91
- Issue:
- 12
- ISSN:
- 0021-8936
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
A linear elastic circular disc is analyzed under a self-equilibrated system of loads applied along its boundary. A distinctive feature of the investigation, conducted using complex variable analysis, is the assumption that the material is incompressible (in its linearized approximation), rendering the governing equations formally identical to those of Stokes flow in viscous fluids. After deriving a general solution to the problem, an isoperimetric constraint is introduced at the boundary to enforce inextensibility. This effect can be physically realized, for example, by attaching an inextensible elastic rod with negligible bending stiffness to the perimeter. Although the combined imposition of material incompressibility and boundary inextensibility theoretically prevents any deformation of the disc, it is shown that the problem still admits non-trivial solutions. This apparent paradox is resolved by recognizing the approximations inherent in the linearized theory, as confirmed by a geometrically nonlinear numerical analysis. Nonetheless, the linear solution retains significance: it may represent a valid stress distribution within a rigid system and can identify critical conditions of interest for design applications.more » « less
-
Consider the elastic scattering of a time-harmonic wave by multiple well-separated rigid particles with smooth boundaries in two dimensions. Instead of using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of the Helmholtz equation. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations is established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving elastic scattering problems with multiple particles.more » « less
-
Abstract We consider an optimal control problem where the state equations are a coupled hyperbolic–elliptic system. This system arises in elastodynamics with piezoelectric effects—the elastic stress tensor is a function of elastic displacement and electric potential. The electric flux acts as the control variable and bound constraints on the control are considered. We develop a complete analysis for the state equations and the control problem. The requisite regularity on the control, to show the well-posedness of the state equations, is enforced using the cost functional. We rigorously derive the first-order necessary and sufficient conditions using adjoint equations and further study their well-posedness. For spatially discrete (time-continuous) problems, we show the convergence of our numerical scheme. Three-dimensional numerical experiments are provided showing convergence properties of a fully discrete method and the practical applicability of our approach.more » « less
-
In this manuscript, we present a multiscale Adaptive Reduced-Order Modeling (AROM) framework to efficiently simulate the response of heterogeneous composite microstructures under interfacial and volumetric damage. This framework builds on the eigendeformation-based reduced-order homogenization model (EHM), which is based on the transformation field analysis (TFA) and operates in the context of computational homogenization with a focus on model order reduction of the microscale problem. EHM pre-computes certain microstructure information by solving a series of linear elastic problems defined over the fully resolved microstructure (i.e., concentration tensors, interaction tensors) and approximates the microscale problem using a much smaller basis spanned over subdomains (also called parts) of the microstructure. Using this reduced basis, and prescribed spatial variation of inelastic response fields over the parts, the microscale problem leads to a set of algebraic equations with part-wise responses as unknowns, instead of node-wise displacements as in finite element analysis. The volumetric and interfacial influence functions are calculated by using the Interface enriched Generalized Finite Element Method (IGFEM) to compute the coefficient tensors, in which the finite element discretization does not need to conform to the material interfaces. AROM takes advantage of pre-computed coefficient tensors associated with the finest ROM and efficiently computes the coefficient tensors of a series of gradually coarsening ROMs. During the multiscale analysis stage, the simulation starts with a coarse ROM which can capture the initial elastic response well. As the loading continues and response in certain parts of the microstructure starts to localize, the analysis adaptively switches to the next level of refined ROM to better capture those local responses. The performance of AROM is evaluated by comparing the results with regular EHM (no adaptive refinement) and IGFEM under different loading conditions and failure modes for various 2D and 3D microstructures. The proposed AROM provides an efficient way to model history-dependent nonlinear responses for composite materials under localized interface failure and phase damage.more » « less
