Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as ...Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.展开更多
Some new systems of exponentially general equations are introduced and investigated, which can be used to study the odd-order, non-positive and nonsymmetric exponentially boundary value problems. Some important and in...Some new systems of exponentially general equations are introduced and investigated, which can be used to study the odd-order, non-positive and nonsymmetric exponentially boundary value problems. Some important and interesting results such as Riesz-Frechet representation theorem, Lax-Milgram lemma and system of absolute values equations can be obtained as special cases. It is shown that the system of exponentially general equations is equivalent to nonlinear optimization problem. The auxiliary principle technique is used to prove the existence of a solution to the system of exponentially general equations. This technique is also used to suggest some new iterative methods for solving the system of the exponentially general equations. The convergence analysis of the proposed methods is analyzed. Ideas and techniques of this paper may stimulate further research.展开更多
In this study,a stable and robust interface-capturing method is developed to resolve inviscid,compressible two-fluid flows with general equation of state(EOS).The governing equations consist of mass conservation equat...In this study,a stable and robust interface-capturing method is developed to resolve inviscid,compressible two-fluid flows with general equation of state(EOS).The governing equations consist of mass conservation equation for each fluid,momentum and energy equations for mixture and an advection equation for volume fraction of one fluid component.Assumption of pressure equilibrium across an interface is used to close the model system.MUSCL-Hancock scheme is extended to construct input states for Riemann problems,whose solutions are calculated using generalized HLLC approximate Riemann solver.Adaptive mesh refinement(AMR)capability is built into hydrodynamic code.The resulting method has some advantages.First,it is very stable and robust,as the advection equation is handled properly.Second,general equation of state can model more materials than simple EOSs such as ideal and stiffened gas EOSs for example.In addition,AMR enables us to properly resolve flow features at disparate scales.Finally,this method is quite simple,time-efficient and easy to implement.展开更多
A new averaged general dynamic equation (GDE) for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and ...A new averaged general dynamic equation (GDE) for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and fluctuating coagulation. The equation is solved with the Taylor-series expansion moment method in a turbulent pipe flow. The experiments are performed. The numerical results of particle size distribu- tion correlate well with the experimental data. The results show that, for a turbulent nanoparticulate flow, a fluctuating coagulation term should be included in the averaged particle GDE. The larger the Schmidt number is and the lower the Reynolds number is, the smaller the value of ratio of particle diameter at the outlet to that at the inlet is. At the outlet, the particle number concentration increases from the near-wall region to the near-center region. The larger the Schmidt number is and the higher the Reynolds num- ber is, the larger the difference in particle number concentration between the near-wall region and near-center region is. Particle polydispersity increases from the near-center region to the near-wall region. The particles with a smaller Schmidt number and the flow with a higher Reynolds number show a higher polydispersity. The degree of particle polydispersity is higher considering fluctuating coagulation than that without considering fluctuating coagulation.展开更多
By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and non...By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and nonuniformity terms don't exist, the multisoliton solutions are found and the corresponding Painleve II type equation for the variable coefficient KdV equation is given.展开更多
The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be construc...The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be constructed by the solutions of original models if their solutions are well known, such as the standard constant coefficient KdV equation and the standard compound KdV--Burgers equation, and so on. Then any one of these variable-coefficient equations can be considered as an original model to obtain new variable-coefficient equations whose solutions can also be known by means of transformation relations between solutions of the resulting new variable-coefficient equations and the original equation.展开更多
In this paper,we study the hyperstability for the general linear equation f(ax+by)=Af(x)+Bf(y)in the setting of complete quasi-2-Banach spaces.We first extend the main fixed point result of Brzdek and Ciepliński(Acta...In this paper,we study the hyperstability for the general linear equation f(ax+by)=Af(x)+Bf(y)in the setting of complete quasi-2-Banach spaces.We first extend the main fixed point result of Brzdek and Ciepliński(Acta Mathematica Scientia,2018,38 B(2):377-390)to quasi-2-Banach spaces by defining an equivalent quasi-2-Banach space.Then we use this result to generalize the main results on the hyperstability for the general linear equation in quasi-2-Banach spaces.Our results improve and generalize many results of literature.展开更多
In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitatio...In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitations by parametric transformation. Bifurcation response equation and transition sets of sub-harmonic resonance with strong nonlinearity of general Van der Pol equation with parametric and external excitation are worked out.Besides, transition sets and bifurcation graphs are drawn to help to analysis the problems theoretically. Conclusions show that the transition sets of general and nonlinear Van der Pol equation with parametric and external excitations are more complex than those of general and nonlinear Van der Pol equation only with parametric excitation, which is helpful for the qualitative and quantitative reference for engineering and science applications.展开更多
In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three p...In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three parameters of standard linear solid, the forced vibration technique of beam is successfully used for PCL and PMMA specimens. The dynamical characteristics of viscoelastic Timoshenko beams, especially the damping properties, are derived from a considerable number of numerical computations. The analyses show that the viscosity of materials has great influence on dynamical characteristics of structures, especially on damping, and the standard linear solid model is the better one for describing the dynamic behavior of high viscous materials.展开更多
In this paper we consider one dimensional mean-field backward stochastic differential equations(BSDEs)under weak assumptions on the coefficient.Unlike[3],the generator of our mean-field BSDEs depends not only on the s...In this paper we consider one dimensional mean-field backward stochastic differential equations(BSDEs)under weak assumptions on the coefficient.Unlike[3],the generator of our mean-field BSDEs depends not only on the solution(Y,Z)but also on the law PY of Y.The first part of the paper is devoted to the existence and uniqueness of solutions in Lp,1<p≤2,where the monotonicity conditions are satisfied.Next,we show that if the generator/is uniformly continuous in(μ,y,z),uniformly with respect to(t,ω) and if the terminal valueξbelongs to Lp(Ω,F,P)with 1<p≤2,the mean-field BSDE has a unique Lp solution.展开更多
In this paper, we study the general structure of evolution equations of the AKNS eigenvalue problem q(x,t), r(x,t) with the spectrum varying asand AV BV CV are all positive or negative power polynomials of where q, r ...In this paper, we study the general structure of evolution equations of the AKNS eigenvalue problem q(x,t), r(x,t) with the spectrum varying asand AV BV CV are all positive or negative power polynomials of where q, r are not limited with any additional conditions at infinity.展开更多
This paper presents the conservation law of Einstein's General Relativity, in whichthe continuity equation of matter and gravitational field is implicitly derived from theequation of motion of matter. Although the...This paper presents the conservation law of Einstein's General Relativity, in whichthe continuity equation of matter and gravitational field is implicitly derived from theequation of motion of matter. Although the obtained energy-momentum tensor is thesame as the Landau- Lifshitz pseudotensor. the physical and conceptual foundation aredifferent. Two alternative methods to obtain the gravitational radiation are proposed inthis paper as well The radiation will be derived from the equation of motion of matterand geodesic equation separately.展开更多
This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the n...This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the nonlinear term is mere cubic.The new framework of convergence analysis consists of two steps.In the first step,by truncating the nonlinear term into a global Lipschitz function,an alternative numerical method is proposed and proved in a rigorous way to be convergent in the discrete L2 norm;followed in the second step,the maximum bound of the numerical solution of the alternative numerical method is obtained by using a lifting technique,as implies that the two numerical methods are the same one.Under our framework of convergence analysis,with neither any restriction on the grid ratio nor any requirement of the small initial value,we establish the error estimate of the proposed conservative Fourier pseudo-spectral method,while previous work requires the certain restriction for the focusing case.The error bound is proved to be of O(h^(r)+t^(2))with grid size h and time step t.In fact,the framework can be used to prove the unconditional convergence of many other Fourier pseudo-spectral methods for solving the nonlinear Schr¨odinger-type equations.Numerical results are conducted to indicate the accuracy and efficiency of the proposed method,and investigate the effect of the nonlinear term and initial data on the blow-up solution.展开更多
The objective of this paper is to study the oscillatory and asymptotic properties of the general mixed type third order neutral difference equation of the form △(aαn△2(xn+ bnxn-τ1+ cnxn+τ2)) + qnx^a n+1...The objective of this paper is to study the oscillatory and asymptotic properties of the general mixed type third order neutral difference equation of the form △(aαn△2(xn+ bnxn-τ1+ cnxn+τ2)) + qnx^a n+1-σ+ pnx^βn+1+=σ2 0,where {an}, {bn}, {cn}, {qn} and {pn} are positive real sequences, both α and β are ratios of odd positive integers, τ1, τ2, σ1 and σ2 are positive integers. We establish some sufficient conditions which ensure all solutions are either oscillatory or converge to zero.展开更多
This paper proposes a discrete-time robust control technique for an uncertain nonlinear system. The uncertainty mainly affects the system dynamics due to mismatched parameter variation which is bounded by a predefined...This paper proposes a discrete-time robust control technique for an uncertain nonlinear system. The uncertainty mainly affects the system dynamics due to mismatched parameter variation which is bounded by a predefined known function. In order to compensate the effect of uncertainty, a robust control input is derived by formulating an equivalent optimal control problem for a virtual nominal system with a modified costfunctional. To derive the stabilizing control law for a mismatched system, this paper introduces another control input named as virtual input. This virtual input is not applied directly to stabilize the uncertain system, rather it is used to define a sufficient condition. To solve the nonlinear optimal control problem, a discretetime general Hamilton-Jacobi-Bellman(DT-GHJB) equation is considered and it is approximated numerically through a neural network(NN) implementation. The approximated solution of DTGHJB is used to compute the suboptimal control input for the virtual system. The suboptimal inputs for the virtual system ensure the asymptotic stability of the closed-loop uncertain system. A numerical example is illustrated with simulation results to prove the efficacy of the proposed control algorithm.展开更多
Quantitative traits often underlie risk for complex diseases. Many studies collect multiple correlated quantitative phenotypes and perform univariate analyses on each of them respectively. However, this strategy may n...Quantitative traits often underlie risk for complex diseases. Many studies collect multiple correlated quantitative phenotypes and perform univariate analyses on each of them respectively. However, this strategy may not be powerful and has limitations to detect plei- otropic genes that may underlie correlated quantitative traits. In addition, testing multiple traits individually will exacerbate perplexing problem of multiple testing. In this study, generalized estimating equation 2 (GEE2) is applied to association mapping of two correlated quantitative traits. We suppose that a quantitative trait locus is located in a chromosome region that exerts pleiotropic effects on multiple quantitative traits. In that region, multiple SNPs are genotyped. Genotypes of these SNPs and the two quantitative traits affected by a causal SNP were simulated under various parameter values: residual correlation coefficient between two traits, causal SNP heritability, minor allele frequency of the causal SNP, extent of linkage disequilibrium with the causal SNP, and the test sample size. By power ana- lytical analyses, it is showed that the bivariate method is generally more powerful than the univariate method. This method is robust and yields false-positive rates close to the pre-set nominal significance level. Our real data analyses attested to the usefulness of the method.展开更多
In this paper,we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations.The state and co-state are approximated by the lowest order Raviart-...In this paper,we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations.The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions.We derive L2 and H−1-error estimates both for the control variable and the state variables.Finally,a numerical example is given to demonstrate the theoretical results.展开更多
The governing equations for heat transfer and fluid flow are often formulated in a general formfor the simplification of discretization and programming,which has achieved great success in thermal science and engineeri...The governing equations for heat transfer and fluid flow are often formulated in a general formfor the simplification of discretization and programming,which has achieved great success in thermal science and engineering.Based on the analysis of the popular general form of governing equations,we found that energy conservation cannot be guaranteed when specific heat capacity is not constant,which may lead to unreliable results.A new concept of generalized density is put forward,based on which a new general form of governing equations is proposed to guarantee energy conservation.A number of calculation examples have been employed to verify validation and feasibility.展开更多
Linear matrix equations are encountered in many systems and control applications.In this paper,we consider the general coupled matrix equations(including the generalized coupled Sylvester matrix equations as a specia...Linear matrix equations are encountered in many systems and control applications.In this paper,we consider the general coupled matrix equations(including the generalized coupled Sylvester matrix equations as a special case)l t=1EstYtFst = Gs,s = 1,2,···,l over the generalized reflexive matrix group(Y1,Y2,···,Yl).We derive an efcient gradient-iterative(GI) algorithm for fnding the generalized reflexive solution group of the general coupled matrix equations.Convergence analysis indicates that the algorithm always converges to the generalized reflexive solution group for any initial generalized reflexive matrix group(Y1(1),Y2(1),···,Yl(1)).Finally,numerical results are presented to test and illustrate the performance of the algorithm in terms of convergence,accuracy as well as the efciency.展开更多
We present a systematic and efficient Chebyshev spectral method using quasiinverse technique to directly solve the second order equation with the homogeneous Robin boundary conditions and the fourth order equation wit...We present a systematic and efficient Chebyshev spectral method using quasiinverse technique to directly solve the second order equation with the homogeneous Robin boundary conditions and the fourth order equation with the first and second boundary conditions.The key to the efficiency of the method is to multiply quasiinverse matrix on both sides of discrete systems,which leads to band structure systems.We can obtain high order accuracy with less computational cost.For multi-dimensional and more complicated linear elliptic PDEs,the advantage of this methodology is obvious.Numerical results indicate that the spectral accuracy is achieved and the proposed method is very efficient for 2-D high order problems.展开更多
文摘Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.
文摘Some new systems of exponentially general equations are introduced and investigated, which can be used to study the odd-order, non-positive and nonsymmetric exponentially boundary value problems. Some important and interesting results such as Riesz-Frechet representation theorem, Lax-Milgram lemma and system of absolute values equations can be obtained as special cases. It is shown that the system of exponentially general equations is equivalent to nonlinear optimization problem. The auxiliary principle technique is used to prove the existence of a solution to the system of exponentially general equations. This technique is also used to suggest some new iterative methods for solving the system of the exponentially general equations. The convergence analysis of the proposed methods is analyzed. Ideas and techniques of this paper may stimulate further research.
文摘In this study,a stable and robust interface-capturing method is developed to resolve inviscid,compressible two-fluid flows with general equation of state(EOS).The governing equations consist of mass conservation equation for each fluid,momentum and energy equations for mixture and an advection equation for volume fraction of one fluid component.Assumption of pressure equilibrium across an interface is used to close the model system.MUSCL-Hancock scheme is extended to construct input states for Riemann problems,whose solutions are calculated using generalized HLLC approximate Riemann solver.Adaptive mesh refinement(AMR)capability is built into hydrodynamic code.The resulting method has some advantages.First,it is very stable and robust,as the advection equation is handled properly.Second,general equation of state can model more materials than simple EOSs such as ideal and stiffened gas EOSs for example.In addition,AMR enables us to properly resolve flow features at disparate scales.Finally,this method is quite simple,time-efficient and easy to implement.
基金Project supported by the National Natural Science Foundation of China(No.11132008)
文摘A new averaged general dynamic equation (GDE) for nanoparticles in the turbulent flow is derived by considering the combined effect of convection, Brownian diffusion, turbulent diffusion, turbulent coagulation, and fluctuating coagulation. The equation is solved with the Taylor-series expansion moment method in a turbulent pipe flow. The experiments are performed. The numerical results of particle size distribu- tion correlate well with the experimental data. The results show that, for a turbulent nanoparticulate flow, a fluctuating coagulation term should be included in the averaged particle GDE. The larger the Schmidt number is and the lower the Reynolds number is, the smaller the value of ratio of particle diameter at the outlet to that at the inlet is. At the outlet, the particle number concentration increases from the near-wall region to the near-center region. The larger the Schmidt number is and the higher the Reynolds num- ber is, the larger the difference in particle number concentration between the near-wall region and near-center region is. Particle polydispersity increases from the near-center region to the near-wall region. The particles with a smaller Schmidt number and the flow with a higher Reynolds number show a higher polydispersity. The degree of particle polydispersity is higher considering fluctuating coagulation than that without considering fluctuating coagulation.
基金Supported by the Develop Programme Foundation of the National Basic research(G1 9990 3 2 80 1 )
文摘By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and nonuniformity terms don't exist, the multisoliton solutions are found and the corresponding Painleve II type equation for the variable coefficient KdV equation is given.
基金Project supported by the National Natural Science Foundation of China (Grant No. 10675065)the Scientific Research Fundof the Education Department of Zhejiang Province of China (Grant No. 20070979)
文摘The extended symmetry approach is used to study the general Korteweg-de Vries-type (KdV-type) equation. Several variable-coefficient equations are obtained. The solutions of these resulting equations can be constructed by the solutions of original models if their solutions are well known, such as the standard constant coefficient KdV equation and the standard compound KdV--Burgers equation, and so on. Then any one of these variable-coefficient equations can be considered as an original model to obtain new variable-coefficient equations whose solutions can also be known by means of transformation relations between solutions of the resulting new variable-coefficient equations and the original equation.
基金AISTDF,DST India for the research grant vide project No.CRD/2018/000017。
文摘In this paper,we study the hyperstability for the general linear equation f(ax+by)=Af(x)+Bf(y)in the setting of complete quasi-2-Banach spaces.We first extend the main fixed point result of Brzdek and Ciepliński(Acta Mathematica Scientia,2018,38 B(2):377-390)to quasi-2-Banach spaces by defining an equivalent quasi-2-Banach space.Then we use this result to generalize the main results on the hyperstability for the general linear equation in quasi-2-Banach spaces.Our results improve and generalize many results of literature.
基金Supported by the National Natural Science Foundation of China(11201118)
文摘In this paper a modified L-P method and multiple scale method are used to solve sub-harmonic resonance solutions of strong and nonlinear resonance of general Van der Pol equation with parametric and external excitations by parametric transformation. Bifurcation response equation and transition sets of sub-harmonic resonance with strong nonlinearity of general Van der Pol equation with parametric and external excitation are worked out.Besides, transition sets and bifurcation graphs are drawn to help to analysis the problems theoretically. Conclusions show that the transition sets of general and nonlinear Van der Pol equation with parametric and external excitations are more complex than those of general and nonlinear Van der Pol equation only with parametric excitation, which is helpful for the qualitative and quantitative reference for engineering and science applications.
文摘In this paper, a governing differential equation of viscoelastic Timoshenko beam including both extension and shear viscosity is developed in the time domain by direct method. To measure the complex moduli and three parameters of standard linear solid, the forced vibration technique of beam is successfully used for PCL and PMMA specimens. The dynamical characteristics of viscoelastic Timoshenko beams, especially the damping properties, are derived from a considerable number of numerical computations. The analyses show that the viscosity of materials has great influence on dynamical characteristics of structures, especially on damping, and the standard linear solid model is the better one for describing the dynamic behavior of high viscous materials.
基金supported in part by the NSFC(11222110,11871037)Shandong Province(JQ201202)+1 种基金NSFC-RS(11661130148,NA150344)111 Project(B12023)。
文摘In this paper we consider one dimensional mean-field backward stochastic differential equations(BSDEs)under weak assumptions on the coefficient.Unlike[3],the generator of our mean-field BSDEs depends not only on the solution(Y,Z)but also on the law PY of Y.The first part of the paper is devoted to the existence and uniqueness of solutions in Lp,1<p≤2,where the monotonicity conditions are satisfied.Next,we show that if the generator/is uniformly continuous in(μ,y,z),uniformly with respect to(t,ω) and if the terminal valueξbelongs to Lp(Ω,F,P)with 1<p≤2,the mean-field BSDE has a unique Lp solution.
基金The Projects Supported by the National Natural Science Foundation of China
文摘In this paper, we study the general structure of evolution equations of the AKNS eigenvalue problem q(x,t), r(x,t) with the spectrum varying asand AV BV CV are all positive or negative power polynomials of where q, r are not limited with any additional conditions at infinity.
文摘This paper presents the conservation law of Einstein's General Relativity, in whichthe continuity equation of matter and gravitational field is implicitly derived from theequation of motion of matter. Although the obtained energy-momentum tensor is thesame as the Landau- Lifshitz pseudotensor. the physical and conceptual foundation aredifferent. Two alternative methods to obtain the gravitational radiation are proposed inthis paper as well The radiation will be derived from the equation of motion of matterand geodesic equation separately.
基金Jialing Wang’s work is supported by the National Natural Science Foundation of China(Grant No.11801277)Tingchun Wang’s work is supported by the National Natural Science Foundation of China(Grant No.11571181)+1 种基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20171454)Qing Lan Project.Yushun Wang’s work is supported by the National Natural Science Foundation of China(Grant Nos.11771213 and 12171245).
文摘This paper aims to build a new framework of convergence analysis of conservative Fourier pseudo-spectral method for the general nonlinear Schr¨odinger equation in two dimensions,which is not restricted that the nonlinear term is mere cubic.The new framework of convergence analysis consists of two steps.In the first step,by truncating the nonlinear term into a global Lipschitz function,an alternative numerical method is proposed and proved in a rigorous way to be convergent in the discrete L2 norm;followed in the second step,the maximum bound of the numerical solution of the alternative numerical method is obtained by using a lifting technique,as implies that the two numerical methods are the same one.Under our framework of convergence analysis,with neither any restriction on the grid ratio nor any requirement of the small initial value,we establish the error estimate of the proposed conservative Fourier pseudo-spectral method,while previous work requires the certain restriction for the focusing case.The error bound is proved to be of O(h^(r)+t^(2))with grid size h and time step t.In fact,the framework can be used to prove the unconditional convergence of many other Fourier pseudo-spectral methods for solving the nonlinear Schr¨odinger-type equations.Numerical results are conducted to indicate the accuracy and efficiency of the proposed method,and investigate the effect of the nonlinear term and initial data on the blow-up solution.
基金Supported by Science Research Foundation of Guangxi Education Board under grant YB2014117
文摘The objective of this paper is to study the oscillatory and asymptotic properties of the general mixed type third order neutral difference equation of the form △(aαn△2(xn+ bnxn-τ1+ cnxn+τ2)) + qnx^a n+1-σ+ pnx^βn+1+=σ2 0,where {an}, {bn}, {cn}, {qn} and {pn} are positive real sequences, both α and β are ratios of odd positive integers, τ1, τ2, σ1 and σ2 are positive integers. We establish some sufficient conditions which ensure all solutions are either oscillatory or converge to zero.
文摘This paper proposes a discrete-time robust control technique for an uncertain nonlinear system. The uncertainty mainly affects the system dynamics due to mismatched parameter variation which is bounded by a predefined known function. In order to compensate the effect of uncertainty, a robust control input is derived by formulating an equivalent optimal control problem for a virtual nominal system with a modified costfunctional. To derive the stabilizing control law for a mismatched system, this paper introduces another control input named as virtual input. This virtual input is not applied directly to stabilize the uncertain system, rather it is used to define a sufficient condition. To solve the nonlinear optimal control problem, a discretetime general Hamilton-Jacobi-Bellman(DT-GHJB) equation is considered and it is approximated numerically through a neural network(NN) implementation. The approximated solution of DTGHJB is used to compute the suboptimal control input for the virtual system. The suboptimal inputs for the virtual system ensure the asymptotic stability of the closed-loop uncertain system. A numerical example is illustrated with simulation results to prove the efficacy of the proposed control algorithm.
基金supported by grants from the Natural Science Foundation of China (No.30600364,30470534,and 30230210)the NSFC-Canadian Institutes of Health Research (CIHR) Joint Health Research Initia-tive Proposal (No.30811120436)+3 种基金the NSFC/RGC Joint Research Scheme (No.30731160618)Shanghai Leading Academic Discipline Project (No.S30501)supported by grants from NIH (No.P50AR055081,R01AG026564,R01AR050496,and R01AR057049)the Dickson/Missouri endowment
文摘Quantitative traits often underlie risk for complex diseases. Many studies collect multiple correlated quantitative phenotypes and perform univariate analyses on each of them respectively. However, this strategy may not be powerful and has limitations to detect plei- otropic genes that may underlie correlated quantitative traits. In addition, testing multiple traits individually will exacerbate perplexing problem of multiple testing. In this study, generalized estimating equation 2 (GEE2) is applied to association mapping of two correlated quantitative traits. We suppose that a quantitative trait locus is located in a chromosome region that exerts pleiotropic effects on multiple quantitative traits. In that region, multiple SNPs are genotyped. Genotypes of these SNPs and the two quantitative traits affected by a causal SNP were simulated under various parameter values: residual correlation coefficient between two traits, causal SNP heritability, minor allele frequency of the causal SNP, extent of linkage disequilibrium with the causal SNP, and the test sample size. By power ana- lytical analyses, it is showed that the bivariate method is generally more powerful than the univariate method. This method is robust and yields false-positive rates close to the pre-set nominal significance level. Our real data analyses attested to the usefulness of the method.
基金supported by National Natural Science Foundation of China(Grant No.11526036)Scientific and Technological Developing Scheme of Jilin Province(Grant No.20160520108JH).
文摘In this paper,we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations.The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions.We derive L2 and H−1-error estimates both for the control variable and the state variables.Finally,a numerical example is given to demonstrate the theoretical results.
基金supported by the National Natural Science Foundation of China(No.51176204 and No.51134006),and the State Key Laboratory of Multiphase Flow in Power Engineering(Xi’an Jiaotong University).
文摘The governing equations for heat transfer and fluid flow are often formulated in a general formfor the simplification of discretization and programming,which has achieved great success in thermal science and engineering.Based on the analysis of the popular general form of governing equations,we found that energy conservation cannot be guaranteed when specific heat capacity is not constant,which may lead to unreliable results.A new concept of generalized density is put forward,based on which a new general form of governing equations is proposed to guarantee energy conservation.A number of calculation examples have been employed to verify validation and feasibility.
文摘Linear matrix equations are encountered in many systems and control applications.In this paper,we consider the general coupled matrix equations(including the generalized coupled Sylvester matrix equations as a special case)l t=1EstYtFst = Gs,s = 1,2,···,l over the generalized reflexive matrix group(Y1,Y2,···,Yl).We derive an efcient gradient-iterative(GI) algorithm for fnding the generalized reflexive solution group of the general coupled matrix equations.Convergence analysis indicates that the algorithm always converges to the generalized reflexive solution group for any initial generalized reflexive matrix group(Y1(1),Y2(1),···,Yl(1)).Finally,numerical results are presented to test and illustrate the performance of the algorithm in terms of convergence,accuracy as well as the efciency.
基金supported by the grants of National Natural Science Foundation of China(No.10731060,10801120)Chinese Universities Scientific Fund No.2010QNA3019.
文摘We present a systematic and efficient Chebyshev spectral method using quasiinverse technique to directly solve the second order equation with the homogeneous Robin boundary conditions and the fourth order equation with the first and second boundary conditions.The key to the efficiency of the method is to multiply quasiinverse matrix on both sides of discrete systems,which leads to band structure systems.We can obtain high order accuracy with less computational cost.For multi-dimensional and more complicated linear elliptic PDEs,the advantage of this methodology is obvious.Numerical results indicate that the spectral accuracy is achieved and the proposed method is very efficient for 2-D high order problems.