A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentration...A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2] [O2], where the rate constant kif is very slow. Thus, the first time steplength Δt1 = 10-6mol·L-1/max R results long in sec. After the first step, the sequences of all the following Δt’s are very short, in μsec. So, the first time steplength Δt1 can be considered as ignition delay time. 4) It is assumed that explosion corresponds to ignition delay time Δt1 t1 = 10 sec., can be considered as explosion limit curve. This curve is compared to the corresponding one obtained by the accurate computational studies of [2].展开更多
In the traditional radar unmanned aerial vehicle(UAV)detection process,coherent integration and micro-Doppler(m-D)parameter estimation are carried out separately.The target discrimination process needs to obtain the p...In the traditional radar unmanned aerial vehicle(UAV)detection process,coherent integration and micro-Doppler(m-D)parameter estimation are carried out separately.The target discrimination process needs to obtain the position information of the target,which will lose energy.In this paper,a long time integration method of radar signal based on rotating target radon Fourier transform(RTRFT)is proposed.This method modifies the distance and frequency terms in the traditional generalized radon Fourier transform(GRFT),and adds the frequency sinusoidal modulation term.Then,based on the cardinality balanced multi-target multi-Bernoulli(CBMeMBer)filter,the position of the target is detected in the high-dimensional space obtained by RTRFT.This method can combine coherent integration and micro-motion parameter estimation.Simulation experiments show that the proposed method can estimate the main translational parameters and rotational micro-motion parameters of the target while detecting the target,and the target detection performance is improved.展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
Time delay and integration （TDI） charge coupled device （CCD） noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the comp...Time delay and integration （TDI） charge coupled device （CCD） noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the complete sources of CCD noise, we study the effects of TDI operation mode on noise, and the relationship between different types of noise and number of the TDI stage. Then we propose a new technique to identify and measure sources of TDI CCD noise employing mathematical statistics theory, where theoretical analysis shows that noise estimated formulation converges well. Finally, we establish a testing platform to carry out experiments, and a standard TDI CCD is calibrated by using the proposed method. The experimental results show that the noise analysis and measurement methods presented in this paper are useful for modeling TDI CCDs.展开更多
The integrated optical true time delay phased array antenna system has the advantages of high bandwidth,small size,low loss and strong antiinterference capability,etc.The high integration of the optically controlled p...The integrated optical true time delay phased array antenna system has the advantages of high bandwidth,small size,low loss and strong antiinterference capability,etc.The high integration of the optically controlled phased array antenna system is a necessary trend for the future development of the phased array,and it is also a major focus and difficulty in the current research of integrated microwave photonics.This paper firstly introduces the basic principle and development history of optical true time delay phased array antenna system based on microwave photonics,and briefly introduces the main implementation methods and integration platform of optical true time delay.Then,the application and development prospect of optical true time delay technology in beam control of phased array antenna system are mainly presented.Finally,according to the current research progress,the possible research directions of integrated optically controlled phased array antenna systems in the future are proposed.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy...Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.展开更多
The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ...The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.展开更多
A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,...A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,and dispersion characteristics are achieved using appropriate values of the parameter.The proposed scheme has second-order accuracy with and without physical damping.Moreover,its stability,accuracy,and dispersion are analyzed.In addition,its performance is demonstrated by the two-dimensional scalar wave problem,the single-degree-of-freedom problem,two degrees-of-freedom spring system,and beam with boundary constraints.The wave propagation problem is solved in the high frequency wave regime to demonstrate the advantage of the proposed scheme.When the proposed scheme is applied to solve the wave problem,more accurate solutions than those of other methods are obtained by using the appropriate value of the parameter.For the single-degree-offreedom system,two degrees-of-freedom system,and the time responses of beam,the proposed scheme can be used effectively owing to its high accuracy and lower computational cost.展开更多
Long-time coherent integration(LTCI)is an effective way for radar maneuvering target detection,but it faces the problem of a large number of search parameters and large amount of calculation.Realizing the simultaneous...Long-time coherent integration(LTCI)is an effective way for radar maneuvering target detection,but it faces the problem of a large number of search parameters and large amount of calculation.Realizing the simultaneous compensation of the range and Doppler migrations in complex clutter back-ground,and at the same time improving the calculation efficiency has become an urgent problem to be solved.The sparse transformation theory is introduced to LTCI in this paper,and a non-parametric searching sparse LTCI(SLTCI)based maneuvering target detection method is proposed.This method performs time reversal(TR)and second-order Keystone transform(SKT)in the range frequency&slow-time data to complete high-order range walk compensation,and achieves the coherent integra-tion of maneuvering target across range and Doppler units via the robust sparse fractional Fourier transform(RSFRFT).It can compensate for the nonlinear range migration caused by high-order motion.S-band and X-band radar data measured in sea clutter background are used to verify the detection performance of the proposed method,which can achieve better detection performance of maneuvering targets with less computational burden compared with several popular integration methods.展开更多
Mason Reset(MR),a groundbreaking invention by Clesson E.Mason in 1930 that later became a part of“the universal approach to process control instrumentation”,is revisited in this paper and is shown to consists of thr...Mason Reset(MR),a groundbreaking invention by Clesson E.Mason in 1930 that later became a part of“the universal approach to process control instrumentation”,is revisited in this paper and is shown to consists of three actions:fast(errorcorrection),medium(negative feedback for expanded proportional band)and slow(reset for zero steady-state error).The focus of the paper is on the reset action,generated from a positive feedback loop,and its underlying principles with profound implications to our understanding and practice of automatic control,both basic and advanced.For example,we note that reset control and integral control,contrary to common belief,differ fundamentally in design principle and in practicality.Such difference comes to a head in the event of integrator windup:while reset windup is a problem of actuator saturation,the integrator windup is a runaway situation due to controller instability.In fact,there is no advantage gained in replacing MR with an integrator.In other words,one should not integrate the error directly as in standard PID,since doing so makes the closed-loop system internally unstable.With MR-based control formulated in this paper,there is no such threat of instability and,therefore,no need for any anti-windup mechanisms.Furthermore,the integral control is made scalable in this framework as a tradeoff between the steady-state accuracy and the controller stability.This leads to a novel MR-based control design,scalable in gain and in time to accommodate various process characteristics and design specifications.Simple in construction and transparent in principle,this MR-based control,as a basic framework of design,is readily deployable in scale.展开更多
Stratigraphic correlations are essential for the fine-scale characterization of reservoirs.However,conventional data-driven methods that rely solely on log data struggle to construct isochronous stratigraphic framewor...Stratigraphic correlations are essential for the fine-scale characterization of reservoirs.However,conventional data-driven methods that rely solely on log data struggle to construct isochronous stratigraphic frameworks for complex sedimentary environments and multi-source geological settings.In response,this study proposed an intelligent,automatic,log-seismic integrated stratigraphic correlation method that incorporates wavelet frequency-division transform(WFT)and dynamic time warping(DTW)(also referred to as the WFT-DTW method).This approach integrates seismic data as constraints into stratigraphic correlations,enabling accurate tracking of the seismic marker horizons through WFT.Under the constraints of framework construction,a DTW algorithm was introduced to correlate sublayer boundaries automatically.The effectiveness of the proposed method was verified through a stratigraphic correlation experiment on the SA0 Formation of the Xingshugang block in the Lasaxing oilfield,the Songliao Basin,China.In this block,the target layer exhibits sublayer thicknesses ranging from 5 m to 8 m,an average sandstone thickness of 2.1 m,and pronounced heterogeneity.The verification using 1760 layers in 160 post-test wells indicates that the WFT-DTW method intelligently compared sublayers in zones with underdeveloped faults and distinct marker horizons.As a result,the posterior correlation of 1682 layers was performed,with a coincidence rate of up to 95.6%.The proposed method can complement manual correlation efforts while also providing valuable technical support for the lithologic and sand body characterization of reservoirs.展开更多
We present a comprehensive extension of the integral of first passage times(IFS)method to investigate the adsorption kinetics of polymers with multiple binding sites on planar surfaces.While effective for single-point...We present a comprehensive extension of the integral of first passage times(IFS)method to investigate the adsorption kinetics of polymers with multiple binding sites on planar surfaces.While effective for single-point adsorption,the original IFS method was limited in capturing the complex kinetics of multi-point adsorption due to inadequate reaction coordinates and theoretical frameworks.Our approach introduces a center-of-mass-based reaction coordinate and a generalized kinetic model that accounts for multi-barrier free energy landscapes characteristic of collective polymer diffusion and binding.This theoretical advancement,implemented using the adaptive bias force method for efficient sampling,enables prediction of adsorption kinetics across timescales from nanoseconds to seconds.Our results demonstrate that adsorption behavior is governed by two key factors:the number of binding monomers primarily controls desorption barriers and long-term stability,while the configuration of pre-adsorbed layers significantly modulates both adsorption and desorption rates.Polymers with three or more binding sites exhibit effectively irreversible adsorption due to exponentially increasing desorption barriers,whereas different adsorbed layer configurations lead to distinct equilibrium coverages and kinetic profiles.This extended IFS framework provides critical insights for designing functional surfaces in nanoscale sensing,macromolecular recognition,and tailored polymeric coatings where precise control over adsorption kinetics is essential.展开更多
The Lagrangian integral time scale(LITS)is a crucial characteristic for investigating the changes in fluid dynamics induced by the chaotic nature,and the finitetime Lyapunov exponent(FTLE)serves as a key measure in th...The Lagrangian integral time scale(LITS)is a crucial characteristic for investigating the changes in fluid dynamics induced by the chaotic nature,and the finitetime Lyapunov exponent(FTLE)serves as a key measure in the analysis of chaos.In this study,a new LITS model with an explicit theoretical basis and broad applicability is proposed based on the FTLE,along with a verification and evaluation criterion grounded in the Lagrangian velocity correlation coefficient.The model is used to cavitating the flow around a Clark-Y hydrofoil,and the LITS is investigated.It leads to the determination of model constants applicable to cavitating flow.The model is evaluated by the Lagrangian velocity correlation coefficient in comparison with other solution methods.All the results show that the LITS model can offer a new perspective and a new approach for studying the changes in fluid dynamics from a Lagrangian viewpoint.展开更多
For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynam...For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynamic equation can be transformed into Hamilton system from Lagrange system on the basis of the original variable. Under the whole state, the problem discussed can be described from a new view, and the equation can be precisely solved by, the time precise integration method established in linear dynamic system. A numerical example shows the effectiveness of the method.展开更多
A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved ...A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.展开更多
A new class of algorithms for trails lent finite element structural dynamical analysis which is amenable to all efficient implementation inl parallel computers (especially Massively Parallel Computers) is proposed. Th...A new class of algorithms for trails lent finite element structural dynamical analysis which is amenable to all efficient implementation inl parallel computers (especially Massively Parallel Computers) is proposed. The suitability of the method for parallel computation stems from the fact that, gived an arbitrary partition of the finite element mesh, each element in the partition can be processed over a time step independently and simultaneously with the rest, and no global equation solving effort is involved. Although the Proposed EBE time integration algorithms are shown to have the structure of an explicit scheme, they are unconditionally stable over a certain range of the algorithmic parameter.展开更多
In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thi...In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thing is to decompose a dynamic system into two sub-systems,a small-scale low-frequency system and a high-frequency system.The former can be analytically and efficiently solved with the mode superposition method,and the latter is dealt with through a time integration method such as the Newmark method.The summation of the responses of these two sub-systems is the responses of the original dynamic system.It is concluded that,with little sacrifice of efficiency,the combination method based on the strategy is more accurate than the combined time integration method,but it has the same accuracy order as that of the combined method.Numerical experiments validate the effectiveness of the proposed strategy.展开更多
Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in trac...Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
文摘A simple algorithm is proposed for step-by-step time integration of stiff ODEs in Chemical Kinetics. No predictor-corrector technique is used within each step of the algorithm. It is assumed that species concentrations less than 10-6 mol·L-1 do not activate any chemical reaction. So, within each step, the time steplength Δt of the algorithm is determined from the fastest reaction rate maxR by the formula Δt = 10-6mol·L-1/max R. All the reversible elementary reactions occur simultaneously;however, by a simple book-keeping technique, the updating of species concentrations, within each step of the algorithm, is performed within each elementary reaction separately. The above proposed simple algorithm for Chemical Kinetics is applied to a simple model for hydrogen combustion with only five reversible elementary reactions (Initiation, Propagation, First and Second Branching, Termination by wall destruction) with six species (H2, O2, H, O, HO, H2O). These five reversible reactions are recommended in the literature as the most significant elementary reactions of hydrogen combustion [1] [2]. Based on the proposed here simple algorithm for Chemical Kinetics, applied to the global mechanism of proposed five reversible elementary reactions for hydrogen combustion, a simple and short computer program has been developed with only about 120 Fortran instructions. By this proposed program, the following are obtained: 1) The total species concentration of hydrogen combustion, starting from the sum of initial reactants concentrations [H2] + [O2], gradually diminishes, due to termination reaction by wall destruction, and tends to the final concentration of the product [H2O], that is to the 2/3 of its initial value, in accordance to the established overall stoichiometric reaction of hydrogen combustion 2H2 + O2 → 2H2O. 2) Time-histories for concentrations of main species H2, O2, H, H2O of hydrogen combustion, in explosion and equilibrium regions, obtained by the proposed program, are compared to corresponding ones obtained by accurate computational studies of [3]. 3) In the first step of the algorithm, the only nonzero species concentrations are those of reactants [H2], [O2]. So, the maximum reaction rate is that of the forward initiation reaction max R = Rif = kif[H2] [O2], where the rate constant kif is very slow. Thus, the first time steplength Δt1 = 10-6mol·L-1/max R results long in sec. After the first step, the sequences of all the following Δt’s are very short, in μsec. So, the first time steplength Δt1 can be considered as ignition delay time. 4) It is assumed that explosion corresponds to ignition delay time Δt1 t1 = 10 sec., can be considered as explosion limit curve. This curve is compared to the corresponding one obtained by the accurate computational studies of [2].
文摘In the traditional radar unmanned aerial vehicle(UAV)detection process,coherent integration and micro-Doppler(m-D)parameter estimation are carried out separately.The target discrimination process needs to obtain the position information of the target,which will lose energy.In this paper,a long time integration method of radar signal based on rotating target radon Fourier transform(RTRFT)is proposed.This method modifies the distance and frequency terms in the traditional generalized radon Fourier transform(GRFT),and adds the frequency sinusoidal modulation term.Then,based on the cardinality balanced multi-target multi-Bernoulli(CBMeMBer)filter,the position of the target is detected in the high-dimensional space obtained by RTRFT.This method can combine coherent integration and micro-motion parameter estimation.Simulation experiments show that the proposed method can estimate the main translational parameters and rotational micro-motion parameters of the target while detecting the target,and the target detection performance is improved.
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
基金Project supported by the National High Technology Research and Development Program of China (Grant No. 2006AA06A208)
文摘Time delay and integration （TDI） charge coupled device （CCD） noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the complete sources of CCD noise, we study the effects of TDI operation mode on noise, and the relationship between different types of noise and number of the TDI stage. Then we propose a new technique to identify and measure sources of TDI CCD noise employing mathematical statistics theory, where theoretical analysis shows that noise estimated formulation converges well. Finally, we establish a testing platform to carry out experiments, and a standard TDI CCD is calibrated by using the proposed method. The experimental results show that the noise analysis and measurement methods presented in this paper are useful for modeling TDI CCDs.
基金supported by Fund of State Key Laboratory of IPOC(BUPT)(No.IPOC2021ZT16),China.
文摘The integrated optical true time delay phased array antenna system has the advantages of high bandwidth,small size,low loss and strong antiinterference capability,etc.The high integration of the optically controlled phased array antenna system is a necessary trend for the future development of the phased array,and it is also a major focus and difficulty in the current research of integrated microwave photonics.This paper firstly introduces the basic principle and development history of optical true time delay phased array antenna system based on microwave photonics,and briefly introduces the main implementation methods and integration platform of optical true time delay.Then,the application and development prospect of optical true time delay technology in beam control of phased array antenna system are mainly presented.Finally,according to the current research progress,the possible research directions of integrated optically controlled phased array antenna systems in the future are proposed.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China （G2000048702, 2003CB716707）the National Science Fund for Distinguished Young Scholars （10025208）+1 种基金 the National Natural Science Foundation of China （Key Program） （10532040） the Research Fund for 0versea Chinese （10228028）.
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
基金National Natural Science Foundation of China under Grant Nos.51639006 and 51725901
文摘Finite element（FE） is a powerful tool and has been applied by investigators to real-time hybrid simulations（RTHSs）. This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time（TET） decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method（CDM）, the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ（λ = 4） method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.
文摘The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.
文摘A direct time integration scheme based on Gauss-Legendre quadrature is proposed to solve problems in linear structural dynamics.The proposed method is a oneparameter non-dissipative scheme.Improved stability,accuracy,and dispersion characteristics are achieved using appropriate values of the parameter.The proposed scheme has second-order accuracy with and without physical damping.Moreover,its stability,accuracy,and dispersion are analyzed.In addition,its performance is demonstrated by the two-dimensional scalar wave problem,the single-degree-of-freedom problem,two degrees-of-freedom spring system,and beam with boundary constraints.The wave propagation problem is solved in the high frequency wave regime to demonstrate the advantage of the proposed scheme.When the proposed scheme is applied to solve the wave problem,more accurate solutions than those of other methods are obtained by using the appropriate value of the parameter.For the single-degree-offreedom system,two degrees-of-freedom system,and the time responses of beam,the proposed scheme can be used effectively owing to its high accuracy and lower computational cost.
基金supported by the National Natural Science Foundation of China(62222120,61871391,U1933135)Shandong Provincial Natural Science Foundation(ZR2021YQ43).
文摘Long-time coherent integration(LTCI)is an effective way for radar maneuvering target detection,but it faces the problem of a large number of search parameters and large amount of calculation.Realizing the simultaneous compensation of the range and Doppler migrations in complex clutter back-ground,and at the same time improving the calculation efficiency has become an urgent problem to be solved.The sparse transformation theory is introduced to LTCI in this paper,and a non-parametric searching sparse LTCI(SLTCI)based maneuvering target detection method is proposed.This method performs time reversal(TR)and second-order Keystone transform(SKT)in the range frequency&slow-time data to complete high-order range walk compensation,and achieves the coherent integra-tion of maneuvering target across range and Doppler units via the robust sparse fractional Fourier transform(RSFRFT).It can compensate for the nonlinear range migration caused by high-order motion.S-band and X-band radar data measured in sea clutter background are used to verify the detection performance of the proposed method,which can achieve better detection performance of maneuvering targets with less computational burden compared with several popular integration methods.
文摘Mason Reset(MR),a groundbreaking invention by Clesson E.Mason in 1930 that later became a part of“the universal approach to process control instrumentation”,is revisited in this paper and is shown to consists of three actions:fast(errorcorrection),medium(negative feedback for expanded proportional band)and slow(reset for zero steady-state error).The focus of the paper is on the reset action,generated from a positive feedback loop,and its underlying principles with profound implications to our understanding and practice of automatic control,both basic and advanced.For example,we note that reset control and integral control,contrary to common belief,differ fundamentally in design principle and in practicality.Such difference comes to a head in the event of integrator windup:while reset windup is a problem of actuator saturation,the integrator windup is a runaway situation due to controller instability.In fact,there is no advantage gained in replacing MR with an integrator.In other words,one should not integrate the error directly as in standard PID,since doing so makes the closed-loop system internally unstable.With MR-based control formulated in this paper,there is no such threat of instability and,therefore,no need for any anti-windup mechanisms.Furthermore,the integral control is made scalable in this framework as a tradeoff between the steady-state accuracy and the controller stability.This leads to a novel MR-based control design,scalable in gain and in time to accommodate various process characteristics and design specifications.Simple in construction and transparent in principle,this MR-based control,as a basic framework of design,is readily deployable in scale.
基金funded by the Major Science and Technology Project of China National Petroleum Corporation(No.2023ZZ22YJ01).
文摘Stratigraphic correlations are essential for the fine-scale characterization of reservoirs.However,conventional data-driven methods that rely solely on log data struggle to construct isochronous stratigraphic frameworks for complex sedimentary environments and multi-source geological settings.In response,this study proposed an intelligent,automatic,log-seismic integrated stratigraphic correlation method that incorporates wavelet frequency-division transform(WFT)and dynamic time warping(DTW)(also referred to as the WFT-DTW method).This approach integrates seismic data as constraints into stratigraphic correlations,enabling accurate tracking of the seismic marker horizons through WFT.Under the constraints of framework construction,a DTW algorithm was introduced to correlate sublayer boundaries automatically.The effectiveness of the proposed method was verified through a stratigraphic correlation experiment on the SA0 Formation of the Xingshugang block in the Lasaxing oilfield,the Songliao Basin,China.In this block,the target layer exhibits sublayer thicknesses ranging from 5 m to 8 m,an average sandstone thickness of 2.1 m,and pronounced heterogeneity.The verification using 1760 layers in 160 post-test wells indicates that the WFT-DTW method intelligently compared sublayers in zones with underdeveloped faults and distinct marker horizons.As a result,the posterior correlation of 1682 layers was performed,with a coincidence rate of up to 95.6%.The proposed method can complement manual correlation efforts while also providing valuable technical support for the lithologic and sand body characterization of reservoirs.
基金financially supported by the National Natural Science Foundation of China(No.12374207)the Natural Science Foundation of Jiangsu Province(No.BK20233001)supported by the Big Data Computing Center of Southeast University。
文摘We present a comprehensive extension of the integral of first passage times(IFS)method to investigate the adsorption kinetics of polymers with multiple binding sites on planar surfaces.While effective for single-point adsorption,the original IFS method was limited in capturing the complex kinetics of multi-point adsorption due to inadequate reaction coordinates and theoretical frameworks.Our approach introduces a center-of-mass-based reaction coordinate and a generalized kinetic model that accounts for multi-barrier free energy landscapes characteristic of collective polymer diffusion and binding.This theoretical advancement,implemented using the adaptive bias force method for efficient sampling,enables prediction of adsorption kinetics across timescales from nanoseconds to seconds.Our results demonstrate that adsorption behavior is governed by two key factors:the number of binding monomers primarily controls desorption barriers and long-term stability,while the configuration of pre-adsorbed layers significantly modulates both adsorption and desorption rates.Polymers with three or more binding sites exhibit effectively irreversible adsorption due to exponentially increasing desorption barriers,whereas different adsorbed layer configurations lead to distinct equilibrium coverages and kinetic profiles.This extended IFS framework provides critical insights for designing functional surfaces in nanoscale sensing,macromolecular recognition,and tailored polymeric coatings where precise control over adsorption kinetics is essential.
基金Project supported by the Key Project of the National Natural Science Foundation of China(No.52336001)the Natural Science Foundation of Zhejiang Province of China(No.LR20E090001)。
文摘The Lagrangian integral time scale(LITS)is a crucial characteristic for investigating the changes in fluid dynamics induced by the chaotic nature,and the finitetime Lyapunov exponent(FTLE)serves as a key measure in the analysis of chaos.In this study,a new LITS model with an explicit theoretical basis and broad applicability is proposed based on the FTLE,along with a verification and evaluation criterion grounded in the Lagrangian velocity correlation coefficient.The model is used to cavitating the flow around a Clark-Y hydrofoil,and the LITS is investigated.It leads to the determination of model constants applicable to cavitating flow.The model is evaluated by the Lagrangian velocity correlation coefficient in comparison with other solution methods.All the results show that the LITS model can offer a new perspective and a new approach for studying the changes in fluid dynamics from a Lagrangian viewpoint.
文摘For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynamic equation can be transformed into Hamilton system from Lagrange system on the basis of the original variable. Under the whole state, the problem discussed can be described from a new view, and the equation can be precisely solved by, the time precise integration method established in linear dynamic system. A numerical example shows the effectiveness of the method.
文摘A method is presented that coordinates the calculation of the displacement, velocity and acceleration of structures within the time-steps of different types of step-by-step integration. The dynamic equation is solved using an energy equation and the calculating data of the original method. The method presented is better than the original method in terms of calculating postulations and is in better conformity with the system's movement. Take the Wilson-θ method as an example. By using the coordination process, the calculation precision has been greatly im proved (reducing the errors by approximately 90% ), and the greater part of overshooting of the calculation result has been eliminated. The study suggests that the mal-coordination of the motion parameters within the time-step is the major factor that contributes to the result errors of step-by-step integration for the dynamic equation.
文摘A new class of algorithms for trails lent finite element structural dynamical analysis which is amenable to all efficient implementation inl parallel computers (especially Massively Parallel Computers) is proposed. The suitability of the method for parallel computation stems from the fact that, gived an arbitrary partition of the finite element mesh, each element in the partition can be processed over a time step independently and simultaneously with the rest, and no global equation solving effort is involved. Although the Proposed EBE time integration algorithms are shown to have the structure of an explicit scheme, they are unconditionally stable over a certain range of the algorithmic parameter.
基金This work was supported by the National Natural Science Foundation of China(Grant Nos.11872090 and 12172023).
文摘In this article,the mode superposition method is combined with a time integration method like the trapezoidal rule to improve solution accuracy for linear dynamic systems.In this combination strategy,the essential thing is to decompose a dynamic system into two sub-systems,a small-scale low-frequency system and a high-frequency system.The former can be analytically and efficiently solved with the mode superposition method,and the latter is dealt with through a time integration method such as the Newmark method.The summation of the responses of these two sub-systems is the responses of the original dynamic system.It is concluded that,with little sacrifice of efficiency,the combination method based on the strategy is more accurate than the combined time integration method,but it has the same accuracy order as that of the combined method.Numerical experiments validate the effectiveness of the proposed strategy.
基金supported by the National Natural Science Foundation of China(Grant Numbers 11872090,11672019,11472035).
文摘Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
