Simulation Of Propagation Research Articles

A smooth particle dynamics method considering material point failure for crack propagation simulation applied in rock mechanics

A smooth particle dynamics method considering material point failure for crack propagation simulation applied in rock mechanics

Full-Vectorial Light Propagation Simulation of Optimized Beams in Scattering Media

Volumetric scattering prevents imaging modalities in biomedical optics from imaging deep inside tissue. The optimization of the incident wavefront has the potential to improve these imaging modalities. To investigate the optimization and light propagation of such beams inside scattering media rigorously, full-vectorial simulations based on solutions of Maxwell’s equations are necessary. In this publication, we present a versatile two-step beam synthesis method to efficiently simulate the scanning and phase optimization of a focused beam inside a static scattering medium. We present four different approaches to the phase optimization of the energy density and the absolute value of the Poynting vector. We find that these quantities have two regions with different, almost exponential decays over depth for a non-optimized beam. Optimization by conjugating the phase of the projected electric field in various directions at the focus shows an improvement below a certain penetration depth. Seeking global solutions to the optimization problems reveals an even better enhancement in the energy density and the absolute value of the Poynting vector in the focus. For Poynting vector optimization, the differences between the presented optimization approaches are more significant than for the energy density. With the presented method, it is possible to efficiently simulate different imaging methods improved by wavefront shaping to investigate their possible penetration depths.

On dynamic crack propagation in a lattice Boltzmann method for elastodynamics in 2D

AbstractIn recent years, the development of lattice Boltzmann methods (LBMs) for solids has gained attention. Fracture mechanics as a viable application for these methods has been presented before, albeit for mode III cracks in the context of an LBM for antiplane shear deformation. The performance of the LBM itself is promising, while the usage of a regular lattice simplifies the modelling of fractures significantly. Recent advancements in LBMs for solids, especially the description of Dirichlet‐ and Neumann‐type boundary conditions, now make it possible to extend the LBM simulation of crack propagation to the plane strain case with modes I and II crack opening, including growth with non‐uniform speed in arbitrary directions. For this, the configurational force acting on a crack tip is utilised. The definition of the moments of the LBM, which are based on the balance laws of continuum mechanics, render the evaluation of macroscopic fields in the configuration straightforward. In this work, the general in‐plane case of dynamic crack propagation is shown and necessary considerations for the implementation are discussed. Lastly, numerical examples showcase the capabilities of the proposed method to model dynamic fractures and establish a proof‐of‐concept.

Spectral-infinite-element simulations of seismic wave propagation in self-gravitating, rotating 3-D Earth models

SUMMARY Although observation of gravity perturbations induced by earthquakes is possible, simulation of seismic wave propagation in a self-gravitating, rotating Earth model with 3-D heterogeneity is challenging due to the numerical complexities associated with the unbounded Poisson/Laplace equation that governs gravity perturbations. Therefore, gravity perturbations are generally omitted, and only the background gravity is taken into account using the so-called Cowling approximation. However, gravity perturbations may be significant for large earthquakes (Mw ≥ 6.0) and long-period responses. In this study, we develop a time-domain solver based on the spectral-infinite-element approach, which combines the spectral element method inside the Earth domain with a mapped-infinite-element method in the infinite space outside. This combination allows us to solve the complete, coupled momentum-gravitational equations in a fully discretized domain while accommodating complex 3-D Earth models. We compute displacement and gravity perturbations considering various Earth models, including Preliminary Reference Earth Model and S40RTS and conduct comprehensive benchmarks of our method against the spherical harmonics normal-mode approach and the direct radial integration method. Our 3-D simulations accommodate topography, bathymetry, rotation, ellipticity and oceans. Results show that our technique is accurate and stable for long simulations. Our method provides a new scope for incorporating earthquake-induced gravity perturbations into source and adjoint tomographic inversions.

Effect of intermittent joint distribution on the mechanical and acoustic behavior of rock masses

The mechanical characteristics and acoustic behavior of rock masses are greatly influenced by stochastic joints. In this study, numerical models of rock masses incorporating intermittent joints with different numbers and dip angles were produced using the finite element method (FEM) with the intrinsic cohesive zone model (ICZM). Then, the uniaxial compressive and wave propagation simulations were performed. The results indicate that the joint number and dip angle can affect the mechanical and acoustic properties of the models. The uniaxial compressive strength (UCS) and wave velocity of rock masses decrease monotonically as the joint number increases. However, the wave velocity grows monotonically as the joint dip angle increases. When the joint dip angle is 45°–60°, the UCS of the rock mass is lower than that of other dip angles. The wave velocity parallel to the joints is greater than that perpendicular to the joints. When the dip angle of joints remains unchanged, the UCS and wave velocity are positively related. When the joint dip angle increases, the variation amplitude of the UCS regarding the wave velocity increases. To reveal the effect of the joint distribution on the velocity, a theoretical model was also proposed. According to the theoretical wave velocity, the change in wave velocity of models with various joint numbers and dip angles was consistent with the simulation results. Furthermore, a theoretical indicator (i.e. fabric tensor) was adopted to analyze the variation of the wave velocity and UCS.

Locating Three-Dimensional Position of Deep-Seated SERS Phantom Lesions in Thick Tissues Using Tomographic Transmission Raman Spectroscopy.

Locating distinct objects within a thick scattering medium remains a long-standing challenge in the fields of materials science, health, and engineering. Transmission Raman spectroscopy (TRS) with the use of surface-enhanced Raman scattering (SERS) nanoparticles has proven to be an effective approach to detect deep-seated lesions inside thick biological tissues. However, it has not yet been proven to spatially locate deep lesions in three dimensions using optical modalities. Herein, we present the concept of tomographic TRS and report its successful use for accurately locating SERS nanoparticles in elongated rod-like thick tissues. Our work starts with theoretical simulations of Raman photon propagation in tissues. We discovered a linear relationship between the Raman spectral peak ratio and propagation distance of Raman photons in tissues, allowing us to predict the location of lesions tagged by SERS NPs. Based on this, we propose a two-step tomographic TRS strategy, which includes axial scanning and ring scanning. We demonstrate the robustness of our approach using ex vivo thick tissue (4.5 cm in thickness) and locate an embedded SERS phantom lesion, with a ring scanning step of 10-30°. We successfully locate multiple SERS phantom lesions in the ex vivo porcine muscle stack with high accuracy (absolute error of <2 mm). Our method is rapid, efficient, and of low cost compared to current tomographic medical imaging techniques. This work advances Raman techniques for three-dimensional positioning and offers new insights toward practical diagnosis applications.

Dynamic and adaptive mesh-based graph neural network framework for simulating displacement and crack fields in phase field models

Fracture is one of the main causes of failure in engineering structures. Phase field methods coupled with adaptive mesh refinement (AMR) techniques have been widely used to model crack propagation due to their ease of implementation and scalability. However, phase field methods can still be computationally demanding making them unfeasible for high-throughput design applications. Machine learning (ML) models such as Graph Neural Networks (GNNs) have shown their ability to emulate complex dynamic problems with speed-ups orders of magnitude faster compared to high-fidelity simulators. In this work, we present a dynamic mesh-based GNN framework for emulating phase field simulations of single-edge crack propagation with AMR for different crack configurations. The developed framework – ADAPTive mesh-based graph neural network (ADAPT-GNN) – exploits the benefits of both ML methods and AMR by describing the graph representation at each time step as the refined mesh itself. Using ADAPT-GNN, we predict the evolution of displacement fields and scalar damage field (or phase field) with good accuracy compared to the conventional phase-field fracture model. We also compute crack stress fields with good accuracy using the predicted displacements and phase field parameter.

Study on the Filter Integrated with High Transmission and Ultra-Wideband Electromagnetic Shielding Function in the View of Surface Plasmon

The generation of a strong electromagnetic weapon requires the electromagnetic protection performance of the photoelectric system, including the high shielding and ultra-wideband protection characteristics against strong electromagnetic attacks and superhigh infrared transmittance. This article aims to excite surface plasmons by etching a subwavelength circular structure array run through a metal film. At the same time, the local field enhancement effect based on surface plasmons is used to tune specific wavelengths of light, achieving the goal of enhancing the optical transmission effect. Ulteriorly, the amplitude and position of the infrared response peak are controlled by optimizing the structural geometric parameters to achieve filtering effects in different wavebands. Furthermore, through the inversion calculation of the transmission spectrum of the circle array structure in the visible light band, the different colors of its surface under different parameters can correspond one-to-one to its performance, so as to intuitively identify different filter types and their performance. Meanwhile, the remaining metal film after etching off the circle array structure region constitutes an ultra-wideband shielding layer as a continuous conductor, achieving a shielding efficiency of more than 45 dB in the radar band (1–18 GHz). It is worth noting that the cascaded theory combines ultra-thin metal films with thick substrates to solve the problem of software running for long periods of time and under high loads during the simulation of light propagation. This scheme greatly reduces the error between simulation and practical application.

Crack Growth Simulations With SMART Method in Adhesively Bonded Joints

Due to the important advantages of adhesive joints, such as their suitability for multi-material designs, their use has been increasing in the last decade. Determining the fracture behavior of structural adhesive bondings is essential for structural durability. In crack propagation analyses, adaptive meshing has drawn considerable attention because of its improvements in terms of complex preprocessing and time management. This paper presents a recently introduced separating morphing and adaptive remeshing technology (SMART) innovative crack growth simulation for adhesively bonded joints, considering static and cyclic cases. For the static case, an R-curve was obtained for the bonding joints of carbon steel and Araldite 2015. For the cyclic case, the Carbon Fiber Reinforced Polymer (CFRP) bonding joints were analyzed under constant-amplitude loading conditions. The crack-propagation rate and the number of cycles were estimated. Crack propagation simulations were validated using experimental data. Acceptable agreement was achieved between the experimental and estimated results.

Peridynamic simulation of heat transfer during quenching of semi-solid plate with occurrence of hot cracks

Quenching plays an important role in the material processing. Water impingement jet quenching involves phenomena such as nucleate boiling, Leidenfrost effect, rewetting, etc. Stronger thermal gradients could induce cracks in the material. Thermal boundary conditions of the process depends on temperature, space, time, and material discontinuity. This study mainly focuses on the thermal modeling of moving water jet quenching of semi-solid hot vertical aluminum plate with latent heat release during the solidification by using Peridynamics (PD). The obtained thermal profile is compared with the infrared camera’s thermal field data to identify the boiling curve parameters on the cracked plate. The PD model with the incorporation of cracking dynamics is giving agreeable results with the experimental thermal field. Modeling quenching heat transfer using PD will enhance the coupling of hot tear origination and propagation simulation, which is quite challenging with other numerical methods. The application of PD in quenching heat transfer with phase change is the first of its kind. The advantages and limitations of the peridynamic approach for this problem is critically evaluated and the simulated results are compared with the experimental measurement of temperature profiles.

3D Wave Propagation Simulations of Mw 6.5+ Earthquakes on the Tacoma Fault, Washington State, Considering the Effects of Topography, a Geotechnical Gradient, and a Fault Damage Zone

ABSTRACT We simulate shaking in Tacoma, Washington, and surrounding areas from Mw 6.5 and 7.0 earthquakes on the Tacoma fault. Ground motions are directly modeled up to 2.5 Hz using kinematic, finite-fault sources; a 3D seismic velocity model considering regional geology; and a model mesh with 30 m sampling at the ground surface. In addition, we explore how adjustments to the seismic velocity model affect predicted shaking over a range of periods. These adjustments include the addition of a region-specific geotechnical gradient, surface topography, and a fault damage zone. We find that the simulated shaking tends to be near estimates from empirical ground-motion models (GMMs). However, long-period (T = 5.0 s) shaking within the Tacoma basin is typically underpredicted by the GMMs. The fit between simulated and GMM-derived short-period (T = 0.5 s) shaking is significantly improved with the addition of the geotechnical gradient. From comparing different Mw 6.5 earthquake scenarios, we also find that the response of the Tacoma basin is sensitive to the azimuth of incoming seismic waves. In adding surface topography to the simulation, we find that average ground motion is similar to that produced from the nontopography model. However, shaking is often amplified at topographic highs and deamplified at topographic lows, and the wavefield undergoes extensive scattering. Adding a fault damage zone has the effect of amplifying short-period shaking adjacent to the fault, while reducing far-field shaking. Intermediate-period shaking is amplified within the Tacoma basin, likely due to enhanced surface-wave generation attributable to the fault damage zone waveguide. When applied in the same model, the topography and fault damage zone adjustments often enhance or reduce the effects of one another, adding further complexity to the wavefield. These results emphasize the importance of improving near-surface velocity model resolution as waveform simulations progress toward higher frequencies.

Development of early warning system for tsunamis accompanied by collapse of Anak Krakatau volcano, Indonesia

Present tsunami warning systems have been specialized for earthquake-generated tsunamis, but rapidly evaluating the tsunamis caused by volcanic eruptions and/or volcanic sector collapses remains a challenge. In this study, we applied a numerical model to the 2018 Anak Krakatau tsunami event, which was generated by the sector collapse, investigated a tsunami prediction skill by the model, and developed a real-time forecasting method based on a pre-computed database for future tsunamis accompanied by such eruption of Anak Krakatau. The database stores spatiotemporal changes in water surface level and flux, which are simulated under various collapse scenarios, for confined areas in the vicinity of potential source. The areas also cover the locations of observation stations that are virtually placed on uninhabited island surrounding the source area. During an actual volcanic tsunami event, a tsunami is expected to be observed at the observation stations. For real-time tsunami forecasting, the most suitable scenarios to reproduce the observed waveforms are searched quickly in the database. The precomputed results under the identified scenarios are further provided as input for rapid tsunami propagation simulation. Therefore, an effective real-time forecasting can be conducted to densely populated coastal areas located at a considerable distance from the source, such as the coasts of Java and Sumatra. The forecasting performance was examined by applying the method for three hypothetical collapse scenarios assuming different sliding directions. We demonstrated that the tsunamis along the coasts were successfully forecasted. Moreover, we showed that the combination of a pre-computed database and the existence of observation stations near the source area was able to produce appropriate tsunami forecasting for the coastal area even in a volcanic event.

Study on sound transmission characteristics and positioning strategy of gas insulated transmission line shell

Gas insulated transmission line (GIL) is widely used for large-capacity power transmission in complex environment, and has the characteristics of long distance and totally enclosed structure. When the insulation fault occurs inside GIL, its metal shell structure makes it more difficult to monitor through optical, electrical and other physical quantities. Considering the rapidity and accuracy of positioning, the propagation of sound waves along the GIL sheath is a commonly used method in GIL on-site PD testing and running state online monitoring for fault positioning. This paper establishes a high-pressure GIL acoustic propagation simulation model based on pipe acoustics theory, analyzes the frequency dispersion and mode distribution of acoustic waves along the GIL shell, determines an appropriate fault location strategy based on the propagation velocity and attenuation of the selected acoustic wave mode, and carries out experiments to verify it. The results of the research show that the low order bending wave F(1,1) propagated along the GIL shell has good positioning characteristics, and the propagation velocity conforms to the dispersion curve distribution. By taking into account the attenuation and delay of the acoustic wave passing through epoxy insulators or expansion joints in the fault localization strategy, the localization accuracy can be improved by 60%. And this positioning strategy was verified to be effective in the GIL insulation fault localization test conducted in the field.

Conflation of modal strain energy, spectral finite element, and machine learning techniques for damage detection in multi-span steel girder bridges with variable sections

A new approach is proposed to detect micro and small damages in multi-span bridges constructed with variable sections by combining modal strain energy and spectral finite element through the use of machine learning techniques. The traditional finite element methods are inefficient in detecting and estimating structures' micro-damages. Therefore, wave propagation-based methods must be applied. Micro-damages can be modeled through the spectral finite element method, widely adopted in wave propagation simulation. Among the frequency-based methods used in damage identification, the modal strain energy is more sensitive, especially in complicated structures. However, this method is difficult to apply to steel girders with variable sections and different support conditions. By combining the modal strain energy and the spectral finite element and applying the support vector regression method, the location and severity of damage in steel girders could be accurately identified to overcome the difficulties of detection in complicated elements. To improve the ability of the support vector regression method, in the damage severity estimation, it was combined with feature selection. To evaluate the proposed method's ability, first, an element of the beam- type, with the characteristics of a spectral finite element, is defined in the OpenSees software. To assess the accuracy of the defined element in detecting micro-damages, a tested numerical model was adopted to measure the accuracy of the provided element. Next, to detect small damages, the model of the steel girder bridge with the defined element was prepared, and the modal properties of the structure were extracted. Due to the change in cross-section along each girder, the modified modal strain energy-based damage index (β) equation was used in the subject model. By applying the damages of different severities in steel girders in single and multiple damage scenarios, the modified damage index β was calculated. With the damage indexes calculated in the different scenarios and applying the machine learning method, the severity of damage was estimated. The obtained results indicate that this proposed method is accurate and satisfactory for detecting and estimating micro and small damages with high accuracy in bridges with variable section girders.

Time-domain immersed-boundary simulation of acoustic propagation between two spherical gas bubblesa).

Acoustic fields and resonance responses from two spherical gas bubbles are investigated using a time-domain simulation. Interior acoustic fields are obtained simultaneously from the simulation of the entire acoustic field propagation with an immersed-boundary method. The linear resonance responses are obtained and discussed for each of the bubbles with respect to the respective interior gas velocities. Also, these are analyzed in the frequency domain in terms of the coupled interactions. Unlike previous numerical and analytical solutions, the method allows for two bubbles of different sizes and shapes to be in contact with each other, which is representative of applicable underwater scattering targets.

Recent progress in computational predictions on fatigue and fracture resistance of steel connection joints

AbstractSteel structures are often subjected to cyclic loads, and sometimes to monotonic overloads, during the technical lifetime, which causes their performance degradation. Incidentally, they undergo progressive collapse, especially due to local failure of connections or joints. To ensure the overall integrity of the steel structures and to make sure the efficient transfer of the load between different components, it is necessary to evaluate the ultimate capacity and fatigue resistance of steel connections. Advances in material modelling and well‐controlled experiments at various scales in recent years have increased confidence in the finite element method as an effective way to predict the ultimate resistance and fatigue life. On examples of two common steel connections: bolted connection and welded connection, the recent progress will be presented in ductile fracture and fatigue crack initiation and propagation simulation.

LANDSLIDE TSUNAMI HAZARD ASSESSMENT: A NUMERICAL MODEL FOR THE SIMULATION OF MULTIPLE LANDSLIDE-INDUCED TSUNAMIS SCENARIOS IN A MONTE CARLO FRAMEWORK

Submarine landslides can pose serious tsunami hazard to coastal communities. However, performing a comprehensive landslide tsunami hazard assessment for a given area is in general difficult in view of the large uncertainties associated with tsunamigenic source parameters. These are often known only approximately, on the basis of estimates of the landslide geometry, slide material properties, and kinematics. Here, we present an efficient model to perform such hazard analyses, based on solving the linear mildslope equation with a time-dependent source term representing the seafloor motion as detailed in Iorio et al. (2021). This approach allows carrying out many computations, for a large number of landslide scenarios, in a Monte Carlo (MC) approach framework, at a reduced computational cost compared to other available methods, while still providing physically accurate simulations of most landslide tsunami generation and propagation.

Simulation of flexural-gravity wave propagation for elastic plates in shallow water using an energy-stable finite difference method with weakly enforced boundary and interface conditions

We introduce an energy stable, high-order-accurate finite difference approximation of the dynamic, pure bending Kirchhoff plate equations for complex geometries and spatially variable properties. We utilize the summation-by-parts (SBP) framework to discretize the biharmonic operator with variable coefficients, with attention given to free and clamped boundary conditions and corner conditions. Energy conservation is established by combining SBP boundary closures with weak enforcement of the boundary and interface conditions using a penalty (simultaneous approximation term, SAT) technique. Then we couple the plate equations to the shallow water equations to study flexural-gravity wave propagation, and prove that the semi-discrete system of equations is self-adjoint. We demonstrate the stability and accuracy properties of the method on curvilinear multiblock grids using the method of manufactured solutions. The method, which we provide in an open-source code, is then used to model ocean wave interactions with the Thwaites Glacier and Pine Island Ice Shelf in the Amundsen Sea off the coast of West Antarctica.

An immersed peridynamics model of fluid-structure interaction accounting for material damage and failure

This paper develops and benchmarks an immersed peridynamics method to simulate the deformation, damage, and failure of hyperelastic materials within a fluid-structure interaction framework. The immersed peridynamics method describes an incompressible structure immersed in a viscous incompressible fluid. It expresses the momentum equation and incompressibility constraint in Eulerian form, and it describes the structural motion and resultant forces in Lagrangian form. Coupling between Eulerian and Lagrangian variables is achieved by integral transforms with Dirac delta function kernels, as in standard immersed boundary methods. The major difference between our approach and conventional immersed boundary methods is that we use peridynamics, instead of classical continuum mechanics, to determine the structural forces. We focus on non-ordinary state-based peridynamic material descriptions that allow us to use a constitutive correspondence framework that can leverage well-characterized nonlinear constitutive models of soft materials. The convergence and accuracy of our approach are compared to both conventional and immersed finite element methods using widely used benchmark problems of nonlinear incompressible elasticity. We demonstrate that the immersed peridynamics method yields comparable accuracy with similar numbers of structural degrees of freedom for several choices of the size of the peridynamic horizon. We also demonstrate that the method can generate grid-converged simulations of fluid-driven material damage growth, crack formation and propagation, and rupture under large deformations.

A hybrid PML formulation for the 2D three-field dynamic poroelastic equations

Simulation of wave propagation in poroelastic half-spaces presents a common challenge in fields like geomechanics and biomechanics, requiring Absorbing Boundary Conditions (ABCs) at the semi-infinite space boundaries. Perfectly Matched Layers (PML) are a popular choice due to their excellent wave absorption properties. However, PML implementation can lead to problems with unknown stresses or strains, time convolutions, or PDE systems with Auxiliary Differential Equations (ADEs), which increases computational complexity and resource consumption.This article introduces a novel hybrid PML formulation for arbitrary poroelastic domains. Instead of using ADEs, this formulation utilizes time-history variables to reduce the number of unknowns and mathematical operations. The modification of the PDEs to account for the PML is limited to the PML domain only, resulting in smaller matrices while maintaining the governing equations in the interior domain and preserving the temporal structure of the problem. The hybrid approach introduces three scalar variables localized within the PML domain.The proposed formulation was tested in three numerical experiments in geophysics using realistic parameters for soft sites with free surfaces. The results were compared with numerical solutions from extended domains and simpler ABCs, such as paraxial approximation, demonstrating the accuracy, efficiency, and precision of the proposed method. The article also discusses the applicability of the method to complex media and its extension to the Multiaxial PML formulation.The codes used for the simulations are available for download from https://github.com/hmella/POROUS-HYBRID-PML.