A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the c...A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.展开更多
This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differen...This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differences methods with back-propagation algorithm for updating the parameters continuously on the basis of recent data. This method can make the neural network model fit the recent characteristic of the time series as close as possible, therefore improves the prediction accuracy. We built models and made predictions for the sunspot series. The prediction results of adaptive modeling method are better than that of non-adaptive modeling methods.展开更多
The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell's equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method i...The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method is suitable to calculate the ground state of the quantum systems,it has been improved to calculate the higher excited states directly.The improvement is based on modifying the iterative process involved in this method to include two procedures.The first is known as cooling steps and the second is known as a heating step.By determining the required length of the cooling iteration steps using suitable excitation energy estimate,and repeating these two procedures using suitable initial guess function for sufficient times.This modified iteration will lead automatically to the desired excited state.In the two dimensional finite rectangular well potential problem both of the suitable excitation energy and the suitable initial guess wave function are calculated analytically using the separation of variables technique.展开更多
In the measurement of the Newtonian gravitational constant G with the time-of-swing method,the influence of the Earth's rotation has been roughly estimated before,which is far beyond the current experimental preci...In the measurement of the Newtonian gravitational constant G with the time-of-swing method,the influence of the Earth's rotation has been roughly estimated before,which is far beyond the current experimental precision.Here,we present a more complete theoretical modeling and assessment process.To figure out this effect,we use the relativistic Lagrangian expression to derive the motion equations of the torsion pendulum.With the correlation method and typical parameters,we estimate that the influence of the Earth's rotation on G measurement is far less than 1 ppm,which may need to be considered in the future high-accuracy experiments of determining the gravitational constant G.展开更多
In contrast to the solutions of applied mathematics to Zeno’s paradoxes, I focus on the concept of motion and show that, by distinguishing two different forms of motion, Zeno’s apparent paradoxes are not paradoxical...In contrast to the solutions of applied mathematics to Zeno’s paradoxes, I focus on the concept of motion and show that, by distinguishing two different forms of motion, Zeno’s apparent paradoxes are not paradoxical at all. Zeno’s paradoxes indirectly prove that distances are not composed of extensionless points and, in general, that a higher dimension cannot be completely composed of lower ones. Conversely, lower dimensions can be understood as special cases of higher dimensions. To illustrate this approach, I consider Cantor’s only apparent proof that the real numbers are uncountable. However, his widely accepted indirect proof has the disadvantage that it depends on whether there is another way to make the real numbers countable. Cantor rightly assumes that there can be no smallest number between 0 and 1, and therefore no beginning of counting. For this reason he arbitrarily lists the real numbers in order to show with his diagonal method that this list can never be complete. The situation is different if we start with the largest number between 0 and 1 (0.999…) and use the method of an inverted triangle, which can be understood as a special fractal form. Here we can construct a vertical and a horizontal stratification with which it is actually possible to construct all real numbers between 0 and 1 without exception. Each column is infinite, and each number in that column is the starting point of a new triangle, while each row is finite. Even in a simple sine curve, we experience finiteness with respect to the y-axis and infinity with respect to the x-axis. The first parts of this article show that Zeno’s assumptions contradict the concept of motion as such, so it is not surprising that this misconstruction leads to contradictions. In the last part, I discuss Cantor’s diagonal method and explain the method of an inverted triangle that is internally structured like a fractal by repeating this inverted triangle at each column. The consequence is that we encounter two very different methods of counting. Vertically it is continuous, horizontally it is discrete. While Frege, Tarski, Cantor, Gödel and the Vienna Circle tried to derive the higher dimension from the lower, a procedure that always leads to new contradictions and antinomies (Tarski, Russell), I take the opposite approach here, in which I derive the lower dimension from the higher. This perspective seems to fail because Tarski, Russell, Wittgenstein, and especially the Vienna Circle have shown that the completeness of the absolute itself is logically contradictory. For this reason, we agree with Hegel in assuming that we can never fully comprehend the Absolute, but only its particular manifestations—otherwise we would be putting ourselves in the place of the Absolute, or even God. Nevertheless, we can understand the Absolute in its particular expressions, as I will show with the modest example of the triangle proof of the combined horizontal and vertical countability of the real numbers, which I developed in rejection of Cantor’s diagonal proof. .展开更多
The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of c...The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.展开更多
This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo's amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^(3+):YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp(α'L),which is previously employed to describe the relationship between echo's efficiency and thickness,should be modified as 1.3 · 0.09 exp(2.4 ·α'L) with the propagation of echo considered.展开更多
In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald...In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.展开更多
In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler ti...In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler time-stepping. Then we apply the abstract framework of to prove its long-time convergence. Finally, a numerical example for solving driven cavity flows is given.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the ...Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to展开更多
The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites...The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites.By now,dual one-way ranging(DOWR)and two-way time transfer(TWTT)are generally integrated in the same software and hardware system to meet the limitations of small satellites in terms of size,weight and power(SWaP)consumption.However,studies show that pseudo-noise regenerative ranging(PNRR)performs better than DOWR if some advanced implementation technologies are employed.Besides,PNRR has no requirement on time synchronization.To apply PNRR to small satellites,and meanwhile,meet the demand for time difference measurement,we propose the round-way time difference measurement,which can be combined with PNRR to form a new integrated system without exceeding the limits of SWaP.The new integrated system can provide distributed small satellite systems with on-orbit high-accuracy and high-precision distance measurement and time difference measurement in real time.Experimental results show that the precision of ranging is about 1.94 cm,and that of time difference measurement is about 78.4 ps,at the signal to noise ratio of 80 dBHz.展开更多
Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, rob...Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, robustness and safe maimer of operation in cases of failure have made them attractive devices for semi-active real-time control in civil, aerospace and automotive applications. Time response characteristic is one of the most important technical performances of MR dampers, and response time directly affects the control frequency, application range and the actual effect of MR dampers. In this study, one kind of finite difference solution for predicting the response time ofmagneto-rheological dampers from "off-state" to "on-state" is put forward. A laminar flow model is used to describe the flow in the MR valve, and a bi-viscous fluid flow model is utilized to describe the relationship of shear stress and shear rate of MR fluid. An explicit difference format is used to diseretize the Novior-Stoks equation, and stability condition of this algorithm is built by Von-Neumann stability criterion. The pressure gradient along the flow duct is solved by a dichotomy algorithm with iteration, and the changing curve of the damping force versus time of MR damper is obtained as well. According to the abovementioned numerical algorithm, the damping forces versus time curves from "off-state" to "on-state" of a cylindrical piston type MR damper are computed. Moreover, the MR damper is installed in a material test system(MTS), the magnetic field in the wire circles of the MR damper is "triggered" when the MR damper is imposed to do a constant speed motion, and the damping force curves are recorded. The comparison between numerical results and experimental results indicates that this finite difference algorithm can be used to estimate the response time delay of MR dampers.展开更多
The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the ...The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.展开更多
Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this ...Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this paper, which is one of the most basic ways for diagnosis of insulation condition. Methods Ultra high frequency(UHF) PD detection method by using internal sensors has been proved efficient, because it may avoid the disturbance of corona, but the sensor installation of this method will be limited by the structure and operation condition of GIS. There are some of electromagnetic (E-M) waves leak from the place of insulation spacer, therefore, the external sensors UHF measuring PD technique is applied, which isn't limited by the operation condition of GIS. Results This paper analyzes propagated electromagnetic (E-M) waves of partial discharge pulse excited by using the finite-difference time-domain (FDTD) method. The signal collected at the outer point is more complex than that of the inner point, and the signals' amplitude of outer is about half of the inner, because it propagates through spacer and insulation slot. Set up UHF PD measuring system. The typical PD in 252kV GIS bus bar was measured using PD detection UHF technique with external sensors. Finally, compare the results of UHF measuring technique using external sensors with the results of FDTD method simulation and the traditional IEC60270 method detection. Conclusion The results of experiment shows that the UHF technique can realize the diagnosis of insulation condition, the results of FDTD method simulation and the result UHF method detection can demonstrate each other, which gives references to further researches and application for UHF PD measuring technique.展开更多
The valuation relation of potential difference with discharging time in Electrical Discharge Machining (EDM) is analyzed theoretically and tested and verified by experiments designed in this paper and the relation bet...The valuation relation of potential difference with discharging time in Electrical Discharge Machining (EDM) is analyzed theoretically and tested and verified by experiments designed in this paper and the relation between potential difference and spark location is induced and analyzed, and proceed by experiments under the condition of onedimension.展开更多
The effect of vibratory stress relief (VSR) is usually evaluated with the indirect method of observing the change of amplitude frequency response characteristics of structures. A new kind of evaluating method of VSR...The effect of vibratory stress relief (VSR) is usually evaluated with the indirect method of observing the change of amplitude frequency response characteristics of structures. A new kind of evaluating method of VSR based on the ultrasonic time-of-arrival method (UTM), which can obtain the residual stress directly through measuring the propagation time of ultrasonic wave in the material, is presented. At first, the principle of the measuring method of residual stress based on UTM is analyzed. Then the measuring system of the method is described, which is in virtue of ultrasonic flaw detector and high-sampling-rate digital oscillograph. And a set of calibration system that contains a piece of standard specimen is also introduced. Experimental results prove the relation between the residual stress and the propagation time of ultrasonic in workpieces. Finally, the measuring and calibration systems are applied in evaluating the effect of VSR. The final test results show that the method is effective.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
We present a numerical study of the long time behavior of approxima- tion solution to the Extended Fisher-Kolmogorov equation with periodic boundary conditions. The unique solvability of numerical solution is shown. I...We present a numerical study of the long time behavior of approxima- tion solution to the Extended Fisher-Kolmogorov equation with periodic boundary conditions. The unique solvability of numerical solution is shown. It is proved that there exists a global attractor of the discrete dynamical system. Furthermore, we obtain the long-time stability and convergence of the difference scheme and the upper semicontinuity d(Ah,τ, .A) → O. Our results show that the difference scheme can effectively simulate the infinite dimensional dynamical systems.展开更多
基金the National Natural Science Foundation of China(No.11701103)the Young Top-notch Talent Program of Guangdong Province of China(No.2017GC010379)+4 种基金the Natural Science Foundation of Guangdong Province of China(No.2022A1515012147)the Project of Science and Technology of Guangzhou of China(No.202102020704)the Opening Project of Guangdong Province Key Laboratory of Computational Science at the Sun Yat-sen University of China(2021023)the Science and Technology Development Fund,Macao SAR(File No.0005/2019/A)the University of Macao of China(File Nos.MYRG2020-00035-FST,MYRG2018-00047-FST).
文摘A finite difference/spectral scheme is proposed for the time fractional Ito equation.The mass conservation and stability of the numerical solution are deduced by the energy method in the L^(2)norm form.To reduce the computation costs,the fast Fourier transform technic is applied to a pair of equivalent coupled differential equations.The effectiveness of the proposed algorithm is verified by the first numerical example.The mass conservation property and stability statement are confirmed by two other numerical examples.
文摘This paper discusses the modeling method of time series with neural network. In order to improve the adaptability of direct multi-step prediction models, this paper proposes a method of combining the temporal differences methods with back-propagation algorithm for updating the parameters continuously on the basis of recent data. This method can make the neural network model fit the recent characteristic of the time series as close as possible, therefore improves the prediction accuracy. We built models and made predictions for the sunspot series. The prediction results of adaptive modeling method are better than that of non-adaptive modeling methods.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
文摘The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell's equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
文摘The higher excited states for two dimensional finite rectangular well potential are calculated numerically,by solving the Schrödinger equation using the finite difference time domain method.Although,this method is suitable to calculate the ground state of the quantum systems,it has been improved to calculate the higher excited states directly.The improvement is based on modifying the iterative process involved in this method to include two procedures.The first is known as cooling steps and the second is known as a heating step.By determining the required length of the cooling iteration steps using suitable excitation energy estimate,and repeating these two procedures using suitable initial guess function for sufficient times.This modified iteration will lead automatically to the desired excited state.In the two dimensional finite rectangular well potential problem both of the suitable excitation energy and the suitable initial guess wave function are calculated analytically using the separation of variables technique.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11575160 and 11805074)the Postdoctoral Science Foundation of China(Grant Nos.2017M620308 and 2018T110750).
文摘In the measurement of the Newtonian gravitational constant G with the time-of-swing method,the influence of the Earth's rotation has been roughly estimated before,which is far beyond the current experimental precision.Here,we present a more complete theoretical modeling and assessment process.To figure out this effect,we use the relativistic Lagrangian expression to derive the motion equations of the torsion pendulum.With the correlation method and typical parameters,we estimate that the influence of the Earth's rotation on G measurement is far less than 1 ppm,which may need to be considered in the future high-accuracy experiments of determining the gravitational constant G.
文摘In contrast to the solutions of applied mathematics to Zeno’s paradoxes, I focus on the concept of motion and show that, by distinguishing two different forms of motion, Zeno’s apparent paradoxes are not paradoxical at all. Zeno’s paradoxes indirectly prove that distances are not composed of extensionless points and, in general, that a higher dimension cannot be completely composed of lower ones. Conversely, lower dimensions can be understood as special cases of higher dimensions. To illustrate this approach, I consider Cantor’s only apparent proof that the real numbers are uncountable. However, his widely accepted indirect proof has the disadvantage that it depends on whether there is another way to make the real numbers countable. Cantor rightly assumes that there can be no smallest number between 0 and 1, and therefore no beginning of counting. For this reason he arbitrarily lists the real numbers in order to show with his diagonal method that this list can never be complete. The situation is different if we start with the largest number between 0 and 1 (0.999…) and use the method of an inverted triangle, which can be understood as a special fractal form. Here we can construct a vertical and a horizontal stratification with which it is actually possible to construct all real numbers between 0 and 1 without exception. Each column is infinite, and each number in that column is the starting point of a new triangle, while each row is finite. Even in a simple sine curve, we experience finiteness with respect to the y-axis and infinity with respect to the x-axis. The first parts of this article show that Zeno’s assumptions contradict the concept of motion as such, so it is not surprising that this misconstruction leads to contradictions. In the last part, I discuss Cantor’s diagonal method and explain the method of an inverted triangle that is internally structured like a fractal by repeating this inverted triangle at each column. The consequence is that we encounter two very different methods of counting. Vertically it is continuous, horizontally it is discrete. While Frege, Tarski, Cantor, Gödel and the Vienna Circle tried to derive the higher dimension from the lower, a procedure that always leads to new contradictions and antinomies (Tarski, Russell), I take the opposite approach here, in which I derive the lower dimension from the higher. This perspective seems to fail because Tarski, Russell, Wittgenstein, and especially the Vienna Circle have shown that the completeness of the absolute itself is logically contradictory. For this reason, we agree with Hegel in assuming that we can never fully comprehend the Absolute, but only its particular manifestations—otherwise we would be putting ourselves in the place of the Absolute, or even God. Nevertheless, we can understand the Absolute in its particular expressions, as I will show with the modest example of the triangle proof of the combined horizontal and vertical countability of the real numbers, which I developed in rejection of Cantor’s diagonal proof. .
文摘The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
文摘This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo's amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^(3+):YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp(α'L),which is previously employed to describe the relationship between echo's efficiency and thickness,should be modified as 1.3 · 0.09 exp(2.4 ·α'L) with the propagation of echo considered.
文摘In this paper, three implicit finite difference methods are developed to solve one dimensional time fractional advection-diffusion equation. The fractional derivative is treated by applying right shifted Grünwald-Letnikov formula of order α ∈(0, 1). We investigate the stability analysis by using von Neumann method with mathematical induction and prove that these three proposed methods are unconditionally stable. Numerical results are presented to demonstrate the effectiveness of the schemes mentioned in this paper.
基金The project supported by Laboratory of Computational Physics,Institute of Applied Physics & Computational Mathematics,T.O.Box 80 0 9,Beijing 1 0 0 0 88
文摘In this paper, we first provide a generalized difference method for the 2-dimensional Navier-Stokes equations by combing the ideas of staggered scheme m and generalized upwind scheme in space, and by backward Euler time-stepping. Then we apply the abstract framework of to prove its long-time convergence. Finally, a numerical example for solving driven cavity flows is given.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
文摘Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference to
基金supported by the National Natural Science Foundation of China(61401389).
文摘The advancement of small satellites is promoting the development of distributed satellite systems,and for the latter,it is essential to coordinate the spatial and temporal relations between mutually visible satellites.By now,dual one-way ranging(DOWR)and two-way time transfer(TWTT)are generally integrated in the same software and hardware system to meet the limitations of small satellites in terms of size,weight and power(SWaP)consumption.However,studies show that pseudo-noise regenerative ranging(PNRR)performs better than DOWR if some advanced implementation technologies are employed.Besides,PNRR has no requirement on time synchronization.To apply PNRR to small satellites,and meanwhile,meet the demand for time difference measurement,we propose the round-way time difference measurement,which can be combined with PNRR to form a new integrated system without exceeding the limits of SWaP.The new integrated system can provide distributed small satellite systems with on-orbit high-accuracy and high-precision distance measurement and time difference measurement in real time.Experimental results show that the precision of ranging is about 1.94 cm,and that of time difference measurement is about 78.4 ps,at the signal to noise ratio of 80 dBHz.
基金supported by National Natural Science Foundation of China (Grant No. 50675106)
文摘Magneto-rhelological(MR) dampers are devices that employ rheological fluids to modify their mechanical properties. Their mechanical simplicity, high dynamic range, lower power requirements, large force capacity, robustness and safe maimer of operation in cases of failure have made them attractive devices for semi-active real-time control in civil, aerospace and automotive applications. Time response characteristic is one of the most important technical performances of MR dampers, and response time directly affects the control frequency, application range and the actual effect of MR dampers. In this study, one kind of finite difference solution for predicting the response time ofmagneto-rheological dampers from "off-state" to "on-state" is put forward. A laminar flow model is used to describe the flow in the MR valve, and a bi-viscous fluid flow model is utilized to describe the relationship of shear stress and shear rate of MR fluid. An explicit difference format is used to diseretize the Novior-Stoks equation, and stability condition of this algorithm is built by Von-Neumann stability criterion. The pressure gradient along the flow duct is solved by a dichotomy algorithm with iteration, and the changing curve of the damping force versus time of MR damper is obtained as well. According to the abovementioned numerical algorithm, the damping forces versus time curves from "off-state" to "on-state" of a cylindrical piston type MR damper are computed. Moreover, the MR damper is installed in a material test system(MTS), the magnetic field in the wire circles of the MR damper is "triggered" when the MR damper is imposed to do a constant speed motion, and the damping force curves are recorded. The comparison between numerical results and experimental results indicates that this finite difference algorithm can be used to estimate the response time delay of MR dampers.
基金Project supported partly by the Open Research Program in State Key Laboratory of Millimeter Waves of China(Grant No.K200802)partly by the National Natural Science Foundation of China(Grant No.60971122)
文摘The transmission coefficients of electromagnetic (EM) waves due to a superconductor-dielectric superlattice are numerically calculated. Shift operator finite difference time domain (SO-FDTD) method is used in the analysis. By using the SO-FDTD method, the transmission spectrum is obtained and its characteristics are investigated for different thicknesses of superconductor layers and dielectric layers, from which a stop band starting from zero frequency can be apparently observed. The relation between this low-frequency stop band and relative temperature, and also the London penetration depth at a superconductor temperature of zero degree are discussed, separately. The low-frequency stop band properties of superconductor-dielectric superlattice thus are well disclosed.
文摘Objective In order to find early latent faults and prevent catastrophic failures, diagnosis of insulation condition by measuring technique of partial discharge(PD) in gas insulated switchgear (GIS) is applied in this paper, which is one of the most basic ways for diagnosis of insulation condition. Methods Ultra high frequency(UHF) PD detection method by using internal sensors has been proved efficient, because it may avoid the disturbance of corona, but the sensor installation of this method will be limited by the structure and operation condition of GIS. There are some of electromagnetic (E-M) waves leak from the place of insulation spacer, therefore, the external sensors UHF measuring PD technique is applied, which isn't limited by the operation condition of GIS. Results This paper analyzes propagated electromagnetic (E-M) waves of partial discharge pulse excited by using the finite-difference time-domain (FDTD) method. The signal collected at the outer point is more complex than that of the inner point, and the signals' amplitude of outer is about half of the inner, because it propagates through spacer and insulation slot. Set up UHF PD measuring system. The typical PD in 252kV GIS bus bar was measured using PD detection UHF technique with external sensors. Finally, compare the results of UHF measuring technique using external sensors with the results of FDTD method simulation and the traditional IEC60270 method detection. Conclusion The results of experiment shows that the UHF technique can realize the diagnosis of insulation condition, the results of FDTD method simulation and the result UHF method detection can demonstrate each other, which gives references to further researches and application for UHF PD measuring technique.
文摘The valuation relation of potential difference with discharging time in Electrical Discharge Machining (EDM) is analyzed theoretically and tested and verified by experiments designed in this paper and the relation between potential difference and spark location is induced and analyzed, and proceed by experiments under the condition of onedimension.
基金This project is supported by National Natural Science Foundation of China(No.50305036).
文摘The effect of vibratory stress relief (VSR) is usually evaluated with the indirect method of observing the change of amplitude frequency response characteristics of structures. A new kind of evaluating method of VSR based on the ultrasonic time-of-arrival method (UTM), which can obtain the residual stress directly through measuring the propagation time of ultrasonic wave in the material, is presented. At first, the principle of the measuring method of residual stress based on UTM is analyzed. Then the measuring system of the method is described, which is in virtue of ultrasonic flaw detector and high-sampling-rate digital oscillograph. And a set of calibration system that contains a piece of standard specimen is also introduced. Experimental results prove the relation between the residual stress and the propagation time of ultrasonic in workpieces. Finally, the measuring and calibration systems are applied in evaluating the effect of VSR. The final test results show that the method is effective.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming (P1)2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金The NSF (10871055) of Chinathe Fundamental Research Funds (HEUCFL20111102)for the Central Universities
文摘We present a numerical study of the long time behavior of approxima- tion solution to the Extended Fisher-Kolmogorov equation with periodic boundary conditions. The unique solvability of numerical solution is shown. It is proved that there exists a global attractor of the discrete dynamical system. Furthermore, we obtain the long-time stability and convergence of the difference scheme and the upper semicontinuity d(Ah,τ, .A) → O. Our results show that the difference scheme can effectively simulate the infinite dimensional dynamical systems.