Reduction of the wave propagation error of a sigma grid based numerical tank using a vertical spacing based on the constant truncation error

2021 ◽  
Vol 239 ◽  
pp. 109741
Author(s):  
Csaba Pakozdi ◽  
Weizhi Wang ◽  
Arun Kamath ◽  
Hans Bihs
2015 ◽  
Vol 8 (6) ◽  
pp. 2359-2369 ◽  
Author(s):  
R. A. Vincent ◽  
A Dudhia ◽  
L. J. Ventress

Abstract. This work presents a new iterative method for optimally selecting a vertical retrieval grid based on the location of the information while accounting for inter-level correlations. Sample atmospheres initially created to parametrise the Radiative Transfer Model for the Television Infrared Observation Satellite Operational Vertical Sounder (RTTOV) forward model are used to compare the presented iterative selection method with two other common approaches, which are using levels of equal vertical spacing and selecting levels based on the cumulative trace of the averaging kernel matrix (AKM). This new method is shown to outperform compared methods for simulated profile retrievals of temperature, H2O, O3, CH4, and CO with the Infrared Atmospheric Sounding Interferometer (IASI). However, the benefits of using the more complicated iterative approach compared to the simpler cumulative trace method are slight and may not justify the added effort for the cases studied, but may be useful in other scenarios where temperature and trace gases have strong vertical gradients with significant estimate sensitivity. Furthermore, comparing retrievals using a globally optimised static grid vs. a locally adapted one shows that a static grid performs nearly as well for retrievals of O3, CH4, and CO. However, developers of temperature and H2O retrieval schemes may at least consider using adaptive or location specific vertical retrieval grids.


2020 ◽  
Vol 110 (2) ◽  
pp. 754-762 ◽  
Author(s):  
Chuan Li ◽  
Jianxin Liu ◽  
Bo Chen ◽  
Ya Sun

ABSTRACT Seismic wavefield modeling based on the wave equation is widely used in understanding and predicting the dynamic and kinematic characteristics of seismic wave propagation through media. This article presents an optimal numerical solution for the seismic acoustic wave equation in a Hamiltonian system based on the third-order symplectic integrator method. The least absolute truncation error analysis method is used to determine the optimal coefficients. The analysis of the third-order symplectic integrator shows that the proposed scheme exhibits high stability and minimal truncation error. To illustrate the accuracy of the algorithm, we compare the numerical solutions generated by the proposed method with the theoretical analysis solution for 2D and 3D seismic wave propagation tests. The results show that the proposed method reduced the phase error to the eighth-order magnitude accuracy relative to the exact solution. These simulations also demonstrated that the proposed third-order symplectic method can minimize numerical dispersion and preserve the waveforms during the simulation. In addition, comparing different central frequencies of the source and grid spaces (90, 60, and 20 m) for simulation of seismic wave propagation in 2D and 3D models using symplectic and nearly analytic discretization methods, we deduce that the suitable grid spaces are roughly equivalent to between one-fourth and one-fifth of the wavelength, which can provide a good compromise between accuracy and computational cost.


Author(s):  
Elliott S. Wise ◽  
James L. B. Robertson ◽  
Ben T. Cox ◽  
Bradley E. Treeby

2015 ◽  
Vol 8 (3) ◽  
pp. 2591-2621
Author(s):  
R. A. Vincent ◽  
A. Dudhia ◽  
L. J. Ventress

Abstract. This work presents a new iterative method for optimally selecting a vertical retrieval grid based on the location of the information content while accounting for inter-level correlations. Sample atmospheres initially created to parametrise the RTTOV forward model are used to compare the presented iterative vertical selection method with two other common approaches, which are using levels of equal vertical spacing and selecting levels based on the cumulative trace of the averaging kernel matrix (AKM). This new method is shown to outperform compared methods for synthesized profile retrievals with IASI of temperature, H2O, O3, CH4, and CO. However, the benefits of using the more complicated iterative approach compared to the simpler method of referencing the cumulative trace of the AKM are slight and may not justify the added effort. Furthermore, comparing retrievals using a globally optimized static grid vs. an atmosphere specific one shows that a static grid is likely appropriate for retrievals of O3, CH4, and CO. However, developers of temperature and H2O retrieval schemes may at least consider using adaptive or location specific vertical retrieval grids.


Author(s):  
J. M. Galbraith ◽  
L. E. Murr ◽  
A. L. Stevens

Uniaxial compression tests and hydrostatic tests at pressures up to 27 kbars have been performed to determine operating slip systems in single crystal and polycrystal1ine beryllium. A recent study has been made of wave propagation in single crystal beryllium by shock loading to selectively activate various slip systems, and this has been followed by a study of wave propagation and spallation in textured, polycrystal1ine beryllium. An alteration in the X-ray diffraction pattern has been noted after shock loading, but this alteration has not yet been correlated with any structural change occurring during shock loading of polycrystal1ine beryllium.This study is being conducted in an effort to characterize the effects of shock loading on textured, polycrystal1ine beryllium. Samples were fabricated from a billet of Kawecki-Berylco hot pressed HP-10 beryllium.


Author(s):  
J. M. Pankratz

It is often desirable in transmission electron microscopy to know the vertical spacing of points of interest within a specimen. However, in order to measure a stereo effect, one must have two pictures of the same area taken from different angles, and one must have also a formula for converting measured differences between corresponding points (parallax) into a height differential.Assume (a) that the impinging beam of electrons can be considered as a plane wave and (b) that the magnification is the same at the top and bottom of the specimen. The first assumption is good when the illuminating system is overfocused. The second assumption (the so-called “perspective error”) is good when the focal length is large (3 x 107Å) in relation to foil thickness (∼103 Å).


Sign in / Sign up

Export Citation Format

Share Document