We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to sol...We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to solve the cutoff botmdary problem as well as to consider the requirement of using the same subdivision grid in joint multifrequency inversion. The proposed method introduces the attenuation factor, and by adjusting it, acoustic waves are sufficiently attenuated in the attenuation layer to minimize the cutoff boundary effect. Based on the law of exponential attenuation, expressions for computing the attenuation factor and the thickness of attenuation layers are derived for different frequencies. In multifrequency-domain FWI, the conjugate gradient method is used to solve equations in the Gauss-Newton algorithm and thus minimize the computation cost in calculating the Hessian matrix. In addition, the effect of initial model selection and frequency combination on FWI is analyzed. Examples using numerical simulations and FWI calculations are used to verify the efficiency of the proposed method.展开更多
In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object a...In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.展开更多
We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high...We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.展开更多
With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) meth...With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.展开更多
An idealized numerical wave flume has been established by finite element method on the bases of Navier Stokes equations through prescribing the appropriate boundary conditions for the open boundary,incident boundary,...An idealized numerical wave flume has been established by finite element method on the bases of Navier Stokes equations through prescribing the appropriate boundary conditions for the open boundary,incident boundary,free surface and solid boundary in this paper.The characteristics of waves propagating over a step have been investigated by this numerical model.The breaker wave height is determined depending on the kinetic criterion.The numerical model is verified by laboratory experiments,and the empirical formula for the damping of wave height due to breaking is also given by experiments.展开更多
Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI me...Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI media). This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer (PML) absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
Spans occur when a pipeline is laid on a rough undulating seabed or when upheaval buckling occurs due to constrained thermal expansion. This not only results in static and dynamic loads on the flowline at span section...Spans occur when a pipeline is laid on a rough undulating seabed or when upheaval buckling occurs due to constrained thermal expansion. This not only results in static and dynamic loads on the flowline at span sections,but also generates vortex induced vibration (VIV),which can lead to fatigue issues. The phenomenon,if not predicted and controlled properly,will negatively affect pipeline integrity,leading to expensive remediation and intervention work. Span analysis can be complicated by:long span lengths,a large number of spans caused by a rough seabed,and multi-span interactions. In addition,the complexity can be more onerous and challenging when soil uncertainty,concrete degradation and unknown residual lay tension are considered in the analysis. This paper describes the latest developments and a'state-of-the-art' finite element analysis program that has been developed to simulate the span response of a flowline under complex boundary and loading conditions. Both VIV and direct wave loading are captured in the analysis and the results are sequentially used for the ultimate limit state (ULS) check and fatigue life calculation.展开更多
Donnell’s thin shell theory and basic equations based on the wave propagation method discussed in detail here, is used to investigate the natural frequencies of thin finite length circular cylindrical shells under va...Donnell’s thin shell theory and basic equations based on the wave propagation method discussed in detail here, is used to investigate the natural frequencies of thin finite length circular cylindrical shells under various boundary conditions. Mode shapes are drawn to explain the circumferential mode number n and axial mode number m, and the natural frequencies are cal-culated numerically and compared with those of FEM (finite element method) to confirm the reliability of the analytical solution. The effects of relevant parameters on natural frequencies are discussed thoroughly. It is shown that for long thin shells the method is simple, accurate and effective.展开更多
Problems of fluid structure interactions are governed by a set of fundamental parameters. This work aims at showing through simple examples the changes in natural vibration frequencies and mode shapes for wall-cavity ...Problems of fluid structure interactions are governed by a set of fundamental parameters. This work aims at showing through simple examples the changes in natural vibration frequencies and mode shapes for wall-cavity systems when the structural rigidity is modified. Numerical results are constructed using ANSYS software with triangular finite elements for both the fluid (2D acoustic elements) and the solid (plane stress) domains. These former results are compared to proposed analytical expressions, showing an alternative benchmark tool for the analyst. Very rigid wall structures imply in frequencies and mode shapes almost identical to those achieved for an acoustic cavity with Neumann boundary condition at the interface. In this case, the wall behaves as rigid and fluid-structure system mode shapes are similar to those achieved for the uncoupled reservoir case.展开更多
In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experiment...In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations (PDEs) of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.展开更多
A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are p...A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are proved. Numerical examples are presented to verify the theoretical results.展开更多
The authors introduce a new Large Eddy Simulation model in a channel, based on the projection on finite element spaces as filtering operation in its variationM form, for a given triangulation (Th)h〉0. The eddy visc...The authors introduce a new Large Eddy Simulation model in a channel, based on the projection on finite element spaces as filtering operation in its variationM form, for a given triangulation (Th)h〉0. The eddy viscosity is expressed in terms of the friction velocity in the boundary layer due to the wall, and is of a standard sub grid-model form outside the boundary layer. The mixing length scale is locally equal to the grid size. The computational domain is the channel without the linear sub-layer of the boundary layer. The no-slip boundary condition (or BC for short) is replaced by a Navier (BC) at the computational wall. Considering the steady state case, the authors show that the variational finite element model they have introduced, has a solution (Vh,Ph)h〉O that converges to a solution of the steady state Navier-Stokes equation with Navier BC.展开更多
This paper aims to look into the determination of effective area-average concentration and dispersion coefficient associated with unsteady flow through a small-diameter tube where a solute undergoes first-order chemic...This paper aims to look into the determination of effective area-average concentration and dispersion coefficient associated with unsteady flow through a small-diameter tube where a solute undergoes first-order chemical reaction both within the fluid and at the boundary. The reaction consists of a reversible component due to phase exchange between the flowing fluid and the wall layer, and an irreversible component due to absorption into the wall. To understand the dispersion, the governing equations along with the reactive boundary conditions are solved numerically using the Finite Difference Method. The resultant equation shows how the dispersion coefficient is influenced by the first-order chemical reaction. The effects of various dimensionless parameters e.g. Da (the Damkohler number), a (phase partitioning number) and F (dimensionless absorption number) on dispersion are discussed. One of the results exposes that the dispersion coefficient may approach its steady-state limit in a short time at a high value of Damkohler number (say Da 〉 10) and a small but nonzero value of absorption rate (say P 〈0.5).展开更多
基金financially supported by the National High Technology Research and Development Program of China(No.2012AA09A20105)the National Science Foundation Network(No.41574127)
文摘We studied finite-element-method-based two-dimensional frequency-domain acoustic FWI under rugged topography conditions. The exponential attenuation boundary condition suitable for rugged topography is proposed to solve the cutoff botmdary problem as well as to consider the requirement of using the same subdivision grid in joint multifrequency inversion. The proposed method introduces the attenuation factor, and by adjusting it, acoustic waves are sufficiently attenuated in the attenuation layer to minimize the cutoff boundary effect. Based on the law of exponential attenuation, expressions for computing the attenuation factor and the thickness of attenuation layers are derived for different frequencies. In multifrequency-domain FWI, the conjugate gradient method is used to solve equations in the Gauss-Newton algorithm and thus minimize the computation cost in calculating the Hessian matrix. In addition, the effect of initial model selection and frequency combination on FWI is analyzed. Examples using numerical simulations and FWI calculations are used to verify the efficiency of the proposed method.
文摘In order to analyze the electrostatic field concerned with electrostatic proximity fuze problem using the available finite analysis software package, the technology to model the problem with a scale reduction object and boundary was presented. The boundary is determined by the maximum distance the sensor can detect. The object model is obtained by multiplying the terms in Poisson's equation with a scale reduction factor and the real value can be reconstructed with the same reverse process after software calculation. Using the finite element analysis program, the simulation value is close to the theoretical value with a little error. The boundary determination and scale reduction method is suitable to modeling the irregular electrostatic field around air targets, such as airplane, missile and so on, which is based on commonly used personal computer (PC). The technology reduces the calculation and storage cost greatly.
基金supported by the National Nature Science Foundation of China(Grant No.U1262208)the Important National Science & Technology Specific Projects(Grant No.2011ZX05019-008)
文摘We apply the newly proposed double absorbing boundary condition(DABC)(Hagstrom et al., 2014) to solve the boundary reflection problem in seismic finite-difference(FD) modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer(PML), the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.
基金The National Natural Science Foundation of China(No.60702027)the Free Research Fund of the National Mobile Communications Research Laboratory of Southeast University (No.2008B07)the National Basic Research Program of China(973 Program)(No.2007CB310603)
文摘With the linear interpolation method, an improved absorbing boundary condition(ABC)is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain (ADI-FDTD) method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.
文摘An idealized numerical wave flume has been established by finite element method on the bases of Navier Stokes equations through prescribing the appropriate boundary conditions for the open boundary,incident boundary,free surface and solid boundary in this paper.The characteristics of waves propagating over a step have been investigated by this numerical model.The breaker wave height is determined depending on the kinetic criterion.The numerical model is verified by laboratory experiments,and the empirical formula for the damping of wave height due to breaking is also given by experiments.
基金supported by the National Natural Science Foundation of China(No.41674118)the national science and technology major project(No.2016ZX05027-002)
文摘Based on the pure quasi-P wave equation in transverse isotropic media with a vertical symmetry axis (VTI media), a quasi-P wave equation is obtained in transverse isotropic media with a tilted symmetry axis (TTI media). This is achieved using projection transformation, which rotates the direction vector in the coordinate system of observation toward the direction vector for the coordinate system in which the z-component is parallel to the symmetry axis of the TTI media. The equation has a simple form, is easily calculated, is not influenced by the pseudo-shear wave, and can be calculated reliably when δ is greater than ε. The finite difference method is used to solve the equation. In addition, a perfectly matched layer (PML) absorbing boundary condition is obtained for the equation. Theoretical analysis and numerical simulation results with forward modeling prove that the equation can accurately simulate a quasi-P wave in TTI medium.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods (IFDM) for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition (ABC) to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
文摘Spans occur when a pipeline is laid on a rough undulating seabed or when upheaval buckling occurs due to constrained thermal expansion. This not only results in static and dynamic loads on the flowline at span sections,but also generates vortex induced vibration (VIV),which can lead to fatigue issues. The phenomenon,if not predicted and controlled properly,will negatively affect pipeline integrity,leading to expensive remediation and intervention work. Span analysis can be complicated by:long span lengths,a large number of spans caused by a rough seabed,and multi-span interactions. In addition,the complexity can be more onerous and challenging when soil uncertainty,concrete degradation and unknown residual lay tension are considered in the analysis. This paper describes the latest developments and a'state-of-the-art' finite element analysis program that has been developed to simulate the span response of a flowline under complex boundary and loading conditions. Both VIV and direct wave loading are captured in the analysis and the results are sequentially used for the ultimate limit state (ULS) check and fatigue life calculation.
文摘Donnell’s thin shell theory and basic equations based on the wave propagation method discussed in detail here, is used to investigate the natural frequencies of thin finite length circular cylindrical shells under various boundary conditions. Mode shapes are drawn to explain the circumferential mode number n and axial mode number m, and the natural frequencies are cal-culated numerically and compared with those of FEM (finite element method) to confirm the reliability of the analytical solution. The effects of relevant parameters on natural frequencies are discussed thoroughly. It is shown that for long thin shells the method is simple, accurate and effective.
文摘Problems of fluid structure interactions are governed by a set of fundamental parameters. This work aims at showing through simple examples the changes in natural vibration frequencies and mode shapes for wall-cavity systems when the structural rigidity is modified. Numerical results are constructed using ANSYS software with triangular finite elements for both the fluid (2D acoustic elements) and the solid (plane stress) domains. These former results are compared to proposed analytical expressions, showing an alternative benchmark tool for the analyst. Very rigid wall structures imply in frequencies and mode shapes almost identical to those achieved for an acoustic cavity with Neumann boundary condition at the interface. In this case, the wall behaves as rigid and fluid-structure system mode shapes are similar to those achieved for the uncoupled reservoir case.
文摘In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations (PDEs) of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.
基金supported by National Science Foundation of USA(Grant No.DMS-1418934)the Sea Poly Project of Beijing Overseas Talents,National Natural Science Foundation of China(Grant Nos.11625101,91430213,11421101,11771338,11671304 and 11401026)+1 种基金Zhejiang Provincial Natural Science Foundation of China Projects(Grant Nos.LY17A010010,LY15A010015 and LY15A010016)Wenzhou Science and Technology Plan Project(Grant No.G20160019)
文摘A posteriori error estimators for the symmetric mixed finite element methods for linear elasticity problems with Dirichlet and mixed boundary conditions are proposed. Reliability and efficiency of the estimators are proved. Numerical examples are presented to verify the theoretical results.
基金Project supported by the Spanish Government and European Union FEDER Grant(No.MTM200907719)
文摘The authors introduce a new Large Eddy Simulation model in a channel, based on the projection on finite element spaces as filtering operation in its variationM form, for a given triangulation (Th)h〉0. The eddy viscosity is expressed in terms of the friction velocity in the boundary layer due to the wall, and is of a standard sub grid-model form outside the boundary layer. The mixing length scale is locally equal to the grid size. The computational domain is the channel without the linear sub-layer of the boundary layer. The no-slip boundary condition (or BC for short) is replaced by a Navier (BC) at the computational wall. Considering the steady state case, the authors show that the variational finite element model they have introduced, has a solution (Vh,Ph)h〉O that converges to a solution of the steady state Navier-Stokes equation with Navier BC.
文摘This paper aims to look into the determination of effective area-average concentration and dispersion coefficient associated with unsteady flow through a small-diameter tube where a solute undergoes first-order chemical reaction both within the fluid and at the boundary. The reaction consists of a reversible component due to phase exchange between the flowing fluid and the wall layer, and an irreversible component due to absorption into the wall. To understand the dispersion, the governing equations along with the reactive boundary conditions are solved numerically using the Finite Difference Method. The resultant equation shows how the dispersion coefficient is influenced by the first-order chemical reaction. The effects of various dimensionless parameters e.g. Da (the Damkohler number), a (phase partitioning number) and F (dimensionless absorption number) on dispersion are discussed. One of the results exposes that the dispersion coefficient may approach its steady-state limit in a short time at a high value of Damkohler number (say Da 〉 10) and a small but nonzero value of absorption rate (say P 〈0.5).