In this paper, wave simulation with the finite difference method for the Helmholtz equation based on the domain decomposition method is investigated. The method solves the problem by iteratively solving subproblems de...
详细信息
In this paper, wave simulation with the finite difference method for the Helmholtz equation based on the domain decomposition method is investigated. The method solves the problem by iteratively solving subproblems defined on smaller subdomains. Two domain decomposition algorithms both for nonoverlapping and overlapping methods are described. More numerical computations including the benchmark Marmousi model show the effectiveness of the proposed algorithms. This method can be expected to be used in the full-waveform inversion in the future.
A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrim...
详细信息
A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrimination of the elastic wave equation are derived. The algorithm of no overlapping domain decomposition method is given. The method solves the elastic wave equation by iteratively solving sub problems defined on smaller sub domains. Numerical computations both for homogeneous and inhomogeneous media show the effectiveness of the proposed method. This method can be used in the full-waveform inversion.
In this paper, wave simulation with the finite difference method for the Helmholtz equation based on the domain decomposition method is investigated. The method solves the problem by iteratively solving subproblems de...
详细信息
In this paper, wave simulation with the finite difference method for the Helmholtz equation based on the domain decomposition method is investigated. The method solves the problem by iteratively solving subproblems defined on smaller subdomains. Two domain decomposition algorithms both for nonoverlapping and overlapping methods are described. More numerical computations including the benchmark Marmousi model show the effectiveness of the proposed algorithms. This method can be expected to be used in the full-waveform inversion in the future.
A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrim...
A new wave simulation technique for the elastic wave equation in the frequency domain based on a no overlapping domain decomposition algorithm is investigated. The boundary conditions and the finite difference discrimination of the elastic wave equation are derived. The algorithm of no overlapping domain decomposition method is given. The method solves the elastic wave equation by iteratively solving sub problems defined on smaller sub domains. Numerical computations both for homogeneous and inhomogeneous media show the effectiveness of the proposed method. This method can be used in the full-waveform inversion.
Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data...
Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.
Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data...
详细信息
Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.
This paper studied subspace properties of the Celis–Dennis–Tapia(CDT)subproblem that arises in some trust-region algorithms for equality constrained opti*** analysis is an extension of that presented by Wang and Yu...
详细信息
This paper studied subspace properties of the Celis–Dennis–Tapia(CDT)subproblem that arises in some trust-region algorithms for equality constrained opti*** analysis is an extension of that presented by Wang and Yuan(***.104:241–269,2006)for the standard trust-region *** suitable conditions,it is shown that the trial step obtained from the CDT subproblem is in the subspace spanned by all the gradient vectors of the objective function and of the constraints computed until the current *** on this observation,a subspace version of the Powell–Yuan trust-region algorithm is proposed for equality constrained optimization problems where the number of constraints is much lower than the number of variables. The convergence analysis is given and numerical results arealso reported.
Difference equations or discrete systems are mathematical models of various fields such as physics, chemistry, biology, and economics and have been subjects of extensive study of both pure mathematicians and applied m...
详细信息
Difference equations or discrete systems are mathematical models of various fields such as physics, chemistry, biology, and economics and have been subjects of extensive study of both pure mathematicians and applied mathematicians. Through its interaction with modern integrable systems, the theory of difference equations is enriched greatly and has been undergoing a rapid development. SIDE-10, the tenth of a series of biennial conferences devoted to Symmetries and Integrability of Difference Equations and related topics, was held during 10-16 June, 2012 at Ningbo, China. It was sponsored and supported by the National Natural science Foundation of China, Ningbo Association of science and Technology, Ningbo University, Academy of Mathematics and Systems science of Chinese Academy of sciences, China University of Mining and Technology （Beijing）, Tsinghua University, and Shanghai University. The conference attracted over 100 participants from more than a dozen of countries. During the conference, 44 contributed talks were arranged and the topics covered by the meeting include
Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlets and outlets at arbitrary orientations with respect to the lat...
详细信息
Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlets and outlets at arbitrary orientations with respect to the lattice, intermediate Reynolds (Re) number, and unsteady flow. Simple bounce-back is one of the most commonly used, simplest, and most computationally efficient boundary conditions, but many others have been proposed. We implement three other methods applicable to complex geometries [Guo, Zheng, and Shi, Phys. Fluids 14, 2007 (2002); Bouzidi, Firdaouss, and Lallemand, Phys. Fluids 13, 3452 (2001); Junk and Yang, Phys. Rev. E 72, 066701 (2005)] in our open-source application hemelb. We use these to simulate Poiseuille and Womersley flows in a cylindrical pipe with an arbitrary orientation at physiologically relevant Re number (1–300) and Womersley (4–12) numbers and steady flow in a curved pipe at relevant Dean number (100–200) and compare the accuracy to analytical solutions. We find that both the Bouzidi-Firdaouss-Lallemand (BFL) and Guo-Zheng-Shi (GZS) methods give second-order convergence in space while simple bounce-back degrades to first order. The BFL method appears to perform better than GZS in unsteady flows and is significantly less computationally expensive. The Junk-Yang method shows poor stability at larger Re number and so cannot be recommended here. The choice of collision operator (lattice Bhatnagar-Gross-Krook vs multiple relaxation time) and velocity set (D3Q15 vs D3Q19 vs D3Q27) does not significantly affect the accuracy in the problems studied.
暂无评论