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 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.展开更多
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.展开更多
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].展开更多
This paper presents an innovative and effective control strategy tailored for a deregulated,diversified energy system involving multiple interconnected area.Each area integrates a unique mix of power generation techno...This paper presents an innovative and effective control strategy tailored for a deregulated,diversified energy system involving multiple interconnected area.Each area integrates a unique mix of power generation technologies:Area 1 combines thermal,hydro,and distributed generation;Area 2 utilizes a blend of thermal units,distributed solar technologies(DST),and hydro power;andThird control area hosts geothermal power station alongside thermal power generation unit and hydropower units.The suggested control system employs a multi-layered approach,featuring a blended methodology utilizing the Tilted Integral Derivative controller(TID)and the Fractional-Order Integral method to enhance performance and stability.The parameters of this hybrid TID-FOI controller are finely tuned using an advanced optimization method known as the Walrus Optimization Algorithm(WaOA).Performance analysis reveals that the combined TID-FOI controller significantly outperforms the TID and PID controllers when comparing their dynamic response across various system configurations.The study also incorporates investigation of redox flow batteries within the broader scope of energy storage applications to assess their impact on system performance.In addition,the research explores the controller’s effectiveness under different power exchange scenarios in a deregulated market,accounting for restrictions on generation ramp rates and governor hysteresis effects in dynamic control.To ensure the reliability and resilience of the presented methodology,the system transitions and develops across a broad range of varying parameters and stochastic load fluctuation.To wrap up,the study offers a pioneering control approach-a hybrid TID-FOI controller optimized via the Walrus Optimization Algorithm(WaOA)-designed for enhanced stability and performance in a complex,three-region hybrid energy system functioning within a deregulated framework.展开更多
Traditional scheduling algorithms for avionics communication have the shortcoming of messages accumulation,the efficiency and reliability of the service can be improved by combining the distributed integrated modular ...Traditional scheduling algorithms for avionics communication have the shortcoming of messages accumulation,the efficiency and reliability of the service can be improved by combining the distributed integrated modular avionics(DIMA)system with a time trigger mechanism.To further improve the utilization of system resources,the static scheduling algorithm of time triggered service is studied.By making the time trigger message schedule dispersedly,the stabilities of both the available time slots for the event triggered messages and the system will be improved.An improved two-dimensional bin packing algorithm is also presented to achieve the above-mentioned purpose with an extra benefit of better delay performance.展开更多
When the chaotic characteristics of manufacturing quality level are studied, it is not practical to use chaotic methods because of the low speed of calculating the correlation integral. The original algorithm used to ...When the chaotic characteristics of manufacturing quality level are studied, it is not practical to use chaotic methods because of the low speed of calculating the correlation integral. The original algorithm used to calculate the correlation integral is studied after a computer hardware upgrade. The result is that calculation of the correlation integral can be sped up only by improving the algorithm. This is accomplished by changing the original algorithm in which a single distance threshold-related correlation integral is obtained from one traversal of all distances between different vectors to a high-efficiency algorithm in which all of the distance threshold-related correlation integrals are obtained from one traversal of all of the distances between different vectors. For a time series with 3000 data points, this high-efficiency algorithm offers a 3.7-fold increase in speed over the original algorithm. Further study of the high-efficiency algorithm leads to the development of a super-high-efficiency algorithm, which is accomplished by changing the original and high-efficiency algorithms, in which the add-one operation of the Heaviside function is executed n times, such that the execution of the add-one operation occurs only once. The super-high-efficiency algorithm results in increases in the calculation speed by up to 109 times compared with the high-efficiency algorithm and by approximately 404 times compared with the original algorithm. The calculation speed of the super-high-efficiency algorithm is suitable for practical use with the chaotic method.展开更多
Vibration equations of time-varying system are transformed to the form which is suitable to precise integration algorithm.Precision analysis and computation efficiency of new algorithm are implemented.The following co...Vibration equations of time-varying system are transformed to the form which is suitable to precise integration algorithm.Precision analysis and computation efficiency of new algorithm are implemented.The following conclusions can be got.Choosing matrixes M,G and K is certainly flexible.We can place left side of nonlinear terms of vibration equations of time-varying system into right side of equations in precise integration algorithms.The key of transformation from vibration equations of time-varying system to first order differential equations is to form matrix H,which should be assured to be nonsingular.With suitable disposal,precision and computation efficiency of precise integration algorithms are greatly larger than those of general methods.展开更多
We propose a theoretical framework,based on the two-component Gross-Pitaevskii equation(GPE),for the investigation of vortex solitons(VSs)in hybrid atomic-molecular Bose-Einstein condensates under the action of the st...We propose a theoretical framework,based on the two-component Gross-Pitaevskii equation(GPE),for the investigation of vortex solitons(VSs)in hybrid atomic-molecular Bose-Einstein condensates under the action of the stimulated Raman-induced photoassociation and square-optical-lattice potential.Stationary solutions of the coupled GPE system are obtained by means of the imaginary-time integration,while the temporal dynamics are simulated using the fourth-order Runge-Kutta algorithm.The analysis reveals stable rhombus-shaped VS shapes with topological charges m=1 and 2 of the atomic component.The stability domains and spatial structure of these VSs are governed by three key parameters:the parametric-coupling strength(χ),atomicmolecular interaction strength(g_(12)),and the optical-lattice potential depth(V_(0)).By varyingχand g_(12),we demonstrate a structural transition where four-core rhombus-shaped VSs evolve into eight-core square-shaped modes,highlighting the nontrivial nonlinear dynamics of the system.This work establishes a connection between interactions of cold atoms and topologically structured matter waves in hybrid quantum systems.展开更多
基金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.
文摘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.
基金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.
文摘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].
文摘This paper presents an innovative and effective control strategy tailored for a deregulated,diversified energy system involving multiple interconnected area.Each area integrates a unique mix of power generation technologies:Area 1 combines thermal,hydro,and distributed generation;Area 2 utilizes a blend of thermal units,distributed solar technologies(DST),and hydro power;andThird control area hosts geothermal power station alongside thermal power generation unit and hydropower units.The suggested control system employs a multi-layered approach,featuring a blended methodology utilizing the Tilted Integral Derivative controller(TID)and the Fractional-Order Integral method to enhance performance and stability.The parameters of this hybrid TID-FOI controller are finely tuned using an advanced optimization method known as the Walrus Optimization Algorithm(WaOA).Performance analysis reveals that the combined TID-FOI controller significantly outperforms the TID and PID controllers when comparing their dynamic response across various system configurations.The study also incorporates investigation of redox flow batteries within the broader scope of energy storage applications to assess their impact on system performance.In addition,the research explores the controller’s effectiveness under different power exchange scenarios in a deregulated market,accounting for restrictions on generation ramp rates and governor hysteresis effects in dynamic control.To ensure the reliability and resilience of the presented methodology,the system transitions and develops across a broad range of varying parameters and stochastic load fluctuation.To wrap up,the study offers a pioneering control approach-a hybrid TID-FOI controller optimized via the Walrus Optimization Algorithm(WaOA)-designed for enhanced stability and performance in a complex,three-region hybrid energy system functioning within a deregulated framework.
基金Supported by Joint Fundation for Equipment Pre-research of Aerospace Science and Technology
文摘Traditional scheduling algorithms for avionics communication have the shortcoming of messages accumulation,the efficiency and reliability of the service can be improved by combining the distributed integrated modular avionics(DIMA)system with a time trigger mechanism.To further improve the utilization of system resources,the static scheduling algorithm of time triggered service is studied.By making the time trigger message schedule dispersedly,the stabilities of both the available time slots for the event triggered messages and the system will be improved.An improved two-dimensional bin packing algorithm is also presented to achieve the above-mentioned purpose with an extra benefit of better delay performance.
文摘When the chaotic characteristics of manufacturing quality level are studied, it is not practical to use chaotic methods because of the low speed of calculating the correlation integral. The original algorithm used to calculate the correlation integral is studied after a computer hardware upgrade. The result is that calculation of the correlation integral can be sped up only by improving the algorithm. This is accomplished by changing the original algorithm in which a single distance threshold-related correlation integral is obtained from one traversal of all distances between different vectors to a high-efficiency algorithm in which all of the distance threshold-related correlation integrals are obtained from one traversal of all of the distances between different vectors. For a time series with 3000 data points, this high-efficiency algorithm offers a 3.7-fold increase in speed over the original algorithm. Further study of the high-efficiency algorithm leads to the development of a super-high-efficiency algorithm, which is accomplished by changing the original and high-efficiency algorithms, in which the add-one operation of the Heaviside function is executed n times, such that the execution of the add-one operation occurs only once. The super-high-efficiency algorithm results in increases in the calculation speed by up to 109 times compared with the high-efficiency algorithm and by approximately 404 times compared with the original algorithm. The calculation speed of the super-high-efficiency algorithm is suitable for practical use with the chaotic method.
基金Project(50078006) supported by the National Natural Science Foundation of China
文摘Vibration equations of time-varying system are transformed to the form which is suitable to precise integration algorithm.Precision analysis and computation efficiency of new algorithm are implemented.The following conclusions can be got.Choosing matrixes M,G and K is certainly flexible.We can place left side of nonlinear terms of vibration equations of time-varying system into right side of equations in precise integration algorithms.The key of transformation from vibration equations of time-varying system to first order differential equations is to form matrix H,which should be assured to be nonsingular.With suitable disposal,precision and computation efficiency of precise integration algorithms are greatly larger than those of general methods.
基金supported by the National Natural Science Foundation of China(Grant No.62275075)the Natural Science Foundation of Hubei Soliton Research Association(Grant No.2025HBSRA09)+1 种基金joint supported by Hubei Provincial Natural Science Foundation and Xianning of China(Grant Nos.2025AFD401 and 2025AFD405)Israel Science Foundation(Grant No.1695/22).
文摘We propose a theoretical framework,based on the two-component Gross-Pitaevskii equation(GPE),for the investigation of vortex solitons(VSs)in hybrid atomic-molecular Bose-Einstein condensates under the action of the stimulated Raman-induced photoassociation and square-optical-lattice potential.Stationary solutions of the coupled GPE system are obtained by means of the imaginary-time integration,while the temporal dynamics are simulated using the fourth-order Runge-Kutta algorithm.The analysis reveals stable rhombus-shaped VS shapes with topological charges m=1 and 2 of the atomic component.The stability domains and spatial structure of these VSs are governed by three key parameters:the parametric-coupling strength(χ),atomicmolecular interaction strength(g_(12)),and the optical-lattice potential depth(V_(0)).By varyingχand g_(12),we demonstrate a structural transition where four-core rhombus-shaped VSs evolve into eight-core square-shaped modes,highlighting the nontrivial nonlinear dynamics of the system.This work establishes a connection between interactions of cold atoms and topologically structured matter waves in hybrid quantum systems.
