Consistency between domain wall oscillation modes and spin wave modes in nanostrips

Investigations on domain wall(DW) and spin wave(SW) modes in a series of nanostrips with different widths and thicknesses have been carried out using micromagnetic simulation. The simulation results show that the frequencies of SW modes and the corresponding DW modes are consistent with each other if they have the same node number along the width direction. This consistency is more pronounced in wide and thin nanostrips, favoring the DW motion driven by SWs.Further analysis of the moving behavior of a DW driven by SWs is also carried out. The average DW speed can reach a larger value of ~ 140 m/s under two different SW sources. We argue that this study is very meaningful for the potential application of DW motion driven by SWs.

A Time-Domain Numerical Simulation for Free Motion Responses of Two Ships Advancing in Head Waves

The constant panel method within the framework of potential flow theory in the time domain is developed for solving the hydrodynamic interactions between two parallel ships with forward speed.When solving problems within a time domain framework,the free water surface needs to simultaneously satisfy both the kinematic and dynamic boundary conditions of the free water surface.This provides conditions for adding artificial damping layers.Using the Runge−Kutta method to solve equations related to time.An upwind differential scheme is used in the present method to deal with the convection terms on the free surface to prevent waves upstream.Through the comparison with the available experimental data and other numerical methods,the present method is proved to have good mesh convergence,and satisfactory results can be obtained.The constant panel method is applied to calculate the hydrodynamic interaction responses of two parallel ships advancing in head waves.Numerical simulations are conducted on the effects of forward speed,different longitudinal and lateral distances on the motion response of two modified Wigley ships in head waves.Then further investigations are conducted on the effects of different ship types on the motion response.

Electronic states of domain walls in commensurate charge density wave ground state and mosaic phase in 1T-TaS_(2)

Domain walls(DWs)in the charge-density-wave(CDW)Mott insulator 1T-TaS_(2)have unique localized states,which play an important role in exploring the electronic properties of the material.However,the electronic states in DWs in 1TTaS_(2)have not been clearly understood,mostly due to the complex structures,phases,and interlayer stacking orders in the DW areas.Here,we explored the electronic states of DWs in the large-area CDW phase and mosaic phase of 1T-TaS_(2)by scanning tunneling spectroscopy.Due to the different densities of DWs,the electronic states of DWs show distinct features in these phases.In the large area CDW phase,both the domain and the DWs(DW1,DW2,DW4)have zero conductance at the Fermi level;while in the mosaic phase,they can be metallic or insulating depending on their environments.In areas with a high density of DWs,some electronic states were observed both on the DWs and within the domains,indicating delocalized states over the whole region.Our work contributes to further understanding of the interplay between CDW and electron correlations in 1T-TaS_(2).

Preserved amplitude migration based on the one way wave equation in the angle domain

Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proportional to reflection coefficients. In recent years, pre-stack depth migration algorithms which preserve amplitudes and based on the one- way wave equation have been developed. Using the method in the shot domain requires a deconvolution imaging condition which produces some instability in areas with complicated structure and dramatic lateral variation in velocity. Depth migration with preserved amplitude based on the angle domain can overcome the instability of the one-way wave migration imaging condition with preserved amplitude. It can also offer provide velocity analysis in the angle domain of common imaging point gathers. In this paper, based on the foundation of the one-way wave continuation operator with preserved amplitude, we realized the preserved amplitude prestack depth migration in the angle domain. Models and real data validate the accuracy of the method.

Multi-scale seismic full waveform inversion in the frequency-domain with a multi-grid method

Although full waveform inversion in the frequency domain can overcome the local minima problem in the time direction, such problem still exists in the space direction because of the media subsurface complexity. Based on the optimal steep descent methods, we present an algorithm which combines the preconditioned bi-conjugated gradient stable method and the multi-grid method to compute the wave propagation and the gradient space. The multiple scale prosperity of the waveform inversion and the multi-grid method can overcome the inverse problems local minima defect and accelerate convergence. The local inhomogeneous three-hole model simulated results and the Marmousi model certify the algorithm effectiveness.

A modified domain reduction method for numerical simulation of wave propagation in localized regions

A modified domain reduction method(MDRM) that introduces damping terms to the original DRM is presented in this paper. To verify the proposed MDRM and compare the computational accuracy of these two methods, a numerical test is designed. The numerical results of the MDRM and DRM are compared using an extended meshed model. The results show that the MDRM significantly improved the computational accuracy of the DRM. Then, the MDRM is compared with two existing conventional methods, namely Liao's transmitting boundary and viscous-spring boundary with Liu's method. The MDRM shows its great advancement in computational accuracy, stability and range of applications. This paper also discusses the influence of boundary location on computational accuracy. It can be concluded that smaller models tend to have larger errors. By introducing two dimensionless parameters, φ_1 and φ_2, the rational distance between the observation point and the MDRM boundary is suggested. When φ_1 >2 or φ_2>13, the relative PGA error can be limited to 5%. In practice, the appropriate model size can be chosen based on these two parameters to achieve desired computational accuracy.

Decoupling coefficients of dilatational wave for Biot's dynamic equation and its Green's functions in frequency domain

Green＇s functions for Blot＇s dynamic equation in the frequency domain can be a highly useful tool for the investigation of dynamic responses of a saturated porous medium. Its applications are found in soil dynamics, seismology, earthquake engineering, rock mechanics, geophysics, and acoustics. However, the mathematical work for deriving it can be daunting. Green＇s functions have been presented utilizing an analogy between the dynamic thermoelasticity and the dynamic poroelasticity in the frequency domain using the u-p formulation. In this work, a special term ＂decoupling coefficient＂ for the decomposition of the fast and slow dilatational waves is proposed and expressed to present a new methodology for deriving the poroelastodynamic Green＇s functions. The correct- ness of the solution is demonstrated by numerically comparing the current solution with Cheng＇s previous solution. The separation of the two waves in the present methodology allows the more accurate evaluation of Green＇s functions, particularly the solution of the slow dilatational wave. This can be advantageous for the numerical implementation of the boundary element method （BEM） and other applications.

A Fully Nonlinear HOBEM with the Domain Decomposition Method for Simulation of Wave Propagation and Diffraction

A higher-order boundary element method(HOBEM) for simulating the fully nonlinear regular wave propagation and diffraction around a fixed vertical circular cylinder is investigated. The domain decomposition method with continuity conditions enforced on the interfaces between the adjacent sub-domains is implemented for reducing the computational cost. By adjusting the algorithm of iterative procedure on the interfaces, four types of coupling strategies are established, that is, Dirchlet/Dirchlet-Neumman/Neumman(D/D-N/N), Dirchlet-Neumman(D-N),Neumman-Dirchlet(N-D) and Mixed Dirchlet-Neumman/Neumman-Dirchlet(Mixed D-N/N-D). Numerical simulations indicate that the domain decomposition methods can provide accurate results compared with that of the single domain method. According to the comparisons of computational efficiency, the D/D-N/N coupling strategy is recommended for the wave propagation problem. As for the wave-body interaction problem, the Mixed D-N/N-D coupling strategy can obtain the highest computational efficiency.

Time-domain wavefield reconstruction inversion

Wave field reconstruction inversion （WRI） is an improved full waveform inversion theory that has been proposed in recent years. WRI method expands the searching space by introducing the wave equation into the objective function and reconstructing the wavefield to update model parameters, thereby improving the computing efficiency and mitigating the influence of the local minimum. However, frequency-domain WRI is difficult to apply to real seismic data because of the high computational memory demand and requirement of time-frequency transformation with additional computational costs. In this paper, wavefield reconstruction inversion theory is extended into the time domain, the augmented wave equation of WRI is derived in the time domain, and the model gradient is modified according to the numerical test with anomalies. The examples of synthetic data illustrate the accuracy of time-domain WRI and the low dependency of WRI on low-frequency information.

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability.

Surface wave attenuation based polarization attributes in time-frequency domain for multicomponent seismic data

In the paper, we propose a surface wave suppression method in time-frequency domain based on the wavelet transform, considering the characteristic difference of polarization attributes, amplitude energy and apparent velocity between the effective signals and strong surface waves. First, we use the proposed method to obtain time-frequency spectra of seismic signals by using the wavelet transform and calculate the instantaneous polarizability at each point based on instantaneous polarization analysis. Then, we separate the surface wave area from the signal area based on the surface-wave apparent velocity and the average energy of the signal. Finally, we combine the polarizability, energy, and frequency characteristic to identify and suppress the signal noise. Model and field data are used to test the proposed filtering method.

Second-Order Wave Diffraction Around 3-D Bodies by A Time-Domain Method

A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.

A NON-OVERLAPPING DOMAIN DECOMPOSITION ALGORITHM BASED ON THE NATURAL BOUNDARY REDUCTION FOR WAVE EQUATIONS IN AN UNBOUNDED DOMAIN

In this paper, a new domain decomposition method based on the natural boundary reduction, which solves wave problems over an unbounded domain, is suggestted. An circular artifcial boundary is introduced. The original unbounded domain is divided into two subdomains, an internal bounded region and external unbounded region outside the artificial boundary. A Dirichlet-Neumann(D-N) alternating iteration algorithm is constructed. We prove that the algorithm is equavilent to preconditional Richardson iteration method. Numerical studies are performed by finite element method. The numerical results show that the convergence rate of the discrete D-N iteration is independent of the fnite element mesh size.

GLOBAL EXISTENCE FOR THE NONHOMOGENEOUS QUASILINEAR WAVE EQUATION WITH A LOCALIZED WEAKLY NONLINEAR DISSIPATION IN EXTERIOR DOMAINS

It is proved that the global existence for the nonhomogeneous quasilinear wave equation with a localized weakly nonlinear dissipation in exterior domains.

Parallel Simulation of 3D Wave Propagation by Domain Decomposition

In order to perform large scale numerical simulation of wave propagation in 3D heterogeneous multiscale viscoelastic media, Finite Difference technique and its parallel implementation based on domain decomposition is used. A couple of typical statements of borehole geophysics are dealt with—sonic log and cross well measurements. Both of them are essentially multiscales, which claims to take into account heterogeneities of very different sizes in order to provide reliable results of simulations. Locally refined spatial grids help us to avoid the use of redundantly tiny grid cells in a target area, but cause some troubles with uniform load of Processor Units involved in computations. We present results of scalability tests together with results of numerical simulations for both statements performed for some realistic models.

Spin waves and transverse domain walls driven by spin waves: Role of damping

Based on the uniform,helical and spiral domain-wall magnetic configurations,the excited spin waves are studied with emphasis on the role of damping.We find that the damping closes the gap of dispersion,and greatly influences the dispersion in the long-wave region for the spin waves of spiral wall and helical structure.For the uniform configuration,the Dzyaloshinskii-Moriya interaction determines the modification of dispersion by the damping.Furthermore,we investigate the interaction between spin waves and a moving spiral domain wall.In the presence of damping,the amplitude of spin wave can increase after running across the wall for small wave numbers.Driving by the spin waves,the wall propagates towards the spin-wave source with an increasing velocity.Unlike the case without damping,the relation between the wall velocity and the spin-wave frequency depends on the position of wall.

Meshless Method with Domain Decomposition for Submerged Porous Breakwaters in Waves

Based on the improved version of the meshless singular boundary method(ISBM)in multi domain(MD),a numerical method is proposed in this paper to study the interaction of submerged permeable breakwaters and regular waves at normal incidence.To account for fluid flow inside the porous breakwaters,the conventional model of Sollitt and Cross for porous media is adopted.Both single and dual trapezoidal breakwaters are examined.The physical problem is formulated in the context of the linear potential wave theory.The domain decomposition method(DDM)is employed,in which the full computational domain is decomposed into separate domains,that is,the fluid domain and the domains of the breakwaters.Respectively,appropriate mixed type boundary and continuity conditions are applied for each subdomain and at the interfaces between domains.The solution is approximated in each subdomain by the ISBM.The discretized algebraic equations are combined,resulting in an overdetermined full system that is solved using a least-square solution procedure.The numerical results are presented in terms of the hydrodynamic quantities of reflection,transmission,and wave-energy dissipation.The relevance of the results of the present numerical procedure is first validated against data of previous studies,and then selected computations are discussed for various structural conditions.The proposed method is demonstrated to be highly accurate and computationally efficient.

Modeling of Surface Waves in a Fluid Saturated Poro-Elastic Medium under Initial Stress Using Time-Space Domain Higher Order Finite Difference Method

In this present context, mathematical modeling of the propagation of surface waves in a fluid saturated poro-elastic medium under the influence of initial stress has been considered using time dependent higher order finite difference method (FDM). We have proved that the accuracy of this finite-difference scheme is 2M when we use 2nd order time domain finite-difference and 2M-th order space domain finite-difference. It also has been shown that the dispersion curves of Love waves are less dispersed for higher order FDM than of lower order FDM. The effect of initial stress, porosity and anisotropy of the layer in the propagation of Love waves has been studied here. The numerical results have been shown graphically. As a particular case, the phase velocity in a non porous elastic solid layer derived in this paper is in perfect agreement with that of Liu et al. (2009).

Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain

In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.

Conformal Wavelet Finite-difference Time-domain Algorithm for Analysis of Milimeter Wave Attenuation on Coplanar Waveguide

Addressed is the calculation of millimeter wave attenuation on coplanar waveguide(CPW). A novel conformal wavelet finite-difference time-domain(CWFDTD) algorithm is proposed with emphasis on its application in calculation of millimeter wave attenuation on CPW, which is the combination of conformal algorithm dealing with the deformed cell with Wavelet-FDTD using multi-resolution analysis(MRA). Derived is the difference formulation for multi-resolution time domain(MRTD) based on Daubechies wavelets, and also given is the stability conditions for wavelet-FDTD algorithm. To validate its accuracy and efficiency, this novel method is applied to calculate the millimeter wave attenuation on lithium niobate CPW. Numerical results demonstrate that this new CWFDTD algorithm has the same accuracy with the conformal finite-difference time-domain(CFDTD) and conformal finite-difference time-domain based on alternating-direction implicit method(ADI-CFDTD), but saves computational time and computer memory.