Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elim...Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elimination cannot resolve the parallel computation of huge data. Therefore, we use the Gaussian elimination with static pivoting （GESP） method for sparse matrix decomposition and multi-source finite-difference modeling. The GESP method does not only improve the computational efficiency but also benefit the distributed parallel computation of matrix decomposition within a single frequency point. We test the proposed method using the classic Marmousi model. Both the single-frequency wave field and time domain seismic section show that the proposed method improves the simulation accuracy and computational efficiency and saves and makes full use of memory. This method can lay the basis for waveform inversion.展开更多
The effects of hurricane forward speed(V)and approach angle(θ)on storm surge are important and a systematic investigation covering possible and continuous ranges of these parameters has not been done before.Here we p...The effects of hurricane forward speed(V)and approach angle(θ)on storm surge are important and a systematic investigation covering possible and continuous ranges of these parameters has not been done before.Here we present such a study with a numerical experiment using the Finite Volume Community Ocean Model(FVCOM).The hurricane track is simplified as a straight line,such that V andθfully define the motion of the hurricane.The maximum surge is contributed by both free waves and a forced storm surge wave moving with the hurricane.Among the free waves,Kelvin-type waves can only propagate in the down-coast direction.Simulations show that those waves can only have a significant positive storm surge when the hurricane velocity has a down-coast component.The optimal values of V andθthat maximize the storm surge in an idealized semi-circular ocean basin are functions of the bathymetry.For a constant bathymetry,the maximum surge occurs when the hurricane approaches the coast from the normal direction when the free wave generation is minimal;for a stepped bathymetry,the maximum surge occurs at a certain acute approach angle which maximizes the duration of persistent wind forcing;a step-like bathymetry with a sloped shelf is similar to the stepped bathymetry,with the added possibility of landfall resonance when the free and forced waves are moving at about the same velocity.For other cases,the storm surge is smaller,given other parameters(hurricane size,maximum wind speed,etc.)unchanged.展开更多
基金supported by China State Key Science and Technology Project on Marine Carbonate Reservoir Characterization (No. 2008ZX05004-006)
文摘Frequency domain wave equation forward modeling is a problem of solving large scale linear sparse systems which is often subject to the limits of computational efficiency and memory storage. Conventional Gaussian elimination cannot resolve the parallel computation of huge data. Therefore, we use the Gaussian elimination with static pivoting （GESP） method for sparse matrix decomposition and multi-source finite-difference modeling. The GESP method does not only improve the computational efficiency but also benefit the distributed parallel computation of matrix decomposition within a single frequency point. We test the proposed method using the classic Marmousi model. Both the single-frequency wave field and time domain seismic section show that the proposed method improves the simulation accuracy and computational efficiency and saves and makes full use of memory. This method can lay the basis for waveform inversion.
基金The National Key R&D Project under contract No.2017YFC1404201
文摘The effects of hurricane forward speed(V)and approach angle(θ)on storm surge are important and a systematic investigation covering possible and continuous ranges of these parameters has not been done before.Here we present such a study with a numerical experiment using the Finite Volume Community Ocean Model(FVCOM).The hurricane track is simplified as a straight line,such that V andθfully define the motion of the hurricane.The maximum surge is contributed by both free waves and a forced storm surge wave moving with the hurricane.Among the free waves,Kelvin-type waves can only propagate in the down-coast direction.Simulations show that those waves can only have a significant positive storm surge when the hurricane velocity has a down-coast component.The optimal values of V andθthat maximize the storm surge in an idealized semi-circular ocean basin are functions of the bathymetry.For a constant bathymetry,the maximum surge occurs when the hurricane approaches the coast from the normal direction when the free wave generation is minimal;for a stepped bathymetry,the maximum surge occurs at a certain acute approach angle which maximizes the duration of persistent wind forcing;a step-like bathymetry with a sloped shelf is similar to the stepped bathymetry,with the added possibility of landfall resonance when the free and forced waves are moving at about the same velocity.For other cases,the storm surge is smaller,given other parameters(hurricane size,maximum wind speed,etc.)unchanged.
