[Objective] The research aimed to analyze the moist potential vorticity (MPV) and conditional symmetric instability (CSI) of a rainstorm in Guangxi in the prior flood season. [Method] Based on the conventional observa...[Objective] The research aimed to analyze the moist potential vorticity (MPV) and conditional symmetric instability (CSI) of a rainstorm in Guangxi in the prior flood season. [Method] Based on the conventional observation data, precipitation data from the automatic station and 1°×1° six-hourly reanalyzed data from NCEP, MPV and CSI in a rainstorm process which occurred in the prior flood season in Guangxi were analyzed. The characteristics of MPV component and role of CSI before and after the rainstorm occurrence were discussed. [Result] The dense belt of isoline in the southeast of MPV1 positive-value area (MPV2 negative-value area) at 925 hPa, the underneath of sloping MPV1 positive-value column (MPV2 negative-value column) on the profile map and the underneath of dense belt of θe isoline had the good corresponding relationships with the falling zone of strong precipitation. After the system developed strongly, the enhancements of gradient and intensity of MPV1 and MPV2 predicted the rainstorm increase. In the beginning period of strong precipitation occurrence, the convective instability at 850-650 hPa co-existed with CSI. The convective instability and vertical movement played the main role. In the middle and latter periods of strong precipitation, the high-altitude weak cold air further invaded southward. It triggered the release of unstable energy, and compelled that the original sloping θe dense belt became steep. It was favorable for the development of sloping vorticity, and the atmosphere gradually turned into the neutral state of convection. CSI and ramp motion played the major role. Then, the rainstorm continued. [Conclusion] The research provided the theory basis for the application of MPV and CSI into the local routine business.展开更多
In the light of the theory on moist potential vorticity (MPV) investigation was undertaken of the 700 hPa vertical (horizontal) component MP1 (MPV2) for the heavy rain event occurring in July 5–6, 1991. Results show ...In the light of the theory on moist potential vorticity (MPV) investigation was undertaken of the 700 hPa vertical (horizontal) component MP1 (MPV2) for the heavy rain event occurring in July 5–6, 1991. Results show that the distribution features of the two components were closely related to the development of a mesoscale cyclone as a rainstorm-causing weather system in the lower troposphere in such a way that the ambient atmosphere of which MPV1 > 0 and MPV2 < 0 with |MPV1| ≥ |MPV2| favored the genesis of conditional symmetric instability (CSI) and that, as indicated by calculations, a CSI sector was really existent in the lower troposphere during the heavy rain happening and contributed greatly to its development.展开更多
Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and...Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.展开更多
Rayleigh-Taylor (R-T) instability is known as the fundamental mechanism of equatorial plasma bubbles (EPBs). However, the sufficient conditions of R-T instability and stability have not yet been derived. In the pr...Rayleigh-Taylor (R-T) instability is known as the fundamental mechanism of equatorial plasma bubbles (EPBs). However, the sufficient conditions of R-T instability and stability have not yet been derived. In the present paper, the sufficient conditions of R-T stability and instability are preliminarily^derived. Linear equations for small perturbation are first obtained from the electron/ion continuity equations, momentum equations, and the current continuity equation in the equatorial ionosphere. The linear equations can be casted as an eigenvalue equation using a normal mode method. The eigenvalue equation is a variable coefficient linear equation that can be solved using a variational approach. With this approach, the sufficient conditions can be obtained as follows: if the minimum systematic eigenvalue is greater than one, the ionosphere is R-T unstable; while if the maximum systematic eigenvalue is less than one, the ionosphere is R-T stable. An approximate numerical method for obtaining the systematic eigenvalues is introduced, and the R-T stable/unstable areas are calculated. Numerical experiments axe designed to validate the sufficient conditions. The results agree with the derived suf- ficient conditions.展开更多
In this paper the analytical expression of free energy expressed by small parameter r of a weakly interacting Fermi gas trapped in weak magnetic field is derived by using 'the maximum approximation' method and the e...In this paper the analytical expression of free energy expressed by small parameter r of a weakly interacting Fermi gas trapped in weak magnetic field is derived by using 'the maximum approximation' method and the ensemble theory. Based on the derived expression, the exact instability conditions of a weakly interacting Fermi gas trapped in weak magnetic field at both high and low temperatures are given. From the instability conditions we get the following two results. (1) At the whole low-temperature extent, whether the interactions are repulsive or attractive with (αn-k 4εF/3) (n and εF denote the particle-number density and the Fermi energy respectively, α = 4πah^2/m, and a is s-wave scattering length) positive, there is a lower-limit magnetic field of instability; in addition, there is an upper-limit magnetic field for the system of attractive interactions with (αn + 4εF/3) negative. (2) At the whole high-temperature extent, the system with repulsive interactions is always stable, but for the system with attractive interactions, the greater the scattering length of attractive interactions |α| is, the stronger the magnetic field is and the larger the particle-number density is, the bigger the possibility of instability in the system will be.展开更多
Observation has clearly shown that natural space plasmas generally possess a pronounced non-Maxwellian high-energy tail distribution that can be well modeled by a kappa distribution. In this study we investigate the p...Observation has clearly shown that natural space plasmas generally possess a pronounced non-Maxwellian high-energy tail distribution that can be well modeled by a kappa distribution. In this study we investigate the proton cyclotron wave instability driven by the temperature anisotropy (T⊥/TH 〉1) of suprathermal protons modeled with a typical kappa distribution in the magnetosheath. It is found that as in the case for a regular bi-Maxwellian, the supratherreal proton temperature anisotropy is subject to the threshold condition of this proton cyclotron instability and the instability threshold condition satisfies a general form T⊥/T|| - 1 = S/β||^α, with a very narrow range of the fitting parameters: 0.40 ≤ α ≤ 0.45, and a relatively sensitive variation 0.27 ≤ S ≤ 0.65, over 0.01 ≤β|| 〈 10. Furthermore, the difference in threshold conditions between the kappa distribution and the bi-Maxwellian distribution is found to be small for a relatively strong growth but becomes relatively obvious for a weak wave growth. The results may provide a deeper insight into the physics of this instability threshold for the proton cyclotron waves.展开更多
In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensio...In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.展开更多
The computational cost of numerical methods in microscopic-scales such as molecular dynamics(MD) is a deterrent factor that limits simulations with a large number of particles. Hence, it is desirable to decrease the c...The computational cost of numerical methods in microscopic-scales such as molecular dynamics(MD) is a deterrent factor that limits simulations with a large number of particles. Hence, it is desirable to decrease the computational cost and run time of simulations, especially for problems with a symmetrical domain. However, in microscopic-scales, implementation of symmetric boundary conditions is not straightforward. Previously, the present authors have successfully used a symmetry boundary condition to solve molecular flows in constant-area channels. The results obtained with this approach agree well with the benchmark cases. Therefore, it has provided us with a sound ground to further explore feasibility of applying symmetric solutions of micro-fluid flows in other geometries such as variable-area ducts. Molecular flows are solved for the whole domain with and without the symmetric boundary condition. Good agreement has been reached between the results of the symmetric solution and the whole domain solution. To investigate robustness of the proposed method, simulations are conducted for different values of affecting parameters including an external force, a flow density, and a domain length. The results indicate that the symmetric solution is also applicable to variable-area ducts such as micro-nozzles.展开更多
Shock induced symmetric compression has been studied in a spherical target. The shock induced interfacial radius will shrink and would reach a minimum point during implosion situation. However, after implosion the pla...Shock induced symmetric compression has been studied in a spherical target. The shock induced interfacial radius will shrink and would reach a minimum point during implosion situation. However, after implosion the plasma tries to expand in blow off/explosion situation and as a result the interfacial radius will increase. Effects of plasma parameters like density and temperature have been studied numerically. It is seen that the density increases many times due to the mass conservation in imploding situation of a compressible shell like ICF. However, temperature will change rapidly due to change of inner density and so would be the pressure of compressible fluid following adiabatic law. Our analytical results agree qualitatively with those of simulation results in spherical geometry and also experimental observations conducted in cylindrical container.展开更多
The asymptotic behaviors for estimators of the drift parameters in the Ornstein-Uhlenbeck process driven by small symmetricα-stable motion are studied in this paper.Based on the discrete observations,the conditional ...The asymptotic behaviors for estimators of the drift parameters in the Ornstein-Uhlenbeck process driven by small symmetricα-stable motion are studied in this paper.Based on the discrete observations,the conditional least squares estimators(CLSEs)of all the parameters involved in the Ornstein–Uhlenbeck process are proposed.We establish the consistency and the asymptotic distributions of our estimators asεgoes to 0 and n goes to∞simultaneously.展开更多
Imposing appropriate numerical boundary conditions at the symmetrical center r=0 is vital when computing compressible fluids with radial symmetry.Extrapolation and other traditional techniques are often employed,but s...Imposing appropriate numerical boundary conditions at the symmetrical center r=0 is vital when computing compressible fluids with radial symmetry.Extrapolation and other traditional techniques are often employed,but spurious numerical oscillations or wall-heating phenomena can occur.In this paper,we emphasize that because of the conservation property,the updating formula of the boundary cell average can coincide with the one for interior cell averages.To achieve second-order accuracy both in time and space,we associate obtaining the inner boundary value at r=0 with the resolution of the corresponding one-sided generalized Riemann problem(GRP).Acoustic approximation is applied in this process.It creates conditions to avoid the singularity of type 1/r and aids in obtaining the value of the singular quantity using L’Hospital’s rule.Several challenging scenarios are tested to demonstrate the effectiveness and robustness of our approach.展开更多
Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two...Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.展开更多
NCEP-NCAR reanalysis data were used to analyze the characteristics and evolution mechanism of convective and symmetric instability before and during a heavy rainfall event that occurred in Beijing on 21 July 2012.Appr...NCEP-NCAR reanalysis data were used to analyze the characteristics and evolution mechanism of convective and symmetric instability before and during a heavy rainfall event that occurred in Beijing on 21 July 2012.Approximately twelve hours before the rainstorm,the atmosphere was mainly dominated by convective instability in the lower level of 900-800 hPa.The strong southwesterly low-level jet conveyed the moist and warm airflow continuously to the area of torrential rain,maintaining and enhancing the unstable energy.When the precipitation occurred,unstable energy was released and the convective instability weakened.Meanwhile,due to the baroclinicity enhancement in the atmosphere,the symmetric instability strengthened,maintaining and promoting the subsequent torrential rain.Deriving the convective instability tendency equation demonstrated that the barotropic component of potential divergence and the advection term played a major role in enhancing the convective instability before the rainstorm.Analysis of the tendency equation of moist potential vorticity showed that the coupled term of vertical vorticity and the baroclinic component of potential divergence was the primary factor influencing the development of symmetric instability during the precipitation.Comparing the effects of these factors on convective instability and symmetric instability showed some correlation.展开更多
From the formulas of the conjugate gradient, a similarity between a symmetric positive definite (SPD) matrix A and a tridiagonal matrix B is obtained. The elements of the matrix B are determined by the parameters of t...From the formulas of the conjugate gradient, a similarity between a symmetric positive definite (SPD) matrix A and a tridiagonal matrix B is obtained. The elements of the matrix B are determined by the parameters of the conjugate gradient. The computation of eigenvalues of A is then reduced to the case of the tridiagonal matrix B. The approximation of extreme eigenvalues of A can be obtained as a 'by-product' in the computation of the conjugate gradient if a computational cost of O(s) arithmetic operations is added, where s is the number of iterations This computational cost is negligible compared with the conjugate gradient. If the matrix A is not SPD, the approximation of the condition number of A can be obtained from the computation of the conjugate gradient on AT A. Numerical results show that this is a convenient and highly efficient method for computing extreme eigenvalues and the condition number of nonsingular matrices.展开更多
文摘[Objective] The research aimed to analyze the moist potential vorticity (MPV) and conditional symmetric instability (CSI) of a rainstorm in Guangxi in the prior flood season. [Method] Based on the conventional observation data, precipitation data from the automatic station and 1°×1° six-hourly reanalyzed data from NCEP, MPV and CSI in a rainstorm process which occurred in the prior flood season in Guangxi were analyzed. The characteristics of MPV component and role of CSI before and after the rainstorm occurrence were discussed. [Result] The dense belt of isoline in the southeast of MPV1 positive-value area (MPV2 negative-value area) at 925 hPa, the underneath of sloping MPV1 positive-value column (MPV2 negative-value column) on the profile map and the underneath of dense belt of θe isoline had the good corresponding relationships with the falling zone of strong precipitation. After the system developed strongly, the enhancements of gradient and intensity of MPV1 and MPV2 predicted the rainstorm increase. In the beginning period of strong precipitation occurrence, the convective instability at 850-650 hPa co-existed with CSI. The convective instability and vertical movement played the main role. In the middle and latter periods of strong precipitation, the high-altitude weak cold air further invaded southward. It triggered the release of unstable energy, and compelled that the original sloping θe dense belt became steep. It was favorable for the development of sloping vorticity, and the atmosphere gradually turned into the neutral state of convection. CSI and ramp motion played the major role. Then, the rainstorm continued. [Conclusion] The research provided the theory basis for the application of MPV and CSI into the local routine business.
文摘In the light of the theory on moist potential vorticity (MPV) investigation was undertaken of the 700 hPa vertical (horizontal) component MP1 (MPV2) for the heavy rain event occurring in July 5–6, 1991. Results show that the distribution features of the two components were closely related to the development of a mesoscale cyclone as a rainstorm-causing weather system in the lower troposphere in such a way that the ambient atmosphere of which MPV1 > 0 and MPV2 < 0 with |MPV1| ≥ |MPV2| favored the genesis of conditional symmetric instability (CSI) and that, as indicated by calculations, a CSI sector was really existent in the lower troposphere during the heavy rain happening and contributed greatly to its development.
基金funded by National Natural Science Foundation of China (No. 12065015)the Hongliu First-level Discipline Construction Project of Lanzhou University of Technology。
文摘Terahertz(THz) radiation can be generated due to the instability of THz plasma waves in field-effect transistors(FETs). In this work, we discuss the instability of THz plasma waves in the channel of FETs with spin and quantum effects under non-ideal boundary conditions. We obtain a linear dispersion relation by using the hydrodynamic equation, Maxwell equation and spin equation. The influence of source capacitance, drain capacitance, spin effects, quantum effects and channel width on the instability of THz plasma waves under the non-ideal boundary conditions is investigated in great detail. The results of numerical simulation show that the THz plasma wave is unstable when the drain capacitance is smaller than the source capacitance;the oscillation frequency with asymmetric boundary conditions is smaller than that under non-ideal boundary conditions;the instability gain of THz plasma waves becomes lower under non-ideal boundary conditions. This finding provides a new idea for finding efficient THz radiation sources and opens up a new mechanism for the development of THz technology.
基金Project supported by the National Natural Science Foundation of China(Nos.41575026 and 41175025)
文摘Rayleigh-Taylor (R-T) instability is known as the fundamental mechanism of equatorial plasma bubbles (EPBs). However, the sufficient conditions of R-T instability and stability have not yet been derived. In the present paper, the sufficient conditions of R-T stability and instability are preliminarily^derived. Linear equations for small perturbation are first obtained from the electron/ion continuity equations, momentum equations, and the current continuity equation in the equatorial ionosphere. The linear equations can be casted as an eigenvalue equation using a normal mode method. The eigenvalue equation is a variable coefficient linear equation that can be solved using a variational approach. With this approach, the sufficient conditions can be obtained as follows: if the minimum systematic eigenvalue is greater than one, the ionosphere is R-T unstable; while if the maximum systematic eigenvalue is less than one, the ionosphere is R-T stable. An approximate numerical method for obtaining the systematic eigenvalues is introduced, and the R-T stable/unstable areas are calculated. Numerical experiments axe designed to validate the sufficient conditions. The results agree with the derived suf- ficient conditions.
文摘In this paper the analytical expression of free energy expressed by small parameter r of a weakly interacting Fermi gas trapped in weak magnetic field is derived by using 'the maximum approximation' method and the ensemble theory. Based on the derived expression, the exact instability conditions of a weakly interacting Fermi gas trapped in weak magnetic field at both high and low temperatures are given. From the instability conditions we get the following two results. (1) At the whole low-temperature extent, whether the interactions are repulsive or attractive with (αn-k 4εF/3) (n and εF denote the particle-number density and the Fermi energy respectively, α = 4πah^2/m, and a is s-wave scattering length) positive, there is a lower-limit magnetic field of instability; in addition, there is an upper-limit magnetic field for the system of attractive interactions with (αn + 4εF/3) negative. (2) At the whole high-temperature extent, the system with repulsive interactions is always stable, but for the system with attractive interactions, the greater the scattering length of attractive interactions |α| is, the stronger the magnetic field is and the larger the particle-number density is, the bigger the possibility of instability in the system will be.
基金supported by the National Natural Science Foundation of China (Nos. 40474064, 40404012) the Scientific Research Foundation for R0CS, SEM+1 种基金 the Scientific Research Fund of Hunan Provincial Science and Technology Department grant 05FJ3045 the Visiting Scholar Foundation of State Key Laboratory of Space Weather, Chinese Academy of Sciences
文摘Observation has clearly shown that natural space plasmas generally possess a pronounced non-Maxwellian high-energy tail distribution that can be well modeled by a kappa distribution. In this study we investigate the proton cyclotron wave instability driven by the temperature anisotropy (T⊥/TH 〉1) of suprathermal protons modeled with a typical kappa distribution in the magnetosheath. It is found that as in the case for a regular bi-Maxwellian, the supratherreal proton temperature anisotropy is subject to the threshold condition of this proton cyclotron instability and the instability threshold condition satisfies a general form T⊥/T|| - 1 = S/β||^α, with a very narrow range of the fitting parameters: 0.40 ≤ α ≤ 0.45, and a relatively sensitive variation 0.27 ≤ S ≤ 0.65, over 0.01 ≤β|| 〈 10. Furthermore, the difference in threshold conditions between the kappa distribution and the bi-Maxwellian distribution is found to be small for a relatively strong growth but becomes relatively obvious for a weak wave growth. The results may provide a deeper insight into the physics of this instability threshold for the proton cyclotron waves.
基金Basic Scientific Research-related Project from Institute of Engineering Mechanics (01180001 and 2007C01)
文摘In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.
文摘The computational cost of numerical methods in microscopic-scales such as molecular dynamics(MD) is a deterrent factor that limits simulations with a large number of particles. Hence, it is desirable to decrease the computational cost and run time of simulations, especially for problems with a symmetrical domain. However, in microscopic-scales, implementation of symmetric boundary conditions is not straightforward. Previously, the present authors have successfully used a symmetry boundary condition to solve molecular flows in constant-area channels. The results obtained with this approach agree well with the benchmark cases. Therefore, it has provided us with a sound ground to further explore feasibility of applying symmetric solutions of micro-fluid flows in other geometries such as variable-area ducts. Molecular flows are solved for the whole domain with and without the symmetric boundary condition. Good agreement has been reached between the results of the symmetric solution and the whole domain solution. To investigate robustness of the proposed method, simulations are conducted for different values of affecting parameters including an external force, a flow density, and a domain length. The results indicate that the symmetric solution is also applicable to variable-area ducts such as micro-nozzles.
文摘Shock induced symmetric compression has been studied in a spherical target. The shock induced interfacial radius will shrink and would reach a minimum point during implosion situation. However, after implosion the plasma tries to expand in blow off/explosion situation and as a result the interfacial radius will increase. Effects of plasma parameters like density and temperature have been studied numerically. It is seen that the density increases many times due to the mass conservation in imploding situation of a compressible shell like ICF. However, temperature will change rapidly due to change of inner density and so would be the pressure of compressible fluid following adiabatic law. Our analytical results agree qualitatively with those of simulation results in spherical geometry and also experimental observations conducted in cylindrical container.
基金Key Natural Science Foundation of Anhui Education Commission,China(No.KJ2017A568)Natural Science Foundation of Anhui Province,China(No.1808085MA02)+1 种基金Quality Engineering Project of Anhui Province,China(No.2019jyxm0476)Quality Engineering Project of Bengbu University,China(No.2018JYXML8)。
文摘The asymptotic behaviors for estimators of the drift parameters in the Ornstein-Uhlenbeck process driven by small symmetricα-stable motion are studied in this paper.Based on the discrete observations,the conditional least squares estimators(CLSEs)of all the parameters involved in the Ornstein–Uhlenbeck process are proposed.We establish the consistency and the asymptotic distributions of our estimators asεgoes to 0 and n goes to∞simultaneously.
基金This work was partially supported by Science Challenge project TZ2016002NSFC with Nos.11771055,11671050,11871113,11871114,12026607,121710493D numerical simulation platform TB14-1 of the China Academy of Engineering Physics.
文摘Imposing appropriate numerical boundary conditions at the symmetrical center r=0 is vital when computing compressible fluids with radial symmetry.Extrapolation and other traditional techniques are often employed,but spurious numerical oscillations or wall-heating phenomena can occur.In this paper,we emphasize that because of the conservation property,the updating formula of the boundary cell average can coincide with the one for interior cell averages.To achieve second-order accuracy both in time and space,we associate obtaining the inner boundary value at r=0 with the resolution of the corresponding one-sided generalized Riemann problem(GRP).Acoustic approximation is applied in this process.It creates conditions to avoid the singularity of type 1/r and aids in obtaining the value of the singular quantity using L’Hospital’s rule.Several challenging scenarios are tested to demonstrate the effectiveness and robustness of our approach.
基金the Key Projects in the National Science & Technology Pillar Program during the Twelfth Five-year Plan Period(Grant No.2015BAK17B01)Science Foundation of Institute of Engineering Mechanics,CEA under Grant No.2014B10+1 种基金Natural Science Foundation of Heilongjiang Province of China under Grant No.LC201403National Natural Science Foundation under Grant No.51378479 and No.51108431
文摘Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41175060 and 41375052)the Key Program of the Chinese Academy of Sciences (Grant No. KZZD-EW-05)+1 种基金the National Basic Research Program of China (Grant No. 2013CB430105)the State Key Laboratory of Severe Weather (LASW),Chinese Academy of Meteorological Sciences (Grant No. 2013LASW-A06)
文摘NCEP-NCAR reanalysis data were used to analyze the characteristics and evolution mechanism of convective and symmetric instability before and during a heavy rainfall event that occurred in Beijing on 21 July 2012.Approximately twelve hours before the rainstorm,the atmosphere was mainly dominated by convective instability in the lower level of 900-800 hPa.The strong southwesterly low-level jet conveyed the moist and warm airflow continuously to the area of torrential rain,maintaining and enhancing the unstable energy.When the precipitation occurred,unstable energy was released and the convective instability weakened.Meanwhile,due to the baroclinicity enhancement in the atmosphere,the symmetric instability strengthened,maintaining and promoting the subsequent torrential rain.Deriving the convective instability tendency equation demonstrated that the barotropic component of potential divergence and the advection term played a major role in enhancing the convective instability before the rainstorm.Analysis of the tendency equation of moist potential vorticity showed that the coupled term of vertical vorticity and the baroclinic component of potential divergence was the primary factor influencing the development of symmetric instability during the precipitation.Comparing the effects of these factors on convective instability and symmetric instability showed some correlation.
文摘From the formulas of the conjugate gradient, a similarity between a symmetric positive definite (SPD) matrix A and a tridiagonal matrix B is obtained. The elements of the matrix B are determined by the parameters of the conjugate gradient. The computation of eigenvalues of A is then reduced to the case of the tridiagonal matrix B. The approximation of extreme eigenvalues of A can be obtained as a 'by-product' in the computation of the conjugate gradient if a computational cost of O(s) arithmetic operations is added, where s is the number of iterations This computational cost is negligible compared with the conjugate gradient. If the matrix A is not SPD, the approximation of the condition number of A can be obtained from the computation of the conjugate gradient on AT A. Numerical results show that this is a convenient and highly efficient method for computing extreme eigenvalues and the condition number of nonsingular matrices.