Problem Of Propagation Research Articles

Prismatic-element SBPML coupled with SBFEM for 3D infinite transient wave problems

In this paper, an enhanced prismatic-element scaled boundary perfectly matched layer (SBPML) is developed, which is a novel time-domain artificial boundary method for 3D infinite wave problems. The SBPML permits the utilization of an artificial boundary with general geometry and can consider planar physical surfaces and interfaces extending to infinity. Moreover, this enhancement enables the seamless integration with surface meshes of arbitrary polygonal faces (such as triangles, quadrangles, pentagons, heptagons, octagons, etc.) at the truncation boundaries. The interior domain can be modelled by polyhedrons formulated by the scaled boundary finite element method (SBFEM). The geometric shape of the SBPML elements in the proposed method, within the local coordinate system, takes the form of a regular prism with regular polygonal top and bottom faces. Consequently, it is referred to as the Prismatic-element SBPML. Local scaled boundary coordinates based on the polygonal shape function are firstly introduced into the truncated infinite domain to describe its general geometry properties and provide the coupled capabilities with arbitrary polygonal face meshes at the truncation boundaries of numerical models. Then, a complex stretching function from PML is applied to radial direction of the local scaled boundary coordinates to map the physical space onto the complex space, resulting in a SBPML domain. Upon spatial discretization, the proposed method is formulated using a 2nd-order mixed displacement-stress unsplit-field formulation, allowing seamless coupling with FEM and SBFEM originally designed for the interior computational domain. Finally, three benchmark examples of wave propagation problems in unbounded domains and two soil-structure interaction problems of engineering problems are presented to demonstrate the accuracy and robustness of the Prismatic-element SBPML.

Duality analysis of interior penalty discontinuous Galerkin methods under minimal regularity and application to the a priori and a posteriori error analysis of Helmholtz problems

We consider interior penalty discontinuous Galerkin discretizations of time-harmonic wave propagation problems modeled by the Helmholtz equation, and derive novel a priori and a posteriori estimates. Our analysis classically relies on duality arguments of Aubin–Nitsche type, and its originality is that it applies under minimal regularity assumptions. The estimates we obtain directly generalize known results for conforming discretizations, namely that the discrete solution is optimal in a suitable energy norm and that the error can be explicitly controlled by a posteriori estimators, provided the mesh is sufficiently fine.

Joint signal-to-signal beat interference mitigation for the field recovery of symmetric carrier-assisted differential detection with low carrier-to-signal power ratio

As a combination of direct detection and coherent detection technologies, self-coherent detection has the advantages of low cost and optical field recovery ability. However, most of the self-coherent detection techniques are limited to single sideband (SSB) signals. Recently, carrier-assisted differential detection (CADD) has been proposed to realize complex-valued double sideband (DSB) signals, but it requires a high carrier-to-signal power ratio (CSPR) to mitigate the signal-to-signal beat interference (SSBI). Later, a more cost-effective symmetric CADD (S-CADD) has been proposed while the required CSPR is still high. In order to alleviate the high requirements of CSPR, we propose a scheme based on the joint of digital pre-distortion (DPD) at transmitter and clipping at receiver to further improve the S-CADD system performance. This joint processing can not only solve the problem of non-uniform distribution of subcarrier signal-to-noise ratio (SNR) caused by non-ideal transfer function, but also the error propagation problem caused by enhanced SSBI under low CSPR. After the validation of the 64 Gbaud 16-ary quadrature amplitude modulation (16-QAM) orthogonal frequency division multiplexing (OFDM) signal transmitted over 80 km standard single mode fiber (SSMF), the CSPR required by the proposed scheme to reach the 20% soft decision-forward error correction (SD-FEC) and 7% hard decision-forward error correction (HD-FEC) can be reduced by 1.3 dB and 2.8 dB, respectively, with a comparison of the conventional S-CADD. The results show the potential of the proposed scheme in the short-reach optical transmissions.

A syntactic features and interactive learning model for aspect-based sentiment analysis

The aspect-based sentiment analysis (ABSA) consists of two subtasks: aspect term extraction (AE) and aspect term sentiment classification (ASC). Previous research on the AE task has not adequately leveraged syntactic information and has overlooked the issue of multi-word aspect terms in text. Current researchers tend to focus on one of the two subtasks, neglecting the connection between the AE and ASC tasks. Moreover, the problem of error propagation easily occurs between two independent subtasks when performing the complete ABSA task. To address these issues, we present a unified ABSA model based on syntactic features and interactive learning. The proposed model is called syntactic interactive learning based aspect term sentiment classification model (SIASC). To overcome the problem of extracting multi-word aspect terms, the model utilizes part-of-speech features, words features, and dependency features as textual information. Meanwhile, we designs a unified ABSA structure based on the end-to-end framework, reducing the impact of error propagation issues. Interaction learning in the model can establish a connection between the AE task and the ASC task. The information from interactive learning contributes to improving the model’s performance on the ASC task. We conducted an extensive array of experiments on the Laptop14, Restaurant14, and Twitter datasets. The experimental results show that the SIASC model achieved average accuracy of 84.11%, 86.65%, and 78.42% on the AE task, respectively. Acquiring average accuracy of 81.35%, 86.71% and 76.56% on the ASC task, respectively. The SIASC model demonstrates superior performance compared to the baseline model.

Taut cables with hanging masses: A metamaterial-like dynamic behavior

We inspect the dynamic behavior of suspended cables presenting a periodic array of scatter elements, consisting of a discrete set of masses that are hanging by means of elastic or rigid connections. By introducing some approximations, we show that the problem in the continuous domain can be brought back to an equivalent discrete problem, whose solutions depend on the variation of an effective mass. We find that the essential spectrum of the problem presents band gaps. By considering boundary conditions, eigenmodes can only be found at frequencies belonging to pass bands for the propagation problem of the unbounded system. This is no more true when a defect of periodicity is inserted, leading to the formation of localized modes. We provide a quantitative validation of these theoretical results by means of comparison with experimental tests.

Full waveform inversion and Lagrange multipliers

SUMMARY Full-waveform inversion (FWI) is an effective method for imaging subsurface properties using sparsely recorded data. It involves solving a wave propagation problem to estimate model parameters that accurately reproduce the data. Recent trends in FWI have seen a renewed interest in extended methodologies, among which source extension methods leveraging reconstructed wavefields to solve penalty or augmented Lagrangian (AL) formulations have emerged as robust algorithms, even for inaccurate initial models. Despite their demonstrated robustness on synthetic data, challenges remain, such as the lack of a clear physical interpretation and reliance on difficult-to-compute least-squares (LS) wavefields. Moreover, the literature lacks a general and through comparison of these source extension methods with each other and with the standard FWI. This paper is divided into three critical parts. In the first, a novel formulation of these methods is explored within a unified Lagrangian framework. This novel perspective permits the introduction of alternative algorithms that use LS multipliers instead of wavefields. These multiplier-oriented variants appear as regularizations of the standard FWI, are suitable to the time domain, offer tangible physical interpretations, and foster enhanced convergence efficiency. The second part of the paper delves into understanding the underlying mechanisms of these techniques. This is achieved by solving the associated non-linear equations using iterative linearization and inverse scattering methods. The paper provides insight into the role and significance of Lagrange multipliers in enhancing the linearization of the equations. It explains how different methods estimate multipliers or make approximations to increase computing efficiency. Additionally, it presents a new physical understanding of the Lagrange multiplier used in the AL method, highlighting how important it is for improving algorithm performance when compared to penalty methods. In the final section, the paper presents numerical examples that compare different methods within a unified iterative algorithm, utilizing benchmark Marmousi and 2004 BP salt models.

A time-domain spectral finite element method for acoustoelasticity: Modeling the effect of mechanical loading on guided wave propagation

Ultrasonic testing techniques such as guided wave-based structural health monitoring aim to evaluate the integrity of a material with sensors and actuators that operate in situ, i.e. while the material is in use. Since ultrasonic wave propagation is sensitive to environmental conditions such as pre-deformation of the structure, the design and performance evaluation of monitoring systems in this context is a complicated task that requires quantitative data and the associated modeling effort. In our work, we propose a set of numerical tools to solve the problem of mechanical wave propagation in materials subjected to pre-deformation. This type of configuration is usually treated in the domain of acoustoelasticity. A relevant modeling approach is to consider two different problems: a quasi-static nonlinear problem for the large displacement field of the structure and a linearized time-domain wave propagation problem. After carefully reviewing the modeling ingredients to represent the configurations of interest, we propose an original combination of numerical tools that leads to a computationally efficient algorithm. More specifically, we use 3D shell elements for the quasi-static nonlinear problem and the time-domain spectral finite element method to numerically solve the wave propagation problem. Our approach can represent any type of material constitutive law, geometry or mechanical solicitation. We present realistic numerical results on 3D cases related to the monitoring of both isotropic and anisotropic materials, illustrating the genericity and efficiency of our method. We also validate our approach by comparing it to experimental data from the literature.

High-Efficiency and High-Precision Ship Detection Algorithm Based on Improved YOLOv8n

With the development of the intelligent vision industry, ship detection and identification technology has gradually become a research hotspot in the field of marine insurance and port logistics. However, due to the interference of rain, haze, waves, light, and other bad weather, the robustness and effectiveness of existing detection algorithms remain a continuous challenge. For this reason, an improved YOLOv8n algorithm is proposed for the detection of ship targets under unforeseen environmental conditions. In the proposed method, the efficient multi-scale attention module (C2f_EMAM) is introduced to integrate the context information of different scales so that the convolutional neural network can generate better pixel-level attention to high-level feature maps. In addition, a fully-concatenate bi-directional feature pyramid network (Concatenate_FBiFPN) is adopted to replace the simple superposition/addition of feature map, which can better solve the problem of feature propagation and information flow in target detection. An improved spatial pyramid pooling fast structure (SPPF2+1) is also designed to emphasize low-level pooling features and reduce the pooling depth to accommodate the information characteristics of the ship. A comparison experiment was conducted between other mainstream methods and our proposed algorithm. Results showed that our proposed algorithm outperformed other models by achieving 99.4% of accuracy, 98.2% of precision, 98.5% of recall, 99.1% of mAP@.5, and 85.4% of mAP@.5:.95 on the SeaShips dataset.

Free propagation of resonant waves in nonlinear dissipative metamaterials

This paper deals with the free propagation problem of resonant and close-to-resonance waves in one-dimensional lattice metamaterials endowed with nonlinearly viscoelastic resonators. The resonators' constitutive and geometric nonlinearities imply a cubic coupling with the lattice. The analytical treatment of the nonlinear wave propagation equations is carried out via a perturbation approach. In particular, after a suitable reformulation of the problem in the Hamiltonian setting, the approach relies on the well-known resonant normal form techniques from Hamiltonian perturbation theory. It is shown how the constructive features of the Lie Series formalism can be exploited in the explicit computation of the approximations of the invariant manifolds. A discussion of the metamaterial dynamic stability, either in the general or in the weak dissipation case, is presented.

Spectrally-tuned compact finite-difference schemes with domain decomposition and applications to numerical relativity

Compact finite-difference (FD) schemes specify derivative approximations implicitly, thus to achieve parallelism with domain-decomposition suitable partitioning of linear systems is required. Consistent order of accuracy, dispersion, and dissipation is crucial to maintain in wave propagation problems such that deformation of the associated spectra of the discretized problems is not too severe. In this work we consider numerically tuning spectral error, at fixed formal order of accuracy to automatically devise new compact FD schemes. Grid convergence tests indicate error reduction of at least an order of magnitude over standard FD. A proposed hybrid matching-communication strategy maintains the aforementioned properties under domain-decomposition. Under evolution of linear wave-propagation problems utilizing exponential integration or explicit Runge-Kutta methods improvement is found to remain robust. A first demonstration that compact FD methods may be applied to the Z4c formulation of numerical relativity is provided where we couple our header-only, templated C++ implementation to the highly performant GR-Athena++ code. Evolving Z4c on test-bed problems shows at least an order in magnitude reduction in phase error compared to FD for propagated metric components. Stable binary-black-hole evolution utilizing compact FD together with improved convergence is also demonstrated.

Local uncertainty quantification for 3-D time-domain full-waveform inversion with ensemble Kalman filters: application to a North Sea OBC data set

SUMMARY Full-waveform inversion (FWI) has emerged as the state-of-the art high resolution seismic imaging technique, both in seismology for global and regional scale imaging and in the industry for exploration purposes. While gaining in popularity, FWI, at an operational level, remains a heavy computational process involving the repeated solution of large-scale 3-D wave propagation problems. For this reason it is a common practice to focus the interpretation of the results on the final estimated model. This is forgetting FWI is an ill-posed inverse problem in a high dimensional space for which the solution is intrinsically non-unique. This is the reason why being able to qualify and quantify the uncertainty attached to a model estimated by FWI is key. To this end, we propose to extend at an operational level the concepts introduced in a previous study related to the coupling between ensemble Kalman filters (EnKFs) and FWI. These concepts had been developed for 2-D frequency-domain FWI. We extend it here to the case of 3-D time-domain FWI, relying on a source subsampling strategy to assimilate progressively the data within the Kalman filter. We apply our strategy to an ocean bottom cable field data set from the North Sea to illustrate its feasibility. We explore the convergence of the filter in terms of number of elements, and extract variance and covariance information showing which part of the model are well constrained and which are not. Analysing the variance helps to gain insight on how well the final estimated model is constrained by the whole FWI workflow. The variance maps appears as the superposition of a smooth trend related to the geometrical spreading and a high resolution trend related to reflectors. Mapping lines of the covariance (or correlation matrix) to the model space helps to gain insight on the local resolution. Through a wave propagation analysis, we are also able to relate variance peaks in the model space to variance peaks in the data space. Compared to other posterior-covariance approximation scheme, our combination between EnKF and FWI is intrinsically scalable, making it a good candidate for exploiting the recent exascale high performance computing machines.

Implicit-explicit time integration for the immersed wave equation

Efficient solution strategies for wave propagation problems with complex geometries are essential in many engineering fields. Immersed boundary methods simplify mesh generation by embedding the domain of interest into an extended domain that is easy to mesh, introducing the challenge of dealing with cells that intersect the domain boundary. We consider the finite cell method that extends the weak form from the physical into the extended domain, multiplied by a small value, to stabilize badly cut cells. Combined with explicit time integration methods, the presence of finite cells with very little support in the physical domain results in tiny critical time step sizes. While the finite cell stabilization limits how small the critical time step size can become, this limit is still restrictive for many applications. Explicit transient analyses commonly use the spectral element method due to its natural way of obtaining diagonal mass matrices through nodal lumping. The resulting method is very efficient since nodal lumping renders the solution of the equation systems trivial while not reducing the accuracy or the critical time step size. The combination of the spectral element and finite cell methods is called the spectral cell method. Unfortunately, a direct application of nodal lumping in the spectral cell method is impossible due to the special quadrature necessary to treat the discontinuous integrand inside the cut cells. The existing approaches to lump the mass matrices of cut cells significantly reduce the accuracy of the approximation.We analyze an implicit-explicit (IMEX) time integration method to exploit the advantages of the nodal lumping scheme for uncut cells on one side and the unconditional stability of implicit time integration schemes for cut cells on the other. In this hybrid, immersed Newmark IMEX approach, we use explicit second-order central differences to integrate the uncut degrees of freedom that lead to a diagonal block in the mass matrix and an implicit trapezoidal Newmark method to integrate the remaining degrees of freedom (those supported by at least one cut cell). The immersed Newmark IMEX approach preserves the high-order convergence rates and the geometric flexibility of the finite cell method while retaining the efficiency of the nodal lumping scheme for uncut cells without compromising the critical time step size. We analyze a simple system of spring-coupled masses to highlight some of the essential characteristics of Newmark IMEX time integration. We then solve the scalar wave equation with immersed boundaries on two- and three-dimensional examples with significant geometric complexity to show that our approach is more efficient than state-of-the-art time integration schemes when comparing accuracy and runtime. While we focus on the finite and spectral cell methods, we expect other immersed methods to benefit equally from Newmark IMEX time integration.

Recursive label attention network for nested named entity recognition

Nested named entity recognition (NER) has been mainly studied by recognizing inner entities first, then outer entities, or vice versa, depending on the nesting level. However, the previous models on nested NER do not deal with the data sparseness problems of deep nested entities (i.e., insufficiency of training data for deeper nested entities) and the error propagation problems between nested entities (i.e., accumulation of errors caused by inner entities or outer entities). To alleviate the data sparseness problems, we propose a recursive label attention network. In contrast to previous models, the proposed model explicitly reflects the nested level and effectively re-uses lower level label information through level-reflected label embeddings. To counteract error propagation in the recursive architecture, we propose an inner entity pretraining strategy in which the proposed model is sequentially trained from lower to higher levels. In experiments, we use task-specific metrics for nested NER, unlike the evaluation metrics used in previous studies. The experimental results show that the proposed model is particularly effective in cases where entities of the same type were nested. In addition, it showed high F1-scores even when named entities are deeply nested. This fact reveals that our proposed model is effective for nest NER.

SH waves in a weakly inhomogeneous half space with a nonlinear thin layer coating

We investigate the self-modulation of Love waves propagating in a nonlinear half-space covered by a nonlinear layer. We assume that the constituent material of the layer is nonlinear, homogeneous, isotropic, compressible, and hyperelastic, whereas for the half-space, it is nonlinear, heterogeneous, compressible and a different hyperelastic material. By employing the nonlinear thin layer approximation, the problem of wave propagation in a layered half-space is reduced to the one for a nonlinear heterogeneous half-space with a modified nonlinear homogeneous boundary condition on the top surface. This new problem is analyzed by a relevant perturbation method, and a nonlinear Schrödinger (NLS) equation defining the self-modulation of waves asymptotically is obtained. The dispersion relation is derived for different heterogeneous properties of the half-space and the thin layer. Then the results of the thin layer approximation are compared with the ones for the finite layer obtained in Teymur et al. (Int J Eng Sci 85:150–162, 2014). The solitary solutions of the derived NLS equation are obtained for selected real material models. It has been discussed how these solutions are influenced by the heterogeneity of the semi-infinite space.

Initiation and growth of fatigue cracks in sheets with U-shaped notches in the first and mixed modes of fracture

Initiation and growth of fatigue cracks in sheets with U-shaped notches in the first and mixed modes of fracture

A hydro-thermo-damage-mechanical fully coupled cohesive phase-field model for fracking in quasi-brittle thermo-poroelastic media

A hydro-thermo-damage-mechanical fully coupled cohesive phase-field method is developed for modelling of fracking in quasi-brittle thermo-poroelastic media. A new sigmoid interpolation function is proposed to characterize the transitional behaviour from porous media to cracks for permeability, specific heat capacity, and thermal conductivity. A fluid continuity equation with the Darcy-Poiseuille law and an energy conservation equation with the Fourier’s law are then used to simulate fluid flow and heat transfer in cracking porous media, respectively. The fluid pressure and temperature of fluids and solids are fully incorporated into the governing equation of the phase-field regularized cohesive zone model to efficiently simulate complicated quasi-brittle hydro-thermal fracking, without the need of remeshing, crack tracking or auxiliary fields. The resultant displacement–pressure-temperature-damage coupled multiphysics system of equations is solved using a staggered Newton–Raphson iterative algorithm within the finite element framework. The method is first validated by five two-dimensional thermally or hydro-thermally induced fracture examples with analytical solutions, experimental data, and published numerical solutions. Three-dimensional fracking and fluid circulation problems in a vertical five-spot wellbore system with complex natural fracture networks, are then simulated as an application to energy harvesting in enhanced geothermal systems (EGS). It is found that the present method is capable of accurately and robustly modelling complicated hydraulic fracture propagation and fluid circulation problems and has the potential to be used for fracking design and optimization of EGS.

Exact domain truncation for the Morse-ingard equations

Morse and Ingard [1] give a coupled system of time-harmonic equations for the temperature and pressure of an excited gas. These equations form a critical aspect of modeling trace gas sensors. Like other wave propagation problems, the computational problem must be closed with suitable far-field boundary conditions. Working in a scattered-field formulation, we adapt a nonlocal boundary condition proposed in [2] for the Helmholtz equation to this coupled system. This boundary condition uses a Green's formula for the true solution on the boundary, giving rise to a nonlocal perturbation of standard transmission boundary conditions. However, the boundary condition is exact and so Galerkin discretization of the resulting problem converges to the restriction of the exact solution to the computational domain. Numerical results demonstrate that accuracy can be obtained on relatively coarse meshes on small computational domains, and the resulting algebraic systems may be solved by GMRES using the local part of the operator as an effective preconditioner. These numerical results taken together combine several advanced techniques, including higher-order finite elements, geometric multigrid in curvilinear geometry, native use of complex arithmetic, and incorporation of nonlocal operators. These are tied together in a high-level simulation using the Firedrake library.

Convolution, deconvolution, the unit hydrograph and flood routing

Convolution equations are used to relate the input and the output of a system such as rainfall and runoff, or inflow and outflow of a river reach. There have been numerous reports of unsatisfactory results from the deconvolution necessary to calculate the connecting transfer function. The cause is that the equations are ill-conditioned, and it is shown here that the fundamental theoretical solution is that of wild oscillations such as has often been found computationally. A spectral method is proposed for numerical solution, where, instead of individual point values, the transfer function is expressed as series of given continuous functions, where the problem is to determine the coefficients of those functions. The resulting equations have been found to be well-conditioned, and solutions obtained were smooth, bounded, and enabled a certain amount of physical interpretation of the transfer function. The method has been applied to several problems, including typical rainfall–runoff ones and flood routing and wave propagation problems, with quite satisfactory results. Another problem for deconvolution is found to be the traditional use of truncated equations. A remedy is only to use later output data points where convolution with input data does not reach back beyond the initial one. For the routing of larger flood events, the linear methods employed were found to be not so accurate. However as they are a first approximation that requires no knowledge of stream geometry or resistance, and as either discharge or water level hydrographs can be used, they may be useful.

An Efficient Method for Solving Problems of Acoustic Scattering on Three-Dimensional Transparent Structures

The article contains a study of methods for solving integral equations in the context of acoustic problems. The methodology considered is applied to describe acoustic wave propagation and scattering. Efficient discretization methods are used together with iterative methods to solve the operator equations, including an apparatus for fast multiplication of the resulting post-discretization Toeplitz matrices by a vector using the fast Fourier transform. The theoretical analysis of the proposed numerical algorithm demonstrates its efficiency in terms of the required number of arithmetic operations and the memory footprint of the computing system. The presented numerical simulation demonstrates the possibility of solving the problem of acoustic wave propagation in transparent media using the proposed methods. A visualization of the obtained solutions for a practical problem with a high level of discretization of the solution volume domain is also presented.

An Efficient Local Artificial Boundary Condition for Infinite Long Finite Element Euler–Bernoulli Beam

To solve the wave propagation problems of the Euler–Bernoulli beam in an unbounded domain effectively and efficiently, a new local artificial boundary condition technology is proposed. It replaces the residual right-hand side of the truncated discrete equation with an equivalent linear algebraic system. First, the equivalent Schrodinger equation is discussed. Its artificial boundary condition is obtained by first rationalizing the Dirichlet-to-Neumann condition in the frequency domain with a Pade approximation and then inverse transforming each Pade term back into the time domain by introducing auxiliary degrees of freedom. Frequency shifting is employed such that it performs better near a prescribed frequency. Then, the artificial boundary condition of the finite element Euler–Bernoulli beam is obtained by simple algebraic manipulations on that of the corresponding Schrodinger equation. This method only makes local changes to the original truncated discrete dynamic system and thus is very efficient and easy to use. The accuracy of the proposed method can be improved by using more Pade terms and a proper shift frequency. The numerical example shows, with only a few additional degrees of freedom, the proposed artificial boundary condition effectively eliminates the spurious reflection. The idea of the proposed method can also be used in other dispersive wave systems.